A&A 393, 453-460 (2002)
DOI: 10.1051/0004-6361:20021021
R. Bottema 1 - J. L. G. Pestaña2 - B. Rothberg3 - R. H. Sanders1
1 - Kapteyn Astronomical Institute, PO Box 800, NL-9700 AV Groningen, The Netherlands
2 - Dpto. de Física, Univ. de Jaén, Virgen de la Cabeza 2, 23071-Jaén, Spain
3 - Institute for Astronomy, 2680 Woodlawn Dr., Honolulu, HI 96822, USA
Received 18 February 2002 / Accepted 9 July 2002
Abstract
Rotation curves for four spiral galaxies with recently
determined Cepheid-based distances are reconsidered
in terms of modified Newtonian dynamics (MOND). For two of the objects,
NGC 2403 and NGC 7331, the rotation curves predicted by MOND
are compatible with the observed curves when these galaxies are
taken to be at the Cepheid distance.
For NGC 3198, the largest distance for which reasonable
agreement is obtained is 10% smaller than the Cepheid-based distance;
i.e., MOND clearly prefers a smaller distance. This conclusion is
unaltered when new near-infrared photometry of NGC 3198 is taken as the
tracer of the stellar mass distribution. For the large Sb spiral, NGC 2841,
MOND requires a distance which is at least 20% larger than the Cepheid-based
distance. However, the discrepancy of the Tully-Fisher and
SNIa distances with the Cepheid determination casts some doubt upon the
Cepheid method in this case.
Key words: galaxies: individual: NGC 3198, NGC 2841 - galaxies: kinematics and dynamics - galaxies: spiral
It is well established that, in the context of Newtonian dynamics, the observable mass in spiral galaxies cannot account for the observed flat rotation curves in the outer regions of galaxies (Bosma 1978; Begeman 1987; van Albada et al. 1985). The standard explanation for this discrepancy is the proposal that galaxies are embedded in an extended dark halo which dominates the gravitational field in the outer regions (Trimble 1987).
An alternative explanation for the discrepancy
is the possibility that dynamics becomes non-Newtonian
in the limit of low accelerations.
The most successful such proposal is Milgrom's (1983)
modified Newtonian dynamics or MOND.
Here the idea is that below a certain acceleration
threshold (
)
the effective gravitational acceleration
approaches
where
is the usual
Newtonian acceleration. This modification yields
asymptotically flat rotation curves of spiral
galaxies and a luminosity - rotation velocity
relationship of the observed form,
,
the Tully-Fisher relation (Tully & Fisher 1977). But apart from
these general aspects the prescription also successfully
predicts the observed form of galaxy rotation curves
from the observed distribution of stars and gas
with reasonable values for the mass-to-light ratio
of the stellar component (Begeman et al. 1991;
Sanders 1996; Sanders & Verheijen 1998; McGaugh & de Blok 1998).
A crucial element of a very
specific prescription like MOND is that the
precise form of the rotation curve is predicted
by the observed mass distribution given the
value of a single universal parameter; in this
case, the critical acceleration a0. Consequently
MOND can in principle be falsified as soon
as there is one galaxy for which the predicted rotation curve
disagrees significantly with the observed curve; although,
in practice, the usual uncertainties inherent in astronomical data
render a definitive falsification problematic in any
individual case.
In Begeman et al. (1991, hereafter
BBS) MOND is applied to
a sample of galaxies for which high
quality H I rotation curves are available.
For a value of a0 equal to 1.21
10-8 cm s-2the rotation curves of the sample could
be reasonably reproduced, the free parameter in each case being the
mass-to-light ratio of the visible disc. Because MOND is an acceleration
dependent modification, this derived value of
depends upon assumed distance scale (
km s-1 Mpc-1 in
this case). Moreover, the quality of an individual fit depends
upon the adopted distance to the galaxy, and, since the relative distances
to these nearby galaxies have not been known to within
an accuracy, typically, of 25%, this has provided
some freedom to adjust the distance in order to improve the MOND fit; i.e.,
distance, within certain limits, can be considered as an
additional second parameter in the fitting procedure. For most of
the galaxies in the sample of BBS, the distance did not have to
be adjusted significantly (<
)
to improve the MOND fits, and the
improvement was not significant. However, one object, NGC 2841, required
a large readjustment: the Hubble law distance to this galaxy is
about 9 Mpc, but MOND clearly prefers a distance which is twice as
large.
Using ground-based and Hubble Space Telescope observations Cepheid distances to 21 inclined galaxies have now been determined as part of the HST key program on the extragalactic distance scale (e.g. Sakai et al. 2000). Three of the galaxies in this Cepheid sample are also in the sample with high quality rotation curves considered by BBS. These are NGC 2403 (Freedman & Madore 1988), NGC 3198 (Kelson et al. 1999) and NGC 7331 (Hughes et al. 1998). For these three galaxies the MOND prescription can now be considered in the context of the Cepheid distance that is generally considered to be the most precise indicator.
NGC 2841 has been discussed as a critical case for MOND by Sanders (1996). For this galaxy, there is also a large discrepancy between the Tully-Fisher distance and the Hubble law distance (for plausible values of the Hubble constant). Moreover, the galaxy was the site of a recent SNIa (1999by). For these reasons this galaxy has been included, subsequently, in the HST program (Macri et al. 2001).
Here we demonstrate
that for two galaxies in the BBS sample the rotation curve predicted by
MOND is consistent with the observed curve when the galaxies
are placed at the Cepheid distance. However,
for NGC 3198 at the Cepheid distance of 13.8
0.6 Mpc,
the shape of the rotation curve predicted by MOND
systematically deviates (by up to 10 km s-1) from the observed curve,
both in the inner and outer regions. The largest distance which can
be compatible with MOND is about 10% lower than the Cepheid-based distance.
This is not particularly problematic because of likely uncertainties
in the Cepheid method and in the determination of a rotation curve
from the observed two-dimensional velocity field.
NGC 2841, however, remains a difficult case for MOND. The minimum
distance which is consistent with MOND is about 17 Mpc whereas
the Cepheid-based distance is 14.1
1.5 Mpc. We discuss the
implications and seriousness of this discrepancy for MOND, or, alternatively,
for the Cepheid method.
The procedure followed when determining a MOND rotation
curve has been described previously (e.g., BBS).
In the context of MOND, the true
gravitational acceleration g, is related to the Newtonian
acceleration
as
| (1) |
| (2) |
![]() |
(3) |
![]() |
(4) |
| v4 = GMa0, | (5) |
The stellar disc may be assumed to be asymptotically thin or have a finite thickness related to the radial scale length of the disc by an empirical rule (van der Kruit & Searle 1981); this makes little difference in the final result. Applying Eqs. (1), (3), and (4), a least squares fit is then made to the observed rotation curve v(r) where the single free parameter of the fit is the mass-to-light ratio of the disc; in cases where there is an indication of a bulge from the light distribution, M/L of the bulge enters as a second parameter.
For the gaseous component a surface density distribution equal to that of the H I is taken, multiplied by a factor 1.3 to account for primordial helium. The gas layer is taken to be infinitesimally thin. The contribution of the gas to the total rotation is fixed, but does depend on the distance to the galaxy.
In principle, the parameter
should be universal and, having
determined its magnitude, one is not allowed to adopt this as a free
parameter. But as noted above, the derived value of
does
depend upon the assumed distance scale. Sanders & Verheijen (1998)
give MOND fits to the rotation curves of 30 spiral galaxies in the
UMa cluster which they assume to be at 15.5 Mpc. The preferred
value of
with this adopted distance is equal to the BBS value
of
cm s-2. However, based upon the Cepheid-based
re-calibrated Tully-Fisher relation (Sakai et al. 2000), Tully & Pierce
(2000) argue that the
distance to UMa should be taken to be 18.6 Mpc.
We have recalibrated the Tully-Fisher law using this same sample of galaxies
but with the three test galaxies (NGC 2403, NGC 3198, NGC 7331) left out
of the fitting. Within the errors, the slope and intercept of the
Tully-Fisher relation are the same as that found by Sakai et al. (2000),
and the distance to Ursa Major is only 1% smaller than that found
Tully & Pierce (2000). In that case the MOND fits to the UMa galaxies imply
that the value of
should be adjusted to
cm s-2. This is also the
preferred value of
from MOND fits to rotation curves of a sample of
nearby dwarf galaxies with distances taken primarily from group membership
(Swaters & Sanders 2002, in preparation).
In Fig. 1 we show the MOND rotation curve
for all four galaxies from the BBS sample with Cepheid-based distance
determinations.
NGC 2841 and NGC 7331 both contain central bulges as evidenced
in the light distribution, and the radial surface brightness
profile has been appropriately decomposed. Here,
is fixed at the
rescaled value of
cm s-2, and the distance is fixed
at the Cepheid-based values as updated and corrected by Freedman et al.
(2001). The free parameters
of the fit are the disc and, in two cases, bulge masses.
The resulting values and
the corresponding mass-to-light ratios are given in Table 1 for the
four galaxies.
![]() |
Figure 1: MOND rotation curves compared to observed H I rotation curves for the four galaxies from the sample of BBS with Cepheid-based distances. The dotted, long-dashed, and short-dashed lines are the Newtonian rotation curves of the stellar disc, bulge, and gaseous components respectively. |
| Open with DEXTER | |
|
|
D | LB |
|
|
|
|
|
|
|
|
|
|
|
|||
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) |
| NGC 2403 | 3.2 |
0.82 | 0.4 | 1.34 |
1.6 | ||
| NGC 2841 | 14.1 |
4.60 | 2.7 | 29.70 |
8.3 | 1.5 | 0.83 |
| NGC 3198 | 13.8 |
2.44 | 1.6 | 2.63 |
1.1 | ||
| NGC 7331 | 14.7 |
5.26 | 1.4 | 13.20 |
2.0 | 5.7 | 1.8 |
|
(2) The Cepheid-based distance from Freedman et al. (2001). (3) The B-band luminosity (in (4) The total gas mass including primordial helium at the Cepheid distance. (5) The total mass of the stellar disc from the MOND fit. (6) The implied mass-to-light ratio of the stellar disc. (7) The total mass of the stellar bulge in those two cases where a bulge is evident. (8) The implied mass-to-light ratio of the stellar bulge. |
For the other two galaxies, there are clear systematic differences between the MOND rotation curve and the observed curves. Basically, the predicted curves have a different shape than the observed curves: for NGC 2841, the predicted curve is significantly higher than observed in the inner regions (by up to 30 km s-1) and comparably lower in the outer regions. For NGC 3198 the differences are in the opposite sense: about 10 km s-1 lower in the inner regions and 10 km s-1 higher in the outer regions. These differences diminish if NGC 2841 is moved further out and if NGC 3198 is moved closer in; i.e., MOND clearly prefers a larger distance to NGC 2841 (as discussed previously by BBS and by Sanders 1996) and a smaller distance to NGC 3198. We now discuss these two cases with respect to the question of whether or not this mismatch can be interpreted as a falsification of MOND. Because the rotation curve of NGC 3198, when taken at the Hubble law distance of 10 Mpc, is very well predicted by MOND, and because the observed curve is thought to be well-determined, this, at first sight, appears to be the more problematic case, and we begin with this object.
This gas-rich spiral galaxy has a generally symmetric H I distribution, and there are no large scale significant warps or distortions of the velocity field. The rotation curve extends to roughly 10 radial scale lengths and is, to first order, flat and featureless (Begeman 1987). For these reasons it has become the classic case of a spiral galaxy evidencing a large mass discrepancy in its outer regions (van Albada et al. 1985). If any theory, such as MOND, fails to predict the rotation curve of this galaxy, then it would be problematic for that theory.
![]() |
Figure 2: MOND rotation curves for NGC 3198 assumed to be at various distances. The Cepheid-based distance is 13.8 Mpc. The dotted and dashed lines are the Newtonian rotation curves of the stellar and gas discs respectively. |
| Open with DEXTER | |
In Fig. 2 we show the MOND rotation curves of
NGC 3198 when the galaxy is assumed to be at distances of
10 Mpc, 12.5 Mpc and 13.8 Mpc. Again, the MOND acceleration
parameter is assumed to be the BBS value rescaled to the new distance
scale, i.e.,
cm s-2.
The closest assumed distance, 10 Mpc, is roughly the Hubble law distance given the radial velocity of NGC 3198 with respect to the local group; it is also the least-square-fit distance if distance is left as a free parameter in the context of MOND. The distance of 13.8 Mpc is the final Cepheid-based distance given by Freedman et al. (2001); and 12.5 Mpc corresponds to the Cepheid distance less 10%. The disc M/L values in the B-band corresponding to the MOND fits at these various distances are also given in Fig. 2.
Here we see that the MOND rotation curve for a distance of 10 Mpc is essentially a perfect fit to the observed curve. At the distance of 12.5 Mpc, the MOND curve is less than a perfect match, but, nowhere that the rotation curve is well measured, does the predicted rotation curve deviate by more than 5 km s-1 from the curve derived from the observed velocity field. This is typically within the difference in the rotation curves derived from the two sides of the galaxies considered separately- a sensible estimate of the uncertainties (the error bars are formal errors determined from the tilted ring fitting procedure).
At the Cepheid distance of 13.8 Mpc, the MOND rotation curve deviates in the same sense but by now up to -10 km s-1 in the inner regions (8-14 kpc) and by +10 km s-1 in the outer regions (30-40 kpc). The reason for the deteriorating fit with increasing assumed distance is the relatively larger contribution of the gaseous component to the rotational velocity. The rotation curve of NGC 3198 in the outer regions (r>20 kpc) is constant at about 150 km s-1. This would imply, in the context of MOND, that essentially the entire mass of the galaxy is enclosed within about 20 kpc, but this is obviously not the case given the significant surface density of neutral gas in the outer regions - contributing more than 50 km s-1 to the Newtonian rotation curve at the last measured point.
At a distance of 12.5 Mpc, the MOND rotation curve appears to be consistent
with the observed curve (within the likely errors of the method for estimating
rotation curves from 21 cm line data). Although this distance is
formally 2
below the Cepheid-based distance, it is unclear if all systematic
effects connected with this method are well-understood.
It has been noted, for example, that for the galaxy NGC 4258
the kinematic water-maser-based distance is also about 10% less than
the Cepheid-based distance (Maoz et al. 1999).
The error budget of the
Cepheid method is probably on the order of 10%.
![]() |
Figure 3: The I-band Tully-Fisher relation for local calibrator galaxies with Cepheid-based distances. The positions of NGC 3198 (diamonds) and NGC 2841 (crosses), assumed to be at various distances, are indicated. |
| Open with DEXTER | |
Sakai et al. (1999) have calibrated the T-F relation using 21 spiral galaxies with known Cepheid distances in five color bands: B, V, R, I, and H. If one places NGC 3198 on the mean B-band relation its distance should be 12.2 Mpc, while for the I-band this distance is 13.3 Mpc. Thus the Tully-Fisher distance is essentially consistent with the maximum MOND distance. Although the MOND rotation curve fit clearly prefers a somewhat smaller distance than the Cepheid-based distance, the idea is in no sense falsified by this well-determined rotation curve.
The I-band Tully-Fisher relation from Sakai et al. is shown in Fig. 3. The open points show the position of NGC 3198 when at a distance of 10.0 Mpc, 12.5 Mpc, and 13.8 Mpc. It is evident that, given the scatter in the observed relation, it is impossible to distinguish between these possibilities although distances of 12.5 to 13.8 Mpc are clearly preferred.
One possible reason for the small deviation of the MOND curve from the observed curve at the Cepheid-based distance is that the r-band photometry is not a precise tracer of the stellar light distribution due to possible contamination by newly-formed stars and dust absorption. For this reason we have also considered recent near-infrared photometry of this galaxy.
![]() |
Figure 4:
The
|
| Open with DEXTER | |
An image of NGC 3198 in the
band has been obtained by
Rothberg et al. (2000) in order to calibrate the near infrared
Tully-Fisher relation.
The observations and initial stages of the data reduction, like
sky-subtraction and flat-fielding are described in that paper.
The detector was 1024
1024 square pixels
of size 1
68
1
68. Consequently the total
image measures 28.7 arcmin along the sides and NGC 3198 which
has a scale-length of approximately 1 arcmin fits completely within
the image leaving ample margins of pure sky around the galaxy.
Rothberg et al. (2000) derived a total brightness of 7.79
magnitudes
which translates to 3.4
1010
for a distance
of 13.8 Mpc.
In Fig. 4 the image of the central regions of NGC 3198 is reproduced. Clearly discernible is a prominent bulge which is much less obvious in images at bluer wavelengths. As a consequence this central bulge region must be enshrouded in an appreciable amount of dust, which explains the reddening going inward. Surrounding the bulge appears to be a ring of spiral arm features with a light depression between the bulge and this ring.
To determine the radial luminosity profile, ellipses have been fitted
to the image which provided the position and orientation of the galaxy
(Fig. 5). As a next step the intensities have been averaged over
elliptic annuli.
In the inner regions the orientations of the annuli were equal to
those determined by the ellipse fit, while for the intermediate and
outer regions a constant position angle and ellipticity was adopted.
The error of each radial intensity value was calculated by quadratically
adding the error generated by sky-level variations and the noise
appropriate for each annulus. The result is shown in Fig. 6.
The radial profile
in the r-band is also plotted in that figure, and one may notice that
the photometry of the disc is of similar shape for the r and
bands.
It is without doubt that NGC 3198 has a bulge or central light concentration. A possible bulge/disc light decomposition is shown in Fig. 7. Here, it is assumed that the mass surface density is exactly proportional to the observed intensity level, and that the light and mass distribution are axisymmetric. For that case the bulge/disc decomposition illustrated in Fig. 7 is essentially a decomposition by eye. Here, it is further assumed that the stellar disc has a central hole with a radius corresponding to that of the light depression and that the bulge extends slightly beyond this radius.
The light depression might well be caused by the presence
of a central bar. The influence of a bar on the radial velocity field
of the gas is suggested in high-resolution H
images of the
galaxy, where characteristic distortions from circular motion are
evident (Corradi et al. 1991).
The light depression
would then be due to a real deficiency of matter near the L4 and L5
Lagrangian points along the minor axis of the bar (Bosma 1978). In that
case, the bar would be oriented nearly parallel to the line-of-sight
and would not be photometrically conspicuous. Moreover, the bar would
affect the derived rotation curve in the inner regions, or, at least,
the interpretation of the rotation curve as a tracer of the radial force
distribution. A bar aligned with the minor axis of the galaxy image would
have the effect of increasing the apparent rotation velocities in the
inner region (Teuben & Sanders 1985); however, this would be significant
only within the inner 30 arcsec (
2 kpc)
and would have little influence upon
the overall shape of the derived 21 cm line rotation curve.
Keeping this caveat in mind, we proceed using the decomposition
depicted in Fig. 7: assuming a spherical bulge and disc with observed ellipticity, the total
luminosity of 3.4
1010
,
is
divided into 8.23
109 and 29.09
109
for the bulge and disc respectively.
Because the scale-length of the disc in
is nearly equal
to the disc scale-length in the optical, it is not to be expected
that the MOND fit will be much different from that for the
r-band photometry. This is the case, as can be seen in Fig. 8
where the MOND rotation
curve again has been determined at the Cepheid-based distance of
13.8 Mpc. Here, except for a spike in the central regions which is
due to the bulge, the predicted rotation curve is essentially
the same as derived from the r-band photometry; that is to say, the
conclusions are unchanged by the near-infrared results.
In Fig. 9 we show the MOND rotation curve for NGC 2841 compared to
the observed curve at various assumed distances: 15.6 Mpc, which
is the 1
upper limit on the Cepheid-based distance,
17 Mpc which is 20% larger than the Cepheid-based distance,
and 23 Mpc which is the MOND-preferred distance. As in Fig. 2 the
M/L values for the disc and bulge are also given in the figure.
![]() |
Figure 5:
The orientation and inclination of ellipses fitted to the
|
| Open with DEXTER | |
![]() |
Figure 6:
Radial light profile in the
|
| Open with DEXTER | |
![]() |
Figure 7: Radial light profile in the central regions and a possible decomposition into a bulge (dotted line) and a disc (dashed line). |
| Open with DEXTER | |
![]() |
Figure 8:
The MOND rotation curve of NGC 3198 where the surface density distribution
of the stellar component is taken to be traced by the
|
| Open with DEXTER | |
![]() |
Figure 9:
MOND rotation curves for NGC 2841 at various distances ranging from
the one-sigma upper limit on the Cepheid-based distance (15.6 Mpc) to
the Tully-Fisher and SNIa distance (23 Mpc). The
|
| Open with DEXTER | |
The rotation curve, as a tracer of the radial force distribution in
this galaxy, is actually not as well-determined
as that of NGC 3198. There is a significant warp in the outer regions
which must be modelled by the tilted-ring technique, and this adds
uncertainty to the derived rotation curve (see comments by Bosma 2002).
None-the-less, it is
clear that, while the Cepheid distance goes in the right direction
(it is significantly larger than the Hubble based distance), it is not
enough to bring the MOND-predicted rotation curve into agreement with
the observed curve. Moreover, not only does the form of the predicted
curve differ systematically from that observed, but it is clear that
the M/L value for the disc is un-naturally large (6.8
) -
larger than that required for the bulge (2.2).
Both of these problems are relieved somewhat if the distance is taken to be
20% larger than the Cepheid-based determination - at 17 Mpc. There are
still systematic deviations in the form of the rotation curve, but these
become large (in the outer regions), only where the gas layer of the
galaxy is observed to be significantly warped. We may take this as
a lower limit on the distance which would be compatible with MOND, although
the disc M/L does remain uncomfortably large (5.9
)
and
that of the bulge rather small (2.7).
Leaving distance as a free parameter in the fit yields a MOND-preferred distance of 23 Mpc, and here we see that the rotation curve fit is perfect with very reasonable implied M/L values for the disc and bulge. This is entirely consistent with the distance implied by the Cepheid-calibrated Tully-Fisher relation as is also shown in Fig. 3. Taking the galaxy to be at the Cepheid distance of 14.1 Mpc, we see that the galaxy lies about one entire magnitude below the mean line of the TF relation. The distance implied by I-band TF relation is 24 Mpc.
There has been a recent supernova in NGC 2841 (SN 1999by), which is
type Ia, i.e., the fundamental extragalactic "standard candle''.
However, if the galaxy is at the Cepheid distance of 14.1 Mpc,
SN 1999by is one of the least luminous supernovae Ia ever observed,
with a peak absolute magnitude of MB = -17.15
0.23. Based upon
an estimate of the decline-rate parameter (
)
Garnivich et al. (2001) argue that this supernova is a peculiar low
luminosity event, and they use this event and several others to
recalibrate the Phillips relation (Phillips et al. 1999) between decline
rate and peak luminosity. However, if we take the Phillips relation
at face value then the peak luminosity of this object would
be MB = -18.3, which would imply that the distance to the galaxy
would be 23.5 Mpc. It is interesting that an earlier SN event in
NGC 2841, SN 1957A, would be, if the galaxy is at the Cepheid distance,
the faintest supernova type Ia ever observed (MB = -16.4). It is
curious that this galaxy only seems to provide sub-luminous
supernovae.
The deviation of the galaxy from the TF relation and the abnormally low peak powers of supernovae, suggest that the Cepheid distance to this object may be substantially too low. It has been argued that the Cepheid method may be adversely affected by blending: the true apparent brightness of Cepheids is enhanced by blending with the light of nearby stars. This would lead to an underestimate of distances based upon the period-luminosity relation, and would affect, in particular, the more distant objects (see Paczynski & Pindor 2001 for a discussion of these points). All we can conclude, at the moment, is that the MOND-preferred distance to NGC 2841 remains significantly larger than the present Cepheid-based distance.
The main conclusions of this paper can be summarized as follows:
It is clear that NGC 2841 remains a critical case for MOND. The discrepancy between Cepheid-based distance and both the TF and SNIa based distances to NGC 2841 suggests that there may be a problem with the derived Cepheid-based distance.
In general, it is evident that accurate distance determinations to nearby galaxies are extremely relevant to the question of the viability of MOND. MOND, as a modification of Newtonian dynamics attached to an acceleration scale, is far more fragile than the dark matter hypothesis in this regard. It would be useful to obtain more Cepheid-based distances to the sample of galaxies with well-observed rotation curves. Particularly useful would be a Cepheid distance estimate to the Ursa Major cluster as many of these galaxies have well-measured rotation curves and near-infrared photometry.
Acknowledgements
We are grateful to M. Milgrom for a critical reading of the manuscript and we thank the referee A. Bosma, for giving useful comments which have improved the paper.