This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

Both major galaxies in the Local Group host planar distributions of co-orbiting satellite galaxies, the Vast Polar Structure (VPOS) of the Milky Way and the Great Plane of Andromeda (GPoA). The $\Lambda$CDM cosmological model did not predict these features. However, according to three recent studies the properties of the GPoA and the flattening of the VPOS are common features among sub-halo based $\Lambda$CDM satellite systems, and the GPoA can be naturally explained by satellites being acquired along cold gas streams. We point out some methodological issues in these studies: either the selection of model satellites is different from that of the observed ones, or an incomplete set of observational constraints has been considered, or the observed satellite distribution is inconsistent with basic assumptions. Once these issues have been addressed, the conclusions are different: features like the VPOS and GPoA are very rare (each with probability $\lesssim 10^{-3}$, and combined probability $< 10^{-5}$) if satellites are selected from a $\Lambda$CDM simulation combined with semi-analytic modelling, and accretion along cold streams is no natural explanation of the GPoA. The origin of planar dwarf galaxy structures remains unexplained in the standard paradigm of galaxy formation.

Since its first formulation in 1983, Milgromian dynamics (MOND) has been very successful in predicting the gravitational potential of galaxies from the distribution of baryons alone, including general scaling relations and detailed rotation curves of large statistical samples of individual galaxies covering a large range of masses and sizes. Most predictions however rely on static models, and only a handful of N-body codes have been developed over the years to investigate the consequences of the Milgromian framework for the dynamics of complex evolving dynamical systems. In this work, we present a new Milgromian N-body code, which is a customized version of the RAMSES code (Teyssier 2002) and thus comes with all its features: it includes particles and gas dynamics, and importantly allows for high spatial resolution of complex systems due to the adaptive mesh refinement (AMR) technique. It further allows the direct comparison between Milgromian simulations and standard Newtonian simulations with dark matter particles. We provide basic tests of this customized code and demonstrate its performance by presenting N-body computations of dark-matter-free spherical equilibrium models as well as dark-matter-free disk galaxies in Milgromian dynamics.

Prompted by the recent successful predictions of the internal dynamics of Andromeda's satellite galaxies (McGaugh & Milgrom 2013a,b), we revisit the classical Milky Way dwarf spheroidal satellites Draco, Sculptor, Sextans, Carina, and Fornax in the framework of Milgromian dynamics (MOND). We use for the first time a Poisson solver with adaptive mesh refinement (AMR) in order to account simultaneously for the gravitational influence of the Milky Way and its satellites. This allows to rigorously model the important external field effect (EFE) of Milgromian dynamics, which can reduce the effective acceleration significantly. We make predictions on the dynamical mass-to-light ratio (Mdyn/L) expected to be measured by an observer who assumes Newtonian dynamics to be valid. We show that Milgromian dynamics predicts typical Mdyn/L ~ 10-50 Msun/Lsun. The results for the most luminous ones, Fornax and Sculptor, agree well with available velocity dispersion data. Moreover, the central power-law slopes of the dynamical masses agree exceedingly well with values inferred observationally from velocity dispersion measurements. The results for Sextans, Carina and Draco are low compared to usually quoted observational estimates, as already pointed out by Angus (2008). For Milgromian dynamics to survive further observational tests in these objects, one would thus need that either (a) previous observational findings based on velocity dispersion measurements have overestimated the dynamical mass due to, e.g., binaries and contaminant outliers, (b) the satellites are not in virial equilibrium due to the Milky Way tidal field, or (c) the specific theory used here does not describe the EFE correctly (e.g., the EFE could be practically negligible in some other theories), or a combination of (a)-(c).

The Local Group timing has been one of the first historical probes of the missing mass problem. Whilst modern cosmological probes indicate that pure baryonic dynamics is not sufficient on the largest scales, nearby galaxies and small galaxy groups persistently obey Milgrom's MOND law, which implies that dynamics at small scales is possibly entirely predicted by the baryons. Here, we investigate the Local Group timing in this context of Milgromian dynamics. Making use of the latest measured proper motions and radial velocities for Andromeda and the Magellanic clouds, we integrate backwards their orbits by making use of the Milgromian two-body equation of motion. We find that, with the currently measured proper motions and radial velocity of M31, MOND would imply that the Milky Way and M31 first moved apart via Hubble expansion after birth, but then necessarily got attracted again by the Milgromian gravitational attraction, and had a past fly-by encounter before coming to their present positions. This encounter would most probably have happened 7 to 11 Gyr ago (0.8<z<3). The absence of a dark matter halo and its associated dynamical friction is necessary for such a close encounter not to trigger a merger. Observational arguments which could exclude or favour such a past encounter would thus be very important in view of falsifying or vindicating Milgromian dynamics on the scale of the Local Group. Interestingly, the closest approach of the encounter is small enough (<55 kpc) to have severe consequences on the disk dynamics, including perhaps thick disk formation, and on the satellite systems of both galaxies. The ages of the satellite galaxies and of the young halo globular clusters, all of which form the vast polar structure around the Milky Way, are consistent with these objects having been born in this encounter.

Polar ring galaxies are ideal objects with which to study the three-dimensional shapes of galactic gravitational potentials since two rotation curves can be measured in two perpendicular planes. Observational studies have uncovered systematically larger rotation velocities in the extended polar rings than in the associated host galaxies. In the dark matter context, this can only be explained through dark halos that are systematically flattened along the polar rings. Here, we point out that these objects can also be used as very effective tests of gravity theories, such as those based on Milgromian dynamics (MOND). We run a set of polar ring models using both Milgromian and Newtonian dynamics to predict the expected shapes of the rotation curves in both planes, varying the total mass of the system, the mass of the ring with respect to the host, as well as the size of the hole at the center of the ring. We find that Milgromian dynamics not only naturally leads to rotation velocities being typically higher in the extended polar rings than in the hosts, as would be the case in Newtonian dynamics without dark matter, but that it also gets the shape and amplitude of velocities correct. Milgromian dynamics thus adequately explains this particular property of polar ring galaxies.