A&A 395, 753-759 (2002)
DOI: 10.1051/0004-6361:20021352

Spatially resolved spectroscopy of Coma cluster early-type galaxies

II. The minor axis dataset[*]

G. Wegner1,[*] - E. M. Corsini2 - R. P. Saglia3,4 - R. Bender3,4 - D. Merkl3,5 - D. Thomas3,4 - J. Thomas3 - D. Mehlert6


1 - Department of Physics and Astronomy, 6127 Wilder Laboratory, Dartmouth College, Hanover, NH 03755-3528, USA
2 - Dipartimento di Astronomia, Università di Padova, Vicolo dell'Osservatorio 2, 35122 Padova, Italy
3 - Universitätssternwarte München, Scheinerstraße 1, 81679 München, Germany
4 - Max-Planck-Institut für extraterrestrische Physik, Postfach 1312, 85741 Garching, Germany
5 - Max-Planck-Institut für Plasmaphysik, Boltzmannstraße 2, 85748 Garching, Germany
6 - Landessternwarte Heidelberg, Königstuhl, 69117 Heidelberg, Germany

Received 8 July 2002 / Accepted 10 September 2002

Abstract
We present minor axis, offset major axis and one diagonal long slit spectra for 10 E and S0 galaxies of the Coma cluster drawn from a magnitude-limited sample studied before. We derive rotation curves, velocity dispersion profiles and the H3 and H4coefficients of the Hermite decomposition of the line of sight velocity distribution. Moreover, we derive the line index profiles of Mg, Fe and H$\beta $ line indices and assess their errors. The data will be used to construct dynamical models of the galaxies and study their stellar populations.

Key words: galaxies: elliptical and lenticular, cD - galaxies: kinematics and dynamics - galaxies: stellar content -
galaxies: abundances - galaxies: formation


   
1 Introduction

This is the second of a series of papers aimed at investigating the stellar populations and the dynamics of early-type galaxies in the Coma Cluster. Spanning about 4 dex in density, the Coma cluster is the ideal place to investigate these galaxy properties as a function of the environmental density. In the first paper of the series (Mehlert et al. 2000, hereafter Paper I) we presented the galaxy sample, its selection properties, the photometric observations and the long-slit major axis spectra for 35 E and S0 galaxies of the Coma Cluster. Here we complement the spectroscopic database with minor axis, offset major axis and one diagonal long slit spectra for 10 objects of the sample. The combined dataset allows the construction of dynamical axisymmetric models of the objects, to study the properties of the dark matter halos of flattened, rotating E and S0 galaxies. This study is complementary to the one presented by Gerhard et al. (2001), focusing on round, non-rotating ellipticals.

The properties of the stellar populations of the central regions of Coma cluster galaxies have been investigated in detail by means of long-slit (Davies et al. 1987; Jørgensen et al. 1995; Lucey et al. 1997) and fibre-fed spectrographs (Jørgensen 1999; Mobasher et al. 2001; Poggianti et al. 2001a, b; Komiyama et al. 2002; Carter et al. 2002). The data presented here and in Paper I allow the study of the stellar population gradients for a restricted number of galaxies and the empirical verification of the aperture corrections (Jørgensen et al. 1995). Morover, possible systematic differences between the disk and bulge components of S0 galaxies can be investigated.

The spectroscopic observations and data reduction are described in Sect. 2. We present the derived kinematics in Sect. 3 and the line indices in Sect. 4. Conclusions are drawn in Sect. 5.


   
Table 1: Log of spectroscopic observations.
Object Run PA Position Offset Single Exp. Time Total Exp. Time Q
    [$^\circ$]   [''] [s] [h]  
(1) (2) (3) (4) (5) (6) (7) (8)
GMP 0144 4 1 MN - $5\times3600+1776$ 5.5 3
GMP 0282 3 135 MN - $3\times3600$ 3.0 2
GMP 0756 3 88 $\parallel$ MJ 5.8 N $6\times3600$ 6.0 3
  3 178 MN - $3\times4200$ 3.5 2
GMP 1176 1 78 MJ - $3\times1800$ 1.5 2
  1 78 $\parallel$ MJ 2.6 N $5\times3600$ 5.0 3
  1 168 MN - $4\times3600$ 4.0 1
GMP 1750 1 150 MN - $3\times3600$ 3.0 1
GMP 1990 1 135 $\parallel$ MJ 3.5 N $6\times3600$ 6.0 2
  1-2 45 MN - $3\times3600$ 3.0 2
GMP 3510 2 79 MN - $2\times4200$ 2.3 2
GMP 3792 1 127 MJ - $3\times1800$ 1.5 2
  4 127 MJ - $3\times1800$ 1.5 3
  1 37 MN - $3\times4200$ 3.5 1
GMP 3818 3 127 DG - $3\times1800$ 1.5 3
GMP 5975 4 113 MN - $3\times3600$ 3.0 3


Column 1: GMP No. from Godwin, Metcalfe & Peach (1983). Column 2: observing run. Column 3: slit position angle measured North through East. Column 4: slit position. MJ = major axis; $\parallel$ MJ = parallel to major axis; MN = minor axis; DG = diagonal axis. Column 5: northward offset of the slit with respect to galaxy center. Column 6: number and exposure time of the single exposures. Column 7: total exposure time. Column 8: estimated quality of the resulting spectrum. 1: very good; 2: good; 3: medium (see Fig. 1).

   
2 Observations and data reduction

2.1 Galaxy sample

All the observed galaxies (Table 1) belong to sample of 35 E and S0 galaxies of the Coma cluster studied in Paper I). For details about their morphological classification and relevant photometric properties (i.e. total magnitude, effective radius, mean surface brightness within effective radius, ellipticity at effective radius and luminosity weighted a4 parameter) the reader is referred to that paper. They were selected from the sample of Paper I as the objects with the most extended and precise major axis kinematics and therefore best suited for dynamical modelling, balancing the number of E and S0 types.

2.2 Long-slit spectroscopy

Long-slit spectroscopic data of the sample galaxies were obtained with the 2.4-m Hiltner telescope of the MDM Observatory at Kitt Peak, Arizona, USA, during four different runs between 1999 and 2001. Details of the instrumental set-up of the observations carried out on 16-20 April 1999 (run 1), 23-26 April 1999 (run 2), 26-30 May 2000 (run 3) and 9-12 February 2001 (run 4) are given in Table 2.
   
Table 2: Instrumental set-up of spectroscopic observations.
Parameter Run 1-2 Run 3-4
Spectrograph Modular Modular
Grating 1200 $\rm gr~mm^{-1}$ 1200 $\rm gr~mm^{-1}$
CCD "Charlotte'' SITe "Echelle'' SITe
Pixel number $1024\times1024$ $2048\times2048$
Pixel size $24\times24$ $\rm\mu m^2$ $24\times24$ $\rm\mu m^2$
Gain 3.2 $e^{-}{\rm ADU}^{-1}$ 2.7 $e^{-}{\rm ADU}^{-1}$
RON 5.5 e- 7.9 e-
Scale $\rm0\hbox{$.\!\!^{\prime\prime}$ }606\;pixel^{-1}$ $\rm0\hbox{$.\!\!^{\prime\prime}$ }606\;pixel^{-1}$
Dispersion $\rm 1.20\;\AA\;pixel^{-1}$ $\rm 1.01\;\AA\;pixel^{-1}$
Slit width $2\hbox{$.\!\!^{\prime\prime}$ }0$ $1\hbox{$.\!\!^{\prime\prime}$ }9$
Wavelength range 4916-5881 Å 4585-6470 Å
Instrumental FWHM 3.20 Å 3.11 Å
Instrumental $\sigma^{\rm a}$ 79 $\rm km~s^{-1}$ 77 $\rm km~s^{-1}$
Seeing FWHM $1\hbox{$.\!\!^{\prime\prime}$ }9$- $3\hbox{$.\!\!^{\prime\prime}$ }0$ $2\hbox{$.\!\!^{\prime\prime}$ }0$- $3\hbox{$.\!\!^{\prime\prime}$ }0$


$^{\rm a}$ Derived at the 5170 Å Mg triplet.


Minor-axis spectra were obtained for all the sample galaxies, while offset spectra with the slit parallel to the major axis were obtained only for GMP 0756, GMP 1176, GMP 1990. GMP 3818 was observed along a diagonal axis. In addition, we took spectra of GMP 1176 and GMP 3792 along their major axes to perform a consistency check with measurements of kinematics and line strength indices of Paper I. The typical integration time of the galaxy spectra was 3600 s. Total integration times and slit position angle of the galaxy spectra as well as the log of the spectroscopic observations are given in Table 1. At the beginning of each exposure the galaxy was centered on the slit using the guiding camera. In each run several spectra of giant stars with spectral type ranging from late-G to early-K were obtained to be used as templates in measuring stellar kinematics and line strength indices. The template stars were selected from Faber et al. (1985) and González (1993). Additionally we observed at least one flux standard star per night to calibrate the flux of the spectra before the line indices were measured. Spectra of the comparison Hg-Ne-Ar-Xe arc lamp were taken before and/or after object exposures to allow an accurate wavelength calibration. The value of the seeing FWHM during the observing run as measured by fitting a two-dimensional Gaussian to the guide star is given in Table 2.

2.3 Basic data reduction

All the spectra were bias subtracted, flatfield corrected, cleaned of cosmic rays, corrected for bad columns and wavelength calibrated using standard MIDAS[*] routines. The flatfield correction was performed by means of both quartz lamp and twilight sky spectra, which were normalized and divided into all the spectra, to correct for pixel-to-pixel sensitivity variations and large-scale illumination patterns due to slit vignetting. Cosmic rays were identified and eliminated by interpolating over as in Bender et al. (1994, BSG94 hereafter). The residual cosmic rays were eliminated by manually editing the spectra. The wavelength calibration was performed by means of the MIDAS package XLONG. Each spectrum was rebinned using the wavelength solution obtained from the corresponding arc-lamp spectrum. We checked that the wavelength rebinning had been done properly by measuring the difference between the measured and predicted wavelengths (Osterbrock et al. 1996) for the brightest night-sky emission lines in the observed spectral range. The resulting accuracy in the wavelength calibration is $\sim$0.08 Å corresponding to $\sim$ $\rm km~s^{-1}$ at 5170 Å. The instrumental resolution was derived as the mean of the Gaussian FWHMs measured for a number of unblended arc-lamp lines which were distributed over the whole spectral range of a wavelength-calibrated spectrum. The mean FWHM of the arc-lamp lines and the corresponding instrumental resolution derived at 5170 Å is given in Table 2. In the galaxy and the stellar spectra the contribution of the sky was determined by interpolating along the outermost 10''-20'' at the two edges of the slit, where the galaxy or stellar light was negligible, and then subtracted. A sky subtraction better than $1\%$ was achieved. All the spectra were corrected for CCD misalignment following BSG94. The spectra obtained for the same galaxy along the same axis were coadded using the center of the stellar continuum as a reference. This allowed to improve the signal-to-noise ratio (S/N) of the resulting two-dimensional spectrum. A one-dimensional spectrum was obtained for each kinematical template star as well as for each flux standard star. The spectra of the kinematical Lick-system templates were deredshifted to laboratory wavelengths.

   
3 Stellar kinematics

We measured the stellar kinematics from the galaxy absorption features present in the wavelength range and centered on the Mg line triplet ( $\lambda\lambda~5164,~5173,~5184$ Å) by applying the Fourier Correlation Quotient method (Bender 1990) as done by BSG94.

The spectra were rebinned along the dispersion direction to a natural logarithmic scale, and along the spatial direction to obtain a nearly constant $S/N \geq 20$ per resolution element. In a few spectra the S/N decreases to $\sim$10 at the outermost radii. The galaxy continuum was removed row-by-row by fitting a fourth to sixth order polynomial. The quality of the final spectrum depends on the resulting S/N. In Fig. 1 we show examples of central spectra covering the 3 quality classes listed in Table 1. The quality parameter is 1 for $S/N\geq100$, 2 for  $50 \leq S/N < 100$, 3 for 20 < S/N < 50.


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{MS2880f1.eps}\end{figure} Figure 1: Examples of central spectra covering the range of quality classes. Relative fluxes have false zero points for viewing convenience.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.3cm,clip]{MS2880f2a.eps}\hspace*{1cm...
....eps}\hspace*{1cm}
\includegraphics[width=8.3cm,clip]{MS2880f2d.eps}\end{figure} Figure 2: Kinematical parameters and line indices measured along the observed axes of the sample galaxies. For each axis the curves are folded around the nucleus and filled dots and asterisks refer to the two sides of the galaxy. The radial profiles of the line-of-sight velocity (v) after the subtraction of systemic velocity, velocity dispersion ($\sigma $), third (H3) and fourth (H4) order coefficient of the Gauss-Hermite decomposition of the LOSVD are shown in the left panels (from top to bottom). The radial profiles of the line indices H$\beta $, [MgFe], <Fe>, Mgb and Mg2 are plotted in the right panels (from top to bottom).
Open with DEXTER


 \begin{figure}
\par\includegraphics[width=8.3cm,clip]{MS2880f2e.eps}\hspace*{1cm...
...ps}\hspace*{1cm}
\includegraphics[width=8.3cm,clip]{MS2880f2h.eps}
\end{figure} Figure 2: continued.
Open with DEXTER


 \begin{figure}
\par\includegraphics[width=7.5cm,clip]{MS2880f2i.eps}\hspace*{1cm...
...ps}\hspace*{1cm}
\includegraphics[width=7.5cm,clip]{MS2880f2l.eps}
\end{figure} Figure 2: continued.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7cm,clip]{MS2880f3a.eps}\hspace*{1cm}
\includegraphics[width=7cm,clip]{MS2880f3b.eps}\end{figure} Figure 3: Kinematical parameters (v, $\sigma $, H3, and H4) and line indices profiles ($H\beta $, [MgFe], <Fe>, Mgb and Mg2) along the major axis of GMP 1176 (left panels) and GMP 3792 (right panels) measured during run 1 (open symbols) and run 4 (asterisks) compared to those obtained by Mehlert et al. (2000, filled symbols). For each run different symbols refer to the different sides of the galaxy.
Open with DEXTER

To measure the stellar kinematics of the sample galaxies we adopted HR 6817 (K1III) as kinematical template for runs 1, 2 and 3 and HR 3427 (K0III) for run 4 and we considered the wavelength range 5115-5541 Å (see Fig. 1) around the Mg lines of the galaxies. We derived for each galaxy spectrum the line-of-sight velocity distribution (LOSVD) along the slit and measured its moments, namely the radial velocity v, the velocity dispersion $\sigma $ and the values of the coefficients H3 and H4. At each radius, they have been derived by fitting the LOSVD with a Gaussian plus third- and fourth-order Gauss-Hermite polynomials ${\cal H}_3$ and ${\cal H}_4$, which describe the asymmetric and symmetric deviations of the LOSVD from a pure Gaussian profile (van der Marel & Franx 1993; Gerhard 1993). Errors on the LOSVD moments were derived from photon statistics and CCD read-out noise, calibrating them by Monte Carlo simulations as done by BSG94. In general, errors are in the range of 3-10  $\rm km~s^{-1}$ for v and $\sigma $, and of 0.01-0.04 for H3 and H4, becoming larger in the outer parts of some galaxies where for $10 \leq S/N < 20$. These errors do not take into account possible systematic effects due to template mismatch or the presence of dust and/or faint emission. The measured stellar kinematics are reported in Table 3 and plotted in Fig. 2.

Figure 3 shows the comparison between the measurements of v, $\sigma $, H3, and H4 along the major axis of GMP 1176 and GMP 3792 obtained here and measurements obtained in Paper I. The value derived from the different datasets are in agreement within the errors. However in the case of GMP 3792 the velocity dispersion measured in run 4 is systematically lower than that measured both in Paper I and in run 1.

   
4 Line indices

We measured the Mg, Fe, and H$\beta $ line strength indices following Worthey et al. (1994) from flux calibrated spectra, as done in Paper I. Spectra were rebinned in the dispersion direction as well as in the radial direction as before. We indicate the average Iron index with <Fe> $\rm = (Fe_{5270}+Fe_{5335})/2$ (Gorgas et al. 1990) and the usual combined Magnesium-Iron index with $\rm [MgFe] =
\sqrt{\rm Mg{\it b}<Fe>}$ (González 1993). We corrected all the measured indices for velocity dispersion broadening and calibrated our measurements to the Lick systems using stars from Faber et al. (1985). No focus correction was applied because atmospheric seeing was the dominant effect during observations (see Mehlert et al. 1998 for details). Errors on indices were derived from photon statistics and CCD read-out noise, and calibrated by means of Monte Carlo simulations. The measured values of H$\beta $, [MgFe], <Fe>, Mgb and Mg2 are listed in Table 4 and plotted in Fig. 2. Figure 3 shows the comparison between the measurements of H$\beta $, [MgFe], <Fe>, Mgb and Mg2 along the major axis of GMP 1176 and GMP 3792 obtained here and measurements obtained in Paper I. The agreement within the error is excellent. The data we obtained in this paper are fully consistent with those of Paper I. Indeed, Fig. 4 shows that the mean central values of $\sigma $, H$\beta $, [MgFe], <Fe>, Mgb and Mg2 measured along the minor axes of the sample galaxies within an aperture of 2'' and those obtained along the corresponding major axes are in agreement within the errors, assessing that the new galaxy dataset is on the Lick system.


  \begin{figure}
\par\includegraphics[angle=270,width=16.6cm,clip]{MS2880f4.eps}\end{figure} Figure 4: Central values of $\sigma $, H$\beta $, [MgFe], <Fe>, Mgb and Mg2 measured along major and minor axes and averaged within an aperture of 2''.
Open with DEXTER

   
5 Conclusions

We presented new radially resolved spectroscopy of 10 E and S0 galaxies of the Coma cluster. We derived rotation curves, velocity dispersion profiles and the H3 and H4 coefficients of the Hermite decomposition of the line of sight velocity distribution along the minor axis, the offset major axis and one diagonal direction. Moreover, we measured the line index profiles of Mg, Fe and H$\beta $line indices. The data complement the existing set along the major axis and have a precision and radial extent sufficient to construct flattened and rotating dynamical models of the galaxies and study their radially resolved stellar populations. Future papers of this series will address these issues.

Acknowledgements
We acknowledge the DFG grant SFB 375.

References

 


Copyright ESO 2002