Can $\Lambda$CDM model reproduce MOND-like behavior?
  It is usually believed that MOND can describe the galactic rotational curves
with only baryonic matter and without any dark matter very well, while the
$\Lambda$CDM model is expected to have difficulty in reproducing MOND-like
behavior. Here, we use EAGLE's data to learn whether the $\Lambda$CDM model can
reproduce MOND-like behavior. EAGLE's simulation result clearly reproduces the
MOND-like behavior for $a_b\gtrapprox 10^{-12}\text{m/s}^2$ at $z=0$, although
the acceleration constant, $a_0$, is a little larger than the observational
data indicate. We find that $a_0$ increases with the redshift in a way
different from what Milgrom proposed ($a_0\propto H$). Therefore, while galaxy
rotation curves can be fitted by MOND's empirical function in the $\Lambda$CDM
model, there is no clear connection between $a_0$ and the Hubble constant. We
also find that $a_0$ at $z\gtrapprox 1$ is well separated from $a_0$ at $z=0$.
Once we have enough galaxies observed at high redshifts, we will be able to
rule out the modified gravity model based on MOND-like empirical function with
a z-independent $a_0$.
