New massive members of Cygnus OB2
Key Words.:massive stars – OB-type – Cygnus OB2 – new members
Context: The Cygnus complex is one of the most powerful star forming regions at a close distance from the Sun (1.4 kpc). Its richest OB association Cygnus OB2 is known to harbor many tens of O-type stars and hundreds of B-type stars, providing a large homogeneous population of OB stars that can be analyzed. Many studies of its massive population have been developed in the last decades, although the total number of OB stars is still incomplete.
Aims:Our aim is to increase the sample of O and B members of Cygnus OB2 and its surroundings by spectroscopically classifying 61 candidates as possible OB-type members of Cygnus OB2, using new intermediate resolution spectroscopy.
Methods:We have obtained intermediate resolution (R5000) spectra for all of the OB-type candidates between 2013 and 2017. We thus performed a spectral classification of the sample using HeI-II and metal lines rates, as well as the Marxist Ghost Buster (MGB) software for O-type stars and the IACOB standards catalog for B-type stars.
Results:From the whole sample of 61 candidates, we have classified 42 stars as new massive OB-type stars, earlier than B3, in Cygnus OB2 and surroundings, including 11 O-type stars. The other candidates are discarded as they display later spectral types inconsistent with membership in the association. We have also obtained visual extinctions for all the new confirmed massive OB members, placing them in a Hertzsprung-Russell Diagram using calibrations for and luminosity. Finally, we have studied the age and extinction distribution of our sample within the region.
Conclusions: We have obtained new blue intermediate-resolution spectra suitable for spectral classification of 61 OB candidates in Cygnus OB2 and surroundings. The confirmation of 42 new OB massive stars (earlier than B3) in the region allows us to increase the young massive population known in the field. We have also confirmed the correlation between age and Galactic longitude previously found in the region. We conclude that many O and early B stars at mag are still undiscovered in Cygnus.
The Cygnus region is the most powerful nearby stellar complex, conspicuous at all wavelengths and very young, with several rich OB associations, numerous young open clusters and tens of compact H II and star formation regions in the field. Hosting the largest number of nearby massive stars and an intense star forming activity (Reipurth & Schneider, 2008), it provides an updated view of the high-mass stellar population in one of the largest groups of young stars in our Galaxy. It is an ideal place to study the process of massive star formation and evolution, individually and in stellar groups, and their interaction with the surroundings.
Its association Cygnus OB2 (d1.4 kpc, Rygl et al., 2012) has received a lot of attention and has been studied at all wavelengths with different spatial coverage since it hosts a high number of early spectral type stars (Walborn et al., 2002). First studies were carried out by Morgan et al. (1954), Schulte (1956, 1958) and Reddish (1968), but it was Massey & Thompson (1991) who developed an extensive survey of the massive population in the association, identifying 120 possible massive star members, 70 of which were classified as OB stars (42 O-type stars). Knödlseder (2000) proposed that this number should be much larger, around 100 O-type stars. In the last few decades, many other studies were carried out in the region updating continually this number (Comerón et al., 2002; Hanson, 2003; Kiminki et al., 2007; Comerón et al., 2008; Negueruela et al., 2008; Comerón & Pasquali, 2012). These surveys have allowed the global study of the massive population in the region, using photometry to place the stars in a Hertzprung-Russell diagram from which the star formation history and mass function of the association were assessed (Wright et al., 2015). In spite of the many photometric and spectroscopic surveys carried out in the region, only a small homogeneous group of early type stars have been spectroscopically analyzed (Herrero et al., 1999, 2002; Negueruela et al., 2008), and few stars have been observed in the UV range (Herrero et al., 2001). There is still a large number of stars that should be explored. The optical extinction of the region is high ( mag, Wright et al., 2015), but not so much as to prevent obtaining spectra of its most massive stars for a rough spectral classification. Therefore, new spectroscopy to search for previously undiscovered massive stars is mandatory to complete the last census of massive O and B-type stars in the association.
One of the most complete spectroscopic surveys in the Cygnus region was developed by Comerón & Pasquali (2012). They performed spectral classification of a magnitude-limited sample ( 16 mag and mag) selected with a homogeneous photometric criterion over a large area that includes Cygnus OB2 and its surroundings, providing a large sample of known and new OB stars, as well as a list of 61 OB candidates for which no spectral data is available and that are pending spectroscopic confirmation.
The main goal of this work is to complete the spectral classification of this latter sample, aiming a later determination of the stellar parameters. Thus, we have obtained intermediate-resolution spectra of all the list candidates, in order to confirm or reject them as true massive OB-type stars. In Sect. 2 we present the observations. In Sect. 3 we describe the spectral classification criteria used in this work, and in Sect. 4 visual extinctions and stellar parameters derived for the new OB-type members. The results are discussed in Sect. 5 where we show the Hertzsprung-Russell diagram (HRD) of the region and the age distribution found across Galactic longitude. Finally, we summarize our conclusions in Sect. 6.
2 Observations and data reduction
The study developed by Comerón & Pasquali (2012) produced a sample of O and early B stars (in a 6 x 4 region centered on Galactic coordinates 79.8 and +0.8 of Cygnus OB2) which were identified using two homogeneous reddening-free criteria
Q_BJK = 0.196 (B-J)-0.981 (J-K)-0.098 ¿ 0
Q_JHK = 0.447 (J-H)-0.894 (H-K)-0.089 ¡ 0 which allowed them to classify 60 new OB stars and produce a list of 61 candidates pending spectroscopic data. They used photometry tabulated in the USNO-B (Zacharias et al., 2010, 2004) and 2MASS all-sky (Skrutskie et al., 2006) catalogs setting limiting magnitudes of mag and mag. By combining these two magnitude cuts with the colors of OB stars, Comerón & Pasquali (2012) set up a selection method sensitive to main-sequence stars earlier than B1 and obscured by mag.
We have obtained new spectra for all the proposed OB candidates, whose location is shown in Fig. 1. The sample has been observed in five different runs between 2013 and 2017. For an accurate spectral classification of OB-type stars we need blue spectra (4000-5000 ) where the diagnostic He I-II and metal lines are located.
|ISIS - 600B||ORM-WHT||2500||Oct 2013||1|
|WYFFOS - H2400B||ORM-WHT||5000||Jul 2014||7|
|ISIS - H2400||ORM-WHT||7500||Jul 2015||2|
|OSIRIS - R2500U/V||ORM-GTC||2500||May 2016||2|
|IDS - R1200B||ORM-INT||5000||Jul 2016||47|
|ISIS - R1200B||ORM-WHT||5000||Apr 2017||2|
The bulk of stars were observed in July 2016, obtaining spectra of 47 candidates. We have chosen the R1200B grating of IDS (EEV10) on the Isaac Newton Telescope (INT) in La Palma, which provides a resolution of 5000 at 4500 . The remaining stars were observed in different runs between 2013 and 2017, at the William Herschel Telescope (WHT) using the AF2/WYFFOS and ISIS instruments, and at the Gran Telescopio CANARIAS (GTC) using the OSIRIS instrument. Out of the whole sample, three of the stars (J20301097+4120088, J20315433+4010067, J20345785+4143543) belong to the Galactic O-Star Spectroscopic Survey (GOSSS) catalog (Maíz Apellániz et al., 2016).
The information on the different runs carried out in this work is shown in Table 1. The final spectra were reduced using the IRAF procedures, with standard routines for bias and flat-field subtraction and also for the wavelength calibration.
3 Spectral classification
|Object||RA (hhmmss)||Dec ( )||Region||SpT||Binary star|
|J20423509+4256364||20 42 35.08||+42 56 36.43||c||14.480||8.304||9.211||O6IIIz|
|J20371773+4156316||20 37 17.73||+41 56 31.57||c||15.760||8.041||9.071||O7V|
|J20345785+4143543||20 34 57.84||+41 43 54.25||a||15.430||7.417||8.447||O7:Ib|
|J20293563+4024315||20 29 35.63||+40 24 31.45||b||12.468||8.268||8.831||O8IIIz|
|J20222481+4013426||20 22 24.81||+40 13 42.55||c||13.620||8.410||9.034||O8II|
|J20261976+3951425||20 26 19.75||+39 51 42.46||c||15.600||8.348||9.351||O8.5IV|
|J20275292+4144067||20 27 52.92||+41 44 06.65||b||13.330||7.277||8.144||O9.5II|
|J20262484+4001413||20 26 24.84||+40 01 41.25||c||13.290||8.330||9.021||O9.2III|
|J20291617+4057372||20 29 16.17||+40 57 37.19||b||15.030||7.899||8.855||O9.7III|
|J20382173+4157069||20 38 21.72||+41 57 06.89||c||15.810||7.682||8.760||O9.7II|
|J20181090+4029063||20 18 10.89||+40 29 06.29||c||14.940||8.399||9.343||O9.7Ib||yes*|
|J20273787+4115468||20 27 37.87||+41 15 46.79||b||14.570||8.263||9.146||B0II|
|J20301097+4120088||20 30 10.97||+41 20 08.82||b||15.690||8.882||9.855||B0:II:|
|J20323968+4050418||20 32 39.68||+40 50 41.83||a||14.410||8.913||9.631||B0II|
|J20395358+4222506||20 39 53.58||+42 22 50.62||c||15.890||5.822||7.345||B0I||yes*|
|J20281176+3840227||20 28 11.75||+38 40 22.73||c||11.805||7.944||8.349||B0Ib|
|J20323882+4058469||20 32 38.82||+40 58 46.85||a||15.430||8.821||9.701||B0Ib|
|J20225451+4023314||20 22 54.50||+40 23 31.39||c||13.107||8.601||9.175||B0Iab|
|J20253320+4048444||20 25 33.19||+40 48 44.38||c||13.112||7.648||8.340||B0Iab|
|J20272099+4121262||20 27 20.99||+41 21 26.15||b||13.830||8.730||9.448||B0.5V||yes|
|HDE229258||20 24 25.51||+39 49 28.30||c||10.235||8.689||8.833||B0.7V|
|J20330526+4143367||20 33 05.26||+41 43 36.74||a||13.940||8.634||9.286||B0.5III|
|J20361806+4228483||20 36 18.06||+42 28 48.30||c||15.650||8.855||9.814||B0.7III|
|J20233816+3938118||20 23 38.16||+39 38 11.84||c||11.293||8.731||8.996||B0.7Ib|
|HD228973||20 20 07.35||+41 07 46.72||c||10.34||7.684||7.922||B1V||yes|
|J20201435+4107155||20 20 14.34||+41 07 15.45||c||12.412||8.240||8.627||B1V|
|J20230290+4133466||20 23 02.90||+41 33 46.59||c||14.780||7.809||8.762||B1V|
|BD+404193||20 29 13.55||+40 41 03.38||b||10.412||8.812||8.945||B1V|
|BD+404208||20 30 49.97||+40 44 18.53||b||10.654||8.664||8.869||B1V|
|J20314341+4100021||20 31 43.40||+41 00 02.07||a||15.940||8.957||9.885||B1V||yes*|
|J20315898+4107314||20 31 58.98||+41 07 31.41||a||15.490||8.832||9.773||B1V||yes*|
|J20330453+3822269||20 33 04.53||+38 22 26.91||c||11.383||8.790||9.021||B1V|
|J20230183+4014029||20 23 01.83||+40 14 02.90||c||13.465||8.060||8.579||B1III|
|J20274925+4017004||20 27 49.25||+40 17 00.42||b||13.460||8.104||8.713||B1III|
|J20315433+4010067||20 31 54.33||+40 10 06.71||b||15.990||8.884||9.742||B1III|
|J20382889+4009566||20 38 28.88||+40 09 56.63||c||9.996||8.816||8.896||B1III|
|J20440752+4107342||20 44 07.51||+41 07 34.18||c||15.470||8.703||9.562||B1III|
|LSII+3797||20 33 35.52||+38 01 36.73||c||11.838||5.900||6.691||B1Ia||yes*|
|J20211924+3936230||20 21 19.24||+39 36 22.98||c||12.732||7.886||8.488||B1Ib||yes*|
|BD+394179||20 25 27.28||+40 24 00.15||c||11.066||6.667||7.210||B1Ib|
|J20290247+4231159||20 29 02.46||+42 31 15.91||c||13.330||8.426||9.090||B1Ib||yes*|
|J20381289+4057169||20 38 12.88||+40 57 16.86||c||12.870||8.675||9.183||B2V|
The accuracy of the spectral classification depends on the effects of spectral resolution as well as the signal-to-noise ratio . The main diagnostic method for O-type stars is the comparison of He II 4542/He I 4471 ratio (Sota et al., 2011). These lines are similar for a O7 type star. For later O types the relative strengths of He II 4542/He I 4387 and He II 4200/He I 4144 are normally used, and represent the main criteria for types O8-B0. Spectral types earlier than O8 were classified using the criteria described by Gray & Corbally (2009). The presence of metal lines in different ionization stage, such as Si III or Mg II, indicates early B-type stars. The relative strength of HeII 4471/MgII 4481 is a useful indicator for B1-late B stars. As secondary indicators we used the criteria described by Gray & Corbally (2009), which were also used to classify stars of spectral types A, F and G.
Regarding the luminosity class, the criteria used for early O-type stars were introduced by Walborn (1971, 1973), taking into account the emission effects in the He II 4686 line and N III 4634-4640-4642. For late O-type stars, we have used the criteria described by Sota et al. (2011) and for B-type stars the criteria described by Gray & Corbally (2009) and the Balmer lines width.
Although the classification for O-type stars was based on the described He and metal lines diagnostic criteria, we have also used for O-type stars the Marxist Ghost Buster (MGB) code developed by Maíz Apellániz et al. (2012) in order to obtain a more accurate result. This tool compares the observed spectra with a grid of O standards (in this work the GOSSS library), allowing us to vary spectral type, luminosity class, velocity and resolution until obtaining a best match. Furthermore, for B-type stars we have used a sample of IACOB standards (Simón-Díaz et al., 2015) to improve their classification.
|Object||RA (hhmmss)||Dec ( )|
|late B-type stars|
|CCDMJ20323+4152AB||20 32 20.81||+41 52 00.78|
|BD+423785a||20 34 15.39||+43 09 35.28|
|J20252497+3934030||20 25 24.96||+39 34 03.02|
|J20285874+4013302||20 28 58.74||+40 13 30.22|
|J20315984+4120354||20 31 59.84||+41 20 35.41|
|J20320734+3828586||20 32 07.34||+38 28 58.62|
|BD+413801||20 33 30.39||+42 04 17.35|
|J20364336+3906145||20 36 43.36||+39 06 14.53|
|CCDM J20429+4311AB||20 42 54.24||+43 10 38.71|
|J20193232+4042447||20 19 32.32||+40 42 44.72|
|J20253116+4005508||20 25 31.16||+40 05 50.82|
|J20275204+4131200||20 27 52.03||+41 31 19.98|
|CCDM J20420+4015||20 42 01.18||+40 14 42.70|
|J20500396+4300118||20 50 03.96||+43 00 11.76|
|J20504551+421012.6||20 50 45.51||+42 10 12.64|
|J20300022+4337553||20 30 00.21||+43 37 55.29|
|J20340430+4136507||20 34 04.29||+41 36 50.67|
|J20432737+4308525||20 43 27.37||+43 08 52.47|
|J20472235+4220523||20 47 22.34||+42 20 52.30|
4 New confirmed OB type stars
The observed spectra have high enough and resolution for a spectral classification of the stars. All the candidate spectra are plotted in Appendix B (Fig. 12) where the main diagnostic lines used are also indicated.
Out of the 61 candidates 42 are OB type stars, earlier than B3, including 11 O-type stars. Two more are late B-type stars. The location of these new confirmed OB stars is shown in Fig. 1 while their names, coordinates, magnitudes and the derived spectral classification are listed in Table 2. We have also included the region in which they are located: (a) the Cygnus OB2 area considered by Wright et al. (2015) which is the youngest core of Cyg OB2 at present; (b) the Cygnus OB2 area considered by Comerón & Pasquali (2012) which is the extended Cyg OB2 area containing older stars, on average, and (c) the surrounding area which includes part of the Cygnus OB9 association and field population. The remaining observed stars are late B and foreground A-F-G stars, which are globally listed in Table 3. The selection criteria success rate obtained by Comerón & Pasquali (2012) along with the number of the confirmed OB type obtained in this work (), support the success at identifying reddened massive OB stars with this method.
We have also detected nine possible or confirmed SB2 binaries in the sample of new OB stars (see Table 2). In some cases, indicated with asterisks, we can only suggest possible binary nature mainly due to noisy spectra. This number represents a of our massive OB sample. We assign to these stars the spectral classification of their primary component. An example is shown in Fig. 2. Sana & Evans (2011) found that at least of the O star population in clusters and OB associations is comprized of spectroscopic binaries, which indicates that it is highly likely that more binaries are undetected in our sample.
We used the colors and the recent extinction law derived by Wright et al. (2015) () to obtain individual extinctions for all the new confirmed massive OB type stars (see Table 4). Most of them have visual extinctions in the range mag, which agrees with previous studies in the region (Comerón & Pasquali, 2012; Wright et al., 2015). But, as in Comerón & Pasquali (2012), this is partly consequence of the imposed and magnitude limits.
In Fig. 3 we present the , diagram of the confirmed O-type stars in Cygnus OB2. We observe that some of them are aligned below the O9.5V reddening vector, which could suggest a small shift in the adopted intrinsic color calibration since most of them are already known dwarf late O-type members. In spite of this we can assess the incompleteness of the O population in the region by taking into account the ratios of stars in different spectral (O3-O5, O5-O7 and O7-O9.5) and color bins.
All the stars of our sample with colors mag should be bright enough to have been detected by our selection criteria. Figure 3 shows that the population is expected to be complete for mag. But for mag we can not see late O-type stars because they are fainter than mag, which is the imposed magnitude limit. Assuming that the ratios of spectral types are independent of extinction, we can conclude we are loosing the fainter or more obscured O-type stars.
The extinction distribution for the 42 new confirmed OB stars earlier than B3 in Cygnus OB2 and boundaries is shown in Fig. 4 (top). Again, the spectral classification for those stars classified as possible or confirmed SB2 binaries is taken from the primary component. Although O-type stars are clearly more obscured, a median value of mag was found for the whole sample. However, this value is also affected by completeness. O-type stars are more obscured on average because we can detect them up to higher foreground extinctions thanks to their intrinsic brightness. In the bottom histogram the same sample is differentiated by location: the whole sample that belongs to Cygnus OB2, Cygnus OB9 and boundaries, and as a sub-sample, those stars located within the Cygnus OB2 area (from Comerón & Pasquali, 2012) for which a median value of mag is derived.
We see an abrupt cut for mag. Stars with such high visual extinctions shall be intrinsically extremely bright to be seen. An O9V star with a visual extinction of mag will have an apparent magnitude mag, beyond the selection criteria limits and probably this is the reason why we do not find stars beyond this value. However, for the star J20395358+4222506 classified as B0I, we have obtained a visual extinction of 11.0 mag. It has a magnitude in the B band of 15.89 mag but in the Ks band of only 5.82 mag which indicates a very bright star. These cases show again that the sample does not provide a complete census. The magnitude-limited sample of Cygnus OB2 is made incomplete due to extinction.
In Fig. 5 we present the spatial location of the new classified OB stars, where each star is color-coded according to its derived visual extinction . The extinction distribution varies smoothly across the region, increasing from the south-west (Cygnus OB9) to the northeast (Cygnus X-North) where the most extinguished star of our sample, J20395358+4222506, is located.
4.2 Stellar parameters
For each classified O-type star and luminosity classes V, III and I, we have used the effective temperature () scale developed by Holgado et al. (in prep.) based on a revision of the O-type standards in the IACOB database. We have preferred it over the SpT- calibration suggested by Martins et al. (2005) since the latter gives too low values for late O-type stars (Simón-Díaz et al., 2014).
For early B-type stars and luminosity class V, we decided to use the spectral type compilation of Nieva (2013). An excellent agreement with the calibration of Holgado et al. is found. For late B-type dwarf stars we have used the calibration of Pecaut & Mamajek (2013). Regarding B stars with luminosity class I, we have used the scale of Markova & Puls (2008), and for luminosity class III we have interpolated between classes I and V. In Fig. 6 are represented the different scales adopted in this work, and the derived values for all the new OB type stars are shown in Appendix A (Table 4). We see that all scales for B-types fit smoothly the Holgado et al. scale. Temperatures were also obtained for those stars classified as SB2 binaries by taking into account the primary component.
Comerón & Pasquali (2012) used the ‘observational’ effective temperature based on the versus spectral type calibration of Martins et al. (2005). For B-type stars, they used the spectral type compilation of Tokunaga (2000) but applying a scaling factor to force the agreement with Martins et al. (2005) temperature scale. In order to be consistent with them we have recalculated temperatures for all of their known OB stars sample in the Cygnus Region using our criteria. New values are shown in Appendix A (Table 5).
Bolometric correction and luminosity
As Comerón & Pasquali (2012), we have adopted a distance modulus of mag, and used the intrinsic magnitudes derived by Martins & Plez (2006) for O stars and those compiled by Tokunaga (2000) for B stars. Thus we could derive absolute magnitudes (see Comerón et al., 2008). To derive luminosities we adopted bolometric corrections (BCs) from Lanz & Hubeny (2003, 2007). There is an excellent agreement with Martins et al. (2005) for temperatures higher than 32000 K, and also with the bolometric corrections from Nieva (2013) for cooler stars (see Nieva, 2013). The derived luminosity values for the new and already known OB stars are shown in Table 4 and Table 5 respectively.
Some of the new confirmed OB type stars (see Fig. 1) are located at the boundaries assigned to Cygnus OB2. However, the limits of the association are not strongly defined. The first surveys in the association assumed a smaller area where the most luminous members are located (Münch & Morgan, 1953; Massey & Thompson, 1991). Then, Knödlseder (2000) provided wider limits showing that extinction was a limiting factor in previous studies. This area was thus extended with the identification of new early-type members (Comerón et al., 2002, 2008; Comerón & Pasquali, 2012; Wright et al., 2009, 2015). If we assume the extension of Cygnus OB2 as the area adopted by Comerón & Pasquali (2012) (1 deg. radius centered on Galactic coordinates 79.8 and +0.8), as well as the most recent census in the area developed by Wright et al. (2015), we can now update the census of O and B-type stars in Cygnus OB2 from 204 to 221, and the number of confirmed O-type members increases from 66 to 70 stars.
5.1 Hertzsprung-Russell diagram
Some studies about the age of Cygnus OB2 have been developed in the last two decades (Massey et al., 2001; Hanson, 2003; Drew et al., 2008; Negueruela et al., 2008; Wright et al., 2010; Comerón & Pasquali, 2012; Comerón et al., 2016). However, different evolutionary stellar models were used. We decided to explore four different stellar models (two families with two initial rotational velocities) to construct the Hertzsprung-Russell Diagrams (HRD) in order to assess uncertainties on the age distribution.
We have placed the new confirmed OB stars in the HRD to study the evolutionary status of the association. For this aim we have considered the stars within the area adopted by Comerón & Pasquali (2012). We have also added the already known OB stars compiled in their work to complete the sample. Two of them, classified as B0 stars (J20325964+4115146 and J20333700+4116113) were recently reclassified as O9.7III and O9.5IV by Maíz Apellániz et al. (2016). We thus decided to assume this last classification for both. We have used the Geneva evolutionary tracks and isochrones with and without rotation as calculated by Ekström et al. (2012) (Fig. 7 top), and for comparison, the Bonn evolutionary tracks and isochrones with and without rotation as calculated by Brott et al. (2011) (Fig. 7 bottom).
In the Geneva case, we found a difference up to 2 Myr in the stellar ages for those most massive members ( 20 M) depending on whether we consider rotation or not. Larger stellar lifetimes are derived from rotating stellar models. A similar result was found by Wright et al. (2015) by comparing both non-rotating and rotating Geneva stellar models. They derived an age range of 1 - 7 Myr for the association, but the results from rotating stellar models suggested an age of 4 - 5 Myr, while non-rotating models hinted at a younger age (2 - 3 Myr). On the contrary, the Bonn stellar models included on this work do not exhibit a large difference in the estimated age of the stars when they include or not stellar rotation. This different behavior is most likely to be ascribed to their different treatment of rotation. Thus, the positions of the isochrones in both Bonn HR diagrams are very similar. Moreover, in the Bonn case we see an extended terminal age main sequence (TAMS), and therefore, all the OB stars are located on the main sequence (MS). In spite of this, all models suggest that most of the stars are in the age range of 1 - 6 Myr, indicating a continuous (but not necessarily constant) star formation activity. However, a combination of different scenarios (rotating and non-rotating stellar models) would help to narrow the possible ages. Slow rotators with initial masses of 30 M or less would give similar ages as more massive stars born with larger rotational velocities.
In view of the results, we can conclude that uncertainties due to rotation and adopted models affect mainly the most massive members, suggesting larger ages of about 1 - 2 Myr for them. This makes the age determination very uncertain and strengthens the need for additional diagnostics, in particular individual rotation velocity measurements, to better constrain them.
5.2 Age distribution across Galactic longitude
Comerón & Pasquali (2012) studied the correlation between ages and Galactic longitudes in Cygnus OB2. They suggest that star formation has proceeded from lower to higher Galactic longitudes, finding most of the old stars located at low Galactic longitudes while the youngest ones lie at higher Galactic longitudes. In order to corroborate it, we have obtained ages by comparing isochrones from the HRDs, and divided the sample of OB stars in three age groups: the young group contains stars in the range of 0 - 5 Myr, the intermediate group contains stars in the range of 5 - 10 Myr, and the old one contains stars with ages Myr.
Figures 8 and 9 show the histograms of the relative frequency of our OB type sample located at low (78.5- 79.5), central (79.5- 80.5) and high (80.5- 81.5) Galactic longitudes in each age group and for each model considered (Geneva (Fig. 8) and Bonn (Fig. 9) rotating and non-rotating models). In all cases, most of the stars at high Galactic longitudes belong to the young-intermediate age group, while those stars located at low Galactic longitudes belong to the old age group. We only obtain significant age differences at using the rotating Geneva models since in this case rotation gives older ages for the most massive stars. Thus, we can see a change in the relative frequencies of old and intermediate stars. On the other hand, we do not see these age differences at using rotating Bonn models because of their lower sensitivity to rotation as noted earlier.
This analysis thus supports the correlation between ages and Galactic longitudes in Cygnus OB2 as previously suggested by Comerón & Pasquali (2012). Massive star formation in Cygnus OB2 seems to have proceeded from lower to higher Galactic longitudes, regardless of the details of the models used.
5.3 The whole Cygnus region
In order to obtain a big picture of the age distribution across the Galactic longitude, we have performed the same analysis as in Cygnus OB2 but now in a wider area which includes Cygnus OB2, Cygnus OB9 and boundaries. We have placed the whole sample of OB stars, including late B-type stars, in an HRD (see Fig. 10). Due to the similar results obtained in Cygnus OB2 by using different models, we have now used isochrones and evolutionary stellar tracks for non-rotating models from Ekström et al. (2012). We derived ages by comparison with isochrones and divided again the sample in three age groups: the young group (0 - 5 Myr), the intermediate group (5 - 10 Myr), and the old group ( Myr). In Fig. 11 we show the spatial distribution of the different age groups, where most of the younger stars (left) are concentrated at higher Galactic longitudes, while the older ones (right) are located at central-lower Galactic longitudes. These results suggest that massive star formation has proceeded from lower to higher Galactic longitudes, from Cygnus OB9 to Cygnus OB2, with a strong peak in the northern part of the association.
We have carried out several observing runs between 2013 and 2017 to obtain new blue intermediate-resolution spectra suitable for spectral classification for the magnitude-limited ( mag, mag) candidate list compiled by Comerón & Pasquali (2012). Out of 61 candidates, we confirm 42 new massive OB-type stars, earlier than B3, including 11 new O-type members. A of this sample results on confirmed or possible SB2 binaries. Two other stars are late-B, seven are A-type, six are F-type, and the remaining four are G-type stars. Therefore, the spectral classification of the homogeneously selected sample is now completed. The selection criteria success rate obtained by Comerón & Pasquali (2012) along with this work is , supporting the success at identifying reddened massive OB stars with the reddening-free parameters based on photometry (see Sect. 2). However, the magnitude cutoff and dust extinction introduce an incompleteness. We loose the faintest and more obscured late O-type members in Cygnus OB2, and we are still far to obtain a complete census of the early-type population in the association.
We have also estimated individual visual extinctions () for the new confirmed OB-type stars using the extinction law derived by Wright et al. (2015). We found a median value of 5.5 mag for the stars in the whole Cygnus region, and as we expect, a large median value of 6.5 mag for those stars within the Cygnus OB2 association.
We have placed the new sample of OB stars in the HR diagram using both rotating and non-rotating models calculated by Ekström et al. (2012) and Brott et al. (2011) in order to assess uncertainties. We have also placed the already known OB type stars compiled by Comerón & Pasquali (2012) for completeness. To this aim, we derived effective temperatures and luminosities using spectral type calibrations for the whole sample. Although coming from different sources, our adopted calibrations fit very well with each other. Uncertainties about rotation and adopted models affect only the most massive members, suggesting larger ages for them of 1-2 Myr. Even so, all models support the previous age spread observed in Cygnus OB2 of 1 - 6 Myr.
In order to check the correlation between age and Galactic longitude found by Comerón & Pasquali (2012) in Cygnus OB2, we have divided the sample of OB stars in different location and age groups by comparing with isochrones: the young group (0 - 5 Myr), the intermediate group (5 - 10 Myr), and the old group ( Myr). Despite the differences in physical process treatment (rotation and physics of the stellar interior) done by the Geneva and Bonn groups, the spatial distribution of each age group is globally similar. Therefore, the result obtained by Comerón & Pasquali (2012) can not be an effect of using the Geneva stellar models without rotation when age dating the Cygnus region stars.
Most of the stars at high Galactic longitudes belong to the young-intermediate age group, while stars at low Galactic longitudes belong to the old group. Assuming a wider area which includes Cygnus OB2, part of Cygnus OB9 and boundaries, we obtain similar results suggesting that massive star formation in the region proceeds from Cygnus OB9 to Cygnus OB2, with a strong peak in the northern part of the association. Therefore the correlation between age and Galactic longitude is confirmed, regardless of the details of the models used.
Acknowledgements.We acknowledge financial support from the Spanish Ministry of Economy and Competitiveness (MINECO) under the grants AYA2012-39364-C02-01, AYA 2015-68012-C2-01 and Severo Ochoa SEV-2015-0548. A.S. also acknowledges support from MINECO through grants AYA2013-40 611-P and AYA2016-75 931-C2-2-P. AP and CBM acknowledge support from the Sonderforschungsbereich SFB 881 ”The Milky Way System” (subproject B5) of the German Research Foundation (DFG). We also thank the WHT and its service programme (SW2017a04), J. MaÃz ApellÃ¡niz for his useful MGB software, and N. Wright for helpful discussion and comments in this paper.
Appendix A Table of stellar parameters
Table 4 shows the derived stellar parameters for the new classified OB stars in a large area which includes Cygnus OB2 and its surroundings. We have also included the derived individual visual extinctions. Table 5 shows the updated temperatures and luminosities for the already known OB stars from Comerón & Pasquali (2012) assuming the same spectral calibrations used in this work (see Sect. 4.2).
|Object||RA (hhmmss)||Dec ( )||SpT||(K)||log ()|
|Object||RA (hhmmss)||Dec ( )||SpT||(K)||log ()|
|Object||RA (hhmmss)||Dec ( )||SpT||(K)||log ()|
Appendix B Candidate spectra
- magnitudes from the USNO-B catalog. and magnitudes from 2MASS catalog. Region a indicates the 1 deg. circular area centered on Cyg OB28 trapezium adopted by Wright et al. (2015) for the Cygnus OB2 association. Region b indicates the 1 deg. radius area adopted by Comerón & Pasquali (2012) for the same Cygnus OB2 association. Region c indicates the surrounding area outside the 1 deg. radius adopted by Comerón & Pasquali (2012). For binary stars asterisks indicate possible SB2 stars, whose spectral types are refered to the primary component.
- Brott, I., de Mink, S. E., Cantiello, M., et al. 2011, A&A, 530, A115
- Comerón, F., & Pasquali, A. 2012, A&A, 543, A101
- Comerón, F., Pasquali, A., Rodighiero, G., et al. 2002, A&A, 389, 874
- Comerón, F., Pasquali, A., Figueras, F., & Torra, J. 2008, A&A, 486, 453
- Comerón, F., Djupvik, A. A., Schneider, N., & Pasquali, A. 2016, A&A, 586, A46
- Crowther, P. A., Lennon, D. J., & Walborn, N. R. 2006, A&A, 446, 279
- Drew, J. E., Greimel, R., Irwin, M. J., & Sale, S. E. 2008, MNRAS, 386, 1761
- Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146
- Gray, R. O., & Corbally, C., J. 2009, Stellar Spectral Classification. Princeton University Press, 2009.
- Hanson, M. M. 2003, ApJ, 597, 957
- Herrero, A., Corral, L. J., Villamariz, M. R., & Martín, E. L. 1999, A&A, 348, 542
- Herrero, A., Puls, J., Corral, L. J., Kudritzki, R. P., & Villamariz, M. R. 2001, A&A, 366, 623
- Herrero, A., Puls, J., & Najarro, F. 2002, A&A, 396, 949
- Humphreys, R. M., & McElroy, D. B. 1984, ApJ, 284, 565
- Kiminki, D. C., Kobulnicky, H. A., Kinemuchi, K., et al. 2007, ApJ, 664, 1102
- Knödlseder, J. 2000, A&A, 360, 539
- Lanz, T., & Hubeny, I. 2003, ApJS, 146, 417
- Lanz, T., & Hubeny, I. 2007, ApJS, 169, 83
- Lejeune, T., & Schaerer, D. 2001, A&A, 366, 538
- Maíz Apellániz, J., Pellerin, A., Barbá, R. H., et al. 2012, Proceedings of a Scientific Meeting in Honor of Anthony F. J. Moffat, 465, 484
- Maíz Apellániz, J., Sota, A., Arias, J. I., et al. 2016, ApJS, 224, 4
- Markova, N., & Puls, J. 2008, A&A, 478, 823
- Martins, F., & Plez, B. 2006, A&A, 457, 637
- Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049
- Massey, P., & Thompson, A. B. 1991, AJ, 101, 1408
- Massey, P., Johnson, K. E., & Degioia-Eastwood, K. 1995, ApJ, 454, 151
- Massey, P., DeGioia-Eastwood, K., & Waterhouse, E. 2001, AJ, 121, 1050
- Morgan, W. W., Johnson, H. L., & Roman, N. G. 1954, PASP, 66, 85
- Münch, L., & Morgan, W. W. 1953, ApJ, 118, 161
- Negueruela, I., Marco, A., Herrero, A., & Clark, J. S. 2008, A&A, 487, 575
- Nieva, M.-F. 2013, A&A, 550, A26
- Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9
- Portegies Zwart, S. F., McMillan, S. L. W., & Gieles, M. 2010, ARA&A, 48, 431
- Reddish, V. C. 1968, The Observatory, 88, 139
- Reipurth, B., & Schneider, N. 2008, Handbook of Star Forming Regions, Volume I, 4, 36
- Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618
- Rygl, K. L. J., Brunthaler, A., Sanna, A., et al. 2012, A&A, 539, A79
- Sana, H., & Evans, C. J. 2011, Active OB Stars: Structure, Evolution, Mass Loss, and Critical Limits, 272, 474
- Schneider, N., Bontemps, S., Simon, R., et al. 2006, A&A, 458, 855
- Schulte, D. H. 1956, ApJ, 124, 530
- Schulte, D. H. 1958, ApJ, 128, 41
- Simón-Díaz, S., & Herrero, A. 2007, A&A, 468, 1063
- Simón-Díaz, S., Herrero, A., Sabín-Sanjulián, C., et al. 2014, A&A, 570, L6
- Simón-Díaz, S., Negueruela, I., Maíz Apellániz, J., et al. 2015, Highlights of Spanish Astrophysics VIII, 576
- Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163
- Sota, A., Maíz Apellániz, J., Walborn, N. R., et al. 2011, ApJS, 193, 24
- Tokunaga, A. T. 2000, Allen’s Astrophysical Quantities, 143
- Walborn, N. R. 1971, ApJS, 23, 257
- Walborn, N. R. 1973, AJ, 78, 1067
- Walborn, N. R., Howarth, I. D., Lennon, D. J., et al. 2002, AJ, 123, 2754
- Wright, N. J., Drake, J. J., & Drew, J. E. 2009, American Astronomical Society Meeting Abstracts #213, 213, 605.10
- Wright, N. J., Drake, J. J., Drew, J. E., & Vink, J. S. 2010, ApJ, 713, 871
- Wright, N. J., Drew, J. E., & Mohr-Smith, M. 2015, MNRAS, 449, 741
- Zacharias, N., Urban, S. E., Zacharias, M. I., et al. 2004, AJ, 127, 3043
- Zacharias, N., Finch, C., Girard, T., et al. 2010, AJ, 139, 2184