Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://sn.sai.msu.ru/~sil/preprints/ic3653.pdf Äàòà èçìåíåíèÿ: Fri Jan 25 01:15:46 2008 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 19:37:48 2012 Êîäèðîâêà:

Mon. Not. R. Astron. Soc. 376, 1033­1046 (2007)

doi:10.1111/j.1365-2966.2007.11549.x

Kinematics and stellar populations of the dwarf elliptical galaxy IC 3653
I. V. Chilingarian,1
1 2

,2,3

P. Prugniel,

2,4

O. K. Sil'chenko1 and V. L. Afanasiev

5

Sternberg Astronomical Institute of the Moscow State University, Universitetsky pr. 13, Moscow 119992, Russia ´ ´ Universite de Lyon, Lyon F-69000, France; Universite Lyon 1, Villeurbanne, F-69622, France; Centre de Recherche Astronomique de Lyon, Observatoire de ´ ´ Lyon, 9 Av. Charles Andre, Saint-Genis Laval, F-69561, France; CNRS, UMR 5574; Ecole Normale Superieure de Lyon, Lyon, France 3 LERMA Observatoire de Paris-Meudon, 61 Av. de l'Observatoire, Paris 75014, France 4 GEPI Observatoire de Paris-Meudon, 5 place Jules Janssen, Meudon 92195, France 5 Special Astrophysical Observatory of the Russian Academy of Sciences, Nizhniy Arkhyz, Karachaevo-Cherkesia 369167, Russia

Accepted 2007 January 24. Received 2007 January 15; in original form 2006 August 3

ABSTRACT

We present the first 3D observations of a diffuse elliptical galaxy (dE). The good quality data (S/N up to 40) reveal the kinematical signature of an embedded stellar disc, reminiscent of what is commonly observed in elliptical galaxies, though similarity of their origins is questionable. Colour map built from Hubble Space Telescope Advanced Camera for Surveys (ACS) images confirms the presence of this disc. Its characteristic scale (about 3 arcsec =250 pc) is about a half of galaxy's effective radius, and its metallicity is 0.1­0.2 dex larger than the underlying population. Fitting the spectra with synthetic single stellar populations (SSP), we found an SSPequivalent age of 5 Gyr and nearly solar metallicity [Fe/H] =-0.06 dex. We checked that these determinations are consistent with those based on Lick indices, but have smaller error bars. The kinematical discovery of a stellar disc in dE gives additional support to an evolutionary link from dwarf irregular galaxies due to stripping of the gas against the intracluster medium. Key words: galaxies: dwarf ­ galaxies: elliptical and lenticular, cD ­ galaxies: evolution ­ galaxies: individual: IC 3653 ­ galaxies: stellar content.

1 INTR ODUCTION Giant elliptical galaxies (E) were believed to be simple objects, rotationally supported and featureless until the kinematical observations of NGC 4697 by Bertola & Capaccioli (1975) made realize that at least some of them were supported by anisotropic velocity dispersions (Binney 1976; Illingworth 1977). In parallel, new imaging and processing techniques revealed the presence of dust lanes (Bertola & Galletta 1978; Sadler & Gerhard 1985) and fine structures (e.g. shells; Malin & Carter 1983) in many Es. Diffuse elliptical galaxies (dE, also called dwarf elliptical or dwarf spheroidal galaxies) followed the same line more recently when high-quality images revealed fine structures, in particular presence of a spectacular and intriguing stellar spiral in IC 3328 (Jerjen, Kalnajs & Binggeli 2000) and IC 783 (Barazza, Binggeli & Jerjen 2002) or broader structures, like embedded discs or bars in a number of dEs (Barazza et al. 2002). Evidence for ubiquity of discs in bright dE galaxies in the Virgo cluster, based on multicolour photometry, has been shown by Lisker, Grebel & Binggeli (2006). A number of recent papers presented also long-slit spectroscopic data (De Rijcke et al. 2001; Geha, Guhathakurta & van der Marel 2002, 2003; Pedraz et al. 2002; Simien & Prugniel 2002; van Zee, Skillman & Haynes 2004a; van Zee, Barton & Skillman 2004b) revealing a diversity

E-mail: chil@sai.msu.su
C

of degree of rotational support and even some complex structures (De Rijcke et al. 2004; Thomas et al. 2006). The stellar population may also be inhomogeneous (Michielsen et al. 2003) and an ionized interstellar medium (ISM) has been detected in some objects (Michielsen et al. 2004). The complexity of elliptical galaxies put severe constraints on the scenario of their formation and evolution. Ellipticals are thought to form at high redshifts from collapse and hierarchical mergers and accidentally lately evolve through major mergers. For what concerns dEs, the main classes of physical mechanisms intervening in the formation and evolution are (i) the feedback of the star formation on the ISM and (ii) the environment. The stellar population of most present dEs may not have formed in large galaxies, because the metallicity, [Fe/H] -0.3 (Geha et al. 2003; van Zee et al. 2004b), is typically smaller than that of large bulges or Es. It is thought to have formed in small galaxies where the supernova-driven winds certainly play an important role in controlling star formation rate and metal enrichment. In the smallest galaxies, the major part of gas and produced metals will be spread to the intergalactic medium. However, duration of star formation should be long enough, because recent studies (Geha et al. 2003; van Zee et al. 2004b) demonstrate that dE galaxies exhibit solar [Mg/Fe], while in the case of short star formation episode they would have been iron-deficient. The environment also obviously plays a major role through three phenomena:

2007 The Authors. Journal compilation

C

2007 RAS


1034

I. V. Chilingarian et al.
most luminous dE in Virgo and has a relatively high surface brightness. It is located 2 7 from the centre of the cluster, i.e. 0.8 Mpc in . projected distance. Its radial velocity 588 ± 4 km s-1 (this work) confirms its membership to the Virgo cluster, the velocity difference from the mean velocity of Virgo (1054 km s-1 , HyperLeda; Paturel et al. 20031 ) is nearly -470 km s-1 . IC 3653 is located some 100 kpc in the projected distance from NGC 4621, a giant elliptical galaxy having a similar radial velocity value (410 km s-1 , HyperLeda). With other low-luminosity Virgo cluster members, in particular IC 809, IC 3652 for which the radial velocities have been measured, they may belong to a physical substructure of Virgo, crossing the cluster at 500 km s-1 . Velocity and velocity dispersion profiles from Simien & Prugniel (2002) show some rotation. Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) archival images from the Virgo cluster ACS survey (Cote et al. 2004) are also available and will be ^´ discussed here. In Section 2, we describe the observations and the data reduction. In Section 3, we derive the age and the metallicity using spectrophotometric indices, and in Section 4 we introduce and use a new method for measuring the internal kinematics and the population parameters by fitting directly the spectra against high spectral resolution population models. The analysis of the ACS images is presented in Section 5. Section 6 is the discussion. 2 SPECTR OSCOPIC OBSER VA TIONS AND D A T A REDUCTION The spectral data we analyse were obtained with the Multi-Pupil Fiber Spectrograph (MPFS) integral-field spectrograph. The MPFS, operated on the 6-m telescope Bolshoi Teleskop Al'tazimutal'nij (BTA) of the Special Astrophysical Observatory of the Russian Academy of Sciences, is a fibre-lens spectrograph with a microlens raster containing 16 â 16 square spatial elements together with 17 additional fibres transmitting the sky background light, taken 4-arcmin away from the object. The size of each element is 1 â 1 arcsec2 . We used the grating 1200 g mm-1 providing the reciprocal dispersion of 0.75 å pixel-1 with an EEV CCD42-40 detector in the spectral range 4100­5650 å. Observations of IC 3653 were made on 2004 May 24 under good atmosphere conditions (seeing FWHM = 1.4 arcsec). The total integration time was 2 h. The spectral resolution R = / , where is the FWHM (full width at half-maximum) resolution (width of the line-spread function), as determined by analysing twilight spectra, varied from R = 1300 to 2200 (between 2.5 and 3.3 å) over the field of view and wavelength. R is lower in the centre of the field and increases towards top and bottom and also in the red end of the wavelength range (Moiseev 2001). The following calibration frames were taken during the observations of IC 3653 with MPFS. (i) BIAS, DARK. (ii) `Etalon': 17 night-sky fibres illuminated by the incandescent bulb. This frames are used to determine positions of spectra on the frame. (iii) `Neon' (arc lines): by exposing the spectral lamp filled with Ar­Ne­He to perform a wavelength calibration. (iv) The internal flat-field lamp. (v) A spectrophotometric standard (Feige 56 for our observations) used to turn the spectra into absolute flux units.
1

(i) the ram pressure stripping against the intracluster medium (Marcolini, Brighenti & D'Ercole 2003), (ii) the tidal harassment due to distant and repeated encounters with other cluster members (Moore, Lake & Katz 1998) and (iii) the collisions involving large gas-rich objects that may disrupt gas clouds from which a dwarf galaxy may born (tidal dwarfs) and evolve into a dE after fading of the young stellar population (Duc & Mirabel 1999). The latter possibility may probably not be the main scenario because dEs are one order of magnitude more numerous than large galaxies which in average have suffered one major collision in their lifetime, while the merger remnants in the local Universe are only accompanied by a few tidal dwarfs candidates which have generally much smaller mass than bright ellipticals (Weilbacher et al. 2000). The first two mechanisms are undoubtedly the major ingredients. Ram pressure stripping is an accepted explanation for the H I depletion of spiral galaxies in clusters and any type of gas-rich galaxies will experience it. It may well explain the relation between the morphology of dwarfs and the density of the environment (Binggeli, Tarenghi & Sandage 1990; van den Bergh 1994): low-density environment is mostly populated by dIrrs, while dEs reside mostly in clusters where ram pressure stripping is expected to be the most efficient. Recently, Conselice et al. (2003) discovered H I in two dEs in the Virgo cluster, which they interpret as transition objects on their way in the morphological transformation from dIrr to dE. Conselice, Gallagher & Wyse (2001) argue that dEs are recently accreted galaxies which morphologically evolve from late-type galaxies when they cross the centre of the cluster. There is yet no decisive test to disentangle the roles of each mechanism. In this paper, we present 3D spectroscopic observations in order to bring further observational constraints. Velocity fields and spatial distribution of the stellar population are most needed to check if counterparts of the observed kinematical substructures can be detected. We are presenting here the first 3D observations of a dE. IC 3653 is a bright dE galaxy belonging to the Virgo cluster (Binggeli, Sandage & Tammann 1985). In Table 1, we summarize its main characteristics. IC 3653 was chosen because it is amongst the
Table 1. General characteristics of IC 3653. Sersic exponent, ´ kinematical and stellar population parameters are obtained in this paper, other properties are taken from HyperLeda and Goldmine (Gavazzi et al. 2003; http://goldmine.mib.infn.it/) data bases, and from Ferrarese et al. (2006). Uncertainties given for age and metallicity correspond to the measurements on co-added spectra. Name Position B Distance modulus A(B) M(B)corr Spatial scale Effective radius, Re B , mag arcsec-2 Ellipticity, S´ ersic exponent, n Heliocentric cz, km s-1 cent , km s-1 Vmax , km s-1 Vmax / t, Gyr (lum. weighted) [Z/H], dex (lum. weighted) IC 3653, VCC1871 J124115.74+112314.0 14.55 31.15 0.13 -16.78 82 pc arcsec-1 6.7 arcsec 550 pc 20.77 0.12 1.2 588 ± 4 80 ± 3 18 ± 2 0.27 ± 0.08 5.2 ± 0.2 -0.06 ± 0.02

http://leda.univ-lyon1.fr/
C

C

2007 The Authors. Journal compilation

2007 RAS, MNRAS 376, 1033­1046


Kinematics and stellar populations of IC 3653
(vi) A standard for Lick indices and radial velocity (HD 137522 and HD 175743), used also to measure instrumental response: asymmetry and width of the line-spread function. (vii) `SunSky': twilight sky spectra for additional corrections of the systematic errors of the dispersion relation and transparency differences over the fibres.
Table 2. Parameters of the `3-points' binning: number of spatial elements, mean AB magnitude, mean AB surface brightness (mag arcsec-2 ) and mean signal-to-noise ratio at 5000 å. Bin P1 P2 P3 N
spax

1035

m(AB) 16.3 15.2 15.7

(AB) 18.7 19.9 20.9

S/N 69 49 21

2.1 Data reduction We use the original IDL software package created and maintained by one of us (VLA) that we modified for this work: error frames are created using photon statistics and then processed through all the stages to have realistic error estimates for the fluxes in the resulting spectrum. The primary reduction process (up to obtaining flux-calibrated data cube) consists of: (i) Bias subtraction, cosmic ray cleaning. Cosmic ray cleaning implies the presence of several frames. Then they are normalized and combined into the cube (x, y, Num). The cube is then analysed in each pixel through `Num' frames. All counts exceeding some level (5 ) are replaced with the robust mean through the column. Then the individual cleaned frames are summed. This procedure assumes that all the frames have the same atmospheric conditions and cleans only the brighter spikes, though it was found sufficient for our purpose. (ii) Creation of the traces of spectra in the `etalon' image. Accuracy of the traces is usually about 0.02 or 0.03 pixel. (iii) Flat-field reduction and diffuse light subtraction. Flat-field is applied to the CCD frames before extracting the spectra. The scattered light model is also constructed and subtracted from the frames during this step. It is made using parts of the frames not covered by spectra and then interpolated with low-order polynomials. (iv) Creation of the traces for every fibre. On this step the traces are determined for each fibre in the microlens block (presently 256 fibres) using the night-sky fibre traces created on the second step and interpolation between them using the tabulated fibre positions. (v) Spectra extraction. Using the fibre traces determined in the previous steps, spectra are extracted from science and calibration frames using fixed-width Gaussian (usually with FWHM = 5 pixel for the present configuration of the spectrograph). The night-sky spectra are also extracted from the science frames. (vi) Creation of dispersion relations. Spectral lines in the arc lines frame are identified and dispersion relations are computed independently for every fibre. (vii) Wavelength rebinning. All the spectra of night sky, object and standard stars are rebinned independently into logarithm of wavelength (as required for the kinematical analysis). The sampling on the CCD varies between 0.65 and 0.85 å (FWHM resolution 3­5 pixel) and we rebinned to a step of 40 km s-1 , i.e. 0.55­0.75 å, corresponding to the mean oversampling factor of 1.2 (we checked that this oversampling was high enough to have no measurable effect on the result of the analysis). (viii) Sky subtraction. The sky spectrum, computed as the median of the 17 night-sky fibres, is rescaled using flat-field to account for the twice larger aperture of the night fibres compared to object fibres. Then it is subtracted from each fibre. (ix) Determination of the spectral sensitivity and absolute flux calibration. Using the spectrophotometric standard star, the ratio between counts and absolute flux is calculated and then approximated with a high-order polynomial function over the whole wavelength
C

9 77 122

range. Finally, the values in the data cube are converted into F [erg cm-2 s-1 å-1 ]. 2.2 Spatial adaptive binning



In our data, the surface brightness, B , changes from 18.5 mag arcsec-2 in the centre down to 21.2 mag arcsec-2 in the outer parts of the field of view. Though the signal-to-noise ratio at the central part is high enough (30) for detailed analysis, the accuracy of kinematics and Lick indices measurements in individual decreases dramatically at the largest radii. To increase the signal-tonoise ratio in the outer parts, without degrading the resolution in the centre, we applied the Voronoi adaptive binning procedure proposed by Cappellari & Copin (2003) for this exact purpose. It uses variable bin size to achieve equal signal-to-noise ratio over the field of view. The Voronoi 2D-binning produces a set of 1D spectra that are then analysed independently. For the kinematical analysis, we will use a target signal-to-noise ratio of 15, and for stellar population analysis we will use 30. Besides we will also use a tessellation of the data set containing only three bins (3-points binning hereafter): central condensation (3 â 3 arcsec2 region around the centre of the galaxy), elongated discy substructure (14 â 7 arcsec2 ) oriented according to kinematics (see Section 4, Fig. 1, illustrating locations of bins and demonstrating spectra integrated in them) with the central region excluded, and the rest of the galaxy. The location of these three bins is shown in Fig. 1, and their characteristics in Table 2. Such a physically stipulated tessellation allows us to gain high signal-to-noise ratios in the bins in order to have high-quality estimations of the stellar population parameters in the regions where they are expected to be homogeneous. 3 SSP A G E AND MET ALLICITY DERIVED FR OM LICK INDICES A classical and effective method of studying stellar population properties exploits diagrams for different pairs of Lick indices (Worthey et al. 1994). A grid of values, corresponding to different ages and metallicities of single stellar population models (instantaneous burst, SSP), is plotted together with the values computed from the observations. A proper choice of the pairs of indices, sensitive to mostly age or metallicity like H and Mgb, allows us to determine SSP-equivalent age and metallicity. We use a grid of models computed with the evolutionary synthesis code: PEGASE.HR (Le Borgne et al. 2004). These models are based on the empirical stellar library ELODIE.3 (Prugniel & Soubiran 2001, 2004) and are therefore bound to the [Mg/Fe] abundance pattern of the solar neighbourhood (see Wheeler, Sneden & Truran 1989). To show that this limitation is not critical for our (low-mass) galaxy, Fig. 2(a) presents the Mgb versus Fe diagram with the models

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046


1036

I. V. Chilingarian et al.

(a) (c)

(b)

Figure 1. (a) The radial velocity field with three bins overplotted (3-points binning, see the text); (b) fit of the `P1' bin, representing the spectrum, `PEGASE.HR' marking the template spectrum (shifted by -0.3 on the flux axis from its real position, and `residuals' showing the difference between the fit and observed spectrum; (c) `P1', `P2' and `P3': co-added spectra for the 3-points binning.

by Thomas, Maraston & Bender (2003) for different [Mg/Fe] ratios overplotted. These data allow us to conclude that IC 3653 has solar [Mg/Fe] abundance ratio with a precision of about 0.05 dex. A critical limitation of Lick indices is their sensitivity to missed/wrong values in the data, for example, due to imperfections of the detector, or uncleared cosmic ray hits. A simple interpolation of the missed values is a poor solution, because if some important detail in the spectrum, e.g. absorption line, is affected, the final measurement of the index will be biased. In addition, the definition of the indices (see equations 1­3 in Worthey et al. 1994) does not allow us to flag or to decrease the weight of low-quality values. Due to a defect of the detector, our data have a 3-pixel wide bad region (hot pixels) in the middle of the blue continuum of Mgb. Therefore, strictly speaking, we could not measure Mgb at all, neither H on a significant part of the field of view.
C

As a workaround, we replaced all the missing or flagged values in the data cube by the corresponding values of the best-fitting model determined as explained in the next section. Therefore, instead of a mathematical interpolation of missing values, we are using model predictions. We measured Lick indices as defined in Trager et al. (1998), and we also derived the combined iron index Fe = 0.72Fe5270 + 0.28Fe5335 , and `abundance-insensitive' [MgFe] = Mg Fe (Thomas et al. 2003). A good intermediate resolution age tracer, H +Mg+Fe125 (Vazdekis & Arimoto 1999) cannot be used, because the required signal-to-noise ratio of about 100 at around = 4340 å cannot be achieved even after co-adding all the spectra in the data cube due to lower S/N in the blue end of the spectral range. The statistical errors on the measurements of Lick indices were computed according to Cardiel et al. (1998). Table 3 presents the measurements of selected Lick indices for the 3-points binning,
2007 The Authors. Journal compilation
C

2007 RAS, MNRAS 376, 1033­1046


Kinematics and stellar populations of IC 3653
Table 3. Measurements of selected Lick indices for the 3-points binning. All values are in å. The first line for each index corresponds to the measurements made on the real spectra and the second on the best-fitting optimal templates (see the text). (Note that the Fe indices Fe4531 and Fe5015 are essentially sensitive to titanium; see Sil'chenko & Shapovalova 1989.) Name Ca4227 G4300 Fe4383 Ca4455 Fe4531 Fe4668 H Fe5015 Mgb Fe5270 Fe5335 Fe5406 Fe [MgFe] Bin 1 1.062 1.096 5.106 4.995 5.794 4.861 1.187 1.336 2.858 3.499 6.070 5.114 1.841 1.908 5.121 5.491 3.575 3.361 3.016 3.059 2.708 2.675 1.687 1.840 2.930 2.952 3.236 3.150 ± 0.080 ± 0.135 ± 0.176 ± 0.089 ± 0.125 ± 0.181 ± 0.067 ± 0.138 ± 0.065 ± 0.072 ± 0.084 ± 0.065 ± 0.075 ± 0.070 Bin 2 0.874 1.020 5.042 4.820 4.940 4.458 1.096 1.235 2.326 3.367 5.776 4.602 1.823 1.966 4.767 5.211 3.550 3.203 2.969 2.886 2.551 2.524 1.667 1.720 2.852 2.785 3.182 2.987 ± 0.190 ± 0.313 ± 0.389 ± 0.186 ± 0.260 ± 0.361 ± 0.122 ± 0.246 ± 0.116 ± 0.132 ± 0.155 ± 0.121 ± 0.138 ± 0.127 Bin 3 0.689 1.092 6.923 4.979 6.251 4.696 0.750 1.313 1.595 3.457 5.685 4.808 1.800 1.878 4.912 5.292 3.771 3.334 3.362 2.963 2.526 2.598 1.542 1.779 3.128 2.861 3.435 3.089 ± 0.700 ± 1.040 ± 1.243 ± 0.561 ± 0.777 ± 1.025 ± 0.304 ± 0.584 ± 0.278 ± 0.309 ± 0.367 ± 0.286 ± 0.325 ± 0.301

1037

Table 4. Comparison of age and metallicity measurements for the 3-points binning obtained with pixel fitting and based on different pairs of Lick indices grids: measured on real spectra and on best-fitting templates. `P1' tfit , Gyr tH ­Mgb tH - Fe tH ­[MgFe] t t t
H ­Mgbmod H - Fe mod H ­[MgFe]mod

`P2' 4.95 6.97 7.28 6.88 ± ± ± ± 0.30 1.47 2.20 1.89

`P3' 4.97 ± 0.70 11.25 ± 6.02 6.08 ± 3.94 6.70 ± 5.11 4.85 ± 6.02 4.23 ± 3.94 4.15 ± 5.11 -0.17 ± 0.05 -0.34 ± 0.24 0.03 ± 0.13 -0.11 ± 0.15 -0.16 ± 0.24 -0.17 ± 0.13 -0.16 ± 0.15

4.93 7.04 7.02 7.11

± ± ± ±

0.20 1.56 1.33 1.65

5.27 ± 1.56 5.15 ± 1.33 5.22 ± 1.65 0.03 ± 0.01 -0.05 ± 0.09 -0.02 ± 0.04 -0.04 ± 0.05 -0.02 ± 0.09 0.02 ± 0.04 0.00 ± 0.05

5.13 ± 1.47 4.30 ± 2.20 4.12 ± 1.89 -0.14 -0.18 -0.10 -0.14 ± ± ± ± 0.02 0.08 0.05 0.06

Zfit , dex Z H ­Mgb Z H - Fe Z H ­[MgFe] Z H ­Mgbmod Z H - Fe mod Z H ­[MgFe]mod

-0.14 ± 0.08 -0.14 ± 0.05 -0.13 ± 0.06

The median value for age using [MgFe]­H pair is 6 ± 2.5 Gyr, without significant spatial variation. [MgFe] and Fe indices are not very age sensitive, thus the age estimations depend mostly on values of H , and are almost equal for all three pairs of indices.

4 STELLAR POPULA TIONS AND INTERN AL KINEMA TICS USING PIXEL FITTING Various methods have been developed to determine the star formation and metal enrichment history (SFH) directly from observed spectra (Ocvirk et al. 2006a,b; De Rijcke et al. 2004; Moultaka et al. 2004). The procedure that we are proposing here, population pixel fitting, is derived from the penalized pixel-fitting method developed by Cappellari & Emsellem (2004) to determine the line-of-sight velocity distribution (LOSVD). The observed spectrum is fitted in pixel space against the population model convolved with a parametric LOSVD. The population model consists of one or several starbursts, each of them parametrized by some of their characteristics, typically age and metallicity for a single burst while the other characteristics, like initial mass function, remain fixed. A single minimization returns the parameters of LOSVD and those of the stellar population. Ideally, we would like to reconstruct SFH, over all the life of the galaxy. This means, disentangle internal kinematics and distribution in the Hertzsprung­Russell diagram from the integratedlight spectrum. This problem has been discussed in several places (e.g. Ocvirk et al. 2006a,b; De Rijcke et al. 2004), it is clearly extremely degenerated and solutions can be found only if a simplified model is fitted. In this paper, we discuss only the simplest case of SSP characterized by two parameters: age and metallicity. We do not discuss complex SFH, because signal-to-noise ratio of our data is not sufficient. The 2 value (without penalization) is computed as follows: 2 =
N

and Fig. 2(c) presents the main index­index pairs that are used to determine age and metallicity. We see almost no population difference among three bins within the precision we reach. The age (see Table 4) is around 6 Gyr, metallicity is about solar for `P1' and slightly subsolar for `P2' and `P3'. The large scatter of the measurements, seen in Fig. 2, results in a spread of age estimations from 4 to 13 Gyr. It is mainly caused by an imperfect cleaning of the spikes and dark pixels but also weak nebular emission may alter these indices: H index might be affected by emission in H , Mgb ­ by [N I] ( = 5199 å) laying in the red continuum region. Though we do not see any significant emission line residuals when we subtract the best-fitting model, we cannot exclude completely this effect. To reduce the scatter, we measured Lick indices on the optimal templates fitted to the data. This approach may produce biased results in case of model mismatch, due for example to inconsistent abundance ratios between the models and the real stellar population. However, we believe it is not the case, since IC 3653 exhibits solar Mg/Fe abundance ratio (see Fig. 2a). We made inversions of the bi-index grids for three combinations of indices: Mgb­H , [MgFe]­H and Fe ­H . The maps shown in Fig. 3 represent interpolated values of the parameters between intensity-weighted centres of the bins. The metallicity distribution shows slight gradient from -0.15 dex at the periphery to +0.10 in the very centre (the average error bar on the metallicity measurements using [MgFe]­H is 0.15 dex).
C

{ Fi - P1 p [Ti (t , Z ) L(v, , h 3 , h 4 ) + P2q ]}2 , Fi2

(1)

where L is LOSVD; Fi and Fi are observed flux and its uncertainty; Ti is the flux from a SSP spectrum, convolved by the

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046


1038

I. V. Chilingarian et al.
(a) (b)

(c)

(d)

Figure 2. The Mgb­ Fe (a), H ­[MgFe] (b), H ­Mgb (c) and H ­ Fe (d) diagrams. On (a) the models from Thomas et al. (2003) are plotted. On (b), (c) and (d) the displayed grid is constructed of the values of the Lick indices for PEGASE.HR synthetic spectra (SSP) for different ages and metallicities. Bold crosses with pointers represent measurements for three regions of the galaxy (see the text), thin crosses are for individual bins for Voronoi tessellation with target S/N = 30.

line-spread function of the spectrograph (LSF; see the next section); P1p and P2q are multiplicative and additive Legendre polynomial of the order of p and q for correcting a continuum; t is the age, Z is the metallicity, v , , h3 and h4 are radial velocity, velocity dispersion and Gauss­Hermite coefficients, respectively (van der Marel & Franx 1993). Normally we used no additive polynomial continuum, and fifth-order multiplicative one, and for IC 3653, which has a low velocity dispersion resulting in insufficient sampling of the LOSVD, we did not fit h3 and h4 . The problem can be partially linearized: in particular, fitting of additive polynomial continuum, and relative contributions of subpopulations Ti is done linearly on each evaluation of the non-linear functional. Thus we end up with 10 free parameters: t, Z, six coefficients for Pmult5 , v and .
C

The main technical part of our method is a non-linear minimization procedure for 2 difference between the observed and modelled spectra. The latter is made by interpolating a grid of SSP spectra computed with PEGASE.HR and degraded to match the LSF of the observed spectrum. This grid has 25 steps in age (10 Myr to 20 Gyr) and 10 steps in metallicity ([Fe/H] from -2.5 to 1.0). Because the minimization procedure requires that the derivatives of the functions are continuous, we used a two-dimensional spline interpolation. The non-linear minimization is made with the MPFIT package (by Craig B. Markwardt, NASA2 ) implementing constrained variant of Levenberg­Marquardt minimization.

2

http://cow.physics.wisc.edu/craigm/idl/fitting.html ~
C

2007 The Authors. Journal compilation

2007 RAS, MNRAS 376, 1033­1046


Kinematics and stellar populations of IC 3653
(a) (b) (c) (d)

1039

Figure 3. Maps for age and metallicity obtained by inverting bi-index grids for indices measured on optimal template spectra: (a) Mgb­H , (b) [MgFe]­H , (c) Fe ­H ; and on real data: (d) Fe ­H .

4.1 Line-spread function of the spectrograph Before comparing a synthetic spectrum to an observation, it is required to transform it as if it was observed with the same spectrograph and setup, i.e. to degrade its resolution to the actual resolution of the observations. Actually the spectral resolution changes both with the position in the field of view and with the wavelength (thus it is not a mere operation of convolving with the LSF). Taking into account these effects is particularly critical when, as it is the case here, the physical velocity dispersion is of the same order or smaller than the instrumental velocity dispersion. The procedure for properly taking into account the LSF goes in two steps. First, determine the LSF as a function of the position in the field and of the wavelength. Secondly, inject this LSF in the grid of SSP. Therefore we made an exhaustive analysis of the LSF of our observations. A previous study of the change of the resolution of the MPFS over the field of view (Moiseev 2001) qualitatively agrees with our results. To measure the LSF change over the field of view we used the spectra of standard stars (HD 135722 and HD 175743) and twilight sky that we analysed with our fitting procedure. The highresolution spectra ( = 0.55 å; R 10 000) for the corresponding stars (the Sun for the twilight spectra) taken from the ELODIE.3 library (Prugniel & Soubiran 2001, 2004) were used as templates. Since these spectra have exactly the same resolution as the PEGASE.HR SSPs, the `relative' LSF that we determined in this way can be directly injected to the grid of SSP to make it consistent with the MPFS observations. We parametrize the LSF using v , , h3 and h4. The whole wavelength range (4100­5650 å) was split into five parts, overlapping by 10 per cent, and the LSF parameters were extracted in each part independently in order to derive the wavelength dependence of the LSF.
C

Finally, to inject the LSF in the grid of SSPs, we applied the following steps to every fibre. (i) Five convolved SSP grids were created using the LSF measured for the five wavelength ranges. (ii) The final grid was generated by linear interpolation at each wavelength point between the five grids of SSPs. For each fibre of the spectrograph, this produces one grid of models matched to the LSF of the observations. 4.2 Results: kinematics, age and metallicity maps We applied Voronoi adaptive binning procedure to our data, setting the target signal-to-noise ratio to 15. The resulting tessellation includes 76 bins with sizes from 1 to 12 pupils. To get a better presentation one may interpolate the computed values of each parameter over the whole field of view using the intensity-weighted centroids of the bins as the nodes. We measured the systemic radial velocity 588 ± 5 km s-1 . The uncertainty includes possible systematic effects not exceeding 4 km s-1 . The maps of radial velocity and velocity dispersions are presented in Figs 4(d)­(f). The galaxy shows significant rotation and highly inclined disc-like structure. The uncertainties of the velocity measurements were estimated using Monte Carlo simulations and confirmed by 2 maps obtained fitting only multiplicative polynomial continuum on a grid of values of age, metallicity and velocity dispersion (see the Appendix for details). They depend on the signal-to-noise ratio and change from 2.5 km s-1 for the S/N = 30 to 8 km s-1 for the S/N = 10. The velocity dispersion distribution shows a gradient from 45­50 km s-1 near the maxima of rotation to 75 km s-1 in the core. Note a sharp peak of the velocity dispersion of 88 km s-1 slightly shifted to the south-west of the photometric core. The uncertainties of the velocity dispersion measurements are 3.8 km s-1 for the S/N = 30, 5.5 km s-1 for the S/N = 20, and 11 km s-1 for the S/N = 10.

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046


1040

I. V. Chilingarian et al.

(a)

(b)

(c)

(d)

(e)

(f )

Figure 4. The generalized view of the kinematical and stellar population data. (a) intensity map with the positions the kinematical profiles are presented for: (1) position of slit in Simien & Prugniel (2002), (2) and (3) positions of major and minor axes of the embedded stellar disc; (b) map of the luminosity-weighted age distribution (in Myr); (c) map of the luminosity-weighted metallicity distribution ([Z/H], dex); (d) radial velocity field, interpolated between the nodes of the Voronoi tessellation; (e) radial velocity field: values exactly correspond to the Voronoi tessellae; (f) velocity dispersion field, interpolated.

The rotation velocity and velocity dispersion profiles are shown in Fig. 5 (top pair). The rotation curve is raising to 20 km s-1 at the edge of the field (i.e. at a radius of 5 arcsec). The previous study of IC 3653 was made using long-slit spectroscopy (Simien & Prugniel 2002) under poor atmosphere conditions (6 arcsec seeing). After the proper degrade of the spatial resolution of the MPFS data, the agreement with these earlier observations, both for radial velocity and velocity dispersion profiles, is excellent (Fig. 5, bottom pair). 5 PHO T OMETR Y AND MORPHOLOGY FR OM A C S IMA GES We have used ACS images from the HST archive, proposal 9401, `The ACS Virgo Cluster Survey' by Patrick Cote. The photometric ^´ analysis of IC 3653 is not included in Cote et al. (2004). We have ^´ converted ACS counts into corresponding ST magnitudes according to the ACS Data Handbook, available online on the web site of Space Telescope Science Institute, and calibrated into AB magnitudes. The derived colours agree perfectly with those given in Ferrarese et al. (2006). The F475W­F850LP colour map (equivalent to SDSS g­z, Fig. 6) reveals an elongated region with a major axis of about 7 arcsec, axis
C

ratio 3.5, and orientation coinciding with the kinematical disclike feature, which is redder than the rest of the map. This colour map was obtained using Voronoi 2D binning technique applied to the F850LP image in order to reach the signal-to-noise ratio of 80 per bin. In the age and metallicity maps derived from the MPFS data, this elongated feature is seen as a metallicity gradient because of the lower spatial resolution due to the adaptive binning. The F475W­ F850LP colour difference between the nucleus and the surrounding region 5-arcsec away from the centre is about 0.05 mag. Using PEGASE.2 templates (Fioc & Rocca-Volmerange 19973 ) and ACS filters transmission curves (Sirianni et al. 2005), we find that this colour difference corresponds to a metallicity difference of about 0.1 dex (assuming an age of 5 Gyr as found from our spectroscopic analysis), a value consistent with the spectroscopic analysis. We have fitted two-dimensional Sersic profile using the GALFIT ´ package (Peng et al. 2002). We can see significant positive residuals, representing the nucleus in the very centre (around 1.5 arcsec in size, central surface brightness ST475 16.25 mag arcsec2 , slightly asymmetric and offcentred with respect to the centre of the Sersic ´
3

http://www2.iap.fr/users/fioc/PEGASE.html
C

2007 The Authors. Journal compilation

2007 RAS, MNRAS 376, 1033­1046


Kinematics and stellar populations of IC 3653

1041

Figure 6. F475W­F850LP colour map (HST ACS data). Isophotes of the F850LP image are overimposed. The position and size of the MPFS field of view are shown.

Figure 5. Top pair of plots: the radial velocity profile for the slice `2' (major axis of the embedded disc) and velocity dispersion profiles for the slices `2' and `3' (major and minor axes of the embedded disc). Bottom pair of plots: comparison of the kinematical profiles (slice `1', major axis of main galactic body) to Simien & Prugniel (2002): radial velocity and velocity dispersion.

profile having n = 1.88, Re = 6.9 arcsec and = 0.11 (Sersic in´ dex, effective radius and ellipticity, respectively; our values coincide with ones from Ferrarese et al. 2006). There are faint large-scale residuals as well, that can be explained by superposition of several components (at least two). Then we have modelled the images by elliptical isophotes with free centre and orientation. We see some isophote twist and ellipticity change in the inner region of the galaxy that may be due to the embedded disc. Main parameters of the model fitted are presented in Fig. 7. The subtraction of the model from the original image does not reveal any internal feature. On the lower right-hand plot in Fig. 7, the F475W light profile is shown with crosses. The solid line represents the best-fitting S´ ersic profile for the whole galaxy excluding only very centre (inner 1 arcsec) with n = 1.9, and dashed line gives the best fitting (n = 1.2) for the periphery of the galaxy outside 7 arcsec, i.e. beyond the disc found in the colour map. We consider this latter fit as the most representative. 6 DISCUSSION Both line-of-sight velocity field extracted from the MPFS data cube and colour map obtained from the HST imagery provide coinciding arguments for a presence of a faint internal corotating stellar disc
C

Figure 7. Photometrical characteristics of IC 3653: position angle (upper left), ellipticity of the isophotes (lower left), discy/boxy parameter a4 (upper right) in two colours and photometric profile in F475W.

embedded within a rotating spheroid. This is the main observational result from our study of IC 3653, which may be regarded to some extend as an edge-on counterpart of IC 3328, the dE with embedded spiral structure found by Jerjen et al. (2000). The disc of IC 3328 is detected to a larger radius than the one of IC 3653: the spiral structure is seen up to 30 arcsec = 2re , while the disc of IC 3653 is detected only to 3 arcsec = 0.5re . The contrast of the spiral in IC 3328 is a few per cent and in IC 3653 the disc is seen from the colour and only as a weak maximum of ellipticity at the radius of its maximum contrast to the bulge. In both galaxies, the discs are essentially stellar (no sign of star formation) and old. Estimating and comparing the mass of the two discs is strongly model-dependent, and beyond the scope of this study. In this section, we will first compare the results of two methods for estimating stellar population parameters: Lick indices and pixel fitting, then the characteristics of IC 3653 to other dE galaxies.

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046


1042

I. V. Chilingarian et al.
the dense cluster environment from late-type disc galaxy (dIrr progenitor). Gas was removed by means of ram pressure stripping and star formation was stopped. This process must have finished at least 5-Gyr ago, otherwise we would have seen a younger population in the galaxy. However, duration of the star formation period must have been longer than 1 Gyr, otherwise deficiency or iron (i.e. Mg/Fe overabundance) would have been observed. Within this scenario, metallicity excess in the disc can be explained by a slightly longer duration of the star formation episode compared to the spheroid. However, we cannot see the difference in the star formation histories (even luminosity-weighted age), because of insufficient resolution on the stellar population ages. Another way to acquire a disc having higher metallicity than the rest of the galaxy is to experience a minor dissipative merger event (De Rijcke et al. 2004). This is rather improbable for a dwarf galaxy, but cannot be completely excluded. In particular, kinematically decoupled cores, recently discovered in dwarf and low-luminosity galaxies (De Rijcke et al. 2004; Geha, Guhathakurta & van der Marel 2005; Prugniel et al. 2005; Thomas et al. 2006) can be explained by minor merger events. The most popular scenario which is usually considered to explain formation of embedded stellar discs in giant early-type galaxies is star formation in situ after infall of cold gas on to existing rotating spheroid, e.g. from a gas-rich companion (cross-fuelling). This scenario has been used by Geha et al. (2005) to explain counterrotating core in NGC 770, a dwarf S0 which is more luminous than IC 3653 (MB = -18.2 mag) and located close to a massive spiral companion NGC 772 (MB = -21.6 mag) in a group. In NGC 127 (MB = -18.0 mag), a close satellite of another giant gas-rich galaxy, NGC 128, we observe the process of cross-fuelling at present (Chilingarian 2006). Group environment, where relative velocities of galaxies are rather low, favours of the interaction processes on large time-scales, such as smooth gas accretion. Our data on IC 3653 cannot provide a decisive choice between those alternatives. However, from a general point of view, dynamically hot environment of the Virgo does not favour the slow accretion of cold gas. IC 3653 is not a member of a subgroup including large galaxies, which can foster gaseous disc, so we believe that in this particular case the scenario of slow accretion is not applicable. At present, the sample of objects where disc-like substructures were searched either from images or integral field spectroscopy is still too small to draw statistical conclusions. However, it is quite probable that the progenitors of the dEs were disc galaxies (pre-dIrr or small spiral galaxies) and that they evolved due to feedback of the star formation and environmental effects. Present dIrr also experienced feedback but kept their gas, so it is unlikely that the feedback alone can remove the gas. Therefore environmental effects are probably the driver of the evolution of dEs, and the discovery of stellar discs in dEs is consistent with this hypothesis. A CKNO WLEDGMENTS We are very grateful to Alexei Moiseev for supporting the observations of IC 3653 at the 6-m telescope. Visits from PP in Russia and IC in France were supported through a CNRS grant. PhD of IC is supported by the INTAS Young Scientist Fellowship (04-83-3618). The dwarf galaxies investigation by IC, OS and VA is supported by the bilateral Flemish­Russian collaboration (project RFBR-0502-19805-MF a). We also appreciate support by IAU for attending the IAU Colloquium 198, where preliminary results of this paper were presented. We would like to thank our colleagues: Francois
C

Finally, we will review the different possible origins of the particular properties of this galaxy.

6.1 Comparison of two techniques for estimating stellar population characteristics In Table 4, we present comparison between SSP-equivalent age and metallicity obtained with the pixel fitting and inversion of bi-index grids for H ­Mgb, H ­ Fe , H ­[MgFe] using index measurements on observed spectra and best-fitting templates. One can notice particularly good agreement between the approaches. The ages derived from Lick indices appear to be slightly older, but the difference is not significant. The best agreement for both ages and metallicities is reached between pixel fitting and measurements of H and the combined [MgFe] index (Thomas et al. 2003). The internal precision of the parameters derived from pixel fitting is better than that from Lick indices by a factor of 3­4, depending on the indices used. This can be explained by the optimized usage of the information by the pixel-fitting procedure. Though it is difficult to assess the reliability of these small error bars, the relative variations of age and metallicity can be trusted, even when the signal-to-noise ratio is as low as 10 per pixel (with MPFS spectral resolution and wavelength coverage).

6.2 Properties and nature of IC 3653 We computed the position of IC 3653 on the Fundamental Plane (FP; Djorgovski & Davis 1987). `Vertical' deviation from FP (d = -8.666 + 0.314e + 1.14 log 0 - log Re ; Guzman, Lucey & Bower ´ 1993) is d = 0.2. This deviation places IC 3653 into the centre of the cloud, representing dE galaxies in De Rijcke et al. (2005; their fig. 2, left-hand panel), and exactly on the theoretical predictions by Chiosi & Carraro (2002) and Yoshii & Arimoto (1987), overplotted on the same figure. The mean age of the stellar population of IC 3653, t = 5 Gyr, coincides with the mean age of dE galaxies in Virgo: tmean = 5 Gyr in Geha et al. (2003); tmean = 5­7 Gyr in van Zee et al. (2004b). By contrast, the metallicity of the main body, Z = -0.1, is slightly higher than in these other galaxies: Z mean = -0.3 in Geha et al. (2003); Z mean = -0.4 in van Zee et al. (2004b). It may be understood as a consequence of the luminosity­metallicity relation if we remember that IC 3653 is more luminous than most of the galaxies in the samples of Geha et al. (2003) and van Zee et al. (2004b). We see that fundamental properties of IC 3653 do not differ from typical dE galaxies, though the effective radius is one of the smallest within the samples of Virgo dEs presented in Simien & Prugniel (2002), Geha et al. (2003) and van Zee et al. (2004a). We derived the B-band mass-to-light ratio of IC 3653 following the method by Richstone & Tremaine (1986) as M /LB = 8.0 ± 1.5(M /LB ) . The stellar mass-to-light ratio of the best-fitting PEGASE.HR model in the `P1' region is (M /LB ) = 3.7 ± 0.4 (M /LB ) . If we trust this simple dynamical model, it indicates that the dark matter content of IC 3653 is similar to recently found in De Rijcke et al. (2006) for the M31 dE companions, i.e. about 50 per cent. Recent theoretical studies based on N-body modelling of the evolution of disc galaxies within a Lambda cold dark matter ( CDM) cluster by Mastropietro et al. (2005) suggest that discs are never completely destroyed in cluster environment, even when the morphological transformation is quite significant. Our discovery of the faint stellar disc in IC 3653 comforts these results. Thus, one of the possible origins of IC 3653 is morphological transformation in

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046


Kinematics and stellar populations of IC 3653
Simien and Sven De Rijcke, and students: Mina Koleva and Nicolas Bavouzet for fruitful discussions. Special thanks to the Large Telescopes Time Allocation Committee or the Russian Academy of Sciences for providing observing time with MPFS.

1043

REFERENCES
Barazza F. D., Binggeli B., Jerjen H., 2002, A&A, 391, 823 Bertola F., Capaccioli M., 1975, ApJ, 200, 439 Bertola F., Galletta, 1978, ApJ, 226, L115 Binggeli B., Sandage A., Tammann G. A., 1985, AJ, 90, 1681 Binggeli B., Tarenghi M., Sandage A., 1990, A&A, 228, 42 Binney J., 1976, MNRAS, 177, 19 Cappellari M., Copin Y., 2003, MNRAS, 342, 345 Cappellari M., Emsellem E., 2004, PASP, 116, 138 Cardiel N., Gorgas J., Cenarro J., Gonzalez J. J., 1998, A&AS, 127, 597 ´ Chilingarian I., 2006, PhD thesis, Moscow State Univ./Univ. Claude Bernard Lyon-1 (astro-ph/0611893) Chiosi C., Carraro G., 2002, MNRAS, 335, 335 Conselice C., Gallagher J., Wyse R., 2001, ApJ, 559, 791 Conselice C., O'Neil K., Gallagher J., Wyse R., 2003, ApJ, 591, 167 C^ e P. et al., 2004, ApJS, 153, 223 ot ´ de Bruyne V., De Rijcke S., Dejonghe H., Zeilinger W. W., 2003, MNRAS, 349, 461 De Rijcke S., Dejonghe H., Zeilinger W. W., Hau G., 2001, ApJ, 559, L21 De Rijcke S., Dejonghe H., Zeilinger W. W., Hau G. K. T., 2004, A&A, 426, 53 De Rijcke S., Michielsen D., Dejonghe H., Zeilinger W. W., Hau G. K. T., 2005, A&A, 438, 491 De Rijcke S., Prugniel Ph., Simien F., Dejonghe H., 2006, MNRAS, 369, 1321 Djorgovski S., Davis M., 1987, ApJ, 313, 59 Duc P.-A., Mirabel I. F., 1999, in Barnes J. E., Sanders D. B., eds, Proc. IAU Symp. 186, Galaxy Interactions at High and Low Redshift. Kluwer, Dordrecht, p. 61 Ferrarese L. et al., 2006, ApJS, 164, 334 Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950 Gavazzi G., Boselli A., Donati A., Franzetti P., Scodeggio M., 2003, A&A, 400, 451 Geha M., Guhathakurta P., van der Marel R. P., 2002, AJ, 124, 3073 Geha M., Guhathakurta P., van der Marel R. P., 2003, AJ, 126, 1794 Geha M., Guhathakurta P., van der Marel R. P., 2005, AJ, 129, 2617 Guzman R., Lucey J. R., Bower R. G., 1993, MNRAS, 265, 731 ´ Held E., Mould J., 1994, AJ, 107, 1307 Illingworth G., 1977, ApJ, 218, L43 Jerjen H., Kalnajs A., Binggeli B., 2000, A&A, 358, 845 Le Borgne D., Rocca-Volmerange B., Prugniel P., Lancon A., Fioc M., ¸ Soubiran C., 2004, A&A, 425, 881 Lisker T., Grebel E., Binggeli B., 2006, AJ, 132, 497 Malin D. F., Carter D., 1983, ApJ, 274, 534 Marcolini A., Brighenti F., D'Ercole A., 2003, MNRAS, 345, 1329 Mastropietro C., Moore B., Mayer L., Debattista V., Piffaretti R., Stadel J., 2005, MNRAS, 364, 607 Michielsen D., De Rijcke S., Dejonghe H., Zeilinger W. W., 2003, ApJ, 597, L21 Michielsen D., De Rijcke S., Zeilinger W. W., Prugniel Ph., Dejonghe H., Roberts S., 2004, MNRAS, 311, 1293 Moiseev A. V., 2001, BSAO, 51, 140 Moore B., Lake G., Katz N., 1998, ApJ, 495, 139 Moultaka J., Pelat D., 2000, MNRAS, 314, 409 Moultaka J., Boisson C., Joly M., Pelat D., 2004, A&A, 420, 459 Ocvirk P., Pichon C., Lancon A., Thiebaut E., 2006a, MNRAS, 365, 46 ¸ ´ Ocvirk P., Pichon C., Lancon A., Thiebaut E., 2006b, MNRAS, 365, 74 ¸ ´ Paturel G., Petit C., Prugniel Ph., Theureau G., Rousseau J., Brouty M., Dubois P., Cambresy L., 2003, A&A, 412, 45 Pedraz S., Gorgas J., Cardiel N., Sanchez-Blazquez P., Guzman R., 2002, ´ ´ ´ MNRAS, 332, L59
C

Peng C., Ho L. C., Impey C. D., Rix H.-W., 2002, AJ, 124, 266 Prugniel P., Soubiran C., 2001, A&A, 369, 1048 Prugniel P., Soubiran C., 2004, preprint (astro-ph/0409214) Prugniel P., Chilingarian I., Sil'chenko O., Afanasiev V., 2005, in Jerjen H., Binggelli B., eds, Proc. IAU Coll. 198, Near-Field Cosmology with Dwarf Elliptical Galaxies. Cambridge Univ. Press, Cambridge, p. 73 Richstone D., Tremaine S., 1986, AJ, 92, 72 Sadler E. M., Gerhard O. E., 1985, MNRAS, 214, 177 Sil'chenko O., Shapovalova A., 1989, SoSAO, 60, 44 Simien F., Prugniel P., 2002, A&A, 384, 371 Sirianni M. et al., 2005, PASP, 117, 1049 Thomas D., Maraston C., Bender R., 2003, MNRAS, 339, 897 Thomas D., Brimioulle F., Bender R., Hopp U., Greggio L., Maraston C., Saglia R. P., 2006, A&A, 445, L19 Trager S. C., Worthey G., Faber S. M., Burstein D., Gonzalez J. J., 1998, ApJS, 116, 1 van den Bergh S., 1994, ApJ, 428, 617 van der Marel R., Franx M., 1993, ApJ, 407, 525 van Zee L., Skillman E. D., Haynes M. P., 2004a, AJ, 128, 121 van Zee L., Barton E. J., Skillman E. D., 2004b, AJ, 128, 2797 Vazdekis A., Arimoto N., 1999, ApJ, 525, 144 Weilbacher P. M., Duc P.-A., Fritze von Alvensleben U., Martin P., Fricke K. J., 2000, A&A, 358, 819 Wheeler J. C., Sneden C., Truran J. W. Jr, 1989, ARA&A, 27, 279 Worthey G., 1994, ApJS, 95, 107 Worthey G., Faber S., Gonzalez J., Burstein D., 1994, ApJS, 94, 687 ´ Yoshii Y., Arimoto N., 1987, A&A, 188, 13

APPENDIX A: VALID A TION AND ERR OR AN AL YSIS In this appendix, we address questions concerning the error analysis, stability of the solutions, and possible biases for relatively old stellar population (about 5 Gyr). A full description of these aspects extended to any instrument and much wider range of parameters will be described in details in the forthcoming paper (Koleva et al., in preparation). Here we give only essential error analysis required for validation of the results presented in this paper and in the forthcoming papers based on MPFS data for dwarf galaxies. A1 Error analysis Error estimates on the stellar population parameters is complicated because of the internal degeneracies and coupling, and the questions of unicity of the solution and local minima are always critical. A complete and detailed description of one of the possible approaches to locate the alternate solutions using for a allied inversion technique is given in Moultaka & Pelat (2000). We performed some Monte Carlo simulations (about 10 000 per spectrum for the 3-points binning) to demonstrate the consistency between the uncertainties on the parameters reported by the minimization procedure and the real error distributions. These simulations have demonstrated that in case of our data set, where there is neither significant model mismatch due to element abundance ratios, nor strong mistake with subtraction of additive terms (diffuse light and night sky), the uncertainties found from the Monte Carlo simulations using scattering of solutions in the multidimensional parameter space coincide with values reported by the minimization procedure being multiplied by 2 values. Deviation of 2 from 1 might be caused either by poor quality of the fit (model mismatch) or by wrong estimations of absolute flux uncertainties in the input data. We inspected the residuals of the fit visually and by averaging them using wide smoothing window (150 pixels) and found no significant deviations between model and observations. Thus we

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046


1044

I. V. Chilingarian et al.
rors, suggesting that scanning of the 2 hyperspace on v is not necessary. In Fig. A1 (upper line), we present the maps of 2 for the 3points binning on the `t­Z' plane. One can see elongated shapes of the minima, corresponding to well known age­metallicity degeneracy. Three plots on the bottom of Fig. A1 represent slices of the 3D 2 space scan (t ­Z ­ ) for the `P1' bin. One can notice that the width of the minimum on `t­Z' plane has decreased due to a correlation between metallicity and velocity dispersion, that is clearly seen on the `Z ­ ' slice. This degeneracy between velocity dispersion and metallicity can be understood as follows: an overestimate of the metallicity in the template increases depth of the absorption lines and can be compensated by higher velocity dispersion. In principle, an important mismatch of the metallicity when analysing the kinematics with usual methods may bias the velocity dispersion and even produce patterns in velocity dispersion profiles, whereas it is in reality a metallicity gradient. A mismatch of 0.1 dex results in a bias of 15 per cent on the velocity dispersion. A2 Stability of solutions Stability of solutions is a crucial point for every method dealing with multiparametric non-linear minimization. We studied the stability with respect to initial guess, wavelength range being used, and degree of the multiplicative polynomial continuum. A2.1 Initial guess We have made several dozens of experiments with different initial guesses in order to inspect the stability of convergence. We found a correct convergence for a wide range of guesses of age, metallicity and velocity dispersions, but for radial velocity the initial guess must be within twice the velocity dispersion from the solution.

conclude that in some cases our estimations of the uncertainty on the absolute fluxes based on the photon statistics may be wrong by 15 per cent, resulting in values of 2 between 0.7 and 1.3. To estimate the errors more accurately, locate possible alternate solutions and search for degeneracies between kinematical and stellar population parameters we study the distribution of 2 in the age­metallicity­velocity dispersion space. The procedure we followed includes the following steps. (i) We chose a grid of values for age, metallicity and velocity dispersion, that was supposed to cover a reasonable region of the parameter space where we could expect to have solutions. In our case the grid was defined as 2 < t < 14 Gyr with a step of 200 Myr, -0.45 < [Fe/H] < 0.40 with a step of 0.01 dex, 30 < < 100 km s-1 with a step of 0.5 km s-1 . (ii) At every node of t­Z grid, we ran the pixel-fitting procedure in order to determine multiplicative polynomial continuum, and to have the best fit for a given SSP. (iii) Later 2 was computed on the grid of values of by fixing all other components of the solution that had been found in the previous step. The reason to scan the dimension with fixed multiplicative polynomial is to save computer time by reducing the number of free parameters, and it is justified by the strong decoupling between the multiplicative polynomial (affecting only low frequencies) and the velocity dispersion (affecting high frequencies). In other words using our procedure we compute a slice of the full 2 by the hypersurface defined as a set of minimal 2 values for multiplicative continuum terms and radial velocity values, and then reproject it on to `t­Z' and `t­Z ­ ' hyperplanes. The result contains two arrays: 2D age­metallicity and 3D age­metallicity­ velocity dispersion. The values of line-of-sight radial velocity obtained during the fitting procedure on the `t­Z' grid are constant within the er-

Figure A1. Two-dimensional maps of 2 distributions for the 3-points binning (upper row), and slices of the three-dimensional 2 distribution (t­Z­ ) for the `P1' bin.
C

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046


Kinematics and stellar populations of IC 3653

1045

Figure A2. Stability of the parameters with respect to the order of the multiplicative polynomial continuum. Different plotting symbols correspond to three bins of the 3-points binning.

Figure A3. Effects of additive terms on the results. Abscissa represents the contribution of the constant level at = 5000 å. Different plotting symbols correspond to each bin of the 3-points binning.

A2.2 Wavelength range We ran two series of experiments: one with > 4700 å, and another one with the full wavelength range, but regions of Balmer lines (H and H ) masked. The first experiment checks the sensitivity to the wavelength range (several strong metallic features are excluded) and the second one removes the best age estimators (Worthey 1994; Vazdekis & Arimoto 1999), so that instability on age may be expected. The first experiment does not reveal any bias, but only an expectable increase of the uncertainties. The second experiment does not neither reveal instability, but though the errors on metallicity and velocity dispersion are not strongly affected, the error on age is multiplied by 3 (therefore the precision become comparable to Lick indices). It is remarkable that even without Balmer lines it is possible to give estimations of age of the stellar population. A2.3 Order of multiplicative polynomial continuum We also explored the stability of the method with respect to the order of the multiplicative polynomial continuum. The results are shown in Fig. A2. One can see that for n > 5, there is significant changes neither of the estimations of kinematical and stellar population parameters, nor of 2 value. Therefore we chose n = 5 for all our data analysis. A3 Possible biases There are several possible sources of systematic errors on the parameters: (i) additive systematics of the flux calibration due to underor oversubtraction of the night sky, (ii) imperfections of the models, one of the most important among those is non-solar abundance ratio. Since IC 3653 exhibits exactly solar [Mg/Fe] (see Fig. 2) everywhere in the field of view, we do not study possible consequences of abundance ratio mismatch in this paper. The accurate subtraction of the night-sky emission is a challenging step of the data reduction for low-surface brightness objects. Basically, night-sky emission consists of continuum emission, that might include scattered solar light as well, and several bright emission lines. Under- or oversubtraction of night sky brings additive component resulting in changing the relative depths of specC

Table A1. Stability of the solutions for the 3-points binning with respect to the wavelength range. First lines for every parameter correspond to > 4700 å, second ones to full range with Balmer lines excluded. P1 v , km s
-1

P2 603.4 ± 1.4 603.1 ± 1.7 67.3 ± 2.2 65.3 ± 2.6 4.728 ± 0.289 4.448 ± 0.403 -0.14 ± 0.02 -0.13 ± 0.02

P3 603.8 ± 3.0 603.7 ± 3.4 52.1 ± 5.0 52.1 ± 5.7 4.629 ± 0.694 4.238 ± 0.930 -0.15 ± 0.05 -0.15 ± 0.05

601.8 ± 1.0 600.9 ± 1.0 70.9 ± 1.6 71.8 ± 1.6 4.855 ± 0.218 4.714 ± 0.235 0.01 ± 0.02 0.03 ± 0.01

, km s t, Gyr Z, dex

-1

tral features, and therefore affects particularly the estimates of the metallicity. We have conducted two series of experiments: (i) adding a constant term or (ii) the spectrum itself smoothed by a 300-pixel window to simulate the diffuse light in the spectrograph. In every series the fraction of the additive term was between -20 and +50 per cent to model over- and undersubtraction. The results for the constant term as a fraction of flux at 5000 å are shown in Fig. A3. One may notice that 2 reaches minimum on slightly negative (oversubtraction) values of the additive term because we did not change flux uncertainties during our experiments. The remarkable result is the stability of age estimations on a wide range of additive components (-25 to 15 per cent). This is an important advantage of the pixelfitting technique over Lick indices, more sensitive to additive continuum. Metallicity and velocity dispersion exhibit the expected behaviour: growth of and fall of Z. Indeed within a range of contribution between -5 and 5 per cent changes are quite small (8 per cent for , and 0.1 dex for Z) though significant. In order to test the consequences of bad sky subtraction, we made two additional experiments: we tried to fit the data, where the sky spectrum was represented by a low-order polynomial continuum and with no sky subtraction at all. We excluded four regions of the spectrum containing bright emission lines: Hg I = 4358 å, 5461 å, [N I] = 5199 å and [O I] = 5577 å. The experiments

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046


1046

I. V. Chilingarian et al.
Table A2. Determination of the kinematical and stellar population parameters for a case of polynomial night-sky model (first line for every parameter) and no sky subtraction (second line for every parameter). P1 v , km s
-1

were made for the 3-points binning of the data, demonstrating the effects for high, intermediate, and low surface brightness (see Table 2). Basically, we found no significant difference for the `P1' and `P2' bins between the parameters for the correct sky subtraction and subtraction of the low-order polynomial model of the sky (see Table A2). `P3' bin gives younger age and higher metallicity, but the estimations are in agreement with the normal sky subtraction within 2 However, as expected, when sky is not subtracted at all we find sizeable bias on age, metallicity, and velocity dispersion estimations for the `P2' bin, and even stronger effect for `P3'. Due to additive continuum metallicities are found to be lower, ages older, and velocity dispersions higher than expected. These experiments demonstrate that for the surface brightness down to B = 20 mag arcsec-2 features of the night-sky spectrum do not affect the results of the pixel-fitting procedure, and very rough sky subtraction is sufficient to obtain the realistic estimations of kinematical and stellar population parameters.

P2 606.0 ± 1.5 605.9 ± 1.6 64.9 ± 2.4 88.3 ± 2.2 4.547 ± 0.310 7.203 ± 0.390 -0.10 ± 0.02 -0.50 ± 0.02

P3 609.4 ± 3.2 610.3 ± 3.7 54.4 ± 5.2 105.5 ± 4.6 3.956 ± 0.731 12.982 ± 1.641 -0.06 ± 0.04 -0.94 ± 0.03

604.3 ± 1.0 604.4 ± 1.0 71.5 ± 1.5 80.5 ± 1.5 4.868 ± 0.210 4.972 ± 0.185 0.04 ± 0.01 -0.08 ± 0.01

, km s t, Gyr Z, dex

-1

A This paper has been typeset from a TEX/L TEX file prepared by the author.

C

2007 The Authors. Journal compilation

C

2007 RAS, MNRAS 376, 1033­1046