A&A 488, 151-159 (2008)
DOI: 10.1051/0004-6361:20079255
J. Borissova1 - V. D. Ivanov2 - M. M. Hanson3 - L. Georgiev4 - D. Minniti5 - R. Kurtev1 - D. Geisler6
1 - Departamento de Física y Astronomía, Facultad
de Ciencias, Universidad de Valparaíso,
Ave. Gran Bretaña 1111, Playa Ancha, Casilla 5030,
Valparaíso, Chile
2 -
European Southern Observatory, Av. Alonso de Córdoba 3107, Casilla 19,
Santiago 19001, Chile
3 -
University of Cincinnati, Cincinnati, Ohio 45221-0011, USA
4 -
Instituto de Astronomía, Universidad Nacional Autónoma de México, Apartado Postal 70-254, CD
Universitaria, CP 04510 Mexico DF, Mexico
5 -
Departmento de Astronomía y Astrofísica, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 782-0436 Macul, Santiago, Chile and Vatican Observatory, Specola Vaticana, V00120 Vatican City State, Italy
6 - Departmento de Astronomía,
Universidad de Concepcion, Casilla 160-C, Concepcion, Chile
Received 16 December 2007 / Accepted 16 April 2008
Abstract
Context. Recent near- and mid-infrared surveys have brought evidence that the Milky Way continues to form massive clusters.
Aims. We carry out a program to determine the basic physical properties of the new massive cluster candidate [DBS2003] 179.
Methods. Medium-resolution K-band spectra and deep near-infrared images of [DBS2003] 179 were used to derive the spectral types of eight member stars, and to estimate the distance and reddening to the cluster.
Results. Seven of ten stars with spectra show emission lines. Comparison with template spectra indicated that they are early O-type stars. The mean radial velocity of the cluster is
km s-1. Knowing the spectral types of the members and the color excesses, we determined extinction
and distance modulus (
mag (
Kpc). The presence of early O-stars and a lack of red supergiants suggests a cluster age of 2-5 Myr. The total cluster mass is approximated to
and it is not yet dynamically relaxed.
Conclusions. The candidate [DBS2003] 179 further increases the family of the massive young clusters in the Galaxy, although it appears less massive than the prototypical starburst clusters.
Key words: Galaxy: open clusters and associations: general - Galaxy: open clusters and associations: individual [DBS2003] 179 - stars: early-type - stars: emission-line, Be
More than 2000 infrared ``hidden'' star cluster candidates have been discovered in the last 5-6 years from all-sky IR surveys such as the 2 MASS (Skrutskie et al. 1997) and the Spitzer Space Telescope Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE, Benjamin et al. 2003) and can be used to test the recent star formation in the Galaxy.
This work is a part of a larger program aimed at characterizing the hidden current star population in the Galaxy (Ivanov et al. 2002, 2005; Borissova et al. 2003, 2005, 2006; Kurtev et al. 2007). Our main goal is to determine if the Milky Way still forms fairly massive star clusters at a rate similar to the rate seen in other normal spiral galaxies (Larsen 2002).
The open cluster candidate [DBS2003] 179 was discovered by Dutra et al. (2003) via visual inspection of 2 MASS images during their survey of embedded clusters in the areas of known radio and optical nebulae. The cluster is located in the direction toward the H II ionized region G347.6+0.2 (Fig. 1). Caswell & Haynes (1987) measured peak brightest temperature of 5 K at 5 Ghz continuum and angular size of 9 arcmin for G347.6+0.2.
![]() |
Figure 1:
Pseudo-true color rendering of the region around [DBS2003] 179 outlined with a yellow circle. Blue, green, and red represent the Spitzer [3.6], [4.5], and [5.8] bands. The field of
view is
![]() |
Open with DEXTER |
Borissova et al. (2005, hearafter Paper I) obtained high quality
infrared (IR) photometry and confirmed that [DBS2003] 179 is a member
of the large family of ``newly'' discovered heavily obscured clusters.
Based on the photometric analysis and on the so called ``10th brightest star method'' (Dutra et al. 2003a), they determined some of the cluster parameters: AV=19, (
m - M)0=13.5 mag, age of 7 Myr, and total mass of
.
Thus, the cluster seems to be young and moderately massive.
However, the photometry alone is not sufficient to derive an accurate
distance to the star cluster and the uncertainty of the 10th brightest star method can easily be as large as (
mag or a factor of
10 in
distance (see the discussion in Paper I). To obtain more accurate
physical parameters for this cluster we embarked on a project
to determine the spectral types of some of the cluster members using
K-band spectroscopic observations.
The next section describes the data. In the third section, we present the spectral classification of the observed stars and their radial velocities. In Sect. 4 the cluster parameters and the total mass of the cluster are determined. The last section is a summary of the results.
We carried out near-IR broad
and narrow band imaging on Aug 2, 2003 with the near-IR imager Persson's Auxiliary Nasmyth Infrared Camera (PANIC) at the 6.5 m Baade telescope at the Las
Campanas Observatory. The instrument uses a
HgCdTe
Hawaii detector array. The scale is 0.125 arcsec
,
giving a total field of view of
arcmin. The Br
filter is centered on 2.165
m and has a width of 0.03
m.
The images are reduced following the ``standard'' procedures for the
IR, the photometry is obtained using the DAPHOT package within
Image Reduction and Analysis Facility (IRAF)
and the typical photometric uncertainties are 0.05-0.06 mag. We refer
the reader to Paper I for more details. Since the goal of the
narrow-band observations was to select the sources with IR excess, no flux
standards were taken and no attempt was made to transform the
Br
instrumental magnitudes to any standard system.
Based on their position in the statistically-decontaminated near-IR
color-magnitude diagram given in Paper I, four probable members of the
[DBS2003] 179 were selected for follow-up spectroscopy. Two of
them were among the brightest objects and the rest are chosen to be the
8th and 10th brightest cluster members. Medium-resolution (R=9000) K-band
spectra were obtained with Infrared Spectrometer And Array Camera (ISAAC) (Moorwood et al. 1998) on the ESO VLT at Paranal, Chile in service mode, in May 2007. We used the
SWS1-SW mode, 0.3 arcsec wide slit, and the grating was centered at
2.134 m. Two radial velocity standards HD 191639 (B1V,
km s-1) and HD 167785(B2V,
km s-1) were also observed with an identical
spectroscopic set up. We observed with two slit positions. The exposure
times were 20 min, the seeing during the observations was less than
1 arcsec. However, some additional stars fell into the long slit
(120 arcsec length). We present spectra of ten stars in total for
the two slit positions (Fig. 2). Their coordinates and
magnitudes are given in Table 1.
The near-IR magnitudes are from Paper I. The Br
instrumental
magnitudes are from this work. The mid-IR magnitudes are from the
Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE;
Benjamin et al. 2003) archive. The ``sat.'' marks a saturated star.
Table 1: Spectroscopically observed stars in [DBS2003] 179.
![]() |
Figure 2:
All stars observed spectroscopically with ISAAC at the
ESO VLT are marked. The field of view is
![]() |
Open with DEXTER |
The spectra were dark-subtracted and flat-fielded using the ESO ISAAC
pipeline. To correct for telluric emission the object was ``nodded''
between two positions on the slit (A and B) so that the background of
position A is recorded while the source is at position B and is
subtracted from the observation when the source is at position A,
and vice-versa. The offset between the A and B positions was chosen
to avoid source overlap. The spectra were extracted using the IRAF
task apall. For each position
along the slit a night sky spectrum was extracted with the same
extraction parameters. The wavelength calibration of the spectra are
performed with the imprinted sky emission lines. The accuracy of
the wavelength calibration is 5-7 km s-1. Telluric
absorption lines were removed using the telluric star Hip 087164
of spectral type B2, observed under identical sky conditions and
using the method described in Hanson et al. (1996, 2005) and Bik et al. (2005, 2006). First, the photospheric lines (Br
and HeI) in the
spectrum of the telluric B standard was divided out.
We did this with the template spectrum of the
B2 star HD 36166 from the spectral atlas of Hanson et al. (2005).
The error on the resulting Br
equivalent width
(EW) of our target star is about 10%.
After the removal of the stellar lines from
the telluric standard star we remove the telluric lines from
the program objects by taking the ratio of the target
spectrum with this telluric spectrum. We did the telluric-line
correction using the IRAF task telluric, which allows
for a shift in wavelength and a scaling in line strength yielding a
more accurate fit. The task uses a cross-correlation procedure to
determine the optimal shift in wavelength and the scaling factor in
line strength. The shifts are usually less than one pixel (1 pixel
corresponds to 17 km s-1).
![]() |
Figure 3:
The J-H versus
![]() |
Open with DEXTER |
Typically, the location of the stars in the color-color diagrams such as
J-H vs.
allow us to separate between cluster members
and field red giants and to estimate their individual extinction
values. The two color diagram of all stars in our field of view is
shown in Fig. 3. The most probable statistically-selected
cluster members taken from Paper 1 are plotted with large
filled circles. They are determined using a comparison field,
with an area equal to that of the cluster.
The solid lines represent the unreddened sequence
of stars of luminosity classes I and V, taken from Koornneef (1983)
and (Schmidt-Kaler 1982). The majority of the cluster members
occupy a well-defined locus in this diagram at
-0.8 and
mag.
The intrinsic
and J-H color of OB stars vary
between -0.05 and -0.2 (Schmidt-Kaler 1982; Koornneef 1983).
Furthermore, the calibrations of Martins & Plez (2006) give the same colors
for the O3I-O9I; the O3III-O9III and the O3V-O9V type stars:
(
J-K)0=-0.21, (
J-H)0=-0.11, (
H-K)0=-0.1. We use these to
calculate the individual color excesses of the spectroscopically-classified stars:
), E(J-H) and
)
are given in Table 2. We calculated the E(B-V) using the color excess ratios of E(J-K),
), and
E(J-H) from Bessell et al. (1998) and Mathis (1990).
The errors of E(B-V) reddening
determination given in Col. 6 are calculated as a standard
deviation of the mean value and quadratically added errors from
photometry and the variance due to different intrinsic estimations.
In the last column of Table 2 the extinction
corrected
magnitude are given. We adopted
E(B-V)=5.2 for the mean reddening of the cluster and thus
,
using the standard reddening low.
Table 2: Individual reddening of the spectroscopically observed stars in [DBS2003] 179.
We performed a preliminary spectral classification of the program stars, using direct comparison of spectra between a star of known spectral and luminosity class and a star which is unknown. The appropriate template spectra of early stars are taken from spectral atlas of Hanson et al. (2005). Some of the spectra from this atlas are taken with the Subaru telescope and are of higher spectral resolution (R=12 000) and are rebinned to R=9000. In general, our K-band spectra can be divided into two different classes: seven objects are emission line stars and in three of them we can identify only absorption lines.
Our spectral setting, centered at 2.134 m, contains the
following lines for the hot stars: Brackett
(4-7)
2.1661
m (Br
); He II
2.188 (7-10);
He I
2.1127 (
S, triplet), and
2.1138 (
P
-
S, singlet). Among the
hottest stars, the C IV triplet (
-
)
at 2.069,
2.078 and 2.083 is also seen. Hanson et al. (1996) identified the
broad emission feature found at 2.1155
m as N III (7-8).
Recently, Geballe et al. (2006) showed that O III 2.116
m
(8-7) transitions feature is an important contributor to the 2.11
m
emission complex.
The spectra of the stars showing emission lines
are shown in Figs. 4 and 5. The Equivalent Widths (EWs) of their
emission lines are listed in Table 3. They are measured on the continuum
normalized spectra, by the IRAF task splot, using the deblending function. The errors of the EW are estimated taking into account the signal-to-noise ratio (see Col. 2 of Table 3), the peak over continuum ratio of the line (see Bik et al. 2005) and the error from the telluric star subtraction and are estimated to be 10-15% in the worst cases. The EW of the emission lines are negative. The symbol ``Em'' given in Col. 1 of the Table 3 indicates that the line is in emission, ``NP'' means that the line is not present in the spectrum.
![]() |
Figure 4:
Left panel: the emission line stars of [DBS2003] 179. The different spectra have been arbitrarily shifted along the flux
![]() |
Open with DEXTER |
The C IV line is at the limit of our setup, where the errors
are large and it is not given strong consideration. The N III/O III
emission complex appears in all these stars, except for Obj6, where it is
not present. The Br
line appears weakly in absorption in Obj1
(EW=0.8 Å), moderately weak in absorption in Obj2, Obj3, and Obj5 (EW=3.4-3.7 Å)
and weakly in emission in Obj8 (EW=-1.4 Å). The Br
line appears strongly in
emission in Obj4 and Obj6, (EW=-37.2 Å and EW=-12.1 Å, respectively).
The He II
2.188 line is seen in absorption in Obj1, 2, 3, 5, and 8.
We derive the following spectral classification of our objects:
Object 1 and Object 2 are quite similar, given the limitation
of the ``edge'' effects, where both the C IV 2.078 m and
He II 2.1885
m lines are compromised. The EW of the
N III/O III line of Obj1 is -2.8 Å which, according to
Hanson et al. (2005) defines the star as later than O3-O4, but not
later than O6, when He I would be much stronger and
C IV line will not be present. The Br
and He II
2.188
m lines are almost equally strong. The extinction corrected
magnitudes of the stars are 9.9 and 10.04 mag. Comparison with
the smoothed to R=9000 Hanson et al. (2005) templates shows, that Obj1 and Obj2 are
similar to CygOB2-8C-O5If and HD46150-O5Vf. Unfortunately, the S/N is not very
high and for the most part in many spectral ranges, dwarfs, giants, and many
supergiants look nearly the same. Our classification for Obj1 and Obj2 is
early O5-O6 spectral type, but it is hard to distinguish between the
luminosity class I and V.
Object 3 and Object 5 show He II 2.188 m stronger than
than in Obj1 and Obj2, and therefore, higher temperature. The stars are
similar to HD66811-O4If and HD46223-O4Vf. Obj3 is slightly
brighter than Obj1 and Obj2, while Obj5 is fainter. Our spectral
classification is O4-O5, luminosity class I or V.
Object 4 shows Br
and He I lines in strong emission, the He II line is at the very end of our setup and there is no C IV emission. Based on this spectrum, it is extremely hard to tell its spectral type. The star is saturated in our H and
images and have no 2 MASS measurements because of crowding. We can estimate approximately its
magnitude, using narrow-band Br
and Spitzer 3.6 mag values (see Table 1). The mean
from these estimations
is 9.2 mag, which corrected for the mean extinction derived from Table 2 gives
KS0=7.4. That is a magnitude and half brighter than Obj8, which is classified as ``evolved'' O4-O5 If supergiant (see below). So, Obj4 is even more evolved from the main sequence and thus probably more massive (assuming co-eval star formation) than Obj8. If what we see on the red part of the spectrum is really HeII 2.189, then the star is similar to WN7 objects (Fig. 4 of Figer et al. 1997, WR131). If, on the contrary, the strong line at 2.189 is artificial, then the star could be a Ofpe/WR9 like HDE 269582 from the Morris et al. (1996) atlas and to IRS34W (Martins et al. 2007, Fig. 2) or Fig. 3 of Crowther & Bohannan (1997). In addition, the wind velocity is low (
km s-1), as suggested by the HeI 2.112
m line, which also supports the Ofpe classification. Even though more observations are necessary to clarify the nature of this interesting object, which is the brightest in our sample. Our tentative classification for Obj4 is Ofpe/WN9.
Object 6 is probably a massive Young Stellar Object (YSO). The Br
is in
strong emission, Mg II at 2.138 and 2.144
m can
be identified. This is the faintest object in the sample and,
according to the color-color diagram, shows IR excess.
Object 8 shows Br
in emission and also strong
He II 2.188
m in absorption. The star is spectroscopically similar to HD 15570, classified as O4If Hanson et al. (2005). It is significantly brighter than Obj1, 2, 3, and 5 (the extinction corrected
mag). For these reasons, we believe the star may be a supergiant.
The spectra of the absorption line stars are shown in Fig. 6.
Object 10 shows He I 2.1112 m, He II 2.188
m
and Br
,
N III/O III might exist weakly and
the He II line is relatively strong. The star
is similar to HD13268-ON8V type. The extinction corrected
mag is 10.79 mag. This allows us to put the spectral type of
Obj10 as O8 V.
It is impossible to classify the stars Obj7 and Obj9. They have
noisy spectra (S/N=25) and the only visible feature is Brin absorption. Note that from the radial velocity analysis (see below)
these stars are unlikely to be cluster members.
![]() |
Figure 5: Left panel: the emission line stars of [DBS2003] 179. Right panel: the template spectra taken from Hanson et al. (1996). The symbols are the same as in Fig. 4. |
Open with DEXTER |
We derived the radial velocities of all ten stars with the
IRAF task fxcor, which uses a cross-correlation Fourier
method. Since the stars show very different characteristics,
several radial velocity standards are used: for the stars Obj1,
Obj2, Obj3 and Obj5 we used HD 15629, an OV5 spectral type star
with
km s-1; Object 8 was
cross-correlated with HD 15570, an O4If spectral type star with
km s-1; Obj6 and Obj 4 are measured
by comparison with DM+493718, an Be spectral type star; all absorption
line stars are compared with HD 191639 (B1V,
km s-1)
and HD 167785(B2V,
km s-1). The
spectra are taken from Hanson et al (2005) and Hanson et al. (1996)
spectral atlases. The HD 15629 have been shifted to zero velocity,
as measured in air using the sophisticated atmospheric modeling
(Repolust et al. 2005). In this study, the high resolution spectra of
Hanson et al. (2005)
were modeled to derive directly absolute values of luminosity and
temperature. The models (synthesized spectra) are made at zero
velocity and thus the modelers shift the observed spectrum to align
them with the models. It is accurate to
10 km s-1.
The radial velocities and their estimated errors are given in Cols. 8 and 9 of Table 3.
The membership of the stars is determined on the basis of the radial velocity
histogram (see Fig. 7). Most of the stars are concentrated
within the radial velocity interval (-100, -50) km s-1. The stars Obj7 and Obj9 have very different radial velocities and
most probably are back- or foreground stars. The photometry also confirm the non-membership
of these stars: the extinction corrected
magnitudes and
colors are (11.05; -0.3)
and (7.06, 1.29) mag (see Tables 1 and 2)
and they are situated far from the determined cluster boundaries.
The Object 4 also has a relatively high radial
velocity; the reason could be the lack of a suitable template
for comparison or the peculiar nature of the star.
The mean radial
velocity value is then calculated from seven
cluster members (Obj4 is not taken in consideration) and is
km s-1.
Table 3: Spectral classification of the observed stars in [DBS2003] 179.
The cluster seems to be compact and dense. To determine its boundaries we performed direct star counting, assuming spherical symmetry. The projected star number density as a function of radius is shown in Fig. 8. The cluster boundary was determined as the radius at which the density profile exceeds twice the standard deviation of the surface density in the surrounding field. This yields a radius of 32 arcsec.
With the newly-adopted cluster radius, we statistically decontaminated the ``cluster''+``field'' color-magnitude diagram by removing as many stars as are present on the ``field'' color-magnitude diagram. The details of the procedure are described in Borissova et al. (2005). We have selected 180 potential cluster members.
To determine the distance to the cluster, we used the K-band spectral
classification of our objects and the intrinsic color and absolute magnitudes
of the main sequence and supergiants stars
given in the latest calibrations of Martins & Plez (2006).
We feel most confident about the spectral type and luminosity class of
Obj10, as an O8V. Likewise, we feel moderately confident about the classification
of the O supergiant, Obj8, showing weak emission in Br.
However, we do not share the same confidence in knowing the luminosity
class for Obj1, 2, 3, and 5. In Table 4 we have calculated
the distance moduli for each object assuming two luminosity classes for
Obj1, 2, 3, and 5 - as supergiants (I) and dwarfs (V).
As can be seen from Table 4
if we assume the luminosity class I for these four objects, the distance
modulus is more than one magnitude larger than this value calculated
from Obj8 and Obj10. All four of these stars have K0 magnitudes
of about 10 and the fact that these stars show both Br
and HeII in absorption, it seems likely they are main sequence stars.
Finally, the star showing the most extreme spectrum, Obj4, also has the largest KS0 magnitude, estimated to be around 8. This object is
likely an extreme supergiant (Of+).
These luminosity designations of dwarf, supergiant, and Of+ are consistent
with their spectra and the individual distance moduli they predict.
Using them, we adopted the average (
mag or
7.9 Kpc. The error is calculated as a standard
deviation of the mean value and quadratically added errors from
photometry and the variance due to unsure spectral classification.
This is a larger
distance than the value of 5 Kpc estimated in
Paper I. The reason for the difference is in the a priori
assumption used in the 10th star method for the spectral type
of the 10th brightest cluster member.
The
versus
(
color-magnitude diagram of the cluster is
plotted in Fig. 9. We adopted
E(B-V)=5.2 and
(
mag.
The unreddened main sequence (Schmidt-Kaler 1982) is shown with the
solid line as well as the isochrones of pre-main sequence stars
(PMS) of 1 and 4 Myr (Siess et al. 2000). The stars with known
spectral type are plotted with large circles and labeled. In
general, as stated above, their absolute
magnitudes are
consistent with luminosity class V (Obj1, Obj2, Obj3, Obj5, and Obj10)
and I (Obj8).
Our photometry is not deep enough to properly sample the pre-main-sequence stars in [DBS2003] 179 and derive an age based on isochrones shown in Fig. 9. However, the presence of O type supergiants and the absence (so far discovered) of WC stars and/or red supergiants implies an upper limit to the age of 4-5 Myr (Maeder & Meynet 2003). As a lower age limit we adopt 2 Myr taking into account the presence of at least one Ofpe/WN9 type star (Obj4) (Crowther & Dessart 1998; Figer et al. 1999, 2002; Blum et al. 2001; Martins et al. 2007). For comparison we also plotted (Fig. 9) the pre-main sequence isochrones from Siess et al. (2000). Thus, the present stellar population of [DBS2003] 179 appears consistent with an age of order 2-5 Myr.
We determine the total stellar mass of the cluster in the following
way. First, assuming the Salpeter initial mass function, masses of O supergiants 50
(Herrero et al. 2002; Hillier et al. 2003) and masses of
O-type dwarfs 30
(Bouret et al. 2003), we calculated a total mass of 250
for the spectroscopically-confirmed cluster members.
Second, using AV=16.6, (
m-M)0=14.5 mag, and integrating
along the 5 Myr Padova isochrone (Girardi et al. 2002), we
obtain a total present-day mass of the observed main sequence cluster members
down to 12
3100
.
Integrating over 2.5 Myr isochrone gives for the same mass interval
3500
.
In summary, we will adopt the mass of all
observed stars of
3550
.
Since no completeness correction has
been applied this value should be considered as lower mass limit of the cluster.
Our photometry is not deep enough and a lot of cluster members are probably not
yet seen. Nevertheless, some estimation of the expected total cluster mass can be
calculated adopting the Salpeter law, normalizing it to the upper three mass bins
(encompassing stars with masses about 10 solar masses), and integrating down to 0.8 solar masses. This yields of
.
The total number of stars determined with this procedure is 1240.
![]() |
Figure 6: Left panel: the absorption line stars of [DBS2003] 179. Right panel: the template spectra taken from Hanson et al. (2005). The symbols are the same as in Fig. 4. |
Open with DEXTER |
![]() |
Figure 7: The radial velocity histogram. |
Open with DEXTER |
![]() |
Figure 8:
Radial surface density (in number of stars per square arcsec) profile of [DBS2003] 179. The bars show the 3![]() |
Open with DEXTER |
To place [DBS2003] 179 in the context of the currently known
Milky Way supermassive/starburst clusters, we compared it with
the Westerlund 1/2, the Orion Nuclear Cluster (ONC), NGC 3603, Arches, and
R136 clusters. The comparison is given in Table 5,
reproduced from Table 9 in Brandner
et al. (2008). The table rows present: age, relaxation time
,
half-mass radius
,
total number of
cluster members
,
and total cluster mass
.
To calculate the relaxation time of the cluster, we use the method described in Sect. 5.2 of Brandner et al. (2008). First, we calculated the velocity dispersion
as:
Table 4: Distance to the cluster.
![]() |
Figure 9:
The
![]() ![]() |
Open with DEXTER |
Table 5: Comparison of [DBS2003] 179 with other young, massive clusters. The compilation data are from Brandner et al. (2008, Table 9). The footnotes list the sources for the data on each individual cluster.
![]() |
Figure 10: Location of the [DBS2003] 179 and the other massive Milky Way clusters. |
Open with DEXTER |
In summary, [DBS2003] 179 is a young, not yet dynamically relaxed, dense and relatively massive cluster. It has at least six O-type stars with masses in excess of 30
and our rough estimate based on the most massive stars and assuming a Salpeter IMF suggests a cluster mass well in excess of 103, approaching 10
.
The location in the Milky Way of the [DBS2003] 179 and the other clusters from Table 5 is shown in Fig. 10.
Already, a number of optical and infrared clusters have been identified with characteristics similar to [DBS2003] 179, a few of which are listed in Table 5. Presumably, all the optically observable clusters of this mass and age (any age) have been identified. But how many more massive young clusters do we expect current near-infrared searches to find? If our Milky Way galaxy is like every other normal spiral galaxy, its global star formation rate would suggest nearly 100 young clusters similar to [DBS2003] 179 should currently exist within it (Weidner et al. 2005). Just how could so many massive young clusters exist without our seeing them? The 2 MASS survey, mostly because of the extinction limitations at the J- and to a lesser extent the H-bands, is not seeing very deeply into our Galaxy. Typical infrared clusters discovered with 2 MASS have distances of between 1 and 2 kpc (Carpenter 2004). Newer, deeper and longer-wavelength surveys look more promising and suggest that they are sampling to several kpc from the Sun (Mercer et al. 2007). However, it is frequently overlooked that most of the star formation in our galaxy occurs deep within the molecular ring at a galactocentric distance of 4-5 kpc. Except at the closest tangent point to this structure, this constitutes a solar distance of more than a few kpc. The deeply buried, very young massive cluster system of W49 may, in fact, exist within this ring. However, no other known massive infrared (least optical) cluster is positively identified with this region of the Galaxy. Surely, we are missing a significant number of massive clusters associated with this galactic structure.
Infrared searches must continue, but care must be given not to interpret the initial results too quickly. Enormous biases in our current search methods are, in fact, leading to a highly skewed view of the massive cluster population presently in the Milky Way. Until these biases are addressed, it will be difficult to fully appreciate the massive cluster population in our galaxy and relate it properly to the same phenomenon being observed in extragalactic super star cluster work.
Acknowledgements
J.B. acknowledges support from Proyecto Fondecyt Regular # 1080086. D.M., J.B. and D.G. are supported by Centro de Astrofísica FONDAP No. 15010003 and the Chilean Centro de Excelencia en Astrofísica y Tecnologías Afines (CATA). Support for RK is provided by Proyecto Fondecyt Regular # 1080154 and DIPUV grant No. 36/2006, Universidad de Valparaiso, Chile. JB acknowledges support from DIPUV grant No. 6/2007, Universidad de Valparaiso, Chile. MMH acknowledges support by the National Science Foundation under grant 0607497 to the University of Cincinnati. We are grateful for helpful comments from the referee.