Issue |
A&A
Volume 517, July 2010
|
|
---|---|---|
Article Number | A2 | |
Number of page(s) | 10 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200913294 | |
Published online | 23 July 2010 |
Near-infrared spectroscopy in NGC 7538![[*]](/icons/foot_motif.png)
E. Puga1 - A. Marín-Franch2 - F. Najarro1 - A. Lenorzer2 - A. Herrero2 - J. A. Acosta Pulido2 - L. A. Chavarría3 - A. Bik4 - D. Figer5 - S. Ramírez Alegría2
1 -
Centro de Astrobiología (CSIC-INTA), Ctra. de Torrejón a Ajalvir km-4, 28850, Torrejón de Ardoz, Madrid, Spain
2 - Instituto de Astrofísica de Canarias, C/ vía Láctea s/n, 38200, La Laguna, Spain
3 - Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
4 - Max-Planck-Institut für Astronomie, Königstuhl 17, 69117, Heidelberg, Germany
5 - Rochester Institute of Technology, 54 Lomb Memorial Drive, Rochester NY 14623, USA
Received 14 September 2009 / Accepted 18 February 2010
Abstract
Aims. The characterisation of the stellar population in
young high-mass star-forming regions allows fundamental cluster
properties like distance and age to be constrained. These are essential
when using high-mass clusters as probes for conducting Galactic
studies.
Methods. NGC 7538 is a star-forming region with an embedded
stellar population unearthed only in the near-infrared (NIR). We
present the first near-infrared spectro-photometric study of the
candidate high-mass stellar content in NGC 7538. We obtained H and K
spectra of 21 sources with both the multi-object and long-slit
modes of LIRIS at the WHT, and complement these data with subarcsecond
photometry of the region using the imaging mode of the same instrument.
Results. We find a wide variety of objects within the studied
stellar population of NGC 7538. Our results discriminate between a
stellar population associated to the H II region,
but not contained within its extent, and several pockets of more recent
star formation. We report the detection of CO bandhead emission toward
several sources, as well as other features indicative of a young
stellar nature. We infer a spectro-photometric distance of
kpc, an age spread in the range 0.5-2.2 Myr and a total mass
for the older population.
Key words: stars: early-type - stars: formation - stars: pre-main sequence - infrared: stars
1 Introduction
Clusters of total mass >10
are the natural habitat of high-mass stars. Whether motivated by the
study of the formation, evolution of feedback effects of high-mass
stars themselves (Bik et al. 2010; Puga et al. 2009; Martins et al. 2009), or as probes to investigate the physics and chemistry evolution of the Galaxy (Messineo et al. 2009; Davies et al. 2009),
NIR spectroscopic studies have been proven powerful at unveiling and
characterising the complete high-mass stellar content within massive
clusters. Although new NIR surveys, such as UKIDSS and VISTA, will
surely uncover new obscured massive clusters, some clusters that are
already known from optical observations harbour hidden high-mass stars
only accessible at NIR wavelengths (e.g. Cyg OB2, Knödlseder 2000).
NGC 7538 (aka Sh 2-158) is a visible H II
region in the Perseus spiral arm and a site of active star formation.
An early detection of several luminous NIR and far-IR sources in the
vicinity of this region (Campbell & Persson 1988) hinted at a rather massive nature. In particular, NGC 7538 IRS 1 is a high-mass (30
)
protostar with a CO outflow, known to power the ultra-compact (UC)H II
region NGC 7538 A. It also has an associated linear methanol
maser structure, which might trace a Keplerian-rotating circumstellar
disk (Pestalozzi et al. 2004).
Recently, 6.7 GHz methanol masers have been detected toward the
nearby objects NGC 7538 IRS 9 and NGC 7538 S,
tracing other young and embedded massive protostars (Pestalozzi et al. 2006).
Although this region has been widely inspected at long wavelengths
(mostly in the submillimetre window) and optical spectroscopy for two
stars has been obtained (Russeil et al. 2007), only a few detailed near-IR photometric studies of the stellar population in NGC 7538 have been conducted (Balog et al. 2004; McCaughrean et al. 1991; Ojha et al. 2004). At subarcsecond resolution, the region breaks down into several areas of interest (see Fig. 1):
(i) the presumed powering cluster, centred on the location of
IRS 6; (ii) the subcluster sitting amidst the bright NIR
reflection nebula (IRS 5); (iii) the cluster located at the
southern rim of the optical H II region
(IRS 4); (iv) the active region around IRS 1-3; (v) an
embedded stellar cluster at the location of NGC 7538 S, only
detected at wavelengths longward of 3.6 m;
and (vi) the IRS 9 reflection nebula.
Surprisingly, the stellar density detected in the NIR is highest at the
southern rim of the dust bubble, which bounds the optical H II region (Balog et al. 2004) between regions (iii) and (iv).
The distance to NGC 7538 has been estimated several times and values between 2.2 and 3.5 kpc have been reported (Israel et al. 1973; Moreno & Chavarria-K. 1986).
![]() |
Figure 1:
Colour composite image of NGC 7538 obtained with LIRIS ( blue: J, green: H, red: |
Open with DEXTER |
In this article, we present NIR spectra of a sample of candidate high-mass stars within NGG 7538 together with subarcsecond NIR photometry. We describe the observations and data reduction in Sect. 2. In Sect. 3 we analyse the spectro-photometric information of a sample of candidate high-mass stars. Section 4 is dedicated to estimating the distance, age, and mass of the cluster powering the H II region. Finally, we conclude in Sect. 5 with a brief summary.
2 Observations and data reduction
2.1 NIR photometry
We conducted
imaging observations of the field around NGC 7538 with LIRIS (Manchado et al. 2004; Acosta Pulido et al. 2003),
mounted at the Cassegrain focus of the William Herschel Telescope in La
Palma. Observations were obtained on 2006 July 21 as part of the
MASGOMAS programme that comprised a total of 45 young clusters in
the Milky Way (see Marín-Franch et al. 2009, for a description of the programme). The average seeing during the imaging observations of this region was 0
60, 0
63, and 0
64 for J, H, and
,
respectively. We followed the procedure described in Marín-Franch et al. (2009)
for the reduction and PSF extraction of the point-like sources in the
field. The saturation limit for the LIRIS photometry corresponds to
mag.
For brighter objects, we complemented the photometry with the 2MASS
point-source catalogue values that had their photometric errors
determined.
2.2 NIR spectroscopy
Follow-up spectroscopic observations of NGC 7538 were obtained
in 2007 September 20 and 2008 August 19 using the multi-object
spectroscopic (MOS) mode of the LIRIS Instrument. The average seeing
during both nights was 0
9. The designed mask comprised 20 slitlets of 0
8 in width and 6
in length. We used the low-resolution grism (
/
700) that covers the H and K band simultaneously and the medium-resolution grisms in the H and K bands with
/
2500 and
/
1700, respectively. The slitlet mask was oriented along a PA= 142
and centred at
13
39
,
.
We obtained the observations following a classic ABBA telescope nodding pattern with an offset of 2
5.
In the case of two slitlets, the spectra appeared contaminated by other
sources that fell close to one of the nodding positions
(position B). Therefore, we extracted those specific spectra on
only half of the stacked frames (at the position A). Observations of
the AOV standard star HD223386 were obtained in the different
configurations in order to remove any telluric features due to the
Earth's atmosphere.
The data were reduced using LIRISDR,
a specific data reduction package developed under the IRAF environment
by the LIRIS team. The routines available for the LIRIS MOS mode make
use of the a-priori mask design information to trace the slitlets and
their limits in the spatial direction in a uniformly illuminated frame.
Special care was taken to make the flat-field correction. The limited
spatial extent of each slitlet may include detector defects at
comparable or larger spatial scales, which are not detected within
individual flat-field slitlets. We circumvented the problem by coadding
the spectral response information of all apertures into a unique
monodimensional response, after shifting by an appropriate offset along
the spectral axis. This response was
modelled using a high order spline function. The pixel-to-pixel
correction was determined for each slitlet from the extracted
flat-field spectrum divided by the modelled response.
![]() |
Figure 2:
(Left panel) CMD
of the stellar content within
the field of view of the left panel. The red squares
represent the main sequence (MS) of dwarf stars without extinction at a
heliocentric distance of 2.7 kpc. The red discontinuous line
represents the reddening vector for an O9V star, with intervals of
AV = 5 mag indicated by crosses. Dark squares represent objects located at a projected distance less than 30
|
Open with DEXTER |
We performed a standard subtraction of consecutive AB pairs to remove the sky background. Individual flatfielding, retracing, resizing, and extraction for each star in every subtracted frame followed. Argon and Xenon arc frames were used to determine the wavelength calibration (rms < 1 Å). The resulting spectra were coadded weighting by the relative number of counts to provide a final spectrum of each star.
The wavelength calibrated spectra were divided by a normalised reference spectrum to remove the telluric absorption features. This normalised reference spectrum was generated from the observed spectra of the A0V standard star, divided by a Vega model spectrum convolved to the correspondent spectral resolution. For this, we constructed an IDL routine that iteratively fits the components due to the instrument response, an atmosphere's transmissivity template, and the Vega model spectrum.
Source #4 was observed with the LIRIS long-slit mode in the ranges 1.92-2.4 m and 1.52-1.78
m
in the nights of 2007 September 21 and 2008 June 27, respectively. The
intermediate-resolution grisms provided spectral resolutions of
/
and
/
.
In this case, we traced the spectrum using the continuum of the source
and applied an optimal weight to extract the flux along the spatial
direction. The standard star HD223386 was observed in the same
configurations to obtain the telluric correction.
3 Results
3.1 Colour-magnitude and colour-colour diagrams
The left and right panels of Fig. 2
show the colour-magnitude diagram (CMD) and the colour-colour diagram
(CCD) of the field around NGC 7538. In the figure, we have only
depicted those objects that have LIRIS counterparts in the H and
bands and in the three channels for the CMD and CCD, respectively.
The 2MASS source IRS 4 appears to be resolved into two components in our LIRIS images (see also Ojha et al. 2004). The brighter source to the southeast is saturated in our
map, while the northwest source #2 is one of the brightest in the
field. Therefore, we used our LIRIS photometry for source #2, while
still displaying the 2MASS integrated photometry of these two
point-like sources in our CMD (IRS 4). IRS 4 is not detected
in the 2MASS J band survey, therefore, it is important to notice that it shows
and
,
pointing to a strong NIR excess.
This paper aims at identifying the high-mass stellar content of a Galactic massive cluster, thus we selected the candidates for spectroscopic follow-up using a NIR photometric criterion. We looked for reddened bright objects in high stellar density regions. The production design specifications of the slitlet mask hinder a complete study of extremely crowded regions. A geometrical limitation is also imposed by the LIRIS multi-object mode. The spectral coverage for each object depends on the position of the corresponding slit relative to the centre of the chip in the dispersion direction. Therefore, the follow-up sources must be selected in a preferred orientation along the sky. Our CCD brings out the strong variations in extinction toward the sources selected for follow-up spectroscopy.
![]() |
Figure 3:
(Left panel) H-band normalised spectra of the 20 sources with MOS (
|
Open with DEXTER |
3.2 Spectral classification
We present the higher resolution spectra of the selected high-mass star candidates in Figs. 3 and 4. The lower resolution spectra are presented in Fig. 5. The inspection of the NIR spectra of the 21 objects observed in NGC 7538 suggests a classification into the six major groups exemplified in Fig. 6.
![]() |
Figure 4:
(Left panel) H-band normalised spectrum of source #4 obtained with long-slit (
|
Open with DEXTER |
Early-type stars:
only 3 objects within our candidate sample exhibit atmospheric spectral features in absorption, but none of them show traces of helium. This proves that sources #35, #25, and #17 are B-type stars. Considering the measured Br-11 and Br
IRS 6 and IRS 5 are very bright sources in the field of NGC 7538 that we did not include in our NIR spectroscopic follow-up. These objects were already classified by Russeil et al. (2007) according to the spectral features detected in their blue spectra. However, the strength of the HeI and HeII lines on these published spectra indicate that, opposite to what the authors state in the article, IRS 6 is of earlier spectral type than IRS 5. To investigate this, we performed spectroscopic observations in September 2009 with ISIS at the William Herschel Telescope, confirming that the spectral classification of Russeil et al. (2007) is swapped for these two particular objects. This revised classification solves the disparity in the estimates of the cluster distance derived from each star, as discussed later in Sect. 4.1.
H2 emission sources:
NIR ro-vibrational emission from H2 is detected toward source #63. This source is located in the vicinity of the active IRS 1-3 region and amidst the strong extended H2 emission reported by Davis et al. (1998) and Bloomer et al. (1998) and attributed to IRS 1/2 (Kraus et al. 2006). The molecular hydrogen emission indeed has a thermal origin, pointing to the presence of shocks produced by outflows. The spectra of sources #60, #11, and #2 also show traces of H2 1-0S(1) emission. However, this spectral feature has a nebular origin for the earlier case or it is not a dominant feature for the other two.CO bandhead emission sources:
objects #6 and #14 show CO first-overtone emission in their K-band spectra. These CO bandheads are emitted by neutral material with temperatures in the range 2000-5000 K and densities of




Featureless sources:
sources #23, #36, and #18 are identified as featureless. The NIR spectra of these objects could resemble those of low-mass Class I objects in Greene & Lada (1996) or even the Herbig Ae stars found by Rodgers (2001).Brackett+Fe emission sources:
two objects (#26 and #12) in the field show a peculiar spectrum with the Brackett series in emission and several FeII lines. Source #12 shows broad Br15-10 features (average


![]() |
Figure 5:
H- and K-band normalised spectra of the 20 sources with MOS (
|
Open with DEXTER |








There is a strong similarity between the H- and K-band spectra of #26 and #12 with some LBVs at low excitation (e.g. Geballe et al. 2000), showing Brackett lines and permitted and forbidden FeII lines in emission.
However, the low luminosity (
)
of these objects, together with the absence of MgII and SiII lines,
goes against any LBV classification. This should be regarded as a
warning against misidentifications in LBV infrared spectroscopic
searches in our Galaxy.
CO bandhead absorption sources:
seven objects within the sample show typical spectral features of cool stars. Our classification of these sources follows the prescriptions of Winston et al. (2009) and Bik et al. (2010) for late-type stars. This comparative method uses primarily the MgI, NaI, and CaI atomic lines as diagnostics for the temperature and the CO first overtone absorption bands as luminosity indicator. This classification carries an uncertainty of about
Despite the small offsets between the nodding positions and the narrow aperture used for extracting the spectra, nebular H II emission appears strongly in the spectra of sources #58 and #41. These objects are located in the line of sight of areas where a strong gradient in the nebular emission is detected in our NIR maps. Therefore, these lines should be considered with caution.
A summary of the photometric and spectroscopic characterisation of NGC 7538's stellar content is compiled in Table 1.
![]() |
Figure 6: Sample of representative spectra in NGC 7538. B-star: B-type star; H2: dominant H2 lines; YSO: young stellar object; FL: featureless spectrum; Br+Fe: Brackett and FeII emission lines; Cool-St: cool star. |
Open with DEXTER |
The spectral type classification allows the estimate of the extinction toward several sources in NGC 7538. The visual extinctions that appear in Col. 8 of Table 1 were calculated from the spectral type classification and corresponding stellar colours, considering the Rieke & Lebofsky (1985) extinction law. The intrinsic colours of O-type dwarfs were taken from Martins & Plez (2006), while those for B dwarfs and giants were extracted from Tokunaga (2000). For PMS stars, we considered the intrinsic colours reported by Kenyon & Hartmann (1995). Absolute magnitudes are taken from Martins & Plez (2006) for O stars.
3.3 Infrared excess
We have complemented our NIR photometry with IRAC/Spitzer mid-infrared (MIR) data for this region (Chavarria et al. private communication) in order to study the infrared excess toward our follow-up sources. Figure 7 shows the spectral energy distribution (SED) of the spectroscopically identified young sources within NGC 7538.
We fitted the LIRIS and IRAC flux distributions with (i) a Kurucz model depicted with a solid line for objects with a spectral type classification of according temperature and surface gravity; (ii) a single black-body distribution represented by a broken line for those sources whose spectral type could not be constrained. Only in the first case were the fluxes dereddened following the Fitzpatrick (1999) parametrisation. The total luminosity for each object was computed by integrating the area below the respective model fit, and it is considered a lower limit for objects that display any excess.
Object #63 shows both NIR and MIR excess, although the 8.0 m
IRAC map appears saturated at this location. The large IRAC colour
[3.6]-[4.5] = 0.9 mag hints at a Class I source. Martins et al. (2009) show that the H2-dominated sources in the RCW 79 and RCW 120 star-forming regions have bright counterparts at 24.0
m. The authors conclude that this large infrared excess comes from the existence of an envelope, implying that these H2 sources are in an early evolutionary phase.
However, the attribute of massive YSO relies on the additional presence of a strong IR excess. #6, #14 and #11 exhibit both NIR and MIR excesses. Owing to its proximity to the bright source IRS 4 and to saturation in the IRAC maps, we can only report the NIR excess of source #2 at subarcsecond spatial resolution. The photospheric features present in sources #6 and #14 seem indicative of a more evolved phase, when the central star has significantly eroded the disk and photospheric features in absorption are present. This idea is supported by their IRAC colour-colour classification: Class II in the case of #6 and Class I for #14 and #11. However, the fits of the SEDs to black-body distributions indicates decreasing temperature, which could be directly related to the different spectral features.
Objects #23, #36, and #18 possess featureless NIR spectra, and
#23 and #18 show both NIR and MIR excess and a Class II IRAC
classification, indicating a very young nature (Greene & Lada 1997). Source #36 is a deeply embedded source (
mag)
that displays the NIR colours
of a reddened naked photosphere without any intrinsic excess. This
hypothesis is supported by its not being detected in the IRAC maps. We
tentatively classify #36 as an early-type O star whose shallow hydrogen
features are not detected at the limited SNR of our NIR spectra.
The SED of #26 (IRS 3) reveals a strong MIR excess that translates into the IRAC colours of a Class I object.
We have spectroscopically classified source #12 as a Herbig Be star. The SED reveals an IR excess that rises up to 8.0 m. Its NIR photometry renders a reddening-free parameter
equivalent to -0.29. Herbig Ae/Be star candidates cover a range -
1.38 < Q < - 0.22 according to Hernández et al. (2005),
reinforcing our classification. Among the identified PMS stars, only
source #58 shows significant mid-IR excess. However, the non detection
of #55 and #39 in the IRAC maps at 5.8 and 8.0
m prevents further conclusions on the IR excess of these sources.
Despite being classified as OB stars, the SEDs of sources IRS 5
and #35 exhibit strong MIR excesses (i.e. the typical colours of
Class I sources). In the case of #35, this effect could be
explained by the presence of a nearby object (1
2 away) that is not resolved in the longer wavelength IRAC channels.
Table 1: Coordinates, LIRIS photometry and spectral classification of the selected objects within NGC 7538.
![]() |
Figure 7: SEDs of young sources in
NGC 7538. Black dots with error bars are photometric measurements
with LIRIS and IRAC. For the sources with a spectral type
classification and an extinction estimate, the SED has been dereddened.
Solid grey lines represent the Kurucz model fit of the respective
spectral type classification, indicated by the
|
Open with DEXTER |
4 Discussion
4.1 Distance, age, and mass
An initial estimate of the distance based on spectro-photometry by Blitz et al. (1982) situated NGC 7538 at an heliocentric distance of 2.8 0.9 kpc.
Meanwhile, kinematical estimates delivered a distance for the same region as far as 3.5 kpc (Israel et al. 1973).
Russeil et al. (2007) determined the spectral type classification of two stars in the H II region, IRS 6 (O9V), and IRS 5 (O3V), inferring spectro-photometric distances of 1.6
0.17 and 4.24
0.29 kpc, respectively. They estimated an average distance of 2.27
0.15 kpc
for the entire Sh 2-158 region.
The most recent estimate of the distance to NGC 7538 has been
obtained by measuring trigonometric parallaxes of methanol masers,
regularly associated with high-mass star forming regions. This study
yielded a distance of
2.65 +0.12-0.11 kpc Moscadelli et al. (2009).
Our revised optical classification of the O-type sources IRS 6 and IRS 5 described in Sect. 3.2
renders a new distance estimate to NGC 7538. Considering a
classification O3V for the IRS 6 object, we derive a
spectro-photometric distance of 2.99 0.5 kpc to the cluster powering the H II region. Likewise, the classification of IRS 5 yields a value of 2.39
0.4 kpc. We report a final spectro-photometric distance of 2.7
0.5 kpc to this region.
Despite the large extinction variations in the region that hampered robust age determination, Balog et al. (2004) find an older generation of stars (4 Myr) and a younger population of faint stars closer to
1 Myr. A fraction of 30% of the young stars in the region around NGC 7538 was reported to exhibit a near-IR excess.
An upper limit limit to the age of the cluster can be established from the most massive classified star member. According to Schaerer & de Koter (1997), a 60
star (corresponding to a O3V star) has a main-sequence lifetime
2.2 Myr
before it starts the giant phase. Assuming a coeval star-formation
event, the oldest stellar content in NGC 7538 must therefore be
younger than 2.2 Myr.
At the other end of the mass spectrum, we can use the information
derived from the PMS stars to constrain the age of the low-mass
population around NGC 7538. To construct a Hertzsprung-Russell
diagram (HRD), we used the spectral type to temperature conversion from
Kenyon & Hartmann (1995), including the overestimates reported by Cohen & Kuhi (1979) as temperature error. The absolute K
magnitudes were then calculated from the dereddened apparent brightness
and assuming a distance of 2.7 kpc. The HRD of the PMS stars
identified in NGC 7538 is shown in Fig. 8. We depict several pre-main-sequence evolutionary tracks in the figure and isochrones taken from Da Rio et al. (2009). The mass of these PMS stars is estimated in the range 2-4
.
The comparison with the PMS isochrones yields an age range between 0.5 and 2 Myr for the PMS identified population.
We thus conclude that these two estimates result in an age range between 0.5 and 2.2 Myr for the powering cluster.
Clearly, other massive young stellar populations are also present at several locations around NGC 7538. They range from stars for which we observe photospheric spectral features and which are detectable in the optical, but also exhibit very strong IR excess (IRS 5 at region (ii)) to objects whose emission is dominated by a circumstellar disk (region (iii)), and finally to nascent stars that present ongoing outflow activity (IRS 1/2 and IRS 9 in regions (iv) and (vi), respectively). These results support the previous idea of a star-forming sequence in the northwest-southeast direction that may be smoother than previously considered. The quantitative analysis of the individual ages of these populations is beyond the scope of this paper.
![]() |
Figure 8: Extinction corrected Hertzsprung-Russell diagram of the identified PMS stars within NGC 7538. The open symbols correspond to PMS stars without an IR excess detected in their LIRIS and IRAC photometry, while the filled ones represent those with detected IR excess. The solid lines represent the pre-main-sequence isochrones, while the dotted lines correspond to evolutionary tracks by Da Rio et al. (2009). |
Open with DEXTER |
By assuming a Salpeter initial mass function (IMF) and extrapolating it down to 0.8
,
we can calculate the stellar mass of a cluster. Nevertheless, this
estimate is only meaningful for a coeval stellar population. This
should be the case for the powering cluster of the H II
region, with IRS 6 as the ionising star. After normalising the
stellar mass distribution by the detection of one O3V star with a mass
uncertainty between 47 and 64
,
we obtain a total mass
for the cluster.
4.2 Spatial distribution
The second most massive, identified star is IRS 5, and although it is located amidst the ionised emission, its infrared excess and proximity to the rim of the molecular cloud suggest that it may not be a member of the central cluster. Region (ii) in fact comprises a younger generation of stars. Among the spectroscopically identified B-type stars, only object #25 is located amidst the ionised emission, while sources #17 and #35 are farther away. In the case of #35, this source is located in the vicinity of the young region IRS 1-3 or (iv) and exhibits infrared excess typical of Class 0/I objects.
In fact, our CMD in the left panel of Fig. 2 reveals a deficit of early B-type candidates for sources located within a radius of 30
from
the O3V ionising star IRS 6 (objects marked by squares). To
further study this sparseness of early B-type stars associated to the
cluster, we identify B-type candidates of compatible
-band brightness and
mag in the CMD. These objects are indicated by diamonds in the left panel of Fig. 2.
Analogously to our previous analysis of spectroscopically identified
B-type stars, a few of the B-type candidates show a mid-IR excess
(Chavarría et al., private communication) suggesting they belong
to younger pockets of star formation. Other B-type candidates exhibit
only photospheric emission and are located to the east and south of
IRS 6, spreading beyond the boundary of the H II region, possibly due to a mass segregation effect.
We sketch a possible picture in which the cavity created by the ionizing stars in the molecular cloud is open in the observer's direction, and part of its stellar population is located in the foreground of the neutral molecular cloud, appearing almost aligned with other embedded younger star-forming regions. This hypothesis is supported by the work of Balog et al. (2004), who find that the stellar density peak of this region at NIR wavelengths is located at the rim of the H II region.
This projection effect is particularly important to bear in mind when isolating the stellar populations of the different pockets of star-formation that may be triggered by the expansion of H II regions.
5 Summary and conclusions
We have reported the spectro-photometric study of a sample of candidate
high-mass stars in NGC 7538. We used the LIRIS instrument to
obtain
photometry, and multi-object and long-slit spectroscopy in the H and K
bands at low and medium resolution. Complementary IRAC photometry of
the sample has aided in the investigation of infrared excesses.
Besides the optically known O stars IRS 6 and IRS 5, whose spectral type classification we revised, we pinpointed a population of B-type stars located either at or beyond the rim of the H II region. In view of this, we confirm that IRS 6 is the main ionising source of Sh 2-158. Considering the non detection of spectral features and the lack of infrared excess, we suspect that object #36 is a deeply embedded, early O-type star in the vicinity of IRS 1-3. However, higher signal-to-noise-ratio observations are needed to ascertain this hypothesis. We spectroscopically identified four G and K stars still contracting to the main sequence, and therefore labelled them as PMS. These stars are considered to be the precursors of A and F stars.
The observations of IRS 3, aka #26, are consistent with its
previous classification as a B-type candidate that has already started
powering a compact H II region in a high-density environment
cm-3 near the protostars IRS 1-2. The outflow activity traced by the [FeII] and H2 molecular emission also support the young nature of this object.
Five of the observed sources have been classified as YSOs, with their spectra dominated either by H2 or CO features. Two stars display featureless NIR spectra and a prominent infrared excess, while object #12 has been identified as a Herbig Be star. In general, all these YSOs are considered to be the precursors of OB stars.
We established a distinction between an older population of naked stellar photospheres and various regions of more recent star formation around NGC 7538. However, the former population spreads spatially beyond the H II region likely because of a projection effect, which makes it difficult to isolate the populations associated to the individual pockets of star formation. We find that most of the YSOs are concentrated either at the southern rim of the H II region or in the area around the infrared sources IRS 1-3. This, together with the known ongoing formation of the protostars in the field, reinforces the idea of a northwest-southeast star-forming sequence. However, the study of the detailed sequence of events and causality between them cannot be addressed with the data presented.
We have estimated a distance of 2.7 0.5 kpc, an age range 0.5-2.2 Myr, and a total mass
for the evolved cluster content within NGC 7538. The variety of
NIR spectral features shown by young objects in this field makes
NGC 7538 an ideal region for studying NIR spectral features of YSO
for a wide range of masses and evolutionary phases in a determined
environment.
We thank the referee for valuable comments that helped improve the contents of this paper. We would also like to thank Dr. Belén López for her input into the initial manuscript.E.P is funded by the Spanish MICINN under the Consolider-Ingenio 2010 Program grant CSD2006-00070: First Science with the GTC (http://www.iac.es/consolider-ingenio-gtc). This work has been partly supported by grants AYA2007-67456-C02/AYA2008-06166-C03.
This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.
References
- Acosta Pulido, J. A., Ballesteros, E., Barreto, M., et al. 2003, The Newsletter of the Isaac Newton Group of Telescopes, 7, 15 [Google Scholar]
- Balog, Z., Kenyon, S. J., Lada, E. A., et al. 2004, AJ, 128, 2942 [NASA ADS] [CrossRef] [Google Scholar]
- Benjamin, R. A., Skillman, E. D., & Smits, D. P. 1999, ApJ, 514, 307 [NASA ADS] [CrossRef] [Google Scholar]
- Bik, A., Kaper, L., & Waters, L. B. F. M. 2006, A&A, 455, 561 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bik, A., Puga, E., Waters, L. B. F. M., et al. 2010, ApJ, 713, 833 [Google Scholar]
- Blitz, L., Fich, M., & Stark, A. A. 1982, ApJS, 49, 183 [NASA ADS] [CrossRef] [Google Scholar]
- Bloomer, J. D., Watson, D. M., Pipher, J. L., et al. 1998, ApJ, 506, 727 [NASA ADS] [CrossRef] [Google Scholar]
- Campbell, B., & Persson, S. E. 1988, AJ, 95, 1185 [NASA ADS] [CrossRef] [Google Scholar]
- Clark, J. S., & Steele, I. A. 2000, A&AS, 141, 65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cohen, M., & Kuhi, L. V. 1979, ApJS, 41, 743 [NASA ADS] [CrossRef] [Google Scholar]
- Da Rio, N., Gouliermis, D. A., & Henning, T. 2009, ApJ, 696, 528 [NASA ADS] [CrossRef] [Google Scholar]
- Davies, B., Origlia, L., Kudritzki, R., et al. 2009, ApJ, 696, 2014 [NASA ADS] [CrossRef] [Google Scholar]
- Davis, C. J., Moriarty-Schieven, G., Eislöffel, J., Hoare, M. G., & Ray, T. P. 1998, AJ, 115, 1118 [NASA ADS] [CrossRef] [Google Scholar]
- Fitzpatrick, E. L. 1999, PASP, 111, 63 [NASA ADS] [CrossRef] [Google Scholar]
- Geballe, T. R., Najarro, F., & Figer, D. F. 2000, ApJ, 530, L97 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Greene, T. P., & Lada, C. J. 1996, AJ, 112, 2184 [NASA ADS] [CrossRef] [Google Scholar]
- Greene, T. P., & Lada, C. J. 1997, AJ, 114, 2157 [NASA ADS] [CrossRef] [Google Scholar]
- Hamann, F., Simon, M., Carr, J. S., et al. 1994, ApJ, 436, 292 [NASA ADS] [CrossRef] [Google Scholar]
- Hanson, M. M., Conti, P. S., & Rieke, M. J. 1996, ApJS, 107, 281 [Google Scholar]
- Hanson, M. M., Rieke, G. H., & Luhman, K. L. 1998, AJ, 116, 1915 [NASA ADS] [CrossRef] [Google Scholar]
- Hernández, J., Calvet, N., Hartmann, L., et al. 2005, AJ, 129, 856 [NASA ADS] [CrossRef] [Google Scholar]
- Israel, F. P., Habing, H. J., & de Jong, T. 1973, A&A, 27, 143 [NASA ADS] [Google Scholar]
- Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 [NASA ADS] [CrossRef] [Google Scholar]
- Knödlseder, J. 2000, A&A, 360, 539 [NASA ADS] [Google Scholar]
- Kraus, S., Balega, Y., Elitzur, M., et al. 2006, A&A, 455, 521 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Manchado, A., Barreto, M., Acosta-Pulido, J., et al. 2004, in SPIE Conf. Ser., 5492, 1094 [Google Scholar]
- Marín-Franch, A., Herrero, A., Lenorzer, A., et al. 2009, A&A, 502, 559 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martins, F., & Plez, B. 2006, A&A, 457, 637 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martins, F., Pomares, M., Deharveng, L., Zavagno, A., & Bouret, J. 2010, A&A, 510, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- McCaughrean, M., Rayner, J., & Zinnecker, H. 1991, Mem. Soc. Astron. Ital., 62, 715 [Google Scholar]
- Messineo, M., Davies, B., Ivanov, V. D., et al. 2009, ApJ, 697, 701 [NASA ADS] [CrossRef] [Google Scholar]
- Meyer, M. R., Calvet, N., & Hillenbrand, L. A. 1997, AJ, 114, 288 [NASA ADS] [CrossRef] [Google Scholar]
- Moreno, M. A., & Chavarria-K., C. 1986, A&A, 161, 130 [NASA ADS] [Google Scholar]
- Moscadelli, L., Reid, M. J., Menten, K. M., et al. 2009, ApJ, 693, 406 [NASA ADS] [CrossRef] [Google Scholar]
- Ojha, D. K., Tamura, M., Nakajima, Y., et al. 2004, ApJ, 616, 1042 [NASA ADS] [CrossRef] [Google Scholar]
- Pestalozzi, M. R., Elitzur, M., Conway, J. E., et al. 2004, ApJ, 603, L113 [NASA ADS] [CrossRef] [Google Scholar]
- Pestalozzi, M. R., Minier, V., Motte, F., et al. 2006, A&A, 448, L57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Puga, E., Hony, S., Neiner, C., et al. 2009, A&A, 503, 107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618 [NASA ADS] [CrossRef] [Google Scholar]
- Rodgers, B. M. 2001, Ph.D. Thesis, University of Washington [Google Scholar]
- Russeil, D., Adami, C., & Georgelin, Y. M. 2007, A&A, 470, 161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schaerer, D., & de Koter, A. 1997, A&A, 322, 598 [NASA ADS] [Google Scholar]
- Scoville, N. Z., Hall, D. N. B., Ridgway, S. T., et al. 1979, ApJ, 232, L121 [NASA ADS] [CrossRef] [Google Scholar]
- Sewilo, M., Churchwell, E., Kurtz, S., Goss, W. M., & Hofner, P. 2004, ApJ, 605, 285 [NASA ADS] [CrossRef] [Google Scholar]
- Tokunaga, A. T. 2000, Infrared Astronomy, ed. Cox, A. N., 143 [Google Scholar]
- Winston, E., Megeath, S. T., Wolk, S. J., et al. 2009, AJ, 137, 4777 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
- ... NGC 7538
- Based on observations made with the WHT operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias.
- ... LIRISDR
- http://www.iac.es/project/LIRIS
All Tables
Table 1: Coordinates, LIRIS photometry and spectral classification of the selected objects within NGC 7538.
All Figures
![]() |
Figure 1:
Colour composite image of NGC 7538 obtained with LIRIS ( blue: J, green: H, red: |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
(Left panel) CMD
of the stellar content within
the field of view of the left panel. The red squares
represent the main sequence (MS) of dwarf stars without extinction at a
heliocentric distance of 2.7 kpc. The red discontinuous line
represents the reddening vector for an O9V star, with intervals of
AV = 5 mag indicated by crosses. Dark squares represent objects located at a projected distance less than 30
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
(Left panel) H-band normalised spectra of the 20 sources with MOS (
|
Open with DEXTER | |
In the text |
![]() |
Figure 4:
(Left panel) H-band normalised spectrum of source #4 obtained with long-slit (
|
Open with DEXTER | |
In the text |
![]() |
Figure 5:
H- and K-band normalised spectra of the 20 sources with MOS (
|
Open with DEXTER | |
In the text |
![]() |
Figure 6: Sample of representative spectra in NGC 7538. B-star: B-type star; H2: dominant H2 lines; YSO: young stellar object; FL: featureless spectrum; Br+Fe: Brackett and FeII emission lines; Cool-St: cool star. |
Open with DEXTER | |
In the text |
![]() |
Figure 7: SEDs of young sources in
NGC 7538. Black dots with error bars are photometric measurements
with LIRIS and IRAC. For the sources with a spectral type
classification and an extinction estimate, the SED has been dereddened.
Solid grey lines represent the Kurucz model fit of the respective
spectral type classification, indicated by the
|
Open with DEXTER | |
In the text |
![]() |
Figure 8: Extinction corrected Hertzsprung-Russell diagram of the identified PMS stars within NGC 7538. The open symbols correspond to PMS stars without an IR excess detected in their LIRIS and IRAC photometry, while the filled ones represent those with detected IR excess. The solid lines represent the pre-main-sequence isochrones, while the dotted lines correspond to evolutionary tracks by Da Rio et al. (2009). |
Open with DEXTER | |
In the text |
Copyright ESO 2010
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.