Issue |
A&A
Volume 497, Number 1, April I 2009
|
|
---|---|---|
Page(s) | 41 - 51 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/200810506 | |
Published online | 18 February 2009 |
Early-type galaxies in the near-infrared: 1.5-2.4
m spectroscopy
,![[*]](/icons/foot_motif.gif)
M. Cesetti1,2 - V. D. Ivanov1 - L. Morelli2 - A. Pizzella2 - L. Buson3 - E. M. Corsini2 - E. Dalla Bontà2 - M. Stiavelli4 - F. Bertola2
1 - European Southern Observatory, Ave. Alonso de Còrdova 3107, Casilla 19, Santiago 19001, Chile
2 - Dipartimento di Astronomia, Università di Padova, vicolo dell'Osservatorio 3, 35122 Padova, Italy
3 - Istituto Nazionale di Astrofisica, Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 35122 Padova, Italy
4 - Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA
Received 2 July 2008 / Accepted 21 August 2008
Abstract
Context. Near-infrared (hereafter NIR) data may provide complementary information to the traditional optical population synthesis analysis of unresolved stellar populations because the spectral energy distribution of the galaxies in the 1-2.5 m range is dominated by different types of stars than at optical wavelengths. Furthermore, NIR data are less affected by dust absorption and hence could constrain the stellar populations in dust-obscured galaxies.
Aims. We attempt to develop observational constraints on the stellar populations of unresolved stellar systems in the NIR.
Methods. To achieve this goal, we require benchmark NIR spectra of ``simple'' early-type galaxies, to test and calibrate the outputs of population synthesis models. We obtained low-resolution (
)
long-slit spectra between 1.5 and 2.4
m for 14 nearby early-type galaxies using SofI at the ESO 3.5-m New Technology Telescope and higher resolution (
)
long-slit spectra, centered on Mg I at
1.51
m for a heterogeneous sample of 5 nearby galaxies observed with ISAAC at Antu, one of the 8.2-m ESO Very Large Telescope.
Results. We defined spectral indices corresponding to CO, Na I, Ca I and Mg I absorption lines and measured their strengths in the galaxies spectrum. We defined a new global NIR metallicity index, suitable for abundance measurements in low-resolution spectra. Finally, we created an average-combined NIR spectrum of an early-type galaxy, from a homogeneous subset of our data set.
Conclusions. The NIR galaxy spectra are similar to each other and some of the line index strengths correlate with both the iron abundance [Fe/H] and optical metal line indices of the galaxies. The data suggest that the NIR metal features, in combination with a hydrogen absorption feature, may break the age-metallicity degeneracy in a similar way the Mg and Fe features in the rest-frame optical wavelength range.
Key words: infrared: galaxies - Galaxy: abundances - galaxies: elliptical and lenticular, cD - galaxies: spiral - galaxies: stellar content
1 Introduction
State-of-the-art space-based instrumentation can resolve the stellar populations only of galaxies in the Local Group and in more distant objects we can usually detect only the tip of the red giant branch stars, which are rarely sufficient for robust population analyses. We must therefore infer the stellar population properties indirectly from the integrated galaxy light. The complex mix of stellar populations present in most galaxies imply that it is usually possible to study only the properties of the most recent generation of stars. However, properties such as stellar kinematics and present-day metal content are affected significantly by the overall star formation history and, due to the well known age-metallicity degeneracy (O'Connell 1986; Worthey 1994; Faber 1973), they often lead to non-unique solutions for the stellar populations. In this respect, the Lick/IDS system of indices pioneered by Burstein et al. (1984) and Faber et al. (1985), and developed further by Trager et al. (1998,2000,2005), has been particularly successful in interpreting the integrated optical light of galaxies.
However, new constraints are necessary to interpret more complicated systems and one possibility is to widen the spectral range towards the near-infrared (NIR) because light in different wavebands is dominated by different populations of stars. The NIR passbands are dominated by light from older and redder stars and therefore offer us the possibility to study other stellar populations than possible with optical spectra alone. In addition, abundance determinations with optical spectroscopy are impossible for heavily reddened evolved stellar populations such as dusty spheroids or some bulge globular clusters hidden by dust (e.g. Stephens & Frogel 2004). NIR spectroscopic observations could overcome these problems because the extinction in the K-band is only one-tenth that in the V-band.
Most of the previous work at NIR was focused on either active galactic nuclei (AGNs) or objects with very strong star formation, including recent surveys of ultra-luminous infrared galaxies (Goldader et al. 1995; Burston et al. 2001; Murphy et al. 2001,1999), luminous infrared galaxies (Reunanen et al. 2007; Goldader et al. 1997), starbursts (Engelbracht 1997; Coziol et al. 2001), Seyfert galaxies (Sosa-Brito et al. 2001; Ivanov et al. 2000; Boisson et al. 2002; Reunanen et al. 2002), LINERs (Alonso-Herrero et al. 2000; Sosa-Brito et al. 2001; Larkin et al. 1998), and interacting galaxies (Vanzi et al. 1998; Vanzi & Rieke 1997).
Relatively few NIR spectroscopic observations exist for ``normal''
galaxies. Only Mannucci et al. (2001) provided low-resolution (
)
template spectra for galaxies of different Hubble types, including
some giant ellipticals. Such data, together with the corresponding
analysis is a necessary first step towards developing a system of
spectral diagnostics in the NIR because well-understood galaxies with
relatively simple star-forming histories allow us to tune the NIR
population synthesis models. Silva et al. (2008, hereafter S08)
studied the stellar populations of eleven early-type galaxies in the
nearby Fornax cluster by means of K-band spectroscopy.
A few prominent NIR features were first studied by Origlia et al. (1993), who demonstrated that they represent a superb set of indicators for constraining the average spectral type and luminosity class of cool, evolved stars. This conclusion was later confirmed by Förster Schreiber (2000) and Ivanov et al. (2004). Furthermore, the same NIR features appear to be promising abundance indicators (Frogel et al. 2001; Stephens & Frogel 2004).
Here we describe two new data sets of high quality NIR spectra of
ellipticals/spirals designed to provide a benchmark for future NIR
studies of unresolved galaxies: (i) low-resolution spectra
covering the range from 1.5 to 2.4 m that include strong
features such as CO, Na I and Ca I that are traditionally
studied, and (ii) moderate resolution spectra around the Mg
I absorption feature at 1.51
m. This is the second strongest
Mg feature in the H- and K-band atmospheric windows (after the
feature at 1.71
m) and at zero redshift is located in a region
of poor atmospheric transmission, making it difficult to observe in
stars. However, at high redshift the observable rest-frame spectrum of
external galaxies increases to wavelengths at which the atmosphere is
more transparent (Ivanov 2001), close to the 1.71
m Mg I, which itself becomes affected by the red edge of the H-band
atmospheric window.
We are only few years from the launch of the James Webb Space Telescope (Gardner et al. 2006), a space-based infrared telescope that will have unprecedented capabilities. Therefore, in the near future the application of NIR spectroscopy to the study of galaxy properties will be limited not by lack of data but by our understanding of spectral features at these wavelengths. Improving the characterization of the NIR indices is a timely step in this direction.
The paper is organized as follows. The sample selection is discussed in Sect. 2. The NIR spectroscopic observations and data reduction are described in Sect. 3. The definition of the new NIR spectral indices and their measurements are given in Sect. 4. Results are discussed in Sect. 5 and summarized in Sect. 6.
2 The sample
Our main sample consists of 14 nearby bright (
14.1 mag),
undisturbed spheroids. They were selected to cover a wide range of
luminosity and velocity dispersion and to have Lick/IDS line-strength
indices available, from the literature. Six of our galaxies are giant
ellipticals (NGC 4472, NGC 4621,
NGC 4649, NGC 4697, NGC 6909,
IC 4296), five are intermediate-size ellipticals
(NGC 4478, NGC 4564, NGC 3818,
NGC 5077, NGC 5576), one is a bright dwarf
elliptical (NGC 3641), and two are bulges of disk galaxies
(NGC 4281, NGC 4594; our spectra are dominated by
their spheroids). The Mg2 index spans the range between 0.21 and
0.34 mag (Bender et al. 1993), in agreement with other Mg2 literature
values, and the central velocity dispersions are between about 130 and
370 km s-1. The galaxy properties are summarized in
Table 1.
Four galaxies in the sample exhibit evidence for weak nuclear activity: NGC 4594 and NGC 5077 were identified as LINERs by Heckman (1980); NGC 4472 is a weak Seyfert 2 (Ho et al. 1997); and IC 4296 shows radio and X-ray emission (Rinn et al. 2005). Finally, NGC 4649 is a close companion of the giant spiral NGC 4647 and it may have undergone some tidal stripping (Das & Jog 1999).
The low resolution used for the observation of our sample is
insufficient to characterize the 1.51 m Mg I feature, so we
selected higher resolution spectra from the ESO Science Archive,
covering the region around the 1.51
m, for a heterogeneous set
of galaxies. Although these observations were acquired for other
purposes, they provide a first glimpse into the behavior of this
feature. The properties of these galaxies are also described in
Table 1.
Table 1: Basic properties of the sample galaxies.
3 Observation and data reduction
3.1 SofI low-resolution spectra
Long-slit (
)
NIR spectroscopy of 14 galaxies was
performed with the SofI (Son of ISAAC) infrared spectrometer
(Moorwood et al. 1998b) at the Nasmyth A focus of the ESO 3.5-m New Technology
Telescope (NTT) at the European Southern Observatory (ESO) in la Silla
(Chile). The detector is a
pixels Hawaii HgCdTe array
with a 18.5-
m pixel size. It is read out in four quadrants and
the average quantum efficiency is
65 %. We used the low
resolution grism (
)
which provided a spectral dispersion of
10.22 Å pixel-1 and a coverage of between 1.53
m and
2.52
m.
Table 2: Log of the spectroscopic observations.
The smaller galaxies were observed in nodding mode, placing them at two different positions along the slit, separated by 150''. In the case of larger galaxies, the telescope pointing was alternated between the galaxy and a nearby clear sky region to obtain a good sky sampling. We randomly dithered the telescope within a 10'' box along the slit about each position. This dithering, also known as jitter (see e.g. SofI handbook or Devillard 1999), helped us to eliminate bad pixel events, and to improve both the pixel sampling and flat-field correction, because, on each image, the object was placed on slightly different positions of the array.
Immediately before or after each galaxy spectrum we observed a B star to measure the atmospheric transmission and the instrument response, always in nodding mode and with a similar airmass to that of the corresponding galaxy. Finally, to calibrate the spatial distortions we also obtained a sequence of stellar spectra placed at 19 equidistant positions along the slit.
The details of the spectroscopic observations are given in Table 2.
The spectra were reduced with the Image Reduction and Analysis
Facility (IRAF). The sky emission was removed by either subtracting the
succeeding sky image from the corresponding object image or by
subtracting, from each object image, the average of the preceding and
the succeeding sky images, respectively, for the two modes of
observation described above. The data were flat-fielded with screen
flat spectra. Geometric distortion correction was applied along both
the spatial and dispersion axes by using an arc frame and the grid of
stellar spectra mentioned above, respectively. Simultaneously, the
frames were wavelength calibrated. We applied the dispersion
corrections, and removed the hot pixels and cosmic ray hits with the
algorithm described in Pych (2004). Then, we aligned the individual
2-dimensional spectra along the spatial direction (i.e. along
the slit), to combine all individual spectra and extract a
1-dimensional spectrum of the galaxy for the apertures listed in
Table 2.
![]() |
Figure 1:
The telluric absorption correction: the wavelength
calibrated 1-dimensional galaxy spectrum (a) is divided by the
wavelength calibrated 1-dimensional observed spectrum of the
standard star (b) and the result is multiplied by the reference of
the same spectral type of the standard star (c) from the library
of Pickles (1998). The product (d) is the corrected galaxy
spectrum. The large scatter at |
Open with DEXTER |
Finally, we removed the telluric absorption. We first created a
sensitivity function by multiplying the observed spectrum of the
standard by the intrinsic spectrum of a star with the corresponding
(or the closest available) spectral type from Pickles (1998). Some
residual emission was observed, indicating that the spectral types of
the standards may be inaccurate. We fitted these artificial features
with a Gaussian, using the IRAF task SPLOT, and subtracted them from
the sensitivity function. This prompted us to disregard the bluest
part of the spectra (
m) in further analysis
because the standards were rich in hydrogen lines that modified the
appearance of the galaxy spectra. We, next, divided the galaxy spectra
by this sensitivity function to restore the overall continuum shape of
the galaxy spectrum. This process, divided into steps, is illustrated
in Fig. 1 and more detailed descriptions can be
found in Hanson et al. (1996) and Maiolino et al. (1996). The final reduced spectra
were corrected to zero redshift (using values of redshift given in Table
1) and are shown in Figs. 2
and 3.
![]() |
Figure 2:
H-band spectra of the sample galaxies, continuum-divided
(i.e. all continua were normalized to 1) and shifted vertically
for display purposes by adding increasing shifts of 0.25. The
position of the main spectral features is shown with dotted lines.
The passbands of the Mg I (1.71 |
Open with DEXTER |
![]() |
Figure 3:
K-band spectra of the sample galaxies, continuum-divided
(i.e. all continua were normalized to 1) and shifted vertically
for display purposes by adding increasing shifts of 0.25. The
positions of the main spectral features are indicated by dotted
lines. The passbands of the Na I (2.21 |
Open with DEXTER |
3.2 ISAAC high-resolution spectra
Five spectra of galaxies were obtained with ISAAC - a NIR imaging
spectrometer (Moorwood et al. 1998a) at Nasmyth B focus of Antu, one of the
8.2-m ESO VLT unit telescopes of ESO Very Large Telescope at Paranal
(Chile). It is equipped with a
Hawaii Rockwell HgCdTe
array with a spatial scale of 0.146'' pixel-1. The long-slit
(
)
medium (
)
resolution spectroscopy
mode was used. The wavelength coverage spans the range from
1.48
m to
1.62
m and the central wavelength
chosen depended on the object redshift. A summary of observations is
given in Table 2.
In a similar way to the SofI observations, telescope pointings ware alternated between the galaxy and a nearby clear region of sky, or if the galaxy dimensions were sufficiently small, then by nodding among two positions along the slit separated by about 20''. No jittering was used. Standards were observed in nodding mode immediately before or after the main galaxy exposure, at a similar airmass. A sequence of stellar spectra placed at equidistant positions along the slit was obtained to calibrate the spacial distortions, as in the case of the SofI data.
The reduction steps were carried out using IRAF in an identical way to that for the main data set, albeit with some complication due to the lower signal-to-noise ratio of the data and longer exposure times. First, the individual images were cleaned from cosmic rays using SPEC_COSM, a cosmic-ray rejection script for long-slit spectra developed by van Dokkum (2001). Second, bad pixels were masked out before the combination because the small number of images limited the ability of the rejection algorithms to remove the bad pixels during combination of the individual images. Finally, the highly variable sky background produced residuals in the sky-subtracted images, which were removed with the IRAF task BACKGROUND before the geometric distortion correction and the wavelength calibration.
4 Definition of the NIR spectral indices
The NIR spectral indices were usually defined with a central bandpass covering the spectral feature of interest, and two other bandpasses at the red and blue sides used to trace a local continuum level by a linear fit to the mean values in both bands, as in Ali et al. (1995), Frogel et al. (2001), and Ivanov et al. (2004). When this was not possible, for example because the feature was at the end of the atmospheric window, then either one continuum band or a fit to the overall shape of the continuum was used, as in Kleinmann & Hall (1986).
Most of the earlier definitions were developed for stars, where the
spectral features are typically only a few km s-1 wide; in
galaxies, however, the intrinsic velocity dispersion broadens them to
a few hundred km s-1. We therefore needed to define new spectral
indices with broader passbands (Table 3,
Figs. 2 and 3). The continuum
passbands were placed in wavelength regions clear of absorption
features, and not (or only weakly) affected by the telluric
absorption. Furthermore, we attempted to find as wide continuum
passbands as possible, to increase the signal-to-noise ratio of the
index measurements. The passband of each spectral feature was defined
to encompass the maximum extent of the feature itself and avoid, if
possible, contamination by other features. For instance, the passband
of our CO (2.30 m) index is wide but still only includes the
12CO feature and no contribution from the 13CO band head.
Table 3: NIR spectral indices.
We measured the equivalent width (EW) of the lines with respect to a
local continuum obtained by fitting a straight line to the ``clear''
parts of the spectrum:
where





Table 4: Equivalent widths of the newly defined NIR spectral indices for the sample galaxies.
Our definitions have intentionally wide bandpasses to ensure that they
are quite insensitive to the varying velocity broadening of spectral
features in the galaxy spectra. To verify this, we measured the
indices for broadened spectra of M and K type giants. We chose stars
in this spectral range because they are the most representative of
early-type galaxy integrated light spectra. (These spectra are a
reasonable match to the observed galaxy spectra.) We computed for them
the correction factor
for various line-of-sight velocity
distribution, as in S08:
![]() |
(2) |
The EW of the feature was measured from the stars convolved with a Gaussian of width


![]() |
Figure 4: Correction factors of the NIR spectral index for different values of velocity dispersion. They were derived by comparing the EW measured in the spectra of the labeled stars before and after broadening them to a particular velocity dispersion. The spectral types of the stars are in a range between KIII and MIII with a spread of [Fe/H] between -0.35 and +0.41. In each panel, the solid line is the 2nd order polynomial fitting the data. |
Open with DEXTER |
5 Discussion and conclusions
5.1 General appearance of the spectra
The spectra of our sample galaxies appear qualitatively similar in
most of the NIR features except for a large spread in the Na I
values. This similarity is unsurprising because we selected mostly
giant ellipticals and spheroids, which are all relatively metal rich,
have no significant recent star formation, and only weak nuclear
activity, if any. This conclusion agrees with Mannucci et al. (2001), who
demonstrated that galaxies of the same Hubble type have almost
identical NIR spectra. On the other hand, the large spread in NaI,
with respect to the observational errors, appears to be real, and
suggests that a variety of star formation and enrichment histories
exist among the galaxies in our sample. Two of the three galaxies with
systematically weaker NaI exhibit strong H,
indicating that the
weak NaI might be related to the presence of younger stellar
populations (the third galaxy lacks optical spectroscopy). The
strongest features are the CO absorption bands in both the H and
K-band atmospheric windows. They originate in K and M stars, as can
be seen from the libraries of stellar spectra described by
Lancon & Rocca-Volmerange (1992), Origlia et al. (1993), Dallier et al. (1996), Meyer et al. (1998),
Wallace & Hinkle (1997), and Förster Schreiber (2000). A number of weaker metal absorption
features are visible as well: Si I at 1.589
m, Mg I at
1.711
m, Na I at 2.206 and 2.209
m, and Ca I at
2.261, 2.263, and 2.266
m. These are also present mostly in
cool stars (i.e. Kleinmann & Hall 1986; Wallace & Hinkle 1997; Förster Schreiber 2000; Meyer et al. 1998; Origlia et al. 1993). As mentioned
above, we only consider the redder features with
1.65
m.
Unsurprisingly, our galaxies show no Fe II at 1.644 m or
H2 lines line emission, usually associated with supernova activity
and supernova related shocks, respectively. There is also no
Br
emission.
![]() |
Figure 5: The equivalent widths of the NIR spectral indices, iron abundance [Fe/H], and metal abundance [Z/H] in the sample galaxies as a function of the central velocity dispersion. The labels correspond to the row numbers in Table 1. The open and filled circles refer to quiescent and active galaxies, respectively. The small and large circles represent galaxies with an age of between 5 and 10 Gyr and >10 Gyr, respectively. Galaxies with no known ages were assumed to be old. The dashed lines corresponds to the relations found by S08. The solid line in each panel represents the best fit linear regression ( y = ax + b) for all the data points except for #4, #12 and #14. The Pearson correlation coefficient (r) and the results of the linear fit are given. Typical error bars are shown. |
Open with DEXTER |
5.2 Behavior of the NIR spectral indices
Evolutionary synthesis models (i.e. Worthey 1994) are the usual
method to interpret the behavior of spectral features in galaxies, but
the simple comparison of these features with the galaxy parameters can
be informative too. For example, the optical indices, Mg2 ,
Fe
,
and
show tight correlations with the central velocity
dispersion (e.g. Mehlert et al. 2003; Morelli et al. 2008; Terlevich et al. 1981; Bernardi et al. 1998), suggesting that
the chemical and dynamical evolution of ellipticals are intertwined.
The EW of our NIR spectral indices, the iron abundance, and
metallicity are plotted versus the central velocity dispersion in
Fig. 5. The loci of the S08 data (see their Fig. 13)
are shown with dashed lines. Among the galaxies, NGC 5077
(#12) and NGC 6909 (#14) possess very weak Na I,
Ca I lines, and their relative errors are significant. The
correlations between the NIR indices and velocity dispersion for the
sample galaxies are plotted in Fig. 5; the fit do not
consider NGC 4281 (#4), NGC 5077 (#12), and
NGC 6909 (#14) because of the reasons explained further on.
The correlation with CO shows a similar slope to that found by S08 but
with a different zero-point. The difference in the index definition
and especially the position of the continuum bandpass can produce
systematic variation of the EW values. This effect seems to be more
evident in the CO band, whose equivalent width is measured using only
an extrapolated continuum blue-wards of the feature. The agreement
with S08 is closer for the Na I versus
relation,
although our data show larger scatter. Without derive a general
conclusion, we found relatively high NaI values compared with stars,
similar was also found by S08. The iron abundance [Fe/H] and total
abundance [Z/H] were only weakly correlated. Active galaxies were
clustered towards the high end of the velocity dispersion
distribution, which is understandable because their velocity
dispersion measurement may be affected by the presence of a black hole
and active nucleus. We note that our highest velocity dispersion
galaxy is a type 2 Seyfert and the CO band may suffer some dilution,
as discussed by Ivanov et al. (2000).
![]() |
Figure 6: The equivalent widths of the NIR indices plotted against each other. The symbols are the same as in Fig. 5. The dashed lines correspond to the relations found by S08. The dotted line corresponds to the relation by Davidge et al. (2008) for giant stars in the solar neighborhood. All of these relations are plotted in their observed range. The solid line represents the linear regression ( y = ax + b) through all data points except for #4, #12, and #14. The Pearson correlation coefficient (r) and the results of the linear fit are given. Typical error bars are shown. |
Open with DEXTER |
The EWs of our NIR spectral indices are plotted versus each other in
Fig. 6 and compared with the relations for cluster
stars and solar neighborhood giants by S08 and Davidge et al. (2008),
respectively. Three galaxies - NGC 4281 (#4),
NGC 5077 (#12) and NGC 6909 (#14) - are above
the relation of S08 for cluster stars, which suggests that they may
have different chemical enrichment histories with respect to the
remainder of the sample. Unfortunately the literature lacks data about
these objects and although NGC 6909 is of the lowest
velocity dispersion in our sample, these galaxies do not appear to be
extremely distinctive in any respect, including in the terms of
optical Mg2 versus
diagram. Further investigation of
these galaxies is necessary. The remaining galaxies populate a locus
that follows similar trends to the galaxies of S08. In the Ca I
versus CO plot we see a correlation with significant scatter, and with
a different slope than in S08. Finally, the Mg I versus CO plot
is dominated by scatter. Galaxies with traces of nuclear activity,
evident in other wavelength ranges than the NIR, do not appear to
differ from the remainder of the sample, and show no emission lines.
Therefore, the contribution of their AGNs to their detected spectrum
is negligible compared with that of the remainder of the galaxy.
![]() |
Figure 7: The equivalent widths of the NIR indices plotted against the Mg2 indices measured by Bender et al. (1993). The symbols are the same as in Fig. 5. The solid line represents the linear regression ( y = ax + b) through all data points except for #4, #12, and #14. The Pearson correlation coefficient (r) and the results of the linear fit are given. Typical error bars are shown. |
Open with DEXTER |
![]() |
Figure 8:
a- d) The equivalent widths of the NaI, CaI, CO and
H |
Open with DEXTER |
![]() |
Figure 9:
The equivalent widths of the newly defined |
Open with DEXTER |
Figure 7 shows the EW of the NIR spectral indices
versus the Mg2 measurements from Bender et al. (1993). Na I, and to a
lesser extent Ca I and CO, show trends with respect to Mg2,
with NGC 4281 (#4), NGC 5077 (#12), and
NGC 6909 (#14) standing out. The correlation of Na I
with both
and Mg2 suggests that the Na I feature is
dominated by the stellar photosphere rather than the interstellar
medium.
Given that [Fe/H] measurements are unavailable for the entire sample
and Mg2 is unrepresentative of the total chemical abundance of a
galaxy because of the varying abundance ratios, we attempted to create
a combined Iron-and--element index similar to that defined by
González (1993) and more recently by Thomas et al. (2003):
![]() |
(3) |
Such a combined indicator is expected to decrease the effect of variation in the

The results are shown in Figs. 8a-c. The NaI index of NGC 5077 (#12) and the CO index of NGC 6909 (#14) deviate from the main loci of the other galaxies.
To address this issue, we plotted the optical indices of these
galaxies versus the combined [MgFe]'' index
(Fig. 8d). This is an analog of the typical plot
(i.e. Worthey 1994) that allows us to disentangle the
age-metallicity degeneracy: the inverse
index is approximately
proportional to age, while the [MgFe]'' is dominated by metal
abundance. The two galaxies exhibit strong
, which suggests that
they are dominated by populations of 3 Gyr or younger (see for
example Fig. 1 in S08). NGC 5077 and NGC 6909 are
also well separated from the remainder of the galaxies in the
versus Na I plot (Fig. 8e); this leads us to the
conclude that the NIR indices can be used to create a similar
diagnostic plot, but measuring the Br
feature requires higher
quality data than presented here.
We investigated the behavior of the NIR indices with respect to the H and K-band magnitude and H-K color, but no clear correlations were found.
5.3 Combined NIR metal index
Ground-based NIR spectroscopic observations are far more time
consuming to acquire than the corresponding optical observations -
higher and variable background and detectors with poorer cosmetics
often require one to sacrifice either resolution or signal-to-noise
ratio to obtain the data in a reasonable amount of time. To alleviate
this problem, we defined a combined spectral index of all the major
K-band metal features:
![]() |
(5) |
Various weighting schemes were attempted to minimize the scatter of the basic relations. However, the simple average yielded the tightest relations (Fig. 9). The carbon and oxygen are




The new index improves the tightness of the correlations. For example,
the Pearson correlation coefficient is 10% higher for the
versus Mg2 relation and
20%
higher for the
versus [MgFe]''
relation, with respect to the respective correlations in which only
one IR index is used. We derived these relations using only galaxies
with known parameters, i.e. excluding the poorly studied
NGC 4281(#4), NGC 5077(#12) and
NGC 6909(#14) which show systematic peculiarities with
respect to the bulk of our sample.
5.4 Average NIR spectrum of the sample galaxies
Studies of composite stellar systems (i.e. galaxies hosting an AGN)
often need to subtract the contribution of the underlying galaxy.
This prompted us to create an average spectrum of the galaxies in our
sample. We used a homogeneous subset of eight galaxies, excluding the
objects with young populations. We also excluded those with spectra of
low signal-to-noise ratio. The average age of the galaxies used to
create the composite spectrum was Gyr (the median was
9.5 Gyr) and the average [Fe/H] was
(the median was
0.16). The composite spectrum is shown in Fig. 10
and the measurements completed for the spectrum are listed in
Table 5, including r.m.s. values per Angstrom. We
included the H-band section because the artifacts caused by the
spectral-type mismatch of the standards are minimized by the averaging
of galaxies observed at different redshifts.
![]() |
Figure 10:
Average H- a) and K-band b) spectrum of sample
galaxies normalized to unity at 22 000 Å (thick line). The
thin lines correspond to the
|
Open with DEXTER |
![]() |
Figure 11:
Spectra of the galaxies observed with ISAAC in the region
of the Mg I (1.51 |
Open with DEXTER |
Table 5: Average spectrum of the sample galaxies observed at NTT in arbitrary units and normalized to unity at 22 000 Å. The regions with zero rms at the beginning and the end are covered only by one or two spectra. The full table is only available in the electronic version at the CDS.
5.5 The Mg I feature at 1.51
m
Ivanov et al. (2004) pointed out the possibility of using the Mg I
feature at 1.51 m as an
-element abundance indicator in
the NIR. For exploratory purposes, we obtained spectra of a small and
diverse set of galaxies (Table 2 and
Fig. 11). To carry out a quantitative analysis, we
defined an index (Table 3), and for five galaxies, we
measured its equivalent width to be 3.3, 5.0, 3.8, 3.9, and 4.5 Å,
respectively, for Mrk 1055, NGC 1144,
NGC 1362, NGC 4472, and NGC 7714. The
typical uncertainty was
0.3 Å. No corrections for velocity
dispersion were applied. Given the diverse nature of these objects, we
refrain from drawing any conclusions, but we note that the relation
between the optical and NIR Mg features is not straightforward because
the only two Mg2 values present in the literature are inversely
proportional to our index measurements for NGC 1362 and
NGC 4472. Nevertheless, these observations prove that it is
feasible to measure the NIR Mg I feature at 1.51
m and
provide a basis for future NIR synthetic spectral modeling.
6 Summary
We have presented a new data set of NIR spectra of ellipticals/spheroids to provide independent constraints on the properties of their stellar populations. They cover a wide range in luminosity and velocity dispersion, and have available Lick/IDS line-strength indices from literature.
- 1.
- The stronger features were measured and showed mild correlations with the optical metal features. NIR versus NIR features showed significant scatter but the diagram of CO versus Na I appears to be a good diagnostic for detecting different chemical enrichment histories.
- 2.
- The NIR metal absorption, in combination with a hydrogen absorption feature may be able to break the age-metallicity degeneracy, but further investigation and higher quality data are required to verify this possibility.
- 3.
- A new combined NIR index was defined that allows one to
interpret lower signal-to-noise data. It is dominated, at about the
2/3 level, by
-elements. We also present an average spectrum of ellipticals/spheroids, constructed from a homogeneous subset of our sample that can be used as a template integrated early-type galaxy spectrum.
- 4.
- The strength of Mg I feature at 1.5
m was measured for the first time in a heterogeneous sample of galaxies, providing a new constraint for future NIR spectral synthesis models.
Acknowledgements
We thank Reynier Peletier for useful discussion and suggestions. We acknowledge grant CPDA068415/06 by Padua University, which provided support for this research. M.C. thanks the European Southern Observatory for support via ESO Studentship program. L.M. is supported by grant (CPDR061795/06) by Padova University. L.M. acknowledge the Pontificia Universidad Catolica de Chile and the European Southern Observatory for hospitality while this paper was in progress. We are grateful to the La Silla and Paranal staff astronomers for their support.
References
- Ali, B., Carr, J. S., Depoy, D. L., Frogel, J. A., & Sellgren, K. 1995, AJ, 110, 2415 [NASA ADS] [CrossRef] (In the text)
- Alonso-Herrero, A., Rieke, M. J., Rieke, G. H., & Shields, J. C. 2000, ApJ, 530, 688 [NASA ADS] [CrossRef]
- Annibali, F., Bressan, A., Rampazzo, R., Zeilinger, W. W., & Danese, L. 2007, A&A, 463, 455 [NASA ADS] [CrossRef] [EDP Sciences]
- Bender, R., Burstein, D., & Faber, S. M. 1993, ApJ, 411, 153 [NASA ADS] [CrossRef] (In the text)
- Bernardi, M., Renzini, A., da Costa, L. N., et al. 1998, ApJ, 508, L143 [NASA ADS] [CrossRef]
- Bernardi, M., Alonso, M. V., da Costa, L. N., et al. 2002, AJ, 123, 2990 [NASA ADS] [CrossRef]
- Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AJ, 90, 1681 [NASA ADS] [CrossRef]
- Boisson, C., Coupé, S., Cuby, J. G., Joly, M., & Ward, M. J. 2002, A&A, 396, 489 [NASA ADS] [CrossRef] [EDP Sciences]
- Burstein, D., Faber, S. M., Gaskell, C. M., & Krumm, N. 1984, ApJ, 287, 586 [NASA ADS] [CrossRef] (In the text)
- Burston, A. J., Ward, M. J., & Davies, R. I. 2001, MNRAS, 326, 403 [NASA ADS] [CrossRef]
- Cenarro, A. J., Gorgas, J., Vazdekis, A., Cardiel, N., & Peletier, R. F. 2003, MNRAS, 339, L12 [NASA ADS] [CrossRef]
- Coziol, R., Doyon, R., & Demers, S. 2001, MNRAS, 325, 1081 [NASA ADS] [CrossRef]
- Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS All Sky Catalog of point sources, The IRSA 2MASS All-Sky Point Source Catalog, NASA/IPAC Infrared Science Archive, http://irsa.ipac.caltech.edu/applications/Gator/
- da Costa, L. N., Willmer, C. N. A., Pellegrini, P. S., et al. 1998, AJ, 116, 1 [NASA ADS] [CrossRef]
- Dalle Ore, C., Faber, S. M., Jesus, J., Stoughton, R., & Burstein, D. 1991, ApJ, 366, 38 [NASA ADS] [CrossRef]
- Dallier, R., Boisson, C., & Joly, M. 1996, A&AS, 116, 239 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Das, M., & Jog, C. J. 1999, ApJ, 527, 600 [NASA ADS] [CrossRef] (In the text)
- Davidge, T. J., Beck, T. L., & McGregor, P. J. 2008, ApJ, 677, 238 [NASA ADS] [CrossRef] (In the text)
- de Vaucouleurs, G., de Vaucouleurs, A., Corwin, Jr., H. G., et al. 1991, Third Reference Catalogue of Bright Galaxies (Berlin: Springer-Verlag)
- Denicoló, G., Terlevich, R., Terlevich, E., et al. 2005, MNRAS, 356, 1440 [NASA ADS] [CrossRef]
- Devillard, N. 1999, in Astronomical Data Analysis Software and Systems VIII, ed. D. M. Mehringer, R. L. Plante, & D. A. Roberts, ASP Conf. Ser., 172, 333 (In the text)
- Engelbracht, C. W. 1997, Ph .D. Thesis, Univ. of Arizona
- Faber, S. M. 1973, ApJ, 179, 731 [NASA ADS] [CrossRef]
- Faber, S. M., Friel, E. D., Burstein, D., & Gaskell, C. M. 1985, ApJS, 57, 711 [NASA ADS] [CrossRef] (In the text)
- Faber, S. M., Trager, S. C., Gonzalez, J. J., & Worthey, G. 1999, Ap&SS, 267, 273 [NASA ADS] [CrossRef]
- Falco, E. E., Kurtz, M. J., Geller, M. J., et al. 1999, PASP, 111, 438 [NASA ADS] [CrossRef]
- Förster Schreiber, N. M. 2000, AJ, 120, 2089 [NASA ADS] [CrossRef] (In the text)
- Frogel, J. A., Stephens, A., Ramírez, S., & DePoy, D. L. 2001, AJ, 122, 1896 [NASA ADS] [CrossRef]
- Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, Space Sci. Rev., 123, 485 [NASA ADS] [CrossRef] (In the text)
- Goldader, J. D., Joseph, R. D., Doyon, R., & Sanders, D. B. 1995, ApJ, 444, 97 [NASA ADS] [CrossRef]
- Goldader, J. D., Joseph, R. D., Doyon, R., & Sanders, D. B. 1997, ApJS, 108, 449 [NASA ADS] [CrossRef]
- González, J. J. 1993, Ph. D. Thesis, Univ. of Califonia, Santa Cruz (In the text)
- Hanson, M. M., Conti, P. S., & Rieke, M. J. 1996, ApJS, 107, 281 [NASA ADS] [CrossRef] (In the text)
- Heckman, T. M. 1980, A&A, 87, 152 [NASA ADS] (In the text)
- Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 112, 315 [NASA ADS] [CrossRef] (In the text)
- Howell, J. H. 2005, AJ, 130, 2065 [NASA ADS] [CrossRef]
- Ivanov, V. D. 2001, Ph. D. Thesis, Univ. of Arizona (In the text)
- Ivanov, V. D., Rieke, G. H., Groppi, C. E., et al. 2000, ApJ, 545, 190 [NASA ADS] [CrossRef]
- Ivanov, V. D., Rieke, M. J., Engelbracht, C. W., et al. 2004, ApJS, 151, 387 [NASA ADS] [CrossRef] (In the text)
- Keel, W. C. 1996, ApJS, 106, 27 [NASA ADS] [CrossRef]
- Kleinmann, S. G., & Hall, D. N. B. 1986, ApJS, 62, 501 [NASA ADS] [CrossRef] (In the text)
- Kobayashi, C., & Arimoto, N. 1999, ApJ, 527, 573 [NASA ADS] [CrossRef]
- Lancon, A., & Rocca-Volmerange, B. 1992, A&AS, 96, 593 [NASA ADS] (In the text)
- Larkin, J. E., Armus, L., Knop, R. A., Soifer, B. T., & Matthews, K. 1998, ApJS, 114, 59 [NASA ADS] [CrossRef]
- Lauberts, A., & Valentijn, E. A. 1989, The Messenger, 56, 31 [NASA ADS]
- Lu, N. Y., Hoffman, G. L., Groff, T., Roos, T., & Lamphier, C. 1993, ApJS, 88, 383 [NASA ADS] [CrossRef]
- Maiolino, R., Rieke, G. H., & Rieke, M. J. 1996, AJ, 111, 537 [NASA ADS] [CrossRef] (In the text)
- Mannucci, F., Basile, F., Poggianti, B. M., et al. 2001, MNRAS, 326, 745 [NASA ADS] [CrossRef] (In the text)
- Mehlert, D., Thomas, D., Saglia, R. P., Bender, R., & Wegner, G. 2003, A&A, 407, 423 [NASA ADS] [CrossRef] [EDP Sciences]
- Meyer, M. R., Edwards, S., Hinkle, K. H., & Strom, S. E. 1998, ApJ, 508, 397 [NASA ADS] [CrossRef] (In the text)
- Moorwood, A., Cuby, J.-G., Biereichel, P., et al. 1998a, The Messenger, 94, 7 [NASA ADS] (In the text)
- Moorwood, A., Cuby, J.-G., & Lidman, C. 1998b, The Messenger, 91, 9 [NASA ADS] (In the text)
- Morelli, L., Pompei, E., Pizzella, A., et al. 2008, MNRAS, in press [arXiv:0806.2988]
- Murphy, Jr., T. W., Soifer, B. T., Matthews, K., Kiger, J. R., & Armus, L. 1999, ApJ, 525, L85 [NASA ADS] [CrossRef]
- Murphy, Jr., T. W., Soifer, B. T., Matthews, K., Armus, L., & Kiger, J. R. 2001, AJ, 121, 97 [NASA ADS] [CrossRef]
- O'Connell, R. W. 1986, PASP, 98, 163 [NASA ADS] [CrossRef]
- Ogando, R. L. C., Maia, M. A. G., Pellegrini, P. S., & da Costa, L. N. 2008, AJ, 135, 2424 [NASA ADS] [CrossRef]
- Oliva, E., Origlia, L., Kotilainen, J. K., & Moorwood, A. F. M. 1995, A&A, 301, 55 [NASA ADS]
- Origlia, L., Moorwood, A. F. M., & Oliva, E. 1993, A&A, 280, 536 [NASA ADS] (In the text)
- Pickles, A. J. 1998, PASP, 110, 863 [NASA ADS] [CrossRef] (In the text)
- Pych, W. 2004, PASP, 116, 148 [NASA ADS] [CrossRef] (In the text)
- Rampazzo, R., Annibali, F., Bressan, A., et al. 2005, A&A, 433, 497 [NASA ADS] [CrossRef] [EDP Sciences]
- Reunanen, J., Kotilainen, J. K., & Prieto, M. A. 2002, MNRAS, 331, 154 [NASA ADS] [CrossRef]
- Reunanen, J., Tacconi-Garman, L. E., & Ivanov, V. D. 2007, MNRAS, 382, 951 [NASA ADS]
- Rinn, A. S., Sambruna, R. M., & Gliozzi, M. 2005, ApJ, 621, 167 [NASA ADS] [CrossRef] (In the text)
- Sánchez-Blázquez, P., Forbes, D. A., Strader, J., Brodie, J., & Proctor, R. 2007, MNRAS, 377, 759 [NASA ADS] [CrossRef]
- Silva, D. R., Kuntschner, H., & Lyubenova, M. 2008, ApJ, 674, 194 [NASA ADS] [CrossRef] (In the text)
- Smith, R. J., Lucey, J. R., Hudson, M. J., Schlegel, D. J., & Davies, R. L. 2000, MNRAS, 313, 469 [NASA ADS] [CrossRef]
- Sosa-Brito, R. M., Tacconi-Garman, L. E., Lehnert, M. D., & Gallimore, J. F. 2001, ApJS, 136, 61 [NASA ADS] [CrossRef]
- Stephens, A. W., & Frogel, J. A. 2004, AJ, 127, 925 [NASA ADS] [CrossRef] (In the text)
- Terlevich, R., Davies, R. L., Faber, S. M., & Burstein, D. 1981, MNRAS, 196, 381 [NASA ADS]
- Thomas, D., Maraston, C., & Bender, R. 2003, MNRAS, 339, 897 [NASA ADS] [CrossRef] (In the text)
- Trager, S. C., Worthey, G., Faber, S. M., Burstein, D., & Gonzalez, J. J. 1998, ApJS, 116, 1 [NASA ADS] [CrossRef]
- Trager, S. C., Faber, S. M., Worthey, G., & González, J. J. 2000, AJ, 120, 165 [NASA ADS] [CrossRef]
- Trager, S. C., Worthey, G., Faber, S. M., & Dressler, A. 2005, MNRAS, 362, 2 [NASA ADS] [CrossRef]
- van Dokkum, P. G. 2001, PASP, 113, 1420 [NASA ADS] [CrossRef] (In the text)
- Vanzi, L., & Rieke, G. H. 1997, ApJ, 479, 694 [NASA ADS] [CrossRef]
- Vanzi, L., Alonso-Herrero, A., & Rieke, G. H. 1998, ApJ, 504, 93 [NASA ADS] [CrossRef]
- Wallace, L., & Hinkle, K. 1997, ApJS, 111, 445 [NASA ADS] [CrossRef] (In the text)
- Wegner, G., Bernardi, M., Willmer, C. N. A., et al. 2003, AJ, 126, 2268 [NASA ADS] [CrossRef]
- Worthey, G. 1994, ApJS, 95, 107 [NASA ADS] [CrossRef]
Footnotes
- ...
spectroscopy - Based on observations made with ESO Telescopes at the La Silla and Paranal Observatory under programmes 69.B-0311 and 077.B-0163.
- ...
- Full Table 5 is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/497/41
- ... (IRAF)
- Image Reduction and Analysis Facility is a general purpose software system for the reduction and analysis of astronomical data; IRAF is written and supported by the IRAF programming group at the National Optical Astronomy Observatories (NOAO).
All Tables
Table 1: Basic properties of the sample galaxies.
Table 2: Log of the spectroscopic observations.
Table 3: NIR spectral indices.
Table 4: Equivalent widths of the newly defined NIR spectral indices for the sample galaxies.
Table 5: Average spectrum of the sample galaxies observed at NTT in arbitrary units and normalized to unity at 22 000 Å. The regions with zero rms at the beginning and the end are covered only by one or two spectra. The full table is only available in the electronic version at the CDS.
All Figures
![]() |
Figure 1:
The telluric absorption correction: the wavelength
calibrated 1-dimensional galaxy spectrum (a) is divided by the
wavelength calibrated 1-dimensional observed spectrum of the
standard star (b) and the result is multiplied by the reference of
the same spectral type of the standard star (c) from the library
of Pickles (1998). The product (d) is the corrected galaxy
spectrum. The large scatter at |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
H-band spectra of the sample galaxies, continuum-divided
(i.e. all continua were normalized to 1) and shifted vertically
for display purposes by adding increasing shifts of 0.25. The
position of the main spectral features is shown with dotted lines.
The passbands of the Mg I (1.71 |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
K-band spectra of the sample galaxies, continuum-divided
(i.e. all continua were normalized to 1) and shifted vertically
for display purposes by adding increasing shifts of 0.25. The
positions of the main spectral features are indicated by dotted
lines. The passbands of the Na I (2.21 |
Open with DEXTER | |
In the text |
![]() |
Figure 4: Correction factors of the NIR spectral index for different values of velocity dispersion. They were derived by comparing the EW measured in the spectra of the labeled stars before and after broadening them to a particular velocity dispersion. The spectral types of the stars are in a range between KIII and MIII with a spread of [Fe/H] between -0.35 and +0.41. In each panel, the solid line is the 2nd order polynomial fitting the data. |
Open with DEXTER | |
In the text |
![]() |
Figure 5: The equivalent widths of the NIR spectral indices, iron abundance [Fe/H], and metal abundance [Z/H] in the sample galaxies as a function of the central velocity dispersion. The labels correspond to the row numbers in Table 1. The open and filled circles refer to quiescent and active galaxies, respectively. The small and large circles represent galaxies with an age of between 5 and 10 Gyr and >10 Gyr, respectively. Galaxies with no known ages were assumed to be old. The dashed lines corresponds to the relations found by S08. The solid line in each panel represents the best fit linear regression ( y = ax + b) for all the data points except for #4, #12 and #14. The Pearson correlation coefficient (r) and the results of the linear fit are given. Typical error bars are shown. |
Open with DEXTER | |
In the text |
![]() |
Figure 6: The equivalent widths of the NIR indices plotted against each other. The symbols are the same as in Fig. 5. The dashed lines correspond to the relations found by S08. The dotted line corresponds to the relation by Davidge et al. (2008) for giant stars in the solar neighborhood. All of these relations are plotted in their observed range. The solid line represents the linear regression ( y = ax + b) through all data points except for #4, #12, and #14. The Pearson correlation coefficient (r) and the results of the linear fit are given. Typical error bars are shown. |
Open with DEXTER | |
In the text |
![]() |
Figure 7: The equivalent widths of the NIR indices plotted against the Mg2 indices measured by Bender et al. (1993). The symbols are the same as in Fig. 5. The solid line represents the linear regression ( y = ax + b) through all data points except for #4, #12, and #14. The Pearson correlation coefficient (r) and the results of the linear fit are given. Typical error bars are shown. |
Open with DEXTER | |
In the text |
![]() |
Figure 8:
a- d) The equivalent widths of the NaI, CaI, CO and
H |
Open with DEXTER | |
In the text |
![]() |
Figure 9:
The equivalent widths of the newly defined |
Open with DEXTER | |
In the text |
![]() |
Figure 10:
Average H- a) and K-band b) spectrum of sample
galaxies normalized to unity at 22 000 Å (thick line). The
thin lines correspond to the
|
Open with DEXTER | |
In the text |
![]() |
Figure 11:
Spectra of the galaxies observed with ISAAC in the region
of the Mg I (1.51 |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.