Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.arcetri.astro.it/science/starform/preprints/mosca_06.pdf Äàòà èçìåíåíèÿ: Tue Jan 17 17:36:32 2012 Äàòà èíäåêñèðîâàíèÿ: Mon Feb 4 14:29:23 2013 Êîäèðîâêà:

Massive star-formation toward G28.87+0.07 (IRAS 18411-0338) investigated by means of maser kinematics and radio to infrared, continuum observations
J. J. Li1,2 , L. Moscadelli2 , R. Cesaroni2 , R. S. Furuya3 , Y. Xu1 , T. Usuda3 , K. M. Menten4 , M. Pestalozzi5 , D. Elia5 , and E. Schisano5 Received ; accepted

1

Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China;

jjli@pmo.ac.cn
2 3

INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy Subaru Telescope, National Astronomical Observatory of Japan, 650 North A'ohoku

Place, Hilo, HI 96720, USA
4 5

Max-Planck-Institut fur Radioastronomie, Auf dem Hugel 69, 53121 Bonn, Germany ¨ ¨ INAF-Istituto Fisica Spazio Interplanetario, Via Fosso del Cavaliere 100, I-00133 Roma,

Italy


­2­ ABSTRACT

We used the Very Long Baseline Array (VLBA) and the European VLBI Network (EVN) to perform phase-referenced VLBI observations of the three most powerful maser transitions associated with the high-mass star-forming region G28.87+0.07: the 22.2 GHz H2 O, 6.7 GHz CH3 OH, and 1.665 GHz OH lines. We also performed VLA observations of the radio continuum emission at 1.3 and 3.6 cm and Subaru observations of the continuum emission at 24.5 µm. Two centimeter continuum sources are detected and one of them (named "HMC") is compact and placed at the center of the observed distribution of H2 O, CH3 OH and OH masers. The bipolar distribution of line-of-sight (l.o.s) velocities and the pattern of the proper motions suggest that the water masers are driven by a (proto)stellar jet interacting with the dense circumstellar gas. The same jet could both excite the centimeter continuum source named "HMC" (interpreted as free-free emission from shocked gas) and power the molecular outflow observed at larger scales ­ although one cannot exclude that the free-free continuum is rather originating from a hypercompact H II region. At 24.5 µm, we identify two ob jects separated along the north-south direction, whose absolute positions agree with those of the two VLA continuum sources. We establish that 90% of the luminosity of the region (2â105 L ) is coming from the radio source "HMC", which confirms the existence of an embedded massive young stellar object (MYSO) exciting the masers and possibly still undergoing heavy accretion from the surrounding envelope.
Subject headings: ISM: individual ob jects (G28.87+0.07) ­ ISM: kinematics and dynamics ­ masers ­ techniques: interferometric


­3­ 1. INTRODUCTION

Understanding the process of high-mass star formation represents a challenge from both a theoretical and observational point of view. By now, three main competing concepts of massive star formation have been discussed in the recent literature (see Zinnecker & Yorke 2007, and references therein): formation through coalescence of lower-mass stars in extremely dense stellar clusters (Bonnell et al. 1998); monolithic collapse of isolated, turbulence supported cores (Mckee & Tan 2003); and competitive accretion where the most massive protostars benefit of the gravitational well created by the lower-mass stars surrounding them, thus boosting the accretion rate (e.g., Bonnell & Bate 2006). The available observational evidence is still inadequate to discern between these different scenarios, mostly because observations of massive young stellar ob jects (MYSOs) are complicated by the large distances of typically several kpc, implying small angular scales, and the fact that such ob jects are born still deeply enshrouded in dense dusty envelopes, opaque to optical and near-infrared wavelength radiation. A possible way to distinguish between different models is to explore the physical conditions and the dynamical processes of the gas close to the protostars. However, present (sub)millimeter interferometers and infrared (IR) telescopes attain angular resolutions 0. 25, corresponding to thousands of astronomical units (AU) at the typical distances of MYSOs of several kpc. This is insufficient to investigate the kinematics of the gas in the immediate surroundings of the forming star, where crucial phenomena such as accretion, rotation, and ejection are likely most prominent. The above mentioned limitation can be bypassed by Very Long Baseline Interferometry (VLBI) observations of maser lines, which attain linear resolutions of a few AU. A few years ago, we started an observational pro ject to study the high-mass star-forming process by comparing interferometric images of thermal lines of molecular


­4­ tracers with multi-epoch VLBI studies in three well known maser species (OH, CH3 OH and H2 O). At present, three sources (IRAS 20126+4104, G16.59-0.05 and G23.01-0.41) have been already analyzed (Moscadelli et al. 2011; Sanna et al. 2010a,b). The results demonstrate that the synergy between VLBI, multi-species maser observations and (sub)mm interferometric observations of thermal molecular lines is crucial to achieve a multi-scale picture of the environment of newly formed massive (proto)stars. In the present article we report on our VLBI maser study of the high-mass star-forming region (HMSFR) G28.87+0.07, and compare our results with available complementary interferometric data of continuum and thermally excited molecular line emission, as well as with sub-arcsec resolution mid-IR images obtained with the Cooled Mid-Infrared Camera and Spectrometer (COMICS) on the Subaru telescope. In Section 2, we provide a review of previous observations toward this region. Section 3 describes the VLBI observations of the maser transitions, the Very Large Array (VLA) observations of the radio continuum emission at 1.3 cm and 3.6 cm, and the mid-IR observations performed with the Subaru/COMICS telescope. In Section 4, we present the distributions of the three types of masers, the water maser proper motions, and the properties of the radio continuum. The results are discussed in Section 5 and the conclusions are drawn in Section 6.

2.

THE HMSFR G28.87+0.07

The HMSFR G28.87+0.07 (alias IRAS 18411-0338) is located at a kinematic distance of 7.4 kpc1 (Codella et al. 1997; Furuya et al. 2008) and is associated with a pc-scale clump
1

Since the ambiguity between the near (6 kpc) and the far (9 kpc) kinematical distance

is not resolved (Sewilo et al. 2004), we adopt the distance to the tangent point.


­5­ with a mass in excess of 103 M and a bolometric luminosity of 2â105 L (L´ ez-Sepulcre op et al. 2009; Faundez et al. 2004). From the observation of NH3 , the systemic velocity of the ´ region (Vsys ) with respect to the local standard of rest (LSR) is found to be 101.0 km s- (Codella et al. 1997). On a large scale, mid-IR observations at 8, 10.5, and 18.1 µm have resolved the IRAS source in two ob jects along the north-south direction (De Buizer et al. 2005; Furuya et al. 2008, Figure 1). Both sources could be associated with an Midcourse Space Experiment (MSX) point source (Egan et al. 2003) and a 2.2 µm source with strong near-infrared excess, lying 3 ­4 to the southeast of the northern, brighter mid-IR ob ject (Testi et al. 1994). Bolometer observations at 1.2 mm (Faundez et al. 2004) have detected a peak ´ of dust continuum emission associated with the two mid-IR sources and the near-IR source. Single-dish maps in the CS(5­4) and C18 O(2­1) lines have revealed an arcmin-scale molecular clump and a monopolar (red-shifted) outflow is seen in the
13 1

CO(2­1) line

(Shirley et al. 2003; L´ ez-Sepulcre et al. 2009). Single-dish surveys have detected emission op in dense gas tracers as CS(7­6), CS(2­1), and CS(1­0) (Plume et al. 1992; Bronfman et al. 1996; Anglada et al. 1996), and in the H110 radio recombination line (Sewilo et al. 2004; Araya et al. 2007), which suggests the presence of a diffuse H II region. In addition, OH maser emission has been found in the main (1.665 and 1.667 GHz) hyperfine structure lines and in the 1.720 GHz satellite line, whereas the 1.612 GHz OH satellite has been observed in absorption (Little & Cesarsky 1982; Caswell & Haynes 1983). Weak CH3 OH and strong H2 O maser emission have been reported at 6.7 GHz (Caswell et al. 1995) and 22.2 GHz (Genzel & Downes 1979), respectively. Interferometric observations at millimeter wavelengths have revealed multiple sources in this region (Furuya et al. 2008, Figure 1) and show that the dust continuum peak coincides with the position of the northern, stronger mid-IR source. The mass (100 M ) from the


­6­ 3 mm dust continuum emission and the rotational temperature (93 K) inferred from the CH3 CN lines suggest that one is dealing with a hot molecular core (HMC). In this study, we will focus on this HMC, as it contains all the tracers of high-mass star formation that we have investigated with our observations, namely the masers and the free-free and mid-IR continuum sources. As shown by Codella et al. (1997) and Furuya et al. (2008), the HMC is seen in NH3 , HCO+ , and CH3 CN line emission. The observations of Codella et al. (1997) were not sensitive enough to detect continuum emission at 1.3 cm, whereas our new images demonstrate the existence of a free-free continuum source inside the HMC (see below). A wide-angle, massive (100 M ) bipolar outflow has been detected in the
12

CO(1­0) line,

oriented perpendicular to a CH3 CN velocity gradient detected in the HMC (Furuya et al. 2008). The absolute positions of the 1.665 GHz OH and 22.2 GHz H2 O masers, derived by Forster & Caswell (1999) using the VLA with an accuracy of 0. 5, indicate that these maser emissions are associated with the HMC. All these features are consistent with the presence of a deeply embedded MYSO(s) inside the HMC, possibly accompanied by lower mass ob jects lying in the neighborhoods of the HMC itself.

3.

OBSERVATION AND DATA REDUCTION 3.1. VLA: 1.3 cm and 3.6 cm Continuum

We observed the HMSFR G28.87+0.07 (tracking center: R.A.(J2000) = 18h 43m 46.24s and Dec.(J2000) = -03 35 30.4 ) with the VLA at X (3.6 cm) and K (1.3 cm) bands in both the C-array (in 2008 April) and A-array configurations (in 2008 November). At 3.6 cm in both configurations and at 1.3 cm in the C-array configuration, the continuum mode of the correlator was used, resulting in an effective bandwidth of 172 MHz. At 1.3 cm in


­7­ the A-array we used mode "4" of the correlator, with a pair of intermediate frequency (IF) bands, each with a width of 3.125 MHz (64 channels) centered on the strongest H2 O maser line and a pair of 25 MHz wide IF bands (8 channels each) sufficiently offset from the maser lines to obtain a measurement of the continuum emission. The two bands centered at the same frequency (detecting both circular polarizations) were averaged to decrease the noise. At X-band, 3C 286 (5.2 Jy) and 3C 48 (3.1 Jy) were used as primary flux calibrators, while 1832-105 (1.4 Jy) was the phase calibrator. For the K-band observations, the primary flux calibrator was 3C 286 (2.5 Jy), the phase calibrator 1832-105 (1.0 Jy), and the bandpass calibrator (for the A-array only) 1733-130 (3.7 Jy). The data were calibrated with the NRAO2 Astronomical Image Processing System (AIPS) software package using standard procedures. Only for the A-array data at 1.3 cm, several cycles of self-calibration were applied to the strongest maser channel, and the resulting phase and amplitude corrections were eventually transferred to all the other line channels and to the K-band continuum data (Menten & Reid 1995). This procedure resulted in a significant (at least a factor 2) improvement in the signal-to-noise ratio (SNR). Natural-weighted maps were made with task IMAGR of AIPS for both the continuum and the line data. Simultaneous observation of the line and continuum emission at K-band made it possible to obtain the relative position of the continuum image with respect to the H2 O maser spots with great precision (0. 1). Finally, we identified the maser spots that appear in both our Very Long Baseline Array (VLBA) and VLA images and, using the accurate absolute position information obtained from the VLBA observations, we re-centered the VLA spots and continuum image. The correction thus applied turns out to
2

The National Radio Astronomy Observatory is operated by Associated Universities, Inc.,

under cooperative agreement with the National Science Foundation.


­8­ be <20 mas in each coordinate. We conclude that the absolute astrometrical precision of the continuum maps presented in this paper is 0. 1.

3.2.

VLBI: Maser Lines

We conducted VLBI observations of the H2 O and CH3 OH masers (at four epochs) and of the OH maser (at a single epoch) toward G28.87+0.07 at 22.2 GHz, 6.7 GHz, and 1.665 GHz, respectively. We give more technical details on the observations in the following subsections. In order to determine the masers' absolute positions, we performed phase-referencing observations in fast switching mode between the maser source and the calibrator J1834-0301. This calibrator has an angular offset from the maser source of 2 4 . and belongs to the list of sources of the VLBA Calibrator Survey 3 (VCS3). Its absolute position is known with a precision better than ±1.5 mas, and its flux measured with the VLBA at S and X bands is 154 and 156 mJy beam-1 , respectively. Five fringe finders (J1642+3948; J1751+0939; J1800+3848; J2101+0341; J2253+1608) were observed for bandpass, single-band delay, and instrumental phase-offset calibration. The data were reduced with AIPS following the VLBI spectral line procedures. For a description of the general data calibration and the criteria used to identify individual masing clouds, derive their parameters (position, intensity, flux and size), and measure their (relative and absolute) proper motions, we refer to the recent paper on VLBI observations of H2 O and CH3 OH masers by Sanna et al. (2010a). For each VLBI epoch, H2 O maser absolute positions have been derived by employing two methods, direct phase-referencing and inverse phase-referencing, and the two procedures always gave consistent results. CH3 OH and OH maser absolute positions have been derived only by using direct phase-referencing. The derived absolute proper motions have been corrected for the apparent proper motion due to the Earth's orbit around the Sun (parallax), the solar motion and the differential Galactic


­9­ rotation between our LSR and that of the maser source. We have adopted a flat Galactic rotation curve (R0 = 8.3 ± 0.23 kpc, 0 = 239 ± 7 km s-1 ) (Brunthaler et al. 2011), and the solar motion of U = 11.1+0..69 , V = 12.24+0..47 , and W = 7.25+0 -0 75 -0 47 -0 et al. (2010), who recently revised the Hipparcos satellite results.
.37 .36

km s

-1

by Sch¨ onrich

3.2.1.

VLBA: 22.2 GHz H2 O Masers

We observed the HMSFR G28.87+0.07 (tracking center: R.A.(J2000) = 18h 43m 46.22s and Dec.(J2000) = -03 35 29.9 ) with the VLBA in the 616 - 523 H2 O transition (rest frequency 22.235079 GHz). The observations (program code: BM244) consisted of 4 epochs: 2006 April 23, 2006 June 30, 2006 September 28, and 2007 January 18. During a run of about 6 h per epoch, we recorded the dual circular polarization with 16 MHz bandwidth centered at an LSR velocity (VLSR ) of 101.0 km s-1 . The data were processed with the VLBA FX correlator in Socorro (New Mexico) using an averaging time of 1 s and 1024 spectral channels. Images were produced with natural weighting, cleaned and restored with an elliptical Gaussian with full width at half maximum (FWHM) size of about 1.3 mas â 0.6 mas at a position angle (PA) of -4 (east of north), with small variation from epoch to epoch (shown in Table 2). The interferometer instantaneous field of view was limited to 2. 7. In each observing epoch the on-source integration time was 2.5 h resulting in an effective rms noise level of the channel maps ( ) in the range 0.006­0.015 Jy beam-1 . The spectral resolution was 0.2 km s-1 .


­ 10 ­ 3.2.2. EVN: 6.7 GHz CH3 OH Masers We observed G28.87+0.07 (tracking center: R.A.(J2000) = 18h 43m 46.22s and Dec.(J2000) = -03 35 29.9 ) with the European VLBI Network (EVN) in the 51 - 60 A
+

CH3 OH transition (rest frequency 6.668519 GHz). Data were taken in 4 epochs (program code: EM061 and EM069) separated by 1 year: 2006 February 28, 2007 March 18, 2008 March 18, and 2009 March 17. In the first two epochs, the antennae involved in the observations were Cambridge, Jodrell2, Effelsberg, Hartebeesthoek, Medicina, Noto, Torun, and Westerbork. Since we found out that the longest baselines involving the Hartebeesthoek antenna (e.g., the Ef-Hh baseline is 8042 km) heavily resolve the maser emission and do not produce fringe-fit solutions, the Hartebeesthoek antenna was replaced with the Onsala antenna in the third epoch. During a run of 6 h, we recorded the dual circular polarization in two IF bands of width of 2 MHz and 16 MHz, both centered at an LSR velocity of 101.0 km s-1 . The 16 MHz band was used to increase the SNR of the weak continuum calibrator. The data were processed with the MKIV correlator at the Joint Institute for VLBI in Europe (JIVE-Dwingeloo, the Netherlands) using an averaging time of 1 s and 1024 spectral channels for each observing bandwidth. Images were produced with natural weighting, cleaned and restored with an elliptical Gaussian with FWHM size of about 12 mas â 6 mas at a PA of 21 , with small variations from epoch to epoch (shown in Table 2). The interferometer instantaneous field of view was limited to 9. 2. Using an on-source integration time of 2.2 h, the effective rms noise level of the channel map ( ) varied in the range 0.005­0.01 Jy beam-1 . The 2 MHz band spectral resolution was 0.09 km s-1 .


­ 11 ­ 3.2.3. VLBA: 1.665 GHz OH Masers

We observed G28.87+0.07 (tracking center: R.A.(J2000) = 18h 43m 46.34s and Dec.(J2000) = -03 35 29.9 ) with the VLBA in the 2
3/2

J = 3/2 OH transition (rest

frequency 1.665401 GHz) on 2007 April 28 (program code: BM244O). During a run of 6 h, we recorded the dual circular polarization with two IF bands of width of 1 MHz and 4 MHz, both centered at VLSR = 103.0 km s-1 . The 4 MHz bandwidth was used to increase the SNR of the weak L-band signal of the continuum calibrator. The data were processed with the VLBA FX correlator in two correlation passes using 1024 and 128 spectral channels for the 1 MHz and 4 MHz bands, respectively. In each correlator pass, the data averaging time was 2 s. Images were produced with natural weighting, cleaned and restored with an elliptical Gaussian with FWHM size of 17.8 mas â 10.9 mas at a PA of 1 . The interferometer instantaneous field of view was limited to 18. 5. The on-source integration time was 1.9 h, resulting in an effective rms noise level in each velocity channel of 0.02 Jy beam-1 . The spectral resolution was 0.2 km s-1 .

3.3.

Subaru/COMICS: 24.5 µm Mid-infrared Continuum

Using the mid-infrared imaging spectrometer (COMICS; Kataza et al. 2000) at the Cassegrain focus of the 8.2 m Subaru Telescope, we carried out imaging observations of the 24.5 µm emission toward G28.87+0.07 on 2008 July 15. For this purpose, we configured the Q24.5 filter, and employed chopping mode for subtracting the sky-background emission. The camera provides a field of view of 42 â 32 with a pixel size of 0. 13. Flux calibration was performed towards four standard sources listed in Cohen et al. (1999): HD146051, HD186791, HD198542, and HD198542. We estimated the overall uncertainty in the flux


­ 12 ­ calibrations to be less than 10%. Astrometric calibration was performed by smoothing the 24.5 µm image to the same angular resolution as the 24 µm MIPSGAL (survey of the inner Galactic plane using MIPS; Carey et al. 2009) image, i.e. 6 , and overlaying the two images until a good match was found. Note that the heavy saturation of the MIPSGAL image does not hinder the comparison between it and the smoothed COMICS image. In practice, the emission in MIPSGAL is basically pointlike and the MIPS PSF presents a well defined circular structure even far from the (saturated) peak. The match between the smoothed COMICS image and the MIPSGAL image can thus be performed successfully on such circular features. Inspection by eye indicates a positional error of less than 1 , comparable to the astrometric accuracy of MIPSGAL.

4.

RESULTS AND ANALYSIS

This section reports our main results from the observation of the radio continuum, individual maser species and infrared emission towards the HMSFR G28.87+0.07. In the following, the term "spot" is used to refer to maser emission in a single velocity channel, whereas the term "feature" indicates a collection of spots emitting at a similar position over contiguous channels (i.e. an individual masing-cloud). In this work, the parameters of maser emission (position, intensity, proper motion, etc..) are derived for maser features.

4.1.

1.3 cm and 3.6 cm Continuum

Our VLA measurements of the 1.3 and 3.6 cm continuum emission towards G28.87+0.07 (in both the A- and C-array configurations) are summarized in Table 1, and shown in Figure 1. At the angular resolution of the VLA C-array, the emission consists of two components, both detected at 1.3 and 3.6 cm. The northern and southern components,


­ 13 ­ separated by 3. 1, are labeled, respectively, "HMC" and "A" in Table 1. The former name is due to the source being associated with the hot molecular core imaged by Furuya et al. (2008). Source "A" is partly resolved with the C-array at 1.3 cm, and is completely resolved out with the A-array at both 1.3 and 3.6 cm. In the following, we focus our attention on the centimeter continuum source "HMC", because this appears to be spatially associated with the observed 22 GHz H2 O, 6.7 GHz CH3 OH, and 1.665 GHz OH masers in G28.87+0.07 (see Figures 1 and 4). While at 1.3 cm the most extended VLA A-array configuration partly resolves the emission of "HMC" (missing about half of the flux measured with the C-array), at 3.6 cm the emission is essentially compact in both configurations (yielding comparable integral fluxes). Note however that the VLA A-array image at 3.6 cm shows a spur protruding to the south. The peak positions of the 1.3 and 3.6 cm emissions coincide within the uncertainties. The absolute position of the A-array 1.3 cm image has been obtained by aligning positions of persistent water maser spots between the VLA and VLBA observations, and should be accurate to 0. 1. Figure 2 shows the spectral energy distribution (SED) at centimeter wavelengths of "HMC". A spectral index () of 0.6 is derived from a linear fit of the VLA C-array 1.3 and 3.6 cm and B-array 6 cm flux densities. The latter flux is an upper limit obtained from the GPS 6 cm Epoch 2 images downloadable from the MAGPIS web site3 , which result to be the most sensitive among the other GPS maps and the relevant map from the CORNISH survey 4 (see Purcell et al. 2008). We also fit the SED at centimeter wavelengths with a model of a homogeneous H II region, with a Lyman continuum rate of 7.9 â 1045 s and a radius of 0. 04.
3 4 -1

http://third.ucllnl.org/gps available on http://www.ast.leeds.ac.uk/Cornish/public


­ 14 ­ We caution that the shape of the SED could be significantly affected by the different sensitivity to extended structures at different wavelengths. For the same VLA configuration, the shorter the wavelength, the larger can be the fraction of emission resolved out by the interferometer. Therefore, one should consider the possibility that the slope between 3.6 cm and 1.3 cm is steeper than 0.6.

4.2.

Masers Results

4.2.1. 22.2 GHz Water Maser Forster & Caswell (1999) observed the HMSFR G28.87+0.07 with the VLA C-array (HPBW of 2. 2 â 1. 6) and identified 13 water maser spots, mainly grouped in two clusters separated by 0. 2 in the northeast-southwest direction. The absolute position of the brightest spot (with a flux density of 61.3 Jy at the VLSR of 106.3 km s-1 ) differs by 0. 1 from that of the corresponding spot in our VLBA observation. The field of view of our VLBA observations covers the whole area within which water maser emission was detected by Forster & Caswell (1999). Figure 3 (upper panel) shows the total-power spectrum of the 22.2 GHz masers toward G28.87+0.07. This profile was obtained by averaging the total-power spectra of all VLBA antennae, after weighting each spectrum with the antenna system temperature (Tsys ). Water maser absolute positions and LSR velocities are plotted in Figure 4. Imaging the range of LSR velocities from 93.2 to 114.3 km s-1 , we detected 168 distinct water maser features distributed over a region of about 0. 3 â 0. 5. Most (92%) of them concentrate in a small area of about 0. 3 â 0. 2 centered on the VLA 1.3 cm continuum peak. The line-of-sight (l.o.s.) velocity distribution is remarkably bipolar, with red-shifted features located to the west of the radio continuum and blue-shifted ones to the east. While


­ 15 ­ the blue-shifted features mainly concentrate in two small (size <20 mas) clusters, the distribution of red-shifted masers consists of one narrow strip elongated east-west (size of 100 mas), together with two diffuse "clouds" of features spread to the north and to the south of the densely populated, narrow strip. It is also interesting to note the presence of two streamlines of water features (enclosed in dashed boxes in Figure 4 upper panel), oriented north-south and extending over a few tenths of an arcsecond, on either side of the VLA 1.3 cm continuum. We took particular care in filtering out spurious emission owing to poor cleaning of the beam side lobes (elongated mostly along the north-south direction) and can confirm that these water maser streamlines are real. The individual feature properties are presented in Table 3. The position offsets are relative to feature #1. The absolute positions of feature #1 at the four observing epochs are given in Table 2. Feature intensities range from 0.1 to 54 Jy. The spread in LSR velocities ranges from 111.99 km s-1 for the most red-shifted feature (#107) to 93.33 km s
-1

for the most blue-shifted one (#155). Only 61 features (36% of the total) persisted over at least 3 epochs, 34 of which lasting

4 epochs. We calculated the geometric center (hereafter "center of motion", identified with label #0 in Table 3, and indicated by a star in Figure 5) of features with a stable spatial and spectral structure, persisting over the 4 observing epochs, and refer our measurement of proper motions to this point. With a time baseline of 9 months, relative (sky-pro jected) velocities of water maser features are derived with a mean uncertainty of 23%. The magnitude of relative proper motions ranges from 1.9±2.8 km s-1 to 58.3±2.6 km s-1 , with a mean value of 24.0 km s-1 . Figure 5 shows the derived water maser proper motions. While the two blue-shifted clusters move to the east and northeast, most velocities of the red-shifted water masers point either to the south or to the west. In particular, looking at the distribution of the maser velocity vectors along the (east-west) narrow strip outlined by


­ 16 ­ the red-shifted features, we note that maser features closer to the "center of motion" move mainly to the south. Getting closer to the western end of the maser strip, the orientation of the maser velocity vectors changes gradually from south to west. The general pattern of water maser relative velocities appears to indicate expansion from a point on the sky close to the "center of motion". Note also that at larger distances from the "center of motion", maser velocities seem to be larger and better collimated along the (north)east-(south)west direction. In particular, we identify two H2 O maser clusters located at the northeast and southwest ends of the maser distribution (named "Jb " and "Jr ", respectively; see Figure 5), which move with high velocities directed close to the northeast-southwest direction. Figure 6 (upper panel) reports the intrinsic absolute proper motions of the water maser features, derived from the observed absolute proper motion of feature #1 corrected for the apparent motion. The latter has been calculated adding the contribution of the parallax, the solar motion with respect to the LSR (Sch¨ onrich et al. 2010), and the differential Galactic rotation, from a flat rotation curve with Galactic constants (R0 = 8.3 ± 0.23 kpc, 0 = 239 ± 7 km s-1 ) as recently determined by Brunthaler et al. (2011). The obtained absolute proper motions differ from the relative proper motions presented in Figure 5 by a vector of amplitude 40 km s-1 pointing to the southwest, which seems to be the dominant component of all proper motions. We think that this difference could be caused by the large systematic error due to the uncertainties in the source distance and solar and Galactic motion, and/or by a peculiar velocity of the maser source with respect to its LSR reference system. Our VLBA water maser observations were not designed to measure the source parallax. With a positional accuracy of 0.1­0.2 milli-arcseconds, and four epochs spanning (only) 9 months, the observed absolute proper motion of feature #1, not corrected for the apparent motion, looks similar to a straight line. From the small deviations from linear motion, we derive an upper limit of 0.2 mas for the amplitude of the sinusoidal parallax signature. This corresponds to a lower limit of 5 kpc for the source distance, in agreement


­ 17 ­ with the adopted kinematical distance of 7.4 kpc. We have checked the reliability of the relative, water maser proper motions by following an alternative approach. As described in Section. 4.2.2, we have also derived the absolute proper motion of the strongest 6.7 GHz methanol maser feature. To date, VLBI observations of 6.7 GHz (Sanna et al. 2010a,b) and 12 GHz (Moscadelli et al. 2003, 2010, 2011) methanol masers toward several sources indicate that these masers usually show small proper motions, with typical velocities of only a few km s-1 . Assuming that also in G28.87+0.07 methanol masers move significantly slower than the water masers, we can use the absolute proper motion of the methanol masers to correct the absolute proper motion of the water masers, and derive the water maser velocities as approximately seen in a reference system comoving with the star. Figure 6 (lower panel) shows the absolute water maser proper motions after subtracting the absolute methanol maser proper motion (evaluated in Section. 4.2.2). Note that, within the measurement errors, the new absolute proper motions are fully consistent with the relative proper motions calculated with respect to the "center of motion" (Figure 5). This result reinforces our assertion that the "center of motion" can define a suitable reference system to calculate water maser velocities. In the following discussion, we use the relative velocities of the water masers to describe the gas kinematics close to the (proto)star.

4.2.2.

6.7 GHz Methanol Maser

Figure 3 (middle panel) shows the total-power spectrum of the 6.7 GHz methanol masers observed on 2008 March 18 toward G28.87+0.07 using the Effelsberg 100-m telescope. The 6.7 GHz methanol maser spectrum, with two narrow features at VLSR of about 97 and 105 km s-1 , resembles that of the water masers (Figure 3, upper panel). Note that Caswell et al. (1995) detected only the emission at VLSR of 105 km s-1 with the Parkes


­ 18 ­ 64-m telescope. We imaged the range of LSR velocities from 95 to 107 km s-1 with the EVN and detected two distinct methanol maser features, both with VLSR close to 105 km s-1 . Feature #1 is detected at the first, third, and fourth epoch, whereas feature #2 is detected at the third and fourth epoch only. Due to adverse atmospherical conditions and telescope failures, no 6.7 GHz maser signal was detected at the second epoch. Individual feature properties are presented in Table 4. Positional offsets are relative to feature #1, and the absolute positions of feature #1 at the three epochs of detection are given in Table 2. The two methanol maser features lie at the western border of the distribution of the red-shifted water masers, 0. 2 to the southwest of the VLA 3.6 cm continuum peak (see Figure 4). Inspecting our images, we could not find any reliable emission at the velocity of 97 km s-1 , in correspondence of the weakest feature detected in the total-power spectrum. The mean rms noise on the channel maps was 20 mJy beam-1 , and we could have missed that emission if falling below the 5 detection threshold. Because of the narrow line width (1 km s-1 ) and the offset in V
LSR

from the systemic velocity (101.0 km s-1 ) of the HMC,

we tend to exclude that the weak methanol feature at VLSR 97 km s-1 is thermal emission from the HMC resolved out by the EVN observations. After correcting for the apparent motion, the absolute proper motion for feature #1 is: V = -27.8 ± 1.0, V = -45.7 ± 4.6 km s-1 . As discussed in Sect. 4.2.1, we have used this vector to correct the absolute proper motions of the water masers. The relative motion of the methanol maser feature #1 with respect to the "center of motion" of the water masers is found to be: v = 1.7 ± 6.7, v = -15.6 ± 8.0 km s-1 . This is consistent with the results from other VLBI observations which show that 6.7 GHz masers move significantly slower than water masers, and justifies employing the methanol maser absolute motion to correct the water maser absolute motions.


­ 19 ­ 4.2.3. 1.665 GHz OH Maser Forster & Caswell (1999) observed the HMSFR G28.87+0.07 with the VLA A-B hybrid configuration (HPBW of 5. 5 â 1. 2) in the left circular polarization (LCP) band and detected eight 1.665 GHz OH maser spots. Their brightest spot had a flux density of 9.85 Jy at the VLSR of 102.8 km s-1 . Figure 3 (lower panel) shows the Tsys -weighted mean of the VLBA antenna total-power spectra of the 1.665 GHz OH maser emission in the LCP and RCP bands. Three strong, narrow features are visible in the LCP band over the LSR velocity range from 102 to 108 km s-1 . Two weak features are detected in the RCP band, one of them at an LSR velocity of 108 km s-1 . These findings are similar to the results obtained by Caswell & Haynes (1983) with the Parkes 64-m telescope. We imaged the range of LSR velocities from 102 to 110 km s-1 , and identified 3 distinct OH maser features in the LCP band. In the RCP band, by imaging maser data after direct phase-referencing to the calibrator J1834-0301, no signal was detected above a 5 level of 0.1 Jy beam-1 . Individual feature properties are presented in Table 5. The positional offsets are relative to feature #1, whose absolute position is given in Table 2. The derived absolute position is offset by 0. 9 from the reference position of Forster & Caswell (1999). Maser intensities range from 0.7 to 1.9 Jy beam-1 . The LSR velocities vary from 102.8 km s
-1

for the most blue-shifted feature (#1), to 105.5 km s-1 for the most red-shifted feature

(#3). The FWHM line width of individual maser features ranges from 0.4 to 1.4 km s-1 . Positions and LSR velocities of hydroxy1 maser features are plotted in Figure 4. The three detected features lie 0. 6 to the southeast of the VLA 3.6 cm continuum source and are distributed along a line extended 0. 3.


­ 20 ­ 4.3. 24.5 µm Infrared Emission

Figure 7 shows an image of the 24.5 µm continuum emission obtained with Subaru/COMICS at an angular resolution of 0. 75. One can see that the structure of the emitting region is quite complex. Beside two components associated with the radio sources "HMC" and "A" with a total flux of 80 Jy, another fainter source (3 Jy) is visible to the south, offset by 4 from "A". The emission from the northern region is clearly resolved and appears to be roughly associated with the two radio ob jects, although one cannot rule out the existence of additional, fainter unresolved sources contributing to the extended halo around the two main sources. In particular, the 24.5 µm emission from "HMC" seems to split into two sub-components in the E-W direction, with comparable flux densities: 32 Jy for the eastern and 43 Jy for the western component. These could trace two distinct ob jects or diffuse emission surrounding a single source at the origin of the "HMC" radio continuum. A third, more intriguing possibility is that one is looking at an "hour-glass" shaped nebulosity centered on the embedded source at the origin of the free-free emission. In the latter scenario, the source lies at the center of a flattened core lying approximately along the line of sight. This dusty core absorbs most of the IR photons emitted at 24 µm towards the observer, while the IR radiation can escape along the minor axis of the core, lying close to the plane of the sky. We have estimated an optical depth of 4 for the HMC at 24 µm, assuming the H2 column density of 4.2 â 1023 cm
-2

obtained by Furuya et al.

(2008) and the opacity adopted by Whitney et al. (2003a) ­ see their Figure 1. Such a value could be sufficiently large to explain the lack of emission towards the embedded star and determine the "hour-glass" shape of the IR emission. Interestingly, the E­W orientation of this putative bipolar nebula is also roughly consistent with the direction of the H2 O maser jet in Figure 5, suggesting that one could be observing a bipolar flow piercing through a


­ 21 ­ dense core, and revealed on different scales in different tracers. However, such a scenario, albeit appealing, is not consistent with the inclination of the maser jet with respect to the plane of the sky. In fact, the blue-shifted emission lies to the E (see Figure 5) and this implies that the eastern lobe should be pointing towards the observer and thus be brighter than the western lobe in the IR. Looking at Figure 7 one can see that this is not the case, being the emission slightly stronger to the west (see above). In conclusion, we consider both the "hour-glass" nebula and the double source hypotheses too speculative on the basis of the observational evidence collected so far, and in the following we make the conservative assumption that the IR emission is associated with "HMC", with a flux of 75 Jy. The southern source "A" has a 24.5 µm flux density of 5 Jy. Also in this case, the association between the IR emission and the radio source is questionable, because the two are offset each other by 1 . However, this is the astrometrical error of the IR image and we thus believe that such an offset is not sufficient evidence against a common origin for the IR and radio emission toward "A". Therefore, we conclude that both emissions can originate from the same ob ject.

4.4.

Sp ectral Energy Distribution

Figure 8 shows the SED of G28.87+0.07 from 3.6 µm to 1.1 mm. Beside our COMICS data, we have used data from the following surveys: Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE; Benjamin et al. 2003), MIPSGAL (Carey et al. 2009), MSX (Egan et al. 2003), IRAS (Neugebauer et al. 1984), the Herschel infrared Galactic Plane Survey (Hi-GAL; Molinari et al. 2010a,b), the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL; Schuller et al. 2009), and the Bolocam Galactic Plane Survey


­ 22 ­ (BGPS; Drosback et al. 2008). For the Herschel/Hi-GAL data, the most recent release (obtained in July 2011 with HIPE version 7.0.0) was used, resulting from data reduction with the RomaGal software (Traficante et al. 2011), where the most relevant image artifacts have been removed. Apart from the IRAS and MSX flux densities, which have been taken from the corresponding Point Source Catalogues (sources IRAS 18411-0338 and G028.8621+00.0657), all the other values have been calculated by integrating the emission inside the dashed contour in Figure 9. The latter was chosen in such a way to include the whole emitting region at all wavelengths. In this way we overcome the problem of different angular resolutions at different wavelengths. In the calculation of the flux the background emission was subtracted assuming that the background flux inside the contour is in all points equal to the typical flux density measured along the contour itself. The flux densities used for the SED are given in Table 6. Figure 9 illustrates the pc-scale structure of the G28.87+0.07 star forming region, imaged with Herschel from 70 to 500 µm and with APEX at 870 µm, and demonstrates that from the mid-IR to the sub-mm the continuum emission is dominated by a compact source surrounded by a weaker halo. We have fitted the SED with the radiative transfer model developed by Robitaille et al. (2007) and Whitney et al. (2003a,b) (using the SED fitting tool available on http://caravan.astro.wisc.edu/protostars/), which assumes a ZAMS star with a circumstellar disk, embedded in an infalling envelope, and allows derivation of a number of physical parameters. While we are aware that this model cannot adequately reproduce a complex situation such as that of the multiple MYSOs in G28.87+0.07, we believe that the fit is sufficiently good to guarantee good estimates of integrated quantities such as the luminosity and the mass of the envelope. The best-fit values are, respectively, 2.3â105 L ,


­ 23 ­ and 3.0â103 M . Note that in the fitting procedure, the IRAS fluxes were considered upper limits due to the limited angular resolution, whereas the MIPSGAL 24 µm flux is in fact a lower limit because the corresponding image is saturated and only the non saturated pixels have been considered in the calculation. Also, the distance was kept fixed to the nominal value of 7.4 kpc. The best-fit value for the visual foreground extinction is 46 mag.

5. 5.1.

DISCUSSION

Nature of the Radio Continuum

In this section, we focus on the cm-continuum source labeled "HMC" in Table 1, found in good correspondence with the positions of the 22 GHz H2 O, 6.7 GHz CH3 OH, and 1.665 GHz OH masers in G28.87+0.07 (see Figure 4). The continuum spectrum (shown in Figure 2) is consistent with free-free emission from an ultracompact (UC) H II region, with the turnover frequency in the range 10­20 GHz, and ionized by a Lyman continuum rate of 7.9â1045 s-1 . This would correspond to a ZAMS star of spectral type B1­B0.5 (Panagia 1973), with a luminosity of L
star

8 â 103 L , and a mass of M

star

10 M . It is worth

mentioning that these values would be different if more recent studies than the seminal article by Panagia were adopted (e.g., Crowther 2005). However, the basic result that one is dealing with an early B-type star with a luminosity of 104 L would remain valid. We note that the fit to the continuum spectrum requires a very small size of the UC H II region, of 0. 04 or 300 AU, which appears in contradiction with the fact that our A-array 1.3 cm observations partly resolve the continuum emission (see Table 1). Since we considered a homogeneous distribution of ionized gas in our fit, one possibility is that the H II region is not homogeneous. The southern spur identified in the A-array 3.6 cm map of the "HMC" source (see Figure 4), could indeed suggest the existence of a


­ 24 ­ density gradient toward the south. One may speculate that being source "A" more extended and perhaps more evolved, the gas in its surroundings is less dense. In this scenario, the north-south streamlines of water masers on either sides of the southern spur, might trace shocks at the interface between the ionized gas of "HMC" (expanding to the south) and the surrounding molecular gas. We also note that the continuum spectrum of "HMC" between 1.3 and 6 cm can be also fitted with a power law S , with = 0.6 (see Figure 2). This spectral index is typical of thermal jets (Reynolds 1986) and provides us with an alternative explanation for the origin of the free-free emission. However, a caveat is in order. As discussed in Section 4.1, the slope of the spectrum between 3.6 cm and 1.3 cm could be significantly affected by the inteferometer resolving out part of the emission at the shorter wavelength. If the spectrum is steeper than estimated from our data, the slope could be still consistent with a thermal jet, but also with that of an H II region with a density gradient (see e.g., Franco et al. 2000). We conclude that it is impossible to decide the nature of the free-free emission on the basis of the information collected so far. Notwithstanding this caveat, as a working hypothesis we prefer the jet interpretation, because this permits to interpret the continuum and maser emission in the same context, as discussed in the following section.

5.2.

The Water Maser Flow

In this section we propose that both the continuum and the water maser emission are associated with a thermal jet impinging against the surrounding quiescent molecular gas. This idea is based on two observational findings. First of all, Figure 5 shows that, at larger distances from the center of the water maser distribution, maser (relative) velocities appear


­ 25 ­ to be larger and better collimated along the northeast-southwest direction. Secondly, the axis of the powerful, molecular outflow mapped at arcsecond scales toward G28.87+0.07 has a northeast-southwest orientation (Furuya et al. 2008, Figure 5). This it is plausible to assume that the two, fast moving, water maser clusters "Jb " and "Jr ", located, respectively, at the northeast and southwest edge of the maser distribution, are excited by the same jet responsible also for driving the large-scale molecular outflow. One could ob ject against the "jet" interpretation, that the whole pattern of water maser velocities cannot be explained only in terms of a collimated outflow. In fact, at smaller distances from the center of the maser distribution, water maser velocities point to the north or south directions, and form large angles with the direction of the putative jet. We propose that the scattered pattern of proper motions presented by the red-shifted water masers can result from the interaction of the jet with dense circumstellar material, which deflects and increase the opening angle of the flow velocities. If the free-free continuum source "HMC" is a jet, the position of the (proto)star powering it may be offset from the free-free continuum peak. If a jet powered by the (proto)star drives the water maser motions, and, in particular, the expansion to the south and west of the narrow strip of red-shifted water masers, the (proto)star could be located at an intermediate position between "HMC" and the red-shifted maser strip. This point is not far (within tens of mas) from the position of the center of motion of the water masers (cross in Figure 5). Considering the large uncertainties in deriving the center of water maser expansion, in the following discussion we can safely take the center of motion as a reasonable estimate of the (proto)star location. In the jet scenario, one can estimate the momentum rate of the jet from the average distance of the H2 O masers from the MYSO and the average maser velocity. Assuming a H2 pre-shock density of 108 cm-3 , as predicted by models (Elitzur et


­ 26 ­ al. 1989), the momentum rate in the water maser jet is given by the expression
2 P = 1.5 â 10-3 V10 R 2 100

(/4 ) M yr-1 km s-1 . Here V10 is the average maser velocity in the average distance of water masers from the YSO in units of

units of 10 km s-1 , R

100

100 AU, and the solid angle of the jet. This expression has been obtained by multiplying the momentum rate per unit surface transferred to the ambient gas (n
H2

m

H2

V 2 ), by R2 ,

under the assumption that the jet is emitted from a source at a distance R from the masers within a beaming angle . To estimate the average distance and velocity of the water masers tracing the jet in G28.87+0.07, we use all the maser features belonging to the two fast moving clusters "Jb " and "Jr " (see Figure 5). The resulting average distance is 122 mas (or 902 AU) and the average speed is 43 km s-1 . Using these values, the jet momentum rate is P = 2.2 (/4 ) M yr-1 km s-1 . For shock induced ionization in a thermal jet, the momentum rate of the jet can be estimated from the continuum flux assuming the emission is optically thin (see Anglada et al. (1996)). One has: F d2 = 103.5 (/4 )P , where F is the measured continuum flux density in mJy, P the jet momentum rate in M yr-1 km s-1 , the jet solid angle in sterad, and d the source distance in kpc. Using the flux of 1.2 mJy measured at 1.3 cm with the VLA C-array for source "HMC", and the distance of 7.4 kpc, one derives P = 2 â 10-2 (/4 )-1 M yr-1 km s-1 . Note that the momentum rates estimated from the maser motion and the continuum flux have a different dependency on the jet solid angle, which can thus be estimated by equating the two momentum rate estimates. We obtain an expression for the jet solid angle: 2.2 (/4 ) = 2 â 10-2 (/4 )-1 . We find = 1.2 sterad, i.e. a jet semi-opening angle j = 35 . Correspondingly, the jet momentum rate is P = 0.2 M yr-1 km s-1 . This value is large enough to justify the momentum rate of several 10-2 M yr-1 km s-1 , estimated by Furuya et al. (2008) for


­ 27 ­ the large-scale molecular outflow. This result supports our hypothesis that the small-scale thermal and H2 O maser jet is sufficiently powerful to feed the large-scale molecular outflow observed by Furuya et al. (2008). We note that the derived parameters for the jet in G28.87+0.07 are comparable to those of the few collimated outflows observed up to now in massive star-forming regions. For comparison, the jets in HH 80­81 (Marti et al. 1998), in IRAS 16547-4247 (Rodr´ iguez et al. 2008), and in IRAS 20126+4104 (Moscadelli et al. 2011) have momentum rates of the order of 10-1 M yr-1 km s 20­30 .
-1

and opening angles of

5.3.

Methanol Maser Environment

The two detected 6.7 GHz methanol maser features are approximately aligned along the same northeast-southwest line defined by the two maser jet clusters "Jb " and "Jr ", the peak of the VLA A-array 1.3 cm, and the center of motion (see Figure 4). Figure 3 shows that the total-power spectra of the 6.7 GHz methanol and 22 GHz water masers are similar, both presenting a strong red-shifted and a much weaker blue-shifted component. The positions of the methanol maser features roughly along the jet axis, and the similarities in the spectral shape of the water and methanol maser emissions, make us speculate that also the 6.7 GHz methanol masers could originate in the (proto)stellar jet, albeit associated with more quiescent gas than that traced by the water masers. Maser VLBI observations toward the sources IRAS 20126+4104 and G23.01-0.41 (Moscadelli et al. 2011; Sanna et al. 2010b) have also indicated that the 6.7 GHz masers might trace entrained gas at the interface with a jet. Figure 4 clearly shows that the blue-shifted lobe of the water maser jet is sampled much more sparsely, as it presents much less features than the red-shifted lobe. This is consistent with the blue-shifted line being significantly weaker than the red-shifted one


­ 28 ­ in the total power spectrum (see Figure 3, upper panel). This fact suggests that the circumstellar material is not isotropically distributed around the (proto)star, but it is mostly concentrated on the side of the red-shifted lobe. That would also account for the stronger emission line in the total-power spectrum of the CH3 OH masers at the red-shifted velocities.

5.4.

OH Maser Environment

The spectrum of the 1.665 GHz maser is dominated by three narrow LCP features (see Figure 3, lower panel), and is very similar to the spectrum of the 1.667 GHz maser transition observed by Caswell & Haynes (1983). These authors find that most OH interstellar masers are highly polarized, either RCP or LCP, with the same predominant polarization for both the 1.665 GHz and 1.667 GHz transition. They propose an explanation based on large-scale velocity and/or magnetic field gradients across the masing region. Figure 1 shows that the 1.665 GHz OH maser features are found significantly offset ( 0. 6) from the compact VLA source "HMC" and the center of water maser activity, and are closer to the more extended source "A". A possible interpretation is that OH masers trace more external, less dense parts of the molecular core surrounding the (proto)star responsible for the excitation of the water and methanol masers. Observations indicate that the density of massive proto-stellar cores decreases from the center outward, with a power-law index of about 1.6 (n r
- 1 .6

) (Beuther et al. 2002a). In G28.87+0.07, the ratio

of the average distance (from the "center of motion" of water masers) of the water ( 0. 1) and OH ( 0. 6) masers, is about 6, implying a decrease in gas density by a factor of about 20. That agrees with maser excitation models, which predict a H2 pre-shock density of 106 ­107 cm 108 cm
-3 -3

for OH 1.665 GHz (Cesaroni & Walmsley 1991; Cragg et al. 2002), and

for H2 O 22 GHz masers (Elitzur et al. 1989), respectively.


­ 29 ­ 5.5. The Nature of Source "HMC"

So far, we have explored the properties of source "HMC" making use of the maser, radio, and IR continuum emission. Now, we want to use our findings to shed light on the nature of the MYSO(s) associated with this source. We have seen that the radio continuum emission and the water maser features suggest the presence of a jet with a large momentum rate, on the order of 10-1 M km s
-1

yr-1 .

Values that high are consistent with the jet being powered by a high-mass (proto)stars, as indicated, e.g., by the relationship between outflow momentum rate and MYSO luminosity obtained by Beuther et al. (2002b). In particular, from their Figure 4b, one can conclude that the MYSO must have a luminosity of at least several 103 L . Alternatively, if the radio emission were originating from an H II region, we have seen (Section 5.1) that the ionizing star should have a luminosity 104 L . Clearly, both the thermal jet and the H II

region hypotheses imply similar luminosity for the associated MYSO. The question we want to investigate here is whether such a luminosity is consistent with the results obtained from the IR data. In Section 4.4 we have found that the total luminosity of the region is on the order of 2 â 105 L . However, this includes all MYSOs in a radius 10 around "HMC", whereas

we need to estimate the contribution of the sole "HMC", to achieve a consistent comparison with the luminosity estimate obtained from the water maser and/or radio continuum. Our 24.5 µm continuum image can be used for this purpose. As shown in Figure 7 , the "HMC" source is dominating the emission at this wavelength, contributing to 90% of the total flux density, and a similar conclusion holds also for the GLIMPSE images. Due to the lack of sub-arcsecond imaging at longer wavelengths, it is difficult to establish whether also the flux density measured close to the peak of the SED is mostly due to source "HMC". Nonetheless, the Hi-GAL/Herschel data suggest that the 70 µm emission peaks closer to


­ 30 ­ source "HMC" than source "A". A word of caution is in order here. The angular resolution of the Hi-GAL data at 70 µm is estimated to be about 10 â 9 (after reduction with the Roma-Gal software), while the astrometrical error is 3. 5 (the latter results from aligning the PACS 70 µm images to the MIPSGAL 24 µm images by matching suitable point sources detected at both wavelengths). Although large, these values still permit to associate the 70 µm emission with source "HMC". In fact, the angular separation between the 3.6 cm continuum source "HMC" and the 70 µm peak is 3. 7, consistent with the positional error, whereas source "A" lies 6. 5 from the 70 µm peak. This can be appreciated in Figure 10, where one sees that the 70 µm emission peaks to the NW of "HMC". Therefore, in the following we assume that the bolometric luminosity estimated from the SED is likely dominated by the "HMC" source and conclude that the MYSO(s) embedded in the HMC may be as luminous as 2 â 105 L . Notwithstanding the large uncertainties of all our luminosity estimates, the value obtained from the IR images exceeds by an order of magnitude that derived from the maser and radio continuum data. This inconsistency can be explained if "HMC" is hosting a multiple stellar system rather than a single star. One can get an idea of how this may affect our estimate of the stellar properties (mass, etc.) by assuming that the luminosity of 2 â 105 L is emitted by a stellar cluster with a Miller & Scalo (1979) initial mass function. One finds that the luminosity of the most massive star may be as small as 3 â 104 L , depending on the statistical distribution of the stellar masses in the cluster. Given the uncertainties, this value is only marginally greater than the luminosity estimate from the radio continuum and maser data. The same method provides us with the total mass of the stellar cluster, which ranges from 5 â 102 M to 2 â 103 M , depending on the statistical distribution of the stars in the cluster. This stellar mass can be divided by the total mass of the associated molecular


­ 31 ­ clump, 3 â 103 M (see Section 4.4), to calculate the star formation efficiency, 17%­67%. Albeit large, values like these are not implausible, as we are considering the part of the molecular cloud where star formation concentrates (Lada & Lada 2003; Bonnell et al. 2011). In conclusion, we believe that the observational evidence collected so far for source "HMC" suggests that one is observing a young star of 20­30 M , still in the main accretion phase and surrounded by a rich stellar cluster. Such a conclusion holds independently of the scenario adopted to explain the free-free emission (thermal jet or weak H II region) and raises a question: why only faint free-free emission is detected from this ob ject? A 3 â 104 L star is expected to ionize the surrounding gas emitting a flux density of 2 mJy (in the optically thin limit), whereas towards source "HMC" less than 1 mJy is measured. A possible explanation is that the H II region is at least partly quenched by accretion from the parental core, as depicted, e.g., in the model by Keto (2002, 2003). Although no evidence of accretion is seen in a typical infall tracer such as HCO+ (Furuya et al. 2008), this fact is inconclusive because the ma jority of the detected HCO+ (1­0) might be part of the outflowing gas, thus hindering the detection of infall. Finally, it is also worth mentioning that according to recent theoretical models (Hosokawa et al. 2010), large accretion rates are bound to bloat the protostellar radius and decrease the corresponding temperature and Lyman continuum flux. This could be an alternative explanation of the lack of a (strong) H II region.

6.

SUMMARY

We observed the HMSFR G28.87+0.07 with the VLBI phase-referencing technique in three powerful maser transitions: 22.2 GHz H2 O, 6.7 GHz CH3 OH, and 1.665 GHz OH. In addition, we also performed VLA observations of the radio continuum emission at


­ 32 ­ 1.3 and 3.6 cm with both the A- and C-array, and Subaru/COMICS observations of the mid-infrared continuum emission. We also made use of data from the Hi-GAL/Herschel, ATLASGAL/APEX, as well as other surveys of the Galactic plane. From all these observations, we draw the following conclusions: 1. The bipolar distribution of l.o.s. velocities and the general pattern of observed relative proper motions indicate that the water masers are tracing expansion with an average velocity of 24 km s-1 from a point on the sky close to the maser "center of motion". While at larger distance from the "center of motion", maser velocities are higher and better collimated along the northeast-southwest direction, at closer distances water masers show smaller speeds and a larger scatter in the velocity orientations. We interpret these facts as the result of the interaction of a (proto)stellar jet with dense circumstellar material. 2. Two continuum sources (named "HMC" and "A"), separated by 3. 1 along the north-south direction, are observed with the VLA at 3.6 and 1.3 cm. The "HMC" source is compact and spatially associated with the observed H2 O, CH3 OH and OH masers. Whether the free-free emission arising from this source is due to an H II region or a thermal ionized jet cannot be unambiguously established with the present data. However, we believe that the jet hypothesis is more consistent with the direction of the water maser motions, which outline expansion in a bipolar flow. Moreover, the orientation and momentum rate of the putative jet are compatible with those of the molecular outflow detected at arcsecond scales. We hence argue that the jet traced by the water masers is also powering the large-scale outflow. 3. The combined information obtained from sub-arcsecond imaging with Subaru/COMICS at 24.5 µm and lower angular resolution images from the HiGAL/Herschel and ATLASGAL/APEX surveys, permit to estimate with good


­ 33 ­ accuracy the luminosity (2 â 105 L ) and gas mass (3 â 103 M ) of the star forming region of G28.87+0.07, and establish that 90% of the luminosity is coming from the radio source "HMC". We conclude that this source must contain multiple stars, with the most massive being at least as luminous as 3 â 104 L . The lack of an associated Hii region indicates that the 20­30 M star must be still undergoing heavy accretion from the surrounding envelope.

We thank the anonymous referee for careful reading of the text and the constructive criticisms that improved the presentation and analysis of our results. This work was supported by the Chinese NSF through grants NSF 11133008, NSF 11073054, NSF 10733030, and NSF 10703010.


­ 34 ­ REFERENCES Anglada, G., Estalella, R., Pastor, J., Rodriguez, L. F., & Haschick, A. D. 1996, ApJ, 463, 205 Araya, E., Hofner, P., Goss, W. M., et al. 2007, ApJS, 170, 152 Beltr´ M. T., Cesaroni, R., Codella, C., et al. 2006, Nature, 443, 427 an, Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953 Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945 Beuther, H., Schilke, P., Sridharan, T. K., et al. 2002, A&A, 383,892 Bonnell, I. A., Bate, M. R., & Zinnecker, H. 1998, MNRAS, 298, 93 Bonnell, I. A., & Bate, M. R. 2006, MNRAS, 370, 488 Bonnell, I. A., Smith, R. J., Clark, P. C., & Bate, M. R. 2011, MNRAS, 410, 2339 Bronfman, L., Nyman, L.-A., & May, J. 1996, A&AS, 115, 81 Brunthaler, A., Reid, M. J., Menten, K. M., et al. 2011, Astron. Nachr., 332, 461 Carey, S. J., Noriega-Crespo, A., Mizuno, D. R., et al. 2009, PASP, 121, 76 Caswell, J. L., & Haynes, R. F. 1983, Aust. J. Phys., 36, 417 Caswell, J. L., Vaile, R. A., Ellingsen, S. P., Whiteoak, J. B., & Norris, R. P. 1995, MNRAS, 272, 96 Cesaroni, R., & Walmsley, C. M. 1991, A&A, 241, 537 Codella, C., Testi, L., & Cesaroni, R. 1997, A&A, 325, 282


­ 35 ­ Cohen, M., Walker, R. G., Carter, B., et al. 1999, AJ, 117, 1864 Cragg, D. M., Sobolev, A. M., & Godfrey, P. D. 2002, MNRAS, 331, 521 Crowther, D. M. 2005, in Massive Star Birth: a Crossroads of Astrophysics, IAU Symposium 227, ed. R. Cesaroni, M. Felli, E. Churchwell, M. Walmsley (Cambridge University Press), 389 De Buizer, J. M., Radomski, J. T., Telesco, C. M., & Pina, R. K. 2005, ApJS, 156, 179 ~ Drosback, M. M., Aguirre, J., Bally, J., et al. 2008, A&AS, 212, 9601 Egan, M. P., Price, S. D., & Kraemer, K. E. 2003, A&AS, 203, 5708 Elitzur, M., Hollenbach, D. J., & McKee, C. F. 1989, ApJ, 346, 983 Faundez, S., Bronfman, L., Garay, G., et al. 2004, A&A, 426, 97 ´ Forster, J. R., & Caswell, J. L. 1999, A&AS, 137, 43 Franco, J., Kurtz, S., Hofner, P., et al. 2000, ApJ, 542, L143 Furuya, R. S., Cesaroni, R., Takahashi, S., et al. 2008, ApJ, 673, 363 Furuya, R. S., Cesaroni, R., Shinnaga, H. 2011, A&A, 525, 72 Genzel, R., & Downes, D. 1979, A&A, 72, 234 Hosokawa, T., Yorke, H. W., Omukai, K. 2010, ApJ, 721, 478 Kataza, H., Okamoto, Y., Takubo, S., et al. 2000, Proc. SPIE, 4008, 1144 Keto, E. 2002, ApJ, 580, 980 Keto, E. 2003, ApJ, 599, 1196


­ 36 ­ Klaassen, P. D., Wilson, C. D., Keto, E. R., et al. 2011, ApJ, 530, 53 Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 Little, L. T., & Cesarsky D.A. 1982, A&A, 112, 49 L´ ez-Sepulcre, A., Codella, C., Cesaroni, R., Marcelino, N., & Walmsley, C. M. 2009, op A&A, 499, 811 Marti, J., Rodriguez, L. F., & Reipurth, B. 1998, ApJ, 502, 337 McKee, C. F., & Tan, J. C. 2003, ApJ, 585, 850 Menten, K. M., & Reid, M. J. 1995, ApJ, 445, L157 Miller, G. E., & Scalo, J. M. 1979, ApJ, 41, 513 Molinari, S., Swinyard, B., Bally, J., et al. 2010a, PASP, 122, 314 Molinari, S., Swinyard, B., Bally, J., et al. 2010b, A&A, 518, 100 Moscadelli, L., Menten, K. M., Walmsley, C. M., & Reid, M. J. 2003, ApJ, 583, 776 Moscadelli, L., Xu, Y., & Chen, X. 2010, ApJ, 526, 66 Moscadelli, L., Cesaroni, R., Rio ja, M. J., Dodson, R., & Reid, M. J. 2011, A&A, 526, 66 Neugebauer, G., Habing, H. J., van Duinen, R., et al. 1984, ApJ, 278, 1 Panagia, N. 1973, AJ, 78, 929 Plume, R., Jaffe, D. T., & Evans, N. J., II. 1992, ApJS, 78, 505 Purcell, C. R., Hoare, M. G., & Diamond, P. 2008, ASP Conference Series, 387, 389 Reynolds, S. P. 1986, ApJ, 304, 713


­ 37 ­ Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328 Rodr´ iguez, L. F., Moran, J. M., Franco-Hern´ andez, R., et al. 2008, AJ, 135, 2370 Sanna, A., Moscadelli, L., Cesaroni, R., et al. 2010a, A&A, 517, 71 Sanna, A., Moscadelli, L., Cesaroni, R., et al. 2010b, A&A, 517, 78 Schuller, F., Menten, K., Contreras, Y., et al. 2009, A&A, 504, 415 Sch¨ onrich, R., Binney, J. J., & Dehnen, W. 2010, MNRAS, 403, 1829 Sewilo, M., Watson, C., Araya, E., et al. 2004, ApJS, 154, 553 Shirley, Y. L., Evans, N. J. II, Young, K. E., Knez, C., & Jaffe, D. T. 2003, ApJS, 149, 375 Testi, L., Felli, M., Persi, P., & Roth, M. 1994, A&A, 288, 634 Traficante, A., Calzoletti, L., Veneziani, M., et al. 2011, MNRAS, 416, 2932 Whitney, B. A., Wood, K., Bjorkman, J. E., & Wolff, M. J. 2003a, ApJ, 591, 1049 Whitney, B. A., Wood, K., Bjorkman, J. E., & Wolff, M. J. 2003b, ApJ, 598, 1079 Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481

A This manuscript was prepared with the AAS L TEX macros v5.2.


­ 38 ­

Fig. 1.-- Maps (gray scale) of the VLA continuum emission toward G28.87+0.07 at 1.3 and 3.6 cm. (a) 3.6 cm emission with the VLA C-array. (b) 1.3 cm emission with the VLA

C-array (c) 3.6 cm emission with the VLA A-array (d) 1.3 cm emission with the VLA Aarray. Squares and circles report the positions of the 1.665 GHz OH and the 6.7 GHz CH3 OH masers, respectively, detected in this paper (see Section. 4.2.2 and 4.2.3). In each panel, the contour levels range from 30% to 90% of the peak emission at multiples of 10%. The restoring beam is shown in the lower, left or right corner of each panel. A rectangle overlaid on the image indicates the field of view expanded in the adjacent panel. The associated labels indicate the names of the continuum sources identified by us (see Section 4.1 and Table 1).


­ 39 ­

Fig. 2.-- Spectral energy distribution of the VLA component "HMC" toward the HMSFR G28.87+0.07. Dots and error bars report the values and the associated errors (20%) of the fluxes at 1.3 and 3.6 cm measured with the VLA C-array (see Table 1). The triangle indicates the upper limit of the flux at 6 cm obtained with the VLA B-array (Purcell et al. 2008). The linear fit to these measurements, indicated by the dotted line, produces a spectral index = 0.6. The solid line shows the best fit of the measured fluxes with a model of a homogeneous H II region, yielding a Lyman continuum rate of 7.9 â 1045 s H II region radius of 0. 04.
-1

and a


­ 40 ­

Fig. 3.--
(T
sys

Total-power spectra of the H2 O, CH3 OH, and OH masers toward G28.87+0.07. Upp er panel: system-temperature Middle

) weighted average of the 22.2 GHz total-power sp ectra of the 9 VLBA antennas observing on 2006 April 23.
sys

panel: Effelsberg total-power spectrum of the 6.7 GHz methanol maser emission on 2008 March 18. Lower panel: T

-weighted

average of the 1.665 GHz total-power spectra of the 9 VLBA antennas observing on 2007 April 28. Continuous (black) and dashed (red) lines are used to distinguish between right (RCP) and left (LCP) circular polarizations, resp ectively. The dotted line crossing the spectra denotes the systemic velocity (V
sys

) inferred from NH3 measurements.


­ 41 ­

Fig. 4.--

Spatial distribution of 22 GHz H2 O (triangles), 6.7 GHz CH3 OH (dots), and 1.665 GHz OH (squares) masers in

G28.87+0.07. Different colors are used to indicate the maser LSR velocities, according to the color scale on the righthand side of the plot. Upp er Panel: global view of the region where the three maser species distribute. The solid contours give the 3.6 cm continuum emission observed with the VLA A-array. Contour levels range from 20% to 90% of the peak emission (0.35 mJy beam
-1

) at multiples of 10%. Dashed rectangles enclose the two north-south streamlines of water masers observed on either

side of the VLA 3.6 cm continuum. A dotted rectangle overlaid on the image indicates the field of view expanded in the right panel. Lower Panel: zo om on the area where most of water masers concentrate. The solid contours give the 1.3 cm continuum emission observed with the VLA A-array. Contour levels range from 40% to 90% of the peak emission (0.62 mJy beam multiples of 10%.
-1

) at


­ 42 ­

Fig. 5.--

Proper motions of 22 GHz H2 O masers (triangles) relative to their "center of motion", and spatial distribution of

6.7 GHz CH3 OH masers (dots) in G28.87+0.07. Different colors are used to indicate the maser LSR velocities, according to the color scale on the righthand side of the plot. The cones indicate the 3-D velocities of the water maser features relative to the "center of motion" (marked with a cross). Points without an asso ciated cone have b een detected only over one or two epo chs and the associated proper motion cannot b e computed or is considered unreliable. The cone op ening angle gives the 1 uncertainty on the proper motion direction. The length of the cone is prop ortional to the velocity, with the amplitude scale indicated in the lower left corner of the figure. Open dotted rectangles mark two clusters of water masers located at the northeast and southwest edge of the distribution (labeled "Jb " and "Jr ", resp ectively), moving fast and close to the northeast-southwest direction. The VLA 1.3 cm continuum emission is plotted with solid contours. Contour levels range from 40% to 90% of the p eak emission (0.62 mJy beam
-1

) at multiples of 10%. Upp er Panel: view of the whole distribution of water maser proper motions. Lower

Panel: zoom on the most crowded area of water maser concentration, to the west of the VLA 1.3 cm continuum.


­ 43 ­

Fig. 6.--

Absolute proper motions of water masers in G28.87+0.07, derived by applying different corrections. Upper

Panel: Correcting for the apparent motion due to the combination of the parallax, the solar motion with respect to the LSR, and the differential Galactic rotation (adopting the flat rotation model indicated in Section 3.2). Lower Panel: Correcting for the absolute prop er motion of feature #1 of the CH3 OH masers. Symbols, colors, and contour levels have the same meaning as in Figure 5. The cones indicate the 3-D velocities of the water maser features. Points without an asso ciated cone have been detected only over one or two epo chs and the associated proper motion cannot be computed or is considered unreliable. The cone op ening angle gives the 1 uncertainty on the prop er motion direction. The length of the cone is proportional to the velocity, with the amplitude scale indicated in the lower left corner of the figure.


­ 44 ­

Fig. 7.-- Continuum emission maps toward G28.87+0.07 at centimeter and infrared wavelengths. The grey image in the background is the 24.5 µm emission from Subaru/COMICS. The white and black contours indicate the centimeter continuum emission at 3.6 and 1.3 cm, respectively. The values of the contours are the same as in Figure 1.


­ 45 ­

Fig. 8.-- Spectral energy distribution of the VLA component "HMC" toward G28.87+0.07 from 3.6 µm to 1.1 mm. Dots indicate real values, whereas triangles with the vertex up and down denote lower and upper limits, respectively. Section 4.4 lists the infrared telescopes used to reconstruct the SED. We have fitted the SED with the radiative transfer model developed by Robitaille et al. (2007) (using the SED fitting tool available on http://caravan.astro.wisc.edu/protostars/).


­ 46 ­

Fig. 9.-- Each panel shows the pc-scale structure of the G28.87+0.07 star-forming region at a different infrared wavelength (indicated in the upper left corner of the panel). Images from 70 to 500 µm are observed with Herschel, and the map at 870 µm is taken with APEX. The color scale corresponds to increasing intensity from dark purple to white, as indicated in the wedge to the right. The observing beam is reported in the lower right corner of each panel. The dashed white circle denotes the area of integration to calculate the G28.87+0.07 flux at different infrared wavelengths.


­ 47 ­

Fig. 10.-- Overlay of the 3.6 cm continuum map (contours; same as in Figure 7) on the Hi-GAL/Herschel image at 70 µm. Note how the far-IR emission peaks close to the free-free source "HMC".


­ 48 ­

Table 1. Radio continuum emission in G28.87+0.07.
Peak p osition Label Array (cm) HPBW (â) PA ()


Image rms (mJy b eam
-1

R.A.(J2000) ) (
h m s

Dec.(J2000) (


F

p eak -1

Fint ) (mJy)

)

)

(mJy b eam

HMC

VLA-A VLA-C VLA-A VLA-C

1.3 1.3 3.6 3.6 1.3 1.3 3.6 3.6

0.122 â 0.093 1.053 â 0.755 0.33 â 0.24 2.62 â 2.11 0.122 â 0.093 1.053 â 0.755 0.33 â 0.24 2.62 â 2.11

-17 -9 -1 5 -17 -9 -1 5

0.06 0.05 0.03 0.02 0.06 0.05 0.03 0.02

18 34 46.23 18 43 46.23 18 43 46.23 18 43 46.25 ··· 18 43 46.30 ··· 18 43 46.29

-03 35 29.8 -03 35 29.8 -03 35 29.9 -03 35 29.8 ··· -03 35 32.8 ··· -03 35 32.7

0.62 1.03 0.35 0.53 <0.18(3 ) 0.39 <0.09(3 ) 0.81

0.64 1.20 0.42 0.55 ··· 0.60 ··· 0.98

A

VLA-A VLA-C VLA-A VLA-C

Note. -- Column 1 gives the source lab el; columns 2 and 3 report the VLA array configuration and the observing band, respectively; columns 4 and 5 the (ma jor and minor) HPBW size and the PA of the observing beam, respectively; column 6 the image rms; columns 7, 8, 9, and 10 the absolute position, the peak intensity, and the integral intensity of the emission peak, respectively. The PA of the beam is defined East of North.


­ 49 ­

Table 2. Positions and Brightness.
Maser Feature epoch R.A. (J2000) ( H2 O 1 2006 Apr 23 2006 Jun 30 2006 Sep 28 2007 Jan 18 CH3 OH 1 2006 Feb 28 2007 Mar 18 2008 Mar 18 2009 Mar 17 OH 1 2007 Apr 28
h m s

Dec. (J2000) (


F

p eak -1

NW beam ) (masâmas @ deg) 1.2â0.4 @ -17 1.3â0.4 @ -16 1.4â0.5 @ -3 1.5â1.0 @ 20 9.9â5.6 @ 5 ··· 12.7â5.2 @ 45 12.6â5.9 @ 13 17.8â10.9 @ 1

)

)

(Jy beam 57.4 50.3 32.0 20.7 0.044 ··· 0.246 0.175 1.94

18 43 46.22467 18 43 46.22465 18 43 46.22455 18 43 46.22448 18 43 46.22079 ··· 18 43 46.22031 18 43 46.22011 18 43 46.26771

-03 35 29.8161 -03 35 29.8173 -03 35 29.8190 -03 35 29.8212 -03 35 29.8379 ··· -03 35 29.8521 -03 35 29.8595 -03 35 30.4697

Note. -- Column 1 reports the maser species; column 2 gives the label number of the feature, as given in Tables 3, 4, and 5; column 3 lists the observing date; columns 4 and 5 report the absolute p osition; column 6 gives the p eak brightness; column 7 reports the (ma jor and minor) HPBW size and PA of the naturally-weighted (NW) b eam. The PA of the b eam is defined East of North.


­ 50 ­ Table 3. Parameters of VLBA 22.2 GHz H2 O maser features in G28.87+0.07.
Feature Number epochs of detection V S ) (mas) (mas) V (km s
-1

LSR -1



V ) (km s
-1

(km s

(Jy)

)

0 1 2 3 4 5 ···

1,2,3,4 1,2,3,4 1,2,3,4 1 1,2,3 2,3,4

103.17 106.01 105.07 104.36 105.27 104.27

··· 53.99 26.07 25.38 19.92 17.62

66.83±0.05 0 37.44±0.05 -2.88±0.05 -0.85±0.05 -3.29±0.05

0.005±0.05 0 -3.02±0.05 2.60±0.05 0.25±0.05 2.79±0.05

0 -16.9±1.0 -23.0±1.2 ··· -2.8±2.2 -29.8±1.7

0 1.4±1.0 -26.4±1.2 ··· 0.2±2.3 0.2±1.8

Note. -- Column 1 gives the feature lab el number; column 2 lists the observing epo chs at which the feature was detected; columns 3 and 4 report the value of the intensity-weighted LSR velocity and flux density of the strongest sp ot, averaged over the observing ep ochs; columns 5 and 6 give the p osition offsets (with the associated errors) along the R.A. and Dec. axes relative to the feature #1, measured at the first epoch of detection; columns 7 and 8 give the comp onents of the relative proper motion (with the associated errors) along the R.A. and Dec. axes, measured with respect to the reference feature #0. This table is published in its entirety in the electronic edition of the Astrophysical Journal. A p ortion is shown here for guidance regarding its form and content.

Table 4. Parameters of EVN 6.7 GHz CH3 OH maser features in G28.87+0.07.
Feature epochs of detection 1 2 1,3,4 3,4 V F ) I
-1

LSR -1

p eak

(mas) 0 -12.41±0.18

(mas) 0 -11.24±0.18

v (km s
-1

v ) (km s
-1

(km s

(Jy beam

)

(Jy) 0.190 0.283

)

104.88 105.21

0.155 0.268

1.7±6.7 ···

-15.6±8.1 ···

Note. -- Column 1 gives the feature label number; column 2 lists the observing ep ochs at which the feature was detected; columns 3 and 4 provide the intensity-weighted LSR velocity and the flux density of the strongest spot, respectively, averaged over the observing epochs; columns 5 and 6 give the p osition offsets (with the asso ciated errors) along the R.A. and Dec. axes, relative to the feature #1; columns 7 and 8 give the components of the relative proper motion (with the asso ciated errors) along the R.A. and Dec. axes, measured with respect to the reference feature #0 (the "center of motion") of water masers.


­ 51 ­

Table 5. Parameters of VLBA 1.665 GHz OH maser features in G28.87+0.07.

Feature

Pol

VLSR (km s-1 )

Fp

eak

I (Jy)

x (mas)

y (mas)

(Jy beam-1 )

1 2 3

L L L

102.9 106.1 105.4

1.94 1.44 0.67

5.26 3.52 1.78

0

0

-233.57±0.11 143.92±0.13 -139.46±0.17 46.91±0.21

Note. -- Column 1 gives the feature label number; column 2 indicates the circular polarization of the maser emission; columns 3, 4, and 5 report the intensity-weighted LSR velocity, and the intensity and the integral flux of the strongest spot, respectively; columns 6 and 7 list the position offsets (with the associated errors) along the R.A. and Dec. axes, relative to the feature #1.


­ 52 ­

Table 6. Flux densities measured towards the G28.87+0.07 region.
(µm) F instrument



(Jy)

3.6 4.5 5.8 8.0 8.3 12 12.1 14.7 21.3 24 24.6 25 60 70 100 160 250 350 500 870 1100
a

0.9 2.1 5.3 9.4 3.3 8.6
a

Spitzer/GLIMPSE Spitzer/GLIMPSE Spitzer/GLIMPSE Spitzer/GLIMPSE MSX IRAS MSX MSX MSX Spitzer/MIPSGAL Subaru/COMICS IRAS IRAS Herschel/Hi-GAL IRAS Herschel/Hi-GAL Herschel/Hi-GAL Herschel/Hi-GAL Herschel/Hi-GAL APEX/ATLASGAL CSO/BGPS

7.1 15.3 34.0 >21 87.4 105
a a

1720

1790 <3180 1350 630 200 53 8.8 2.5

all IRAS fluxes are assumed to be upp er limits in Figure 8


­ 53 ­ 7. ONLINE-ONLY MATERIALS


­ 54 ­

Table 7. Parameters of VLBA 22.2 GHz H2 O maser features in G28.87+0.07. (Table 3 in the print version.)

Feature Number

epochs of detection

V

LSR -1

S )



(mas)

(mas)

V (km s
-1

V ) (km s
-1

(km s

(Jy)

)

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29

1,2,3,4 1,2,3,4 1,2,3,4 1 1,2,3 2,3,4 2 2 1,2,3,4 1,2,3,4 1,2,3,4 1 1,2,3,4 1,2,3 2 1,2,3,4 1,2,3,4 2 2,3,4 1,2,3,4 2 1,2,3,4 1,2 1,2,3,4 2 2 1,2,3 1,2,3,4 2 1,2,3,4

103.17 106.01 105.07 104.36 105.27 104.27 105.06 105.14 105.21 105.45 105.25 103.92 104.80 104.00 105.35 106.34 106.73 105.12 104.94 105.87 104.70 105.31 104.23 105.42 105.11 105.10 95.16 105.50 104.36 103.91

··· 53.99 26.07 25.38 19.92 17.62 15.08 14.36 12.16 11.21 9.40 9.39 9.30 9.15 8.57 8.22 7.53 7.33 7.23 7.22 7.02 6.90 6.82 6.57 6.48 6.45 6.20 6.16 5.79 5.57

66.83±0.05 0 37.44±0.05 -2.88±0.05 -0.85±0.05 -3.29±0.05 2.36±0.05 14.86±0.05 36.29±0.05 -1.71±0.05 21.35±0.05 -3.79±0.05 29.94±0.05 -5.85±0.05 -1.00±0.05 56.56±0.05 -1.84±0.05 3.31 ±0.05 29.24±0.05 34.93±0.05 -11.90±0.05 29.01±0.05 36.56±0.05 27.79±0.05 -0.06±0.05 15.54±0.05 202.71±0.05 21.32±0.05 -1.96±0.05 -6.62±0.05

0.005±0.05 0 -3.02±0.05 2.60±0.05 0.25±0.05 2.79±0.05 -3.68±0.05 -146.72±0.05 -3.33±0.05 0.79±0.05 -2.95±0.05 3.17±0.05 -5.50±0.05 4.22±0.05 24.53±0.05 -3.51±0.05 0.91±0.05 -46.40±0.05 -165.12±0.05 -5.68±0.05 192.10±0.05 -5.21±0.06 -3.49±0.05 -4.98±0.06 -42.00±0.05 -104.30±0.05 50.48±0.05 21.04±0.05 2.10±0.05 4.52±0.05

0 -16.9±1.0 -23.0±1.2 ··· -2.8±2.2 -29.8±1.7 ··· ··· 0.5±1.3 -26.3±1.5 0.9±1.5 ··· 8.2±1.5 -21.5±2.2 ··· -21.2±1.4 -25.3±1.1 ··· ··· 21.1±1.6 ··· 7.6±1.5 ··· 8.3±1.6 ··· ··· 55.3±1.8 11.5±1.9 ··· -42.5±1.3

0 1.4±1.0 -26.4±1.2 ··· 0.2±2.3 0.2±1.8 ··· ··· -10.0±1.3 7.0±1.5 -5.7±1.5 ··· -6.9±1.6 -13.5±2.3 ··· 3.0±1.5 6.4±1.2 ··· ··· -25.0±1.6 ··· -0.7±1.6 ··· -12.6±1.7 ··· ··· 18.4±1.9 -11.6±1.9 ··· -4.2±1.3


­ 55 ­

Table 7--Continued
Feature Number epochs of detection V S )


LSR -1

(mas)

(mas)

V (km s
-1

V ) (km s
-1

(km s

(Jy)

)

30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59

2,3,4 1,2,3,4 1,2,4 1 1,2,3,4 3 1,2,3,4 1,2,3 2 3 1,2,3,4 2 1,2 2 1 1,2,3,4 2 1,2,3,4 1,2,3 1,2,3,4 1,2,3,4 1,2 1 2,3,4 1,2,3,4 2,3 4 3,4 2,3,4 3

103.78 105.13 104.22 104.35 106.26 104.71 105.46 105.15 105.54 104.72 108.62 104.93 104.71 105.12 104.23 106.20 105.52 104.67 104.04 103.85 105.44 105.06 105.41 105.55 106.82 106.46 105.35 105.95 106.37 105.39

5.30 5.28 5.23 5.23 5.21 5.17 4.86 4.78 4.72 4.57 4.46 4.37 4.32 4.17 4.06 4.03 3.96 3.95 3.91 3.84 3.71 3.67 3.60 3.24 3.22 3.13 3.07 3.02 2.93 2.90

-2.48±0.05 31.91±0.05 -2.89±0.05 -8.60±0.05 21.33±0.05 14.22±0.05 20.48±0.05 28.10±0.05 36.95±0.05 -10.01±0.05 2.52±0.05 36.76±0.05 -5.05±0.05 25.61±0.05 -6.09±0.05 61.07±0.05 33.94±0.05 -2.36±0.05 -4.62±0.05 -6.13±0.05 21.97±0.05 38.10±0.05 -0.82±0.05 35.42±0.05 23.89±0.05 13.73±0.05 3.11±0.05 -4.34±0.05 4.53±0.05 25.72±0.05

2.28±0.05 -9.86±0.05 -40.09±0.05 8.75±0.05 -4.99±0.05 -239.02±0.05 -4.65±0.05 11.48±0.05 13.71±0.05 121.21±0.05 -2.55±0.05 -27.47±0.05 46.44±0.05 0.01±0.05 30.67±0.06 -7.08±0.05 43.85±0.05 26.39±0.05 3.58±0.05 47.09±0.05 -20.76±0.05 -20.04±0.05 -23.65±0.05 -29.60±0.05 -3.16±0.05 -123.75±0.05 -47.13±0.06 5.38±0.06 -2.16±0.05 -25.59±0.05

4.2±1.9 -5.2±1.8 -50.6±2.7 ··· -1.6±1.6 ··· 4.1±1.7 1.9±3.1 ··· ··· -7.9±1.2 ··· ··· ··· ··· 6.2±1.7 ··· -50.8±1.8 -41.3±2.5 -31.9±1.7 1.2±1.9 ··· ··· 10.5±2.7 -3.5±1.4 ··· ··· ··· -38.8±2.4 ···

-23.6±2.0 -26.4±2.0 -8.5±2.7 ··· -6.3±1.6 ··· -3.7±1.7 -8.7±3.2 ··· ··· -6.2±1.2 ··· ··· ··· ··· -5.8±1.7 ··· 13.5±1.8 7.8±2.7 10.5±1.7 1.5±2.1 ··· ··· -22.8±2.7 -9.8±1.4 ··· ··· ··· -21.0±2.4 ···


­ 56 ­

Table 7--Continued
Feature Number epochs of detection V S )


LSR -1

(mas)

(mas)

V (km s
-1

V ) (km s
-1

(km s

(Jy)

)

60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89

2 1 4 1,2,3,4 3 1,2 4 1,2,3 1,2,3 1,2,3,4 2,3 2 4 3 1,2,3,4 2 3 1,2,3,4 3 3 3 3 2,3 1,2,3,4 3 4 1 3 3 4

105.01 104.90 105.34 105.46 105.83 105.30 103.22 105.21 106.13 106.85 105.04 106.54 105.71 105.30 106.74 106.40 106.37 104.43 106.05 106.54 106.27 105.94 106.59 94.90 106.57 105.12 105.73 106.06 106.37 105.42

2.86 2.81 2.77 2.72 2.63 2.51 2.40 2.34 2.30 2.28 2.26 2.21 2.15 2.08 1.97 1.95 1.91 1.87 1.87 1.82 1.82 1.81 1.80 1.80 1.72 1.67 1.66 1.66 1.64 1.61

0.23±0.05 -1.77±0.05 2.65±0.05 37.26±0.05 35.05±0.05 57.47±0.05 -5.42±0.05 38.23±0.05 20.32±0.05 3.34±0.05 40.16±0.05 26.23±0.05 0.42±0.06 6.51±0.05 22.28±0.05 56.00±0.05 20.16±0.05 23.93±0.05 -4.11±0.05 8.25±0.05 20.56±0.05 -1.45±0.05 3.66±0.05 138.95±0.05 23.19±0.05 38.38±0.06 22.67±0.05 2.01±0.05 20.18±0.05 29.71±0.06

1.02±0.05 2.26±0.05 -5.13±0.06 20.82±0.05 -48.30±0.06 -3.62±0.05 3.83±0.05 -3.01±0.05 -47.78±0.05 -2.71±0.05 -5.89±0.05 -8.07±0.05 -42.97±0.06 -1.41±0.05 -2.91±0.05 -3.39±0.05 -23.71±0.06 -4.51±0.05 71.74±0.06 35.04±0.05 -65.42±0.05 43.02±0.06 -45.51±0.05 67.32±0.05 -46.20±0.05 33.81±0.06 -4.96±0.05 -1.62±0.06 1.15±0.05 -98.11±0.06

··· ··· ··· -9.8±1.8 ··· ··· ··· -14.0±3.5 -10.4±3.1 -3.2±1.4 ··· ··· ··· ··· -5.8±1.7 ··· ··· -9.5±1.8 ··· ··· ··· ··· ··· 6.7±1.2 ··· ··· ··· ··· ··· ···

··· ··· ··· -12.8±1.8 ··· ··· ··· 12.2±3.7 2.4±3.3 -3.8±1.5 ··· ··· ··· ··· 0.8±1.7 ··· ··· -2.7±2.0 ··· ··· ··· ··· ··· 14.6±1.2 ··· ··· ··· ··· ··· ···


­ 57 ­

Table 7--Continued
Feature Number epochs of detection V S )


LSR -1

(mas)

(mas)

V (km s
-1

V ) (km s
-1

(km s

(Jy)

)

90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119

1,2,3,4 1 1,2 4 1,2,3 2,3,4 4 1,2,3,4 3 2,3 4 2 1 2,3 3 4 4 2,3 1,2,3,4 2,3 3 1 1,2,3 4 3 4 4 1,2,3,4 1,2 4

95.86 103.89 104.29 105.12 104.78 107.30 105.31 106.59 103.95 106.78 104.94 106.96 103.88 106.11 104.07 105.11 105.50 111.99 106.79 108.68 106.57 106.80 106.82 106.91 104.67 103.55 103.92 103.66 94.32 103.64

1.60 1.58 1.52 1.49 1.48 1.45 1.39 1.37 1.35 1.30 1.28 1.25 1.25 1.24 1.12 1.11 1.11 1.10 1.09 1.03 0.97 0.94 0.86 0.83 0.81 0.79 0.73 0.72 0.69 0.68

220.77±0.05 0.44±0.05 23.52±0.05 31.87±0.06 7.67±0.05 -3.86±0.05 29.71±0.06 5.76±0.05 -6.91±0.05 46.15±0.05 30.01±0.06 52.39±0.05 -6.46±0.05 -0.24±0.05 -0.25±0.05 33.36±0.06 37.49±0.06 71.41±0.05 49.43±0.05 100.94±0.05 1.57±0.05 1.44±0.05 48.34±0.05 -1.46±0.06 23.33±0.05 -11.43±0.05 -7.58±0.06 179.75±0.05 139.93±0.05 0.34±0.06

55.75±0.06 -2.15±0.05 -4.52±0.05 -34.73±0.06 -1.02±0.05 1.76±0.06 -55.31±0.06 -2.52±0.05 -38.75±0.06 0.89±0.05 -31.50±0.06 -8.39±0.05 7.37±0.07 3.29±0.06 -44.67±0.06 2.56±0.06 -29.30±0.06 0.79±0.06 -3.39±0.05 20.38±0.06 -46.93±0.05 -1.20±0.05 -3.64±0.05 -42.03±0.06 -4.89±0.06 10.12±0.06 8.39±0.06 -103.20±0.05 66.86±0.05 -114.27±0.06

53.6±1.4 ··· ··· ··· -13.0±3.9 -27.6±2.3 ··· -22.6±1.8 ··· ··· ··· ··· ··· ··· ··· ··· ··· ··· 6.3±1.7 ··· ··· ··· 2.5±4.0 ··· ··· ··· ··· 18.9±2.2 ··· ···

17.8±1.5 ··· ··· ··· -12.2±4.1 7.0±2.4 ··· -9.8±1.9 ··· ··· ··· ··· ··· ··· ··· ··· ··· ··· -14.1±1.7 ··· ··· ··· -9.2±4.2 ··· ··· ··· ··· 1.3±2.4 ··· ···


­ 58 ­

Table 7--Continued
Feature Number epochs of detection V S )


LSR -1

(mas)

(mas)

V (km s
-1

V ) (km s
-1

(km s

(Jy)

)

120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149

3,4 1,2,3 1 1,2,3 3,4 2 1 1,2,3,4 4 1,2,3,4 3 2,3,4 2 1,2,3 1,2 1,2,3,4 1,2,3,4 1,2,3 3,4 4 3 3 1,2 1,2,3 1,2,3 4 1,2,3 1,2,3 1,2 4

107.47 107.35 107.18 94.48 107.11 95.24 107.01 94.75 108.82 93.64 103.41 107.11 95.94 93.60 99.32 93.62 95.47 109.59 99.16 93.97 103.39 107.37 95.14 107.23 108.39 95.90 94.16 111.87 99.12 107.41

0.66 0.61 0.57 0.57 0.53 0.49 0.48 0.41 0.38 0.38 0.38 0.37 0.36 0.36 0.34 0.34 0.34 0.33 0.31 0.31 0.30 0.28 0.27 0.25 0.25 0.25 0.23 0.22 0.22 0.21

167.34±0.05 50.20±0.05 -3.75±0.05 214.01±0.05 -6.35±0.05 202.75±0.05 -5.96±0.05 206.85±0.05 -1.59±0.06 140.63±0.05 -6.86±0.06 84.03±0.05 138.45±0.05 213.24±0.05 121.98±0.05 205.00±0.05 204.63±0.06 98.29±0.05 122.29±0.05 214.80±0.06 58.58±0.06 65.53±0.05 202.97±0.05 166.00±0.05 45.92±0.05 218.01±0.06 136.39±0.05 171.96±0.05 122.43±0.05 -8.25±0.06

-261.03±0.06 -3.27±0.05 0.95±0.05 54.68±0.06 6.41±0.07 51.60±0.05 6.03±0.05 51.60±0.06 3.18±0.06 66.57±0.05 74.45±0.07 -0.53±0.06 67.75±0.06 55.03±0.11 91.75±0.06 51.40±0.06 51.38±0.06 11.15±0.06 91.99±0.06 55.26±0.06 -3.38±0.07 5.02±0.06 49.48±0.05 -257.48±0.06 -6.11±0.05 61.98±0.07 74.76±0.06 -42.98±0.05 91.64±0.05 7.57±0.07

··· -0.6±3.7 ··· 48.1±2.6 ··· ··· ··· 43.1±1.5 ··· 5.9±1.6 ··· -0.8±2.7 ··· 46.1±2.9 ··· 42.0±1.7 52.1±1.8 -12.3±3.0 ··· ··· ··· ··· ··· ··· 2.7±3.7 ··· 6.9±1.9 20.3±2.3 ··· ···

··· -9.6±4.0 ··· 18.7±2.8 ··· ··· ··· 14.0±1.6 ··· 14.3±1.6 ··· -15.5±2.9 ··· 21.3±3.1 ··· 11.4±1.7 17.2±2.1 -12.6±3.4 ··· ··· ··· ··· ··· ··· -8.6±4.4 ··· 7.9±2.2 -6.5±2.6 ··· ···


­ 59 ­

Table 7--Continued
Feature Number epochs of detection V S )


LSR -1

(mas)

(mas)

V (km s
-1

V ) (km s
-1

(km s

(Jy)

)

150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168

4 2 4 1,2 1,2,3 2,3,4 1 4 1,2,3 1 2 3 2 4 2 3 3 2 3

103.23 95.20 108.66 102.57 95.22 93.33 99.90 103.02 96.62 95.29 107.64 108.71 107.42 94.87 94.98 108.47 95.84 98.34 107.50

0.21 0.19 0.19 0.16 0.16 0.16 0.14 0.14 0.13 0.13 0.13 0.13 0.12 0.12 0.12 0.12 0.11 0.11 0.10

-6.00±0.07 135.49±0.05 2.04±0.07 100.89±0.05 202.79±0.06 205.03±0.06 99.55±0.05 -2.92±0.07 217.56±0.05 202.24±0.06 182.83±0.05 2.47±0.06 170.05±0.06 135.39±0.06 201.51±0.06 82.90±0.06 225.25±0.06 128.73±0.06 -1.99±0.07

-22.28±0.07 77.14±0.07 -4.78±0.07 8.93±0.05 56.68±0.05 58.09±0.06 9.10±0.06 -0.52±0.08 54.15±0.06 47.80±0.06 -174.07±0.07 -4.01±0.05 -263.08±0.08 73.40±0.08 43.74±0.08 1.27±0.09 7.28±0.10 86.69±0.09 2.06±0.06

··· ··· ··· ··· 46.4±3.0 27.8±4.5 ··· ··· 51.8±2.0 ··· ··· ··· ··· ··· ··· ··· ··· ··· ···

··· ··· ··· ··· 1.8±3.7 -1.8±4.5 ··· ··· 18.0±2.4 ··· ··· ··· ··· ··· ··· ··· ··· ··· ···

Note. -- Table 3 in printer. Column 1 gives the feature label number; column 2 lists the observing epo chs at which the feature was detected; columns 3 and 4 report the value of the intensity-weighted LSR velocity and flux density of the strongest spot, averaged over the observing ep ochs; columns 5 and 6 give the p osition offsets (with the asso ciated errors) along the R.A. and Dec. axes relative to the feature #1, measured at the first epo ch of detection; columns 7 and 8 give the components of the relative prop er motion (with the asso ciated errors) along the R.A. and Dec. axes, measured with respect to the reference feature #0.