# A global view on star formation: the GLOSTAR Galactic plane survey

Gisela Ortiz-León, Karl Menten, Andreas Brunthaler, Timea Csengeri, James Urquhart, Friedrich Wyrowski, Yan Gong, Michael Rugel, Sergio Dzib, Aiyuan Yang, et al.

► To cite this version:

Gisela Ortiz-León, Karl Menten, Andreas Brunthaler, Timea Csengeri, James Urquhart, et al. A global view on star formation: the GLOSTAR Galactic plane survey. Astronomy and Astrophysics - A&A, 2021, 651, pp.A87. 10.1051/0004-6361/202140817. insu-03439597

HAL Id: insu-03439597

https://insu.hal.science/insu-03439597

Submitted on 22 Nov 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A global view on star formation: the GLOSTAR Galactic plane survey

## III. 6.7 GHz methanol maser survey in Cygnus X

Gisela N. Ortiz-León<fnref n="1" />, Karl M. Menten<fnref n="1" />, Andreas Brunthaler<fnref n="1" />, Timea Csengeri<fnref n="2" />, James S. Urquhart<fnref n="3" />,
Friedrich Wyrowski<fnref n="1" />, Yan Gong<fnref n="1" />, Michael R. Riegel<fnref n="1" />, Sergio A. Dzib<fnref n="1" />, Aiyuan Yang<fnref n="1" />, Hans Nguyen<fnref n="1" />,
William D. Cotton<fnref n="4" />, Sac Nicté X. Medina<fnref n="1" />, Rohit Dokara<fnref n="1" />, Carsten König<fnref n="1" />, Henrik Beuther<fnref n="5" />,
Jagadheep D. Pandian<fnref n="6" />, Wolfgang Reich<fnref n="1" />, and Nirupam Roy<fnref n="7" />

Received 15 March 2021 / Accepted 14 May 2021

## ABSTRACT

The Cygnus X complex is covered by the Global View of Star Formation in the Milky Way (GLOSTAR) survey, an unbiased radio-wavelength Galactic plane survey, in 4–8 GHz continuum radiation and several spectral lines. The GLOSTAR survey observed the 6.7 GHz transition of methanol (CH₃OH), an exclusive tracer of high-mass young stellar objects. Using the Very Large Array in both the B and D configurations, we observed an area in Cygnus X of 7°×3° in size and simultaneously covered the methanol line and the continuum, allowing cross-registration. We detected thirteen sources with Class II methanol maser emission and one source with methanol absorption. Two methanol maser sources are newly detected; in addition, we found four new velocity components associated with known masers. Five masers are concentrated in the DR21 ridge and W75N. We determined the characteristics of the detected masers and investigated the association with infrared, (sub)millimeter, and radio continuum emission. All maser sources are associated with (sub)millimeter dust continuum emission, which is consistent with the picture of masers tracing regions in an active stage of star formation. On the other hand, only five masers (38 ± 17%) have radio continuum counterparts seen with GLOSTAR within ~1", testifying to their youth. Comparing the distributions of the bolometric luminosity and the luminosity-to-mass ratio of cores that host 6.7 GHz methanol masers with those of the full core population, we identified lower limits $L_{\text{Bol}} \sim 200 L_{\odot}$ and $L_{\text{Bol}}/M_{\text{core}} \sim 1 L_{\odot} M_{\odot}^{-1}$ for a dust source to host maser emission.

**Key words.** masers – ISM: molecules – techniques: interferometric – radio lines: ISM – radio continuum: ISM – stars: formation

## 1. Introduction

The Cygnus X complex is a nearby (1.4 kpc<fnref n="1" />; Rygl et al. 2012), ~10° wide region rich in molecular clouds that is undergoing active star formation and hosts a number of OB associations that testify to the prodigious star formation that has occurred over the past few million years. The vast richness of the present star formation activity in this region is revealed by a multitude of observations of centimeter-wavelength free-free radio emission, (sub)millimeter molecular line emission, dust continuum emission, and infrared-wavelength imaging (see Reipurth & Schneider 2008 for a review). Numerous (~150) massive dense cores and young stellar objects have been identified (Motte et al. 2007; Beerer et al. 2010; Bontemps et al. 2010; Cao et al. 2019),

as have around ten ultracompact (UC) HII regions (Roy et al. 2011; Cao et al. 2019) and hundreds of more developed compact and extended HII regions (Wendker et al. 1991). Products of past star formation manifest themselves as OB associations, with the 3–5 Myr old Cygnus OB2 association being the most prominent (Wright et al. 2010), as well as dozens of evolved objects, including Wolf-Rayet and carbon stars, planetary nebulae, supernova remnants, and the famous X-ray binary Cygnus X 1 (Uyaniker et al. 2001; Kraemer et al. 2010). Figure 1 of Wendker (1984) gives a schematic overview of Cygnus X and its rich composition.

Due to its proximity (being one of the nearest massive star-forming regions) and large gas reservoir, Cygnus X serves as an excellent laboratory for studying the formation process of high-mass stars. The region extends roughly between 76° and 83° in Galactic longitude and -2° and 2° in latitude (see Fig. 1 in Schneider et al. 2006). With a total mass of $3-4 \times 10^6 M_{\odot}$ (Schneider et al. 2006), this complex is one of the most massive giant molecular clouds (GMCs) in the extended solar neighborhood.

<footnote n="1">¹ Max Planck Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
e-mail: gortiz@mpifr-bonn.mpg.de</footnote>

<footnote n="2">² Laboratoire d'astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France</footnote>

<footnote n="3">³ Centre for Astrophysics and Planetary Science, University of Kent, Canterbury, CT2 7NH, UK</footnote>

<footnote n="4">⁴ National Radio Astronomy Observatory (NRAO), 520 Edgemont Road, Charlottesville, VA 22903, USA</footnote>

<footnote n="5">⁵ Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany</footnote>

<footnote n="6">⁶ Department of Earth Space Sciences, Indian Institute of Space Science and Technology, Trivandrum 695547, India</footnote>

<footnote n="7">⁷ Department of Physics, Indian Institute of Science, Bengaluru 560012, India</footnote>

<footnote n="1">¹ This distance comes from the average of the individual distances to several star-forming regions in Cygnus X North. AFGL 2591, which is projected in the direction of Cygnus X South, was found to be more distant, at 3.3 kpc. Whether Cygnus X is a single connected region is still under debate.</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Locations of methanol sources detected with the VLA (blue circles), overlaid on a *Herschel* SPIRE 500 µm continuum image of the Cygnus X region (grayscale). The mosaic was generated by Cao et al. (2019) from data obtained as part of the HOBYS (Motte et al. 2010; Hennemann et al. 2014) and Hi-GAL projects (Molinari et al. 2010). The region outlined in green represents the area mapped with GLOSTAR. A few well-known star-forming regions and radio continuum sources (Downes & Rinehart 1966) are labeled in yellow.</figcaption></figure>

The $S_1-6_0 A^+$ line of methanol at 6.7 GHz was first detected in the interstellar medium by Menten (1991b). The observed complex line shapes, high intensities, and association with star-forming regions immediately indicated maser emission – in fact, after the 22.2 GHz transition from $H_2O$, it is the second strongest and second most frequently found interstellar maser line (e.g., Caswell et al. 1995).

The comprehensive Methanol Multibeam (MMB) survey in this transition that covered the Galactic plane from $l=186^\circ$ to $l=60^\circ$ found 972 maser sources (see Green et al. 2009, 2010, 2012, 2017; Caswell et al. 2010, 2011). Various studies have shown that 6.7 GHz methanol masers are exclusively associated with high-mass star formation (Minier et al. 2003; Bourke et al. 2005; Ellingsen 2006; Pandian et al. 2008, 2010; Xu et al. 2008). Most of these masers were found to have no radio continuum counterpart (at the few mJy sensitivity level of previous observational efforts) by Walsh et al. (1998), who first suggested that most methanol masers are associated with massive young (proto)stellar objects (MYSOs) in their earliest evolutionary phases (i.e., before they excite a UC HII region), which has been widely confirmed. The comprehensive studies of Urquhart et al. (2013, 2015) and Billington et al. (2019) found that 99% of the 958 methanol masers identified by the MMB in the Galactic plane are also associated with submillimeter continuum radiation from dust detected by the 870 µm APEX Telescope Survey of the Galaxy (ATLASGAL; Schuller et al. 2009; Csengeri et al. 2014; Urquhart et al. 2014) and the 850 µm JCMT Plane Survey (Moore et al. 2015; Eden et al. 2017) conducted with the 15 meter *James Clerk Maxwell Telescope (JCMT)*. Given that 99% of sources have been associated with dust, it is probably safe to assume that all methanol masers are thus associated with compact dust clumps that are massive enough to harbor at least one high-mass (proto)star, along with its associated lower-mass (proto)stellar cluster members. We note further that the 6.7 GHz line is a Class II methanol maser transition and thus, by definition, is radiatively pumped (Batrla et al. 1987; Menten 1991a,b). These masers thus naturally require a source emitting strongly in the mid-infrared regime in their close vicinity

(Sobolev & Deguchi 1994; Cragg et al. 2005), which gives a natural explanation for the above-discussed maser–MYSO association.

The Global View of Star Formation in the Milky Way (GLOSTAR) project is a large radio survey aimed at mapping a significant part of the Galactic mid-plane to characterize the properties of sites with star formation in the earliest evolutionary stages. It takes advantage of the upgraded wideband capabilities of the *Karl G. Jansky Very Large Array (VLA)*, which provides a vastly increased bandwidth and better sensitivity. An overview of the survey is given by Brunthaler et al. (2021), while the first catalog of continuum radio sources detected in our pilot region ($l=28^\circ$ to $36^\circ$ with $b=\pm1^\circ$) is presented in Medina et al. (2019).

Here we report on new VLA observations of 6.7 GHz methanol masers that were carried out in the course of our GLOSTAR coverage of the Cygnus X region, from $l=76^\circ$ to $83^\circ$ and $b=-1^\circ$ to $+2^\circ$ (see Fig. 1 for survey coverage). Primarily, this paper presents the properties of the detected methanol masers and their association with radio continuum emission. Since the maser positions are accurate at the subarcsecond level, we establish their associations with other tracers of star formation to investigate the relationship between methanol masers and massive star formation in the Cygnus X complex. In Sect. 2 we describe the observations and data analysis. Section 3 presents the method for source extraction and the detected methanol masers. In Sect. 4 we discuss the properties of the masers and their association with radio continuum, (sub)millimeter, and infrared emission. We also investigate the properties of the dust cores with which the masers are associated and comment on the possible nature of each maser source. Finally, in Sect. 5 we discuss the relationship between maser emission and dust continuum source properties.

## 2. Observations and data reduction

The observations were carried out as part of the GLOSTAR survey, which is described in detail in Brunthaler et al. (2021). Here we discuss the Cygnus X region portion of the survey.
<table><caption>**Table 1.** Summary of the VLA observations.</caption><thead><tr><th>Observing date<br>D-conf.</th><th>Observing date<br>B-conf.</th><th>Galactic coverage</th><th>Beam size <sup>(a)</sup>; PA<br>(arcsec × arcsec; deg)</th></tr></thead><tbody><tr><td>2014-08-16</td><td>2015-05-02</td><td>76° &lt; l &lt; 77°; -1.0° &lt; b &lt; +0.5°</td><td>14 × 12; -28°</td></tr><tr><td>2014-08-25</td><td>2015-05-08</td><td>76° &lt; l &lt; 77°; +0.5° &lt; b &lt; +2.0°</td><td>17 × 11; -43°</td></tr><tr><td>2014-08-27</td><td>2015-05-07</td><td>77° &lt; l &lt; 78°; -1.0° &lt; b &lt; +0.5°</td><td>15 × 12; -54°</td></tr><tr><td>2014-08-24</td><td>2015-05-06</td><td>77° &lt; l &lt; 78°; +0.5° &lt; b &lt; +2.0°</td><td>18 × 11; -45°</td></tr><tr><td>2014-08-22</td><td>2015-04-25</td><td>78° &lt; l &lt; 79°; -1.0° &lt; b &lt; +0.5°</td><td>14 × 11; -46°</td></tr><tr><td>2014-08-07</td><td>2015-04-04</td><td>78° &lt; l &lt; 79°; +0.5° &lt; b &lt; +2.0°</td><td>16 × 12; -27°</td></tr><tr><td>2014-08-16</td><td>2015-04-24</td><td>79° &lt; l &lt; 80°; -1.0° &lt; b &lt; +0.5°</td><td>15 × 12; -46°</td></tr><tr><td>2014-08-15</td><td>2015-04-26</td><td>79° &lt; l &lt; 80°; +0.5° &lt; b &lt; +2.0°</td><td>15 × 12; -42°</td></tr><tr><td>2014-08-11</td><td>2015-05-14</td><td>80° &lt; l &lt; 81°; -1.0° &lt; b &lt; +0.5°</td><td>15 × 11; -29°</td></tr><tr><td>2014-08-17</td><td>2015-05-03</td><td>80° &lt; l &lt; 81°; +0.5° &lt; b &lt; +2.0°</td><td>13 × 11; -27°</td></tr><tr><td>2014-08-06</td><td>2015-05-08</td><td>81° &lt; l &lt; 82°; -1.0° &lt; b &lt; +0.5°</td><td>16 × 12; -29°</td></tr><tr><td>2014-08-18</td><td>2015-05-10</td><td>81° &lt; l &lt; 82°; +0.5° &lt; b &lt; +2.0°</td><td>14 × 11; -26°</td></tr><tr><td>2014-08-26</td><td>2015-05-09</td><td>82° &lt; l &lt; 83°; -1.0° &lt; b &lt; +0.5°</td><td>14 × 11; -41°</td></tr><tr><td>2014-08-12</td><td>2015-05-09</td><td>82° &lt; l &lt; 83°; +0.5° &lt; b &lt; +2.0°</td><td>14 × 11; -38°</td></tr></tbody></table>

**Notes.** <sup>(a)</sup>These values are for the methanol dirty cubes constructed for the D-array configuration data.

To characterize the radio emission and various tracers of star formation, we carried out simultaneous observations of C-band continuum emission (covering 4.2–5.2 and 6.5–7.5 GHz), the methanol (CH$_3$OH) line at 6.7 GHz, seven radio recombination lines (H96α, H98α, H99α, H110α, H112α, H113α, and H114α), and formaldehyde (H$_2$CO) transition at 4.8 GHz. Observations were taken in both the D and B configurations to simultaneously achieve the best sensitivity for extended emission and excellent angular resolution (see Table 1 for details). The D- and B-configuration data provide angular resolutions of 15" × 11" and 1" × 5 × 1"0, respectively. This corresponds to approximate physical scales of 0.1 pc and 1700 au, respectively, at a distance of 1.4 kpc. Two 1 GHz wide base bands in full polarization mode were registered for the continuum. These were centered at 4.7 and 6.9 GHz and consisted of eight 128 MHz wide spectral windows each. The methanol maser line, which has a rest frequency of 6668.5192 MHz (Breckenridge & Kukolich 1995), was observed with 8 MHz of bandwidth and 2048 channels, resulting in a channel spacing of 0.176 km s$^{-1}$, and a total velocity coverage of 360 km s$^{-1}$.

Our observations, which covered almost the whole Cygnus X complex (i.e., an area of 7° × 3°), required a total of 28 epochs, which were scheduled as part of GLOSTAR under program ID 14A-420. Each epoch observed a 1° × 1.5° strip (see Table 1), typically using 532 single pointings. Two 11-second scans were spent on each pointing for a total (on-source) integration time of ~15 s. The spectral line data from each epoch were calibrated independently within the Common Astronomy Software Applications (CASA) package using version v4.6.0 and a customized version of the VLA pipeline$^2$. The continuum subtraction was performed on the visibility data using the CASA task uvcontsub, where we excluded channels expected to contain emission lines in order to fit the continuum. Specifically, we excluded the local standard of rest (LSR) velocity range from approximately -30 to +20 km s$^{-1}$ for the methanol line, which covers the velocity range of the complex (Schneider et al. 2006).

The calibrated visibilities were imaged using the CLEAN algorithm as implemented in the package tclean within CASA

version v5.4 with a pixel size of 2"5 and 0"2 for D- and B-configuration data, respectively, and a spectral resolution of 0.18 km s$^{-1}$. Each 1° × 1.5° stripe requires 2800 × 2800 pixels in D-configuration maps, while B-configuration maps require about 100 times more pixels. Due to limitations of the computer resources required to perform the deconvolution of such a large map over a large number of spectral channels, we only constructed dirty D-configuration images for the methanol line and performed the maser detection search on these dirty cubes (this is described in more detail in Sect. 3). Dirty cubes of the methanol line were produced for the velocity range from -200 to +100 km s$^{-1}$. Maps of the B-configuration data were only constructed, in a subsequent step, for smaller regions (3'4 × 3'4) centered at the positions of the detected masers determined in the D-configuration cubes and using the data from all neighboring pointings. The synthesized beams are on average 15" × 11" at position angle (PA) = -49° and 1" × 5 × 1"0 at PA = -33° in the D and B configuration, respectively. The 1σ rms noise measured in channels free of maser emission is, on average, 0.028 Jy beam$^{-1}$ per 0.18 km s$^{-1}$ channel for both configurations.

To complement our study, we also analyzed GLOSTAR continuum images toward sites with maser emission. A full description of the GLOSTAR continuum data calibration and imaging is given in Brunthaler et al. (2021), while the full analysis of continuum images of Cygnus X will be presented in a forthcoming paper. Here, we briefly discuss the imaging strategy. The calibration and imaging of the continuum data was performed with the Obit package (Cotton 2008). The 2 GHz bandwidth was first rearranged into nine frequency subbands, which were used to image each pointing individually. Then, for each frequency subband the pointings were combined into large individual mosaics to cover the entire observed area. Finally, we combined the different frequencies to obtain the image at the reference frequency, which has circular beams of 19" and 1.5" in the D and B configuration, respectively. Continuum and methanol line maps from Effelsberg observations have also been obtained as part of the GLOSTAR survey (Brunthaler et al. 2021, Rugel et al., in prep.) We note that continuum images were constructed for Effelsberg data, the VLA D configuration, the VLA B configuration, a combination of the VLA D and B (D+B) configurations, and a combination of the VLA D configuration and

<footnote>$^2$ https://science.nrao.edu/facilities/vla/data-processing/pipeline</footnote>
Effelsberg observations. The central frequency of these images is 5.8 GHz. Here, we only use B-configuration continuum maps to study the region of the investigated methanol maser positions and D+B maps of the region around DR21 (see Sect. 4.5). Methanol line data from Effelsberg were also inspected to look for flux variations in the VLA-detected masers (Sect. 4.4). The noise in the continuum images is not uniform, but rather varies across the mapped region, and can be high around strong sources with complex or extended emission. We locally measured the noise in regions close to the maser locations, resulting in $1\sigma$ values in the range from 0.056 to 0.43 mJy beam$^{-1}$ for B-configuration images. For the D configuration, the $1\sigma$ rms noise ranges from 0.10 to 2.6 mJy beam$^{-1}$. The higher values measured in D-configuration data are due to bright extended emission, which is present across the Cygnus X region, and are resolved out by the array in the B configuration. The highest local rms noise occurs around the strong radio source, DR21, a compact HII region.

### 3. Source detection

We searched for methanol line sources in the D-configuration dirty cubes using the source extraction code introduced in Brunthaler et al. (2021). Briefly, this code scans through the full data cubes to locate brightness peaks with a signal-to-noise ratio (S/N) above a fixed value; here we chose a threshold of S/N=4, which was found to be a good balance between detecting weak sources and dismissing artifacts. For the calculation of the S/N, the rms noise was measured in a small spatial box (typically with a size of 50 × 50 pixel). Starting at the bottom-left corner of the image, a pixel within this box is considered to represent a potential source if the S/N of the emission remains above the S/N threshold in at least three consecutive spectral channels. After scanning the bandwidth for this spatial box, the search was repeated in the next adjacent spatial box until the full image (~1600 × 2300 pixel) was scanned. A few additional filters were applied to the candidate detection list in order to remove duplicate detections (both in position and velocity) and spurious detections caused by, for example, sidelobes. That is, we checked if the candidate detection was within 200 pixels of any previous detection. We kept the detection as a potential one when it appeared in the same channel range as the previous detection and had a higher S/N (in this case, we discarded the previous detection) or when it appeared in a different channel range regardless of its S/N. A full description of the source finding algorithm will be given in an upcoming paper (Nguyen et al., in prep.).

Once a list of potential detections was produced by the extraction code, we imaged a region of 512 × 512 pixel (1024 × 1024 pixel) or ~21′ × 21′ (3.4 × 3.4) for the D configuration (B configuration) centered on the recorded brightness peak position. Additionally, we used the visibility data from all VLA pointings (observed at any date) that lay within the 21′ × 21′ region. The cubes were CLEANed to the expected rms by setting a threshold of 38 mJy. This threshold corresponds to the expected $1\sigma$ rms for a single field. The final sensitivities, however, are higher (cf. Sect. 2) because each field was covered by six neighboring fields. D- and B-configuration data were imaged separately using a velocity spacing of 0.18 km s$^{-1}$. The spectral extent of these cubes covers LSR velocities from -100 to +50 km s$^{-1}$. The CLEANed cubes produced for all potential source candidates were first scanned with our source extraction code and then inspected visually to search for detections. Finally, we manually collected a list of channels that showed spectral peaks with a signal above $4\sigma$. From this list we obtained

the source positions and flux densities at individual channels by performing two-dimensional Gaussian fits to the brightness distribution with the CASA task *imfit*. For the Gaussian fits we selected a box, typically with spatial extents of 14 × 14 pixel and 20 × 20 pixel in the D and B configuration, respectively (≈1.3–1.7 times the beam), centered on the position of the brightest pixel.

The statistical error in maser position is given by the astrometric uncertainty, $\theta_{res}/(2 \times S/N)$, where $\theta_{res}$ is the full width at half maximum (FWHM) size of the restoring beam and S/N is that of the source (Thompson et al. 2017). The D- and B-configuration maps of methanol have an average beam size of 13′′ and 1′′2, respectively. For a maser detected at S/N = 10, the formal (statistical) precision in position is ≈0.′7 and ≈0.′06 for the D and B configurations, respectively. However, additional errors – due to, for instance, uncompensated atmospheric delays that arise when applying the phase corrections determined for the calibrator sources to the Cygnus X data – likely worsen the position accuracy. In order to estimate realistic position uncertainties, we compared our VLA positions with the positions determined from observations with the European Very-long-baseline interferometry (VLBI) Network (EVN; Rygl et al. 2012; see Appendix A for details). For the B configuration, the differences have rms values of 0.′15 and 0.′12 and mean values of 0.′09 and 0.′05 in RA and Dec, respectively. Thus, additional position errors of ±0.′18 and ±0.′13, respectively, need to be added to the formal errors. Similarly, from the differences between the VLA D configuration and EVN positions, we obtained systematic errors of ±0.′4 and ±0.′3 in RA and Dec, respectively, for the D configuration.

### 4. Results

#### 4.1. Detected methanol masers

We detected methanol maser emission toward 13 different locations and one case of methanol absorption (toward the continuum emission of the compact HII region DR21); two of these maser sources (G76.8437+0.1233 and G78.9884+0.2211) are newly detected. The spatial distribution of these methanol sources is shown in Fig. 1, overlaid on a mosaic image obtained at 500 µm with the Spectral and Photometric Imaging REceiver (SPIRE) aboard the Herschel Space Observatory (Cao et al. 2019) from data taken as part of the Herschel imaging survey of OB Young Stellar objects (HOBYS; Motte et al. 2010; Hennemann et al. 2014) and the Herschel infrared Galactic Plane Survey (Hi-Gal; Molinari et al. 2010). Remarkably, several of the masers (five in total) are associated with the ridge of dense molecular material that extends from DR21 to W75N, which is the region with the highest degree of star formation activity in Cygnus X and home to numerous high-mass stars.

Table 2 lists the names, equatorial coordinates, peak LSR velocity, and peak and integrated fluxes of the main velocity components (or “features”) of the methanol sources detected in the D-configuration maps. Maser properties of the features as measured in the B-configuration cubes are shown in Table 3, while properties and maps of the distribution of all detected maser spots above $4\sigma$ are provided in the appendix (Table A.1 and Fig. A.3, respectively). Here, “spot” refers to emission detected in a single velocity channel and “feature” to emission observed in contiguous velocity channels at nearly the same position. The features consist of several contributing spots.

Class II methanol masers are, in general, time variable, and the 8–9 month time span between our D- and B-configuration
<table><caption>**Table 2.** Properties of methanol sources from D-configuration maps.</caption><thead><tr><th>Maser ID (1)</th><th>Name (2)</th><th>α/Δα (h:m:s)/(arcsec) (3)</th><th>δ/Δδ (°:′′)/(arcsec) (4)</th><th>V<sub>LSR</sub> (km s<sup>-1</sup>) (5)</th><th>S<sub>v, Peak</sub> (Jy beam<sup>-1</sup>) (6)</th><th>S<sub>v, Int.</sub> (Jy) (7)</th><th>Common name (8)</th></tr></thead><tbody><tr><td>1</td><td>G76.0932+0.1580</td><td>20:23:23.683<br>0.46<br>1.52</td><td>+37:35:36.17<br>0.70<br>-1.32</td><td>4.84<br>-6.50<br>6.46</td><td>0.54 ± 0.02<br>0.46 ± 0.02<br>0.41 ± 0.02</td><td>0.55 ± 0.03<br>0.53 ± 0.04<br>0.49 ± 0.04</td><td>—</td></tr><tr><td>2</td><td>G76.8437+0.1233</td><td>20:25:43.758</td><td>+38:11:13.26</td><td>-5.42</td><td>4.61 ± 0.04</td><td>4.70 ± 0.08</td><td>—</td></tr><tr><td>3</td><td>G78.8870+0.7087</td><td>20:29:24.947</td><td>+40:11:19.81</td><td>-7.04</td><td>0.77 ± 0.02</td><td>0.80 ± 0.04</td><td>AFGL2591</td></tr><tr><td>4</td><td>G78.9690+0.5410</td><td>20:30:22.725</td><td>+40:09:23.90</td><td>4.84</td><td>1.16 ± 0.03</td><td>1.22 ± 0.06</td><td>IRAS 20286+3959</td></tr><tr><td>5</td><td>G78.9884+0.2211</td><td>20:31:47.286<br>-0.34<br>-0.72<br>-0.00<br>1.66</td><td>+39:59:00.52<br>-0.64<br>1.28<br>-0.06<br>-2.84</td><td>-68.60<br>-66.26<br>-55.46<br>-54.92<br>-60.32</td><td>1.87 ± 0.03<br>0.29 ± 0.02<br>0.21 ± 0.02<br>0.20 ± 0.02<br>0.14 ± 0.02</td><td>1.97 ± 0.06<br>0.32 ± 0.04<br>0.20 ± 0.03<br>0.18 ± 0.03<br>0.15 ± 0.03</td><td>—</td></tr><tr><td>6</td><td>G79.7358+0.9904</td><td>20:30:50.699<br>0.13<br>0.09</td><td>+41:02:27.43<br>-0.10<br>-0.07</td><td>-5.60<br>-3.98<br>-3.08</td><td>23.21 ± 0.10<br>5.03 ± 0.03<br>4.28 ± 0.04</td><td>23.42 ± 0.17<br>5.04 ± 0.06<br>4.38 ± 0.06</td><td>IRAS 20290+4052</td></tr><tr><td>7</td><td>G80.8617+0.3834</td><td>20:37:00.991<br>0.04<br>-1.09<br>-0.12<br>-1.96<br>0.52</td><td>+41:34:56.18<br>0.29<br>-0.58<br>-0.60<br>-0.05<br>0.21</td><td>-4.16<br>-2.00<br>-11.00<br>-12.26<br>1.06<br>-1.10</td><td>11.79 ± 0.07<br>4.18 ± 0.03<br>0.90 ± 0.02<br>0.49 ± 0.01<br>0.35 ± 0.02<br>0.30 ± 0.02</td><td>11.89 ± 0.12<br>4.21 ± 0.06<br>0.98 ± 0.04<br>0.51 ± 0.03<br>0.36 ± 0.03<br>0.33 ± 0.04</td><td>DR20</td></tr><tr><td>8<sup>(a)</sup></td><td>G81.6790+0.5378</td><td>20:39:01.165</td><td>+42:19:33.34</td><td>-1.10</td><td>-0.70 ± 0.03</td><td>-1.04 ± 0.07</td><td>DR21</td></tr><tr><td>9</td><td>G81.7219+0.5711</td><td>20:39:01.066<br>0.11<br>-0.29<br>-1.88</td><td>+42:22:48.89<br>-0.02<br>-0.23<br>-0.78</td><td>-2.72<br>-3.08<br>-3.80<br>8.98</td><td>5.11 ± 0.04<br>3.79 ± 0.03<br>1.43 ± 0.03<br>0.23 ± 0.03</td><td>5.09 ± 0.06<br>3.73 ± 0.06<br>1.50 ± 0.05<br>0.38 ± 0.07</td><td>DR21(OH)</td></tr><tr><td>10</td><td>G81.7444+0.5910</td><td>20:39:00.369<br>0.19</td><td>+42:24:37.00<br>-0.00</td><td>4.48<br>3.76</td><td>11.72 ± 0.21<br>8.68 ± 0.10</td><td>12.32 ± 0.39<br>9.02 ± 0.18</td><td>W75S-FIR1 /<br>DR21 B<sup>(b)</sup></td></tr><tr><td>11</td><td>G81.7523+0.5908</td><td>20:39:01.996<br>-0.06<br>0.03<br>0.24<br>0.27</td><td>+42:24:59.13<br>-0.07<br>-0.06<br>-0.32<br>-0.17</td><td>-8.66<br>-5.78<br>-6.86<br>-2.90<br>1.11</td><td>14.06 ± 0.07<br>7.09 ± 0.05<br>4.88 ± 0.05<br>1.80 ± 0.05<br>1.11 ± 0.04</td><td>14.37 ± 0.13<br>7.30 ± 0.09<br>5.03 ± 0.10<br>1.92 ± 0.09<br>1.14 ± 0.08</td><td>W75S-FIR2 /<br>DR21 A<sup>(b)</sup></td></tr><tr><td>12</td><td>G81.7655+0.5972</td><td>20:39:02.945</td><td>+42:25:50.93</td><td>-1.28</td><td>3.18 ± 0.17</td><td>3.82 ± 0.33</td><td>CygX-N53</td></tr><tr><td>13</td><td>G81.8713+0.7807</td><td>20:38:36.435<br>-0.16<br>-0.13<br>-0.27<br>-0.26<br>-0.12<br>-0.14</td><td>+42:37:34.70<br>0.39<br>0.46<br>0.34<br>0.38<br>0.45<br>-0.43</td><td>7.18<br>4.66<br>4.12<br>5.74<br>5.20<br>3.40<br>9.34</td><td>336.19 ± 1.65<br>286.99 ± 1.50<br>235.08 ± 1.17<br>120.86 ± 0.58<br>109.74 ± 0.55<br>92.67 ± 0.40<br>37.03 ± 0.16</td><td>334.63 ± 2.89<br>288.37 ± 2.65<br>236.85 ± 2.06<br>121.43 ± 1.02<br>110.23 ± 0.97<br>94.04 ± 0.71<br>37.43 ± 0.29</td><td>W75N(B)</td></tr><tr><td>14</td><td>G82.3079+0.7296</td><td>20:40:16.642</td><td>+42:56:29.12</td><td>10.42</td><td>20.67 ± 0.07</td><td>20.78 ± 0.13</td><td>—</td></tr></tbody></table>

**Notes.** Column 2 gives the GLOSTAR source name, constructed from Galactic coordinates; Cols. 3 and 4 are equatorial coordinates. We give the absolute position of the maser component with maximum intensity. For the other components, we list position offsets relative to the position of the strongest component. The position uncertainties are ±0′′4 and ±0′′3, in RA and Dec, respectively (see Sect. 3). For each feature, Col. 5 gives the LSR radial velocity of the peak. Columns 6 and 7 are the peak and integrated fluxes at the channel corresponding to the peak of the feature. Column 8 gives the common source names used in the literature. For masers 6, 7, 9, 10, 11, and 13, Rygl et al. (2012) reported astrometric VLBI observations. <sup>(a)</sup>Seen in absorption. <sup>(b)</sup>Alternative nomenclature used by Rygl et al. (2012).

observations (see Table 1) allows us to check the variability in our sources. Individual maser components can spread by up to several arcseconds on the sky (see Fig. A.3 and also Hu et al. 2016), which is comparable to our B-configuration synthesized beam; as such, to allow for a meaningful comparison, we convolved the B-configuration cubes to the D-configuration beam and then, for every source, extracted each spectrum at the same pixel.

In Fig. 2, the D- and B-configuration spectra are compared with one another. We note that the noise in the B-configuration images convolved with the D-configuration beam is higher by a factor of ≈10. However, this is expected since the noise, in brightness temperature units, decreases with the square root of the beam area, and the brightness temperature is converted to flux density by multiplying with the beam area. This results in a linear increase in the noise with the beam size. Because of this,
<table><caption>**Table 3.** Properties of methanol sources from B-configuration maps.</caption><thead><tr><th>Maser ID (1)</th><th>Name (2)</th><th>α/Δα (h:m:s)/(arcsec) (3)</th><th>δ/Δδ (°:′:′′)/(arcsec) (4)</th><th>V<sub>LSR</sub> (km s<sup>-1</sup>) (5)</th><th>S<sub>v, Peak</sub> (Jy beam<sup>-1</sup>) (6)</th><th>S<sub>v, Int.</sub> (Jy) (7)</th><th>Common name (8)</th></tr></thead><tbody><tr><td>1</td><td>G76.0932+0.1580</td><td>20:23:23.7140<br>-0.103<br>-0.281</td><td>+37:35:35.558<br>0.055<br>-0.005</td><td>6.64<br>4.84<br>-6.50</td><td>0.57 ± 0.01<br>0.38 ± 0.01<br>0.17 ± 0.01</td><td>0.59 ± 0.02<br>0.48 ± 0.01<br>0.18 ± 0.01</td><td>—</td></tr><tr><td>2</td><td>G76.8437+0.1233</td><td>20:25:43.7657</td><td>+38:11:12.734</td><td>-5.42</td><td>0.17 ± 0.00</td><td>0.17 ± 0.01</td><td>—</td></tr><tr><td>3</td><td>G78.8870+0.7087</td><td>20:29:24.9436</td><td>+40:11:19.626</td><td>-7.04</td><td>0.10 ± 0.01</td><td>0.19 ± 0.03</td><td>AFGL2591</td></tr><tr><td>4</td><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>4.84</td><td>1.25 ± 0.01</td><td>1.28 ± 0.02</td><td>IRAS 20286+3959</td></tr><tr><td>5</td><td>G78.9884+0.2211</td><td>20:31:47.3082<br>-0.028</td><td>+39:58:59.903<br>-0.103</td><td>-68.96<br>-66.26</td><td>0.25 ± 0.01<br>0.21 ± 0.01</td><td>0.30 ± 0.01<br>0.30 ± 0.02</td><td>—</td></tr><tr><td>6</td><td>G79.7358+0.9904</td><td>20:30:50.6680<br>0.055<br>0.045</td><td>+41:02:27.403<br>-0.010<br>-0.039</td><td>-5.42<br>-3.98<br>-3.26</td><td>22.96 ± 0.21<br>5.88 ± 0.09<br>4.97 ± 0.06</td><td>23.81 ± 0.39<br>6.11 ± 0.15<br>5.07 ± 0.12</td><td>IRAS 20290+4052</td></tr><tr><td>7</td><td>G80.8617+0.3834</td><td>20:37:00.9571<br>-0.029</td><td>+41:34:55.606<br>0.331</td><td>-3.98<br>-2.00</td><td>4.61 ± 0.11<br>1.56 ± 0.06</td><td>5.02 ± 0.21<br>1.77 ± 0.12</td><td>DR20</td></tr><tr><td>9</td><td>G81.7219+0.5711</td><td>20:39:01.0501<br>0.202<br>-0.325<br>-1.163</td><td>+42:22:49.124<br>0.048<br>-0.133<br>-0.226</td><td>-2.72<br>-3.08<br>-3.80<br>8.98</td><td>4.07 ± 0.05<br>2.72 ± 0.05<br>1.30 ± 0.03<br>1.22 ± 0.03</td><td>4.30 ± 0.10<br>3.12 ± 0.10<br>1.43 ± 0.05<br>1.25 ± 0.06</td><td>DR21(OH)</td></tr><tr><td>10</td><td>G81.7444+0.5910</td><td>20:39:00.3722<br>0.015</td><td>+42:24:37.089<br>0.020</td><td>4.48<br>3.58</td><td>17.18 ± 0.12<br>6.58 ± 0.08</td><td>17.88 ± 0.22<br>6.61 ± 0.14</td><td>W75S-FIR1</td></tr><tr><td>11</td><td>G81.7523+0.5908</td><td>20:39:01.9870<br>-0.014<br>-0.036<br>0.123<br>0.074</td><td>+42:24:59.261<br>-0.140<br>-0.107<br>-0.300<br>-0.363</td><td>-8.66<br>-5.78<br>-6.86<br>-2.36<br>-2.90</td><td>15.19 ± 0.10<br>7.93 ± 0.07<br>4.29 ± 0.04<br>1.33 ± 0.03<br>1.52 ± 0.03</td><td>15.67 ± 0.19<br>8.12 ± 0.13<br>4.37 ± 0.08<br>1.45 ± 0.06<br>1.70 ± 0.07</td><td>W75S-FIR2</td></tr><tr><td>12</td><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>-1.28</td><td>3.90 ± 0.05</td><td>4.09 ± 0.09</td><td>CygX-N53</td></tr><tr><td>13</td><td>G81.8713+0.7807</td><td>20:38:36.4232<br>-0.176<br>-0.170<br>-0.261<br>-0.241<br>-0.100<br>-0.122</td><td>+42:37:34.756<br>0.425<br>0.444<br>0.346<br>0.352<br>0.480<br>-0.424</td><td>7.18<br>4.66<br>4.48<br>5.74<br>5.20<br>3.40<br>9.34</td><td>336.13 ± 0.96<br>324.83 ± 1.15<br>238.51 ± 1.07<br>126.12 ± 0.76<br>114.47 ± 0.74<br>88.83 ± 0.80<br>31.67 ± 0.61</td><td>346.76 ± 1.73<br>336.46 ± 2.08<br>246.98 ± 1.95<br>128.52 ± 1.37<br>116.80 ± 1.33<br>90.64 ± 1.43<br>35.80 ± 1.16</td><td>W75N(B)</td></tr><tr><td>14</td><td>G82.3079+0.7296</td><td>20:40:16.6487</td><td>+42:56:29.288</td><td>10.42</td><td>23.82 ± 0.18</td><td>24.19 ± 0.31</td><td>—</td></tr></tbody></table>

**Notes.** Column 2 gives source Galactic coordinates; Cols. 3 and 4 are equatorial coordinates. We give the absolute position of the maser component with maximum intensity. For the other components, we list position offsets relative to the position of the strongest component. The position uncertainties are ±0.′′18 and ±0.′′13, in RA and Dec, respectively (see Sect. 3). For each feature, Col. 5 gives the LSR radial velocity of the peak. Columns 6 and 7 are the peak and integrated fluxes at the channel corresponding to the peak of the feature. Column 8 gives the common source names used in the literature.

we do not show the spectra from the convolved data cubes for the weakest lines (those with $S_{\text{peak}} < 1$ Jy beam$^{-1}$) in Fig. 2. At our 0.18 km s$^{-1}$ channel spacing, even narrow single velocity features are generally spectrally resolved. In a few cases, extremely narrow features are unresolved or only partially resolved. In the case of strong features, this leads to observable “Gibbs ringing,” which appears most pronounced in the G82.3079+0.7296 spectrum between 7 and 10 km s$^{-1}$ but can also be seen in a few of the other maser spectra. The vertical lines show the systemic velocities of the molecular gas in the line of sight, which were taken from the literature. The LSR velocity of all but one of the maser features ranges from -12 to +10 km s$^{-1}$. The exception is source G78.9884+0.2211, which has $V_{\text{LSR}} \approx -70$ km s$^{-1}$. In all cases, methanol-maser velocities are observed to be close, within $\le 7$ km s$^{-1}$, to the systemic velocities.

We see small differences in the maser emission between D- and B-configuration data. In general, fewer maser features

(i.e., emission observed in contiguous velocity channels at nearly the same position) are seen in the B-configuration images. The peak fluxes are also weaker in the B-configuration data for five sources; for these sources we see variations in the range of 20–90%, while the velocity profiles do not change significantly (Fig. 2). This variability is further discussed for individual sources in Sect. 4.4.

Maser luminosities were estimated using

$$L_{\text{maser}} = 4\pi D^2 S_{\text{int}} f/c, \qquad (1)$$

where $S_{\text{int}}$, in units of Jy km s$^{-1}$, is the maser flux integrated over the velocity range $\Delta V$ (i.e., across the maser line with emission above $4\sigma$), $f = 6668.5192$ MHz is the rest frequency of the $S_1-6_0 A^+$ methanol line, $c$ is the speed of light, and $D$ is the distance to the source. The term $f/c = \Delta f/\Delta v$ takes the change in frequency, $\Delta f$, into account, resulting from a change in velocity
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Observed spectra of the 13 sources with methanol maser emission and the one with methanol absorption. For every source, the top and bottom panels represent D- and B-configuration data, respectively. These spectra were extracted at the peak pixel from the data cubes (see text). The vertical lines show the systemic LSR velocities of the dense molecular gas of the star-forming regions that host the masers, which were taken from the literature (see Table 6).</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>Fig. 2. continued.</figcaption></figure>

of $\Delta v = 0.18$ km s$^{-1}$. The velocity ranges as well as the resulting values of the maser luminosity are given in Table 4. For all sources (except AFGL 2591) we used a distance of $D = 1.4$ kpc, which is the average distance of the Cygnus X complex measured from trigonometric parallaxes of several maser sources (Rygl et al. 2012; see also Dzib et al. 2013, who obtain a parallax measurement to a radio continuum source in the Cygnus OB2 association). On the other hand, source AFGL 2591 is much farther away, at 3.33 kpc (Rygl et al. 2012), based also on maser parallax measurements. Figure 3 shows the velocity range of the maser emission as a function of integrated maser luminosity. As expected, the maser velocity range (which can be used as a proxy for the line profile complexity) increases with maser luminosity. A similar relation was previously noticed by, for example, Breen et al. (2011) and Billington et al. (2019). This may simply be the result of the fact that for a maser whose strongest component has a low flux density, any additional even weaker velocity components are below the S/N cutoff.

## 4.2. Association with millimeter-to-submillimeter emission

To compare the 6.7 GHz methanol maser activity with other tracers of star formation, we first compared the maser positions

with those of the submillimeter continuum sources that pinpoint the high density gas. For this, we used JCMT Submillimeter Common-User Bolometer Array 2 (SCUBA-2) maps of Cygnus X at 450 and 850 micron (Cao et al. 2019). We found that all detected masers have associated 850 µm emission (Fig. 4) with a comparable beam size of 13'' (cf. Table 1 in Cao et al. 2019). Thus, all our methanol masers are associated with dust continuum emission. This is consistent with the results of Urquhart et al. (2015) and Billington et al. (2019), who found an almost ubiquitous association (99%) between the 6.7 GHz methanol maser and dust continuum sources. In order to derive the properties (mass, temperature, and luminosities) of these dust cores, we fit the spectral energy distribution (SED) of the dust emission with a modified blackbody model,

$$F_v = \frac{\kappa_v B_v(T) M}{D^2}, \qquad (2)$$

where $B_v(T)$ is the Planck function for the temperature, $T$, of the dust and $M$ is the core mass. Here, we adopted a dust mass opacity of $\kappa_v = \kappa_0(v/v_0)^\beta$, where $\kappa_0 = 0.1 \text{ cm}^2 \text{ g}^{-1}$, $v_0 = 1$ THz, and $\beta = 2$. For the construction of the SED (see Fig. B.1), we used data from Herschel obtained as part of the HOBYS (Motte et al. 2010; Hennemann et al. 2014) and Hi-GAL projects
<table><caption>**Table 4.** Properties of the cores associated with methanol masers.</caption><thead><tr><th>GLOSTAR name (1)</th><th>Common name (2)</th><th>ΔV<sub>D</sub> (km s<sup>-1</sup>) (3)</th><th>ΔV<sub>B</sub> (km s<sup>-1</sup>) (4)</th><th>L<sub>maser, D</sub> (L<sub>☉</sub>) (5)</th><th>L<sub>maser, B</sub> (L<sub>☉</sub>) (6)</th><th>Size (′′) (7)</th><th>M<sub>core</sub> (M<sub>☉</sub>) (8)</th><th>T<sub>core</sub> (K) (9)</th><th>L<sub>FIR</sub> (L<sub>☉</sub>) (10)</th><th>L<sub>Bol</sub> (L<sub>☉</sub>) (11)</th></tr></thead><tbody><tr><td>G76.0932+0.1580</td><td>–</td><td>1.3</td><td>1.6</td><td>5.8 × 10<sup>-9</sup></td><td>5.9 × 10<sup>-9</sup></td><td>18</td><td>27 ± 3</td><td>20.4 ± 0.9</td><td>189</td><td>227</td></tr><tr><td>G76.8437+0.1233</td><td>–</td><td>2.2</td><td>0.5</td><td>3.8 × 10<sup>-8</sup></td><td>1.0 × 10<sup>-9</sup></td><td>25</td><td>12 ± 2</td><td>16.0 ± 0.7</td><td>19</td><td>26</td></tr><tr><td>G78.8870+0.7087</td><td>AFGL2591</td><td>0.7</td><td>0.9</td><td>2.3 × 10<sup>-8</sup></td><td>2.4 × 10<sup>-8</sup></td><td>24</td><td>937 ± 72</td><td>30.2 ± 1.9</td><td>67 750</td><td>66 262</td></tr><tr><td>G78.9690+0.5410</td><td>IRAS 20286+3959</td><td>1.4</td><td>1.4</td><td>1.3 × 10<sup>-8</sup></td><td>1.4 × 10<sup>-8</sup></td><td>22</td><td>72 ± 9</td><td>13.0 ± 0.4</td><td>33</td><td>42</td></tr><tr><td>G78.9884+0.2211<sup>(a)</sup></td><td>–</td><td>3.6</td><td>0.7</td><td>1.8 × 10<sup>-8</sup></td><td>2.9 × 10<sup>-9</sup></td><td>25</td><td>8 ± 1</td><td>21.3 ± 1.2</td><td>69</td><td>176</td></tr><tr><td>G79.7358+0.9904</td><td>IRAS 20290+4052</td><td>3.4</td><td>3.2</td><td>2.4 × 10<sup>-7</sup></td><td>2.4 × 10<sup>-7</sup></td><td>19</td><td>42 ± 4</td><td>17.7 ± 0.5</td><td>122</td><td>233</td></tr><tr><td>G80.8617+0.3834</td><td>DR20</td><td>4.7</td><td>1.4</td><td>1.0 × 10<sup>-7</sup></td><td>3.7 × 10<sup>-8</sup></td><td>23</td><td>122 ± 11</td><td>20.9 ± 0.8</td><td>972</td><td>1328</td></tr><tr><td>G81.7219+0.5711</td><td>DR21(OH)</td><td>2.3</td><td>2.5</td><td>5.2 × 10<sup>-8</sup></td><td>4.9 × 10<sup>-8</sup></td><td>16</td><td>930 ± 140</td><td>20.9 ± 1.3</td><td>7344</td><td>7299</td></tr><tr><td>G81.7444+0.5910</td><td>W75S-FIR1</td><td>2.2</td><td>2.2</td><td>1.3 × 10<sup>-7</sup></td><td>1.3 × 10<sup>-7</sup></td><td>11</td><td>74 ± 5</td><td>21.9 ± 0.6</td><td>776</td><td>1487</td></tr><tr><td>G81.7523+0.5908</td><td>W75S-FIR2</td><td>5.4</td><td>5.4</td><td>1.8 × 10<sup>-7</sup></td><td>2.0 × 10<sup>-7</sup></td><td>9</td><td>298 ± 34</td><td>16.4 ± 0.6</td><td>552</td><td>955</td></tr><tr><td>G81.7655+0.5972</td><td>CygX-N53</td><td>1.3</td><td>1.6</td><td>2.3 × 10<sup>-8</sup></td><td>2.2 × 10<sup>-8</sup></td><td>11</td><td>325 ± 27</td><td>14.5 ± 0.3</td><td>294</td><td>221</td></tr><tr><td>G81.8713+0.7807</td><td>W75N(B)</td><td>10.3</td><td>7.2</td><td>5.8 × 10<sup>-6</sup></td><td>5.8 × 10<sup>-6</sup></td><td>–</td><td>1020 ± 91</td><td>27.3 ± 1.3</td><td>40 588</td><td>32 930</td></tr><tr><td>G82.3079+0.7296</td><td>–</td><td>2.3</td><td>1.3</td><td>1.2 × 10<sup>-7</sup></td><td>1.3 × 10<sup>-7</sup></td><td>16</td><td>24 ± 3</td><td>19.7 ± 0.8</td><td>139</td><td>183</td></tr></tbody></table>

**Notes.** Columns 1 and 2 report the maser name, which is based on Galactic coordinates, and its common label; Cols. 3 and 4 are the velocity ranges of the maser emission in D- and B-configurations; Cols. 5 and 6 are the maser integrated luminosities; Col. 7 is the deconvolved size at 160 µm band, where the dust is optically thin. It corresponds to the intensity-weighted diameter given by $\theta_D = 2 \sum_i r_i A_i / \sum_i A_i$, where $r_i$ is the distance of the ith pixel to the center of the source and $A_i$ is the pixel intensity (Purcell et al. 2013); Cols. 8–11 are the mass, temperature, integrated FIR luminosity, and bolometric luminosity of the associated dust cores. $^{(a)}$The mass, FIR luminosity, and bolometric luminosity of this core must be rescaled to 347 $M_☉$, 3093 $L_☉$, and 6237 $L_☉$, respectively, for a distance of 9.35 kpc (see Sect. 4.4).

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Velocity range of detected maser emission in Cygnus X as a function of maser integrated luminosity. The colors indicate D- and B-array configurations (blue and orange, respectively).</figcaption></figure>

(Molinari et al. 2010) in the Photodetector Array Camera and Spectrometer (PACS) 70 and 160 µm and the Spectral and Photometric Imaging Receiver (SPIRE) 250, 350, and 500 µm bands<fnref n="3" />, as well as JCMT SCUBA-2 data at 450 and 850 µm (Cao et al. 2019). In addition, data at 1.1 mm are available from the Bolocam instrument installed at the 10.4 m Caltech Submillimeter Observatory, which have a resolution of 33" and were taken as part of the Bolocam Galactic Plane Survey (BGPS; Aguirre et al. 2011; Rosolowsky et al. 2010). At a higher angular resolution of 11", 1.2 mm data are also available from the Max-Planck-Millimeter-Bolometer (MAMBO/MAMBO-2) instruments installed at the Institut de Radioastronomie Millimétrique (IRAM) telescope

(Motte et al. 2007). Not all dust cores associated with the methanol sources were covered by the observations of Motte et al. (2007). Thus, we used 1.2 mm IRAM data for eight sources and 1.1 mm Bolocam data for six sources. Figure 4 displays the images of (sub)millimeter continuum data for the GLOSTAR 6.7 GHz methanol masers. The images are 3' × 3' in size and are centered at the positions of the masers.

In the SED we also show infrared data from the Infrared Array Camera (IRAC) 3.6, 4.5, 5.8, and 8 µm and the Multi-band Imaging Photometer for *Spitzer* (MIPS) 24 µm bands taken as part of the *Spitzer* Legacy Survey of the Cygnus X Complex (Beerer et al. 2010). IRAC aperture photometry and MIPS point spread function (PSF) photometry were retrieved from the NASA/IPAC Infrared Science Archive (IRSA4). The corresponding images can be seen in Fig. 5, where we present both three-color infrared images and single band images of 0.7 × 0.7 in size around the location of the masers.

The mass and luminosity of five maser-associated cores are not available in the literature, while for the other eight these properties have been measured, for example by Cao et al. (2019). We thus derived these properties from the SED fitting not only for the five missing cores, but – for consistency and to avoid systematic deviations introduced by different extraction techniques – for all cores that have associated methanol maser emission detected by GLOSTAR. When fitting the SED, we note that the point at 70 µm shows, in almost all cases, significant excess emission to the fit since emission at this wavelength traces warm dust that cannot be reproduced with a single-temperature model (Cao et al. 2019; König et al. 2017). The 70 and 24 µm data points can be reproduced if we add a second, warmer component to the SED fit (see Fig. B.1). We thus carried out a two-component fit for sources for which 70 and 24 µm data are both available. We fit Eq. (2) to 24–70 µm and 160–1000 µm data for the warm and cold component, respectively.

<footnote n="3">³ Data were retrieved from the Canadian Astronomy Data Centre.</footnote>

<footnote n="4">⁴ https://irsa.ipac.caltech.edu/data/SPITZER/Cygnus-X/</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Dust continuum emission images of DR20 from: Herschel PACS (70 and 160 µm) and SPIRE (250, 350, and 500 µm); JCMT SCUBA-2 (450 and 850 µm); and either IRAM MAMBO/MAMBO-2 (1.2 mm) or CSO Bolocam (1.1 mm) instruments. The color scale is in units of Jy beam⁻¹. The VLA position of the detected methanol maser is indicated by the black cross. The galactic coordinate name of the GLOSTAR source is indicated at the top of the *bottom-left panel*. The beams are shown at the bottom-left corner of each panel. Saturated pixels are shown in gray. The same figure for all sources is provided in Appendix A.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Infrared images (color scale and white contours) of the environment of the VLA-detected methanol masers toward DR20. **Top-left panel:** three-color map constructed using *Spitzer* 3.6 µm (blue), 4.5 µm (green), and 8 µm (red) images. **Other five panels:** infrared emission for each of the *Spitzer* bands. The nth white contour is at (√2)ⁿ × Smax × p, where Smax is the maximum flux shown in the color bar (MJy sr⁻¹) for each panel, n=0, 1, 2, ..., and p is equal to 30%. The green crosses mark the positions of the maser features seen in B-array configuration maps. The Galactic coordinate name of the GLOSTAR source is indicated at the top of the *top-left panel*. The beams are shown at the bottom-left corner of each panel. Saturated pixels are shown in gray. The same figure for all sources is provided in Appendix A.</figcaption></figure>

The fit was done with the Python routine *curve_fit* from the *scipy* package. The fluxes from all images that are saturated at the location of the masers (cf. Figs. A.1 and A.2) are not shown in the SED and were thus not used in the fitting. Specifically, saturation occurs for G78.8870+0.7087⁵ (24 and

250 µm), G80.8617+0.3834 (24 µm), G81.7219+0.5711 (250 and 350 µm), G81.7444+0.5910 (24 µm), G81.7523+0.5908 (24 µm), and G81.8713+0.7807 (24, 160, 250, and 350 µm). In the case of G78.8870+0.7087 and G81.8713+0.7807, for which we have a limited number of data points due to saturation, the 70 µm point was included in the fit of the cold component, and the data were weighted with their 20% flux errors. The IRAC data points between 3.6 and 8 µm – likely originating from a third, hot and inner component – were not used in the SED fitting.

<footnote n="5">⁵ The image of this source has issues in the 3.6, 4.5, 5.8, and 8 µm bands as well. Thus, the fluxes from these bands are not shown in the SED.</footnote>
Finally, far-infrared (FIR) luminosities were obtained by integrating the SED of the cold component over wavelength,

$$L_{\text{FIR}} = 4\pi D^2 \int_{0}^{+\infty} S_{\lambda} d\lambda. \qquad (3)$$

To obtain bolometric luminosities from the SED, we should also have integrated the emission from the near- and mid-infrared; however, these data are not available for all our sources (cf. Fig. A.2). Instead, we estimated the bolometric luminosities of the cores by using the empirical relation between luminosity and 70 µm flux found by Dunham et al. (2008),

$$L_{\text{bol}} = 3.3 \times 10^8 \left( \left( \frac{\nu F_{70 \mu\text{m}}}{\text{ergs cm}^{-2} \text{s}^{-1}} \right) \left( \frac{D}{140 \text{ pc}} \right)^2 \right)^{0.94} L_{\odot}. \qquad (4)$$

The derived parameters for the cold components (masses, FIR luminosities, and bolometric luminosities) are listed in Table 4, and the SED and model fit are shown in Fig. B.1. We found that the two methods give luminosities that agree within $1\sigma$ in eleven sources and within $3\sigma$ in two sources. We now compare our results with those previously published in the literature.

For the core extraction, Cao et al. (2019) selected regions with column density $N_{H_2} \ge 3.5 \times 10^{22} \text{cm}^{-2}$, which, for typical deconvolved sizes of 0.1 pc, corresponds to core masses above 23 $M_{\odot}$. Additional criteria were used by Cao et al. (2019) to reject cores with poorly extracted fluxes and/or large relative errors in mass and temperature estimation from the SED fitting, resulting in a minimum mass of 36 $M_{\odot}$ in their catalog. We see that some methanol maser-associated cores not reported previously have masses close to the threshold of 23 $M_{\odot}$ adopted by Cao et al. (2019), which explains why they were not reported by these authors. Although we follow a similar approach to Cao et al. (2019) to estimate core properties, here the fluxes were extracted individually at each wavelength using the radio flux plug-in in DS9⁶. The integrated source flux was obtained over an aperture that delimits emission within a 30% contour. The local background flux was estimated over an aperture farther out from the source and then subtracted from the source flux. We note that fluxes extracted in this way do not show large deviations from the model fit, unlike some fluxes extracted with the getsources package used by Cao et al. (2019), especially at 450 µm. The (sub)millimeter emission in the DR21 ridge is complex and extended, which probably affects the performance of background subtraction by getsources. As a result, the masses derived by us agree with those from Cao et al. (2019) within ~55%.

When comparing the estimated masses with the masses obtained by Motte et al. (2007) for the eight overlapping cores from 1.2 mm data alone, after correcting for the different distance and dust mass opacity used by these authors, we see a slightly better agreement, with a relative difference of ~50% on average. The properties derived above for the dust cores that are associated with a methanol maser source are used in Sect. 5 to investigate the relationship between methanol masers and massive star-forming dust cores.

### 4.3. Comparison with other maser surveys

Pestalozzi et al. (2005) published a catalog of all known 6.7 GHz methanol masers from both targeted and unbiased surveys in

the Galactic plane. Their catalog contains six masers spread throughout the region observed by GLOSTAR in Cygnus X. Out of these six masers, four have a GLOSTAR counterpart within 1', which is the average of the uncertainties in the maser coordinates of Pestalozzi's catalog. The sources from Pestalozzi et al. (2005) without counterparts are G78.62+0.98 and G80.85+0.43. Source G78.62+0.98 has a quoted peak flux of 3 Jy ($\sigma=1$ Jy) at $V_{LSR} = -39.0 \text{ km s}^{-1}$. This velocity is very different from the velocities of our masers detected with GLOSTAR. The source G80.85+0.43 has a quoted peak flux of 4.3 Jy at $V_{LSR} = -4.1 \text{ km s}^{-1}$. This source is 2.9 away from a strong maser detected by GLOSTAR (G80.8617+0.3834/DR20), which shows a peak flux of 12 Jy ($\sigma=0.12$ Jy) at $V_{LSR} = -4.16 \text{ km s}^{-1}$ and is not reported by Pestalozzi et al. (2005). The two masers G80.85+0.43 and G80.8617+0.3834/DR20 could correspond to the same source given that the Pestalozzi et al. (2005) catalog collects positions obtained with a variety of single dish observations and hence different positional accuracy. They report an uncertainty of up to 3' for their Galactic coordinates that would suggest that these two sources are the same. We created cleaned images centered at the position of sources G80.85+0.43 and G78.62+0.98 and did not see emission. These masers are not detected in our Effelsberg maps either.

Another targeted survey of this methanol maser line was conducted by Hu et al. (2016) with the VLA in C configuration. They selected all known methanol masers visible from the northern hemisphere and with peak fluxes above 2 Jy. The eight maser sources listed in this catalog within $76^{\circ} < l < 83^{\circ}$ and $-1^{\circ} < b < +2^{\circ}$ were also detected in our survey. More recently, Yang et al. (2019) performed a new targeted survey in the whole Galaxy with the 65 m Shanghai Tianma Radio Telescope. They observed targets from the all-sky Wide Field Infrared Survey Explorer (WISE) point source catalog that could be associated with 6.7 GHz methanol masers based on their WISE magnitudes. We detected their three new masers and the five already known masers with WISE bright emission listed in Yang et al. (2019). The only maser from this catalog not detected by us, likely due to variability, is G81.794+0.911. However, we do not see this maser in recent Effelsberg data either. In summary, out of thirteen GLOSTAR masers, four were previously reported by both Pestalozzi et al. (2005) and Hu et al. (2016), five by both Hu et al. (2016) and Yang et al. (2019), three others by Hu et al. (2016), and three others by Yang et al. (2019). Two maser sources (G76.8437+0.1233 and G78.9884+0.2211) are newly detected by our observations. These two sources display high flux variations between the two observed epochs, which may explain why they were missed by previous surveys. Names, positions, and fluxes of all velocity components detected in the surveys by Hu et al. (2016) and Yang et al. (2019) are listed for reference in Table C.1.

### 4.4. Comments on individual sources and their possible nature

In this section we give a short description of each 6.7 GHz methanol source detected in Cygnus X (here, sources are referred to by their GLOSTAR identifier, GLL.111+/-BB.bbbb, and also by their most common names if they are known sources). In order to discuss the association between these sources and their radio continuum emission, we show in Fig. 6 the positions of maser components determined from the B-configuration images, overlaid on GLOSTAR radio continuum emission that is also from B-configuration data. Maser positions measured from D-configuration images are less precise as compared with B-configuration positions, and thus they do not provide additional

<footnote n="6">⁶ https://www.extragalactic.info/~mjh/radio-flux.html</footnote>
<figure><img src="image_8.png" /></figure>

**Fig. 6.** Positions of maser features (filled circles) overlaid on continuum emission (blue scale and blue contours) for B-array data. The masers are color coded by LSR velocity (color bar). When only one velocity component is detected, the velocity is given at the bottom of the panel. Names and positions of radio continuum sources identified in previous observations are also marked (black stars). The *n*th blue contour is at (√2)*n* × *S*<sub>max</sub> × *p*, where *S*<sub>max</sub> is the maximum continuum flux shown in the blue bar for each panel, *n*=0, 1, 2 ..., and *p* is equal to 10%. The magenta rectangle in panel 1 depicts the approximate extent of a 1.3 mm continuum source detected with NOEMA (Gieser et al. 2019). The gray grid shows the Galactic coordinate system. The beam size of the radio continuum maps is 1′′5 × 1′′5. The maser position errors are ≈0′′2 (see Sect. 3).
<table><caption>**Table 5.** Measured properties of radio continuum sources detected in GLOSTAR B-array maps close to the methanol line sources.</caption><thead><tr><th>Name<br>(1)</th><th><i>l</i><br>(°)<br>(2)</th><th><i>b</i><br>(°)<br>(3)</th><th>α<br>(J2000)<br>(4)</th><th>δ<br>(J2000)<br>(5)</th><th><i>S</i><sub>Peak</sub><sup>(a)</sup><br>(mJy beam<sup>-1</sup>)<br>(6)</th><th>SNR<br>(7)</th><th>Offset<sup>(b)</sup><br>(′′)<br>(8)</th></tr></thead><tbody><tr><td>G78.8870+0.7087</td><td>78.886737</td><td>0.708870</td><td>20h29m24.85s</td><td>+40d11m19.7s</td><td>1.69 ± 0.37</td><td>4.5</td><td>1.2</td></tr><tr><td>G80.8617+0.3834</td><td>80.861469</td><td>0.383462</td><td>20h37m00.91s</td><td>+41d34m56.0s</td><td>0.79 ± 0.13</td><td>6.0</td><td>0.6</td></tr><tr><td>G81.6790+0.5378<sup>(c)</sup></td><td>81.678758</td><td>0.537054</td><td>20h39m01.30s</td><td>+42d19m31.5s</td><td>1210 ± 4<sup>(d)</sup></td><td>333</td><td>2.4</td></tr><tr><td>G81.7219+0.5711</td><td>81.721609</td><td>0.571246</td><td>20h39m00.95s</td><td>+42d22m48.8s</td><td>3.05 ± 0.40</td><td>7.6</td><td>0.3</td></tr><tr><td>G81.7444+0.5910</td><td>81.744319</td><td>0.591176</td><td>20h39m00.30s</td><td>+42d24m37.4s</td><td>1.69 ± 0.22</td><td>7.8</td><td>0.8</td></tr><tr><td>G81.8713+0.7807</td><td>81.871218</td><td>0.780755</td><td>20h38m36.39s</td><td>+42d37m34.9s</td><td>4.44 ± 0.08</td><td>58.2</td><td>0.4</td></tr><tr><td>G76.0932+0.1580</td><td>—</td><td>—</td><td>—</td><td>—</td><td>&lt;0.16</td><td>—</td><td>—</td></tr><tr><td>G76.8437+0.1233</td><td>—</td><td>—</td><td>—</td><td>—</td><td>&lt;0.18</td><td>—</td><td>—</td></tr><tr><td>G78.9690+0.5410</td><td>—</td><td>—</td><td>—</td><td>—</td><td>&lt;0.20</td><td>—</td><td>—</td></tr><tr><td>G78.9884+0.2211</td><td>—</td><td>—</td><td>—</td><td>—</td><td>&lt;0.33</td><td>—</td><td>—</td></tr><tr><td>G79.7358+0.9904</td><td>—</td><td>—</td><td>—</td><td>—</td><td>&lt;0.22</td><td>—</td><td>—</td></tr><tr><td>G81.7523+0.5908</td><td>—</td><td>—</td><td>—</td><td>—</td><td>&lt;0.63</td><td>—</td><td>—</td></tr><tr><td>G81.7655+0.5972</td><td>—</td><td>—</td><td>—</td><td>—</td><td>&lt;0.45</td><td>—</td><td>—</td></tr><tr><td>G82.3079+0.7296</td><td>—</td><td>—</td><td>—</td><td>—</td><td>&lt;0.19</td><td>—</td><td>—</td></tr></tbody></table>

**Notes.** The beam size of the continuum maps is $1.5 \times 1.5$°. $^{(a)}$A $3\sigma$ upper limit is given for the cases where no radio continuum emission is detected around the position of methanol masers. $^{(b)}$Offset in position of the radio continuum source from the methanol maser intensity peak. $^{(c)}$Methanol line is seen in absorption. $^{(d)}$This flux was measured in the map made for the combination of D- and B-configuration data.

information on the distribution of maser spots. For this reason, we do not show D-configuration positions in Fig. 6.

A clear association between maser (colored circles) and continuum emission (blue scale) and within a radius of ~1" is seen in five cases. The coordinates, peak fluxes, and the S/N of the detected radio continuum sources as measured with BLOBCAT (Hales et al. 2012; Medina et al. 2019) are given in Table 5. A $3\sigma$ upper limit is given in the table when no radio continuum emission is detected.

**G78.8870+0.7087/AFGL2591.** This maser is close (~1"; cf. Table 5) to the VLA3 hot core (Campbell 1984; Johnston et al. 2013; Gieser et al. 2019) in the well-known high-mass star-forming region AFGL2591. VLA3 is also detected in the GLOSTAR B-configuration map (see panel 1 of Fig. 6). A spectral index of $0.5 \pm 0.02$ between 8.3 and 43 GHz was estimated by Johnston et al. (2013) for VLA3. Extended bright infrared emission at 3.6–24 µm is also present in the environment of the maser (see Figs. A.2 and 7). However, this infrared source is saturated in the *Spitzer* images. Motte et al. (2007) identified a dust core in their 1.2 mm observations whose peak position is closest to VLA3, and this dust core was recognized as a UC HII region. More recently, Gieser et al. (2019) resolved the millimeter emission toward VLA3 with higher angular observations taken with the NOrthern Extended Millimeter Array (NOEMA) and revealed an extended, almost spherically symmetric envelope at ~2" scales. We note that both the maser and the radio continuum emission lie within the extension of this source, which is approximately depicted by the magenta circle in panel 1 of Fig. 6. The integrated flux of this maser in the peak channel changes from $0.80 \pm 0.04$ Jy in the D-configuration image to $0.29 \pm 0.04$ Jy in the B-configuration image due to potential variability. The Effelsberg data, however, show an integrated flux of ~0.7 Jy, which is consistent with the flux seen in the D-configuration map. Its 3.3 kpc maser parallax distance places AFGL 2591, like Cygnus X, in the local arm (see, e.g., Fig. 1 of Reid et al. 2019).

**G80.8617+0.3834/DR20.** We detected two additional features in the D-configuration spectrum, at −12.3 and 1 km s$^{-1}$,

which were not seen in the spectrum from Hu et al. (2016) and Yang et al. (2019). We note that these two components and the component at −11 km s$^{-1}$ are not detected in the B-configuration data, likely due to variability. We see a weak source ($S_{\text{peak}}=0.79$ mJy beam$^{-1}$) in the B-array radio continuum map coincident with the location of the maser (panel 2 of Fig. 6). The infrared environment around the maser shows slightly extended 3.6–8 µm emission (Fig. A.2). The source at 24 µm is saturated. This object is recognized as a high-luminosity infrared protostellar core (Motte et al. 2007).

**G81.7219+0.5711/DR21(OH).** Compared to the spectrum from Hu et al. (2016), we detected an additional weak feature at 9 km s$^{-1}$. A bright infrared peak is seen at 3.6–24 µm close to the location of the maser (Fig. A.2); however, its position does not exactly coincide with the maser intensity peak, showing a small offset of ~1". This maser is associated with the well-known high-mass star-forming region DR21(OH). Its position coincides with the main core in DR21(OH), which is a massive infrared-quiet (<10 Jy at 21 µm) protostar in an early evolutionary stage (Motte et al. 2007). We see a weak source in the B-array radio continuum map close to the location of the masers (panel 3 of Fig. 6). The strongest maser component is offset by 1"1 with respect to the peak position of the continuum source; however, the component at $V_{\text{LSR}}=9$ km s$^{-1}$ coincides in position with the continuum source. The continuum source has a peak flux density of $S_{\text{peak}}=3$ mJy beam$^{-1}$. In a previous work, Araya et al. (2009) conducted radio continuum observations at 3.6, 1.3, and 0.7 cm toward DR21(OH) at subarcsecond angular resolution. They detected a cluster of radio sources; their source MM1-NW ($S_{\text{peak}}=1.14 \pm 0.02$ mJy beam$^{-1}$ at 8.46 GHz) is the closest to the source we detect in our GLOSTAR B-array continuum data. Source MM1-NW is also coincident with the maser at $V_{\text{LSR}}=9$ km s$^{-1}$, while source R2 is the closest to the masers at smaller velocities (see panel 3 of Fig. 6). Araya et al. (2009) suggested that the two brightest sources detected toward the molecular core (MM1-NW and MM1-SE) trace radio emission from shock-ionized gas in a jet, while the radio emission from source R2 was proposed to be gyrosynchrotron radiation
<figure><img src="image_9.png" /><figcaption>**Fig. 7.** MSX infrared images (color scale and contours) of the environment around the location of G78.8870+0.7087 and G81.8713+0.7807. The *Spitzer* and WISE images of these two sources are saturated in at least one band. The *first panel* in each row shows a three-color map constructed using MSX 8.28 µm (blue), 12.13 µm (green), and 14.65 µm (red) images. *Other four panels:* infrared emission for each of the four MSX bands. The *n*th white contour is at (√2)ⁿ × *S*<sub>max</sub> × *p*, where *S*<sub>max</sub> is the maximum flux shown in the color bar (W m⁻² sr⁻¹) for each panel, *n* = 0, 2, 4 ..., and *p* is equal to 3%. The green crosses mark the positions of the masers as measured in D-array configuration maps. The Galactic names of the GLOSTAR sources are indicated at the bottom of the second image in each row.</figcaption></figure>

from a low-mass star. We cannot resolve these individual contin-
uum components in the GLOSTAR data, which prevents us from
estimating spectral indices of the individual radio sources to con-
firm or reject the conclusion by Araya et al. (2009). Based on
their higher angular resolution images, Araya et al. (2009) mea-
sured a spectral index of 0.8 ± 0.2 between 8.4 and 43 GHz for
both MM1-NW and MM1-SE, and 0.1 ± 0.2 for R2. Thus, given
the correspondence in the positions measured by these authors
with the maser positions, it is tempting to suggest that MM1-NW
and R2 are UC HII regions.

G81.7444+0.5910/W75S-FIR1. We see a weak source in the B-array radio continuum map very close (at 0′′8) to the location of the maser position (panel 4 of Fig. 6). The continuum source has a peak flux density of $S_{\text{peak}} = 1.7$ mJy. The maser has associated compact emission at 3.6–8 µm (Fig. A.2). At 24 µm the source is saturated, but we see some extended emission surrounding the maser. The source with which this maser is associated is the high-luminosity ($\gtrsim 10^3 L_{\odot}$) infrared object W75S FIR1 (Harvey et al. 1986), a well-known massive embedded protostar (Chandler et al. 1993; Motte et al. 2007).

G81.8713+0.7807/W75N(B). We see that the masers trace a velocity gradient of a few km s⁻¹ in the north-south direction across a scale of ~2". We detect continuum emission close to this maser in both D-array and B-array VLA maps. The maser position is within 0′′4 of a radio continuum peak ($S_{\text{peak}} = 4.4$ mJy) seen in the B-array map (panel 5 of Fig. 6). The radio continuum peak corresponds to VLA 1, a source that was found to have a jet-like appearance at 1.3 cm (Torrelles et al. 1997). Recently, Rodríguez-Kamenetzky et al. (2020) conducted multifrequency VLA observations and concluded that VLA1 consists of a thermal radio jet surrounded by a hyper-compact (HC) HII region. The infrared environment toward the maser shows three peaks in

the 3.6–8 µm maps (Fig. A.2). At 24 µm, the source is saturated. The maser seems to be located between two infrared peaks seen in the 5.8 µm map, but its position coincides with an infrared peak we see at 21.3 µm in the Midcourse Space Experiment (MSX) map (see the bottom panel of Fig. 7). The maser is found in the star-forming region W75N(B), which is mainly composed of two UC HII regions (Hunter et al. 1994), namely W75N(Ba) and W75N(Bb) (also known as VLA1 and VLA3, respectively, in the nomenclature of Torrelles et al. 1997).

The following maser sources do not have counterparts in
either our VLA radio continuum maps or previous radio con-
tinuum surveys (within a matching radius of 10").

G76.0932+0.1580. Yang et al. (2019) reported on the veloc-
ity components detected at 4.8 and 6.4 km s⁻¹ (see Table C.1).
Here, we detect an additional component at −6.5 km s⁻¹ in the
D-array VLA spectrum. The maser is associated with a compact
weak source seen at 4.5–24 µm (Fig. A.2). The infrared source
has been classified as a protostar candidate by Kryukova et al.
(2014) and shows a rising infrared SED.

G76.8437+0.1233. This source is a newly detected 6.7 GHz maser. The maser is associated with compact bright emission at 4.5–24 µm (Fig. 5). It also appears in the catalog of protostar candidates by Kryukova et al. (2014) as a rising SED protostar. One single strong velocity component with an integrated flux in the peak channel of 4.7 ± 0.08 Jy is detected in the D-configuration image taken in August 2014. The flux declines to 0.17 ± 0.01 Jy in the B-configuration image taken in May 2015. The maser is undetected in Effelsberg observations taken in August 2019 with a 3σ upper limit of 0.36 Jy beam⁻¹. Such variability in the 6.7 GHz emission may explain why the maser has not been reported in previous surveys of methanol masers.
G78.9690+0.5410/IRAS 20286+3959. This maser was recently detected by Yang et al. (2019), with a peak velocity and peak flux very consistent with our VLA detection. The maser is associated with slightly extended emission at 3.6–24 µm (Fig. A.2). It has also been classified as a rising infrared SED protostar (Kryukova et al. 2014).

G78.9884+0.2211. This 6.7 GHz maser has not been reported before. The maser LSR velocity ($V_{\text{peak}} = -69$ km s⁻¹) is completely different from the velocities of the other masers detected in Cygnus X, which are within −15 to 15 km s⁻¹; therefore, this maser is probably at a different distance. The maser has five velocity components, with the strongest component having an integrated flux density at the peak channel of 1.97 ± 0.06 Jy in the D-configuration image. This declines to 0.32 ± 0.02 Jy in the B-configuration image, while the Effelsberg spectrum shows only one velocity component with an integrated flux of ~0.6 Jy. The second strongest component does not change in flux between the two VLA epochs, while the other three velocity components are not seen in the B-configuration map nor in the Effelsberg data.

We inspected CO (J = 1–0) data<fnref n="7" /> taken with the 45 m Nobeyama telescope by Yamagishi et al. (2018). The spectrum extracted from the maser position shows two peaks, one at $V_{\text{LSR}} = 9.2$ km s⁻¹ and one at $V_{\text{LSR}} = -62.8$ km s⁻¹, the latter of which is close to the methanol maser velocity. The emission at this velocity is distributed across a region of ≈2', while emission around $V_{\text{LSR}} = +9.2$ km s⁻¹ is more diffuse and extended. We calculated the distance to the $V_{\text{LSR}} = -62.8$ km s⁻¹ component using the Bayesian distance estimator developed by Reid et al. (2016)<fnref n="8" />. The estimator yields a distance of 9.35 ± 0.53 kpc with an integrated probability of 100% association with the Outer Arm. At a similar Galactic longitude there are two maser stars in the Outer Arm with parallaxes that are consistent with this distance (G075.29+01.32 and G073.65+00.19; Sanna et al. 2011; Reid et al. 2019). Certainly, the CO-traced $V_{\text{LSR}} = -62.8$ km s⁻¹ cloud associated with the methanol maser and likely also the $V_{\text{LSR}} = +9.2$ km s⁻¹ cloud are not part of the Cygnus X region, but are rather unrelated clouds lying on the same line of sight by chance.

Assuming that G78.9884+0.2211 is at 9.35 kpc, the mass, FIR luminosity, and bolometric luminosity must be rescaled to 347 $M_{\odot}$, 3093 $L_{\odot}$, and, 6237 $L_{\odot}$, respectively. Given the uncertainty in the distance, we do not include this maser source in the analysis presented in Sect. 5.

Wendker (1970) reported on radio continuum emission at 2.695 GHz, but at >8" from the maser location. There are several bright infrared peaks seen at 3.6–24 µm in the vicinity of the maser surrounded by extended emission (Fig. A.2). The maser position coincides with one infrared peak. This infrared source also shows a rising SED (Kryukova et al. 2014).

G79.7358+0.9904/IRAS 20290+4052. This maser is associated with bright, somewhat extended 3.6–24 µm emission (Fig. A.2). Kryukova et al. (2014) classified the infrared source as a protostar candidate with a rising SED. The infrared source was also recognized early on as a UC HII region based on their FIR colors (Bronfman et al. 1996; van der Walt et al. 1996), which meet the criteria for UC HII regions proposed by Wood & Churchwell (1989). However, as already pointed out, no radio continuum emission has been detected toward this source, which suggests that this object may be a precursor of a UC HII region.

G81.7523+0.5908/W75S-FIR2. The infrared environment toward this maser shows bright and slightly extended 3.6–8 µm emission (Fig. A.2). The source is saturated at 24 µm. The maser position coincides with the location of object W75S FIR2 (Harvey et al. 1986), a high-luminosity infrared protostellar core (Motte et al. 2007).

G81.7655+0.5972/CygX-N53. The emission seen in the D-array spectrum close to 7 km s⁻¹ is from sidelobes by the strong nearby maser G81.8713+0.7807/W75N(B). The 3.6–8 µm emission seen at the location of the maser (Fig. A.2) is remarkably weak (measured flux density at the location of the maser is 1–11 mJy beam⁻¹, respectively). The associated molecular core has been recognized as a massive infrared-quiet protostellar core (Motte et al. 2007). Observations at 1.3 mm with the IRAM Plateau de Bure interferometer revealed that the core is broken up into four smaller fragments (Bontemps et al. 2010). The maser position coincides with the strongest fragment, CygX-N53 MM1, which is a self-gravitating protostellar object and a good candidate to be a precursor of an OB star (Bontemps et al. 2010). Among the sample of 17 massive infrared-quiet protostellar cores identified by Motte et al. (2007), CygX-N53 and DR21(OH) are the only ones that show maser emission. However, the infrared emission from CygX-N53 is about an order of magnitude weaker than that of DR21(OH). The sample of Motte et al. (2007) comprises other infrared-quiet protostars ($S < 10$ Jy at 21 µm) with similar properties (mass, size, density) as CygX-N53, also located in the DR21 ridge; thus, it is unclear which conditions favor the onset of maser emission in CygX-N53 and not in the other cores that have similar physical properties.

G82.3079+0.7296. This maser was recently detected by Yang et al. (2019). They reported a peak flux of 58.4 Jy, which is more than two times higher than the flux we detected with the VLA. The maser is associated with extended 3.6–24 µm emission (Fig. A.2). The infrared source has been classified as a rising SED protostar (Kryukova et al. 2014).

### 4.5. Methanol absorption toward DR21

The region surrounding the prominent thermal radio source DR21 (Downes & Rinehart 1966) hosts a compact (about 30" × 30") and massive molecular core, which is located near the southern end of a long molecular cloud ridge oriented in the north-south direction and referred to as the DR21 ridge (see, e.g., Reipurth & Schneider 2008; Hennemann et al. 2012). This core harbors a group of five compact HII regions that are surrounded by a diffuse halo of ionized gas with a cometary morphology, with the most prominent compact component, DR21 M, at the head of the “comet” (Harris 1973; Roelfsema et al. 1989; Cyganowski et al. 2003).

The DR21 M HII region itself is surrounded by a dense molecular envelope with a photodissociation region (PDR) as an interface (Ossenkopf et al. 2010). The more compact and denser envelopes of UC HII regions often host Class II CH₃OH and OH masers, with W3(OH) the archetypal example (Menten et al. 1992). In contrast, the more developed envelope of DR21 M does not host any masers, but exhibits absorption features of OH in the hyperfine transitions of the rotational ground and excited states against the strong continuum emission of the compact HII region (see Jones et al. 1994, and references therein). As shown in Fig. 8, we detect absorption in the 6.7 GHz methanol line with a shape similar to that of the 1667 MHz OH line (taken from Koley et al. 2021). While the OH line extends to much

<footnote n="7">⁷ https://cygnus45.github.io/</footnote>

<footnote n="8">⁸ http://bessel vlbi astrometry org / bayesian</footnote>
<figure><img src="image_10.png" /><figcaption>**Fig. 8.** Absorption in the 6.7 GHz methanol line (lower panel) and the 1667 MHz OH line (upper panel) observed toward DR21 M. The dashed line indicates the systemic LSR velocity.</figcaption></figure>

lower velocities than that of CH₃OH, this may be a result of the
poorer S/N of the latter. Both lines show their deepest absorp-
tion at about VLSR = -1 km s⁻¹, close to the systemic velocity
of DR21 M, ≈2–3 km s⁻¹. Figure 9 shows that the methanol
absorption is seen only against the strongest continuum from
DR21 M, although this may also be caused by our limited S/N
toward regions with weaker extended continuum emission.

We note that the OH transition, but not the CH₃OH tran-
sition, also shows absorption at velocities between +6 and
+12 km s⁻¹. This represents low density material from a sec-
ond extended molecular cloud, sometimes called the “W75 N
9 km s⁻¹ cloud”, that covers the part of the W75 region that is
kinematically distinct from the cloud, which is also extended,
at the LSR velocity of DR21 (Dickel et al. 1978). Emission or
absorption from this cloud is only observed from transitions with
a low critical density, that is, ground-state or near ground-state
lines (see Koley et al. 2021).

We searched the D-configuration data cubes for other sources of methanol absorption, specifically toward the DR21 ridge, but did not detect absorption at any other location. In particular, we did not find absorption toward the Class I methanol maser sources in the vicinity of DR21(OH) and DR 21 (see Fig. 10). As explained by Menten (1991b) and Leurini et al. (2016), the same (purely collisional) pumping process that inverts the energy levels of the $J_{K=0} - (J-1)_{K=1} A^+$ Class I maser lines of A-type CH$_3$OH ($J \ge 7$) in molecular outflows also “anti-inverts” those of the 6.7 GHz $5_1-6_0 A^+$ line, causing enhanced absorption (overcooling) over more extended regions; this has actually been observed with single dish telescopes (Menten 1991b)$^9$. Apparently, in our data of the DR21 region, absorption signals confined to compact Class I maser spots are not strong enough to be

<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Methanol absorption (magenta contours, integrated from $V_{LSR} = -13.9$ to 0.5 km s⁻¹) seen in D-array data, overlaid on extended continuum emission toward DR21 (blue scale and contours) as seen in our VLA D+B array configuration map. The *n*th contour is at $(\sqrt{2})^n \times S_{max} \times p$, where: $S_{max} = 1.2$ Jy beam⁻¹, *n* = 1, 3, 5 ..., and *p* is equal to 0.4% for the radio continuum; and *n* = 0, 1, 2 ..., $S_{max} = -3.2$ Jy beam⁻¹ km s⁻¹, and *p* = 16% for the methanol absorption. The synthesized beams (4'' × 4'' for the VLA D+B map and 14''×4×10'' for PA = -84.1° for the D-configuration map) are shown at the bottom corner. The white stars mark the positions of Class I methanol masers observed in the 95 GHz $8_1-7_0 A^+$ line (Plambeck & Menten 1990). The orange stars are 44 GHz $7_0-6_1 A^+$ methanol masers (Kurtz et al. 2004).</figcaption></figure>

detectable at our sensitivity. In Appendix D we discuss con-
straints on the methanol abundance in the DR21 M molecular
core derived from our absorption spectrum.

**5. Discussion**

5.1. *Association with radio continuum emission*

One important aspect to establish an evolutionary sequence
model of methanol masers in star-forming regions is the coin-
cidence of masers and UC HII regions. Walsh et al. (1998)
proposed that methanol masers are present before the forma-
tion of a UC HII region around a massive star and persist until
the destruction of methanol as the UC HII region evolves. Our
high angular resolution observations can probe, for the first
time, the coincidence of masers and UC HII regions (by means
of their radio continuum emission) in the Cygnus X complex
down to scales of ~1700 au. While the GLOSTAR observations
are not ultra deep (our sensitivity limit is 0.06 mJy beam⁻¹),
they are more sensitive than most previous efforts and, notably,
have resulted in the detection of weak (a few mJy) continuum
emission toward five of our maser sources.

We found that five masers without a radio continuum coun-
terpart are associated with infrared sources that show a rising
SED between 4.5 and 24 µm (Kryukova et al. 2014), sug-
gesting that they host high-mass protostellar objects. Since

<footnote n="9">⁹ Analogously, for *E*-type CH₃OH, the *J*<sub>*k*=-1</sub> - (*J* - 1)><sub>*k*=0</sub> *E* lines (*J* ≥ 4) are masing, whereas the 12.1 GHz 2<sub>0</sub>-3<sub>-1</sub> *E* line is overcooled.</footnote>
<figure><img src="image_12.png" /></figure>

**Fig. 10.** Locations of 6.7 GHz Class II methanol masers detected with the VLA (blue circles) toward the dense molecular ridge extending from DR 21 M in the south to DR 21 (OH) (and beyond), overlaid on a 70 µm dust continuum image obtained with Herschel (red scale; Hennemann et al. 2012). Gray contours represent the GLOSTAR combined B+D-array continuum emission. The magenta contours show the velocity-integrated absorption of the 6.7 GHz methanol line. The nth contour is at ($\sqrt{2}$)ⁿ × $S_{max}$ × p, where: $S_{max}$ = 1.2 Jy beam⁻¹, n = 1, 3, 5 ..., and p is equal to 0.4% for the continuum; and n=0, 1, 2 ..., $S_{max}$ = -3.2 Jy beam⁻¹ km s⁻¹, and p = 16% for the methanol absorption. The white stars and green Xs mark, respectively, the positions of Class I methanol masers observed in the 95 GHz 8₁–7₀ A⁺ line (Plambeck & Menten 1990) and the 18 cm OH masers associated with DR21(OH) and W75S-FIR 1 (Norris et al. 1982; Argon et al. 2000). The orange stars represent 44 GHz Class I methanol masers (Kurtz et al. 2004).

they do not have associated radio continuum, these sources
may be in a younger evolutionary state compared to the
UC HII region (pre-UC HII). Three masers are associated with
known UC HII regions identified previously in the literature
(AFGL2591, Motte et al. 2007; IRAS 20290+4052, Bronfman
et al. 1996; and W75N(B), Hunter et al. 1994). However, we
found that G79.7358+0.9904/IRAS 20290+4052 does not show
radio continuum emission at a 3σ level of 0.22 mJy; thus, it
could be in a pre-UC HII state. DR21 is also associated with
a UC HII region; however, it does not have maser emission but
rather absorption in the methanol line. Three other masers are
found to be coincident (within ~1") with sources of compact
radio continuum emission; therefore, they are likely UC HII or
HC HII regions. We have indicated in Col. 4 of Table 6 whether
the masers are associated with a pre-UC HII or a UC or HC HII
region. In total, five dust cores host both a UC or HC HII region
and maser emission<fnref n="10" />, and all of them have $L_{Bol} > 1000 L_\odot$. This

<figure><img src="image_13.png" /></figure>

**Fig. 11.** Maser integrated luminosity against continuum flux density for the sample of masers detected by GLOSTAR in B configuration. The arrows indicate upper limits of continuum flux density.

implies an association rate between masers and HII regions of
38 ± 17%, which is consistent within the error with that found in
other studies (Walsh et al. 1998, 2003; Beuther et al. 2002; Hill
et al. 2005; Hu et al. 2016; Billington et al. 2019). However, we
should note here that our sample is much smaller than the sam-
ples studied in these previous works, which is reflected by the
large uncertainty in the association rate. As discussed in Sect. 1,
many of the MYSOs that power a methanol maser could be in
the HC HII stage. Because of their small angular size (≤0.03 pc;
Kurtz 2005), such sources are expected to have weak radio emis-
sion (~20 mJy at 1.4 kpc; Thompson et al. 2016) that, because of
their high electron density (≥10⁶ cm⁻³; Kurtz 2005), would be
optically thick up to high radio frequency, implying a spectral
index of ≈2. Unfortunately, the low flux levels of the contin-
uum sources (which causes these sources to remain undetected
at the individual subbands across the bandwidth; cf. Sect. 2)
and the very limited frequency range covered (of only 2 GHz)
preclude a meaningful determination of their spectral indices,
which would have allowed us to unambiguously classify them as
HC HII regions.

In Fig. 11 we show the 6.7 GHz maser integrated luminosity against continuum flux density for the five methanol masers that are associated with compact radio continuum emission as well as the 3σ upper limits for the non-detections. We find no significant correlation between these values (a Spearman correlation test yields a coefficient of 0.67 with a p-value = 0.22). If confirmed with more significant statistics, for instance by the GLOSTAR observations from the full Galactic plane, this would suggest that the 6.7 GHz maser phenomenon is not directly linked to the ionizing radiation.

5.2. *Mass and luminosity correlation*

We estimated the properties of the dust cores that are associ-
ated with a methanol source (Sect. 4.2). In order to study the
relationship between methanol masers and massive star-forming
dust cores, in this section we investigate the physical proper-
ties of dust cores with and without associated masers as well
as the correlations between the methanol masers and dust core
properties.

The upper panel of Fig. 12 shows the distribution of bolo-
metric luminosity for the sample of Cao et al. (2019) and the

<footnote n="10">¹⁰ We found that the number of expected background sources inside our matching radius is well below one source (see Appendix E). This implies a very low probability for the association between masers and continuum on the line of sight being just by chance.</footnote>
<table><caption>**Table 6.** Properties of the detected methanol line sources.</caption><thead><tr><th>ID</th><th>GLOSTAR name</th><th>SIMBAD or Common name</th><th>Possible nature</th><th>1.2 mm core<sup>(b)</sup></th><th>6-GHz continuum</th><th>V<sub>LSR</sub><sup>(a)</sup></th><th>V<sub>LSR</sub> Ref.</th></tr><tr><th>(1)</th><th>(2)</th><th>(3)</th><th>(4)</th><th>(5)</th><th>(6)</th><th>(7)</th><th>(8)</th></tr></thead><tbody><tr><td>1</td><td>G76.0932+0.1580</td><td>J202323.73+373535.16<sup>(c)</sup></td><td>pre-UC HII</td><td>—</td><td>—</td><td>0.7</td><td>1</td></tr><tr><td>2</td><td>G76.8437+0.1233</td><td>J202543.79+381112.98<sup>(c)</sup></td><td>pre-UC HII</td><td>—</td><td>—</td><td>-4.4</td><td>1</td></tr><tr><td>3</td><td>G78.8870+0.7087</td><td>AFGL2591</td><td>UC HII region</td><td>S26</td><td>yes</td><td>-5.7</td><td>2</td></tr><tr><td>4</td><td>G78.9690+0.5410</td><td>IRAS 20286+3959</td><td>pre-UC HII</td><td>—</td><td>—</td><td>7.0</td><td>3</td></tr><tr><td>5</td><td>G78.9884+0.2211</td><td>J203147.25+395900.33<sup>(c)</sup></td><td>pre-UC HII</td><td>—</td><td>—</td><td>-62.8</td><td>4</td></tr><tr><td>6</td><td>G79.7358+0.9904</td><td>IRAS 20290+4052</td><td>pre-UC HII</td><td>—</td><td>—</td><td>-1.4</td><td>5</td></tr><tr><td>7</td><td>G80.8617+0.3834</td><td>DR20, IRAS 20352+4124</td><td>UC or HC HII region</td><td>N14</td><td>yes</td><td>-1.5</td><td>6</td></tr><tr><td>8<sup>(d)</sup></td><td>G81.6790+0.5378</td><td>DR21</td><td>UC HII region</td><td>N46</td><td>yes</td><td>-3.0</td><td>7</td></tr><tr><td>9</td><td>G81.7219+0.5711</td><td>DR21(OH)</td><td>UC HII region</td><td>N44</td><td>yes</td><td>-2.6</td><td>3</td></tr><tr><td>10</td><td>G81.7444+0.5910</td><td>W75S-FIR1, DR21 B</td><td>UC or HC HII region</td><td>N43</td><td>yes</td><td>-3.1</td><td>7</td></tr><tr><td>11</td><td>G81.7523+0.5908</td><td>W75S-FIR2, DR21 A</td><td>High-luminosity IR protostellar core</td><td>N51</td><td>—</td><td>-3.3</td><td>7</td></tr><tr><td>12</td><td>G81.7655+0.5972</td><td>CygX-N53</td><td>Massive IR-quiet protostellar core</td><td>N53</td><td>—</td><td>-4.4</td><td>8</td></tr><tr><td>13</td><td>G81.8713+0.7807</td><td>W75N(B)</td><td>UC HII region</td><td>N30</td><td>yes</td><td>12.0</td><td>9</td></tr><tr><td>14</td><td>G82.3079+0.7296</td><td>J204016.75+425629.49<sup>(c)</sup></td><td>pre-UC HII</td><td>—</td><td>—</td><td>11.4</td><td>1</td></tr></tbody></table>

**Notes.** (a) Systemic LSR velocity of molecular gas. (b) From Motte et al. (2007). (c) From Kryukova et al. (2014). (d) Seen in absorption.

**References.** 1. Schneider et al. (2006); 2. van der Tak et al. (1999); 3. Keown et al. (2019); 4. Yamagishi et al. (2018); 5. Bronfman et al. (1996); 6. Pipenbrink & Wendker (1988); 7. This work; 8. Duarte-Cabral et al. (2014); 9. Gottschalk et al. (2012).

maser associated cores. Although the sample of cores with associated masers is small, it does seem to cover a narrower range in luminosity than the range occupied by the full sample of massive dense cores; however, it still covers a range of a few orders of magnitude (from ~$10^{1.4}$ to ~$10^5 L_\odot$). A similar result was reported by Billington et al. (2019) for the ATLASGAL clumps. They found that maser associated clumps have, on average, higher luminosities when compared to their full sample. Figure 12 also suggests a luminosity limit for maser emission of ~200 $L_\odot$, which is lower than the limit of $10^3 L_\odot$ estimated by Bourke et al. (2005) and more similar to the limit found in Paulson & Pandian (2020) and Jones et al. (2020). The bottom panel of Fig. 12 shows the distribution of the bolometric luminosity-to-core-mass ratio, $L_{\text{Bol}}/M_{\text{core}}$, for the sample of Cao et al. (2019) and the maser associated cores found with GLOSTAR. There seems to be a threshold of this parameter for maser emission since the minimum $L_{\text{Bol}}/M_{\text{core}}$ of a source that is associated with maser emission is ~1 $L_\odot M_\odot^{-1}$. A similar value was suggested by Billington et al. (2019).

Figure 13 shows the correlations between methanol masers and dust core properties. The left and middle panels present maser integrated luminosity as a function of core mass and bolometric luminosity, respectively. We see in these figures that if we exclude the maser with the strongest luminosity, the correlation between these properties is nearly nonexistent. However, this may be due to the limited size of our sample. Billington et al. (2019) found a weak correlation between clump bolometric luminosity and maser luminosity for a much larger sample (958 dust clumps that host methanol masers). In another study, Paulson & Pandian (2020) also found a weak correlation between clump and maser luminosity for their sample of 320 MMB masers. These authors also noted that the relationship between clump and maser luminosity for the low-luminosity ($L < 10^{-6} L_\odot$) maser population is consistent with no correlation. If this non-correlation is confirmed for a larger sample, it would suggest that, for low-luminosity masers, the maser luminosity is probably not mainly driven by core luminosity but by other factors, as suggested by Paulson & Pandian (2020), such as gas density, gas temperature,

<figure><img src="image_14.png" /><figcaption>**Fig. 12.** Distribution of core bolometric luminosities (top) and the $L_{\text{Bol}}/M_{\text{core}}$ ratio (bottom). In both panels the full sample of dust cores by Cao et al. (2019) is shown in blue. The distribution of cores with associated masers detected by GLOSTAR is shown in red.</figcaption></figure>

and methanol fractional abundance, which are not explored in this paper. We note that the two new detected masers are associated with dust cores that have the lowest mass, $M_{\text{core}}$, and also have the lowest luminosity, $L_{\text{maser}, B}$ (cf. Table 4). While previous maser studies toward Cygnus X have selected the brightest
<figure><img src="image_15.png" /><figcaption>**Fig. 13.** Maser integrated luminosity measured in B-configuration GLOSTAR data as a function of core mass (**left**) and bolometric luminosity (**middle**). Right panel: maser integrated luminosity against the $L_{\text{Bol}}/M_{\text{core}}$ ratio.</figcaption></figure>

sources (e.g., Hu et al. 2016; Yang et al. 2019), GLOSTAR is an unbiased survey and is also more sensitive than most of these previous efforts, allowing the detection of more new masers at the lower end of the maser luminosity distribution function.

The luminosity-to-mass ($L/M$) ratio has been proposed to be an evolutionary indicator during the formation process of massive clumps (e.g., Molinari et al. 2008, see, however, Ma et al. 2013). The $L/M$ ratio increases with age as a result of the envelope mass falling into the forming star. Previous works have found a weak correlation between the $L/M$ ratio of maser associated clumps and maser luminosity (e.g., Billington et al. 2019), suggesting that higher-luminosity masers are associated with more evolved clumps. We investigated this relation for the maser host cores in Cygnus X. The right panel of Fig. 13 shows maser integrated luminosity against the $L_{\text{Bol}}/M_{\text{core}}$ ratio. The correlation between these quantities is not significant (the Spearman correlation coefficient is 0.3 with $p$-value = 0.3), which can be attributed to the large scatter of the measurements and the small size of our sample. However, we note that the masers associated with UC or HC HII regions (cf. Table 6) occupy the region in this diagram with high $L_{\text{Bol}}/M_{\text{core}}$ ratios, as expected for more evolved objects; specifically, they all lie above $L_{\text{Bol}}/M_{\text{core}} \approx 8 L_{\odot} M_{\odot}^{-1}$.

In summary, the results above suggest that there is a core bolometric luminosity threshold for the maser emission of ~200 $L_{\odot}$, which is in agreement with that found in recent works that have studied a significant sample of maser host clumps in the Galactic plane but excluded the Cygnus X complex. For the minimum core mass, we found a value of ≈10 $M_{\odot}$ in our sample (cf. Table 4), which is likely sufficient to form a single high-mass star. In a previous study, Billington et al. (2019) reported a minimum FWHM mass (mass within 50% of the 850 or 870 µm contour) of 17 $M_{\odot}$ for their distance-limited sample of cores that host methanol masers, while the minimum mass reported by Paulson & Pandian (2020) is 11 $M_{\odot}$. The distribution of the luminosity-to-mass ratio of dust cores that host maser emission is narrower than the distribution of the full core population, with methanol host cores having ratios between ~1 and $10^2 L_{\odot} M_{\odot}^{-1}$. This is consistent with that expected from the evolutionary sequence model of Breen et al. (2010) and Billington et al. (2019), who suggest that the maser emission starts at a $L/M$ ratio of $10^{0.6} L_{\odot} M_{\odot}^{-1}$, followed by the maser emission and continuum emission from the HII region being simultaneously detectable at higher values of $L/M$, leading to the final decline of the maser emission due to the disruptive effects

of the HII region at a $L/M$ value of $10^{2.2} L_{\odot} M_{\odot}^{-1}$. The threshold of the luminosity-to-mass ratio of ~1 $L_{\odot} M_{\odot}^{-1}$ is also consistent with what has been measured for a larger sample of star-forming regions in the Galactic plane (Billington et al. 2019, 2020).

## 6. Conclusions

We have conducted the first unbiased survey of Class II methanol masers at 6.7 GHz in the Cygnus X complex. A total of 13 masers were detected in the observed 7° × 3° area, for which we derive fluxes and integrated luminosities, as well as absolute positions of the maser features. In addition, methanol absorption is detected toward DR21 against the compact continuum from the HII regions inside the DR21 core. We have established the associations between maser sources and radio continuum, (sub)millimeter, and infrared emission to within ≈1". All masers are associated with dust continuum cores, and all of them are close to or coincident with near-infrared or FIR peaks, as expected from the maser-MYSO relationship. Only a fraction of them (38 ± 17%) have radio continuum detected by our survey or in previous observations. These sources are more likely to be UC HII or HC HII regions.

By fitting the SED to the dust continuum emission, we have derived the properties of the cores that host maser emission and compared their distributions with those of the full population of dust cores in the Cygnus X complex identified in the literature. We found thresholds for maser emission in luminosity and the luminosity-to-mass ratio of ~200 $L_{\odot}$ and ~1 $L_{\odot} M_{\odot}^{-1}$, respectively. These values are in agreement with recent works that have investigated a larger sample of methanol maser hosts in the Galactic plane.

*Acknowledgements.* The authors are grateful to the anonymous referee, whose comments helped to improve this paper. G.N.O.-L. would like to thank Antonio Hernandez-Gómez for a detailed reading of an early version of the manuscript and valuable suggestions. We thank Yue Cao and Keping Qiu for providing mosaics maps of (sub)millimeter emission. We also thank Nicola Schneider for sharing her IRAM CO data. G.N.O.-L. acknowledges support from the von Humboldt Stiftung. T.Cs. has received financial support from the French State in the framework of the IdEx Université de Bordeaux Investments for the future Program. H.B. acknowledges support from the European Research Council under the European Community's Horizon 2020 framework program (2014–2020) via the ERC Consolidator grant “From Cloud to Star Formation (CSF)” (project number 648505). H.B. further acknowledges support from the Deutsche Forschungsgemeinschaft (DFG) via Sonderforschungsbereich (SFB) 881 “The Milky Way System” (sub-project B1). The National Radio Astronomy Observatory (NRAO) is operated by Associated Universities Inc., under cooperative agreement with the National Science Foundation.
## References

Aguirre, J. E., Ginsburg, A. G., Dunham, M. K., et al. 2011, ApJS, 192, 4

Anglada, G., Villuendas, E., Estalella, R., et al. 1998, AJ, 116, 2953

Araya, E. D., Kurtz, S., Hofner, P., & Linz, H. 2009, ApJ, 698, 1321

Argon, A. L., Reid, M. J., & Menten, K. M. 2000, ApJS, 129, 159

Batrla, W., Matthews, H. E., Menten, K. M., & Walmsley, C. M. 1987, Nature, 326, 49

Beerer, I. M., Koenig, X. P., Hora, J. L., et al. 2010, ApJ, 720, 679

Beuther, H., Walsh, A., Schilke, P., et al. 2002, A&A, 390, 289

Billington, S. J., Urquhart, J. S., König, C., et al. 2019, MNRAS, 490, 2779

Billington, S. J., Urquhart, J. S., König, C., et al. 2020, MNRAS, 499, 2744

Bonfand, M., Belloche, A., Menten, K. M., Garrod, R. T., & Müller, H. S. P. 2017, A&A, 604, A60

Bontemps, S., Motte, F., Csengeri, T., & Schneider, N. 2010, A&A, 524, A18

Bourke, T. L., Hyland, A. R., & Robinson, G. 2005, ApJ, 625, 883

Breckenridge, S. M., & Kukolich, S. G. 1995, ApJ, 438, 504

Breen, S. L., Ellingsen, S. P., Caswell, J. L., & Lewis, B. E. 2010, MNRAS, 401, 2219

Breen, S. L., Ellingsen, S. P., Caswell, J. L., et al. 2011, ApJ, 733, 80

Bronfman, L., Nyman, L. A., & May, J. 1996, A&AS, 115, 81

Brunthaler, A., Menten, K. M., Dzib, S. A., et al. 2021, A&A, 651, A85 (Paper I)

Campbell, B. 1984, ApJ, 287, 334

Cao, Y., Qiu, K., Zhang, Q., et al. 2019, ApJS, 241, 1

Caswell, J. L., Vaile, R. A., Ellingsen, S. P., Whiteoak, J. B., & Norris, R. P. 1995, MNRAS, 272, 96

Caswell, J. L., Fuller, G. A., Green, J. A., et al. 2010, MNRAS, 404, 1029

Caswell, J. L., Fuller, G. A., Green, J. A., et al. 2011, MNRAS, 417, 1964

Chandler, C. J., Gear, W. K., & Chini, R. 1993, MNRAS, 260, 337

Cotton, W. D. 2008, PASP, 120, 439

Cragg, D. M., Sobolev, A. M., & Godfrey, P. D. 2005, MNRAS, 360, 533

Csengeri, T., Urquhart, J. S., Schuller, F., et al. 2014, A&A, 565, A75

Cyganowski, C. J., Reid, M. J., Fish, V. L., & Ho, P. T. P. 2003, ApJ, 596, 344

Dickel, J. R., Dickel, H. R., & Wilson, W. J. 1978, ApJ, 223, 840

Downes, D., & Rinehart, R. 1966, ApJ, 144, 937

Duarte-Cabral, A., Bontemps, S., Motte, F., et al. 2014, A&A, 570, A1

Dunham, M. M., Crapsi, A., II, N. J. E., et al. 2008, ApJS, 179, 249

Dzib, S. A., Rodríguez, L. F., Loinard, L., et al. 2013, ApJ, 763, 139

Eden, D. J., Moore, T. J. T., Plume, R., et al. 2017, MNRAS, 469, 2163

Ellingsen, S. P. 2006, ApJ, 638, 241

Gieser, C., Semenov, D., Beuther, H., et al. 2019, A&A, 631, A142

Gottschalk, M., Kothes, R., Matthews, H. E., Landecker, T. L., & Dent, W. R. F. 2012, A&A, 541, A79

Green, J. A., Caswell, J. L., Fuller, G. A., et al. 2009, MNRAS, 392, 783

Green, J. A., Caswell, J. L., Fuller, G. A., et al. 2010, MNRAS, 409, 913

Green, J. A., Caswell, J. L., Fuller, G. A., et al. 2012, MNRAS, 420, 3108

Green, J. A., Breen, S. L., Fuller, G. A., et al. 2017, MNRAS, 469, 1383

Hales, C. A., Murphy, T. J. T., Curran, J. R., et al. 2012, MNRAS, 425, 979

Harris, S. 1973, MNRAS, 162, 5P

Harvey, P. M., Joy, M., Lester, D. F., & Wilking, B. A. 1986, ApJ, 300, 737

Hennemann, M., Motte, F., Schneider, N., et al. 2012, A&A, 543, L3

Hennemann, M., Motte, F., & Schneider, N. 2014, in The Labyrinth of Star Formation (Berlin: Springer), 36, 271

Hill, T., Burton, M. G., Minier, V., et al. 2005, MNRAS, 363, 405

Hu, B., Menten, K. M., Wu, Y., et al. 2016, ApJ, 833, 18

Hunter, T. R., Taylor, G. B., Felli, M., & Tofani, G. 1994, A&A, 284, 215

Johnston, K. G., Shepherd, D. S., Robitaille, T. P., & Wood, K. 2013, A&A, 551, A43

Jones, K. N., Field, D., Gray, M. D., & Walker, R. N. F. 1994, A&A, 288, 581

Jones, B. M., Fuller, G. A., Breen, S. L., et al. 2020, MNRAS, 493, 2015

Kauffmann, J., Bertoldi, F., Bourke, T. L., Evans, II, N. J., & Lee, C. W. 2008, A&A, 487, 993

Keown, J., Di Francesco, J., Rosolowsky, E., et al. 2019, ApJ, 884, 4

Koley, A., Roy, N., Menten, K. M., et al. 2021, MNRAS, 501, 4825

König, C., Urquhart, J. S., Csengeri, T., et al. 2017, A&A, 599, A139

König, C., Urquhart, J. S., Wyrowski, F., Colombo, D., & Menten, K. M. 2021, A&A, 645, A113

Kraemer, K. E., Hora, J. L., Egan, M. P., et al. 2010, AJ, 139, 2319

Kryukova, E., Megeath, S. T., Hora, J. L., et al. 2014, AJ, 148, 11

Kurtz, S. 2005, IAU Symp., 227, 111

Kurtz, S., Hofner, P., & Álvarez, C. V. 2004, ApJS, 155, 149

Leurini, S., Menten, K. M., & Walmsley, C. M. 2016, A&A, 592, A31

Ma, B., Tan, J. C., & Barnes, P. J. 2013, ApJ, 779, 79

Medina, S. N. X., Urquhart, J. S., Dzib, S. A., et al. 2019, A&A, 627, A175

Menten, K. M. 1991a, ASP Conf. Ser., 16, 119

Menten, K. M. 1991b, ApJ, 380, L75

Menten, K. M., Walmsley, C. M., Henkel, C., & Wilson, T. L. 1988, A&A, 198, 253

Menten, K. M., Reid, M. J., Pratap, P., Moran, J. M., & Wilson, T. L. 1992, ApJ, 401, L39

Minier, V., Ellingsen, S. P., Norris, R. P., & Booth, R. S. 2003, A&A, 403, 1095

Molet, J., Brouillet, N., Nony, T., et al. 2019, A&A, 626, A132

Molinari, S., Pezzuto, S., Cesaroni, R., et al. 2008, A&A, 481, 345

Molinari, S., Swinyard, B., Bally, J., et al. 2010, PASP, 122, 314

Moore, T. J. T., Plume, R., Thompson, M. A., et al. 2015, MNRAS, 453, 4264

Motte, F., Bontemps, S., Schilke, P., et al. 2007, A&A, 476, 1243

Motte, F., Zavagno, A., Bontemps, S., et al. 2010, A&A, 518, L77

Müller, H. S. P., Schläder, F., Stutzki, J., & Winnewisser, G. 2005, J. Mol. Struc., 742, 215

Norris, R. P., Booth, R. S., Diamond, P. J., & Porter, N. D. 1982, MNRAS, 201, 191

Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943

Ossenkopf, V., Röllig, M., Simon, R., et al. 2010, A&A, 518, L79

Pandian, J. D., Leurini, S., Menten, K. M., Belloche, A., & Goldsmith, P. F. 2008, A&A, 489, 1175

Pandian, J. D., Momjian, E., Xu, Y., Menten, K. M., & Goldsmith, P. F. 2010, A&A, 522, A8

Paulson, S. T., & Pandian, J. D. 2020, MNRAS, 492, 1335

Pestalozzi, M. R., Minier, V., & Booth, R. S. 2005, A&A, 432, 737

Pipenbrink, A., & Wendker, H. J. 1988, A&A, 191, 313

Plambeck, R. L., & Menten, K. M. 1990, ApJ, 364, 555

Price, S. D., Egan, M. P., Carey, S. J., Mizuno, D. R., & Kuchar, T. A. 2001, AJ, 121, 2819

Purcell, C. R., Hoare, M. G., Cotton, W. D., et al. 2013, ApJS, 205, 1

Reid, M. J., Dame, T. M., Menten, K. M., & Brunthaler, A. 2016, ApJ, 823, 77

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2019, ApJ, 885, 131

Reipurth, B., & Schneider, N. 2008, ASP Monograph Pub., 4, 36

Rodríguez-Kamenetzky, A., Carrasco-González, C., Torrelles, J. M., et al. 2020, MNRAS, 496, 3128

Roelfsema, P. R., Goss, W. M., & Geballe, T. R. 1989, A&A, 222, 247

Rosolowsky, E., Dunham, M. K., Ginsburg, A., et al. 2010, ApJS, 188, 123

Roy, A., Ade, P. A. R., Bock, J. J., et al. 2011, ApJ, 727, 114

Rygl, K. L. J., Brunthaler, A., Sanna, A., et al. 2012, A&A, 539, A79

Sanna, A., Reid, M. J., Dame, T. M., et al. 2011, ApJ, 745, 82

Schneider, N., Bontemps, S., Simon, R., et al. 2006, A&A, 458, 855

Schuller, F., Menten, K. M., Contreras, Y., et al. 2009, A&A, 504, 415

Sobolev, A. M., & Deguchi, S. 1994, A&A, 291, 569

Thompson, M., Goedhart, S., Goedhart, S., et al. 2016, in MeerKAT Science: On the Pathway to the SKA, 15

Thompson, A. R., Moran, J. M., & Swenson, George W., J. 2017, Interferometry and Synthesis in Radio Astronomy, 3rd edn. (Hoboken: Wiley)

Torrelles, J. M., Gómez, J. F., Rodríguez, L. F., et al. 1997, ApJ, 489, 744

Urquhart, J. S., Moore, T. J. T., Schuller, F., et al. 2013, MNRAS, 431, 1752

Urquhart, J. S., Csengeri, T., Wyrowski, F., et al. 2014, A&A, 568, A41

Urquhart, J. S., Moore, T. J. T., Menten, K. M., et al. 2015, MNRAS, 446, 3461

Uyaniker, B., Fürst, E., Reich, W., Aschenbach, B., & Wielebinski, R. 2001, A&A, 371, 675

van der Tak, F. F. S., van Dishoeck, E. F., Evans, Neal J., I., Bakker, E. J., & Blake, G. A. 1999, ApJ, 522, 991

van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, A&A, 468, 627

van der Walt, D. J., Retief, S. J. P., Gaylard, M. J., & MacLeod, G. C. 1996, MNRAS, 282, 1085

Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, G. 1998, MNRAS, 301, 640

Walsh, A. J., Macdonald, G. H., Alvey, N. D. S., Burton, M. G., & Lee, J. K. 2003, A&A, 410, 597

Wendker, H. J. 1970, A&A, 4, 378

Wendker, H. J. 1984, A&AS, 58, 291

Wendker, H. J., Higgs, L. A., & Landecker, T. L. 1991, A&A, 241, 551

Wood, D. O. S., & Churchwell, E. 1989, ApJ, 340, 265

Wright, N. J., Drake, J. J., Drew, J. E., & Vink, J. S. 2010, ApJ, 713, 871

Xu, Y., Li, J. J., Hachisuka, K., et al. 2008, A&A, 485, 729

Yamagishi, M., Nishimura, A., Fujita, S., et al. 2018, ApJS, 235, 9

Yang, K., Chen, X., Shen, Z.-Q., et al. 2019, ApJS, 241, 18
## Appendix A: Supplementary figures and tables

In this appendix, Figs. A.1 and A.2 show (sub)millimeter dust emission and *Spitzer* infrared emission in a region of $3' \times 3'$ and $0.7' \times 0.7'$, respectively, around maser positions.

Table A.1 lists the catalog of maser spots detected above $4\sigma$ in B-configuration data.

Figure A.3 displays the distribution of all methanol maser spots measured in the B-configuration images. For five of our sources, the 6.7 GHz methanol maser emission was imaged by Rygl et al. (2012) with the EVN with a resolution of a few milliarcseconds (mas), much better than our $1''$. A comparison between the EVN positions, which have absolute position uncertainties better than 1 mas, and even our higher resolution B-array VLA data is not straightforward for several reasons:

generally, in the VLA images several maser spots remain spatially unresolved, and our positions represent averages. Moreover, for several sources the EVN data show maser emission only over a fraction of the velocity range of the VLA data, leaving the determination of an emission centroid biased. Nevertheless, we note that for most sources the EVN-determined maser spots lie within the distributions of the VLA B-array spots (within tens of mas). For two sources, G81.7219+0.5711 [DR21(OH)] and G81.8713+0.78-7 [W75N(B)], we find offsets between the EVN and VLA positions of $\approx 0.''2$ and $\approx 0.''15$, respectively. For DR21(OH), the significance of the offset is difficult to assess as many fewer spots (a total of three) appear in the EVN data than in the VLA data over a much narrower velocity range. Still, the offsets are smaller than or comparable to the position uncertainty discussed in Sect. 3.

<table><caption>**Table A.1.** Maser spot catalog for B-configuration data.</caption><thead><tr><th>Name<br>(1)</th><th>Ref. RA<br>(h:m:s)<br>(2)</th><th>Ref. Dec<br>(°:′:′′)<br>(3)</th><th>Δα<br>(arcsec)<br>(4)</th><th>Δδ<br>(arcsec)<br>(5)</th><th>V<sub>LSR</sub><br>(km s<sup>-1</sup>)<br>(6)</th><th>S<sub>v, Peak</sub><br>(Jy beam<sup>-1</sup>)<br>(7)</th><th>S<sub>v, Int.</sub><br>(Jy)<br>(8)</th></tr></thead><tbody><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>0.000</td><td>0.000</td><td>6.64</td><td>0.57 ± 0.01</td><td>0.59 ± 0.02</td></tr><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>-0.089</td><td>-0.024</td><td>6.46</td><td>0.38 ± 0.01</td><td>0.44 ± 0.02</td></tr><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>-0.103</td><td>0.055</td><td>4.84</td><td>0.38 ± 0.01</td><td>0.48 ± 0.01</td></tr><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>-0.281</td><td>-0.005</td><td>-6.50</td><td>0.17 ± 0.01</td><td>0.18 ± 0.01</td></tr><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>-0.067</td><td>-0.132</td><td>5.56</td><td>0.15 ± 0.01</td><td>0.12 ± 0.01</td></tr><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>-0.202</td><td>-0.126</td><td>-6.32</td><td>0.14 ± 0.01</td><td>0.12 ± 0.01</td></tr><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>0.015</td><td>-0.005</td><td>4.66</td><td>0.15 ± 0.00</td><td>0.18 ± 0.01</td></tr><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>-0.143</td><td>0.135</td><td>5.02</td><td>0.10 ± 0.01</td><td>0.12 ± 0.01</td></tr><tr><td>G76.0932+0.1580</td><td>20:23:23.7140</td><td>+37:35:35.558</td><td>-0.172</td><td>-0.183</td><td>5.74</td><td>0.11 ± 0.01</td><td>0.11 ± 0.01</td></tr><tr><td>G76.8437+0.1233</td><td>20:25:43.7657</td><td>+38:11:12.734</td><td>0.000</td><td>0.000</td><td>-5.42</td><td>0.17 ± 0.00</td><td>0.17 ± 0.01</td></tr><tr><td>G76.8437+0.1233</td><td>20:25:43.7657</td><td>+38:11:12.734</td><td>0.082</td><td>0.079</td><td>-5.24</td><td>0.15 ± 0.00</td><td>0.18 ± 0.01</td></tr><tr><td>G76.8437+0.1233</td><td>20:25:43.7657</td><td>+38:11:12.734</td><td>0.110</td><td>0.156</td><td>-5.60</td><td>0.07 ± 0.00</td><td>0.06 ± 0.01</td></tr><tr><td>G78.8870+0.7087</td><td>20:29:24.9436</td><td>+40:11:19.626</td><td>0.000</td><td>0.000</td><td>-7.04</td><td>0.10 ± 0.01</td><td>0.19 ± 0.03</td></tr><tr><td>G78.8870+0.7087</td><td>20:29:24.9436</td><td>+40:11:19.626</td><td>-0.054</td><td>-0.160</td><td>-7.22</td><td>0.09 ± 0.01</td><td>0.11 ± 0.01</td></tr><tr><td>G78.8870+0.7087</td><td>20:29:24.9436</td><td>+40:11:19.626</td><td>0.170</td><td>-0.404</td><td>-7.40</td><td>0.09 ± 0.01</td><td>0.09 ± 0.01</td></tr><tr><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>0.000</td><td>0.000</td><td>4.84</td><td>1.25 ± 0.01</td><td>1.28 ± 0.02</td></tr><tr><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>-0.008</td><td>-0.063</td><td>4.66</td><td>1.07 ± 0.00</td><td>1.08 ± 0.01</td></tr><tr><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>0.035</td><td>-0.014</td><td>4.48</td><td>0.95 ± 0.01</td><td>1.02 ± 0.01</td></tr><tr><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>0.063</td><td>-0.019</td><td>5.02</td><td>0.88 ± 0.01</td><td>0.94 ± 0.01</td></tr><tr><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>0.025</td><td>-0.131</td><td>4.30</td><td>0.57 ± 0.01</td><td>0.61 ± 0.01</td></tr><tr><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>0.021</td><td>0.038</td><td>5.20</td><td>0.21 ± 0.01</td><td>0.21 ± 0.01</td></tr><tr><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>-0.098</td><td>-0.019</td><td>4.12</td><td>0.12 ± 0.00</td><td>0.10 ± 0.00</td></tr><tr><td>G78.9690+0.5410</td><td>20:30:22.6762</td><td>+40:09:23.348</td><td>0.181</td><td>-0.099</td><td>5.38</td><td>0.10 ± 0.00</td><td>0.09 ± 0.01</td></tr><tr><td>G78.9884+0.2211</td><td>20:31:47.3082</td><td>+39:58:59.903</td><td>0.000</td><td>0.000</td><td>-68.96</td><td>0.25 ± 0.01</td><td>0.30 ± 0.01</td></tr><tr><td>G78.9884+0.2211</td><td>20:31:47.3082</td><td>+39:58:59.903</td><td>-0.224</td><td>0.096</td><td>-66.44</td><td>0.20 ± 0.00</td><td>0.16 ± 0.01</td></tr><tr><td>G78.9884+0.2211</td><td>20:31:47.3082</td><td>+39:58:59.903</td><td>-0.028</td><td>-0.103</td><td>-66.26</td><td>0.21 ± 0.01</td><td>0.30 ± 0.02</td></tr><tr><td>G78.9884+0.2211</td><td>20:31:47.3082</td><td>+39:58:59.903</td><td>-0.564</td><td>0.122</td><td>-68.78</td><td>0.22 ± 0.01</td><td>0.34 ± 0.03</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.000</td><td>0.000</td><td>-5.42</td><td>22.96 ± 0.21</td><td>23.81 ± 0.39</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.029</td><td>-0.001</td><td>-5.60</td><td>21.87 ± 0.24</td><td>22.88 ± 0.44</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.028</td><td>-0.009</td><td>-5.78</td><td>11.84 ± 0.16</td><td>12.70 ± 0.30</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.007</td><td>-0.016</td><td>-5.24</td><td>7.38 ± 0.09</td><td>7.85 ± 0.17</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.047</td><td>-0.001</td><td>-5.96</td><td>5.62 ± 0.08</td><td>5.79 ± 0.15</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.055</td><td>-0.010</td><td>-3.98</td><td>5.88 ± 0.09</td><td>6.11 ± 0.15</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.045</td><td>-0.039</td><td>-3.26</td><td>4.97 ± 0.06</td><td>5.07 ± 0.12</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.074</td><td>-0.022</td><td>-3.08</td><td>4.62 ± 0.09</td><td>4.81 ± 0.17</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.010</td><td>-0.015</td><td>-3.80</td><td>2.86 ± 0.06</td><td>2.88 ± 0.11</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.088</td><td>-0.026</td><td>-3.62</td><td>2.34 ± 0.05</td><td>2.41 ± 0.10</td></tr></tbody></table>

**Notes.** The position uncertainties are $\pm 0.''18$ and $\pm 0.''13$ in RA and Dec, respectively (see Sect. 3).
<table><caption>Table A.1. continued.</caption><thead><tr><th>Name<br>(1)</th><th>Ref. RA<br>(h:m:s)<br>(2)</th><th>Ref. Dec<br>(°:′′)<br>(3)</th><th>Δα<br>(arcsec)<br>(4)</th><th>Δδ<br>(arcsec)<br>(5)</th><th>V<sub>LSR</sub><br>(km s<sup>-1</sup>)<br>(6)</th><th>S<sub>v, Peak</sub><br>(Jy beam<sup>-1</sup>)<br>(7)</th><th>S<sub>v, Int.</sub><br>(Jy)<br>(8)</th></tr></thead><tbody><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.013</td><td>0.004</td><td>-6.14</td><td>2.26 ± 0.04</td><td>2.35 ± 0.08</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.180</td><td>0.021</td><td>-4.88</td><td>0.95 ± 0.04</td><td>1.03 ± 0.07</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.018</td><td>0.020</td><td>-6.32</td><td>0.90 ± 0.04</td><td>0.93 ± 0.07</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.062</td><td>0.002</td><td>-4.70</td><td>0.86 ± 0.03</td><td>0.82 ± 0.05</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>0.012</td><td>-0.054</td><td>-4.16</td><td>0.68 ± 0.06</td><td>0.74 ± 0.10</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.132</td><td>-0.111</td><td>-5.06</td><td>0.39 ± 0.05</td><td>0.55 ± 0.10</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.250</td><td>-0.031</td><td>-4.52</td><td>0.37 ± 0.02</td><td>0.33 ± 0.03</td></tr><tr><td>G79.7358+0.9904</td><td>20:30:50.6680</td><td>+41:02:27.403</td><td>-0.070</td><td>0.017</td><td>-4.34</td><td>0.14 ± 0.01</td><td>0.10 ± 0.01</td></tr><tr><td>G80.8617+0.3834</td><td>20:37:00.9571</td><td>+41:34:55.606</td><td>0.000</td><td>0.000</td><td>-3.98</td><td>4.61 ± 0.11</td><td>5.02 ± 0.21</td></tr><tr><td>G80.8617+0.3834</td><td>20:37:00.9571</td><td>+41:34:55.606</td><td>-0.003</td><td>0.001</td><td>-4.16</td><td>4.10 ± 0.12</td><td>4.55 ± 0.24</td></tr><tr><td>G80.8617+0.3834</td><td>20:37:00.9571</td><td>+41:34:55.606</td><td>-0.029</td><td>0.331</td><td>-2.00</td><td>1.56 ± 0.06</td><td>1.77 ± 0.12</td></tr><tr><td>G80.8617+0.3834</td><td>20:37:00.9571</td><td>+41:34:55.606</td><td>0.153</td><td>-0.035</td><td>-4.34</td><td>1.24 ± 0.07</td><td>1.42 ± 0.14</td></tr><tr><td>G80.8617+0.3834</td><td>20:37:00.9571</td><td>+41:34:55.606</td><td>-0.000</td><td>0.320</td><td>-2.18</td><td>0.89 ± 0.06</td><td>1.02 ± 0.11</td></tr><tr><td>G80.8617+0.3834</td><td>20:37:00.9571</td><td>+41:34:55.606</td><td>0.206</td><td>-0.093</td><td>-4.52</td><td>0.61 ± 0.06</td><td>0.81 ± 0.13</td></tr><tr><td>G80.8617+0.3834</td><td>20:37:00.9571</td><td>+41:34:55.606</td><td>0.603</td><td>0.044</td><td>-3.44</td><td>0.14 ± 0.04</td><td>0.61 ± 0.20</td></tr><tr><td>G80.8617+0.3834</td><td>20:37:00.9571</td><td>+41:34:55.606</td><td>-0.095</td><td>0.204</td><td>-3.62</td><td>0.23 ± 0.02</td><td>0.26 ± 0.04</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>0.000</td><td>0.000</td><td>-2.72</td><td>4.07 ± 0.05</td><td>4.30 ± 0.10</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>0.202</td><td>0.048</td><td>-3.08</td><td>2.72 ± 0.05</td><td>3.12 ± 0.10</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>0.216</td><td>0.044</td><td>-3.26</td><td>2.19 ± 0.04</td><td>2.49 ± 0.07</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>0.098</td><td>0.034</td><td>-2.90</td><td>2.36 ± 0.04</td><td>2.54 ± 0.07</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-0.325</td><td>-0.133</td><td>-3.80</td><td>1.30 ± 0.03</td><td>1.43 ± 0.05</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-0.023</td><td>0.025</td><td>-2.54</td><td>1.60 ± 0.04</td><td>1.79 ± 0.07</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-1.163</td><td>-0.226</td><td>8.98</td><td>1.22 ± 0.03</td><td>1.25 ± 0.06</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>0.187</td><td>0.073</td><td>-3.44</td><td>0.89 ± 0.02</td><td>1.03 ± 0.05</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-0.433</td><td>-0.126</td><td>-3.98</td><td>0.88 ± 0.02</td><td>0.87 ± 0.04</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-1.240</td><td>-0.223</td><td>8.80</td><td>0.68 ± 0.03</td><td>0.70 ± 0.05</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-0.320</td><td>-0.141</td><td>-3.62</td><td>0.55 ± 0.04</td><td>0.63 ± 0.07</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-1.305</td><td>-0.011</td><td>9.16</td><td>0.18 ± 0.02</td><td>0.33 ± 0.06</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-0.467</td><td>-0.105</td><td>-4.16</td><td>0.15 ± 0.01</td><td>0.12 ± 0.01</td></tr><tr><td>G81.7219+0.5711</td><td>20:39:01.0501</td><td>+42:22:49.124</td><td>-1.415</td><td>-0.087</td><td>9.34</td><td>0.13 ± 0.01</td><td>0.13 ± 0.01</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>0.000</td><td>0.000</td><td>4.48</td><td>17.18 ± 0.12</td><td>17.88 ± 0.22</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>0.002</td><td>0.010</td><td>4.66</td><td>11.57 ± 0.09</td><td>11.89 ± 0.17</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>0.015</td><td>0.020</td><td>3.58</td><td>6.58 ± 0.08</td><td>6.61 ± 0.14</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>0.006</td><td>0.002</td><td>3.76</td><td>5.37 ± 0.07</td><td>5.42 ± 0.12</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>-0.025</td><td>0.004</td><td>3.94</td><td>2.89 ± 0.05</td><td>3.00 ± 0.09</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>-0.045</td><td>0.015</td><td>3.40</td><td>2.52 ± 0.04</td><td>2.56 ± 0.08</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>-0.027</td><td>-0.039</td><td>4.30</td><td>2.00 ± 0.03</td><td>1.92 ± 0.06</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>-0.001</td><td>0.044</td><td>4.12</td><td>1.72 ± 0.04</td><td>1.73 ± 0.07</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>-0.074</td><td>0.049</td><td>3.22</td><td>1.50 ± 0.03</td><td>1.54 ± 0.06</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>-0.030</td><td>0.009</td><td>3.04</td><td>1.20 ± 0.03</td><td>1.25 ± 0.05</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>-0.005</td><td>-0.089</td><td>5.02</td><td>0.27 ± 0.03</td><td>0.35 ± 0.05</td></tr><tr><td>G81.7444+0.5910</td><td>20:39:00.3722</td><td>+42:24:37.089</td><td>0.073</td><td>-0.040</td><td>5.38</td><td>0.17 ± 0.01</td><td>0.12 ± 0.01</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.000</td><td>0.000</td><td>-8.66</td><td>15.19 ± 0.10</td><td>15.67 ± 0.19</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.030</td><td>0.010</td><td>-8.48</td><td>11.28 ± 0.09</td><td>11.62 ± 0.17</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.014</td><td>-0.140</td><td>-5.78</td><td>7.93 ± 0.07</td><td>8.12 ± 0.13</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.018</td><td>-0.141</td><td>-5.96</td><td>5.47 ± 0.06</td><td>5.63 ± 0.11</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.036</td><td>-0.107</td><td>-6.86</td><td>4.29 ± 0.04</td><td>4.37 ± 0.08</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.005</td><td>0.003</td><td>-8.30</td><td>5.21 ± 0.06</td><td>5.52 ± 0.10</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.020</td><td>0.021</td><td>-8.12</td><td>4.44 ± 0.04</td><td>4.65 ± 0.08</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.053</td><td>-0.127</td><td>-5.60</td><td>2.52 ± 0.04</td><td>2.65 ± 0.07</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.009</td><td>-0.110</td><td>-6.68</td><td>2.56 ± 0.03</td><td>2.65 ± 0.05</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.022</td><td>-0.111</td><td>-6.32</td><td>1.88 ± 0.03</td><td>1.97 ± 0.06</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.002</td><td>-0.094</td><td>-7.04</td><td>1.97 ± 0.03</td><td>2.06 ± 0.05</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.038</td><td>-0.125</td><td>-5.42</td><td>1.72 ± 0.03</td><td>1.83 ± 0.05</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.020</td><td>-0.020</td><td>-8.84</td><td>1.81 ± 0.04</td><td>1.91 ± 0.07</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.020</td><td>-0.112</td><td>-6.50</td><td>1.58 ± 0.03</td><td>1.69 ± 0.06</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.123</td><td>-0.300</td><td>-2.36</td><td>1.33 ± 0.03</td><td>1.45 ± 0.06</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.023</td><td>0.001</td><td>-9.02</td><td>1.15 ± 0.02</td><td>1.17 ± 0.04</td></tr></tbody></table>
<table><caption>Table A.1. continued.</caption><thead><tr><th>Name<br>(1)</th><th>Ref. RA<br>(h:m:s)<br>(2)</th><th>Ref. Dec<br>(°:′:′′)<br>(3)</th><th>Δα<br>(arcsec)<br>(4)</th><th>Δδ<br>(arcsec)<br>(5)</th><th>V<sub>LSR</sub><br>(km s<sup>-1</sup>)<br>(6)</th><th>S<sub>v, Peak</sub><br>(Jy beam<sup>-1</sup>)<br>(7)</th><th>S<sub>v, Int.</sub><br>(Jy)<br>(8)</th></tr></thead><tbody><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.074</td><td>-0.363</td><td>-2.90</td><td>1.52 ± 0.03</td><td>1.70 ± 0.07</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.037</td><td>-0.018</td><td>-9.20</td><td>0.77 ± 0.02</td><td>0.89 ± 0.04</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.035</td><td>-0.124</td><td>-5.24</td><td>0.66 ± 0.02</td><td>0.70 ± 0.04</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.014</td><td>-0.374</td><td>-3.26</td><td>0.99 ± 0.04</td><td>1.19 ± 0.09</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.055</td><td>-0.090</td><td>-7.22</td><td>0.50 ± 0.02</td><td>0.56 ± 0.04</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.094</td><td>-0.431</td><td>-3.08</td><td>0.73 ± 0.04</td><td>0.86 ± 0.07</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.110</td><td>-0.365</td><td>-2.72</td><td>0.78 ± 0.05</td><td>1.03 ± 0.10</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.022</td><td>-0.188</td><td>-2.18</td><td>0.27 ± 0.02</td><td>0.34 ± 0.04</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.113</td><td>-0.051</td><td>-6.14</td><td>0.28 ± 0.02</td><td>0.34 ± 0.04</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>0.075</td><td>0.009</td><td>-2.54</td><td>0.22 ± 0.04</td><td>0.52 ± 0.11</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.106</td><td>0.155</td><td>-9.38</td><td>0.12 ± 0.02</td><td>0.34 ± 0.08</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.013</td><td>-0.173</td><td>-5.06</td><td>0.14 ± 0.02</td><td>0.19 ± 0.03</td></tr><tr><td>G81.7523+0.5908</td><td>20:39:01.9870</td><td>+42:24:59.261</td><td>-0.062</td><td>0.090</td><td>-7.40</td><td>0.10 ± 0.00</td><td>0.15 ± 0.01</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>0.059</td><td>0.013</td><td>-7.76</td><td>0.09 ± 0.00</td><td>0.11 ± 0.01</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>0.000</td><td>0.000</td><td>-1.28</td><td>3.90 ± 0.05</td><td>4.09 ± 0.09</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>-0.037</td><td>0.017</td><td>-1.46</td><td>2.38 ± 0.03</td><td>2.49 ± 0.06</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>-0.009</td><td>-0.003</td><td>-1.10</td><td>1.32 ± 0.03</td><td>1.49 ± 0.06</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>0.015</td><td>0.061</td><td>-1.64</td><td>0.45 ± 0.04</td><td>0.70 ± 0.09</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>0.171</td><td>-0.057</td><td>-1.82</td><td>0.16 ± 0.02</td><td>0.18 ± 0.03</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>-0.010</td><td>-0.140</td><td>-2.00</td><td>0.13 ± 0.01</td><td>0.16 ± 0.03</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>-0.359</td><td>-0.088</td><td>-0.74</td><td>0.09 ± 0.01</td><td>0.07 ± 0.01</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>-0.085</td><td>0.039</td><td>-0.92</td><td>0.08 ± 0.01</td><td>0.07 ± 0.01</td></tr><tr><td>G81.7655+0.5972</td><td>20:39:02.9329</td><td>+42:25:50.954</td><td>0.154</td><td>-0.025</td><td>-2.18</td><td>0.08 ± 0.01</td><td>0.06 ± 0.01</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>0.000</td><td>0.000</td><td>7.18</td><td>336.13 ± 0.96</td><td>346.76 ± 1.73</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.176</td><td>0.425</td><td>4.66</td><td>324.83 ± 1.15</td><td>336.46 ± 2.08</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.158</td><td>0.491</td><td>4.12</td><td>231.68 ± 0.89</td><td>238.50 ± 1.61</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.170</td><td>0.444</td><td>4.48</td><td>238.51 ± 1.07</td><td>246.98 ± 1.95</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.021</td><td>0.018</td><td>7.36</td><td>199.40 ± 0.95</td><td>206.21 ± 1.71</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.124</td><td>0.470</td><td>3.94</td><td>181.99 ± 0.88</td><td>188.35 ± 1.59</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.261</td><td>0.346</td><td>5.74</td><td>126.12 ± 0.76</td><td>128.52 ± 1.37</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.241</td><td>0.352</td><td>5.20</td><td>114.47 ± 0.74</td><td>116.80 ± 1.33</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.100</td><td>0.480</td><td>3.40</td><td>88.83 ± 0.80</td><td>90.64 ± 1.43</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.233</td><td>0.349</td><td>5.02</td><td>89.76 ± 0.83</td><td>91.29 ± 1.49</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.101</td><td>0.474</td><td>3.58</td><td>64.76 ± 0.74</td><td>66.49 ± 1.33</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.022</td><td>-0.138</td><td>6.82</td><td>41.99 ± 0.58</td><td>45.36 ± 1.08</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.122</td><td>-0.424</td><td>9.34</td><td>31.67 ± 0.61</td><td>35.80 ± 1.16</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.309</td><td>0.381</td><td>5.56</td><td>35.79 ± 0.69</td><td>37.69 ± 1.26</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.315</td><td>0.352</td><td>5.92</td><td>30.37 ± 0.61</td><td>31.77 ± 1.12</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.193</td><td>0.414</td><td>4.84</td><td>45.54 ± 0.76</td><td>47.59 ± 1.38</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.113</td><td>-0.435</td><td>9.52</td><td>16.86 ± 0.52</td><td>20.40 ± 1.04</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.229</td><td>0.347</td><td>5.38</td><td>25.80 ± 0.64</td><td>28.06 ± 1.19</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.199</td><td>0.444</td><td>3.76</td><td>20.97 ± 0.46</td><td>22.90 ± 0.87</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.076</td><td>0.000</td><td>7.72</td><td>6.08 ± 0.37</td><td>8.11 ± 0.79</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.130</td><td>0.442</td><td>3.22</td><td>10.31 ± 0.44</td><td>12.30 ± 0.88</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.172</td><td>-0.060</td><td>6.46</td><td>2.47 ± 0.26</td><td>4.53 ± 0.70</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>0.044</td><td>-0.287</td><td>9.16</td><td>2.89 ± 0.32</td><td>6.33 ± 0.97</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.025</td><td>-0.349</td><td>8.80</td><td>0.93 ± 0.12</td><td>2.03 ± 0.36</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.100</td><td>-0.261</td><td>7.00</td><td>7.90 ± 0.62</td><td>9.60 ± 1.22</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.099</td><td>-0.366</td><td>8.98</td><td>2.15 ± 0.24</td><td>4.27 ± 0.69</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.174</td><td>0.410</td><td>2.86</td><td>0.79 ± 0.10</td><td>1.88 ± 0.31</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>0.130</td><td>-0.086</td><td>9.70</td><td>0.90 ± 0.10</td><td>2.19 ± 0.35</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.566</td><td>0.371</td><td>6.28</td><td>1.00 ± 0.10</td><td>1.46 ± 0.24</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.426</td><td>0.400</td><td>4.30</td><td>7.97 ± 0.52</td><td>9.50 ± 1.04</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.387</td><td>-0.150</td><td>6.64</td><td>1.72 ± 0.25</td><td>3.76 ± 0.77</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.021</td><td>0.108</td><td>8.08</td><td>0.60 ± 0.09</td><td>1.17 ± 0.26</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.019</td><td>-0.017</td><td>6.10</td><td>0.85 ± 0.15</td><td>1.95 ± 0.48</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.053</td><td>-0.448</td><td>8.62</td><td>0.52 ± 0.06</td><td>0.92 ± 0.17</td></tr></tbody></table>
<table><caption>Table A.1. continued.</caption><thead><tr><th>Name<br>(1)</th><th>Ref. RA<br>(h:m:s)<br>(2)</th><th>Ref. Dec<br>(°:′:′′)<br>(3)</th><th>Δα<br>(arcsec)<br>(4)</th><th>Δδ<br>(arcsec)<br>(5)</th><th>V<sub>LSR</sub><br>(km s<sup>-1</sup>)<br>(6)</th><th>S<sub>v, Peak</sub><br>(Jy beam<sup>-1</sup>)<br>(7)</th><th>S<sub>v, Int.</sub><br>(Jy)<br>(8)</th></tr></thead><tbody><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>0.113</td><td>-0.042</td><td>8.44</td><td>0.27 ± 0.05</td><td>0.71 ± 0.18</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>0.081</td><td>-0.498</td><td>9.88</td><td>0.27 ± 0.04</td><td>0.48 ± 0.10</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>0.127</td><td>0.048</td><td>10.42</td><td>0.27 ± 0.04</td><td>0.69 ± 0.14</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.213</td><td>0.344</td><td>2.50</td><td>0.26 ± 0.03</td><td>0.64 ± 0.11</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>0.066</td><td>-0.038</td><td>10.06</td><td>0.23 ± 0.03</td><td>0.67 ± 0.12</td></tr><tr><td>G81.8713+0.7807</td><td>20:38:36.4232</td><td>+42:37:34.756</td><td>-0.574</td><td>0.603</td><td>2.68</td><td>0.14 ± 0.03</td><td>0.28 ± 0.07</td></tr><tr><td>G82.3079+0.7296</td><td>20:40:16.6487</td><td>+42:56:29.288</td><td>0.000</td><td>0.000</td><td>10.42</td><td>23.82 ± 0.18</td><td>24.19 ± 0.31</td></tr><tr><td>G82.3079+0.7296</td><td>20:40:16.6487</td><td>+42:56:29.288</td><td>0.027</td><td>0.006</td><td>10.24</td><td>15.28 ± 0.13</td><td>15.85 ± 0.23</td></tr><tr><td>G82.3079+0.7296</td><td>20:40:16.6487</td><td>+42:56:29.288</td><td>0.006</td><td>0.001</td><td>10.78</td><td>6.38 ± 0.07</td><td>6.46 ± 0.12</td></tr><tr><td>G82.3079+0.7296</td><td>20:40:16.6487</td><td>+42:56:29.288</td><td>-0.079</td><td>-0.048</td><td>10.60</td><td>4.47 ± 0.05</td><td>4.54 ± 0.09</td></tr><tr><td>G82.3079+0.7296</td><td>20:40:16.6487</td><td>+42:56:29.288</td><td>0.072</td><td>0.006</td><td>9.88</td><td>1.27 ± 0.03</td><td>1.30 ± 0.05</td></tr><tr><td>G82.3079+0.7296</td><td>20:40:16.6487</td><td>+42:56:29.288</td><td>-0.042</td><td>-0.021</td><td>11.14</td><td>1.04 ± 0.04</td><td>1.27 ± 0.07</td></tr><tr><td>G82.3079+0.7296</td><td>20:40:16.6487</td><td>+42:56:29.288</td><td>0.212</td><td>-0.053</td><td>9.52</td><td>0.24 ± 0.03</td><td>0.39 ± 0.08</td></tr></tbody></table>
<figure><img src="image_16.png" /><figcaption>**Fig. A.1.** Dust continuum emission images from: Herschel PACS (70 and 160 µm) and SPIRE (250, 350, and 500 µm); JCMT SCUBA-2 (450 and 850 µm); and either IRAM MAMBO/MAMBO-2 (1.2 mm) or CSO Bolocam (1.1 mm) instruments. The color scale is in units of Jy beam⁻¹. The location of the methanol masers detected with the VLA are indicated by the black crosses. The Galactic names of the GLOSTAR sources are indicated at the top of the bottom-left panels. The beams are shown at the bottom-left corner of each panel. Saturated pixels are shown in gray.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>Fig. A.1. continued.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>Fig. A.1. continued.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>Fig. A.1. continued.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>Fig. A.1. continued.</figcaption></figure>
<figure><img src="image_21.png" /></figure>

**Fig. A.2.** Infrared images (color scale and contours) of the environment around the location of the VLA-detected methanol masers. The top-left panel shows a three-color map constructed using Spitzer 3.6 µm (blue), 4.5 µm (green), and 8 µm (red) images. The other five panels show infrared emission for each of the Spitzer bands. The *n*th white contour is at (√2)ⁿ × S<sub>max</sub> × *p*, where S<sub>max</sub> is the maximum flux shown in the color bar (MJy sr⁻¹) for each panel, *n*=0, 1, 2 ..., and *p* is equal to 30%. The green crosses mark the positions of the maser features seen in B-array configuration maps. The Galactic names of the GLOSTAR sources are indicated at the top of the top-left panels. The beams are shown at the bottom-left corner of each panel. Saturated pixels are shown in gray.
<figure><img src="image_22.png" /><figcaption>Fig. A.2. continued.</figcaption></figure>
<figure><img src="image_23.png" /><figcaption>Fig. A.2. continued.</figcaption></figure>
<figure><img src="image_24.png" /></figure>

Fig. A.2. continued.
<figure><img src="image_25.png" /></figure>

Fig. A.2. continued.
<figure><img src="image_26.png" /></figure>

**Fig. A.3.** Distribution of all the maser spots (filled circles) contributing to the maser features for B-array data. The background image shows B-configuration continuum emission (blue scale and contours). The spots are color coded by LSR velocity (color bar). The *n*th blue contour is at ($\sqrt{2}$)ⁿ × *S*<sub>max</sub> × *p*, where *S*<sub>max</sub> is the maximum continuum flux shown in the blue bar for each panel, *n*=0, 1, 2 ..., and *p* is equal to 10%. Positions of radio continuum sources identified in previous observations are marked with black stars along with their names. Maser spots detected with the EVN by Rygl et al. (2012) are shown as filled squares. The gray grid shows the Galactic coordinate system. The beam size of the radio continuum maps is 1′′5 × 1′′5. The maser position errors are ≈0′′2 (see Sect. 3).
<figure><img src="image_27.png" /><figcaption>Fig. A.3. continued.</figcaption></figure>
<figure><img src="image_28.png" /><figcaption>Fig. A.3. continued.</figcaption></figure>
## Appendix B: SED fitting

Figure B.1 displays the SED fitting of the dust cores that host maser emission.

<figure><img src="image_29.png" /></figure>

**Fig. B.1.** SED fitting for the mass and dust temperature derivation of the cores with associated maser emission. A model consisting of a warm and a cold component (solid line) is fit when the 24 and 70 µm points are available. Otherwise, the data are fit with a single gray-body model. The two dotted lines present the two components separately. The IRAC 3.6–8 µm points – likely originating from a third, hot and inner component – were not used in the fitting.
<figure><img src="image_30.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
## Appendix C: Previously detected masers

Methanol masers at 6.7 GHz detected in previous surveys by Hu et al. (2016) and Yang et al. (2019) are listed in Table C.1.

<table><caption>**Table C.1.** Methanol masers toward Cygnus X reported in previous surveys.</caption><thead><tr><th>RA</th><th>Dec</th><th>Name</th><th>V<sub>LSR</sub><br>(km s<sup>-1</sup>)</th><th>S<sub>v, Peak</sub><br>(Jy beam<sup>-1</sup>)</th></tr></thead><tbody><tr><td>20 29 24.937</td><td>+40 11 19.29</td><td>G078.886+0.708</td><td>-7.02</td><td>0.54</td></tr><tr><td rowspan="3">20 30 50.685</td><td rowspan="3">+41 02 27.50</td><td rowspan="3">G079.735+0.990</td><td>-5.44</td><td>23.17</td></tr><tr><td>-3.86</td><td>5.00</td></tr><tr><td>-3.16</td><td>1.91</td></tr><tr><td rowspan="3">20 37 00.980</td><td rowspan="3">+41 34 55.05</td><td rowspan="3">G080.861+0.383</td><td>-11.06</td><td>2.53</td></tr><tr><td>-4.04</td><td>8.52</td></tr><tr><td>-1.93</td><td>1.17</td></tr><tr><td>20 39 00.879</td><td>+42 22 48.93</td><td>G081.721+0.571</td><td>-3.51</td><td>0.86</td></tr><tr><td rowspan="2">20 39 00.378</td><td rowspan="2">+42 24 36.92</td><td rowspan="2">G081.744+0.590</td><td>-2.63</td><td>2.96</td></tr><tr><td>3.86</td><td>3.82</td></tr><tr><td rowspan="3">20 39 01.573</td><td rowspan="3">+42 24 59.09</td><td rowspan="3">G081.752+0.590</td><td>4.57</td><td>5.53</td></tr><tr><td>-8.60</td><td>8.35</td></tr><tr><td>-6.67</td><td>3.43</td></tr><tr><td rowspan="3">20 39 02.274</td><td rowspan="3">+42 25 50.81</td><td rowspan="3">G081.765+0.597</td><td>-5.80</td><td>3.87</td></tr><tr><td>-2.81</td><td>1.51</td></tr><tr><td>-2.28</td><td>0.90</td></tr><tr><td rowspan="8">20 38 36.423</td><td rowspan="8">+42 37 35.01</td><td rowspan="8">G081.871+0.780</td><td>-1.23</td><td>2.54</td></tr><tr><td>3.51</td><td>42.91</td></tr><tr><td>4.04</td><td>99.94</td></tr><tr><td>4.74</td><td>166.73</td></tr><tr><td>5.27</td><td>76.09</td></tr><tr><td>5.80</td><td>67.27</td></tr><tr><td>7.20</td><td>174.93</td></tr><tr><td>9.48</td><td>17.48</td></tr><tr><th></th><th></th><th></th><th>V<sub>LSR</sub>, ΔV<br>(km s<sup>-1</sup>)</th><th>S<sub>v, Peak</sub><br>(Jy)</th></tr><tr><td>20 23 23.65</td><td>+37 35 34.3</td><td>G76.093+0.158</td><td>4.92 (4.27, 6.90)</td><td>0.49</td></tr><tr><td>20 29 20.35</td><td>+40 11 37.5</td><td>G78.882+0.723</td><td>-6.94 (-7.41, -6.51)</td><td>0.95</td></tr><tr><td>20 30 22.77</td><td>+40 09 23.2</td><td>G78.969+0.541</td><td>4.74 (4.01, 5.48)</td><td>1.24</td></tr><tr><td>20 30 50.70</td><td>+41 02 28.8</td><td>G79.736+0.991</td><td>-5.51 (-6.63, -2.84)</td><td>19.4</td></tr><tr><td>20 37 01.02</td><td>+41 34 56.9</td><td>G80.862+0.383</td><td>-4.01 (-4.57, -1.77)</td><td>9.67</td></tr><tr><td>20 39 02.01</td><td>+42 24 59.3</td><td>G81.752+0.591</td><td>-8.57 (-9.39, 5.04)</td><td>15.4</td></tr><tr><td>20 37 47.39</td><td>+42 38 39.0</td><td>G81.794+0.911<sup>(a)</sup></td><td>7.19 (-4.73, 7.54)</td><td>0.59</td></tr><tr><td>20 38 36.67</td><td>+42 37 30.6</td><td>G81.871+0.779</td><td>4.53 (2.55, 9.87)</td><td>225.7</td></tr><tr><td>20 40 16.72</td><td>+42 56 28.6</td><td>G82.308+0.729</td><td>10.3 (10.0, 11.3)</td><td>58.4</td></tr></tbody></table>

**Notes.** The top table lists maser components reported in the catalog of Hu et al. (2016) and gives the LSR velocity of the component and the flux density. The bottom table lists masers in the catalog of Yang et al. (2019) and gives LSR velocity of peak emission, the LSR velocity interval of the maser emission in parentheses, and peak flux density. (a) Not detected in GLOSTAR.

## Appendix D: Methanol absorption toward DR21

Based on our absorption spectrum, we now constrain the methanol abundance in DR21 M. In local thermodynamic equilibrium (LTE), the column density in a spectral line's lower energy level is given by the line's observed optical depth, $\tau$, integrated over the line profile,

$$N_1 = \frac{8\pi k g_1}{c^3 h g_u} T_{ex} \int \tau dv, \qquad (D.1)$$

from which the total column density of the molecule, $N(\text{CH}_3\text{OH})$, can be calculated:

$$N = \frac{N_1}{g_1} Q(T_{\text{rot}}) e^{E_1/kT}. \qquad (D.2)$$

Here, $k$, $h$, and $c$ are the Planck constant, the Boltzmann constant, and the speed of light, respectively, and $g_1$ and $g_u$ are the degeneracies of the line's lower and upper energy levels, respectively ($g_1 = 13$ and $g_u = 11$); $Q$ is the partition function, and $E_1$ is the energy above ground state of the lower energy level (48.7 K).
<figure><img src="image_31.png" /><figcaption>**Fig. D.1.** Dust spectral energy distribution of DR21 M. The best fit SED is marked by the blue line and is the sum of the simultaneously fit gray body (right) and blackbody (left) components (green dashed lines).</figcaption></figure>

Integrating our absorption spectrum from -17 to +2 km s⁻¹, we determine an integrated optical depth of 0.47 km s⁻¹.

By definition, under LTE, the excitation temperature and rotation temperature of the line are the same and are also equal to the kinetic temperature, $T_{ex} = T_{rot} \equiv T_{kin}$. A question then arises as to the appropriate temperature for the calculation of the methanol column density.

Modeling a variety of radio-wavelength OH absorption lines, Jones et al. (1994) derived an H₂ density of (1.8 ± 0.7) × 10⁷ cm⁻³ and a kinetic temperature of 175 K for the part of DR21 M's dense molecular core from which the OH-bearing gas originates, namely the dense PDR interface with the compact HII region (see Sect. 4.5). These numbers might not be representative of the more extended condensation that gives rise to the methanol absorption.

Therefore, we characterized the molecular gas associated with DR21 M by fitting its dust SED (Fig. D.1) obtained through aperture photometry of mid-infrared to submillimeter continuum maps to derive its temperature, source size, H₂ column, and volume density.

Using the same procedures as described in König et al. (2021), we determined the flux densities from Herschel/Hi-GAL (Molinari et al. 2010) and MSX (Price et al. 2001) continuum maps in nine bands using an aperture-annulus approach. As the 160 µm PACS band and the D-configuration maps of methanol have similar beam widths (i.e., 13′′5 and 13′′, respectively), we estimated the source size (i.e., the FWHM₁₆₀) through a two-dimensional Gaussian fit to the emission in the 160 µm PACS band centered on the peak of the emission. The same position, aperture, and annulus sizes were subsequently used for all bands to measure the source and background flux, from which the background corrected source flux was calculated. The resulting SED was then fit with a two-component model consisting of a blackbody for the warm emission originating from the evolving star and a gray body representing the cold dust envelope, yielding the dust temperature, $T_{dust}$.

With the ratio of the measured flux density, $F_{350}$, to the gray body's intensity, $B_{350}(T_{dust})$, reflecting the dust emissivity, we calculated the H₂ column density according to

$$N_{\mathrm{H}_2} = \frac{F_{350}}{B_{350}(T_{\mathrm{dust}})} \cdot \frac{\gamma}{\kappa_{350}} \cdot \frac{1}{\Omega_{\mathrm{app}} \cdot \mu_{\mathrm{H}_2} \cdot m_{\mathrm{H}}} \quad (\mathrm{D.3})$$

where $\gamma = 100$ (the gas-to-dust ratio), $\kappa_{350} = 1.1 \text{ cm}^2\text{g}^{-1}$ (the dust opacity at 350 µm calculated as the mean of the dust models from Ossenkopf & Henning 1994), $\Omega_{src}$ is the source

<table><caption>**Table D.1.** Summary of the physical properties of DR21 M.</caption><thead><tr><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>Apparent source size FWHM<sub>160</sub>:</td><td>37′′8</td></tr><tr><td>Linear source radius r<sub>src</sub>:</td><td>0.128 pc</td></tr><tr><td>Dust temperature T<sub>dust</sub>:</td><td>36.0 ± 2.1 K</td></tr><tr><td>Average H<sub>2</sub> column density N<sub>H<sub>2</sub></sub>:</td><td>(1.7 ± 0.4) × 10<sup>23</sup> cm<sup>-2</sup></td></tr><tr><td>Average H<sub>2</sub> volume density n<sub>H<sub>2</sub></sub>:</td><td>(2.2 ± 0.7) × 10<sup>5</sup> cm<sup>-3</sup></td></tr></tbody></table>

solid angle, $\mu_{H_2} = 2.8$ is the mean molecular weight of the interstellar medium with respect to a hydrogen molecule (Kauffmann et al. 2008), and $m_H$ is the mass of a hydrogen atom.

Using a distance of $d = 1.4$ kpc to DR21 and taking into account the fact that the measured source size FWHM is convolved with the beam, we determined the linear source size:

$$r_{\mathrm{src}} = d \cdot \tan \left( \sqrt{\mathrm{FWHM}_{160}^2 - \theta_{160}^2} \right) \quad (\mathrm{D.4})$$

with $\theta_{70} = 13′′51$ the beamwidth of the 160 µm band.

Using the linear source size, $r_{src}$, we determined the H₂ volume density:

$$ \begin{aligned} n_{\mathrm{H}_2} &= \frac{M}{V} \cdot \frac{1}{\mu_{\mathrm{H}_2} \cdot m_{\mathrm{H}}} \\ &= \frac{3}{4\pi \cdot r_{\mathrm{src}}^3} \cdot \frac{d^2 \cdot F_{350}}{B_{350}(T_{\mathrm{dust}})} \cdot \frac{\gamma}{\kappa_{350}} \cdot \frac{1}{\mu_{\mathrm{H}_2} \cdot m_{\mathrm{H}}} \end{aligned} \quad (\mathrm{D.5}) $$

We summarize the physical properties of DR21 M in Table D.1.

At the derived density, the dust and gas temperatures, $T_{dust}$ and $T_{kin}$, respectively, are well coupled (i.e., equal to 36 K). If we assume that the CH₃OH energy levels are in LTE at this temperature, we can derive the total CH₃OH column density using Eqs. (D.1) and (D.2), to calculate a total methanol column density, $N(\mathrm{CH}_3\mathrm{OH})$, of $2.8 \times 10^{17} \text{ cm}^{-2}$ and a CH₃OH abundance of $1.8 \times 10^{-6}$ relative to H₂ averaged over the DR21 M molecular core, for which we have interpolated the value for the partition function, $Q(36K)$, as 977 from the values tabulated on a CH₃OH data web page<fnref n="11" /> of the Cologne Database for Molecular Spectroscopy (CDMS; Müller et al. 2005).

While CH₃OH abundances of this magnitude (or even higher) have been determined for compact hot ($T \sim 150$ K) molecular cores that often surround newly formed high-mass stars (see, e.g., Menten et al. 1988; Bonfand et al. 2017; Molet et al. 2019), $1.8 \times 10^{-6}$ appears to be a high value for the more extended, cooler, and more developed DR21 M molecular clump. We note that the much higher temperature (175 K) invoked from modeling the OH radio absorption lines characterizes a more limited hotter region that represents the dense PDR interface of the molecular clump with the compact HII region.

As discussed in Sect. 4.5, for certain ranges of density and temperature, the $S_1 \to 6_0 A^+$ transition can be “overcooled” (i.e., show deeper “enhanced” absorption) than expected under LTE, even against the cosmic microwave background radiation: Pandian et al. (2008) observed this line in absorption toward the hot corinos<fnref n="12" /> NGC 1333-IRAS 4A and 4B, which do not show

<footnote n="11">¹¹ https://cdms.astro.uni-koeln.de/cgi-bin/cdmsinfo?file=e032504.cat</footnote>

<footnote n="12">¹² Hot corinos are the equivalent of hot cores surrounding low-mass protostellar objects or low-mass young stellar objects.</footnote>
radio continuum emission. Were this the case, the CH3OH abundance presented above would be an overestimate. To investigate this, we used the non-LTE radiative transfer program RADEX (van der Tak et al. 2007) to study the line’s excitation. For the gas densities and temperatures determined by our modeling, we find that the line is indeed predicted to show enhanced absorption, that is, its excitation temperature is lower than the temperature of the cosmic microwave background. This means that the methanol abundance derived above represents an upper limit to the true value. Meaningful determinations of N(CH3OH) for DR21 M require constraints provided by data from thermally excited methanol lines.

## Appendix E: Expected number of background sources

Following Anglada et al. (1998), we can estimate the number of expected background sources, Nbg, inside a field of diameter θF as

$$N_{bg} = 1.4 \left\{ 1 - \exp\left[-0.0066 \left(\frac{\theta_F}{\text{arcmin}}\right)^2 \left(\frac{\nu}{5 \text{ GHz}}\right)^2\right] \right\} \times \left(\frac{S_0}{\text{mJy}}\right)^{-0.75} \left(\frac{\nu}{5 \text{ GHz}}\right)^{-2.52}, \quad (E.1)$$

where S0 is the detectable flux density threshold and ν the observing frequency. In our observations, ν = 5.8 GHz and S0 = 3 × rms = 0.18 mJy. For field sizes of θF = 1, 0.5, and 0.017 arcmin, we expect Nbg = 0.03, 0.008, and 8.6 × 10⁻⁶.