Issue 
A&A
Volume 636, April 2020



Article Number  A56  
Number of page(s)  10  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201936964  
Published online  17 April 2020 
A first attempt to differentiate between modified gravity and modified inertia with galaxy rotation curves
^{1}
CP 3Origins, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark
email: petersen@cp3.sdu.dk
^{2}
European Southern Observatory, KarlSchwarschildStrasse 2, 85748 Garching bei Munchen, Germany
email: LelliF@cardiff.ac.uk
^{3}
School of Physics and Astronomy, Cardiff University, Queens Buildings, The Parade, Cardiff CF24 3AA, UK
Received:
21
October
2019
Accepted:
19
December
2019
The phenomenology of modified Newtonian dynamics (MOND) on galaxy scales may point to more fundamental theories of either modified gravity (MG) or modified inertia (MI). In this paper, we test the applicability of the global deepMOND parameter Q which is predicted to vary at the 10% level between MG and MI theories. Using mockobserved analytical models of disk galaxies, we investigate several observational uncertainties, establish a set of quality requirements for actual galaxies, and derive systematic corrections in the determination of Q. Implementing our quality requirements to the SPARC database yields 15 galaxies, which are close enough to the deepMOND regime as well as having rotation curves that are sufficiently extended and sampled. For these galaxies, the average and median values of Q seem to favor MG theories, albeit both MG and MI predictions are in agreement with the data within 1.5σ. Improved precision in the determination of Q can be obtained by measuring extended and finelysampled rotation curves for a significant sample of extremely lowsurfacebrightness galaxies.
Key words: galaxies: kinematics and dynamics / galaxies: structure / dark matter
© ESO 2020
1. Introduction
The missing mass problem is established by a host of astronomical observations, including the dynamical behavior of galaxy clusters (Zwicky 1933; Clowe et al. 2006), the rotation curves of disk galaxies (Rubin et al. 1978; Bosma 1978; van Albada et al. 1985), and the properties of the cosmic microwave background (Planck Collaboration XIII 2016; Komatsu et al. 2003). Solutions to the missing mass problem have been proposed in the form of either some unobserved mass (dark matter) or a modification to the classical laws of dynamics. The former option has been vastly explored and led to the standard cosmological model Λ cold dark matter (ΛCDM), which provides a good description of the Universe on large scales but has some problems on galaxy scales (Bullock & BoylanKolchin 2017). The latter option led to the development of modified Newtonian dynamics (MOND; Milgrom 1983), in which the classical laws of dynamics are modified at small accelerations (see e.g., Sanders & McGaugh 2002; Famaey & McGaugh 2012; Milgrom 2014 for reviews of MOND). MOND is primarily motivated by the close relationship between the baryonic matter and the observed dynamical behavior of galaxies (Tully & Fisher 1977; Faber & Jackson 1976; McGaugh et al. 2000, 2016; Lelli et al. 2017) as well as the existence of a characteristic acceleration scale () below which the dark matter effect appears (Sanders 1990; McGaugh 2004; McGaugh et al. 2016).
The MOND phenomenology can be interpreted either as a modification of inertia (MI) or a modification of gravity (MG). In the context of the nonrelativistic action of the theory, MI and MG differ in that the former alters the kinetic term, whereas the latter alters the potential term (Milgrom 2002, 2005). The first fullfledged MG theory was proposed by Bekenstein & Milgrom (1984) and it is commonly known as AQUAL due to the aquadratic form of the Lagrangian. Another MG theory was proposed by Milgrom (2009) and is dubbed quasilinear MOND (QMOND). On the other hand MI theories has proven to be more difficult to build. Milgrom (1994), however, demonstrated that MI theories must be strongly nonlocal and that for purely circular orbits, MI leads to the algebraic MOND equation
where g_{obs} is the observed kinematic acceleration, g_{N} is the classical Newtonian gravitational acceleration due to the baryons, μ is an arbitrary interpolation function that reproduces the Newtonian regime at large accelerations, that is μ → 1 for g_{obs}/a_{0} → ∞, and asymptotic flat rotation curves at low accelerations, that is for g_{obs}/a_{0} → 0. In the MG context, instead, Eq. (1) is only valid for highly symmetric mass distributions, such as in the case of spherical geometry (Bekenstein & Milgrom 1984).
Both MI and MG have struggled on larger scales, for example in accounting for the entire gravitational anomaly in galaxy clusters (Sanders 2003) and cosmological observations (Skordis et al. 2005; Dodelson & Liguori 2006; Dodelson 2011). However, the two theories remain a popular source of inspiration for model building (e.g., see Chashchina et al. 2017; Edmonds et al. 2017; Dai & Lu 2017; Cai et al. 2018; Berezhiani et al. 2018; Costa et al. 2019) due to their ability to provide an intuitive explanation for galactic data. Recently MI has been debated in regards to galactic scales, with McGaugh et al. (2016) and Li et al. (2018), presenting evidence in favor of MI and Petersen & Frandsen (2017), and Frandsen & Petersen (2018), Petersen (2019) presenting evidence in opposition with MI.
In principle, one could distinguish between MG and MI by performing detailed rotationcurve fits under the two different prescriptions and compare the corresponding residuals. This approach, however, relies on the knowledge of the galaxy distance and disk inclination, which often dominate the error budget (e.g., Lelli et al. 2017; Li et al. 2018) and may well exceed the expected difference of ∼10% expected between MG and MI (Brada & Milgrom 1995).
In this paper, we investigate the possibility of differentiating between MI and MG theories using a parameter proposed by Milgrom (2012)
where v_{obs} is the observed rotation velocity, Σ_{bar} is the baryonic mass surface density (gas plus stars), m_{bar} is the total baryonic mass, and v_{∞} is the asymptotically flat velocity obtained at r → ∞. Milgrom (2012) shows that for diskonly galaxies everywhere in the deep MONDian limit (DML, maximum acceleration ≪a_{0}), MI predicts Q^{(MI)} ≈ 0.73 ± 0.01 whereas MG predicts . Thus, there is potential for differentiating between the two model classes by measuring Q in low density, low acceleration galaxies. The result Q^{(MG)} = 2/3 has general validity and relies only on the form of the deepMOND virial relation, which is the same for MG theories like AQUAL and QMOND (Zhao & Famaey 2010). The result for MI, instead, is not general but depend on the adopted mass distribution. Milgrom (2012), however, finds that Q^{(MI)} is very close to 0.73 for a large set of realistic mass density profiles for disk galaxies. A major advantage of the Q parameter is that uncertainties due to galaxy distance and disk inclination cancel out due to the normalization factors m_{bar} and v_{∞}. The drawback, however, is that both v_{∞} and the integral are defined for r → ∞. Thus, one needs to quantify the systematic effects introduced by measuring Q in finitesize galaxies with finite spatial resolution.
In Sect. 2, we create a set of mock galaxies to quantify systematic effects in the estimate of Q. In Sect. 3, we present a first attempt at measuring Q for galaxies in the SPARC database. Finally, in Sect. 4, we draw our conclusions and discuss future prospects.
2. Mock galaxies
In the following, we create a set of mock galaxies in order to quantify a number of systematic uncertainties in the derivation of Q: (1) The effect of deviating from the theoretical deep MOND limit (DML), (2) the effect of the finitesize of galaxy disks and (3) the effect of the finite spatial resolution. We will show that these uncertainties can be accounted for by using the following equation:
where Q^{(m)} and Q^{(c)} are the measured and corrected Q values, respectively. ΔQ^{(DML)}, ΔQ^{(range)} and ΔQ^{(res.)} are correction factors that can be applied as long as max[g_{obs}] ≲ 0.4a_{0} (acceleration scale requirement), max[r] ≳ 3r_{d} (sampled range requirement) and a spacing between sampled points smaller than 1.3r_{d} (resolution requirement).
For the mock data, let the gas and stellar mass densities be represented by a single disk profile characterized by a central disk mass density and a disk scale length. The mock data are generated by adding Gaussian noise around each point in the (r_{d}, Σ_{d})plane of data from the SPARC database and subsequently removing galaxies with negative values. On top of this – to make sure that the deep MONDian regime is sampled well enough – 500 points are randomly generated at decreasing acceleration scales. Figure 1 illustrates both the original SPARC data (red) and the mock data (black) in the (r_{d}, Σ_{d})plane. For the analysis of mock galaxies, the Kuzmin disk is considered as a case study. For this mass distribution, MG has an exact analytic solution (Brada & Milgrom 1995), so the exact difference between MI and MG can be extrapolated beyond the deep MONDian regime. The surface mass density of the Kuzmin disk is given by
Fig. 1.
Distribution of real (red) and mock (black) galaxies in the r_{d} − Σ_{d} plane. The blue line denotes the scale of MOND in terms of the characteristic surface density . 
where Σ_{0} denotes the central surface density and r_{d} denotes the disk scale length. Additional mass distributions are considered in Appendix 2, using the BradaMilgrom approximation for MG (Eq. (25) in Brada & Milgrom 1995) which is valid for AQUALlike theories. All these mass distributions give similar results as the Kuzmin disk.
In order to asses when the DML is reached, an interpolation function must be specified for both MI and MG in the BekensteinMilgrom formulation. For example, for a step function, the DML is reached as soon as the maximum acceleration is less than a_{0}, but this is clearly an unrealistic situation. Thus, we use the (inverse) interpolation function corresponding to the radial acceleration relation (RAR, McGaugh et al. 2016);
2.1. The deep MONDian limit (DML)
To gauge when the DML is reached, Q is calculated from the mock data within bins of maximum total acceleration (max(g_{obs})) for MI and MG. This is shown in the left panel of Fig. 2, from which it is clear that the expected DML value of Q is reached for max(g_{obs})≲0.01a_{0}. These are extremely low accelerations: known disk galaxies are never entirely in such a deep MOND regime. However, the predictions for Q for MI and MG drift in a similar way as max(g_{obs}) increases. As long as the drifts of MI and MG are approximately equal, the measured Q can – in the context of differentiating between MI and MG – be rescaled such that galaxies above the deep MONDian regime yield “correct” DML values of Q. Denote the drifts in MI and MG viz
Fig. 2.
Left panel: Q calculated for both MI and MG for a Kuzmin disk using mock data binned into acceleration bins. The gray lines denote the deep MONDian limits for MI and MG, respectively. Right panel: ΔQ^{MI} − ΔQ^{MG} as a function of acceleration scale. 
where and denotes the exact values of Q. As long as ΔQ^{MI} − ΔQ^{MG} ≃ 0, Q can be rescaled to the deep MONDian value with either ΔQ^{MI} or ΔQ^{MG}. The right panel of Fig. 2 illustrates the difference ΔQ^{MI} − ΔQ^{MG}. ΔQ^{MI} − ΔQ^{MG}≲0.01 ≃ const for max_{r}(g_{obs})≲0.4a_{0}, meaning that the measured Q^{m} can be corrected with ΔQ^{(DML)} = ΔQ^{MI} ≃ ΔQ^{MG} for galaxies which fulfill max(g_{obs})≲0.4a_{0}. Table 1 lists ΔQ^{(DML)} for different acceleration bins.
Acceleration scale corrections.
2.2. Range of sampling
To gauge the effect of the finite extension of galaxy rotation curves and mass surface density profiles, Q is computed at decreasing multiples of the disk scale length (r_{cut} = N_{range}r_{d}). The left panel of Fig. 3 illustrates the arithmetic mean of Q of the mock data for which max[g_{obs}]≤0.01a_{0} as a function of decreasing N_{range}. Q is shown for both MI and MG calculated in two different ways: (i) with v_{∞} fixed at r = 500 kpc to isolate the effect of neglecting part of the integral in the numerator of Q and (ii) v_{∞} approximated by v_{tot}(r_{cut}) such that it provides a better representation of the observational situation. From the left panel of Fig. 3 it is clear that all curves approach the theoretical values at increasing r_{cut}. At small ranges Q decreases for both approximations of v_{∞}. However, with the moving approximation of v_{∞}, Q decreases significantly less than the fixed one (for both MI and MG). This is a lucky occurrence indicating that the underestimate of Q from cutting down the integral is somewhat counterbalanced by the underestimate of v_{∞} from the finite size of galaxy disks. Indeed, for galaxies with max[g_{obs}]< 0.01a_{0}, the rotation curve is still slowly rising up to 500 kpc, so v_{obs}(r_{cut}) < v_{∞}. Hence, the measured value of Q is closer to the theoretical value in the actual observational situation, since we can only estimate v_{flat} near r_{cut}.
Fig. 3.
Left panel: arithmetic mean of Q for all mock galaxies with max[g_{obs}]≤0.01a_{0} as a function of radial range in units of r_{d}. Black and blue lines represent MI and MG, respectively, with v_{∞} ≃ v_{flat} = v_{tot}(500 kpc). Cyan and brown lines represent MI and MG, respectively, with v_{∞} ≃ v_{flat} = v_{tot}(r_{cut}), where r_{cut} denotes the radii at which the integrals are cut. The inset is a magnification of N_{range} = [0, 10]. Right panel: ΔQ^{MI} − ΔQ^{MG} as a function of radial range. The error bars are smaller than the points. 
From the left panel of Fig. 3 it is clear that MI and MG drift similarly as the range is decreased. As long as the drifts of MI and MG are approximately equal, the measured Q can be corrected similarly to the case of the acceleration scale in the previous section. The right panel of Fig. 3 shows ΔQ^{MI} − ΔQ^{MG} in the context of varying range. ΔQ^{MI} − ΔQ^{MG}≲0.01 ≃ const for N_{range} ≳ 3, meaning that the measured Q can be corrected with ΔQ^{(range)} = ΔQ^{MI} ≃ ΔQ^{MG} for galaxies which fulfill max[r] ≳ 3r_{d}. Table 2 lists ΔQ^{(range)} for different ranges. The approach to the theoretical limit for increasing range is slow for the Kuzmin disk. For this reason ΔQ^{(range)} = 0 is only reached for N_{range} > 55 in this case. The corrections in Table 2 isolate the range effect from the acceleration effect. To mimick the observational situation, however, we have repeated the same exercise considering galaxies with max(g_{obs}) < 0.15a_{0}, applying the corrections in Table 1. We find that the corresponding corrections due to N_{range} are very similar to those given in Table 2.
Range corrections.
2.3. Sampling resolution
In order to gauge the effect of the spatial resolution of the observations, Q is calculated as a function of the spacing between sampled points in units of r_{d}. The sampling for each galaxy is performed in the range r_{j} ∈ [0.01, 200]r_{d}. In order to compute Q from a discrete set of points, the integral of Eq. (2) is discretized viz
with
where v_{f} ≃ v_{∞} being the asymptotically flat velocity which is determined as the arithmetic mean between the chain of points that are within 5% of the arithmetic mean of the two outermost points (similarly to Lelli et al. 2016). If the third outermost point is not within 5% of the arithmetic mean of the two outermost points, the galaxy is not assigned an asymptotically flat velocity.
The left panel of Fig. 4 shows the arithmetic mean of Q calculated via Eq. (7) – as predicted by MI (black) and MG (blue) – of the mock data for which max(g_{obs})≤0.01a_{0} as a function of the spacing between the sampled points in units of r_{d}. As expected, Q approaches the theoretical value for extremely small spacings (r < 0.1r_{d}) and progressively increases for larger spacing between the sampled points.
Fig. 4.
Left panel: arithmetic mean of Q for all mock galaxies with max[g_{obs}]≤0.01a_{0} as a function of spacing between sampled points. Black and blue lines represent MI and MG, respectively. Right panel: ΔQ^{MI} − ΔQ^{MG} as a function of spacing between sampled points. 
From the right panel of Fig. 4 it is clear that the increase for MI and MG are similar, meaning that – similarly to the acceleration scale and range cases – the measured Q can be rescaled such that galaxies with spacings larger than ∼0.5r_{d} yield Q approximately equal to that calculated with infinitely small spacing. ΔQ^{MI} − ΔQ^{MG}≲0.01 ≃ const for spacing ≳1.3r_{d}, meaning that the measured Q can be corrected with ΔQ^{(res.)} = ΔQ^{MI} ≃ ΔQ^{MG} for galaxies with spacing ≲1.3r_{d}. Table 3 lists ΔQ^{(res.)} for different spacing between sampled points.
Spacing corrections.
3. Application to SPARC galaxies
We will now use actual galaxies from the SPARC database (Lelli et al. 2016) to calculate Q. The SPARC database consists of 175 rotationally supported galaxies spanning broad dynamic ranges in stellar mass, surface brightness and gas fraction. The database provide observed rotational velocities (v_{obs}), along with the associated uncertainties (δv_{obs}), as well as distance (D) and inclination (i) measurements for each galaxy. The database also provides the rotational velocities due to the baryonic components, which are computed using the Spitzer [3.6] surface brightness profile for the stars and the HI surface density profile for the gas. Finally, the database provides the central surface brightness (μ_{d}) and exponential scale length (r_{d}) of the stellar disk (see Fig. 1). In this paper, however, we estimated the scale length of the baryonic disk (gas plus stars) by fitting an exponential profile to the total surface mass density profile of both stars and gas (see Appendix B).
Following Lelli et al. (2017), 12 galaxies are discarded from the analysis based on a quality flag (see Lelli et al. 2016). We note that in this work we keep faceon galaxies (i < 30°) because Q does not depend on i.
The baryonic velocity is computed viz (recall we are only considering diskonly galaxies)
where v_{g} is the gas contribution, v_{d} is the stellar contribution, and Υ_{d} is the stellar masstolight ratio. In line with McGaugh & Schombert (2014) and Lelli et al. (2016) is taken with a 25% uncertainty.
Q is discretized as shown in Eq. (7), however with Σ_{d} → μ_{d}Υ^{d} + 1.33Σ_{HI}. Requiring that there be a measured gas profile (approximated by the HI profile) and a well defined value for v_{f} on top of the accelerations scale, range and resolution requirements yields 15 galaxies from the SPARC database. In Appendix B, we present the rotation curves for these galaxies as well as a table detailing Q^{(m)}, Q^{(c)}, the acceleration scale, the range and the spacing for each galaxy. Figure 5 shows Q as a function of the acceleration scale (max[g_{obs}]), the range and the resolution. The left column of the figure shows Q^{(m)} whereas the right column shows Q^{(c)}. The averages are calculated viz
Fig. 5.
Q for individual galaxies as a function of max[g_{obs}] (first row), (second row) and spacing (third row). These 15 galaxies are deep enough in the MOND regime and have sufficiently extended and finely sampled rotation curves (see Sect. 2). The gray line shows the predicted Q from MOND modified gravity, while the gray band shows the ranges of values expected in MOND modified inertia. 
where N = ∑_{j ∈ gal}1, with gal denoting the relevant galaxies, is the number of measurements from a galaxy. The median of Q from individual galaxies can be written
From all three rows there is a slight tendency of the average Q decreasing as a the requirements are fulfilled to a higher degree. Let ⟨Q^{(m)}⟩ denote a list of the arithmetic mean and the median of Q^{m}, respectively. From the galaxies shown in Fig. 5 ⟨Q^{m}⟩ is given by
For the corrected Q
Comparing Eqs. (12) and (13) it is clear that (i) the arithmetic mean and median yield similar results and (ii) the sum of corrections have a very little impact on the two measures of the average Q, meaning that the corrections approximately cancel out.
Both the mean and median values of Q^{m} and Q^{c} are consistent with the MOND predictions within 1.5 sigma. With the current level of accuracy, we cannot distinguish between MG and MI, albeit there is a slight preference for the lower MG value. We also note that individual galaxies can significantly deviate from the predicted Q values (at more than 5σ), but this is likely due to the fact that the errorbars on a single object cannot possibility account for all the systematic uncertainties. A larger galaxy sample would have allowed us to determine whether Q converges toward a characteristic value and to rigorously estimate the corresponding error. We expect that a sample of ∼100−200 galaxies satisfying our quality criteria should be sufficient to carry out this experiment. This may be possible in the near future thanks to large HI surveys with the Square Kilometer Array (SKA) and its pathfinders (e.g., Duffy et al. 2012).
The effect of the stellar masstolight ratio on the measured Q is relatively minor in these 15 galaxies. If we assume rather than we obtain a slightly different set of galaxies with different corrections since the value of r_{d} of each galaxy slightly changes. Then Eqs. (12) and (13) become
These values are consistent with the previous ones within the errors, but they are systematically higher by ∼0.02−0.04 and become closer to the predicted value from MI.
In Li et al. (2018) the RAR is fitted to individual galaxies yielding Υ^{d} values that are maximally favorable for MI. The Q values corresponding to Eqs. (12) and (13) are in this case
The Q values are very close to the case of . Intriguingly, even when we adopt M/L values that are maximally favorable for MI, the resulting Q values remain closer to the predictions of MG.
Lastly we note that the analysis of Q conducted in this paper is closely related to the fractions of accelerations investigated in Frandsen & Petersen (2018) and Petersen (2019) since both quantities involve ratios between velocities and benefit from the cancellation of systematic uncertainties (inclination angle, distance and partially the mass to light ratio). However, despite the apparent similarities, the finer details in calculating Q and the fractions of accelerations in Frandsen & Petersen (2018) and Petersen (2019) lead to a different analysis and in the end the results cannot be compared.
4. Summary and conclusions
In this paper, we present a first attempt to differentiate between MOND modified inertia (MI) and MOND modified gravity (MG) using galactic rotation curve data. In particular, we investigate defined in Milgrom (2012), for which MI predicts Q = 0.73 ± 0.01 and MG predicts for diskonly galaxies everywhere in the deep MONDian regime. We compare theoretical predictions to data from the SPARC database. To do so, we thoroughly discuss the impact of systematic uncertainties in the value of Q via investigating a set of mock galaxies. Specifically we consider the impact on Q of the acceleration scale (max[g_{obs}]), the spacing between sampled points (resolution) and the range of sampling (max(r)). We find that the systematic uncertainties in Q can be approximately accounted for as long as max[g_{obs}] ≲ 0.4a_{0} (acceleration scale requirement), max[r] ≳ 3r_{d} (sampled range requirement) and the spacing between points is ≲1.3r_{d} (resolution requirement). Imposing these criteria on the SPARC database leaves 15 galaxies. Before correcting for systematic effects, the arithmetic mean and median of Q is given by ⟨Q^{(m)}⟩={0.63 ± 0.09, 0.64 ± 0.07}, respectively. After correction ⟨Q^{(c)}⟩ ≈ {0.63 ± 0.08, 0.65 ± 0.06}. From ⟨Q^{(m)}⟩ and ⟨Q^{(c)}⟩ several things can be noted; (i) the arithmetic mean and median yield similar predictions, (ii) these measurements line up closely with the predictions of MG both before and after correction, although the prediction of MI is still within 1.5σ.
Future HI surveys with the SKA and its pathfinders are expected to provide HI rotation curves for thousands of galaxies, which will allow us to achive more solid estimates of Q and push down the statistical errors.
Acknowledgments
We thank the referee for a constructive report that improved the clarity of our paper. JP also thanks ESO for the hospitality that made this collaboration possible as well as the partial funding from The Council For Independent Research, grant number DFF 610800623. The CP3Origins center is partially funded by the Danish National Research Foundation, grant number DNRF90.
References
 Bekenstein, J., & Milgrom, M. 1984, ApJ, 286, 7 [Google Scholar]
 Berezhiani, L., Famaey, B., & Khoury, J. 2018, JCAP, 2018 [Google Scholar]
 Bosma, A. 1978, PhD Thesis, Groningen Univ. [Google Scholar]
 Brada, R., & Milgrom, M. 1995, MNRAS, 276, 453 [Google Scholar]
 Bullock, J. S., & BoylanKolchin, M. 2017, Annu. Rev. Astron. Astrophys., 55, 343 [Google Scholar]
 Cai, R.G., Liu, T.B., & Wang, S.J. 2018, Phys. Rev. D, 97, 023027 [NASA ADS] [CrossRef] [Google Scholar]
 Chashchina, O., Foot, R., & Silagadze, Z. 2017, Phys. Rev. D, 95, 023009 [NASA ADS] [CrossRef] [Google Scholar]
 Clowe, D., Bradac, M., Gonzalez, A. H., et al. 2006, ApJ, 648, L109 [NASA ADS] [CrossRef] [Google Scholar]
 Costa, R., Franzmann, G., & Pereira, J. P. 2019, ArXiv eprints [arXiv:1904.07321] [Google Scholar]
 Dai, D.C., & Lu, C. 2017, Phys. Rev. D, 96, 124016 [NASA ADS] [CrossRef] [Google Scholar]
 Dodelson, S. 2011, Int. J. Mod. Phys. D, 20, 2749 [NASA ADS] [CrossRef] [Google Scholar]
 Dodelson, S., & Liguori, M. 2006, Phys. Rev. Lett., 97, 231301 [NASA ADS] [CrossRef] [Google Scholar]
 Duffy, A. R., Moss, A., & StaveleySmith, L. 2012, PASA, 29, 202 [NASA ADS] [CrossRef] [Google Scholar]
 Edmonds, D., Farrah, D., Minic, D., Ng, Y. J., & Takeuchi, T. 2017, Int. J. Mod. Phys. D, 27, 1830001 [NASA ADS] [CrossRef] [Google Scholar]
 Faber, S. M., & Jackson, R. E. 1976, ApJ, 204, 668 [NASA ADS] [CrossRef] [Google Scholar]
 Famaey, B., & McGaugh, S. 2012, Liv. Rev. Rel., 15, 10 [Google Scholar]
 Frandsen, M. T., & Petersen, J. 2018, ArXiv eprints [arXiv:1805.10706] [Google Scholar]
 Komatsu, E., Kogut, A., Nolta, M. R., et al. 2003, ApJS, 148, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Lelli, F., McGaugh, S. S., & Schombert, J. M. 2016, AJ, 152, 157 [Google Scholar]
 Lelli, F., McGaugh, S. S., Schombert, J. M., & Pawlowski, M. S. 2017, ApJ, 836, 152 [Google Scholar]
 Li, P., Lelli, F., McGaugh, S., & Schombert, J. 2018, A&A, 615, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 McGaugh, S. S. 2004, ApJ, 609, 652 [NASA ADS] [CrossRef] [Google Scholar]
 McGaugh, S., & Schombert, J. 2014, AJ, 148 [Google Scholar]
 McGaugh, S. S., Schombert, J. M., Bothun, G. D., & de Blok, W. J. G. 2000, ApJ, 533, L99 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 McGaugh, S., Lelli, F., & Schombert, J. 2016, Phys. Rev. Lett., 117, 201101 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 1983, ApJ, 270, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 1994, Ann. Phys., 229, 384 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 2002, New Astron. Rev., 46, 741 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 2005, EAS Publ. Ser., 20, 217 [CrossRef] [Google Scholar]
 Milgrom, M. 2009, Phys. Rev. D, 80, 123536 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 2012, Phys. Rev. Lett., 109, 251103 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 2014, Scholarpedia, 9, 31410 [NASA ADS] [CrossRef] [Google Scholar]
 Petersen, J. 2019, ArXiv eprints [arXiv:1906.09798] [Google Scholar]
 Petersen, J., & Frandsen, M. T. 2017, ArXiv eprints [arXiv:1710.03096] [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rubin, Jr., V. C., Ford, W. K., & Thonnard, N. 1978, ApJ, 225, L107 [NASA ADS] [CrossRef] [Google Scholar]
 Sanders, R. H. 1990, A&ARv, 2, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Sanders, R. H. 2003, MNRAS, 342, 901 [Google Scholar]
 Sanders, R. H., & McGaugh, S. S. 2002, Annu. Rev. Astron. Astrophys., 40, 263 [Google Scholar]
 Skordis, C., Mota, D. F., Ferreira, P. G., & Boehm, C. 2005, Phys. Rev. Lett., 96, 011301 [Google Scholar]
 Tully, R. B., & Fisher, J. R. 1977, A&A, 54, 661 [NASA ADS] [Google Scholar]
 van Albada, T. S., Bahcall, J. N., Begeman, K., & Sancisi, R. 1985, ApJ, 295, 305 [NASA ADS] [CrossRef] [Google Scholar]
 Zhao, H., & Famaey, B. 2010, Phys. Rev. D, 81, 087304 [NASA ADS] [CrossRef] [Google Scholar]
 Zwicky, F. 1933, Helv. Phys. Acta, 6, 110 [Gen. Rel. Grav., 41, 207 (2009)] [NASA ADS] [Google Scholar]
Appendix A: Mock galaxies for different mass distributions
Here we consider the approach to the DML for the different mass distributions considered in Milgrom (2012);
Exponential disk:
Double exponential disk:
Finite galaxy disk:
with
being hypergeometric functions. For the MG calculation, we adopt the approximation from Brada & Milgrom (1995). The results are shown in Fig. A.1.
Fig. A.1.
Q calculated from MOND modified inertia and MOND modified gravity in the Brada & Milgrom (1995) approximation from the mock data binned into acceleration bins. Each panel represents a different mass distribution (denoted on the panel). See Eqs. (A.1) and (A.2) for details. 
Appendix B: Galaxy set
Here we show the mass models for the 15 galaxies (see Fig. B.1) that satisfy our quality criteria: max(g_{obs}) < 0.4a_{0}, max(r) > 3r_{d}, and rotationcurve spacing smaller than 1.3r_{d}. The properties of these 15 galaxies are given in Table B.1.
Fig. B.1.
Rotation curves of galaxies in the SPARC database after imposing max[g_{obs}]≤0.4a_{0}, and spacing ≤1.3r_{d} as well as requiring that there be a gas profile and a welldefined approximation of v_{∞}. Black is v_{obs}, blue is v_{bar}, cyan is v_{g}, red is v_{d} and brown denotes the 1σ region for v_{f}. The brown band illustrates the points used to compute v_{f}. 
Galaxy set.
All Tables
All Figures
Fig. 1.
Distribution of real (red) and mock (black) galaxies in the r_{d} − Σ_{d} plane. The blue line denotes the scale of MOND in terms of the characteristic surface density . 

In the text 
Fig. 2.
Left panel: Q calculated for both MI and MG for a Kuzmin disk using mock data binned into acceleration bins. The gray lines denote the deep MONDian limits for MI and MG, respectively. Right panel: ΔQ^{MI} − ΔQ^{MG} as a function of acceleration scale. 

In the text 
Fig. 3.
Left panel: arithmetic mean of Q for all mock galaxies with max[g_{obs}]≤0.01a_{0} as a function of radial range in units of r_{d}. Black and blue lines represent MI and MG, respectively, with v_{∞} ≃ v_{flat} = v_{tot}(500 kpc). Cyan and brown lines represent MI and MG, respectively, with v_{∞} ≃ v_{flat} = v_{tot}(r_{cut}), where r_{cut} denotes the radii at which the integrals are cut. The inset is a magnification of N_{range} = [0, 10]. Right panel: ΔQ^{MI} − ΔQ^{MG} as a function of radial range. The error bars are smaller than the points. 

In the text 
Fig. 4.
Left panel: arithmetic mean of Q for all mock galaxies with max[g_{obs}]≤0.01a_{0} as a function of spacing between sampled points. Black and blue lines represent MI and MG, respectively. Right panel: ΔQ^{MI} − ΔQ^{MG} as a function of spacing between sampled points. 

In the text 
Fig. 5.
Q for individual galaxies as a function of max[g_{obs}] (first row), (second row) and spacing (third row). These 15 galaxies are deep enough in the MOND regime and have sufficiently extended and finely sampled rotation curves (see Sect. 2). The gray line shows the predicted Q from MOND modified gravity, while the gray band shows the ranges of values expected in MOND modified inertia. 

In the text 
Fig. A.1.
Q calculated from MOND modified inertia and MOND modified gravity in the Brada & Milgrom (1995) approximation from the mock data binned into acceleration bins. Each panel represents a different mass distribution (denoted on the panel). See Eqs. (A.1) and (A.2) for details. 

In the text 
Fig. B.1.
Rotation curves of galaxies in the SPARC database after imposing max[g_{obs}]≤0.4a_{0}, and spacing ≤1.3r_{d} as well as requiring that there be a gas profile and a welldefined approximation of v_{∞}. Black is v_{obs}, blue is v_{bar}, cyan is v_{g}, red is v_{d} and brown denotes the 1σ region for v_{f}. The brown band illustrates the points used to compute v_{f}. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.