A&A 407, 905-912 (2003)
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
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 (BC) 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) 2.13 m, the red star is indicated by an arrow. c) 5300 Å, the two components of N88 are labelled A and B. d) The H2 emitting region with its main N-S filament and peak of intensity, is indicated by an arrow. The size of each field corresponds to or pc.|
|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/Hratio in N88A was found to be 4.2, corresponding to a mean extinction AV of 1.14 mag. In N88B the H/Hratio 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,
stars surrounded by an HII region can be derived from the flux of Hphotons
(Devereux & Scowen 1994)
|Figure 3: E-W rectified composite spectra of the ionizing source embedded in the core N88A, in a 4450-4750 wavelength range obtained in 1991 (grey line) and 1988 (black line). The spectrum of the component B is also presented (dotted line). The continuum level is normalized to unity.|
|Open with DEXTER|
|Figure 4: E-W matched red and blue CCD spectra obtained in 1997 (long exposures) through an aperture. The intensity is in units of 10-16 ergs cm-2 Å-1.|
|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
as derived from
includes the medium empty of electrons.
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
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.
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.