A&A 407, 105-120 (2003)
DOI: 10.1051/0004-6361:20030808
N. G. Guseva 1 - P. Papaderos 2 - Y. I. Izotov 1 - R. F. Green 3 - K. J. Fricke 2 - T. X. Thuan4 - K. G. Noeske2
1 - Main Astronomical Observatory,
Ukrainian National Academy of Sciences,
Zabolotnoho 27, Kyiv 03680, Ukraine
2 -
Universitäts-Sternwarte, Geismarlandstraße 11,
37083 Göttingen, Germany
3 -
National Optical Astronomy Observatory,
Tucson, AZ 85726, USA
4 -
Astronomy Department, University of Virginia,
Charlottesville, VA 22903, USA
Received 21 November 2002 / Accepted 16 May 2003
Abstract
We present a detailed optical spectroscopic and B, V, I, H
photometric study of the metal-deficient cometary blue compact dwarf (BCD) galaxy
SBS 1415+437. We derive an oxygen abundance
and
(
)
in the two brightest H II
regions,
among the lowest in BCDs. The helium mass fractions in these regions
are
and
.
Four techniques based on the equivalent widths of the hydrogen emission and
absorption lines, the spectral energy distribution and
the colours of the galaxy are used to put
constraints on the age of the stellar population in the low-surface-brightness
(LSB) component of the galaxy,
assuming two limiting cases of star formation (SF), the
case of an instantaneous burst and that of a continuous SF with
a constant or a variable star formation rate
(SFR). The spectroscopic and photometric data for different regions of the LSB
component are well reproduced by a young stellar population with an
age t
250 Myr, assuming a small extinction in the range
AV = 0-0.6 mag.
Assuming no extinction, we find that the upper limit for the mass of the old
stellar population, formed between 2.5 Gyr and 10 Gyr, is
not greater than
(1/20-1) of that of the stellar population formed
during the last
250 Myr. Depending on the region considered, this also
implies that the SFR in the most recent SF period
must be 20 to 1000 times greater than the SFR at ages
2.5 Gyr.
We compare the photometric and spectroscopic properties of SBS 1415+437 with
those of a sample of 26 low-metallicity dwarf irregular and BCD galaxies.
We show that there is a clear trend for
the stellar LSB component of lower-metallicity galaxies to be
bluer. This trend cannot be explained only by
metallicity effects. There must be also a change in the age of the
stellar populations. The most metal-deficient galaxies have also smaller
luminosity-weighted ages.
Key words: galaxies: abundances - galaxies: dwarf - galaxies: evolution - galaxies: compact - galaxies: starburst - galaxies: stellar content - galaxies: individual: SBS 1415+437
Since its discovery as a metal-deficient
blue compact dwarf (BCD) galaxy (Thuan et al. 1995),
SBS 1415+437 (CG 389)
has been considered as a probable nearby young dwarf galaxy. Situated at a
distance D = 11.4 Mpc it was classified by Thuan et al.
(1999) as a cometary BCD with a very bright supergiant H II region
at the SW tip of the galaxy.
From 4 m Kitt Peak National Observatory (KPNO) telescope spectra, Thuan et al.
(1995) first derived
an oxygen abundance of
in SBS 1415+437 placing the
galaxy among the most metal-deficient BCDs known. Later, Izotov & Thuan
(1998,1999) derived from the same spectrum
using five-level atom models for abundance determination
instead of the three-level atom model used by Thuan et al. (1995).
Thuan et al. (1999), using
Multiple Mirror Telescope (MMT) and Hubble
Space Telescope (HST) FOS observations, derived
and
,
respectively. The high brightness and
low metallicity of the H II region in SBS 1415+437 make this galaxy
one of the best objects for helium abundance determination. Izotov & Thuan
(1998) and Thuan et al. (1999) derived respectively a helium mass
fraction
and
for it, close to
the primordial helium mass fraction of
by
Izotov et al. (1999).
Thuan et al. (1999) have discussed the evolutionary status of
SBS 1415+437, using ground-based MMT spectroscopic and HST/WFPC2
photometric data. Based on the (V-I) vs. I colour-magnitude diagrams (CMD)
and spectral energy distributions (SED) in the optical range, they concluded
that SBS 1415+437 is a truly young galaxy that did not start to form stars
until 100 Myr ago. However, the V and I images used by Thuan et al.
(1999) were not deep enough for the detection of old red giant
branch (RGB) stars in the CMD.
Furthermore, they considered an instantaneous burst model which gives
only a lower limit to the age of the stellar population in SBS 1415+437.
In this paper we combine new spectroscopic and photometric data with previous observations to derive elemental abundances and to better constrain the age of the stellar population in SBS 1415+437. For the latter task we use four different techniques of age determination and consider different star formation (SF) histories. The paper is organized as follows. In Sect. 2 we describe the observations and data reduction. The photometric properties of SBS 1415+437 are described in Sect. 3. We derive in Sect. 4 the chemical abundances in the two brightest H II regions. In Sect. 5 we discuss the properties of the stellar populations in SBS 1415+437 and compare them with those in other low-metallicity dwarf galaxies. Finally, Sect. 6 summarises the main conclusions of this study.
Narrow-band images of SBS 1415+437 in the
H
line at
6563 Å
through a passband with a full width at half maximum
(
)
of 74 Å, and
in the adjacent continuum at
6477 Å through a passband with
Å were obtained with the Kitt Peak
2.1 m telescope
on April 22, 1999 during a photometric night.
The telescope was equipped with a Tektronix
CCD detector operating at a gain of 3 e- ADU-1, giving an
instrumental scale of 0
305 pixel-1 and field of view of 5
.
The total exposures of 50 min in the H
line and 40 min
in the
adjacent continuum bluewards of H
were split up into 5 and 4 subexposures, slightly offset with respect to each other for removal
of cosmic particle hits and bad pixels. The point spread function
has a
of 2
2.
Bias and flat-field frames were obtained during the same night.
The standard stars Feige 34 and HZ 44 were observed
in both filters during the same night at several airmasses for
absolute flux calibration.
Another broad-band B image (15 min) of
SBS 1415+437 was obtained on March 9, 1997 under
photometric conditions, with the CAFOS focal reducer
attached to the 2.2 m telescope of the German-Spanish Astronomical
Center, Calar Alto, Spain.
CAFOS was equipped with a SITe 2048
2048 CCD operating at a gain
of 2.3 e- ADU-1, with a read-out noise of <3 counts (rms).
With a focal ratio of f/4.4, the instrumental scale was 0
53 pixel-1
and field of view
15
.
The seeing during
the observations was 3
1 (
).
Standard stars from Christian et al. (1985) were observed at several
airmasses for calibration.
Standard reduction steps, including bias subtraction, flat-field
correction, removal of cosmic particle hits
and absolute flux calibration were carried out using
IRAF
and MIDAS
.
The ground-based photometric data were supplemented by the HST/WFPC2 V (F569W) and I (F791W) images described by Thuan et al. (1999).
![]() |
Figure 1: HST V image of SBS 1415+437. The long-slit positions during the two observations with the MMT and 4 m Mayall telescope are labeled "1'' and "2'', respectively. North is up and east is to the left. Regions e1-e3 with hydrogen emission lines in the spectra and regions a1-a4 with hydrogen absorption lines are labeled. |
Open with DEXTER |
![]() |
Figure 2:
Continuum-subtracted H![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Spectroscopic observations were carried out on June 18, 1999, at the
Kitt Peak 4 m
Mayall telescope with the Ritchey-Chrétien spectrograph and a
T2KB
CCD detector. The 2
300
slit was centered on the brightest H II region e1
(slit 2 in Fig. 1) with position angle
so as to include the second brightest H II region
e2 to the SW of region e1.
We used the KPC-10A grating in first order and a GG 375 order separation filter.
The spatial scale along the slit
was 0
69 pixel-1 and the spectral resolution
7 Å (
).
The spectra were obtained at an airmass 1.27. The total
exposure time of 60 min was broken up into 3 subexposures.
No correction for atmospheric refraction was made because of
the small airmass during the observations.
Two Kitt Peak
spectrophotometric standard stars were observed for flux calibration.
For wavelength calibration, He-Ne-Ar comparison spectra were obtained
after each exposure.
The data reduction was performed with the IRAF software package. This includes bias-subtraction, flat-field correction, cosmic-ray removal, wavelength calibration, night sky background subtraction, correction for atmospheric extinction and absolute flux calibration of the two-dimensional spectrum.
For abundance determination, one-dimensional spectra of regions
e1 and e2
were extracted within
apertures of 2
4
6 and 2
4
0,
respectively.
In addition, we extracted spectra of the low-surface-brightness (LSB)
regions a3 and a4 with strong hydrogen
Balmer absorption lines,
and of region e3 with H
and H
in emission (Fig. 1).
We also used the two-dimensional MMT spectrum obtained by Thuan et al. (1999)
with the slit oriented at
(slit 1 in Fig. 1).
We extracted one-dimensional spectra of the LSB regions a1 and
a2 within apertures of 1
5
3
4 and
,
respectively (Fig. 1).
These spectra show strong hydrogen Balmer absorption lines.
The selected LSB regions are listed in Tables 3 and 4 with their positions and spatial extents. Origins are set at the center of the brightest region e1 (Figs. 1, 2). The spectra of the LSB regions are used to study the stellar populations and to constrain the age of the oldest stars which contribute to the light of these regions.
It is seen from the continuum-subtracted H
images (Fig. 2)
that star-forming activity in SBS 1415+437 is primarily occurring
in regions e1 and e2, with some additional faint H
emission present in region e4.
However, the available narrow-band data are not deep enough for tracing
faint H
emission in other regions of the LSB component, such as in
regions a1, a3 and e3 where H
and H
have
been detected spectroscopically (Table 3).
The latter fact suggests that some low-level SF may be present at
various locations within the LSB component.
The photometric properties of the stellar LSB component of SBS 1415+437
were first investigated by Thuan et al. (1999) using
HST/WFPC2 V and I images.
These authors fit the surface brightness profiles (SBPs)
of SBS 1415+437 with an exponential distribution
in the radius range
with a scale length
(
0.3 kpc).
However, their study was limited by the fact
that the HST images do not include the outermost NE part of the LSB
component.
Furthermore, the SBPs in Thuan et al. (1999) reach only
a surface brightness level
mag arcsec-2.
It is known, however, that the star-forming component may contribute
to the optical BCD emission to fainter surface brightness levels
(see e.g. Papaderos et al. 2002 and references therein).
A comparatively young stellar population has
been observed in the inner part of the LSB component of other cometary BCDs,
several hundred pc away from the brightest H II region
(Noeske et al. 2000; Guseva et al. 2001; Fricke et al. 2001).
Therefore deeper images are needed to study the outer parts of
the LSB component in SBS 1415+437.
Note that the SBPs by Thuan et al. (1999) show
in the outermost part (
,
or
)
a steeper
exponential intensity decrease than the one observed at intermediate intensity levels.
This slope change, not discussed in Thuan et al. (1999), is found
independently by us in the HST/WFPC2 SBPs derived with
the method iv of Papaderos et al. (2002) and
ellipse fitting to the visible
part of the LSB component. While
the latter method extends surface photometry out to
larger radii (
), it is subject to large uncertainties
because the NE part of the LSB component with
mag arcsec-2 (
)
lies outside the HST/WFPC2
field of view.
Additionally, the outermost LSB isophotes show considerable
deviations from ellipticity.
![]() |
Figure 3:
a)
Surface brightness profiles (SBPs) of SBS 1415+437 in
V and I (filled and open circles, respectively), derived from
HST/WFPC2 data. A linear fit to the V SBP in the radius interval
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
To study the surface brightness distribution at large radii we
use the ground-based Calar Alto B image. Despite the poor
spatial resolution this image allows us to study the
entire LSB component out to its Holmberg radius.
The change in the exponential slope for
is confirmed from the ground-based B SBP.
At large radii, however, this SBP reveals a flatter, outer exponential
regime with a scale length
fairly comparable to that previously
obtained at intermediate intensity levels from HST data
(upper thick-grey line in Fig. 3a).
From fitting an exponential model to the B SBP for
(lower thick-grey line in Fig. 3a) we obtain a central surface brightness
B mag arcsec-2 and a scale length
kpc.
Note, however, that the inner exponential profile studied by Thuan et al. (1999) is 0.3 mag brighter than the outer one, which suggests that
more than 1/4 of the emission associated with this profile originates
from the part of the LSB component
between regions e1 and a2.
The present data provide no compelling evidence for a large age difference
between the stellar population which dominates within the inner exponential regime
discussed in Thuan et al. (1999) and that responsible for the outermost LSB
emission (i.e. for R*>16
).
The (V-I) profile reveals a strong colour
gradient (
mag kpc-1; inner solid-grey
line in Fig. 3b) within the inner 5
,
or
roughly the V band effective radius
(
)
of SBS 1415+437. At larger radii, however,
linear fits to the (V-I) profile yield, depending on whether they
are error-weighted or not, a gradient not exceeding 0.1 and 0.2
mag kpc-1, respectively. Comparable values are also found from subtraction
of the exponential fits in V and I (Fig. 3a) in the
radius range
(Fig. 3b,
solid-grey line at
intermediate radii).
![]() |
Figure 4: The KPNO 4 m telescope spectra of the brightest H II regions e1 and e2 with the identified emission lines. The lower spectra in a) and b) are the observed spectra downscaled by factors of 50 and 30, respectively. |
Open with DEXTER |
Table 1:
Observed (F()) and extinction-corrected
(I(
)) fluxes and equivalent widths (EW) of emission lines
in the H II regions e1 and e2.
Table 2: Element abundances in regions e1 and e2.
In this section we derive the elemental abundances of regions e1 and e2 using the Kitt Peak 4 m telescope observations. Their spectra with strong emission lines are shown in Fig. 4.
The observed (F()) and extinction-corrected
(I(
)) emission line fluxes relative to the H
emission line
fluxes, their equivalent widths EW, the extinction coefficients
C(H
), the observed fluxes of the H
emission line, and the
equivalent widths of the hydrogen absorption lines
for regions e1 and e2 are shown in Table 1.
Despite the differences in aperture
(
for the Kitt Peak 4 m data,
and
for the MMT data from Thuan et al. 1999),
the relative fluxes of the emission lines for region e1
are in agreement within the errors with those derived by Thuan et al. (1999).
The physical conditions and heavy element abundances in regions e1
and e2 were derived
following Izotov et al. (1994,1997a) and Thuan et al. (1995).
The electron temperatures (O III),
(S III),
(O II) for the high-, intermediate- and low-ionization
regions, respectively, the electron number
densities
(S II), ionization correction factors (ICF), and
ionic and total heavy element abundances are shown in Table 2 for
both regions.
The oxygen abundance
(
/20)
and heavy element abundance ratios for region e1 are in good agreement
with those derived by
Thuan et al. (1999). The oxygen abundance
and heavy element abundance ratios in region e2
are consistent with those for region e1 within the errors.
In Fig. 5 we compare the heavy element abundance
ratios in the two brightest regions of
SBS 1415+437 with data for a sample of low-metallicity BCDs.
The Ne/O, S/O, Ar/O and [O/Fe] abundance ratios for the comparison sample are taken from
Izotov
Thuan (1999), while the Cl/O abundance ratios are collected from
Izotov
Thuan (1998) and Izotov et al. (1997a).
The heavy element abundance ratios for regions e1 (large squares)
and e2 (small squares) are in good agreement with those for other BCDs.
Note, that the Cl/O ratio does not show any significant increase with increasing oxygen
abundance. This conclusion is strengthened by the observations of
Esteban et al. (1998, 1999a, 1999b) who derived
log(Cl/O) in the range from -3.28 to -3.47 for
high-metallicity H II regions
in Orion, M17 and M8 with
,
8.50 and 8.60, respectively. For comparison,
is derived for the Sun (Anders & Grevesse 1989).
![]() |
Figure 5: Comparison of the elemental abundance ratios, obtained for the brightest H II regions e1 (large squares) and e2 (small squares) with data for other BCDs (open circles). |
Open with DEXTER |
The high brightness of regions e1 and e2
allows for a reliable determination of the 4He abundance. Nine He I
emission lines are detected in the spectrum of region e1 (Table 1). Two
of them, He I 3889 and
4713, are blended with other
lines. Six He I lines are detected in region e2.
The five brightest He I
3889,
4471,
5876,
6678,
7065 emission lines are used to correct
their fluxes for collisional and fluorescent enhancement. This is
done by minimizing the deviations of the corrected He I line
flux ratios from the recombination ratios, through varying the electron number
density in the He+ zone and the optical depth in the He I
3889
emission line. The flux of this line was preliminarily corrected for the
contribution of the H I
3889 emission line, according to
prescriptions of Izotov et al. (1994,1997a).
Helium abundances He+/H+, derived from the corrected He I
4471,
5876,
6678 line fluxes and their weighted
mean are listed in Table 2. The abundance He+2/H+ is
added to He+/H+ for region e1, as He II
4686 is present in its
spectrum. Note the lower He abundance derived from the He I
4471 flux which is most likely due to significant
underlying stellar He I
4471 absorption.
The effect of underlying absorption for the other He I emission
lines used in the He abundance determination is much smaller, as they
have much larger equivalent widths compared to the He I
4471 emission line (Table 1).
The mean 4He mass
fractions
and
in regions e1 and
e2 (Table 2) are
consistent with the values derived for SBS 1415+437 by Izotov & Thuan (1998)
and Thuan et al. (1999). They are also consistent with the primordial
4He mass fraction
,
derived by extrapolating the
Y vs. O/H linear regression to O/H = 0 (Izotov & Thuan 1998), or to
derived from spectroscopic observations of
the two most metal-deficient BCDs known,
I Zw 18 and SBS 0335-052 (Izotov et al. 1999).
We consider next the spectroscopic and photometric properties of
the LSB regions labeled a1, a2 (slit
position 1 in Fig. 1) and a3, e3, a4 (slit position 2), to constrain the age of the stellar populations
contributing to the light in those regions.
H
and H
emission lines are present in regions a1,
a3 and e3 while H
and H
absorption
lines are detected in the spectra of all regions except for region e3.
This
allows us to derive the age of the stellar population using four
methods, based on: (1) the time evolution of equivalent widths (EW)
of hydrogen emission lines, (2) the time evolution of EWs of hydrogen
absorption lines, (3) the comparison of the observed
and theoretical spectral energy distributions, and (4) the broad-band
colours.
The requirement of consistency of the ages determined from the
reddening-insensitive methods 1 and 2 and from the
reddening-sensitive methods 3 and 4 allows to simultaneously derive
the extinction coefficient and constrain the SF history (Guseva et al.
2001,2003a,2003b).
We measured the fluxes and equivalent widths of the
H
and H
emission lines and the H
and H
absorption
lines in the spectra of the LSB regions, and list them in
Tables 3 and 4.
Because the H
emission line is narrower than the absorption line in
these regions and does not fill the absorption component, its flux was
measured using the continuum level at the bottom of the absorption line.
This level was chosen by visually interpolating from the absorption line
wings to the center of the line.
Table 3:
Fluxes, equivalent widths of the H
and H
emission
lines and the extinction coefficients C(H
)
in the LSB regions.
Table 4:
Equivalent widths of the H
and
H
absorption lines in the LSB regions.
The extinction coefficient C(H)
in those regions is derived
from the H
/H
flux ratio.
We adopt the theoretical recombination H
/H
flux ratio of
2.8, which is typical for hot low-metallicity H II regions.
No correction for the absorption line equivalent widths has been made.
The extinction coefficients C(H
)
are shown in Table 3.
Hydrogen absorption lines are seen in the spectra of all
regions labeled in Fig. 1 except for the
brightest H II regions e1 and e2
and the LSB region e3.
The blue part of the spectrum of region a1 with hydrogen absorption
and emission lines is shown in Fig. 6.
Table 4 lists the equivalent widths with their errors of
the H
and H
absorption lines measured
in the wavelength intervals or "windows'' of Bica & Alloin (1986).
The errors include
the errors in the fitting of line profiles with Gaussians and the noise
dispersion in the
continuum. A careful placement of the continuum level is very important
for deriving accurate EWs.
For this purpose, we choose points in the spectrum free of nebular and stellar lines,
which were then fitted by cubic splines. The uncertainties
were estimated from several different measurements of the
equivalent widths of hydrogen absorption lines with independent continuum
fittings. They are of the same order as the errors in
Table 4.
![]() |
Figure 6: The blue part of the MMT spectrum of region a1 with labeled emission and absorption lines. |
Open with DEXTER |
The calibration of the age of stellar populations using
the equivalent widths
of the H
and H
nebular emission lines,
those of the H
and H
stellar absorption lines and the spectral energy distributions is
discussed in
detail in Guseva et al. (2001,2003a,2003b). Here
we only briefly describe these calibrations.
The temporal evolution of the H
and H
emission
line equivalent widths depends on the star formation history.
We consider here the two limiting cases of instantaneous burst
and continuous star formation models.
The equivalent widths for the instantaneous burst model
with a heavy element mass fraction
/20
are calculated
using the galactic evolution code PEGASE.2 (Fioc & Rocca-Volmerange
1997). The dependence of the H
emission line equivalent
width on time is shown in Fig. 6a of Guseva et al. (2003b)
by the thick solid line.
These models are appropriate for regions e1 and e2 with strong
emission lines.
The equivalent widths of Balmer emission lines in region e1
(EW(H
Å and EW(H
Å)
and region e2 (EW(H
Å and EW(H
Å) correspond to an instantaneous burst
age of 4 Myr.
However, for the LSB regions, models with continuous
star formation are more appropriate. For these models we adopt
a constant star formation rate (SFR)
within the time interval from
when star formation starts to
when it stops. Time is zero now and increases to the past.
The equivalent widths of hydrogen emission
lines and SEDs for a set of instantaneous burst models (Fioc &
Rocca-Volmerange 1997) are used to calculate the temporal evolution of
EWs for continuous SF with a constant
SFR. The temporal dependence of the equivalent widths of the
H
emission line is shown in Fig. 6a of Guseva et al.
(2003b) for different
and
.
Another way to derive the age of a stellar population is to use
the relation between the H
and H
absorption line
equivalent widths and age, derived
by González Delgado et al. (1999).
Their instantaneous burst models predict a steady increase of the equivalent
widths with age from 1 Myr to 1 Gyr. However,
they did not extend the calculations for ages
1 Gyr when the equivalent
widths of the absorption lines decrease with age (Bica & Alloin
1986). Hence, each value of the hydrogen absorption
line equivalent width corresponds to two values of the age,
1 Gyr
and
1 Gyr. This ambiguity can be resolved with the use of other age
constraints discussed in this paper.
Furthermore, the models by González Delgado et al. (1999)
probably overestimate the equivalent widths of the absorption lines at ages
1 Gyr (Guseva et al. 2003b).
Therefore, in the age range from 1 Myr to 16.5 Gyr instead of the
calibration by González Delgado et al. (1999)
we use an empirical calibration of the hydrogen absorption line
equivalent widths versus age by Bica & Alloin (1986).
This calibration is based on the integrated spectra of 63 open and globular
stellar clusters with
known ages, metallicities and reddenings
which can be used as templates for
stellar populations formed in an instantaneous burst.
For consistency we use the same wavelength intervals or "windows''
for H
and H
flux measurements as Bica & Alloin (1986)
(
4318-4364 Å and
4082-4124 Å, respectively).
![]() |
Figure 7:
Spectra of region a2 on which
synthetic SEDs are superposed.
Synthetic SEDs are calculated for stellar populations forming continuously
during one a) or two ( b)- g)) time intervals. In the case
of two intervals the SFR is constant
within each interval but varies from one interval to another one by a
factor b = SFR(young)/SFR(old).
The spectra in left panel are corrected for interstellar
extinction with C(H![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 8:
a) V surface brightness distribution
along the slit oriented at position angle PA = 22![]() ![]() ![]() ![]() |
Open with DEXTER |
Table 5: Predicted equivalent widths of emission and absorption lines, extinction coefficients and (V-I) colours from models of continuous star formation.
The behaviour of the empirical H
absorption line equivalent width
with the age for an instantaneous burst (Bica & Alloin
1986) is shown in Fig. 6b of Guseva et al. (2003b)
by the thick solid line.
The temporal evolution of the H
and H
absorption line equivalent widths in the case of continuous SF
is calculated similarly to that of the H
and H
emission line equivalent widths described in Sect. 5.1.1.
More specifically, we use the empirical equivalent widths of hydrogen
absorption lines (Bica & Alloin 1986)
and SEDs for instantaneous bursts
(Fioc & Rocca-Volmerange 1997)
to calculate the temporal evolution of EWs in the case of continuous
SF with constant SFR. The results are shown in
Fig. 6b of Guseva et al. (2003b) for
SF with different
and
.
The shape of the spectrum reflects the properties of the
stellar population. However, it is also dependent on the reddening.
A precise determination of the extinction can be done only for the
two brightest H II regions e1 and e2 which possess
many strong hydrogen emission lines (Table 1).
We derived an extinction coefficient C(H
in these regions.
In the LSB regions a1, a3 and
e3, only H
and H
emission lines are present.
The extinction coefficients obtained from the H
/H
flux ratio in
these regions are small (Table 3). However, they are more uncertain as
compared to the ones in regions e1 and e2 because of the weakness
of the emission lines and significant contribution of the stellar
absorption lines.
H
and H
emission lines are not detected in the other LSB regions.
Therefore, the observed SED cannot directly give information on the age,
but should be used together with the methods discussed in Sects. 5.1.1
and 5.1.2 for simultaneous determination of the age and
interstellar extinction.
We used the galactic evolution code PEGASE.2 (Fioc &
Rocca-Volmerange 1997) to produce a grid of theoretical SEDs
for an instantaneous burst of star formation with ages ranging between
0 and 10 Gyr, and a heavy element mass fraction of Z = /20.
The SEDs for continuous SF in the time interval between
ago and
ago are derived by integration of instantaneous burst SEDs.
Table 6: (V-I) colours of the extended LSB components in some irregular and BCD galaxies.
In this section we derive self-consistently the ages of the stellar populations
in the LSB regions using:
1) the equivalent widths of emission lines, 2) the equivalent widths of
absorption lines, 3) the SEDs and 4) the colours.
For this we adopt a continuous
SF scenario with constant or variable SFR. In the latter case we
consider a simplified scenario with two time intervals of SF, which we
call young and old, with
different SFRs. To quantify the difference in SFRs we introduce the parameter
(young)/SFR(old). Because of the noisy spectrum of region
e3 and the significant contamination by nebular emission we were not able
to measure equivalent widths of absorption lines
in that region. Therefore, we exclude it from the analysis. For the
remaining LSB regions, the model predictions providing the best fits to the observed
SEDs are shown in Table 5. For each region we show the
age range of one or two SF episodes, the parameter b, the mass ratio
/
of the young-to-old populations,
the model equivalent widths of hydrogen lines,
the extinction coefficients C(H
)
derived from the best
match between observed and calculated SEDs,
and the intrinsic and reddened (V-I) colours (see Sect. 5.2.1).
We use the
relation E
(V-I) = 0.891 C(H
)
(e.g.,
Aller 1984) to correct for reddening.
Our young and old stellar populations in Table 5 include stars
with ages not older than 250 Myr (
)
and not younger than 160 Myr
(
), respectively. These definitions differ from those conventionally used.
In fact, the old stellar population in our case includes not only several
Gyr old stars but also intermediate-age stars with age
1 Gyr.
Negative EWs denote absorption lines,
positive EWs refer to emission lines. Models with highlighted EWs are
those in which the equivalent
widths are not reproduced well despite a good fit of the observed SEDs
(compare model EWs in Table 5 with observed EWs in
Tables 3 and 4).
![]() |
Figure 9:
Dependence of the (V-I) colours of the LSB component a)
and logarithm of the derived ages b)
on the oxygen abundance for the dwarf irregular and BCD galaxies from our sample.
Dashed lines in a) show theoretical dependences of
the (V-I) colour on the oxygen abundance for an instantaneous burst calculated
using the galactic evolution code PEGASE.2 (Fioc ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
We show in Fig. 7
the observed
and predicted SEDs in the outer LSB region a2 (MMT
observations) with no emission features.
The H
and H
absorption
lines in this region as well as in the other outer
region a4 (4 m telescope observations)
are not contaminated by nebular emission from young stellar populations.
Hence,
the absorption line equivalent widths in regions
a2 and a4 are more accurate than in other LSB regions
(Table 4).
We consider several SF histories and vary extinction in the LSB regions to put constraints on the age of their stellar populations. Figure 7 shows the spectra of region a2 on which are superposed the synthetic SEDs with different SF histories which best fit it.
We have assumed two extinction coefficients:
C(H
,
derived
by Thuan et al. (1999), and corresponding to
,
and C(H
.
The model SEDs, adopting C(H
are
shown in Fig. 7, left panel.
The observed properties of region a2 can be reproduced by a
single young stellar
population with age between 12 Myr and 160 Myr, if a small amount of extinction
is assumed (Fig. 7a).
There is no need to invoke an older stellar population which,
if present, has to be much smaller in mass than the young stellar population.
We note that there is no ongoing SF in this region.
The most recent SF in region a2
stopped at time
12 Myr ago, otherwise
its spectrum would have shown
a detectable H
emission line (Fig. 6a in
Guseva et al. 2003b).
Next we consider SF scenarios with a varying SFR which include older stellar populations with an age of up to 2 Gyr (Fig. 7b-7c) and 10 Gyr (Fig. 7d). But even in these cases, the young population completely dominates the light and mass of region a2 (Table 5).
An upper limit to the age of region a2 can be obtained by assuming
no extinction, i.e. C(H)
= 0 (Fig. 7,
right panel).
Then the observed SED cannot be reproduced by a
synthetic SED with stars forming between 12 Myr and 160 Myr ago, as the
latter is too blue. However, by varying the parameter b, the observed
SED can be reproduced by SEDs of stellar populations with other SF scenarios.
In Figs. 7e-g, we show the best
fits labeled by the adopted value of b.
We exclude the models shown in
Figs. 7c and f because
they do not fit
the observed equivalent widths of H
and H
.
The remaining two
models with C(H
satisfy all observational constraints.
In these models, the SFRs for the young population are respectively
7 times (Fig. 7e) and 17 times
(Fig. 7g) larger than SFRs for the old population.
These ratios are, however,
significantly smaller than the corresponding parameters b in the case with
C(H
.
The relative mass fraction of the young stellar
population is therefore smaller in the extinction-free case,
being
/
(Fig. 7e)
and 0.3 (Fig. 7g),
respectively. Therefore, the presence of a 2 Gyr old or even 10 Gyr old
population is not excluded in region a2, which would however not
dominate the optical emission. If the extinction is non-zero, which is
likely the case, then there is no need to invoke a significant old population
to explain the observed properties of region a2. It is likely that
for other regions, a1 and a3 (Table 5), the extinction
is even larger.
Similar consideration for region a4 shows that the mass fraction of the old stellar population is small even in the absence of extinction. In any case, the old stellar population, if present does not contribute significantly to the luminosity. All the observed properties of this region can be reproduced with a young stellar population formed between 12 Myr and 160 Myr ago. However, the presence of an older stellar population cannot be excluded in that case, but one needs to increase the parameter b to match the observations (Table 5).
If some extinction is present in region a4, then our age
upper limits will decrease.
There is some hint that the extinction may increase with
decreasing distance to the brightest H II region e1. Indeed,
by considering models which include also a 10 Gyr old population and
C(H)
= 0,
we find from Table 5 that the closer a region is
to the brightest part of the galaxy
(region a1 compared to region a2, region a3 compared to
region a4), the smaller is the relative mass
/
of the young stellar population. This is in contrast to the trend
found in BCDs
where the relative mass of the young stellar population decreases
outwards. Hence, we conclude, that the extinction in SBS 1415+437 is
larger in the inner brighter regions.
In summary, our spectroscopic analysis of the LSB regions
favors a relatively young luminosity-weighted age of the stars
populating those regions.
A model with a constant SFR continuing
from 10 Gyr ago until now is definitely excluded.
An old population 2 Gyr could be present only in models with
very specific SF scenarios, with a very low SF activity during
the first 2-10 Gyr period, a very high star formation rate during the last
(10-200) Myr, and a quiescent period in between. If, however, low SF
activity has occurred in the period between
200 Myr and 2 Gyr ago,
then there is no need to invoke a stellar population with age
2 Gyr, and all the spectroscopic properties of the LSB regions
can be explained by only young and intermediate-age
stellar populations.
From the HST images, we derived V and I surface brightness and colour distributions for the regions covered by the spectroscopic observations with the MMT (Figs. 8a and b) and the KPNO 4 m telescope (Figs. 8c and d). Note the trend for the (V-I) colour of the LSB component to decrease with increasing distance from region e1, especially for the KPNO 4 m data (Fig. 8d). This is again suggestive of larger extinction in the brighter regions. We compare the observed colours with predictions from our population synthesis modeling. The results of this comparison are shown in Figs. 8b and d. The predicted colours, obtained from convolving the theoretical SEDs with the appropriate filter bandpasses, are shown by different symbols. The transmission curves for the Johnson V and Cousins I bands are taken from Bessell (1990). The zero points are from Bessell et al. (1998).
Since the contribution of ionized gas emission to the total brightness of
region e1 is significant,
the theoretical SED for this region has been constructed
using a 4 Myr old stellar population SED for a heavy element mass fraction
to which the gaseous continuum SED and the
observed emission lines were added
(see Guseva et al. 2001).
For comparison with the observed colour at
,
we reddened
the predicted colour of region
e1 by an amount corresponding to C(H
(Table 3 in Thuan et al. 1999),
and at
we
adopted C(H
(Table 1).
The observed colour of region e1 is very blue, (V-
(Figs. 8b and d), and cannot be
reproduced by a 4 Myr old stellar population alone ((V-
).
Strong gaseous continuum and line emission need to be added
(Table 1).
On the other hand, in all the LSB regions the contribution of the gaseous
emission to the total flux is negligible.
In Figs. 8b and d we show by different symbols the
colours of the LSB regions with various star formation
history (Table 5).
The total colour (stellar plus gaseous emission) of region e1
is shown by a filled squares.
The agreement between the (V-I) colours
and those derived from the spectral energy
distributions is very good.
However, some uncertainties in the colours may be introduced by
the uncertainties in the reddening curves and the ratios of total to selective
extinction
R=A(V)/E(B-V). We use R = 3.2 by Aller (1984), based
on the reddening curve by Seaton (1979).
Schlegel et al. (1998) give a slightly different R = 3.315.
Nevertheless, because the extinction in the studied regions is small, the errors
introduced by the use of different reddening curves and Rare less than 2% in the V band and negligible in the I band.
The (B-V) colour of 0.2 mag, derived from the
SEDs of the LSB regions, is also consistent with the observed value,
derived from the B and V SBPs, which are shown in Fig. 3a.
Our comprehensive studies of selected galaxies
with oxygen abundances
and
blue LSB components
(SBS 0335-052, Izotov et al. 1997b, Papaderos et al. 1998;
I Zw 18, Izotov et al. 2001b, Papaderos et al. 2002;
SBS 0940+544, Guseva et al. 2001;
Tol 1214-277, Fricke et al. 2001;
Tol 65, Papaderos et al. 1999;
SBS 1129+576, Guseva et al. 2003a;
HS 1442+4250, Guseva et al. 2003b;
SBS 1415+437, this paper) have led us to the conclusion that these galaxies might
be young. This is in contrast
to the large age estimates for some well-studied higher-metallicity
irregular and BCD galaxies, such as VII Zw 403 and
UGCA 290 (Schulte-Ladbeck et al. 1998; Crone et al. 2002).
To investigate this apparent inconsistency we therefore
compare the (V-I) colours of the LSB components ((V-I)
)
of the galaxies from our
sample with those in galaxies where a large age was derived from
colour-magnitude diagrams (CMD).
In Table 6, we show the oxygen abundances and (V-I) colours of the
LSB components in several dwarf irregular and BCD galaxies.
The telescopes used for
photometric and spectroscopic observations are also given in the table.
Ground-based photometric data in V and I were obtained with the 2.1 m
KPNO telescope by Y. Izotov and R. Green. The HST photometric data
were retrieved from the archive of the Space Telescope Science Institute
(STScI).
We included in the sample only galaxies with low internal and/or foreground
interstellar extinction to avoid uncertainties introduced by the
correction of the (V-I) colour for reddening.
Thus, we do not include the nearby star-forming galaxy NCG 1569
studied with HST, for example.
The chemical abundances for the majority of galaxies are
obtained from spectroscopic observations of
their H II regions. An exception
is I Zw 18C, the faint component of the BCD I Zw 18, where
no emission lines of heavy elements
were detected.
Therefore, for I Zw 18C we adopt the oxygen abundance
derived for the bright main body of I Zw 18.
(V-I) colours of the galaxies in Table 6 are corrected for
reddening in our Galaxy, with E(V-I)
taken from the NASA/IPAC Extragalactic Database (NED).
The dependence of the (V-I)
colours
on oxygen abundance for 26 selected galaxies
is shown in Fig. 9a.
Thin dashed lines represent theoretical dependences
for an instantaneous burst
in the age range from log
to log
= 10.0
(
in yr).
The models for continuous star formation are shown
by thick solid lines.
These models are calculated for a constant SFR which started at time
,
with log
between 7.5 and 10.0 and continuing until now
(
). The blueing of the (V-I) colour with decreasing
oxygen abundance is in agreement with model predictions. However, this blueing
trend for galaxies with
is too
steep to be explained only by metallicity effects.
The blue colours of the LSB components of these low-metallicity
BCDs
are also not due to ionized gas emission, as the latter
dominates the outer parts of only two BCDs: I Zw 18
(Izotov et al. 2001b; Papaderos et al. 2002) and SBS 0335-052 (Izotov et al. 1997b;
Papaderos et al. 1998). In other galaxies,
including I Zw 18C, the LSB component emission has mainly a stellar origin.
The steep trend cannot be explained by
reddening effects because the interstellar extinction derived
from the spectroscopic observations is small for the galaxies
shown in Fig. 9a.
Therefore, it is likely that the blueing is mainly due to a change in
the age of the stellar populations, and that
low-metallicity galaxies are younger than high-metallicity ones.
All objects with
and (V-I)
mag
are in the range of colours predicted
for ages
1-2 Gyr by instantaneous and continuous
models of star formation. On the other hand, galaxies with
(V-I)
mag are likely to be older,
with ages
10 Gyr. Yi (2003)
has analysed the uncertainties in the synthetic
integrated (V-I) colours caused by uncertainties in the stellar
evolutionary models, population synthesis techniques and stellar atmosphere
models.
He has shown that, despite all the uncertainties,
a (V-
mag can be
attributed to an intermediate-age stellar population with an age not larger
than 2 Gyr, in the case of an instantaneous burst of star formation,
in agreement with the age derived here from the PEGASE.2 models
(dashed lines in Fig. 9a).
The smaller luminosity-weighted age of galaxies with blue
LSB components in Fig. 9a is supported by the detailed
analysis of the spectroscopic and photometric properties of some of these
galaxies, as demonstrated in this series of papers.
Izotov & Thuan (1999) have suggested that the oxygen abundance may be a good age indicator. That this appears to be the case is shown in Fig. 9b where a clear trend of increasing age of the LSB stellar population with increasing oxygen abundance is seen for 5 star-forming galaxies with available data. The ages for 4 galaxies (I Zw 18, Izotov et al. 2001b; SBS 0940+544, Guseva et al. 2001; HS 1442+4250, Guseva et al. 2003b; SBS 1415+437, this paper) are determined by the four methods described before. The age for SBS 0335-052 (Papaderos et al. 1998) is derived using only colours and SEDs. This is because its LSB component is embedded into the H II region produced by the young central clusters and no absorption lines are present in the spectra.
The results of
a detailed photometric and spectroscopic study
of the metal-deficient blue compact dwarf galaxy SBS 1415+437
are presented. H
images and spectra in the optical range have been obtained
with the Kitt Peak 2.1 m and 4 m telescopes, respectively.
A B image has been obtained with the 2.2 m Calar Alto telescope.
These data are
supplemented by HST/WFPC2 V and I images and MMT spectra from
Thuan et al. (1999).
The main conclusions of this study can be summarized as follows:
Assuming no extinction, we find that the upper limit to the
mass of the old
stellar population in SBS 1415+437, formed between 2.5 Gyr and 10 Gyr, is
not greater than (1/20-1) of the mass of the stellar population formed
during the last
250 Myr. Depending on the region considered, this also
implies that the star formation rate during the
most recent star formation
in SBS 1415+437 must be 20 to 1000 times greater than the SFR at ages
2.5 Gyr.
Acknowledgements
N.G.G. has been supported by DFG grant 436 UKR 17/2/02 and Y.I.I. acknowledges the Göttingen Academy of Sciences for a Gauss professorship. N.G.G. and Y.I.I. have been partially supported by Swiss SCOPE 7UKPJ62178 grant. They are grateful for the hospitality of the Göttingen Observatory. Y.I.I. and T.X.T. have been partially supported by NSF grant AST-02-05785. Research by P.P. and K.J.F. has been supported by the Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) under grant 50 OR 9907 7. K.G.N. acknowledges the support from the Deutsche Forschungsgemeinschaft (DFG) grants FR 325/50-1 and FR 325/50-2. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.