C. J. Evans1 -
S. J. Smartt2 -
J.-K. Lee2 -
D. J. Lennon1 -
A. Kaufer3 -
P. L. Dufton2 -
C. Trundle4 -
A. Herrero4,5 -
S. Sim
n-Díaz4 -
A. de Koter6 -
W.-R. Hamann7 -
M. A. Hendry8 -
I. Hunter2 -
M. J. Irwin8 -
A. J. Korn9 -
R.-P. Kudritzki10 -
N. Langer11 -
M. R. Mokiem6 -
F. Najarro12 -
A. W. A. Pauldrach13 -
N. Przybilla14 -
J. Puls13 -
R. S. I. Ryans2 -
M. A. Urbaneja10 -
K. A. Venn15 -
M. R. Villamariz4
1 - The Isaac Newton Group of Telescopes, Apartado de Correos 321, 38700
Santa Cruz de La Palma, Canary Islands, Spain
2 -
The Department of Pure and Applied Physics, The Queen's University of Belfast,
Belfast BT7 1NN, Northern Ireland, UK
3 -
European Southern Observatory, Alonso de Cordova 3107, Santiago 19, Chile
4 -
Instituto de Astrofísica de Canarias, 38200 La Laguna, Tenerife, Spain
5 -
Dept. de Astrofísica, Universidad de La Laguna, Avda. Astrofísico Francisco S
nchez, s/n, 38071 La Laguna, Spain
6 -
Astronomical Institute Anton Pannekoek, University of Amsterdam,
Kruislaan 403, 1098 SJ Amsterdam, The Netherlands
7 -
Lehrstuhl Astrophysik der Universit
t Potsdam, Am Neuen Palais 10,
14469 Potsdam, Germany
8 -
Institute of Astronomy, University of Cambridge,
Madingley Road, Cambridge, CB3 0HA, UK
9 -
Uppsala Astronomical Observatory, Box 515, 75120 Uppsala, Sweden
10 -
Institute for Astronomy, 2680 Woodlawn Drive,
Honolulu, HI 96822, USA
11 -
Astronomical Institute, Utrecht University, Princetonplein 5,
3584 CC Utrecht, The Netherlands
12 -
Instituto de Estructura de la Materia, CSIC,
C/ Serrano 121, 28006 Madrid, Spain
13 -
Universit
ts-Sternwarte M
nchen,
Scheinerstr. 1, 81679 M
nchen, Germany
14 -
Dr. Remeis-Sternwarte Bamberg, Sternwartstr. 7, 96049 Bamberg, Germany
15 -
Macalester College, 1600 Grand Avenue, Saint Paul, MN 55105, USA
Received 29 November 2004 / Accepted 14 March 2005
Abstract
We introduce a new survey of massive stars in the Galaxy and the
Magellanic Clouds using the Fibre Large Array Multi-Element
Spectrograph (FLAMES) instrument at the Very Large Telescope (VLT).
Here we present observations of 269 Galactic stars with the
FLAMES-Giraffe Spectrograph (
), in fields centered on
the open clusters NGC 3293, NGC 4755 and NGC 6611. These data are
supplemented by a further 50 targets observed with the Fibre-Fed
Extended Range Optical Spectrograph (FEROS, R = 48 000). Following
a description of our scientific motivations and target selection
criteria, the data reduction methods are described; of critical importance
the FLAMES reduction pipeline is found to yield spectra that are in
excellent agreement with less automated methods. Spectral
classifications and radial velocity measurements are presented for
each star, with particular attention paid to morphological
peculiarities and evidence of binarity. These observations represent
a significant increase in the known spectral content of NGC 3293
and NGC 4755, and will serve as standards against which our subsequent
FLAMES observations in the Magellanic Clouds will be compared.
Key words: stars: early-type - stars: fundamental parameters - open clusters and associations: individual: NGC 3293 - open clusters and associations: individual: NGC 4755 - open clusters and associations: individual: NGC 6611
The Fibre Large Array Multi-Element Spectrograph (FLAMES) at the Very Large Telescope (VLT) provides multi-object spectroscopy over a corrected 25' diameter field-of-view. As part of a European Southern Observatory (ESO) Large Programme, comprising over 100 h of VLT time, we have observed in excess of 50 O-type stars and 500 B-type stars, in a total of seven clusters distributed over the Galaxy and the Magellanic Clouds. This unprecedented survey of massive stars, and its homogeneous spectroscopic analysis using state-of-the-art model atmospheres, will be the subject of a series of papers. This first paper gives an overview of the scientific goals of the project, discusses the data reduction techniques employed for the survey, and presents the data for the three Galactic clusters.
A diverse range of topics depend critically on our understanding of
the evolution of massive stars. Through their stellar winds and as
supernova explosions they are the main source of kinetic energy into
the interstellar medium of their host galaxies, and they are
responsible for the continuous build-up of elements heavier than
helium ("metals'') from successive generations of star formation
(e.g. Shapley et al. 2004). Indeed, one of the key challenges is to
fully understand the role of metallicity on stellar evolution, to
place better constraints on extremely low-metallicity stellar models
(e.g. Kudritzki 2002). Theoretical simulations of star formation in
the early Universe have suggested that masses larger than
100
are strongly favoured (Abel et al. 2002; Bromm et al. 2002), and
these "zero-metallicity'' (or Population III) stars are considered to
be prime candidates for the re-ionization of the Universe
(e.g. Haehnelt et al. 2001; Wyithe & Loeb 2003) that may have happened at
redshifts greater than six (Becker et al. 2001).
Furthermore, with ongoing studies of the nature of supernova progenitors
(e.g. Smartt et al. 2004), and following recent discoveries that some of
the most spectacular supernovae appear to be related to
-ray bursts
(e.g. Galama et al. 1998; Hjorth et al. 2003; Bloom et al. 2002), understanding the evolution of
massive stars will gain even greater importance in the near future.
However, fundamental problems remain as to how the evolution of
massive stars depends on the physics of mass-loss and rotation, and
the role of metallicity. The FLAMES survey of massive stars aims at a
comprehensive study of these processes. Some of the main scientific
issues that we will address with the project are:
In the current paper we discuss the methods used for target selection and detail the data reduction techniques employed for the survey. We also present spectral classifications, radial velocities and cross-identifications with previous catalogues for our targets in the three Galactic clusters: NGC 3293, NGC 4755 (aka The Jewel Box), and NGC 6611 (aka The Eagle Nebula; Messier 16). The data presented here also include observations of the brightest cluster members with the Fibre-Fed, Extended Range Optical Spectrograph (FEROS) at the 2.2-m Max Planck Gesellschaft (MPG)/ESO telescope. A thorough analysis of the rotational velocities of the Galactic sample will be given elsewhere (Lee et al. in preparation) and similar catalogues for the Magellanic Cloud targets will be presented in a further publication (Evans et al. in preparation).
Table 1: Galactic and Magellanic Cloud clusters observed with FLAMES. The earliest main-sequence object observed in each cluster is given as an indication of their relative ages.
Table 2: Summary of the relevant EIS pre-FLAMES observations of Galactic fields used for target selection and astrometry.
Astrometry and photometry for our target fields were acquired from observations for the ESO Imaging Survey (EIS) that used the Wide Field Imager (WFI) at the 2.2-m MPG/ESO telescope (Momany et al. 2001). The relevant B (B/99, ESO#842) and V (V/89, ESO#843) science and calibration images were taken from the ESO archive and reduced using a modified version of the Isaac Newton Telescope-Wide Field Camera (INT-WFC) data reduction pipeline (Irwin & Lewis 2001). Multiple exposures of the fields (aimed at shallow and deeper magnitudes) were combined to produce the final photometry; the observations are summarised in Table 2.
Photometric transformations for the WFI data to the Johnson-Cousins
system were determined for each cluster from published photometry,
using visually matched stars - in the absence of observations of
standard fields this "bootstrap'' approach was both necessary and
sufficient for our requirements. The brightest targets in each
cluster were saturated, even in the short WFI exposures,
although the astrometry of these sources was more than adequate for
positioning of the FLAMES fibres. The photometric transformations were
determined using those targets that displayed a linear response i.e.,
V > 10.75
in NGC 3293 and V > 11.5
in
NGC 4755 and NGC 6611.
In Figs. 1 and 2 we show fits to the colour terms in V and B; published photometry was taken from Turner et al. (1980, NGC 3293), Sanner et al. (2001, NGC 4755), and
Hillenbrand et al. (1993, NGC 6611). The transformation equations found (from 51 stars) for
NGC 3293 were:
| (1) |
| (2) |
| (3) |
| (4) |
| (5) |
| (6) |
![]() |
Figure 1:
Comparison of
|
| Open with DEXTER | |
![]() |
Figure 2:
Comparison of
|
| Open with DEXTER | |
Following transformation to the standard system, the observed stars in each cluster were assigned running numbers on the basis of the V-band magnitudes, with "3293-001'' the brightest star observed in NGC 3293 and so on. This is not strictly the case for the stars in NGC 6611, following refinement of the calibrations (using an enlarged sample of published values) after the assignment of the running numbers.
In the process of determining the photometric transformations we noticed that
two stars in NGC 4755, namely IV-11 and IV-17 in the numbering system of
Arp & van Sant (1958), are listed by Sanner et al. (2001) as
having
14
.
These are significantly discrepant from
the (albeit photographic) values given by Arp & van Sant (V
= 11.35 and 9.76
respectively). Photographic (IV-11: V
= 11.41
)
and photoelectric (IV-17: V = 9.96
)
results from Dachs & Kaiser (1984) roughly tally with the Arp & van Sant values,
and the WFI observations confirm that they are much brighter than 14th mag. We suggest that these may be misidentifications by Sanner et al. (2001) rather than significant, intrinsic variability.
More curious is the case of Walker 442 in NGC 6611, with Thé et al. (1990)
reporting V = 8.26. The photographic observations by
Walker (1961, V = 15.09m) and the value from the
WFI data (V = 14.99
)
are both much fainter for this star.
Though the star is not explicitly cross-referenced, Hillenbrand et al. (1993) also
give V = 15.00 for a star at
18
19
0.86
,
-13
48'47.2''(J2000), almost certainly W442. It is hard to reconcile these
observations with the Thé et al. (1990) value which, if real, would make
W442 one of the visually brightest stars in the cluster. The sources
of this discrepancy are not clear, but we note that the SIMBAD
database reports the Thé et al. results.
The primary element of FLAMES is the Giraffe spectrograph. Our
programme employed the mode of Giraffe in which it is fed by 132 so-called "Medusa'' fibres. Input catalogues of possible targets were compiled for each cluster from the WFI photometry and astrometry; these were then used to allocate the Medusa fibres with the FLAMES Fibre Positioner Observation Support Software (FPOSS). The philosophy behind compilation of Giraffe targets was to select bright blue stars, with cuts of
(B - V)0 < 0.00
(corresponding to the threshold between B and A spectral types) and to a faint limit such that the
main sequence was sampled down to
B5 (corresponding to
).
An important consideration for application of these selection criteria
in the Galactic clusters is an estimate of the extinction and
reddening. Estimates of E(B - V) were incorporated into our
selection routines, adopting
E(B - V) = 0.30
for NGC 3293
(e.g., Turner et al. 1980) and 0.40
for NGC 4755 (e.g., de Waard et al. 1984).
The treatment of reddening in NGC 6611 was somewhat more complicated;
the extinction is
large and strongly variable
across the cluster and a value of
E(B - V) = 1.20
was
adopted, which is at the upper end of those found by Hillenbrand et al. (1993).
This value ensured inclusion of the majority of the early-type objects
(necessary for a reasonable density of targets across the FLAMES
field-of-view), but contamination by (less-reddened) later-type stars
can then become a problem. Indeed, as demonstrated by the spectral
types present in our final sample (see Table 11), selection of
solely early-type objects in such a cluster is extremely difficult
without previous spectroscopic information.
Following compilation of the input catalogues, FPOSS was used to assign the Medusa fibres. In addition to selection effects such as avoiding fibre collisions and the physical size of the fibre buttons, the final spectral samples are influenced by some additional external constraints. Where published spectral types were available (and the object was of interest) a higher weight was assigned to the target when using FPOSS. Similarly, known Be-type stars were (in general) deliberately omitted from the input catalogues used with FPOSS; this is primarily relevant to configuration of our Magellanic Cloud fields but is mentioned here to note that this survey is not suitable for statistical studies of the incidence of Be-type stars.
The brightest stars in each cluster were deliberately omitted from the FLAMES observations as they were too bright for the dynamic range of the instrument; such stars would simply saturate during the exposure times needed for the fainter B-type main-sequence stars. These stars were observed separately using FEROS at the 2.2-m MPG/ESO telescope (see Sect. 2.3).
The colour-magnitude diagram for all the stars in the FLAMES
field-of-view centered on NGC 3293 is shown in Fig. 3.
Magnitudes and colours for stars fainter than V = 10.75
are
from our WFI photometry (corrected to the Johnson system); with
photometry for the brighter stars from Turner et al. (1980) and
Feinstein & Marraco (1980). The faint limit of our targets (
13.25
)
is marked by the dotted line. The colour-magnitude diagrams for
NGC 4755 and 6611 are shown in Figs. 4 and 5, with the photometry for the brightest stars from Sanner et al. (2001) and Hillenbrand et al. (1993) respectively.
![]() |
Figure 3: Colour-magnitude diagram for all objects in the FLAMES field-of-view for NGC 3293. The marked targets were observed with FLAMES-Giraffe (solid circles) or with FEROS (solid triangles). |
| Open with DEXTER | |
![]() |
Figure 4: Colour-magnitude diagram for all objects in the FLAMES field-of-view for NGC 4755. The marked targets were observed with FLAMES-Giraffe (solid circles) or with FEROS (solid triangles). |
| Open with DEXTER | |
![]() |
Figure 5: Colour-magnitude diagram for all objects in the FLAMES field-of-view for NGC 6611. The marked targets were observed with FLAMES-Giraffe (solid circles) or with FEROS (solid triangles). |
| Open with DEXTER | |
The brightest targets in NGC 3293 and NGC 4755 were observed using FEROS at the ESO/MPG 2.2-m telescope on 2004 July 4-6. We also sought observations of remaining targets of interest in the NGC 6611 field and, to increase the number of Galactic O-type stars in the survey, we included three stars classified as O-type by Hillenbrand et al. (1993) that were not in the FLAMES field-of-view (6611-001, 6611-005, and 6611-045; W412, HD 168504, and W584 respectively).
FEROS is a fixed configuration instrument (with R = 48 000), giving a wide wavelength coverage of 3600-9200 Å in one exposure. The FEROS data were reduced using the reduction pipeline that runs under the MIDAS environment (Kaufer et al. 1999). In the majority of cases the signal-to-noise ratio of the FEROS data is in excess of 100 per resolution element. Targets observed with FEROS are marked in Figs. 3-5 as solid triangles.
Table 3: Summary of VLT-FLAMES observations of Galactic clusters. The exposure times are for each wavelength setting.
Each of the three Milky Way clusters was observed with one FLAMES
fibre configuration, centered on the cluster core. Six of the
standard high-resolution Giraffe settings were used: HR02 (with a
central wavelength of
3958 Å), HR03 (
4124), HR04 (
4297), HR05 (
4471), HR06 (
4656), and HR14 (
6515). The spectral range for each setting is approximately 200 Å, i.e. there is sufficient overlap between settings to give a continuous
spectrum of 3850-4755 Å, with additional coverage of
6380-6620 Å, primarily to observe the H
Balmer line. The
observations are summarized in Table 3, with details of
the individual Galactic targets given in Tables 9-11.
The Galactic FLAMES data were reduced using IRAF
(v2.11). Standard routines were used for bias
subtraction, fibre extraction and wavelength calibration. To remove
the flat-field response from the science data, each extracted spectrum
was divided through by the normalised flat-field spectrum (taken as
part of the daytime calibrations) from the same fibre. Further
manipulation of the data was undertaken using the STARLINK
package DIPSO.
The Galactic targets are relatively bright and sky subtraction was not
a critical issue, however such methods were fully explored for these
data as it will be more important for the fields observed in the Magellanic
Clouds. A master sky spectrum was created from merging the available
sky spectra, rejecting those contaminated with a significantly larger
than average Balmer line emission component or those contaminated by
arc lines (see Sect. 3.1). In a cluster such as
NGC 6611 there is considerable variation in the nebular emission
across very small scales (let alone across the FLAMES field-of-view)
and the subtraction of such features will not be 100
accurate; an
inherent drawback of multi-fibre spectrscopy. The logic behind
rejecting those sky spectra with significantly large nebular emission
is driven by our desire to minimise over-subtraction in the line
cores, although this proved unavoidable in some of the NGC 6611
targets. The extracted flat-field fibres were used to determine
relative throughputs (compared to the median number of counts) for
each fibre; the master sky spectrum was then appropriately scaled
prior to subtraction from the science spectrum.
The exact resolving power of the grating varies with the different central wavelength settings. The mean FWHM of the arc-lines in the central science fibre at each setting (averaged over the three observed fields) is summarized in Table 4. The typical signal-to-noise of the data ranges from 100-150 per resolution element.
Table 4:
Summary of the mean FWHM of the arc lines and effective
resolving power, R, at each Giraffe central wavelength setting, ![]()
.
For classification purposes the spectra for each region were rectified automatically (using pre-defined continuum regions) and then merged (using DIPSO routines developed by I. D. Howarth for echelle/fibre reductions) to give a continuous spectrum of 3850-4755 Å, with additional coverage of 6380-6620 Å.
For programmes that demand high-accuracy radial velocity measurements
there are five Medusa fibres assigned for simultaneous arc calibrations. These are distributed evenly across the CCD and (in the absence of information to the contrary at the time) were left enabled for the Galactic observations. Inspection of these arc spectra showed
no significant wavelength shifts (for a given central wavelength) in
comparison to the daytime calibrations.
However, these checks did
highlight the overspill of the simultaneous arc spectra into those from adjacent
fibres on the CCD. Several of our science spectra contain arc-line contamination
from this overspill, although the contamination is generally weak and
in the majority of cases limited to the red (HR14,
6515) wavelength
setting.
During the course of this survey, the Giraffe Base-Line Data Reduction
Software (girBLDRS) underwent significant development at Observatoire
de Gen
ve (see Blecha et al. 2003, for full details); this
pipeline has been used for reduction of the Magellanic Cloud
observations (Evans et al. in preparation).
Prior to reduction of the Magellanic Cloud data, we re-reduced the observations in NGC 3293 using girBLDRS (v.1.10) for comparison with the IRAF reductions. In brief, we used the calibration data observed at the time of our science frames to update those released with girBLDRS, and then used the "extract'' pipe. The pipeline default is an optimal extraction of each fibre, however this option relies on correction of the localisation of the fibres on the CCD. The master localisation is obtained using the flat-field frames, but the correction employs the simultaneous arc spectra. In anticipation of the fact that the simultaneous arcs were disabled for the majority of our FLAMES observations in the Magellanic Clouds (because of the overspill discussed in Sect. 3.1), the simpler summed extraction method was used for our comparisons. After extraction with girBLDRS, a combined sky spectrum was subtracted from each target using DIPSO, employing the same methods as those used in the IRAF reductions. These re-reduced data were rectified using the same scripts as before and the final spectra are in excellent agreement.
We have conventional longslit spectroscopy of the majority of our
O-type targets in NGC 6611 (see Table 11). These were
observed as part of separate programmes at the 4.2-m William Hershel
Telescope (WHT) on 2003 June 15-17 and October 16-17, using the
Intermediate-dispersion Spectroscopic and Imaging System (ISIS) with
1200B and 1200R gratings. Observations were made at three central
wavelength settings in both arms of ISIS, giving continuous spectral
coverage of 3800-5100 Å, with additional coverage of
6200-6800 Å, at
7000. The basic properties of the ISIS observations are summarized in Table 5, in which they are compared with those from FLAMES-Giraffe and FEROS.
Table 5: Overview of the wavelength coverage and resolution obtained from the different instrumentation employed in the current study.
In general we find excellent agreement between the FLAMES-Giraffe and
ISIS spectra (aside from detection of binarity, see Sect. 6). Indeed, minor differences in the final spectra arise in the far-blue region, not because of problems with the FLAMES spectra but from residual features from the ISIS dichroic used to
simultaneously observe in both blue and red wavelength regions. Of
particular note is that the equivalent widths of important diagnostic
lines (such as He I
4471, He II
4686 and H
)
agree,
adding confidence to the multi-fibre FLAMES data.
In Fig. 6 we compare the H
Balmer line profiles
for the six stars in NGC 6611 observed with both FLAMES-Giraffe and
WHT-ISIS; the FLAMES spectra have been degraded to the effective
resolution of the WHT-ISIS data to ensure a more meaningful
comparison. As one would expect from multi-fibre observations, the
nebular subtraction is not perfect in the FLAMES spectra. However, the
nebular component is easily resolved and the wings of the profiles are
in good agreement with the longslit spectroscopy, and will still offer
strong constraints on the physical parameters. The differences in
Fig. 6 for 6611-011 and 6611-014 are accounted for by
binarity, further discussed in Sect. 6.
![]() |
Figure 6:
Comparison of sky-subtracted FLAMES-Giraffe (solid line)
and WHT-ISIS (dot-dashed line) H |
| Open with DEXTER | |
Table 6: Total numbers of stars in each cluster by spectral type (incorporating both the FLAMES and FEROS observations).
The FLAMES-Giraffe spectra of our Galactic targets were classified in the MK system (Morgan et al. 1943) by visual inspection. The Giraffe spectra have a much greater resolution than those traditionally used in spectral classification, so the data were degraded to an effective resolution of 1.5 Å and then classified according to published morphological precepts. Work is underway by other groups (e.g. Bagnulo et al. 2003) to compile high-resolution observations of spectral standards and our survey will significantly add to such efforts. Furthermore, some of our stars have robust spectral types in the literature and so also serve as internal standards for the current dataset. The classification criteria applied to each spectrum are summarized below, points of interest and spectral peculiarities are also discussed. Following initial classification by CJE, the data were reinspected to ensure internal consistency and, if necessary, minor revisions of spectral type were made. A subset of these data were independently classified by DJL to provide an external check on the spectral types; agreement to within one spectral subtype was found in all cases. The distribution of observed stars by spectral type is summarised in Table 6.
The O-type spectra were classified using the digital atlas of
Walborn & Fitzpatrick (1990) as the principal reference; the primary classification
criterion is the ratio of He I to He II lines. A total of 13 O-type
stars were observed, all of which are in NGC 6611. The standard
suffixes are used to indicate the intensity of N III emission and the
behaviour of He II
4686 (see Walborn & Fitzpatrick for further
details). FLAMES-Giraffe spectra of three O-type stars in NGC 6611
are shown in Fig. 7; note that 6611-080 is one of the
faintest targets in the Galactic sample and yet the data are of
sufficient quality to permit detailed analysis in a future study.
![]() |
Figure 7:
Illustrative FLAMES O-type spectra in NGC 6611. For display purposes the
spectra have been smoothed and rebinned to a resolution of 1.5 Å FWHM and are offset by 0.75 continuum units. A section of the longslit ISIS spectrum of 6611-014 is also shown, illustrating its binarity (see Sect. 6). The sky spectrum has been
subtracted from the blue-region data but, to highlight the nebular
contribution, it has not been removed in the red. The spectral lines
identified in 6611-080 are, from left to right by species,
H I |
| Open with DEXTER | |
The observed B-type spectra were classified from comparisons with the standards in the Walborn & Fitzpatrick (1990) atlas, which extends as late as B3 for dwarfs. These standards were complemented by unpublished observations of B-type dwarfs using the Intermediate Dispersion Spectrograph (IDS), at the 2.5-m Isaac Newton Telescope (INT) in 1999 December, and with FEROS in 2004 July. These additional standards are summarized in Table 7.
Some of the Giraffe spectra are assigned the non-unique spectral type of B6-7 (see Tables 9 and 10). Examples of B6 and B7-type supergiants were presented by Lennon et al. (1992), with subtle variations seen in the metallic spectra as one progresses from B5 through to B8. The FLAMES spectra however are largely of less-luminous objects and such subtleties are not seen because of the intrinsically weaker metallic lines in non-supergiants. Additionally, rotational velocities are typically larger in non-supergiant B-type stars, which further confuses the issue. As a result, here we classify those spectra that are clearly between the observed B5 and B8 standards as "B6-7''.
Table 7: Summary of B-type classification standards in addition to those from Walborn & Fitzpatrick (1990).
The width of the lines in the B-type spectra offer
evidence of a wide range of rotational velocities, some are
particularly sharp-lined, whilst others are significantly broadened. For now
we limit ourselves to simply adding the conventional "n'' suffix to
those B-type spectra which display the most significant broadening, employed
here when the He I
4471 and Mg II
4481 lines are strongly blended.
A thorough analysis of the distribution of rotational velocities will be
presented elsewhere (Lee et al. in preparation).
Luminosity types were assigned to the B-type spectra using
measurements of the H
equivalent widths and the calibrations
of Balona & Crampton (1974). In some cases the measured widths are mid-way
between the published values for class III and V objects, resulting in
an adopted class of "III-V''.
Balona & Crampton (1974) largely confirmed the earlier results of Petrie (1965) for B-type spectra. Their (photographic) work remains the most useful and reliable source of such calibrations, although some less wide-ranging studies (primarly concerned with absolute magnitudes of supergiants, of use in consideration of extragalactic distances) have made use of more modern detectors (e.g. Hill et al. 1986). We envisage that with completion of our survey of Galactic and Magellanic Cloud targets we will be in a position to revisit these relations with an unprecedented (and fully homogenous) dataset. Indeed, two of the clusters observed by Balona & Crampton (1974) were NGC 3293 and 4755; pending a new study we do not explore these issues further in the current work.
Digital observations of A-type spectra were presented by Evans & Howarth (2003).
Their principal classification criterion uses the intensity of the
Ca K line compared to the H
/Ca H blend. Although their
work was principally motivated by studies of A-type supergiants, the
adopted scheme was found to be largely insensitive to luminosity
effects (e.g. their Fig. 4). However, reliance on the calcium line
ratio is not without problems (see discussion by Evans & Howarth),
and with the high quality of the FLAMES data we are also able to give
consideration to the metallic features in the spectra, which increase
as one progresses to later-types. Thus, the observed A-type stars
were classified using the criteria of Evans & Howarth (2003), with secondary
consideration to the intensity of the metal-lines. Luminosity types
were allocated using the H
equivalent-width criteria from
Evans et al. (2004b)
As a result of the significant and variable reddening in the direction
of NGC 6611, our target selection methods yielded numerous (22 of 73)
targets with spectral types of F0 or later. Without prior
spectroscopy to make informed decisions in target selection, this is
not a surprising result. It is not our intention to analyse these
stars in the future, but they are included here for completeness. The
classification bins applied to these later-type spectra are largely
those described by Evans & Howarth (2003), with additional reference to the
useful on-line digital spectral atlas of
R. O. Gray
, for those in the region G5-K2.
Useful as a diagnostic of stellar winds and Be-type stars, when of
particular note the morphology of the H
profiles of our targets is
summarized in the final column of Tables 9-11. The majority of the sample display H
in absorption, with many of these also showing superimposed (narrow) nebular
emission. Such stars would be best described as "abs+em'' but are
not explicitly noted in the tables - exceptions to this are
3293-001 and 4755-003 in which the emission is broader and more likely
stellar in origin. The notation employed to describe the other
H
profiles is: "broad em'' for broad emission; "twin'' for
twin-peaked emission (usually synonomous with disk-like features);
"abs+twin'' for those spectra in which twin-peaked emission is seen
superimposed on an absorption profile; and "P Cyg'' for stars showing a P Cygni emission profile.
A small number of emission line objects were observed in our Galactic
fields, the majority of which can be categorized as normal Be-type
spectra. Three spectra (namely 3293-022, 3293-045, and 4755-014) have
strong twin-peaked H
emission profiles; five (3293-027,
3293-040, 4755-038, 4755-050 and 4755-057) have twin-peaked emission
visible in the wings of the H
lines, superimposed on
absorption profiles; and four (3293-011, 4755-018, 6611-010 and
6611-028) display broad, single-peaked emission features. These have
all been classified as Be or Ae stars. Broad H
emission is
also seen in 6611-022, which is discussed further in Sect. 4.6.
Further caution regarding the problems of nebular subtraction is
prompted by the significantly large number of Be/Ae objects in
NGC 6611 reported by Hillenbrand et al. (1993). From observations with a
slitless, grism spectrograph, Herbig & Dahm (2001) attribute most of the
emission features as due to nebular contamination,
rather than from intrinsic stellar phenomena. For instance, the
H
emission in four targets reported as Be-type by
Hillenbrand et al. is best described here as "absorption with narrow
(presumably nebular) emission'' (6611-020, 064, 066 and 068). Each of
these is reported by Herbig & Dahm as having H
absorption
profiles (albeit their data are at relatively low resolution); these
results reinforce our cautious approach to employing the Be-notation
A number of peculiarities were noted in the process of classification, which are now discussed in turn:
![]() |
Figure 8:
FLAMES-Giraffe spectra of particular interest. The spectral lines identified in 6611-006 are, from left to right by species, He II |
| Open with DEXTER | |
Strong interstellar Ca II H and K lines at
3933 and
3968 are observed towards all three clusters, although blended
with stellar features for the cooler targets. For the two older clusters
(NGC 3293 and NGC 4755) the absorption is not resolved into separate
components and probably arises mainly in the warm diffuse interstellar
medium; no other interstellar lines are apparent although
it should be noted that the spectral region containing the Na I
doublet was not observed.
The interstellar spectra for the NGC 6611 targets are far richer, likely reflecting the substantial circumstellar and interstellar material associated with this young cluster. For example, Thé et al. (1990), Hillenbrand et al. (1993) and Yadav & Sagar (2001) have all discussed both the variable extinction across the cluster and the non-standard reddening law. Furthermore, Welsh (1984) found evidence for high velocity interstellar components in IUE spectra. He interpreted these as arising from a shock-front (that is driven by stellar winds) interacting with the ambient neutral interstellar gas. In the NGC 6611 FLAMES-Giraffe spectra, the Ca II lines show significant spatial variations - with a single strong component (which could be unresolved multiple components), complemented by up to four additional high velocity components. The strong component probably again reflects the warm diffuse interstellar medium, whilst the other features are likely to be spatially associated with the cluster.
Many sightlines show molecular absorption lines (with significant
variations in the line strengths) from CH (
4300), CH+ (
4232), and CN (
3875), implying the presence of a cold molecular cloud. In some targets the Ca I
4227 line is also
observed, further evidence of a dense, optically-thick cloud.
Finally, He I absorption at
3889 (from the 2s
3p
transition) is observed towards some objects, but with large
spatial variations. This was first observed by Wilson (1937) and arises in an ionized plasma in which the helium ions are recombining to a metastable triplet state (e.g. Oudmaijer et al. 1997).
The source of the ionization in this instance is presumably the strong X-ray
continuum emitted by the O-type stars in the cluster.
In general there is good agreement between our classifications and those extant in the literature, with many of the fainter stars now receiving more refined classifications, as one would expect from the new high-quality data. In Table 12 we give previous spectral types for our target stars. Published classifications are included from: Bidelman (1954); Bosch et al. (1999); Conti & Leep (1974); Feast (1958, 1963); Hernandez (1960); Hillenbrand et al. (1993); Hiltner et al. (1969); Hiltner & Morgan (1969); Hoffleit (1956); Morgan et al. (1955); Schild (1970); Turner et al. (1980); Walker (1961); Walborn (1972, 1973, 1976, 1982). We omit the spectral types in NGC 6611 from de Winter et al. (1997) as they are largely based on photometric methods, and are thus "low weight'' classifications.
Many of the targets in NGC 3293 (92 of 126 stars) and NGC 4755 (86 of 108 stars) do not have previous classifications from the sources compiled in Table 8. Even in NGC 6611, in which the spectral content has been well explored in the past (most significantly by Hillenbrand et al. 1993), we provide classifications for 31 stars without previous types in the literature. Our main science drivers do not require spectroscopic completeness but, to be thorough, we note that due to observational constraints (primarily the fact that only one field configuration was observed) a small number of B-type stars with published spectral types from Hillenbrand et al. (1993) are not included in the current survey.
Radial velocities (
)
of our targets were measured from
the Doppler shift of the H
,
H
,
H
,
and H8
(H
) Balmer lines. This relatively simple approach (compared to
using multiple metallic lines) was adopted to ensure internal
consistency over the wide range of observed spectral types. The line
centres were measured by Gaussian fits to the inner region of the
profile, or by manual estimates of the position of minimum intensity.
Only the inner part of the Balmer lines (
3 Å of the centre)
were used for the Gaussian fits because of the effect of linear Stark
broadening on the wings; the whole profile is poorly fit by a
Gaussian function but, in general, the inner part is successfully matched.
This approach worked well for symmetric profiles and was also
particularly effective when weak nebular emission features were
present in the line centres, because the fits are relatively
unaffected. In cases where the lines were very sharp or asymmetric
this approach was not successful and the line centres were determined
from manual measurements. In the double-lined binary 4755-024, two Gaussian profiles were fitted, with the quoted value being that of the primary.
Heliocentric radial velocities are included in Tables 9-11. The values are plotted for the stars in each cluster in the left-column of Fig. 9;
the upper panels show the mean velocity for each star, with the
standard deviation (
)
of the individual measurements
shown in the lower panels. The radial velocity distribution for each
cluster is shown in the right-column plots of Fig. 9.
There are a number of specific comments regarding the radial velocity
measurements. In the O- and B-type spectra the interstellar Ca H
3967 line is easily resolved and was ignored when fitting the H
line. In contrast, in the A-type spectra the stellar
component of the Ca H line becomes increasingly strong, and the H
line was excluded from calculation of the mean velocity. Similarly, as discussed in Sect. 4.7, some stars in
NGC 6611 display sharp and strong He I absorption close to the
centre of the H8 line; this line was not included in calculation of the
mean velocities for any of the NGC 6611 stars. For the G- and K-type
spectra in NGC 6611 the quoted velocities are the means of manual
measurements of several Fe I lines (![]()
4071.7, 4528.6, 4602.9)
and the Cr I
4254.3 line.
The mean velocity and its associated dispersion for each cluster is given in Table 8. In calculation of the mean velocities, their dispersion around the mean, and the mean standard deviation of the measurements within each stellar spectrum, it is appropriate to exclude cluster non-members. In NGC 3293 and NGC 4755 we consider stars with spectral types of B8 or later as probable non-members. Given the imposed faint magnitude limit applied in target selection this is reasonable - indeed any B8 or later-type stars observed are often classified as giants and are therefore likely to be non-members of these relatively young clusters (or at least non-coeval with the observed upper main-seqeunce). Following similar arguments in NGC 6611, we do not include stars of type B5 or later.
Table 8:
The mean radial velocity of the clusters
,
with their dispersions
,
and the mean error in the individual measurements
.
The number of stars, n, used for each cluster (after removal of outliers) is also given.
![]() |
Figure 9:
Distribution of stellar radial velocities. The left-column plots show the
radial velocity of each star, with the 1-sigma scatter shown in the
lower panel for each cluster. The right-column plots are the binned
histograms for each cluster, with a 10
|
| Open with DEXTER | |
The mean of the standard deviations of the individual radial
velocities was calculated to give an estimate of the expected
uncertainty (referred to hereafter as
), and this can
be directly compared to the velocity dispersion measurements in the
clusters to determine if we are measuring statistically significant
velocity dispersions. These values (i.e.
), the mean
radial velocity of each cluster (
)
and the dispersion (
)
are also affected by the stars with variable radial velocities, some of which are the outliers in the plots in Fig. 9. To exclude these we employed a consistent method. The
mean (
)
and standard deviation (
)
of the velocity
distributions were calculated, and any star lying more than 2
from
was rejected and the process repeated. This iteration converged quickly and resulted in clear outliers being removed. The same method was used to determine the
intrinsic mean error in our analysis
in each
cluster. The values are given in Table 8, together with
the final number of stars used, n. For each of the clusters the
measured velocity dispersion is less than 2
,
i.e., we
do not see any strong evidence for detecting real velocity dispersions. The measurement errors could possibly be reduced using more sophisticated velocity measurement techniques, such as
cross-correlation with model spectra, however the rationale for that is beyond the scope of this paper. Note that the mean standard deviation of the individual measurements
ranges from 3-5
,
approximately the one-pixel resolution of the spectra.
Lastly, significant nebular emission or broad, asymmetric profiles made precise location of the line centre difficult in a number of stars; these are flagged with the usual ":'' identifier in the tables to indicate a greater uncertainty (
20
)
than for the other values.
Figure 9 shows that there is a tendency for the
bright cluster members in NGC 4755 to have a smaller approach velocity
than the fainter ones. This was suggested by Feast (1963), although
he saw this effect when he combined data from three clusters. We see
no evidence in NGC 3293 and NGC 6611 for the effect, however in
NGC 4755 the
values for the bright sample (star numbers
20, or V< 10) is indeed significantly different from the faint sample (star numbers > 20, or V> 10), see Table 8. Whether this is physically meaningful remains unclear.
The repeated observations of the NGC 6611 field offer an opportunity
to check for single-lined spectroscopic binaries. We find evidence of
two stars with companions from this method, 6611-030 and 6611-068.
Similarly, the July 18 spectrum of 6611-007 appears to show a
relatively sharp-lined early B-type star, however, the July 15
observation reveals it as a double-lined spectroscopic binary. A B0.5 V spectral type is assigned to one component because of weak He II
4686 absorption; the other component is also early B-type, approximately B1.
The two FLAMES observations of 6611-014 display some small profile
asymmetries but otherwise appear largely unremarkable. However,
inspection of the ISIS spectrum of this star reveals it to be a
double-lined binary. From the FLAMES spectra we had initially
classified this object as O6 V((f(+))), the suffix indicating the
N III emission and evidence of Si IV emission at
4116. From
the ISIS spectrum the He I
4471 line of the hotter component
appears weaker than for an O6 and we assign a spectral type of O5 V((f(+))). The cooler component appears to be that of a late O-type with He I
4471 stronger than the components of the He II lines at
4200 and 4541. This confirms the suggestion of Bosch et al. (1999) that 6611-014
(W175) comprises two O-type stars (note that a third, relatively
low-mass component was reported by Duchêne et al. 2001).
Two further stars in NGC 6611 were identified as binaries by
Bosch et al. (1999), namely 6611-001 (W412) and 6611-003 (W197).
They found that the former is a short period (
4 days)
single-lined binary, but the FEROS and ISIS spectra are in
good agreement and thus offer no additional insight. In the case of 6611-003 (initially classified here as O7) Bosch et al. noted the presence of weak Si III absorption that is inconsistent with a mid O-type classification, suggesting a cooler companion. Indeed, shifts
of up to 100
are seen in the Si III lines between the ISIS and FEROS spectra, similarly with the C III blend at
4650. The intensity of the
4200 and
4542 He II lines is invariant, but the
4686 line is weaker in the ISIS spectrum, reinforcing the suggestion by Bosch et al. of an approximately B0-type secondary - such a star would contribute to
4686, but have less of an impact on the other He II lines. The presence of such a secondary
could justify revision of the classification of the primary to an earlier type of O6.5, or even O6, depending on the contribution to the He I spectrum by the secondary; therefore we adopt the type employed by Bosch et al. of O6-7 V((f)).
For the remaining 7 stars studied by Bosch et al., our radial
velocities are in excellent agreement. However, the comparison of the
H
Balmer line profiles (see Fig. 6) highlights
the disagreement between the FLAMES and ISIS spectra of 6611-011 (W314); no differences were revealed between the two FLAMES observations. Closer inspection of the blue-region FLAMES data
reveals slightly peculiar profiles for a number of the He I lines,
namely ![]()
4026, 4121, 4388 and 4471, with possibly weak line-doubling seen in the latter two. Both Trumpler (from Walker 1961) and Bosch et al. also suggested possible secondary components. The exact nature of this system is still not entirely clear but the
new data here, in particular the H
profiles in Fig. 6, provide further evidence of a possible companion. That the lines are seen to vary slightly between different epochs may
also account for the slightly earlier classification assigned here to
the FLAMES-Giraffe spectrum of O9 V (cf. B0 V from Bosch et al. 1999; Hillenbrand et al. 1993).
We have only one observational epoch for NGC 3293 and NGC 4755,
but inspection of the FLAMES spectra revealed a double-lined
binary, namely 4755-024. The principal component appears as a
relatively broad-lined B2.5 dwarf, whilst the relative ratio of the
He I
4471 line to the Mg II
4481 in the secondary spectrum
suggests a spectral type of approximately B5 for the companion.
In comparison to published values, the radial velocities in Sect. 5 also provide an opportunity to detect potential binaries in NGC 3293 and NGC 4755. Feast (1958) published velocities for a number of stars in NGC 3293, 20 of which were observed with FEROS or
FLAMES. In general there is reasonable agreement between the two sets of results, but a few stars have significantly discrepant (
20
)
velocities. In NGC 3293 these stars are 3293-008 (F007), 3293-013 (F005), and 3293-014 (F019), all of which were identified by Feast as having variable velocities. Furthermore, Feast suggested that 3293-029 (F028) displayed evidence of line doubling and that it was probably a binary. The FLAMES-Giraffe
spectrum of this star displays broad diffuse lines, with some very mild asymmetry in the weaker He I lines such as
4144 and
4388. Its radial velocity (8
)
is different to Feast's
singular measurement (-13
), but his value was flagged as uncertain and the true nature of this star is still not clear.
In NGC 4755 our results are in good agreement with those of
Hernandez (1960). Comparison with the results from Feast (1963) again
highlights a number of potential binaries, namely, 4755-005 (II-23),
4755-007 (III-05), 4755-012 (III-07), 4755-013 (IV-18) and 4755-021 (I), where the aliases are again taken from Arp & van Sant (1958). Interestingly Feast presented many measurements for 4755-005, commenting that it is a double-lined binary; our derived velocity of -55
is different to those tabulated by Feast, and the spectrum does not appear to contain a secondary component.
We have introduced our FLAMES survey of early-type stars in Galactic and Magellanic Cloud open clusters, in particular explaining the philosophy employed in target selection in the Galactic fields. For one of the clusters, NGC 3293, we have reduced the FLAMES data using the pipeline reduction software (girBLDRS) and also using more manual methods in IRAF, finding excellent agreement in the final rectified spectra. Similarly, longslit ISIS spectra of stars in NGC 6611 are also in excellent agreement with the FLAMES data, providing further confidence in the multi-fibre observations.
The spectra were classified, with discussion of specific peculiarities, potential Be-type stars and evidence of binarity. In particular, with the benefit of the high-resolution FLAMES spectra we have detected interstellar He I toward some of the NGC 6611 stars, indicative of some highly ionized plasma in the region of the cluster. Stellar radial velocities are given for the majority of our targets, which help to highlight cluster non-members and further potential binaries.
The distributions of spectral types for the three Galactic clusters are summarized in Table 6; in total we have observed 319 Galactic stars. These observations represent a significant increase in the known spectral content of NGC 3293 and NGC 4755, and will serve as comparison stars for our forthcoming data in the Magellanic Clouds.
Acknowledgements
We are extremely grateful to Francesca Primas and the staff at both Paranal and La Silla for their invaluable assistance with the observational programme, and to Andre Blecha and Gilles Simond at the Observatoire de Genve for their helpful responses to our girBLDRS queries. We also thank Nolan Walborn for numerous insightful discussions, Damian Christian and Peter van Hoof for their input on the identification of the non-stellar helium absorption, the referee for their constructive comments, and Cameron Reed for his input on the correct identification of CPD-59
4553. CJE acknowledges financial support from the UK Particle Physics and Astronomy Research Council (PPARC) under grant PPA/G/S/2001/00131. AH, SSD, and MRV thank the Spanish Ministerio de Educacion y Ciencia under project PNAYA 2001-0436. FN acknowledges support through grants AYA2003-02785E and ESP2002-01627. This paper is based in part on observations made with the Isaac Newton and William Herschel telescopes, 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.
To further complicate matters, when the coordinates from the CD catalogue (epoch 1875, available from the Centre de Donn
es astronomiques de Strasbourg, CDS) for
CD-57
3346 are precessed to epoch 1900, they do not tally with those published by Morgan et al. (1955) for MCW 1181. A similar check for CD-57
3340 (MCW 1179) finds agreement in the positions, suggesting possible typographical or numerical errors in the case of MCW 1181.
The spectra of MCW 1179 (CD-57
3340) and MCW 1182 (CD-57
3348) were both classified by Morgan et al. (1955) as B0.5 V. These stars were matched by Feast (1958) as CPD-57
3507 (his #14) and as CPD-57
3500 (#16), which are 3293-019 and 3293-010
respectively in the current study - although as a result of the identical classifications it is not clear which CPD star refers to which MCW target. Again use of the SIMBAD database highlights
apparent mismatches. The SIMBAD information for MCW 1179 refers to Feast #10, and for MCW 1182 to Feast #22 (the very star to which Feast appears to have attributed the spectral type for MCW 1181). As for MCW 1181, we were again frustrated in our search for the explicit
source for the cross-references employed by SIMBAD.
Star I-05 from Arp & van Sant (1958) was identified as CPD-59
4553
by Hernandez (1960). However, SIMBAD matches CPD-59
4553 as
star R from Feast (1963)
and precession of the coordinates in the CPD catalogue leads
to similar conclusions. Perry et al. (1976) reported V = 8.86 for
CPD-59
4553; a likely explanation is that stars
R and S were spatially unresolved in their study, i.e. the composite
object was the original entry in the CPD catalogue and that the cross-identification
of Hernández was incorrect.
Table 9:
NGC 3293: observational Parameters of Target Stars.
The identifications in column two are from: F, Feast (1958);
T, Turner et al. (1980); HM, Herbst & Miller (1982). WFI photometry is quoted for V > 10.75
,
brighter than this there were saturation problems and values are from Turner et al.
and Feinstein & Marraco (1980). Radial velocities (
)
are given in
.
Instrument codes refer to FEROS (F) and Giraffe (G).
Table 10:
NGC 4755: observational Parameters of Target Stars. The identifications in column two are those of Arp & van Sant (1958); uncertain matches from the finding charts are denoted by "
''. WFI photometry is quoted for V > 11.5
,
brighter than this there were saturation problems and values are from Sanner et al. (2001). Radial velocities (
)
are given in
.
Instrument codes refer to FEROS (F) and Giraffe (G).
Table 11:
NGC 6611: observational Parameters of Target Stars. Cross-identifications in column two are from W, Walker (1961), K, Kamp (1974), and T, Tucholke et al. (1986). WFI photometry is quoted for V > 11.47
,
brighter than this there were saturation problems and values are from Hillenbrand et al. (1993). Radial velocities (
)
are given in
.
Note that stars 6611-001, 6611-005 and 6611-045 are actually outside of the FLAMES field-of-view. Photometry for 6611-040 and 6611-045 is also from Hillenbrand et al. (1993). Instrument codes refer to FEROS (F), Giraffe (G), and ISIS (I).
Table 12: Comparison of current classifications with published spectral types. Sources of spectral types are: B54 (Bidelman 1954); BMN99 (Bosch et al. 1999); CL74 (Conti & Leep 1974); F58 (Feast 1958); F63 (Feast 1963); H60 (Bidelman, from Hernandez 1960); HMSM (Hillenbrand et al. 1993); HGS69 (Hiltner et al. 1969); HM69 (Hiltner & Morgan 1969); H56 (Hoffleit 1956); MCW55 (Morgan et al. 1955); S70 (Schild 1970); T80 (Turner et al. 1980); TrW61 (Trumpler, from Walker 1961); W61 (Walker 1961); W72 (Walborn 1972); W73 (Walborn 1973); W76 (Walborn 1976); and W82 (Walborn 1982).