Issue |
A&A
Volume 557, September 2013
|
|
---|---|---|
Article Number | A51 | |
Number of page(s) | 8 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/201321694 | |
Published online | 27 August 2013 |
CEN 34 – high-mass YSO in M 17 or background post-AGB star?⋆
1 Purple Mountain Observatory, Chinese Academy of Sciences, 2 West Beijing Road, 210008 Nanjing, PR China
e-mail: zwchen@pmo.ac.cn
2 Astronomisches Institut, Ruhr-Universität Bochum, Universitätsstrasse 150, 44801 Bochum, Germany
3 University of Chinese Academy of Sciences, 100039 Beijing, PR China
4 Instituto de Astronomía, Universidad Católica del Norte, Avenida Angamos 0610, Casilla 1280 Antofagasta, Chile
Received:
14
April
2013
Accepted:
10
July
2013
We investigate the proposed high-mass young stellar object (YSO) candidate CEN 34, thought to be associated with the star-forming region M 17. Its optical to near-infrared (550−2500 nm) spectrum reveals several photospheric absorption features, such as Hα, the Ca ii triplet, and the CO bandhead, but lacks emission lines. The spectral features in the range 8375−8770 Å are used to constrain an effective temperature Teff = 5250 ± 250 K (early-/mid-G) and a log g = 2.0 ± 0.3 (supergiant). The spectral energy distribution (SED) displays a faint infrared excess that resembles that of a high-mass YSO or an evolved star of intermediate mass. Moreover, the observed temperature and surface gravity are identical for high-mass YSOs and evolved stars. The radial velocity of CEN 34 relative to the local standard of rest (VLSR) as obtained from various photospheric lines is of the order of −60 km s-1 and thus distinct from the +25 km s-1 found for several OB stars in the cluster and for the associated molecular cloud. The SED modeling yields 10-4 M⊙ of circumstellar material, which contributes only a tiny fraction to the total visual extinction (11 mag). The distance of CEN 34 is between 2.0 kpc and 4.5 kpc. In the case of a YSO, a dynamical ejection process is proposed to explain the VLSR difference between CEN 34 and M 17. Additionally, to match the temperature and luminosity, we speculate that CEN 34 had accumulated the bulk of its mass with an accretion rate >4 × 10-3M⊙/yr over a very short time span (~103 yrs), and it is currently undergoing a phase of gravitational contraction without any further mass gain. However, all the aforementioned characteristics of CEN 34 are compatible with an evolved star of 5−7 M⊙ and an age of 50−100 Myr, so it is most likely a background post-AGB star with a distance between 2.0 kpc and 4.5 kpc. We consider the latter classification as the more likely interpretation. Further discrimination of the two possible scenarios should come from the stricter confinement of CEN 34’s distance.
Key words: stars: late-type / stars: AGB and post-AGB / stars: formation / stars: pre-main sequence
© ESO, 2013
1. Introduction
Before reaching the main sequence, young stars accumulate the bulk of their mass during the so-called pre-main sequence (PMS) evolutionary phase. This phase can last as long as a few Myr for low- and intermediate-mass stars ( < 8 M⊙), while it can be as short as ~104−105 yr for high-mass stars (>8 M⊙; Zinnecker & Yorke 2007). During the majority of their PMS phase, young stars are surrounded and obscured by significant amounts of molecular gas and dust; the younger the forming stars, the more deeply embedded they are in their natal environments, and thus, the more difficult their detection and characterization through observations at optical and near-infrared (NIR) wavelength (e.g., Nielbock et al. 2007). The detection of high-mass young stellar objects (YSOs) is rare, because these objects are deeply embedded, but equally hampering is the extremely short PMS evolutionary phase, which makes it statistically more difficult to find them. Furthermore, the search for and the identification of high-mass YSOs is hampered by confusion with much more evolved objects, showing similar or even identical observational signatures. For instance, if infrared (IR) excess is considered as the sole criterion for source selection, high-mass YSO candidates and post-asymptotic giant branch (post-AGB) stars might be easily mixed up, although they are totally distinct in evolutionary phase. In both cases the observed IR excess originates in absorption of stellar photons by circumstellar dust and subsequent re-emission as thermal radiation at longer wavelengths. Alternatively, both the forming high-mass YSOs and the evolved stars have larger stellar radii and correspondingly cooler effective temperatures than the main sequence stars. This also prevents distinguishing high-mass YSOs from evolved stars solely from a spectroscopic point-of-view.
Fig. 1 The projection of CEN 34 (red cross) on the sky overlaid on the JHKs pseudo-color image of M 17 taken by the InfraRed Survey Facility (IRSF) telescope. The ionizing sources of M 17-the O4-type binary, and the hypercompact H ii region-M 17-UC1 are also marked. |
As one of the brightest H ii regions in the Galaxy, M 17 shows evidence of multiple epochs of massive star formation (Jiang et al. 2002; Povich et al. 2009; Chen et al. 2012). It features two prominent bars, the northern H ii bar (N-bar) and the southwestern H ii bar (S-bar), which are ionized by a couple of OB-type stars located in the center of the H ii region (Chen et al. 2012). A number of intermediate- and high-mass YSOs have been discovered in the S-bar (Chini et al. 2004; Nürnberger et al. 2007; Nielbock et al. 2007; Ochsendorf et al. 2011). In addition, the S-bar harbors several high-mass YSO candidates, which display rising spectral energy distributions (SEDs) towards IR wavelengths but still await confirmation of their preliminary classification (Nielbock et al. 2001). Among these high-mass YSO candidates, CEN 34 (also known as B 358) has its own confusing peculiarities. It is very bright in the near IR (K(2MASS) = 7.88 mag) and rather faint at optical wavelengths (B 22.8 mag, V 18.8 mag), which is indicative of high extinction. Its K-band spectrum shows the CO bandheads longwards of 2.29 μm in absorption (Hanson et al. 1997; Hoffmeister et al. 2006), reminiscent of late-type Class II/III YSOs (Casali & Eiroa 1996; Hoffmeister et al. 2006). However, the same feature is typical of late-type stars, as well as of evolved red stars (e.g., Wallace & Hinkle 1997). Therefore, the true nature of CEN 34 is unclear.
To constrain the effective temperature and luminosity class of CEN 34, which is crucial to better understanding its nature, we performed new spectroscopic observations with high sensitivity and large wavelength coverage (from the optical to the NIR), using the medium-resolution spectrograph X-shooter on ESO’s Very Large Telescope (VLT). In this paper we report on our findings with respect to CEN 34’s spectral classification, its radial velocity in comparison to that of other cluster members, its extinction and its most likely distance.
Fig. 2 X-shooter spectrum of CEN 34, covering the visible and near-IR wavelength range. Prominent photospheric absorption lines are marked. Gray shaded regions mark the major telluric absorption features. |
2. Observations and data reduction
The VLT/X-shooter observations of CEN 34 (, J2000) were performed on August 12, 2009 during the first science verification run of X-shooter, which is the first of the second-generation instruments at the VLT (D’Odorico et al. 2006; Vernet et al. 2011). The average optical seeing and airmass during the observations are and 1.376, respectively. All observations were acquired by nodding the star on the slit using an ABBA sequence, allowing for proper sky subtraction and bad pixel removal. On each nodding pointing, the integration time is 685 s for the UVB arm (300−600 nm), 570 s for the VIS arm (550−1000 nm), and 33 s for the NIR arm (1000−2500 nm). In total, the on-source integration time is 45 min for the UVB arm, 38 min for the VIS arm, and 2.2 min for the NIR arm. The data taken in the UVB arm were not satisfying because the faintness of CEN 34 provided an insufficient signal-to-noise (S/N). For the VIS and the NIR arms, a slit width was used to guarantee R = 8800 (Δλ = 0.09 nm) and R = 5600 (Δλ = 0.3 nm), respectively. The S/N is ~30 in the VIS arm and ~20 in the NIR arm. The standard procedures of data reduction were applied using the X-shooter pipeline version 1.3.7 (Goldoni 2011), including bias-, dark-, and flat-field correction, wavelength calibration, order tracing, and flux calibration. Regarding the wavelength calibration, the transformation from pixel space to (λ, s) space was performed by analyzing a nine-pinhole Th-Ar lamp frame with a list of lines and a table containing the expected positions of the lines in the frame. The accuracy of wavelength calibration using the pipeline is better than ~2 km s-1 over the entire wavelength range of X-shooter. The standard star EG 274 (a white dwarf) was used to calibrate the fluxes for the spectra of the three individual arms. Meanwhile, the telluric absorption correction was conducted for the spectrum of the NIR arm by using the telluric standard stars HD 100858.
3. Results
The relative location of CEN 34 with respect to the H ii region is shown in the finding chart (see Fig. 1). The most massive stars, the O 4-type binary (Hoffmeister et al. 2008), are located on the left edge of the finding chart, meanwhile the giant molecular cloud M 17-SW is located on the right edge.
The flux-calibrated spectrum of CEN 34 from the visual to the NIR is shown in Fig. 2. The optical and NIR parts of the spectrum are quite noisy due to the faintness at optical wavelengths and the short exposure time at NIR wavelengths, respectively. Only absorption lines are detected, such as Hα, Ca ii triplet, and CO bandheads in the K-band. The analysis based on the Ca ii triplet and CO bandheads are presented in the following sections.
3.1. Spectral type and luminosity class
Our spectral classification of CEN 34 relies on the wavelength range 8375–8770 Å, which avoids strong telluric absorption. It also offers a rich combination of lines relevant for astrophysical diagnostics and spectral classification (Gray & Corbally 2009). First, neutral Fe and Ti lines appear from late-F stars onwards, and grow gradually in intensity with decreasing temperature until late-K stars. Second, the Paschen series (e.g., Paschen 14 line 8598 Å, hereafter P14) disappear for main sequence stars later than G0, but they can still be seen for mid-G supergiants due to the luminosity effect (Munari & Tomasella 1999). Third, neutral Ti and Fe lines, the Ca ii triplet (8498 Å, 8542 Å, 8662 Å), as well as the lines of the Paschen series can provide luminosity discrimination by comparing their absolute strengths.
Fig. 3 X-shooter spectrum of CEN 34 (emphasized in red) in the wavelength region 8375–8770 Å. The two features apparently seen in emission in the CEN 34 spectrum at wavelengths ~8590 Å and ~8730 Å are artifacts with unnatural line widths. For comparison, the spectra of several reference stars are displayed; they are adopted from the literature as follows: F8 III (HIP 74975), G1 III (HIP 117503), and G1 II (HIP 47908) from Valdes et al. (2004); G0 Ib (HIP 96481) and G0 III (HIP 5454) from Andrillat et al. (1995); G2 Ib (HIP 109074), G2 III (HIP 86731), G5 III (HIP 100524), and G5 Ib (HIP 107348) from Carquillat et al. (1997). All reference spectra have a resolution around 1 Å, very close to that of our CEN 34 spectrum. All spectra are normalized and are displayed by applying a corresponding shift along the Y-axis for better visual comparison. |
Line parameters.
Figure 3 shows the wavelength range 8375–8770 Å of the X-shooter spectrum of CEN 34 in comparison to the spectra of several reference stars of luminosity classes I (supergiants); II (bright giants); and III (normal giants). Prominent absorption lines – like the Ca ii triplet, neutral metal lines (Fe i, Ti i, and O i), and the P14 line – are marked, and their derived parameters listed in Table 1. The appearance of Fe i, Ti i, and Paschen lines in the spectrum of CEN 34 indicates a rough temperature class of late-F to mid-G. The absolute strengths of these absorption lines are greater than those for giants, suggesting a luminosity class II or I. In addition, Ti i and Paschen lines stronger than in the spectrum of the F8 III reference star can be due to the brighter luminosity class (II or I) and/or later temperature class.
Fig. 4 CEN 34’s X-shooter spectrum (black solid) around the Ca ii 8498 Å, P14, and Ca ii 8662 Å lines. The best fitting grid (red) and other nearby grids are shown as dashed lines. The synthetic spectral library is based on the Kurucz codes, and has an uniform dispersion of 1 Å/pix (for more details see Munari et al. 2005). |
3.2. Effective temperature and surface gravity
Effective temperature and gravity of CEN 34 are further constrained by comparing the observed spectrum to the library of synthetic spectra published by Munari et al. (2005) and assuming solar metallicity. The metallicity-dependent Ca ii triplet is used as the primary discriminator because of its line strength and best S/N ratio. In addition, to rule out any metallicity effect, the metallicity independent P14 line is used as secondary discriminator.
We construct a grid of synthetic spectra by covering the parameter space , . On the basis of the spectral type of CEN 34, we suggest a rotation velocity of 10 km s-1 (the typical value for early-G supergiants; De Medeiros et al. 2002) during the fitting procedure. The best fitting synthetic spectrum should first match the Ca ii triplet and then match the P14 line; fitting results are checked by visual inspection alone. Figure 4 shows the best fitting grid and other nearby grids. Among the five different grids, we notice that the line strengths of the Ca ii lines are very sensitive to log g, while not varying much with Teff changing; in contrast, the P14 line is sensitive to both log g and Teff. Only the grid of Teff = 5250 K, log g = 2.0, provides an acceptable match for the two Ca ii lines, as well as the P14 line. Because the S/Ns of the two Ca ii lines are much higher than that of the P14 line, log g is better fitted than Teff. Therefore, the uncertainty of log g is estimated to be 0.3, smaller than the grid step, but the uncertainty estimate of Teff is 250 K, similar to the grid step. Through this approach, Teff is determined as 5250 ± 250 K and log g as 2.0 ± 0.3.
3.3. CO bandhead absorption
The CO bandhead absorption longwards of 2.29 μm is typical of late-type stars, whose outer atmospheric layers have the correct temperatures (~1000−3000 K) to produce such features. The CO bandhead absorption generally start to appear in the spectra of late-G dwarfs, but their appearance could shift to earlier spectral types for supergiants (e.g., Wallace & Hinkle 1997). However, low-mass Class II/III YSOs may also exhibit CO bandheads in absorption, because the effective temperature is not high enough to dissociate CO molecules, and the circumstellar veiling is less than that of Class I YSOs (Casali & Eiroa 1996; Hoffmeister et al. 2006).
Figure 5 shows the CO bandhead absorption for CEN 34, as well as a reference mid-G supergiant. Since the CO bandhead absorption is strongly related to the Teff and to the surface gravity, González-Fernández & Negueruela (2012, and reference therein) find a good correlation between the equivalent width of the CO bandhead in absorption (hereafter COEW) and spectral types for giants and supergiants: COEW linearly increases with the decreasing Teff, and is larger for supergiants than for giants of the same Teff. Therefore the CO bandhead can be used to verify the above spectral type and luminosity class of CEN 34. Following Fig. 2 in González-Fernández & Negueruela (2012), the COEW ≈ 11 Å of CEN 34 suggests an early-K giant or mid-G supergiant, and the latter agrees well with the previously classified spectral type and luminosity class, if the involved uncertainty is considered.
Based on the spectral classification and synthetic spectral fitting presented above, CEN 34 has a spectral type of early-/mid-G (Teff = 5250 ± 250 K) and log g = 2.0 ± 0.3, that is compatible with a supergiant.
Fig. 5 Normalized X-shooter spectrum of CEN 34 (black) in the wavelength range 2.28–2.39 μm, covering the CO bandhead absorption. For comparison, the reference spectrum of the mid-G supergiant HIP 107348 (SpT G5 Ib; Wallace & Hinkle 1997) is shown in red. The spectrum of CEN 34 has been smoothed to the spectral resolution (R 3000) of the reference spectrum. |
3.4. Radial velocity relative to LSR (VLSR)
Usually, newly formed stars conserve some kinematic imprint of their natal molecular cloud; i.e., stars having originated in the same cloud are expected to display similar radial velocities (Herbig 1977). Therefore, VLSR measurements are an excellent tool for investigating the membership of stars in star-forming regions or stellar clusters.
In the case of M 17, comparison can be made to the young O- and B-type, whose membership is well established. For many of them, reliable spectral classification has been derived by Hoffmeister et al. (2008), primarily on the basis of low-resolution NIR spectroscopy with ISAAC complemented by higher resolution (R = 6600) optical-red GIRAFFE spectra covering the wavelength range 8200−9400 Å. While the spectra of the B-type stars are featureless in the 8000−9000 Å range, accurate reference VLSR values can be derived for 11 O-type stars with strong Ar I absorption lines at 8620 Å. In addition, the X-shooter spectrum of the B7-type PMS star CEN 24/B 275 also serves as VLSR reference (Ochsendorf et al. 2011). Its Balmer series yield VLSR 24.7 km s-1). In total, we have calculated the VLSR for 12 well-established OB-type members of M 17, and their mean VLSR is +29.4 km s-1, with a scatter of 7.1 km s-1. Furthermore, the prevalent VLSR of the associated molecular clouds (M17-SW and M17-North) as obtained from radio and millimeter observations (Wilson et al. 2003; Xu et al. 2011) is +20 km s-1.
To determine the VLSR of CEN 34, we use the Ca ii triplet and the CO first overtone in the X-shooter spectrum. The lines of the triplet are fitted with Gaussian absorption profiles to derive their central wavelengths (see Table 1). Assuming R = 8800, the velocity resolution of the X-shooter VIS arm is around 34 km s-1. Indeed, the central wavelengths fitted by Gaussian absorption profiles have errors below the velocity resolution. By testing different fitting parameters, the wavelength deviations are about 0.2 Å, corresponding to velocity resolution of 7 km s-1, i.e., 20% of the velocity resolution. Comparison to the rest wavelengths in air of the Ca ii triplet (see Table 1) then yields a VLSR of −57 ± 6 km s-1. To obtain an alternative estimate of its VLSR from the CO first overtone, we cross-correlated CEN 34’s spectrum with the high-resolution (R = 100 000) vacuum rest-frame spectrum of Arcturus (see Fig. 6). For more details on this approach, the reader is referred to Figer et al. (2003). The spectrum of Arcturus was first resampled to match the lower spectral resolution (R = 5500) of the X-shooter NIR arm and then shifted in velocity by −64 km s-1 to overlap with the spectrum of CEN 34. The uncertainties involved in the velocity shift are dominated by the lower spectral resolution of the X-shooter spectrum and as before assumed to be on the order of 20% of the velocity resolution, or 11 km s-1. Obviously, the two independent VLSR estimates for CEN 34 agree reasonably well within the errors. Thus, CEN 34 has a VLSR that is significantly distinct from those of the OB stars and from those of the associated molecular clouds. This suggests that CEN 34 is not a member of M 17, although one could argue that the high-mass stars in M 17 did conserve their original velocity, while stars with lower mass have been kicked out by tidal interactions.
The VLSR of CEN 34 differs from that of any other molecular cloud within the inner Galaxy along the same line-of-sight (l = 15°), by adopting the Galactic rotation model of Reid et al. (2009). A spiral arm lying beyond the Outer Arm in the first Galactic quadrant ~23 kpc from the Galactic center shows an LSR velocity of −15 km s-1 at l = 15° (Dame & Thaddeus 2011). Therefore, it is not possible to estimate a kinematic distance for CEN 34 by its VLSR.
Fig. 6 Normalized CO first overtone bandhead of CEN 34 (black solid) and those of Arcturus (red and blue dotted lines). The red-dotted line presents the rest frame spectrum of Arcturus, while the blue-dotted one is the spectrum shifted by the velocity corresponding to the result from the cross correlation with respect to CEN 34’s spectrum. The rest frame spectrum is shifted by arbitrary amounts along the vertical axis for presentation purposes. |
Multi-wavelength broadband photometry of CEN 34.
3.5. Extinction and luminosity
In Table 2, we collected the broadband photometric results for CEN 34 available from both the literature and the catalogs of recent infrared surveys (DENIS, 2MASS, and Spitzer GLIMPSE) to construct the SED of CEN 34 (see Fig. 7). The SED shows moderate infrared excess, most likely due to the presence of circumstellar material in an envelope and/or a disk. If CEN 34 is an evolved star the circumstellar envelope is a more plausible assumption than a circumstellar disk (Fujii et al. 2002).
We used the simulated annealing approach described in Liu et al. (2012) to model the observed SED and to estimate the mass of circumstellar material. We explored the scenario of a spherical envelope that is commonly adopted to mimic the dust shell surrounding a post-AGB star. The corresponding dust-density distribution is assumed to be a power law ρenv ∝ r−α with α = 2. Our best model fit is shown in Fig. 7. The observed SED at λ ≲ 4 μm, i.e. dominated by the stellar radiation, is reproduced by a star with Teff = 5260 K reddened by AV = 11 mag. The SED longward of 4 μm can be reproduced by a spherical envelope of 9 × 10-5 M⊙, yielding an optical depth τv = 0.29 for the circumstellar material along the line-of-sight. This corresponds to AV = 0.3 mag and is only a small fraction of the total extinction of AV = 11 mag mentioned above. Considering that Hoffmeister et al. (2008) derived a foreground extinction of 2 mag toward M 17, the largest fraction of extinction originates inside the H ii region and behind it, in case CEN 34 is not a cluster member. Therefore, the distance of CEN 34 is at least greater than 2 kpc, the distance of M 17. With this lower limit of CEN 34’s distance, the lower limit of the luminosity is determined as 1600 L⊙. On the other hand, the derived temperature and surface gravity also constrain the absolute magnitude of CEN 34. Using the calibration from (Straizys & Kuriliene 1981) we obtain MV = −4 mag, which transforms into a photometric distance of 2.3 kpc. However, this distance estimate has a large error due to the uncertainty of log g.
Fig. 7 Observed SED of CEN 34 generated from the photometric data listed in Table 2 (diamonds), together with the best-fitting model SED (dashed line) based on a star+envelope system of Menv = 9 × 10-5M⊙. After dereddening (Av = 11 mag), the extinction free model SED (solid line) can be compared to the dereddened fluxes of the reference post-AGB star IRAS 23304+6147 (triangles; Fujii et al. 2002). In addition, the best-fitting stellar blackbody (Teff = 5260 K) is shown as a dotted line. |
Fig. 8 Position of CEN 34 in the Hertzsprung-Russell diagram (HRD) with a well determined effective temperature (5250 ± 250 K) and loosely constrained luminosity (1600 L⊙ ≲ L∗ ≲ 8100 L⊙). Left panel: the position of CEN 34 denoted by the full dot with error bar is compared with the model of massive star formation via disk accretion (Hosokawa et al. 2010). The luminosity of 1600 L⊙ adopted here is the lower limit (denoted by the upward arrow) when the distance of CEN 34 is assumed to be 2 kpc. Right panel: the position of CEN 34 (the rectangle representing the uncertainties of Teff and L∗) is placed against the isochrones of post-main-sequence evolution (Marigo et al. 2008). The lower and upper isochrones (roughly the 100 Myr and 50 Myr ones, respectively) that enclose the position of CEN 34 are highlighted in red, and the corresponding masses are highlighted in blue. |
4. Discussion
In Fig. 8 we outline the location of CEN 34 in the Hertzsprung-Russell diagram (HRD) based on effective temperature and luminosity ranges constrained in Sects. 3.2 and 3.5. Comparison is made to two distinct sets of theoretical evolutionary tracks and isochrones: while in the lefthand panel, we compare to high-mass star formation models of Hosokawa et al. (2010), in the righthand panel we overplot post-main sequence isochrones of Marigo et al. (2008). Obviously both high-mass YSOs and evolved stars are compatible with the observations of CEN 34.
As demonstrated by the numerical simulations of Hosokawa et al. (2010), very young and accreting high-mass YSOs tend to have rather low Teff and large radius, similar to those obtained for CEN 34 in Sect. 3.2. However, if we assume that CEN 34 is associated with the M 17 region and located at a distance of 2 kpc, its Teff ≈ 5250 K and luminosity L∗ ≳ 1600 L⊙ place CEN 34 in an area of the HRD where the numerical models suggest accretion rates in excess of 4 × 10-3M⊙/yr. Nevertheless, accretion rates may vary strongly with the evolution of high-mass YSOs. As a consequence, CEN 34 does not need to maintain such a high accretion rate throughout its entire formation process. Such extremely high accretion rates are usually expected to correlate with the presence of a very massive (~1 M⊙) envelope/disk surrounding the central high-mass YSO, which obviously contradicts to the amount of circumstellar material actually estimated for CEN 34 (~10-4 M⊙; see Sect. 3.5). This discrepancy could be circumvented by a scenario in which the high-mass YSO has already accumulated most of its circumstellar material in a very short time span, so that the circumstellar material becomes optically thin. We speculate that CEN 34 might be in a phase of gravitational contraction, during which it will not gain more material from its surroundings and will move to the left in the HRD.
To solve the VLSR inconsistency between CEN 34 and M 17, one might argue that CEN 34 has been ejected from M 17. Recently, Gvaramadze et al. (2013) have reported that even O-type stars can be ejected from a massive cluster because of a dynamical few-body encounter in the cluster’s core. If we propose that CEN 34 is dynamically kicked o the cluster center, then the time that CEN 34 needs to move to the current position relative to the cluster center can be used to test our speculation. The distance from CEN 34 to the cluster center is about 1 pc, which is the relative distance projected to the sky plane. Assuming 1 pc to the relative distance perpendicular to the sky plane, the VLSR difference of ~80 km s-1 yields a time of ~1.2 × 104 yrs, which is needed for CEN 34 to move to the current position. This traveling time is comparable to the age of CEN 34 if it had accumulated the bulk of its mass with accretion rate >4 × 10-3M⊙/yr within a period of ~103 yrs. A byproduct of this ejection scenario is to place CEN 34 in an isolated environment where feeding through its parental cloud cannot be continued.
Alternatively, besides the fact that the spectrum of CEN 34 shows very similar characteristics to those of post-AGB stars, which typically are late-type supergiants with spectral absorption features as discussed in Sects. 3.1–3.3 (Szczerba et al. 2007), the overall shape of CEN 34’s SED is reminiscent of post-AGB stars. To foster the hypothesis of CEN 34 being a post-AGB star, we compared its dereddened SED (Fig. 7) to that of a reference post-AGB star of spectral type G2 Ia (IRAS 23304+6147) from Fujii et al. (2002), located at a distance (d~4.5 kpc). IRAS 23304+6147 has mid-infrared excess that is slightly stronger than that of CEN 34, most likely due to larger amount of circumstellar material. In view of the similar spectral types and luminosity classes of CEN 34 and IRAS 23304+6147, their comparable optical and NIR fluxes imply a comparable distance, too. Because the optical and NIR fluxes of CEN 34 are still somewhat higher than those of IRAS 23304+6147, the distance of IRAS 23304+6147 (i.e., 4.5 kpc) serves as upper limit for CEN 34’s distance. Considering the most likely range of CEN 34’s distance (2 kpc ≲ d ≲ 4.5 kpc), the corresponding luminosity range is easily derived (1600 L⊙ ≲ L∗ ≲ 8100 L⊙). By comparing this luminosity range and the effective temperature to post-main sequence isochrones of Marigo et al. (2008) (see the right panel of Fig. 8), we constrain the ranges of stellar mass and age to 5−7 M⊙ and 100−50 Myr, respectively. This would suggest that CEN 34 is an evolved intermediate-mass star. While post-AGB stars typically yield prominent far-infrared emissions (e.g., Suárez et al. 2006; Szczerba et al. 2007), far-infrared flux measurements of CEN 34 (e.g., from IRAS or AKARI) are highly contaminated by the strong dust continuum of the M 17 S-bar. For this reason, the 60 μm and 100 μm fluxes listed for CEN 34 in the IRAS Point Source Reject Catalog (released in 1989) are unreliable, because of the large IRAS beam size (~1 arcmin at 60 μm); although AKARI observed with higher spatial resolution and better sensitivities than IRAS, unfortunately the AKARI/FIS All-Sky Survey Point Source Catalogue (Yamamura et al. 2010) does not return an entry matching the position of CEN 34.
5. Summary and conclusions
Originally considered a high-mass YSO candidate in the M 17 star-forming region, CEN 34 has been reclassified as an early-/mid-G supergiant on the basis of the new X-shooter spectrum. The SED of CEN 34 resembles that of either a high-mass YSO or an evolved intermediate-mass star (likely a post-AGB star) with a visual extinction AV = 11 mag that it is dominated by interstellar extinction. If CEN 34 is a high-mass YSO, its derive properties (Teff, log g, and VLSR) can be explained by a scenario where the star has been ejected from the M 17 cluster. In this case CEN 34 had accumulated the bulk of its mass during the first ~103 years, and it is currently gravitationally contracting with only little circumstellar material left. If, on the other hand, CEN 34 is an evolved star of intermediate mass, its luminosity is in the range 1600 L⊙ ≲ L∗ ≲ 8100 L⊙, corresponding to a distance range of 2 kpc ≲ d ≲ 4.5 kpc. Comparison to the isochrones of post-main sequence evolution then suggests that CEN 34 has a mass of 5−7 M⊙ and an age of 50−100 Myr, which are reminiscent of a post-AGB star.
In conclusion, the ambiguities concerning infrared excess and spectral types between high-mass YSOs and evolved intermediate-mass stars prevent a unique classification of CEN 34. The arguments for interpreting CEN 34 as a high-mass YSO are not as straightforward as in the case of a post-AGB star, and they are based on a number of assumptions. For the time being, we prefer the interpretation of a post-AGB star and hope that a future, more accurate distance estimate will settle the issue.
Acknowledgments
We acknowledge the anonymous referee whose comments helped to clarify the paper at several points. The ESO Paranal staff is acknowledged for obtaining the X-shooter spectrum of CEN 34 in Service Mode. M.F. acknowledges the support by the NSFC through grants 11203081. This publication used data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research made use of the SIMBAD database, operated at the CDS, Strasbourg, France. This research made use of the NASA/IPAC Infrared Science Archive (IRSA) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research made use of NASA’s Astrophysics Data System Bibliographic Services.
References
- Andrillat, Y., Jaschek, C., & Jaschek, M. 1995, A&AS, 112, 475 [NASA ADS] [Google Scholar]
- Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953 [NASA ADS] [CrossRef] [Google Scholar]
- Carquillat, M. J., Jaschek, C., Jaschek, M., & Ginestet, N. 1997, A&AS, 123, 5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Casali, M. M., & Eiroa, C. 1996, A&A, 306, 427 [NASA ADS] [Google Scholar]
- Chen, Z., Jiang, Z., Wang, Y., et al. 2012, PASJ, 64, 110 [NASA ADS] [Google Scholar]
- Chini, R., & Krügel, E. 1985, A&A, 146, 175 [NASA ADS] [Google Scholar]
- Chini, R., Hoffmeister, V., Kimeswenger, S., et al. 2004, Nature, 429, 155 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Dame, T. M., & Thaddeus, P. 2011, ApJ, 734, L24 [NASA ADS] [CrossRef] [Google Scholar]
- De Medeiros, J. R., Udry, S., Burki, G., & Mayor, M. 2002, A&A, 395, 97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- D’Odorico, S., Dekker, H., Mazzoleni, R., et al. 2006, in SPIE Conf. Ser., 6269 [Google Scholar]
- Epchtein, N., de Batz, B., Capoani, L., et al. 1997, The Messenger, 87, 27 [NASA ADS] [Google Scholar]
- Figer, D. F., Gilmore, D., Kim, S. S., et al. 2003, ApJ, 599, 1139 [NASA ADS] [CrossRef] [Google Scholar]
- Fujii, T., Nakada, Y., & Parthasarathy, M. 2002, A&A, 385, 884 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Goldoni, P. 2011, Astron. Nachr., 332, 227 [NASA ADS] [CrossRef] [Google Scholar]
- González-Fernández, C., & Negueruela, I. 2012, A&A, 539, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gray, R. O., & Corbally, J. C. 2009, Stellar Spectral Classification (Princeton University Press) [Google Scholar]
- Gvaramadze, V. V., Kniazev, A. Y., Chené, A.-N., & Schnurr, O. 2013, MNRAS, 430, L20 [NASA ADS] [Google Scholar]
- Hanson, M. M., Howarth, I. D., & Conti, P. S. 1997, ApJ, 489, 698 [NASA ADS] [CrossRef] [Google Scholar]
- Herbig, G. H. 1977, ApJ, 214, 747 [NASA ADS] [CrossRef] [Google Scholar]
- Hoffmeister, V. H. 2008, Ph.D. Thesis, Ruhr-Universität Bochum [Google Scholar]
- Hoffmeister, V. H., Chini, R., Scheyda, C. M., et al. 2006, A&A, 457, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hoffmeister, V. H., Chini, R., Scheyda, C. M., et al. 2008, ApJ, 686, 310 [NASA ADS] [CrossRef] [Google Scholar]
- Hosokawa, T., Yorke, H. W., & Omukai, K. 2010, ApJ, 721, 478 [NASA ADS] [CrossRef] [Google Scholar]
- Jiang, Z., Yao, Y., Yang, J., et al. 2002, ApJ, 577, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, Y., Madlener, D., Wolf, S., Wang, H., & Ruge, J. P. 2012, A&A, 546, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marigo, P., Girardi, L., Bressan, A., et al. 2008, A&A, 482, 883 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Munari, U., & Tomasella, L. 1999, A&AS, 137, 521 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Munari, U., Sordo, R., Castelli, F., & Zwitter, T. 2005, A&A, 442, 1127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nielbock, M., Chini, R., Jütte, M., & Manthey, E. 2001, A&A, 377, 273 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nielbock, M., Chini, R., Hoffmeister, V. H., et al. 2007, ApJ, 656, L81 [NASA ADS] [CrossRef] [Google Scholar]
- Nürnberger, D. E. A., Chini, R., Eisenhauer, F., et al. 2007, A&A, 465, 931 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ochsendorf, B. B., Ellerbroek, L. E., Chini, R., et al. 2011, A&A, 536, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Povich, M. S., Churchwell, E., Bieging, J. H., et al. 2009, ApJ, 696, 1278 [NASA ADS] [CrossRef] [Google Scholar]
- Reader, J., Corliss, C. H., Wiese, W. L., & Martin, G. A. 1980, Wavelengths and transition probabilities for atoms and atomic ions: Part 1. Wavelengths, part 2. Transition probabilities [Google Scholar]
- Reid, M. J., Menten, K. M., Zheng, X. W., et al. 2009, ApJ, 700, 137 [NASA ADS] [CrossRef] [Google Scholar]
- Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
- Straizys, V., & Kuriliene, G. 1981, Ap&SS, 80, 353 [NASA ADS] [CrossRef] [Google Scholar]
- Suárez, O., García-Lario, P., Manchado, A., et al. 2006, A&A, 458, 173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Szczerba, R., Siódmiak, N., Stasińska, G., & Borkowski, J. 2007, A&A, 469, 799 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Valdes, F., Gupta, R., Rose, J. A., Singh, H. P., & Bell, D. J. 2004, ApJS, 152, 251 [NASA ADS] [CrossRef] [Google Scholar]
- Vernet, J., Dekker, H., D’Odorico, S., et al. 2011, A&A, 536, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wallace, L., & Hinkle, K. 1997, ApJS, 111, 445 [NASA ADS] [CrossRef] [Google Scholar]
- Wilson, T. L., Hanson, M. M., & Muders, D. 2003, ApJ, 590, 895 [NASA ADS] [CrossRef] [Google Scholar]
- Xu, Y., Moscadelli, L., Reid, M. J., et al. 2011, ApJ, 733, 25 [NASA ADS] [CrossRef] [Google Scholar]
- Yamamura, I., Makiuti, S., Ikeda, N., et al. 2010, VizieR Online Data Catalog: II/298 [Google Scholar]
- Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 The projection of CEN 34 (red cross) on the sky overlaid on the JHKs pseudo-color image of M 17 taken by the InfraRed Survey Facility (IRSF) telescope. The ionizing sources of M 17-the O4-type binary, and the hypercompact H ii region-M 17-UC1 are also marked. |
|
In the text |
Fig. 2 X-shooter spectrum of CEN 34, covering the visible and near-IR wavelength range. Prominent photospheric absorption lines are marked. Gray shaded regions mark the major telluric absorption features. |
|
In the text |
Fig. 3 X-shooter spectrum of CEN 34 (emphasized in red) in the wavelength region 8375–8770 Å. The two features apparently seen in emission in the CEN 34 spectrum at wavelengths ~8590 Å and ~8730 Å are artifacts with unnatural line widths. For comparison, the spectra of several reference stars are displayed; they are adopted from the literature as follows: F8 III (HIP 74975), G1 III (HIP 117503), and G1 II (HIP 47908) from Valdes et al. (2004); G0 Ib (HIP 96481) and G0 III (HIP 5454) from Andrillat et al. (1995); G2 Ib (HIP 109074), G2 III (HIP 86731), G5 III (HIP 100524), and G5 Ib (HIP 107348) from Carquillat et al. (1997). All reference spectra have a resolution around 1 Å, very close to that of our CEN 34 spectrum. All spectra are normalized and are displayed by applying a corresponding shift along the Y-axis for better visual comparison. |
|
In the text |
Fig. 4 CEN 34’s X-shooter spectrum (black solid) around the Ca ii 8498 Å, P14, and Ca ii 8662 Å lines. The best fitting grid (red) and other nearby grids are shown as dashed lines. The synthetic spectral library is based on the Kurucz codes, and has an uniform dispersion of 1 Å/pix (for more details see Munari et al. 2005). |
|
In the text |
Fig. 5 Normalized X-shooter spectrum of CEN 34 (black) in the wavelength range 2.28–2.39 μm, covering the CO bandhead absorption. For comparison, the reference spectrum of the mid-G supergiant HIP 107348 (SpT G5 Ib; Wallace & Hinkle 1997) is shown in red. The spectrum of CEN 34 has been smoothed to the spectral resolution (R 3000) of the reference spectrum. |
|
In the text |
Fig. 6 Normalized CO first overtone bandhead of CEN 34 (black solid) and those of Arcturus (red and blue dotted lines). The red-dotted line presents the rest frame spectrum of Arcturus, while the blue-dotted one is the spectrum shifted by the velocity corresponding to the result from the cross correlation with respect to CEN 34’s spectrum. The rest frame spectrum is shifted by arbitrary amounts along the vertical axis for presentation purposes. |
|
In the text |
Fig. 7 Observed SED of CEN 34 generated from the photometric data listed in Table 2 (diamonds), together with the best-fitting model SED (dashed line) based on a star+envelope system of Menv = 9 × 10-5M⊙. After dereddening (Av = 11 mag), the extinction free model SED (solid line) can be compared to the dereddened fluxes of the reference post-AGB star IRAS 23304+6147 (triangles; Fujii et al. 2002). In addition, the best-fitting stellar blackbody (Teff = 5260 K) is shown as a dotted line. |
|
In the text |
Fig. 8 Position of CEN 34 in the Hertzsprung-Russell diagram (HRD) with a well determined effective temperature (5250 ± 250 K) and loosely constrained luminosity (1600 L⊙ ≲ L∗ ≲ 8100 L⊙). Left panel: the position of CEN 34 denoted by the full dot with error bar is compared with the model of massive star formation via disk accretion (Hosokawa et al. 2010). The luminosity of 1600 L⊙ adopted here is the lower limit (denoted by the upward arrow) when the distance of CEN 34 is assumed to be 2 kpc. Right panel: the position of CEN 34 (the rectangle representing the uncertainties of Teff and L∗) is placed against the isochrones of post-main-sequence evolution (Marigo et al. 2008). The lower and upper isochrones (roughly the 100 Myr and 50 Myr ones, respectively) that enclose the position of CEN 34 are highlighted in red, and the corresponding masses are highlighted in blue. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.