Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.mso.anu.edu.au/~kcf/surya/bars/Perez05 Äàòà èçìåíåíèÿ: Tue May 31 21:04:22 2005 Äàòà èíäåêñèðîâàíèÿ: Sun Apr 10 06:04:34 2016 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: ï ï ï ï ï ï ï ï

Astronomy & Astrophysics manuscript no. aa2703-05 (DOI: will be inserted by hand later)

Dark matter in the inner parts of barred galaxies: The data
I. PÈrez1,3 , I. MÀrquez2 , K. Freeman3 , and R. Fux3,
1

,

4

2 3 4

Kapteyn Astronomical Institute, University of Groningen, Postbus 800, Groningen 9700 AV, The Netherlands e-mail: isa@astro.rug.nl Instituto de AstrofÌsica de AndalucÌa (CSIC), Apartado 3004, 18080 Granada, Spain RSAA, Mount Stromlo Observatory, Cotter road, Weston Creek, ACT 2611, Australia Geneva Observatory, Ch. des Maillettes 51, 1290 Sauverny, Switzerland

Received 14 January 2005 / Accepted 11 April 2005
Abstract. This paper presents surface photometry ( B, V, I , J, H, K ) and H rotation curves of 27 isolated spiral galaxies. The final goal is to obtain the mass distribution of a sample of isolated spiral galaxies in order to model their gas kinematics. This is then compared to the observed rotation curve, to determine the necessity of a dark halo in the inner parts (PÈrez et al. 2004). The azimuthally averaged radial surface brightness profiles and the integrated magnitudes obtained from ellipse fitting are given for each of the sample galaxies. The ellipse fitting technique applied to the light distribution also allowed us to obtain the size of the bar, and the inclination and position angle of the outer isophotes that allow the galaxy deprojection. Using these profiles, 1-D disk-bulge decomposition was performed to obtain the disk scale-length and the bulge effective radius for the different bands. Through the fitting of a parametric function to the observed rotation curve, the maximum rotational velocity and the corresponding radius was obtained. The correlation between the bulge and disk parameters is in agreement with previous studies (de Jong 1996a; MÀrquez & Moles 1999; Baggett et al. 1998). Regarding the Kormendy relation (Kormendy 1977), in agreement with de Jong, no correlation between the bulge effective radius and its surface brightness is found, possibly due to the small range of bulge magnitudes covered. We find a smaller scatter in the structural relations when compared to non-isolated samples in agreement with MÀrquez & Moles (1999). Finally, a correlation between the disk scale-length and the bar size is observed, possibly reflecting the rapid growth of a bar. Key words. galaxies: fundamental parameters ­ galaxies: structure ­ galaxies: photometry ­ galaxies: kinematics and dynamics

1. Introduction
The coupling between the massive dark halo and the visible component in spiral galaxies leads to a degenerate problem. The velocity field of an axisymmetric galaxy does not allow to uniquely disentangle the contribution of both components. The "maximum disk" assumption requires that the stellar disk accounts for most of the rotational support of the galaxy in the inner parts. There is no consensus yet on whether maximum disks are preferred (Freeman 1992; Courteau & Rix 1999; Salucci & Persic 1999). The match to the overall shape of the rotation curve is often met just as well by pure halo or pure disk models. It is not even clear whether the Milky Way favours a maximum disk or not (Lewis & Freeman 1989; Fux & Martinet 1994; Sackett 1997; van der Kruit 1999). To test whether the luminous mass in the inner parts of spiral galaxies can account for their observed gas kinematBased on observations obtained at the European Southern Observatory (La Silla), Las Campanas Observatories (Chile) and Siding Spring Observatory (RSAA, ANU, Australia). Figure 11 and Tables 3­8 are only available in electronic form at http://www.edpsciences.org

ics or whether a more axisymmetric dark matter component is required, the inner gas dynamics of a sample of barred spiral galaxies has been modelled by running a 3-D composite N -body/hydrocode on the potential derived directly from the light distribution. This approach takes into account the nonaxisymmetric structures, such as bars, that are very important in the dynamics of the inner region. These results are presented in a companion paper; PÈrez et al. (2004), hereafter Paper I. This paper presents the initial data set and sample characterization of a sample of isolated spiral galaxies from which the modelled galaxies presented in Paper I were chosen. This paper also presents the data analysis required to prepare the galaxies for the modelling of their stellar mass distribution prior to the dynamical simulations. Since the final aim is to investigate the distribution of luminous matter and test whether an additional dark component is needed to explain the observed kinematics, a good tracer of the stellar mass distribution is needed. Near infrared (NIR) observations have been chosen because they are good tracers of the underlying stellar population. The light in the NIR bands is dominated by the main sequence low mass stars, although there could also be an important contribution from giant stars.


2

I. PÈrez et al.: Dark matter in the inner parts of barred galaxies: The data

Furthermore, the NIR is less affected by internal extinction than the optical bands. Long slit observations were chosen as kinematic tracers. Integral field spectrographs offer an advantageous way of obtaining the kinematics of galaxies. However, long slit spectrographs have several features that make long-slit observations an interesting instrument for kinematic studies; it has more wavelength and spatial coverage than most 2-D spectrographs. The wavelength coverage is very useful to characterise the internal extinction of the galaxies and the large spatial coverage is useful to reach the regions well outside the bar region. Furthermore, very high signal-to-noise can be achieved for all the covered wavelengths. Different slit positions can help to understand the modelled 2-D velocity field. To obtain reliable mass-to-light ratios through population synthesis, is necessary to obtain photometry in different bands. The available stellar libraries are more complete for the optical bands so together with the NIR observations, further B, V and I band observations were obtained. Extinction information can be obtained from the Balmer line ratios which can help to determine the opacity of the disk as a function of radius. This information can be obtained from the spectroscopic observations. The colour maps can also help to trace the dust absorption: the profiles are expected to redden with dust absorption. Since the galaxies in the sample are too faint in HI to obtain high resolution gas maps, the 2-D dust maps will help tracing the gas and will be compared to the simulated gas distribution (Paper I). The selection of the sample will be described in Sect. 2. In Sect. 3 the observations and data reduction of the galaxies of the sample is explained. The methodology to prepare the images for the SPH modelling is presented in Sect. 4. The methodology and the correlation between the different parameters is presented in Sect. 5. Finally, the conclusions are given in Sect. 6

(see Fig. 1): all the isolated spiral galaxies in the volume given by cz 6000 km s-1 with galactic longitude(b) 40 , to avoid being affected much by the extinction of our galaxy, and those with inclination between 30 and 70 . These inclination angles were chosen to ensure that good measurement of the rotation curve was done; also, not too inclined galaxies were chosen to avoid internal extinction and to simplify the deprojection process. A preliminary list was created using the Third Reference Catalogue of bright galaxies (RC3), de Vaucouleurs et al. (1991), for the southern hemisphere. 2100 galaxies were chosen to be southern ( < 15 ) and volume confined (see above). From these, 293 were chosen to be isolated, using the criteria described above. However, non-classified objects could have still been found close to these galaxies. Therefore, all candidates were checked for companions with unknown redshifts using the Digital Sky Survey (DSS). Good quality data was obtained for 27 galaxies out of the final sample of 35 isolated galaxies. Up to this date, 5 barred galaxies have been modelled using hydrodynamical simulations (Paper I).

2.1. Characterisation of the sample
Figure 2 shows the number of galaxies as a function of morphological type. The morphological types were taken from the RC3 catalog and are based on the B-band; no attempt has been made to compare the morphologies in the different observed bands (apart from the presence of bars). Table 3 presents the main properties of the sample galaxies. From the sample of 35 galaxies, 18 galaxies (51%) present a bar (RC3 classification) in the optical band. From the final 27 galaxies, 8 have no information about the presence of a bar and 3 are classified as non-barred. After the ellipse fitting analysis (Sect. 5.3) and visual check, from the 8 galaxies with no bar information 3 showed clearly in the NIR or the I -band a bar morphology (UGC 10130, UGC 4861 and UGC 1553). From the 3 galaxies previously classified as non-barred, two showed bar morphologies in the NIR (NGC 6014 and ESO417-G006). With respect to AGN presence, 4 (11%) galaxies present an active nucleus. These numbers are found in larger samples. Due to the small number statistics these numbers are just indicative and are only useful to characterise the sample. One of the galaxies, NGC 6014, has been reclassified as a Sy2 galaxy using the Osterbrok diagnosis diagrams (Osterbrock 1989), only the [OII] emission line is not covered in the spectral range observed. The final sample of 27 galaxies does not differ in distribution of parameters to that of the initial 35 galaxies (all the AGN, except NGC 6014, are part of the final sample).

2. Sample selection
The properties of an initially unperturbed galaxy can be significantly modified by gravitational interaction with other galaxies. Such modifications, showing up either in the morphology or the dynamical properties, can be relatively mild (forming floculent like structures or tidal bridges) or so strong that they could lead to the complete disruption of the original system to form a different one, like an elliptical galaxy from the merging of two spiral galaxies (Barnes & Hernquist 1996). To avoid environmental effects where dark halos could be merged or overlapped it is necessary to study a well defined sample of isolated galaxies. It is, of course, not a simple question to define an isolated galaxy. Here, we define an operational criteria to ensure that possible perturbation effects have been erased, assuming that the typical time scale for the perturbation effect decay is no longer than a few times 109 years. For more detailed explanation about the selection process and isolation criteria see MÀrquez et al. (1999, 2000). We defined isolated systems as those without a companion within a projected radius of 0.5 Mpc or within the volume given by (cz) 500 km s-1 . In this way, one ensures that no encounter occurred between the galaxies in a time scale of the order of 109 yr. The initial sample was chosen as follows

3. Observations and reduction

3.1. Long-slit obser vations
The rotation curves were ob scopic (H) observations alon ing the Cassegrain Boller & 1.52 m telescope at La Silla Beam Spectrograph (DBS) on tained from optical spectrog several position angles usChivens spectrograph on the Observatory and the Double the 2.3 m telescope at Siding


I. PÈrez et al.: Dark matter in the inner parts of barred galaxies: The data

3

RC3: 16088 Spiral galaxies

Selection Filter

3473 Volume confined galaxies

cz < 6000km/s

2100 Southern galaxies

< 15

r<0.5MPc z< 600km/s 293 Candidates to be isolated

Barred 16%

*

Non-Barred* 84%

SBa/b SBb/c SBc/d SBd/m 13% 33% 21% 24%

S0/0+ Sa/b 12% 17%

Sb/c 32%

Sd/dm Irr 13%

35 Isolated galaxies After checking the DSS
Fig. 1. Flow diagram showing the selection criteria.

Spring Observatory (SSO). At La Silla, the detector was a 2K â 2K Loral-Lesser CCD with a pixel scale of 0.82 arcsec. The achived dispersion was 1.0 å pixel-1 using grating #26, with a wavelength coverage of 4900­5700 å. The detector at SSO was a SITe 1752 â 532 CCD. The gratings used were the 1200B and the 1200R for the blue and red arms, respectively. This set-up gives a dispersion of 1.13 å pixel-1 for the blue arm and 1.09 å pixel-1 for the red arm in the wavelength interval from 4315­5283 å and 6020­6976 å ; respectively, giving

a dispersion resolution 45 km s-1 . The slit was in both cases visually aligned with the galaxy nucleus at major and minor axis (RC3) and the slit width was set to 1.5 arcsec on the sky, giving a projected width of 2.2 pixels. The average seing was 1.2 arcsec. A comparison lamp exposure was obtained before and after each exposure to wavelength calibrate the frames. Spectrophotometric standards were observed with a slit of 10 arcsec. Total integration times for the galaxies varied from 1


4

I. PÈrez et al.: Dark matter in the inner parts of barred galaxies: The data Table 1. Calibration coefficients of the optical observations. Band B V I B V I B V I Zero-point ­21.19 ± 0.050 ­22.44 ± 0.190 ­21.67 ± 0.040 ­22.19 ± 0.050 ­22.25 ± 0.040 ­21.65 ± 0.044 ­22.40 ± 0.037 ­22.70 ± 0.068 ­22.14 ± 0.140 Colour coeff. June 1999 LCO 0.060 ± 0.020 ­0.080 ± 0.010 0.014 ± 0.010 January 1999 LCO 0.030 ± 0.011 ­0.050 ± 0.010 0.010 ± 0.010 June 2001 SSO ­0.037 ± 0.004 0.008 ± 0.008 ­0.025 ± 0.017 Airmass coeff. ­0.250 ± 0.030 ­0.230 ± 0.150 ­0.070 ± 0.030 ­0.210 ± 0.040 ­0.090 ± 0.030 0.008 ± 0.035 ­0.214 ± 0.032 ­0.111 ± 0.059 ­0.188 ± 0.122

Fig. 2. Morphological type distribution for the final sample of 27 galaxies.

to 3 h. In the case of the DBS the exposures had to be shorter than 20 minutes due to flexure of the instrument. The observations were split in three or more exposures per galaxy to eliminate cosmic rays from the frames. The observing log is presented in Table 4. All the spectra obtained were reduced using standard IRAF routines. Overscan and bias was subtracted from all the frames. No dark subtraction was done due to the low dark current of the chips used (1 e- pix-1 h-1 ). At the 1.52 m telescope, dome and sky flats were acquired in order to correct properly for slit width variation. Flatfielding correction was achieved to the 2% level. Geometrical corrections were applied to the frames in order to correct for any misalignment. The relatively small, compared to the slit size, size of the galaxies allowed for a correct sky subtraction for each frame prior to combining. A first order polynomial was fit along the spatial direction and then subtracted from the frames. In most of the cases, all the sky lines were well subtracted, but for some of the brightest lines some residuals persisted and had to be masked manually. The data reduction for the DBS spectra was carried out in a similar way, but due to the double beam nature of the spectrograph, care had to be taken to deconvolve the dichroic using spectrophotometric standards. Ionised gas rotation curves were obtained by fitting a Gaussian to the H and the [NII] (6583.4 å) emission lines using the MIDAS package ALICE. These lines were clearly detected in all the galaxies observed. The gas curves derived independently from H and [NII] are in good agreement at all radii. The H rotation curves are presented in Fig. 11.

servations at SSO were taken with the direct imager and the SITe 2K detector with a pixel scale of 0.6 arcsec pixel-1 . B, V and I filters were used, three images were taken per object per band, except that two frames were obtained in some cases. Total exposure times are shown in Table 5. Three standard fields (Landolt 1992) were acquired per night. A series of twilight sky flats were obtained each night. The average seeing was 1.2 arcsec for the LCO runs and 1.3 arcsec for the SSO runs. The optical data was reduced using standard IRAF routines. Bias and overscan subtraction was performed. The bias frames at the Swope telescope showed in the last columns in the CCD reading direction an increase of 1% in the number of counts, so the bias frame was subtracted from the observations. The bias levels were very stable during the runs. Calibration for the photometric nights1 was done using APPHOT on the Landolt fields and then using PHOTCAL to fit the zero-point magnitudes, the colour, and the extinction coefficients of the form: b = B + b1 + b2 ( B - V ) + b3 X v = V + v1 + v2 ( B - V ) + v3 X i = I + i1 + i2 ( B - V ) + i3 X where B, V , I are the standard magnitudes, b, v, i the instrumental magnitudes per second, X the airmass and b/v/i j the unknown transformation coefficients. In the cases where the Thuan-Gunn i filter were used, a transformation was applied from Thuan-Gunn to Cousins system using the transformation equation by Schombert et al. (1990). The resulting parameters are presented in Table 1.

3.3. Near-infrared imaging
NIR imaging data was obtained at Las Campanas Observatory, using the DuPont 100 inch telescope and the Swope 40 inch telescope; and also, the 2.3 m telescope at Siding Spring Observatory. The NIR observations at the DuPont were obtained with the IR camera with a Rockwell NICMOS3 HgCdTe 256 â 256 array with a pixel scale of 0.42 arcsec pixel-1 (1.8 arcmin â 1.8 arcmin field of view). The IR observations
1 The non photometric observations at LCO were repeated afterward during photometric conditions at SSO.

3.2. Optical imaging
Optical imaging observations were performed at the Swope 40 inch telescope at Las Campanas Observatory (LCO) and the 40 inch telescope at SSO. The observations at the Swope telescope were taken with the SITe 3 chip with a pixel scale of 0.43 arcsec pixel-1 and a field of view of 12 arcmin. The ob-


I. PÈrez et al.: Dark matter in the inner parts of barred galaxies: The data

5

at the Swope telescope were taken with the IR camera with a NICMOS3 HgCdTe 256 â 256 with a pixel scale of 0.6 arcsec pixel-1 (2.5 arcmin â 2.5 arcmin field of view). CASPIR, the imager at the 2.3 m telescope at SSO, uses a Santa Barbara Research Center 256 â 256 InSb detector array, with a focal plane scale of 0.5 arcsec pixel-1 . The filters used were J , H , Kn for the imaging with CASPIR and Jshort , H , Kshort for the LCO runs. The NICMOS3 detector becomes nonlinear when the total counts (sky + object) exceed 17 000 ADU. For this reason and the sky variability in the NIR, we exposed for 60­ 120 s in J 30 s in H and 15 s in K . Each object was observed at several positions on the array. The average seeing was 1 arcsec for the LCO runs and 1.2 arcsec for the SSO runs. Total exposure times and the observing log are presented in Table 6. For the NIR imaging obtained at LCO only the smallest angular size galaxies of the sample were observed. Therefore no offset to different sky positions was made and the sky was simply calculated by averaging a number of the stacked images for a given object with a clipping rejection algorithm. The CASPIR images were recorded at several positions on the array with a dither pattern and larger offsets were applied to obtain sky frames with the same exposure time as the object frames. The dithering between frames was typically 20 arcsec and the offset for the sky frames was typically of 2.5 arcmin. Sky frames were obtained after each second object exposure, beginning and ending each exposure sequence with a sky exposure. Each single exposure time was the same as for the Las Campanas runs. The NIR data was reduced using standard IRAF routines and eclipse (IR data reduction package developed at ESO). All the frames had to be corrected for instrumental effects such as nonlinearity, dark current and pixel-to-pixel variation. The observations at LCO were carefully made in order to avoid the nonlinear regime of the array, so no correction for nonlinear deviations was needed. The response of the CASPIR detector has a quadratic nonlinearity which must be corrected before dark correction and flatfielding is done. CASPIR data is linearised by subtracting a bias frame and then applying the quadratic correction. To better avoid the thermal emission of the telescope, series of on-lamp and off-lamp dome flats were used to flat-field the images in the case of the CASPIR data and good flat-fielding (around 1%) was achieved at LCO using twilight flats. The NIR background has contributions from OH airglow in the J , H and K bands, moonlight, especially in the J and H -band, and from thermal emission from the telescope and sky in the K band. To construct a sky image for the LCO data, we selected a number of neighbouring images taken within ±5 min, to ensure a sky variation not larger than 1% to 2%, from the given image; then, the images are combined with a clipping rejection algorithm, and finally, subtracted from the given image. For the CASPIR data we created a sky by subdividing each object observation sequence in 10 min observing blocks, creating the sky frame from the off-source images and then subtracting this sky from the on-source images of the subset. When no satisfactory subtraction was achieved the nearest off-source frames were used. Only in a few cases could the sky

Table 2. Calibration coefficients of the NIR observations. Band J H K Zero-point 21.55 ± 0.028 21.15 ± 0.020 21.11 ± 0.036 Color coeff. Du-Pont LCO ­0.069 ± 0.08 ­0.096 ± 0.01 0.094 ± 0.11 March 1999, Swope LCO ­0.082 ± 0.015 ­0.013 ± 0.061 0.030 ± 0.0256 May 2000, SSO ­0.034 ± 0.053 ­0.015 ± 0.090 0.092 ± 0.073 August 2000, SSO ­0.027 ± 0.062 ­0.013 ± 0.020 0.075 ± 0.064 Airmass coeff. ­0.094 ± 0.023 ­0.056 ± 0.025 ­0.075 ± 0.030

J H K J H K J H K

19.89 ± 0.032 19.45 ± 0.054 19.05 ± 0.04 21.35 ± 0.083 21.27 ± 0.104 20.02 ± 0.15 21.42 ± 0.098 21.15 ± 0.068 20.41 ± 0.069

­0.054 ± 0.020 ­0.155 ± 0.008 ­0.095 ± 0.040 ­0.054 ± 0.056 ­0.152 ± 0.094 ­0.20 ± 0.12 ­0.031 ± 0.042 ­0.08 ± 0.03 ­0.27 ± 0.243

values be obtained from the edges of the galaxy frames for the CASPIR data (eso416-G012, IC 2822 and NGC 7246). For flux calibration a number of standard stars from the HST standards list Persson et al. (1998) and the IRIS standard star list Carter (1995) were used. These lists provide stars fainter than 8th magnitude to avoid saturation of the detector during the acquisition process. 6 to 7 standard star observations were done per night, whenever possible repeating each standard up to 3 times. When comparing different photometric systems, transformation equations can be derived comparing common standard stars. To transform from the MSSSO system to the CTIO system the equations derived by McGregor (1994) were used. The difference between Kn and Kshort is minimal, giving a very similar band pass. The transformation between K and Kn was determined from model stellar atmosphere energy distributions and measured filter transmission curves (McGregor 1995, CASPIR manual): Kn = K - 0.022( J - K ). The NIR calibration coefficients for the different runs are presented in Table 2.

4. Preparing the images for the modelling

4.1. Star removal techniques
In order to model the gas flows using SPH, a smooth mass distribution map has to be input in the code. From the observed light distribution one can derive a mass model. However, the foreground stars in the field of the galaxies have to be removed so that the light from these sources does not affect the mass distribution and the SPH modelling. Furthermore, the ellipse fitting applied to the galaxies would also be affected by the light of these foreground stars. For the non-saturated stars, a point spread function (PSF) was calculated for each frame. The stars in each frame were found by using a detection algorithm which searches for the local density maxima with a given FWHM and


6

I. PÈrez et al.: Dark matter in the inner parts of barred galaxies: The data

Fig. 3. NGC 5728 I -band before and after star removal.

a peak intensity greater than certain value. Then, photometry is done on the found stars; and finally, the scaled PSF centered on each object is subtracted. In some cases, the residuals had to be cleaned using a linear interpolation algorithm. In the case of saturated stars, a model for each of the stars was created using a King profile for the halo and the star itself for the center. In most of the cases, the diffraction spikes were not cleanly removed and linear interpolation had to be done. For the finding stars algorithm, photometry and PSF matching, routines within DAOPHOT in IRAF were used. In the case where small background galaxies appeared, a box enclosing the galaxy was replaced by the pixel values of a box of the same size located at the same radial position at 180. Codes in IDL were written for the linear interpolation and the saturated star modelling. An example of the cleaning procedure is shown in Fig. 3.

inner regions of the disk. To combine both images first we have to create images of the same dimensions. Then, both frames are registered using common stars and making use of the GEOMAP task within IRAF to create the transformation function. GEOMAP determines higher order coordinate transformations (coordinates, rotation, pixel scale). GREGISTER transforms images from one coordinate system to another, applying the transformation map created by GEOMAP; then resampling was done using linear interpolation. Once the images have been centered and have the same pixel scale, the PSF of the two images are matched by convolving with a kernel computed directly from objects in the image with worst seeing. If f is the reference image, h is the best seeing image and k is the PSF matching function then by the convolution theorem; f = k h,

4.2. Combining H-band and I-band images
To model the stellar mass distribution, we would ideally prefer to use a photometric band that is not affected by dust. The bands that offer a good compromise between representing the old stellar population and being less affected by dust, the NIR bands, are nevertheless very time consuming to map due to the small size of the NIR arrays, to the fact that they are not as sensitive as the optical CCDs to obtain deep large images, and to the sky brightness limitation for NIR imaging. But, are galaxy disks optically thick all throughout their disks or are they optically thin in the outer regions and opaque or moderately optically thick in their inner parts?. Recent studies by Xilouris et al. (1999) present evidence showing that disks are optically thick in their inner parts and transparent in the outer parts in the I -band, for the inclinations encountered in this work (between 30 and 70 ). Therefore for our sample it is a good approach to combine NIR data for the inner parts of galaxies with I -band data for the outer parts. In order to obtain the mass distribution from the light down to the fainter isophotes, the I -band has been used for the outer optically thin regions while the H -band has been used for the

f~ ~ k= ~ h where denotes the Fourier transform. Once the images have been matched, a mask is created to isolate in both images the regions that will be combined. These regions are chosen using the extinction values obtained from the spectroscopy and choosing the H -band region to extend into the optically thin regions ( 1) where the observed colours do not change significantly. The noise in the H -band images is higher than in the I -band images. The SPH code deals with the different noise levels in the different parts of the images. Before the images were finally combined, they were also deprojected. This deprojection is done assuming a flat disk and then simple geometrical stretching is done. An example is shown in Figs. 4 and 5, where the masked images of NGC 7267 and the respective central line profiles are presented. One can note the fact that the H -band seems to be more symmetric and smoother; this is more evident if one looks at the combined image and the I -band image, due to the higher dust sensitivity in the I -band compared to the H -band.


I. PÈrez et al.: Dark matter in the inner parts of barred galaxies: The data

7

Fig. 4. On the parts. On the the outer part then combine

left the I -band image of NGC 7267 masked in the inner right, the H -band image of the same galaxy masked in s. Both masks are used on their respective images and into a single light distribution.

4

and 10.42 ± 0.06 and 10.32 ± 0.01 for NGC 7483). Two of the sample galaxies (NGC 6012 and NGC 5728) have integrated K and H magnitudes from previous works (MÀrquez & Moles 1999; de Jong 1996a; Glass & Moorwood 1985). These results are also in good agreement with the results obtained in this study. I -band contours, surface brightness profiles in all available bands, rotation curves and a plot of the ellipticity and the position angle are presented in Fig. 11. At the end of the section the methodology and results for the extraction of the parameters of the non-axisymmetric structures will be presented. In the inner parts the ellipticity and the position angle (PA) obtained from the ellipse fitting technique will be the indicators of the existence of non-axisymmetric structures such as bars and the position and size of these structures (MÀrquez et al. 1999, 2000). The velocity used in Sect. 2 is the asymptotic velocity obtained from fitting a parametric rotation curve to the observed one. The formula representing a galaxy by a sequence of flattened spheroids by Brandt (1960) and by Bettoni & Galetta (1997) was adopted, namely; v(r, ) = 3V
max

3

2

r/rmax 1 + 2(r/rmax )

3/2

P(i, PA - )

1 0 100 200 300 400

4

3

where, Vmax is the maximum velocity and the radius rmax is the radius at which this velocity is reached, i is the inclination, PA is the position angle of the line of nodes, is the observed position angle and P(i, PA - ) is the projection factor. Kannappan et al. (2002) derived a velocity value by fitting the empirical function: v = v0 + v
c

2

( 1 + x) (1 + x )1/



1 0 100 200 300 400

Fig. 5. Profiles of a cut along the major axis of NGC 7267 for the combined profile (top panel) and the I -band profile (lower panel). Notice the smooth transition in the combined image from the H -band region to the I -band region.

5. General galaxy parameters
To quantitatively describe the sample, some global parameters were extracted from the data. The total magnitude in the different bands were derived by integrating the flux within a certain radius corresponding to the point where the isophotal magnitude is 25 mag arcsec-2 in the B-band (see Table 7). For the galaxies for which the isophotes do not reach 25 mag arcsec-2 in the B-band, the surface brightness was extrapolated to the required magnitude using the exponential disk parameters obtained from the surface brightness profiles (Sect. 5.1). The total magnitudes have been checked both internally (measurements of the same galaxy at different runs and telescopes) and externally (measurements of the same galaxy by different authors). For the internal comparison the difference in H magnitude obtained at the 40 inch telescope at LCO and the magnitud