text:
# High-resolution rotation curves of Low Surface Brightness galaxiesbased on observations at the Observatoire de Haute Provence

## 1 Introduction

Despite much effort, it is still unclear to what extent rotation curves can give clues about the distribution of the visible and dark matter in bright spiral galaxies (e.g. Bosma 1999; Sellwood 1999). It is thought, however, that low surface brightness (LSB) galaxies and dwarf galaxies are dark matter dominated, and that therefore the analysis of their rotation curves can yield directly important information about the properties and distribution of their associated dark matter halos (de Blok & McGaugh 1997; Verheijen 1997; Swaters 1999).

This implication has far-reaching consequences, as early results by e.g. Moore (1994) already showed. Predictions from cold dark matter (CDM) simulations were found to disagree with observations of rotation curves of several dwarf galaxies; the data indicating much less cuspy distributions of matter than the simulations. It was thought at the time that this problem could be solved once the effect of feedback due to star formation was understood. However, the behaviour of the rotation curves of low surface brightness (LSB) galaxies is rather similar, and their low star formation rate at present and in the past indicate that star formation in such galaxies might never have been important enough to modify their structure drastically.

Early work on LSB galaxies used 21-cm Hi line work to determine the rotation curves (e.g. van der Hulst et al. 1993; de Blok et al. 1996), and as such the results are likely to suffer from modest angular resolution effects (commonly called beam smearing). Even though these can be partly modelled, experience shows that direct measurements are preferable (cf. discussions in Bosma 1978; Begeman 1987); in particular, supplementary data in the optical emission lines, such as H$`\alpha `$ and \[Nii\] are always useful (cf. Van der Kruit & Bosma 1978; Rubin et al. 1989; Corradi et al. 1991; Swaters et al. 2000). Thus a whole industry has sprung up to combine optical and 21-cm line rotation curves of all sorts of gas rich galaxies.

For the particular problem of the dark matter distribution in LSB galaxies Swaters et al. (2000) presented supplementary H$`\alpha `$ data for five LSB galaxies already observed in Hi by de Blok et al. (1996), and concluded that the influence of beam smearing on the Hi curves was severe enough to question earlier conclusions regarding dark matter content and rotation curve shape of LSB galaxies. McGaugh et al. (2001) and de Blok et al. (2001a) (hereafter dBMR) reanalysed these data, and concluded that the discrepancy between H$`\alpha `$ and Hi data is only really significant for one of these five galaxies. Conclusions regarding the dark matter content and the shape of the dark matter distribution in LSB galaxies are thus not affected.

In a comparison of pseudo-isothermal and CDM halo models using high-resolution rotation curves of a sample of a further 29 LSB galaxies dBMR show that the so-called “universal” CDM halo-profile as parameterised in Navarro et al. (1996) is not a good description of the data: the rotation curves generally show linear solid-body rise in the inner parts, which is inconsistent with the steeper rise necessary for the CDM rotation curves. Rather, the rotation curves prefer a pseudo-isothermal (i.e. core-dominated) halo model. de Blok et al. (2001b) furthermore showed that *all* rotation curves of LSB galaxies measured so far are consistent with a pseudo-isothermal core model. Analyses that show that these rotation curves confirm the CDM NFW halo model (e.g. van den Bosch & Swaters 2001) can be traced back to the fact that at lower resolutions the NFW and pseudo-isothermal models look sufficiently similar and the errors are large enough that a NFW model can be made to fit the data. The higher-resolution data presented in de Blok et al. (2001b) now seem to have settled the observational side of the debate in favour of core-dominated LSB galaxy halos. These results are independent of the stellar mass-to-light ratios one assumes in constructing these mass models.

It would therefore seem that current models of structure formation and galaxy evolution need to take into account the fact that most late-type galaxies have a constant density dark matter *core* rather than a cusp (see also results by Borriello & Salucci 2001; Salucci & Burkert 2000; Salucci 2001 for HSB galaxies).

As the signature of the core is clearest at small radii it is important to find galaxies where these inner parts are well-resolved and well-sampled. In this paper we thus supplement the collection of high-resolution rotation curves of LSB galaxies from de Blok et al. (2001a) with curves for an additional 26 galaxies, of which 12 are entirely new and 14 were already used in the analysis by de Blok et al. (2001b). The 12 new galaxies in our sample have been specifically chosen to have small distances so that we can easily verify the discrepancies between the NFW model and the observations.

This paper is organised as follows: we describe our data in Section 2, and present the rotation curves in Section 3. Section 4 describes individual galaxies. In Section 5 we describe the derivation of the final rotation curves. The mass models are then presented in Section 6. In Section 7 we discuss these models. Section 8 digresses into the consequences that systematic observational effects may or may not have on the data. In Section 9 we discuss the mass densities profiles inferred from the rotation curves. Section 10 summarises the paper.

## 2 Observations and Data reduction

### 2.1 Sample selection

We constructed a sample of galaxies to observe in H$`\alpha `$ from lists of LSB galaxies already observed in the 21-cm Hi line by van der Hulst et al. (1993) and de Blok et al. (1996). A representative sample of nearby dwarf galaxies from the theses of Swaters (1999) and Stil (1999) was also included (see also Swaters et al. 2002a, b). The sample was supplemented by two well-known nearby dwarf galaxies (NGC 1560 and NGC 100) and a galaxy (UGC 711) from the Flat Galaxy Catalog (Karachentsev et al. 1993). Hi velocity fields were available for all galaxies (except U711 and N100). Since we observed with a long-slit spectrograph, we determined for each galaxy the position angle of the major axis of the velocity field, using the position angle of the major axis of the light distribution as an additional check. For U711 and N100 we only used the position angle of the light distribution: since these are edge-on galaxies these angles are well determined.

Table 1 lists parameters of these galaxies for which we managed to obtain useful, high S/N data that will be used in our mass modelling. The last column in Table 1 gives references to the data presented in Table 1 not derived in this paper. Images and surface brightness profiles can also be found in these references. See also Sect. 2.4. Table 2 gives a brief summary of those galaxies that we did observe, but which for various reasons were deemed to be of insufficient quality to warrant mass modelling. Our main sample therefore consists of the galaxies listed in Table 1.

### 2.2 Optical spectra

The emission line observations have been made with the 193-cm telescope at the Observatoire de Haute Provence and its long-slit Carelec spectrograph from 2-9 Jan 2000 and 22-27 Feb 2001. The spectrograph has been described in detail by Lemaître et al. (1991). As a detector, we used a EEV 2048 x 1024 CCD chip. The slit length on the sky was about 5.5′, and we used a 2.0″ slit-width. This works out to a pixel size of 0.6″ $`\times `$ 45.7 km s-1. Spectral resolution was 54 km s-1 FWHM. Typical exposure times were 1 hour per spectrum, preceded and followed by calibration spectra from a Neon arc lamp. The typical seeing during the observations was $`2\mathrm{}`$. Care was taken to align the slit with the optical centre of the galaxy. This was done with off-set pointing from nearby stars. In many cases the central part of the galaxy was visible in the guiding camera, and the accuracy of the off-set procedure could be verified. For a few of these galaxies we repeated the off-set procedure a few times, and found that the galaxy acquisition was indeed repeatable, with an error of less than 1′′. See Sect. 8 for a full discussion of the effects of (mis-)pointing.

### 2.3 Data reduction

The data were corrected for instrumental effects (bias, flat-field, etc.) in the standard manner. The spectra were then further reduced with the Figaro package: the calibration spectra were averaged, and for each of the spectra a wavelength solution was found and the wavelength calibration applied. Cosmic rays were corrected with an automated method, which was verified interactively.

The data were then further processed using the emsao programme of the rvsao package (Kurtz & Mink 1998) within the iraf environment. Emission lines were traced and their wavelengths determined using known night sky lines. Spectral shifts were then converted into barycentric radial velocities. The spectra were averaged every 3 pixels in spatial direction (1.8″). Some spectra with low S/N were averaged every 9 pixels (5.4″) as indicated in Table 1.

### 2.4 Data collected from the literature

Since our aim is to make mass models of our galaxies, we collected supplementary data from the literature, basically optical surface photometry and Hi data. The references are listed in the notes to Table 1 (see also Swaters et al. 2002a, b). For 22 out of the 26 galaxies studied here optical $`R`$-band ($`B`$-band where noted in Table 1) is available, as well as Hi imaging and Hi and H$`\alpha `$ rotation curves. Of these 22 galaxies, 5 are (close to) edge-on and because of the uncertainties involved in de-projecting surface brightness profiles of edge-on galaxies we will not consider their surface photometry here. See also Sect. 6.1.

For U628 we only have $`R`$-band photometry and H$`\alpha `$ velocity data. For two galaxies (U711 and N100) we only have H$`\alpha `$ data. For U10310 we only have Hi and H$`\alpha `$ data. For these four galaxies we can construct only limited mass models.

Generally we adopted the distances given in the source papers (corrected to $`H_0=75`$ km s-1 Mpc-1), except for DDO 47 where we used the revised distance given in Walter & Brinks (2001), and for U711 where we used the recession velocity. Inclinations were also adopted from the source papers. These inclinations are for the major part kinematically derived and as such only have a few degrees error at most. It should be emphasised that the shape of the rotation curve and hence applicability of any particular model does not depend on inclination, and thus an uncertainty of a few degrees does not affect our main conclusions. The scale lengths and surface brightnesses were are also taken from the references listed in Table 1. We refer to these papers for precise details. Here it suffices to say that they were derived using a simple straight-line fit to the exponential part of the profile (most of the sample galaxies are dominated by an exponential disk).

## 3 Results

### 3.1 Raw HI and H$`\alpha `$ data

Our principal results are collected in Figs. 1 to 5. Fig. 1 shows the raw velocity curves derived from the H$`\alpha `$ data for the galaxies listed in Table 1. Fig. 2 shows the raw data of the galaxies listed in Table 2. These latter data will not be analysed further.

Fig. 3 presents overlays of the H$`\alpha `$ data on a major axis position-velocity diagram from the Hi data. Such a presentation has the advantage of immediately showing any discrepancies between the two datasets, if they exist. The Hi data are taken from the literature, with references given in Table 1. Also indicated is the Hi rotation curve as derived from the position-velocity diagram. These overlays will be discussed in more detail in Sect. 4.

### 3.2 Symmetrised Rotation curves

In most cases the continuum emission from the central parts of the galaxies was visible in the spectra. This was used to determine the central position in the spectra and the systemic velocities. In a few cases (DDO185, U10310, U5999) the centre was not seen and we varied the central position and systemic velocity over a small range until we found the values that gave the maximum symmetry. The systemic velocities based on the H$`\alpha `$ data are given in Table 1. A comparison with the systemic velocities derived from the HI data (also listed in Table 1 and taken from the source papers) shows a good correspondence with the H$`\alpha `$ velocities. The absolute difference is 5 km s-1 or less, except in 6 cases where the difference lies between 10 and 14 km s-1. In 4 of the 6 cases the HI data were obtained using a velocity resolution of over 20 km s-1, implying that the comparison with H$`\alpha `$ is still very favourable. Only for DDO 47 and DDO 52 do we find a difference of $`1.5`$ HI channel widths. Given the limited spatial and velocity resolution of the HI data it is not clear where this difference originates.

Figure 4 show the folded and symmetrised rotation curves, corrected for inclination.

Figure 5 compares the folded H$`\alpha `$ rotation curves with the Hi rotation curves as presented in the literature (see Table 1 for references). It is clear that *grosso modo* there is reasonable agreement. A small number of galaxies were obviously affected by beam-smearing, just as for some galaxies the H$`\alpha `$ rotation curve shows signs of non-circular motion.

## 4 Description of individual galaxies

In this section we describe the various rotation curves of each galaxy in some detail. We will compare the Hi and the H$`\alpha `$ data, as well as the Hi and H$`\alpha `$ rotation curves as presented in Figs. 1 and 3-5. For some galaxies we make a comparison with independent data from the literature.

F563-1: The H$`\alpha `$ data show reasonable agreement with the Hi data, showing that beamsmearing is not a major problem. At the receding side the H$`\alpha `$ data disagree slightly with the Hi curve. An independent observation by dBMR of this galaxy is available and agrees with our observation (see Fig. 6); any systematic differences between the data sets is less than the uncertainties in individual data points. Obtaining high-resolution H$`\alpha `$ rotation curves of LSB galaxies is thus a repeatable exercise.

U5005: The H$`\alpha `$ data agree with the Hi position-velocity diagram, but the resulting high-resolution curve rises less steeply than the Hi curve. The Hi data have only a low spatial resolution of $`40\mathrm{}`$ (van der Hulst et al. 1993) and uncertain corrections for beam smearing are thus unavoidable.

U1230: The H$`\alpha `$ data are of high-quality, and shows that beamsmearing has affected the $`40^{\prime \prime }`$-resolution Hi curve (van der Hulst et al. 1993). The high-resolution curve rises more steeply, though still with a solid-body signature.

U5750: The Hi data are poor, but in agreement with the H$`\alpha `$ data. The linearly rising Hi curve also shows reasonable agreement with the H$`\alpha `$ data. For this galaxy an independent observation of dBMR is available and shows excellent agreement (Fig. 6), again showing that systematic effects due to telescope pointing etc. are negligible.

U731/DDO9: the H$`\alpha `$ is weak at the receding side and superimposed on a night sky line. This makes it virtually impossible to say anything about the H$`\alpha `$ velocities between $`+0.5^{}`$ and $`+2^{}`$. The H$`\alpha `$ is stronger at the approaching side, and agrees well with the ridge in the Hi $`pV`$ diagram. The H$`\alpha `$ curve rises more steeply than the Hi curve, which turns over too quickly. U731 is a “kinematically lopsided” galaxy (Swaters et al. 1999), and the Hi rotation curve of the approaching side is different from that of the receding side. The H$`\alpha `$ data agrees well with the steep approaching side of the Hi position-velocity diagram. The H$`\alpha `$ *curve* is however significantly steeper than the Hi *curve*. This discrepancy cannot be explained by lopsidedness effects, which are much smaller than the difference shown here. The Hi curve extends over $`8`$ beams in radius, and this galaxy thus demonstrates that even for well-resolved galaxies corrections for beam-smearing are ambiguous, and do not always yield the correct outcome.

U1281, U3137, U3371/DDO39, U4325/N2552: For these four galaxies there is excellent agreement between the Hi position velocity data and the H$`\alpha `$ curve. There is some disagreement between the Hi and H$`\alpha `$ rotation curves for U4325 and U3371. In these cases the Hi curves rise more steeply than the H$`\alpha `$ data. The agreement between the *raw* Hi and H$`\alpha `$ data suggests that the cause of this disagreement is likely an overcorrection for beamsmearing in the Hi data. While suited for correcting for large-scale resolution effects, beam smearing corrections thus give non-unique solutions on small scales — a rotation curve that looks like it has linear rise, can in reality be a steeply rising rotation curve (as the beam-smearing corrections in this case usually assume), *or* it can be an intrinsically linearly rising curve, as the H$`\alpha `$ data in this case show.

U4173: the H$`\alpha `$ data are rather poor in quality, but do seem to agree with the Hi data.

N100/U231: good agreement between the data sets. Given the edge-on orientation of the galaxy, it reaffirms the transparency of this small galaxy, cf. Bosma et al. (1992).

N1560/U3060: excellent agreement between Hi and H$`\alpha `$. This curve is a combination of three spectra covering the entire extent of the optical disk of N1560. The H$`\alpha `$ curve shows clear signs of non-circular motions in the disk of N1560 (high-velocity gas near star forming regions, etc.).

N2366/U3851/DDO42: The H$`\alpha `$ rotation-curve is dominated at the approaching side by the effects of non-circular motions. This position corresponds to a large star forming region (Hunter et al. 2001), that is disturbing the dynamics on that side of the galaxy. The unaffected raw H$`\alpha `$ data show excellent agreement with the Hi position-velocity diagram, the H$`\alpha `$ rotation curve, however, rises steeper than the Hi curve.

N3274/U5721, N4395/U7524: In general there is excellent agreement between the various raw data sets. The H$`\alpha `$ curves rise a bit faster than the Hi curves.

N4455/U7603: The H$`\alpha `$ agrees well with the Hi position velocity diagram. The Hi rotation curve seems to over-estimate the rotation velocity in the outer parts by a significant amount. For this galaxy we found a $`10`$ km s-1 difference between the systemic velocities of the Hi and H$`\alpha `$ curves. The velocity resolution of the Hi data is, however, 25 $`\text{km}\text{s}^1`$, so this difference may not be significant.

N5023/U8286: Excellent agreement between the data sets, except between $`50\mathrm{}`$ and $`100\mathrm{}`$ where the Hi curve overestimates the velocity. The good agreement in the inner parts of this edge-on galaxy confirm again that these late-type galaxies are virtually transparent.

U10310/Arp2: The H$`\alpha `$ data of this galaxy is rather poorer in quality than the previous couple of data sets. Nevertheless the data sets are in good agreement, except perhaps the Hi rotation curve which seems to overestimate the true rotation velocity.

DDO47/U3974: The H$`\alpha `$ curve is rather poor in quality, but shows reasonable agreement with the Hi curve.

DDO52/U4426: Rather poor H$`\alpha `$ data but a reasonable match with the Hi curve, except in the outer parts, where beam smearing effects presumable cause a small discrepancy.

DDO64/U5272: Good H$`\alpha `$ data show that the Hi curve suffers from slight beam smearing effects. Non-circular motions presumably cause the difference between approaching and receding sides at $`r30\mathrm{}`$.

DDO185/U8837: The H$`\alpha `$ data are of good quality, but do not correspond very well with the Hi position-velocity overlay nor the Hi curve. The most likely explanation is non-circular motions in especially the approaching side, where the DSS image shows what is likely a large star formation complex. The receding side does show a linear rise. We do find an H$`\alpha `$ systemic velocity that differs $`15`$ km s-1 from the Hi systemic velocity. Alternatively, as this is one of the galaxies where the continuum emission of the central part could not be detected, we may have chosen the wrong systemic velocity. This curve is therefore rather uncertain.

DDO189/U9211: The H$`\alpha `$ data are of average quality. The various data sets are all consistent with each other.

In summary the comparison between Hi and H$`\alpha `$ shows that the Hi curves in this sample tend to suffer from two systematic effects. The low-resolution curves clearly suffer from beam-smearing which underestimates the rotation velocities. This is apparent in e.g. U1230. A second effect is apparent for a small number of galaxies where the *raw* Hi and H$`\alpha `$ data are an excellent match, but where the Hi *curve* overestimates the true rotation velocity. Because of the good match between the raw data sets, this systematic effect must have been introduced as a result of beam-smearing corrections; U4325 is a good example.

This comparison underscores the need for this kind of high-resolution H$`\alpha `$ data even in moderately well-resolved galaxies. Beam-smearing corrections are useful in deriving the general properties of rotation curves. For example, they can tell us what the general shape of the original curve was, or how likely it is that a particular rotation curve shape can give rise to the observed data. When it comes to determining an accurate slope or detailed shape for an individual galaxy, beam smearing correction are frequently ambiguous.

## 5 Deriving hybrid rotation curves

In many cases the Hi data extend further out than the H$`\alpha `$ data. Due to the flattening of the Hi curve towards larger radii, resolution effects become less important, and to get a more extended rotation curve it therefore makes sense to combine the Hi and H$`\alpha `$ data. We have defined hybrid Hi/H$`\alpha `$ curves, always giving preference to H$`\alpha `$ data where available. In practice this meant that we simply adopt the H$`\alpha `$ rotation velocities interior to a radius $`r_{\mathrm{H}\alpha }`$, usually corresponding to the outermost point of the H$`\alpha `$ rotation curve, and the Hi velocities outside that radius. Table 3 lists these radii, along with supplementary notes.

### 5.1 Asymmetric Drift

In order to combine the Hi and H$`\alpha `$ rotation curves one should in principle correct for the pressure gradients in the gas. This correction for asymmetric drift is given by

$$v_c^2=v_\varphi ^2\sigma ^2\left(\frac{\mathrm{ln}\rho }{\mathrm{ln}R}+\frac{\mathrm{ln}\sigma ^2}{\mathrm{ln}R}\right),$$

where $`v_c`$ is the true rotation velocity, $`v_\varphi `$ the observed gas rotation velocity, $`\sigma `$ the velocity dispersion in the gas, and $`\rho `$ the volume density. The velocity dispersion in late-type galaxies is found to lie between 7 and 10 km s-1 and is to first order constant with radius. Assuming a constant scale height we then find that the correction for asymmetric drift is typically $`3`$ km s-1 or less. Given other uncertainties due to non-circular motions, inclination, etc. the asymmetric drift corrections are negligible and we have therefore not applied this correction.

### 5.2 Deriving smooth rotation curves

Some of the main assumptions that are made when deriving mass models from rotation curves is that there is symmetry, that all mass is on circular orbits and there is continuity with radius. The raw data still show some scatter and non-circular motions, which can occasionally result in virtual or ambiguous rotation velocities. In order to estimate the (smooth) radial run of the potential we want to represent the data by a curve whose velocity and velocity-derivative change smoothly, retaining small-scale details. As we are combining Hi and H$`\alpha `$ points we also want the density of the different data points to be approximately equal, in order to prevent the densely sampled H$`\alpha `$ curves from dominating the fits. We have used a local regression method to derive the smooth curve (Loader 1999). This method fits polynomial pieces *locally* to subsets of data but does not impose any *global* functional description (see dBMR for a more extensive description). As this process inevitably introduces a small degree of correlation between adjacent data points, the smooth (H$`\alpha `$) curves were resampled every 6′′.

The final errorbar in these resampled data points consists of two components: (i) the average measurement uncertainty, derived from the average weighted measurement error in each 6′′ interval; (ii) an error estimate due to large-scale asymmetries and non-circular motions. We used the difference between the weighted mean velocity and the smooth curve in each interval as an estimate. These two uncertainties were added quadratically giving a conservative estimate for the uncertainty.

A general problem with H$`\alpha `$ rotation curves is that the formal errorbars of well-sampled, high S/N data points become unrealistically small (sometimes less than 1 km s-1). These errorbars no longer have any physical significance but simply indicate that the Gaussian fit to the profile was well-determined. A similar problem also exists with Hi curves where the tilted-ring fitting programs produce error bars that give the uncertainty in the fit and *not* the actual uncertainty due to profile width etc. Realistically, observational and physical uncertainties make it difficult to determine a meaningful rotation velocity with an accuracy more than a few km s-1. As the final step in our derivation of the rotation curves we have imposed a minimum error of 4 km s-1 (prior to inclination correction) on each data point.

We compare the resampled smooth curves with the raw data in Fig. 7. To avoid cluttering the plot we represent the smooth curves as a line. The full curves, complete with data points, are also shown in Fig. 8. It is evident that the method does not introduce any systematic effects that affect the shape of the curve. For an extensive discussion see dBMR. The curves shown in Fig. 7 will be used as the basis for our mass modelling.

## 6 Mass models

In order to find the signature of the dark matter halo one needs to decompose the observed rotation curve in a number of separate dynamical components, described below.

### 6.1 Stellar and gas components

For the stellar disk the $`R`$-band photometry presented in de Blok et al. (1995) and Swaters (1999) (see also Swaters et al. 2002a) was used. For NGC 1560 we used the $`B`$-band photometry presented in Broeils (1992). For DDO47 and DDO64 we used the radial surface brightness profiles from Stil (1999). These profiles are presented as relative to the sky-level without absolute calibration and we therefore used the surface photometry presented in Hopp & Schulte-Ladbeck (1991) and Makarova (1999) to put the profiles on an absolute $`R`$-band surface brightness scale. (Since the Stil photometry is in the $`I`$-band, this assumes that there are no strong $`RI`$ gradients.) See Table 1 for more details.

The rotation curve of the disk was computed using the rotmod task in gipsy. The disk was assumed to have a vertical sech2 distribution with a scale height $`z_0=h/6`$ (van der Kruit & Searle 1981). The rotation curves of the stellar component were resampled at the same radii as the smooth curves. We assume $`M/L_{}`$ constant with radius, which is a reasonable assumption for the range of variation expected from plausible stellar populations (de Jong 1996).

For the gas component we used the Hi surface densities from the references given in Table 1. They were scaled by a factor of 1.4 to take into account the contribution of helium and metals. Their rotation curve was derived assuming the gas is distributed in a thin disk. For a few nearby galaxies we had to slightly smooth the Hi radial profile, mainly because small-scale structure at small radii cause unrealistically large fluctuations in the gas rotation curve and hence the potential. At large radii this small-scale structure is averaged out so the Hi profile there is a good description of the global Hi structure.

### 6.2 Dark halo

After subtraction of the above components from the observed rotation curve, any residuals are usually taken to indicate the presence of a dark matter halo. The halos are usually parametrised by fitting a model to these residuals. The quality of the fit tells us something about the applicability of any particular model. Here we restrict ourselves to two well-known models: the pseudo-isothermal halo, used in most of the classic rotation curve studies, and the CDM NFW halo.

#### 6.2.1 Pseudo-isothermal halo

We assume a spherical pseudo-isothermal halo with a density profile

$$\rho _{ISO}(R)=\rho _0\left[1+\left(\frac{R}{R_C}\right)^2\right]^1,$$
(1)
where $`\rho _0`$ is the central density of the halo, and $`R_C`$ the core radius of the halo. This density profile results in the rotation curve

$$V_{\mathrm{halo}}(R)=\sqrt{4\pi G\rho _0R_C^2\left[1\frac{R_C}{R}\mathrm{arctan}\left(\frac{R}{R_C}\right)\right]}.$$
(2)
The asymptotic velocity of the halo, $`V_{\mathrm{}}`$, is given by

$$V_{\mathrm{}}=\sqrt{4\pi G\rho _0R_C^2}.$$
(3)
To characterise this halo only two out of the three parameters $`(\rho _0,R_C,V_{\mathrm{}})`$ are needed, as equation (3) determines the value of the third parameter.

#### 6.2.2 NFW halo

The NFW mass density distribution (Navarro et al. 1996, 1997) takes the form

$$\rho _{NFW}(R)=\frac{\rho _i}{\left(R/R_s\right)\left(1+R/R_s\right)^2}$$
(4)
where $`R_s`$ is the characteristic radius of the halo and $`\rho _i`$ is related to the density of the universe at the time of collapse. These parameters are not independent and are set by the cosmology. They are usually expressed in a slightly different form by the concentration parameter $`c=R_{200}/R_s`$, and the rotation velocity $`V_{200}`$ at radius $`R_{200}`$. The latter is the radius where the density contrast exceeds 200, roughly the virial radius (Navarro et al. 1996). The rotation curve of the NFW halo then is given by

$$V(R)=V_{200}\left[\frac{\mathrm{ln}(1+cx)cx/(1+cx)}{x[\mathrm{ln}(1+c)c/(1+c)]}\right]^{1/2},$$
(5)
where $`x=R/R_{200}`$.

### 6.3 Stellar Mass-to-light ratios

One of largest uncertainties in any mass model is the precise value of $`M/L_{}`$. Though broad trends in $`M/L_{}`$ have been measured and modelled (Bottema 1997; Bell & De Jong 2000), the precise value from galaxy to galaxy is unknown, and depends on extinction, star formation history, initial mass function, etc. Because of this problem we present disk-halo decompositions using four different assumptions for $`(M/L)_{}`$ for those galaxies with photometry available. For the galaxies without photometry only the first (minimum-disk) model is presented.

1. *minimum disk.* This ignores all visible components and assumes that the observed rotation curve is entirely due to dark matter. This gives a strong limit on the amount of dark matter in a galaxy. This is what is usually called “minimum disk” in the CDM analyses.
2. *minimum disk + gas.* This is the classical definition of minimum disk mostly used in observational rotation curve papers (and usually called “minimum disk” there). The contribution of the gas is taken into account, but $`M/L_{}`$ is assumed to be zero.
3. *constant $`M/L_{}`$.* Here $`M/L_{}`$ is fixed to a constant and representative value for all galaxies. In this case we have chosen $`M/L_{}(R)=1.4`$ which is at the “blue” or light-weight end of the range of plausible values derived (under reasonable assumptions for IMF and star formation history) from population synthesis models using the observed colours as constraints. See dBMR for a further motivation of this value. By using this “lightweight” value we give the halo models maximum opportunity to fit the data.
4. *maximum disk.* Here the rotation curve of the stellar component is scaled to the maximum value allowed by the smooth curve, without exceeding it (van Albada & Sancisi 1986). In practice the stellar curve is scaled until the inner points match those of the smooth curve, but avoiding a hollow halo.

Each of the rotation curves was fitted using the gipsy task rotmas. The program determines the best-fitting combination of $`R_C`$ and $`V_{\mathrm{}}`$ (for the pseudo-isothermal halo) or $`c`$ and $`V_{200}`$ (for the NFW halo), using a least squares fitting routine.

## 7 Modelling Results

The results of the modelling are presented in Fig. 8 and Tables 4 and 5. Fig. 8 displays the models for the four $`M/L_{}`$ assumptions, or, if no photometry was available, only the minimum disk model. The NFW and pseudo-isothermal (“ISO”) halo fits are shown side by side. The values of the halo parameters are given in Table 4 for the NFW model and Table 5 for the pseudo-isothermal halo model. Both Figure and Tables also list the reduced $`\chi ^2`$ goodness-of-fit values. We have also done a simple $`\chi ^2`$ test to compute the chance $`p`$ that the data and the models could originate from the same distribution. We regard $`p>0.95`$ as a good fit, indicating that the data and the model match. Values of $`p<0.05`$ generally indicate that the wrong model is used to describe the data. Derived values for $`p`$ are shown in Fig. 8.

For the NFW fits the fitting routine often preferred very small or negative $`c`$ values. These are obviously unphysical, and we have restricted the parameter space for $`c`$ to values $`0.1`$. When the fit preferred a value less than 0.1 we set $`c`$ simply equal to 0.1, and found the best fitting value for $`V_{200}`$. This is indicated in Table 4 by giving the values in italics.

### 7.1 Comments on the mass modelling of individual galaxies

F563-1: The pseudo-isothermal halo fits are comparable in quality to the NFW fits. We can compare our results with the independent fits by dBMR. For the minimum disk pseudo-isothermal halo we find $`R_C=2.0\pm 0.2`$ and $`\rho _0=70.4\pm 13.1`$. dBMR find $`R_C=1.7\pm 0.2`$ and $`\rho _0=91.9\pm 21.6`$. These values thus agree within their $`1\sigma `$ errors. For the NFW models our values of $`c=7.8\pm 1.4`$ and $`V_{200}=106.8\pm 10.3`$ compare well with the dBMR values of $`c=10.7\pm 1.2`$ and $`V_{200}=93.1\pm 4.3`$. The $`c`$ values agree at the 1.1$`\sigma `$ level, whereas the $`V_{200}`$ values agree within the $`1\sigma `$ errorbars. Two completely independent analyses of different observations on different telescopes are thus able to reproduce the same parameters, showing that systematic effects due to telescope pointing, smoothing etc. are small.

U628: The pseudo-isothermal and NFW fits are comparable. One should note that for this galaxy no Hi data are available, and that the “minimum disk” and “minimum disk+gas” fits are identical. Also no gas rotation curve is included in the “constant $`M/L_{}`$” and “maximum disk” cases. The gas rotation curve does however usually have only a small impact on the fit parameters, and the values presented here will be close to the true ones.

U731: The pseudo-isothermal models and NFW fits are comparable except that the latter overestimate the rotation velocity in the inner parts. Both models have trouble fitting the sharp kink at $`r0.8`$ kpc. This kink is artificial and caused by a lack of data between 0.8 and 4.2 kpc due to the presence of a strong night-sky line at these velocities.

U1230: For both models maximum disk is the best fit. This comes at the cost, however, of a high $`M/L_{}(R)=6.1`$ which is hard to explain in term of reasonable population synthesis $`M/L_{}`$ values, which usually find values between $`0.5`$ and $`2`$ for LSB galaxies (van den Hoek et al. 2000).

U4173: This rotation curve is virtually identical to the Hi curve as the H$`\alpha `$ data do not add much extra constraint. Both models fit equally well. The Hi surface density profile given in Swaters (1999) gives rise to large fluctuations in the inner gas rotation curve and was therefore smoothed slightly as explained above.

U4325: Good fits for the pseudo-isothermal models, bad ones for the NFW models.

U5005: Both sets of fits are comparable. The data are however not of high enough quality to constrain either model.

U5750: This galaxy can again be compared with independent observations by dBMR. We first compare the fitting parameters. For the isothermal model we find $`R_C=5.0\pm 0.9`$ and $`\rho _0=7.9\pm 1.6`$. dBMR find $`R_C=4.3\pm 0.4`$ and $`\rho _0=10.6\pm 1.0`$ (again minimum disk). These agree within the $`1\sigma `$ level. For the NFW model we find $`c=1.9\pm 2.1`$ and $`V_{200}=145.7\pm 122.9`$, compared to their $`c=2.6\pm 1.5`$ and $`V_{200}=123.1\pm 58.8`$, again in good agreement. The large uncertainties in the NFW parameters indicate that the rotation curve of U5750 cannot be adequately described by that model. This is mainly caused by the very linear rise of the inner rotation curve. For the NFW model the only constraint is imposed by the two Hi points at the flat part of the curve. Without these the fit would diverge to $`c=0`$ and $`V_{200}=\mathrm{}`$. Clearly, higher-resolution Hi observations of the *outer* parts are needed further fill in the flat part of the curve and bring down the uncertainties of the fit parameters. This would however not change the conclusion that this galaxy is very hard to fit with the NFW model. Again we see that independent observations are capable of reproducing the same set of parameters at the $`1\sigma `$ level.

N1560: Due to the large number of data points and the small scale structure present in the rotation curve, the $`\chi ^2`$ values are high. The pseudo-isothermal and NFW models seem to perform equally well, though there is still a slight tendency of the NFW models to overestimate the inner parts. Whereas the problem with rotation curves is usually that there is not enough resolution, NGC 1560 illustrates the opposite: the non-circular motions in the H$`\alpha `$ data are large. For the pseudo-isothermal model maximum disk fits best, but again with a value $`M/L_{}(B)=4.9`$ which is much higher than one would find on the basis of reasonable star formation histories and stellar initial mass functions.

DDO185: This is a very linear curve that is reasonably fitted by pseudo-isothermal models, but poorly by NFW models. The curve is however only of average quality, and better observations are needed to constrain the models further.

DDO189: This curve is equally well fitted by pseudo-isothermal and NFW models.

DDO64: This curve shows a bump in the inner parts, that is very well fitted by the stellar disk in the maximum disk model, perhaps indicating that the inner part of this galaxy is dominated by stars. Similar bumps have been found in other highly resolved observations of dwarf galaxies (e.g. Blais-Ouellette et al. 2001). The pseudo-isothermal model is a better fit.

N2366: This curve is a challenge for both ISO and NFW models. The very linear rise in the inner parts rapidly changes in a flat part at larger radii. It is hard to reproduce such a sharp change in slope. It is possible that the Hi curve on which the outer points are based underestimates the rotation velocity as a comparison with the Hi position-velocity diagram perhaps suggests. Another possible explanation would be that non-circular motions due to the bar-like structure in the centre affect the curve. In their study of NGC 2366 Hunter et al. (2001) find some weak evidence for this from the Hi velocity field, but they also note that the rotation curve does not seem to be affected. Clearly there is room for additional high-resolution velocity measurements.

N3274, N4395, N4455: The NFW models fit as well as the ISO models. Parts of the rotation curve of N4395 are dominated by effects of star formation

U1281, U3137, U711, N100, DDO52, N5023, IC2233 and U10310: For these eight galaxies only minimum disk models are presented. For U711, N100, DDO52 and U10310 this is because no surface photometry is available. U1281, U3137, N100, IC2233, and N5023 are edge-on galaxies and deriving their face-on surface brightness profiles depends on various assumptions about the properties of their disk which would introduce additional uncertainties in the models.

### 7.2 Goodness of fit

The $`\chi ^2`$ values of the pseudo-isothermal models are generally smaller than those of the NFW models. This is illustrated in Fig. 9 where we plot both sets of $`\chi ^2`$ values against each other for the four different $`M/L_{}`$ scenarios. In all cases the NFW models are either comparable or worse than the isothermal models. We emphasise this point again in Table 6. Here we list the number of galaxies with good ($`p>0.95`$) and bad ($`p>0.05`$) fits for the two models under the four $`M/L_{}`$ assumptions.

We note that the $`\chi ^2`$ values given here should not be regarded as absolute. As described in Sect. 5, the combination of H$`\alpha `$ and Hi data, the symmetrization of the curves, the definition of the subsequent error and the imposition of a minimum error result in a conservative estimate of the uncertainty. The errorbars we give here are likely an overestimate of the true uncertainty. This at least partly explains the very low $`\chi ^2`$ values found for some curves.

The bottom panel in Fig. 9 compares the $`\chi ^2`$ values for the minimum disk case, but measured separately for the inner and outer halves of the rotation curves, respectively. The largest discrepancies between pseudo-isothermal and NFW are found in the inner parts, as expected, but it is clear that also in the outer parts pseudo-isothermal models generally provide better fits. It is remarkable that despite the large errorbars, there is still an obvious preference for the pseudo-isothermal model. This conclusion is independent of any over- or underestimate of the true error: a less conservative estimate would even strengthen the trends observed here.

The pseudo-isothermal model has no bad ($`p<0.05`$) fits, in contrast with NFW. The absolute number of good fits is a factor 2 larger for the pseudo-isothermal model. The proportion of “average” fits ($`0.05p0.95`$) is larger in the NFW case. The conclusion is that the pseudo-isothermal halo gives a better description of the data.

### 7.3 NFW halos

The mass and concentration of (numerically simulated) CDM halos depend on the cosmological assumptions that are used as input to the simulations. NFW halos follow a specific relation between $`c`$ and $`V_{200}`$ (Navarro et al. 1997). Independent simulations have also given a reasonable idea of the scatter one expects at fixed halo mass (Jing 2000; Bullock et al. 2001), even though they do not agree on the details. In Fig. 10 we compare our results with these predictions.

Shown as the full line is the prediction for a $`\mathrm{\Lambda }`$ dominated flat universe ($`\mathrm{\Omega }_m=0.3`$, $`\mathrm{\Omega }_\mathrm{\Lambda }=0.7`$). There is a general trend of increasing $`c`$ values towards lower $`V_{200}`$. There is some disagreement regarding the theoretical scatter in the $`c`$-values. Bullock et al. (2001) quotes a scatter in $`\mathrm{log}c`$ of 0.18, while Jing (2000) derives a smaller scatter in $`\mathrm{log}c`$ of 0.08. Here we adopt the larger value from Bullock et al. (2001) in order to give the CDM models as much leeway as possible. The dotted lines in Fig. 10 indicate the 1$`\sigma `$ (grey area) and $`2\sigma `$ scatter where $`\sigma =\sigma (\mathrm{log}c)=0.18`$. The small dots represent the values for other LSB and dwarf galaxies derived in dBMR. Our data show a similar behaviour, though the locus of our data seems to be better defined than that of the dBMR data. The spread in $`c`$ larger than the models predict, which is best illustrated by the excess of low-$`c`$ points.

## 8 Pointing effects

In an analysis of this nature it is crucial that the slit is aligned with the center of the galaxy. Any off-set from center will lead to a rotation curve that is less steep than a properly centered rotation curve unless the entire galaxy is a pure solid body rotator. It is therefore theoretically conceivable that systematic effects may have caused us to underestimate the inner slope of the rotation curve, thus mistaking the incorrectly lowered slope of the NFW profile for that of an isothermal halo.

Here we present several arguments why this is not the case, and why the observationally derived slopes are very close to the true slopes.

### 8.1 Evidence from observations

As already mentioned in Sect. 2, we took considerable care when acquiring the galaxy at the telescope. Most of the nearby galaxies in our sample were visible in the guiding-camera, and they were used to test our off-set procedure. This involved re-acquiring a galaxy several times in order to test the stability of the procedure, as well as moving back and forth to the acquisition star to test for repeatability. We found the stability of the system to be better than $`1^{\prime \prime }`$. Of course, for the faintest galaxies, we rely entirely on the position given for the center, which we took from Swaters (1999) or Stil (1999) or the NED database.

There are no systematic differences between the rotation curves acquired by us at the 1.93m telescope at OHP, by McGaugh et al. (2001) at the Kitt-Peak 4m telescope, and the Las Campanas du Pont 2.5 telescope, by Swaters et al. (2000) at the 200” Hale Telescope at Mt. Palomar and by Pickering et al. (1998) at the MMT. This indicates that either all or none of the data sets suffer from systematic effects due to mispointing. A number of galaxies have been observed multiple times by different observers, and similarly show no systematic differences. This is illustrated in Fig. 11 where the raw data of four LSB galaxy rotation curves are compared. In one case there are three independent data sets available, all agreeing with each other. It is unlikely that independent observers miss the center of the galaxy by the same amount in the same direction every time.

If the true rotation curves were really bona-fide NFW curves we would expect that with the increased resolution of the H$`\alpha `$ data, at least some of them should start to look more like NFW curves. After all, the combined samples in this paper and McGaugh et al. (2001) contain a few dozen rotation curves. Even in the presence of (hypothetical) pointing effects, some should have hit the nucleus and have shown the characteristic NFW shape. This has not happened. The situation is quite the reverse: the best resolved rotation curves show the flattest cores (see Sect. 9 and de Blok et al. 2001b).

This leads to an interesting contradiction: if indeed we imagine that a pointing offset leads to a systematic under estimate of the slope, we ought to see a trend of rotation curve shape with distance. This is not seen. A pointing offset of a given number of arcseconds corresponds to a different physical scale in each galaxy, depending on its distance. Our sample and that of McGaugh et al. (2001) span a range in distance of a factor of $`200`$. If we imagine that a pointing offset of, say, 0.5′′ causes us to mistake an NFW halo for a pseudo-isothermal one in the distant galaxies in our sample, then this implies that in order to mislead ourselves in a similar way, we ought to be making an error of almost 1.5 in our pointing when observing the nearest galaxies. Though this is an extreme example, it is highly improbably that systematic effects in target acquisition at the telescope depend on the distance of the target.

Assuming then that we managed in the large majority of the cases to indeed home in on the center of the light distribution of the galaxy, one could still assume that the dynamical center of the galaxy does not coincide with the center of the light distribution. Centering on the light distribution would thus not give us the true rotation curve. This has of course just the same effect of mis-pointing, but now the fault lies with the galaxy, rather than the observer.

To test this one needs high-resolution velocity fields, to which one can then fit a tilted ring model. The tilted ring model has the central position of the rings as a free parameter, and will thus immediately show whether there is a significant discrepancy between the optical and the dynamical center or not. This is precisely what Blais-Ouellette et al. (2001) have done. They have derived H$`\alpha `$ rotation curves for a number of late-type galaxies using full 2D Fabry-Perot velocity fields. They find the same slowly-rising rotation curves, and conclude that NFW models are incompatible with their data, which prefer a core-dominated halo or a model with a shallow inner density slope.

This is the same conclusion as reached by Bolatto et al. (2001), who obtained high-resolution $`5^{\prime \prime }`$ CO observations of the nearby dwarf galaxy NGC 4605. They derived a rotation curve from a 3D data cube, again using a tilted ring model, and find a slowly rising rotation curve incompatible with NFW.

A number of galaxies in our sample is edge-on, and it is conceivable that, despite contrary evidence, these data are affected by optical depth or projection effects. These edge-on galaxies are however only a small fraction of the galaxies investigated. For example, of the 56 galaxies presented here and in de Blok et al. (2001a), only 16 have inclinations larger than 70. Though there are no systematic differences between the $`30`$ per cent with $`i70^{}`$ and the $`70`$ per cent with $`i<70^{}`$, one could easily disregard the high-inclination galaxies without affecting any of the conclusions. There are no significant trends of $`c`$ or slope with inclination.

In summary, there is no evidence that the data set presented in this paper suffers from significant systematic effects, nor is it inconsistent with any other relevant data set in the literature. Pointing effects play no significant role, and the observed rotation curves represent the overall dynamics of our galaxies well.

### 8.2 Evidence from models

To quantify the effects of mispointing we have created model NFW velocity fields and have derived rotation curves for various pointing offsets. We consider model velocity fields assuming NFW halos with $`i=60^{}`$, $`V_{200}=100`$ km s-1 and $`c=4,9,20`$, respectively. Each of these three velocity fields is “observed” with a 1′′ wide slit, offset by $`d=0^{\prime \prime },2^{\prime \prime },5^{\prime \prime },10^{\prime \prime }`$ from the center, parallel to the major axis. We assume a distance of 10 Mpc, which is the average distance of the large majority of galaxies in our sample (disregarding the 4 galaxies with $`D>40`$ Mpc.)

The top panel in Fig. 12 shows the rotation curves, for each combination of $`c`$ and $`d`$. The most important question we want to answer here is under what conditions an NFW velocity field can produce a rotation curve that mimics a pseudo-isothermal curve.

The three bottom panels in Fig. 12 show the derived offset NFW curves, with the best fitting pseudo-isothermal models overplotted. Also shown are the differences between both types of curves. It is clear that in almost all cases the shape of the pseudo-isothermal halo curve is distinctly different from that of the NFW curve. The only cases where one might mistake an offset NFW curve for a pseudo-isothermal curve are those for ($`c=9`$; $`d=10^{\prime \prime }`$) and ($`c=20`$; $`d=5^{\prime \prime },10^{\prime \prime }`$). Pointing errors of this magnitude are simply not present in the data. If we accept an uncertainty of at most a few tenths of an arcsecond, it is easy to see that pointing effects are only likely to affect galaxies at distances $`100`$ Mpc.

We can also regard the NFW curves in Fig. 12 as best fits to the plotted pseudo-isothermal curves. In this case we can make a direct comparison with the fits plotted in Fig. 8, and see that the characteristic overprediction of the inner part of the curve is also present in the data. Except for the three cases mentioned above, Fig. 12 shows that in all cases the difference between NFW and pseudo-isothermal is most pronounced in the inner 2 kpc with residuals between $`5`$ and $`20`$ km s-1. The majority of galaxies in Fig. 8 does show such residuals, meaning that even in the presence of (hypothetical) modest pointing offsets, NFW curves would look different from the curves observed and shown in Fig. 8.

In summary, for the galaxies in this sample the effect of missing the centers of the galaxies would not be strong enough to masquerade NFW curves as pseudo-isothermal curves, except in the case of very (unrealistically) large offsets ($`5^{\prime \prime }`$). This conclusion, combined with the conclusion from Sect. 8.1 that any pointing errors must be below the $`1^{\prime \prime }`$ level, shows that the observed rotation curves must be close to the true rotation curves. Systematic observational effects cannot hide NFW halos in LSB galaxies.

## 9 Mass profiles

### 9.1 Introduction

In de Blok et al. (2001b) (hereafter dBMBR) the rotation curve sample presented in dBMR and part of the sample (the Jan00 galaxies) presented in this work were used to calculate the mass-density profiles that give rise to the observed rotation curves. As the different halo models predict different shapes of the mass distribution this gives one a direct test of the applicability of a particular model.

The mass density profile is derived by assuming a spherical mass distribution. Then, from $`^2\mathrm{\Phi }=4\pi G\rho `$ and $`\mathrm{\Phi }=GM/r`$ one can derive the mass density $`\rho (r)`$:

$$4\pi G\rho (r)=2\frac{v}{r}\frac{v}{r}+\left(\frac{v}{r}\right)^2.$$
(6)
In the above $`v`$ is the rotation velocity and $`r`$ is the radius.

This inversion is only valid if the contribution by the gas and stars is negligible, i.e. one implicitly assumes a situation (close to) minimum disk. It is well-established that this is a reasonable assumption for LSB (dwarf) galaxies. A minimum disk also provides an upper limit on the steepness of the slopes of the halo mass-density profiles: inclusion of gas and stars will necessarily tend to flatten the slopes.

dBMBR find that the shape of the minimum-disk halo profiles can usually be characterised by two power-law components of the form $`\rho r^\alpha `$. The outer slope has a value close to $`\alpha =2`$, whereas the inner slope is more shallow. The models tested here make distinct predictions regarding the value of the inner slope: the pseudo-isothermal halo predicts $`\alpha =0`$, the NFW halo predicts $`\alpha =1`$, and more recent CDM simulations (Moore et al. 1999) find even steeper values of $`\alpha =1.5`$.

dBMBR find that the distribution of the inner slopes $`\alpha `$ is asymmetric, with a well-defined peak at $`\alpha =0.2\pm 0.2`$, inconsistent with CDM predictions. This distribution has a broad wing however, extending to values of $`\alpha <1.8`$, seemingly consistent with the steep slopes demanded by CDM. dBMBR show that the inner slope one derives depends on the radius of the innermost sampled point $`r_{\mathrm{in}}`$ of the rotation curve. Larger values of $`r_{\mathrm{in}}`$ sample the range of radii where the NFW profile and the pseudo-isothermal profile have similar slopes (which does not imply they have similar *shapes*!). Small values of $`r_{\mathrm{in}}`$ probe the region where the predicted slopes are distinctly different. The galaxies with small values for $`r_{\mathrm{in}}`$ indeed show the clearest evidence for the presence of a constant-density core (i.e. $`\alpha 0`$). dBMBR conclude that all data are consistent with LSB galaxies having core-dominated halos with core-radii of a few kpc.

The implication is that galaxies with large values of $`r_{\mathrm{in}}`$ and consequently with steep slopes, should show much shallower slopes once the resolution is increased (and $`r_{\mathrm{in}}`$ decreased). The current data set puts us in an excellent position to test this prediction. For a large number of galaxies we now have high-resolution rotation curves, as well as independent lower-resolution curves from Swaters (1999). We should thus see a systematic decrease in inner slope when moving from the pure Hi data to the high-resolution data.

### 9.2 Results

#### 9.2.1 Minimum disk profiles

Fig. 13 plots the derived mass profiles for our complete sample. It is clear that most of the galaxies from our Feb01 run, which were specifically chosen to be nearby objects, are also characterised by an almost flat inner core with a radius of a few kpc, in contrast with the steep $`\alpha =1.5`$ power-law slope demanded by CDM, or even the inner $`\alpha =1.0`$ power-law slope of the NFW profile.

Table 7 lists the values of the slopes derived using the method described in dBMBR. In short, after determining the “break-radius” where the slope changes most rapidly, we determined the slope of the inner component using a weighted least-squares fit. The uncertainty was determined by re-measuring the slope twice, once by including the first data point outside the break-radius, and once by omitting the data point at the break-radius. The maximum difference between these two values and the original slope was adopted as the uncertainty. Following dBMBR we re-plot the values of the inner slope against the value of $`r_{\mathrm{in}}`$, this time also including the Feb01 data (see Fig. 14). The new galaxies are consistent with a core, not a cusp, and the galaxies with small values of $`r_{\mathrm{in}}`$ have flat mass density slopes.

#### 9.2.2 The case of N3274

The only exception seems to be N3274 which shows a slope close to $`1`$ even though its value for $`r_{\mathrm{in}}`$ is quite small. The pseudo-isothermal halo fit shows that this galaxy has the smallest (minimum disk) core-radius and highest central density of the entire sample. Though seemingly consistent with the NFW profile, it is also consistent with a pseudo-isothermal halo profile with a small core radius. Furthermore, it has by far the highest surface brightness and smallest scale-length of the entire sample, and in this case it is very likely that the minimum disk assumption breaks down and the stellar component *is* important in the inner parts.

#### 9.2.3 Constant $`M/L_{}`$ profiles

We tested this by re-deriving the slopes for the constant $`M/L_{}=1.4`$ case (see dBMR for a justification of this value). The halo rotation curve was found by quadratically subtracting the gas-rotation curve and the rotation curve of the stars. This treatment is likely to be too naive, as in a non-minimum disk case one expects the disk to influence the dark matter distribution to some degree (though perhaps not for LSB galaxies). A full treatment of this problem is beyond the scope of this paper, and the naive procedure suffices to illustrate the main point: as the stellar mass-to-light ratio is increased the inner slope of the halo mass-density profile becomes flatter.

The slopes for the constant $`M/L_{}`$ case are listed in Table 7. In Fig. 15 we again plot the derived slopes against the inner radii for the galaxies in our sample where a constant $`M/L_{}`$ model was available. A comparison with Fig. 14 shows that the data points have all moved up, as expected. In a number of cases we found rather large positive slopes, which would imply that these galaxies have hollow halos. This is rather improbable, and a more realistic explanation is that the value $`M/L_{}=1.4`$ is an overestimate of the true $`M/L_{}`$.

N3274 has a slope of $`0.66`$ in the constant $`M/L_{}`$ case, which is consistent with the slope one would expect for a halo with a core-radius of a few tenths of a kpc. Though this galaxy obviously cannot be used to prove or disprove either model, it is clear that galaxies which show steep inner slopes are likely to be of high surface brightness with inner regions dominated by stars. In order to unambiguously prove that a galaxy has a bona-fide NFW halo and is inconsistent with the pseudo-isothermal halo model, it needs to have a steep inner slope, *and* a small value of $`r_{\mathrm{in}}`$ *and* a low surface brightness. Such galaxies have not been found yet.

#### 9.2.4 Resolution and slope

One of the conclusions of dBMBR is that as the resolution of rotation curves ) is increased (i.e. $`r_{\mathrm{in}}`$ is decreased), the derived inner slope should become shallower. We test this by comparing the slopes and inner radii for a number of curves for which we have high-resolution curves and medium-resolution ($`15^{\prime \prime }`$) Hi curves. We determined the slopes of both the high-resolution curves and the Hi curves.

In Fig. 16 we plot the inversions of the high-resolution hybrid curves and the lower-resolution Hi curves. The slopes as derived from these curves are shown in the Figure. Those derived from the Hi curves tend to be steeper than those from the equivalent high-resolution curves. The values of the slopes in Fig. 16 are slightly different from those given in Fig. 13. This is because the Hi profiles generally did not extend as far inward as the hybrid profiles, and we had to choose a larger upper fitting radius. In Fig. 17 we plot the positions that the high-resolution and Hi curves occupy in the $`\alpha r_{\mathrm{in}}`$ diagram. The general trend is indeed that as the resolution is increased, the galaxies move towards smaller slopes, more or less along the pseudo-isothermal halo model curves.

It is interesting to compare our analysis with that of van den Bosch & Swaters (2001) (hereafter vdBS), who investigated the Hi curves of a number of galaxies for which we now have high-resolution data. Their conclusion was that the majority of the galaxies they investigated was consistent with a $`\mathrm{\Lambda }`$CDM scenario, but they also noted that it is difficult to distinguish between a cusp and a core model at the typical resolution of the Hi data. It is clear that this was because the lower spatial and velocity resolution makes it possible to force a NFW fit on the data, even though it is not always of the shape preferred by the data. We will discuss the galaxies they consider in turn.

U731: The Hi curve differs substantially from the high-resolution curve. Though our NFW fit gives similar parameters to the vdBS fit ($`c18`$ versus their $`c16`$), our curve provides a much inferior fit. This is mainly caused by the more linear inner curve. The inner slope of the high-resolution mass density profile shows a large change away from the optimum CDM value.

U3371: The high-resolution curve is more linear and rises less steeply than the Hi curve. The slope of the inner curve again shows a large change: the Hi value is $`\alpha 1`$ and consistent with CDM, the new value is $`\alpha 0`$ and consistent with a constant-density core. We were unable to derive a sensible NFW model fit.

U4325: This galaxy shows the same behaviour as U3371: the slope derived from the Hi data is optimum for CDM, whereas the slope derived from the high-resolution data is consistent with a constant density core.

N4395 (U7524): This is the best-resolved Hi curve in the vdBS sample, but here we also see the same trends noted above: the slope of the Hi curve is steeper than the slope of the high-resolution curve. The new curve also differs markedly from the Hi curve. The difference between NFW and pseudo-isothermal is most marked in the inner parts. The increase in resolution has made this galaxy less consistent with CDM.

N4455 (U7603): vdBS note that for this galaxy no meaningful CDM fit could be derived. We find identical slopes for both sets of data. As a note of interest, we note that the NFW fit to the Hi curve presented in vdBS shows the same systematic residuals as the H$`\alpha `$ LSB rotation curves presented in dBMR.

DDO189 (U9211): We find similar slopes as vdBS and also the models are of similar quality. The new data does not prefer one model over the other.

Of the 6 galaxies in this small subsample vdBS find that 5 are consistent with CDM, with one dubious case. The increased resolution suggests that only one of these 6 galaxies (DDO189) is still consistent with CDM (this is also the galaxy with the poorest high-resolution data in the subsample). For two galaxies CDM is perhaps inconsistent with the data (U731 and N4455), whereas for the remaining three the data are clearly inconsistent with the CDM models. In summary, we conclude that the slopes of rotation curves are best described by a pseudo-isothermal halo model.

## 10 Conclusions

We have presented high-resolution $`H\alpha `$/Hi rotation curves of a sample of 26 LSB and dwarf galaxies. We have fitted mass-models to these rotation curves assuming both a pseudo-isothermal (core-dominated) halo and a CDM NFW (cusp-dominated) halo. We find that the pseudo-isothermal halos generally provide better fits, though the difference is maybe not as pronounced as in dBMR, which can be traced back to the fact that our sample contains more galaxies of higher surface brightness than the dBMR sample.

We find more galaxies with low concentration parameters than predicted by numerical CDM simulations. An analysis of the mass-density profiles of the halos, as derived from minimum-disk rotation curves, shows that the galaxies in our sample are dominated by more-or-less constant density cores. As shown in dBMBR the few galaxies that show slopes consistent with CDM are usually the ones that are not well-resolved, so that one traces the edges of the constant-density cores rather than the cores themselves. We have illustrated this explicitly by comparing the mass-density profiles of a number of galaxies for which we have both high-resolution $`H\alpha `$ curves and medium-resolution Hi curves. In general the slopes derived from the high-resolution curves are less steep than those from the Hi curves.

In conclusion, our high resolution rotation data on nearby dwarfs and LSB galaxies show that the halos of late-type LSB and dwarf galaxies are dominated by kpc-sized constant density cores inconsistent with the predictions of cuspy dark matter halos in cosmological numerical simulations.

## Acknowledgements

EdB would like to thank ANSTO for their financial support which made attending the February 2001 observing run possible. AB thanks the Programme National Galaxies for financial support of his observing runs at the Observatoire de Haute Provence. We thank Stacy McGaugh for useful discussions.