A&A 407, 905-912 (2003)
DOI: 10.1051/0004-6361:20030974
G. Testor1, -
J. L. Lemaire2,
-
D. Field3,
1 - LUTH, UMR 8102 du CNRS, Observatoire de Paris, 92195 Meudon, France
2 -
LERMA, UMR 8112 du CNRS, Observatoire de Paris, 92195 Meudon, France
and Université de Cergy-Pontoise, 95031 Cergy Cedex, France
3 -
Department of Physics and Astronomy, Århus University, 8000
Århus C, Denmark
Received 1 April 2003 / Accepted 12 June 2003
Abstract
We present optical, near-infrared (NIR) images and spectra of the
compact, dusty
and complex HII region N88A in the SMC. The data were obtained mainly in the
visible with the Boller Chivens (B
C) spectrograph at the ESO 1.5 m telescope and recently with the
infrared
spectrometer and array camera (ISAAC) at the VLT, around 2
m.
The NIR observations identified a region emitting in H2, both
in and around the two major components of N88. Using long-slit spectra
obtained in the 3600-10 100 Å wavelength range we have re-determined the key physical
parameters of electron temperature and densities. The spectral type of the
complex exciting source of N88A
has been classified as an O6-O8 V, using He and H lines in the visible and NIR, with
the HeI 1.7
m/2.11
m line ratio in the K band spectrum, suggesting that
N88A is effectively ionized by a source of type <O7 V. However, from reddening correction and the temperature
determined by photoionization model calculations, the ionizing source should also contain a very early type
O star.
A comparison with a recent photometric
study shows that the V magnitude of the exciting source of N88A has previously been
overestimated. The chemical abundances of He, O, N, Ne, S, and Ar, from emission-line intensities,
are also derived. Abundancies are consistent with the average abundances for SMC HII regions.
Key words: galaxies: Magellanic Clouds - ISM: individual objects: N88A (SMC) - ISM: HII regions - ISM: abundances - infrared: ISM - stars: early-type
The recent discovery of the heavily reddened compact HII region N26A-B
(Testor 2001) in the SMC has strengthened interest in continuing a detailed investigation of
N88A, a similar neighbour region extensively studied since its discovery by Testor & Pakull
(1985). Like N26A-B, N88A has a close component B, only 5
distant and less
excited. N88A is classified as a member of a rare class of compact HII regions (
pc in H
)
(Testor 2002), generally associated with more extended
HII regions detected by Henize (1956). N88A is strongly excited and is among the brightest
objects in the Magellanic Clouds (MCs).
The first highly excited compact HII region in the MCs, now generically named
HEBs, was discovered two decades ago (Heydari & Testor 1982) in a study of individual
regions in a low metallicity environment. In spite of their small angular size, as the
clouds are sufficiently close, HEBs are the best candidates for the study of individual
star forming regions, both in detail and globally. HEBs should be the final stage of ultra
compact HII regions where new massive stars are forming (Churchwell 1990). Up to now,
eight HEBs have been found in the MCs. Save for N26 A-B, their characteristics
are reviewed by Heydari et al. (2002). Israel & Koorneef (1988, 1991) have shown that the
N88 region has a strong IR continuum with additional emission from shocked H2.
Wilcots (1994a,b) has presented radio continuum imagery,
optical imagery and spectroscopy as well as UBV photometry.
N88A was intensively studied by Kurt et al. (1999) using UV, optical and HST
observations both in imaging and spectroscopy. Kurt et al. described this region, formed from
a dusty molecular cloud, as a very young HII region, cut by a dark band and
excited by at least two stars.
They noted that the reddening, especially in the UV, must be different from
the SMC extinction and presented an accurate determination of elemental abundances. These were
found to be similar to values in other recent studies, except that C was found
overabundant. Subsequently, thanks to the refurbished HST, Heydari et al. (1999, hereafter
HM99) give a detailed description of the complex morphology of N88A, especially of the
bright and dusty central region where the extinction is estimated to have a peak
AV of
3.5 mag
at the position of the dark lane. In a region of extent
0.3
,
west of the dark lane, HM99 found three faint exciting stars arriving on the zero age
main sequence that could belong to the so called Vz luminosity class (Walborn
Parker 1992).
However, despite such detailed studies, little can be said about the
massive stars in N88A except on the basis of the H
flux. The faint main
exciting star remains unidentified. Up to now, no spectrum of the ionizing cluster has been
published. Vermeij & van der Hulst (2002) have revisited the physical structure of the region. They suggest that abundance determinations suffer from a number of
uncertainties in line fluxes mainly due to the extinction correction, a point which we confirm here.
Massey (2002) presents an accurate UBV photometry of the region, except for the A
and B components of N88.
It is evident from these numerous studies of N88 that many fundamental quantities await accurate determination. In the present paper we work towards this goal, presenting results, accumulated over several years, involving optical and NIR spectroscopic and imaging data, largely for the bright region N88A. The paper is organized as follows: Sect. 2 gives a brief description of old and new data as well as the photometry of N88A and its nearby environment. Section 3 briefly outlines the morphology of the region and Sect. 4 the extinction and the spectral classification of the ionizing source of N88A, obtained from Lyman photons and from spectroscopy. In addition physical parameters such as electron density, temperature and mass of ionized gaz are given. Section 5 describes the chemical abundances followed by a summary of results in Sect. 6.
Table 1 lists all the observations of long-slit optical and NIR spectra
obtained during four observing runs at the ESO 1.5 m telescope with the BC spectrograph, as
follows:
- In December 1988 a dispersion of 114 Å mm-1, giving a spectral coverage from 3400 to 5400 Å and a FWHM of 4.8 Å.
- In October 1991 a high resolution
CCD Ford at
per pixel corresponding to
on the sky and a grating with a dispersion of 114 Å mm-1, giving a large spectral coverage
from 3600 to 7400 Å and a FWHM of 4.8 Å.
Table 1:
Journal of the spectrographic observations conducted with the BC
spectrograph.
- In January and December 1995 the same
CCD Ford but with a
holographic grating, giving a reciprocal dispersion of 30 Å mm-1. In the range
4000-5200 Å, this grating yielded a
spectral resolution
1.1 Å.
- In December 1997 the Loral
m chip
with two different gratings of dispersion 114 Å mm-1 from 3600 to 7400 Å and of 116 Å mm-1 from 6500 to 10 100 Å giving a FWHM of
7.3 Å.
Both short and long exposures were taken because of the problem of overexposure
in the very strong [O III]
5007 and H
lines. The 2
and
4
long-slit widths were oriented E-W.
Spectroscopic data were reduced mainly with the Midas software package LONG available at Paris-Meudon Observatory. Corrections such as flat-fielding and sky substraction were applied in the standard manner. The spectra were wavelength calibrated with a helium-argon arc. The flux was calibrated using observations of the spectrophotometric standard stars Feige 24 and 110 (Oke 1974).
![]() |
Figure 1:
Extracted images of monochromatic CCD frames in the optical and NIR continuum a) K-band image with stars labelled. b)
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The Infrared Spectrometer and Array Camera (ISAAC) installed on the ESO-VLT UT1
has been used for imaging in the infrared.
Details of the data and reduction methods will be described in a subsequent (Paper II).
Images of N88A were obtained using the Ks broad band filter (Fig. 1a) and two
narrow-band filters at 2.13
and at 2.19
,
both with 1.3% bandwidth, where data
for the former filter are shown in Fig. 1b. For the bright objects, the background continuum is dominant at both
wavelengths. Below, we refer to these bright images in the continuum as
cont
and cont
respectively. In addition,
in the ionized region a Helium line may appear at 2.11
m, in the blue wing
of the 2.13
m filter. On the faint objects in the 2.13
m filter, the 1-0 S(1) H2
line at 2.121
m may in some cases be seen in addition to the continuum (see Sect. 3.1).
The 2.19
m filter, which does not contain any H2 lines, shows only the
continuum emission.
In the visible, a continuum image (Fig. 1c), taken through a narrow-band filter
(
Å) centered at 5300 Å (cont
), was also used to
compare the optical and NIR continuum. This cont
image
was obtained from a CCD (
chip,
47),
attached to the Danish 1.5 m telescope at la Silla. Table 2 summarizes
the log of CCD frames. The images were reduced using MIDAS.
Table 2: Log of CCD images.
From the K band and cont5300 images the instrumental
magnitudes of a few stars were computed
using DAOPHOT/ALLSTAR (Stetson 1987) in a field of
centered on N88A (Fig. 1a). The FWHM of N88A is 2.3
and the point spread function
.
Therefore the integrated magnitude of N88A was derived using circular apertures.
Photometric calibration in the K band was obtained using the reference
star HD 7644 of K
(van der Bliek et al. 1996), whereas for
the cont5300 we used the star of coordinates [1.24 04.89-73 09 15.0
(J2000)] and V = 14.70 found in Massey (2002).
Table 3 shows magnitudes of stars in the region, taken from literature values,
and the present cont
and K magnitudes (Cols. 15 and 16). Assigned spectral types
are also shown. Table 3 also contains, as available, the UBV
photometry obtained by Massey (2002) (Cols. 4-6) and the IJHK photometry found
in the infrared data base of Denis (Cioni et al. 2000) (Cols 7-9) and 2MASS
(IPAC/UMass, 2000) (Cols. 10-12). In Cols. 13 and 14 are shown V and
Strömgren y magnitudes from Wilcots (1994b) and HM99. The magnitudes derived
by Wilcots and HM99 compared with Massey's photometry seem overestimated by at least one
magnitude. Table 3 represents an extract of a larger VLT-ISAAC field for which the
complete K band photometry will be published in Paper II.
Table 3:
N88A and a few neighbouring stars: values of the
magnitudes U, B, V, I, J, K, y (Strömgren 5475 Å) found in the
literature and values of the cont5300 and K magnitudes from the present
work. Identification corresponds
to the numbering in Fig. 1a.
H2 emission has been detected in N88A. Detection was achieved by substracting
the image in the 2.19 m filter, which allows the passage only of continuum
radiation, from the image in the 2.13
m filter (see Sect. 2.2). Data were corrected
for the different responses of the instrument
through the filters, through observations of stars in the field. The
resulting image is shown in Fig. 1d. The very precise removal of all stars in
this image indicates accurate continuum subtraction. The strong emission directly over
N88A has an unknown contribution from the HeI line at 2.11
m. The H2
contribution will be discussed in Paper II. However around N88A there is highly extended
emission, due purely to H2 in three roughly equally spaced arms. The length of each of these arms is 6 to 7 pc.
The extent of the region containing H2 emission is similar to the extent of H2 emitting
regions in highly luminous massive star forming regions in the Galaxy, for
example W51 (e.g. Hodapp & Davis 2002), which has previously been likened to
extremely active star-forming regions in other galaxies e.g. the 30 Doradus in
the LMC (Hodapp & Davis 2002). Further details of the H2 emitting region
around N88A will be reported in Paper II, including K-band spectra of the brightest
zones.
Figure 1a shows, 8
to the north of N88A, a star labelled 33 which is
brighter in the infrared (Figs. 1a,b) than in the visible (Fig. 1c).
Using a distance modulus
m - M = 19.09 (Jacoby et al. 1990)
and the tables of Schmidt-Kaler (1982)
and Johnson
Morgan (1953), we derive for star 33 an absolute magnitude MV of -2.17 and a V-K of 2.31 (Table 3, Col. 17). From these parameters this
star may be classified as a red supergiant (RSG) of type G9-K1 II.
Using the
slit aperture, the H
/H
ratio in N88A was found to be 4.2, corresponding to a mean extinction
AV of 1.14 mag. In N88B the H
/H
ratio is 3.39 (
mag). With a 2
aperture the H
/H
ratio is 4.65 (
AV = 1.42 mag).
Each ionizing source of A and B is of type O (Table 3). However the N88A/N88B
ratio derived from cuts oriented N88A-N88B (Fig. 1c) is 45 (4.13 mag)
in the cont
or cont
images and 5.5 (1.85 mag) in the
cont5300 image, as illustrated in Fig. 2.
The magnitude excess of A in the cont
image is estimated to be
2.3 mag, giving for N88A a mean AV of 3.44 mag.
This strong mean AV found from NIR observations could be interpreted as
dense internal dust associated with
the central core of N88A and should have a local extinction
3.44 mag. Our estimated
extinction strengthens the argument for the strong AV of at least 6 mag
estimated by HM99 at the position of the faintest star of the ionizing source.
In connection with nebular reddening,
it is clear in the cont
image (Fig. 1c) that the intensity ratio
of the N88A to the N88B components is lower than in the
cont
image (Fig. 1b). Thus N88A is differentially more reddened than
N88B, consistent with the greater contribution of longer wavelength emission in N88A.
![]() |
Figure 2: Brightness distribution of the 5300 Å (dashed line) and NIR continuum (solid line) images along the direction A-B. |
Open with DEXTER |
HST observations of HM99 have shown that the core A contains, in a
diameter of 0
3, two stars
1 and
2 and a third fainter
star of
20 mag, the latter more dust embedded, just to the east of the brightest
core A. These three stars of integrated magnitude
M5470 = 17.5 mag should be the
sources responsible for the ionization of N88A.
As our V magnitudes of the stars compared with Wilcots (1994b) and HM99
appear one magnitude brighter, the integrated V magnitude of the ionizing
source should be
mag rather than
obtained by HM99 above
the nebular continuum. Corrected for the reddening (
AV = 3.70 mag, Sect. 4.1), the V magnitude is expected to be 12.8 mag. Using Schmidt-Kaler's tables and a distance modulus of 19.09, such a magnitude suggests that the ionizing source could be a very early star
of type O3.
The Lyman continuum photon rate in photons per second,
,
produced by
stars surrounded by an HII region can be derived from the flux of H
photons
(Devereux & Scowen 1994)
![]() |
(1) |
![]() |
Figure 3:
E-W rectified composite spectra of the ionizing source embedded in the core N88A, in a 4450-4750 ![]() |
Open with DEXTER |
![]() |
Figure 4:
E-W matched red and blue CCD spectra obtained in 1997 (long exposures) through an
![]() |
Open with DEXTER |
Table 4:
Forbidden and permitted SMC N88A
B emission line intensities
normalized to H
.
To determine the spectral type of the faint embedded ionizing source, we
used a spectrum resulting from an average of five east-west 1800 s exposures. A
single 1D-spectrum was extracted
by summing a pixel range along the 2
slit width, corresponding
to 4
.
This 2
aperture was centered on N88A, which is of
in H
.
The two components A and B are spatially resolved and the contribution of
B is negligible.
In Fig. 3 the E-W spectrum of the component A, in the range 4450-4750 Å, clearly
shows the HeII absorption lines
4541 and
4686.
It is not possible to discriminate between absorption and emission in the He I
4471 line, even if this
line is analyzed with the higher spectral resolution spectrum
Å obtained in 1995 (Table 1).
Comparing the composite spectrum shown in Fig. 3 with that in the atlas of
Walborn & Fitzpatrick (1990), we estimate that the exciting source of N88A has an equivalent spectral type of O6-O8 V.
The spectral type classification is strengthened by the NIR spectrum (Paper II), where the ratio
of 2.11 m He I to Br
.
According to Hanson et al. (2002), this
ratio is consistent with an ionizing source
of spectral type
O7 V.
Data in Fig. 3 also show the spectrum of N88B and in addition allow the
classification of the ionizing star of N88B, which appears to be a late O star.
Physical parameters were computed using the 1996 version of Abelion
developed by Stasinska. This program uses a five-level atom code and the atomic
data are described by Stasinska
Leitherer (1996). The emission-line
intensities (Fig.4) were corrected for reddening using the observed H
/H
ratios and Seaton's law (Seaton 1979) for N88A, whereas for the fainter
component N88B we used the H
/H
ratio. The reddening-uncorrected
lines normalized to H
are listed
in Table 4 with the corresponding signal to noise ratio (S/N) derived from the method
given by Newberry (1994).
Four electron temperatures were derived, associated with
different ionization regions in the nebula. These temperatures were
computed from the following line intensity ratios: [O III]
4363/
(4959+5007),
[N II]
5755/
(6548+6584),
[O II]
3727/
(7319+7330) and [S III]
6312/
9532.
Two electron densities were determined from the [S II]
6717/
6731 and
[Cl III]
5517/
5537 ratios.
In the determination of the physical parameters obtained from S and N lines,
only the [O III]
4363/
(4959+5007) and [S II]
6717/
6731 ratios were used.
The mean electronic densities,
,
and the mean electron
temperatures
[O III], are shown in the last five rows of Table 4 as well as temperatures and
densities obtained with other ratios, for comparison.
Density diagnostics such as [SII] line ratios measure the density of
the matter ([S II]) while the mean electron density
(Pereira 2000),
as derived from
,
includes the medium empty of electrons.
![]() |
(2) |
![]() |
(3) |
![]() |
(4) |
N88A is an HII region well adapted to computation of abundances. It is bright and its angular size is such that the Strömgren sphere fills the slit of the spectrograph.
Table 5 lists the ionic abundances of O+, O++, S+ and S++ for component A. The total abundances of Ne, N, S and Ar have been calculated using
ionization corrector factors (ICFs)
described by Vermej & van der Hulst (2002, Eqs. (6)-(8)) with
,
using the derived
[O III] and
[S II] values for the
nebula. Results are presented in Table 6. For comparison, Table 6 also contains
the latest
average HII abundances in the SMC derived by Reyes (1999) and Garnett (1999) as well as the N88A abundances derived from the studies of Kurt et al. (1999) and
Vermej (2002).
Small discrepancies are noticed between the abundances derived from 1991
and 1997 observations as well as between the mean abundances given
by Reyes (1999) and Garnett (1999). Despite our small error bars, it is
difficult to maintain that the scatter of results is physically real. Except
for S that appears underabundant,
all our variations are within the errors given by most of the different
authors. Our detection of a weak emission feature at [CI]
9850 (Table 4)
should, according to Rudy et al. (2001),
strengthen the conclusion that N88A is carbon rich in nature, as established by
Kurt et al. (1999).
Vilchez & Pagel (1988) defined a radiation softness parameter
from the ionic ratio
![]() |
(5) |
Based on optical spectra and NIR data:
Table 5: N88A Ionic abundances: 10-8 X/H.
- An H2 emitting region was detected in the direction of N88A, composed of 3 main filaments, with the strongest oriented N-S and located a few arcsecs to the east of N88A.
Table 6:
Element abundances in N88 A & B in the form (12+log[N(X)/N(H)]), derived from spectra obtained with 2
and 4
slit widths.
- The integrated magnitude of the ionizing source of N88A was found to be one V mag brighter than the value previously estimated.
- The first composite spectrum of N88A is presented, showing the HeII
absorption lines. This spectrum allows classification of the complex ionizing source as O6-O8 V.
This result is in agreement with the exciting
source stellar type derived from the H
luminosity and the HeI/Br
ratio
but does not allow us to determine with certainty if the source contains a very
early O star, as suspected from the strong extinction
detected from NIR observations and the
parameter.
- The presence of the red super-giant (33)
is consistent with the OB star population towards N88.
- We calculate, using accurate electron temperatures, electron densities and empirical ICFs, new elemental abundances of He, O, Ne, S and Ar. These values are found to be consistent with the mean average of the SMC. However our values present small discrepancies both with average SMC abundances and with earlier estimates of N88A abundances, which themselves display differences. The likely cause of discrepancies lies in uncertainties in reddening corrections.
- We confirm the carbon rich nature of N88A (Kurt et al. 1999), through the
detection of the [CI] 9850 emission line.
In conclusion we note that the nature of the exciting sources of N88A and B could be more precisely determined using high spatial resolution, NIR imaging and spectroscopy, for example with the adaptive optics capability which has recently become available on UT4 (NAOS-CONICA). In Paper II, we extend the present study to consider in detail the nature of the H2 emitting region in N88.
Acknowledgements
We thank the referee, Manuel Peimbert, for a careful reading of the manuscript and several suggestions, which significantly improved the paper. D.F would like to acknowledge the support of the Aarhus Centre for Atomic Physics (ACAP), funded by the Danish Basic Science Foundation. We should also like to thank the Director and Staff of ESO especially for making possible the IR observations reported in this paper.