Maser Activity of Organic Molecules toward Sgr B2(N)

At centimeter wavelengths, single-dish observations have suggested that the Sagittarius (Sgr) B2 molecular cloud at the Galactic Center hosts weak maser emission from several organic molecules, including CH2NH, HNCNH, and HCOOCH3. However, the lack of spatial distribution information on these new maser species has prevented us from assessing the excitation conditions of the maser emission as well as their pumping mechanisms. Here, we present a mapping study toward Sgr B2 north (N) to locate the region where the complex maser emission originates. We report the first detection of the Class I methanol (CH3OH) maser at 84 GHz and the first interferometric map of the methanimine (CH2NH) maser at 5.29 GHz toward this region. In addition, we present a tool for modeling and fitting the unsaturated molecular maser signals with non-LTE radiative transfer models and Bayesian analysis using the Markov Chain Monte Carlo approach. These enable us to quantitatively assess the observed spectral profiles. The results suggest a two-chain-clump model for explaining the intense CH3OH Class I maser emission toward a region with low continuum background radiation. By comparing the spatial origin and extent of maser emission from several molecular species, we find that the 5.29 GHz CH2NH maser has a close spatial relationship with the 84 GHz CH3OH Class I masers. This relationship serves as observational evidence to suggest a similar collisional pumping mechanism for these maser transitions.


INTRODUCTION
Maser emission from molecular species has provided critical insights into the kinetic properties of young stars and has even enabled monitoring of their accretion processes (Hunter et al. 2018, and references therein).Energetic events (e.g., outflows and collimated jets) and infrared radiation surrounding massive protostars provide the energy to overpopulate the upper energy levels of specific molecular transitions and amplify the molecular emission intensities (Moscadelli et al. 2016;Sanna et al. 2016).Common species exhibiting maser emission include, but are not limited to, OH, SiO, H 2 O, H 2 CO, NH 3 , and CH 3 OH.An empirical taxonomy based on pumping mechanisms classifies molecular masers into Corresponding author: Ci Xue cixue@mit.edutwo categories (Batrla et al. 1987;Cyganowski et al. 2009); H 2 O and Class I CH 3 OH masers are most often excited via collisions and tend to reside in shocked material associated with outflows and expanding H II regions (e.g.Baudry et al. 1974;Elitzur & Fuqua 1989;Voronkov et al. 2014;Towner et al. 2017), while OH and Class II CH 3 OH masers are excited via radiation and are generally close to young stellar objects and luminous infrared sources (e.g., Sobolev & Deguchi 1994;Sobolev et al. 1997;Cragg et al. 2002Cragg et al. , 2005)).The most extensively studied Class I CH 3 OH masers include the 4 −1 − 3 0 E transition at 36 GHz and the 7 0 − 6 1 A transition at 44 GHz, while Class II CH 3 OH masers include the 5 1 − 6 0 A transition at 6.7 GHz and the 2 0 − 3 −1 E transition at 12.2 GHz (e.g.Voronkov et al. 2014;Song et al. 2022).
By correlating the maser emission from various species with a sample of star-forming regions at different stages, statistical analyses have also suggested that masers may help to trace the evolutionary phases of massive stars (e.g. Figure 6 in Breen et al. 2010).For example, Forster & Caswell (1989) found that isolated H 2 O masers tend to appear in sources at earlier evolutionary stages than OH masers based on a survey toward 74 star-forming regions.In addition, Class II CH 3 OH masers, which appear exclusively in massive star-forming regions (Breen et al. 2013;Xu et al. 2008), typically have a significant overlapping timeline with H 2 O masers and arise before the OH masers and ultra-compact (UC) H II regions turn on (Walsh et al. 2003;Szymczak et al. 2005).
However, within the picture of massive star formation, the position for Class I CH 3 OH masers in an evolutionary sequence is less clear.Initial comparisons of the mid-infrared properties of the sample of CH 3 OH maser sources suggested that Class I masers are signposts of the earliest stage of massive star formation, even prior to H 2 O masers (Ellingsen 2006).In contrast, detailed mapping of multiple Class I CH 3 OH maser transitions toward a luminous young stellar object, IRAS 16547-4247, led Voronkov et al. (2006) to argue that Class I CH 3 OH masers can appear earlier than the Class II masers but also last long enough to coexist with the OH masers.Subsequent observations toward the massive protocluster G18.67+0.03further suggested that the Class I CH 3 OH masers might associate with both young sites and relatively evolved UC H II regions within the same protocluster (Cyganowski et al. 2012).In any case, high angular resolution imaging is essential for determining which protostar(s) within a protocluster currently power these masers (Towner et al. 2017).
Recently, new maser species, such as HNCO and CH 2 NH, have been discovered to exhibit compelling roles in inferring fine-scale structures and gas dynamics in the maser region.In particular, the HNCO maser was found to be associated with accretion flows, which helped to build an arm-spiral-structure model of one particular high-mass young stellar object (Chen et al. 2020); and the CH 2 NH maser toward galactic nuclei helps to probe the inner nuclear processes in extragalactic environments (Gorski et al. 2021).The pumping mechanism of these new maser species, however, is less constrained, as is their relation to other common maser species.Observations of these species in maser-active regions that contain both known collisionally and radiatively pumped masers are therefore critical as a first step in understanding their pumping mechanisms.
The high-mass star-forming region Sgr B2 North (N) displays rich maser activity from a diverse set of molecules, ranging from common maser species to rare complex molecules (Morita et al. 1992;McGrath et al. 2004;Hoffman et al. 2007;Qiao et al. 2018;Fujisawa et al. 2014).The single-dish PRebiotic Interstellar MOlecular Survey (PRIMOS) project has identified numerous species exhibiting maser activity at centimeter wavelengths.In fact, the discovery of interstellar carbodiimide (HNCNH) was only achievable due to its maser emission (McGuire et al. 2012).Radiative transfer calculations have also revealed that most emission lines of both methyl formate (HCOOCH 3 ) and CH 2 NH detected in PRIMOS are weak masers amplifying the background continuum radiation of Sgr B2(N) at cm wavelengths (Faure et al. 2014(Faure et al. , 2018)).However, the large ∼2.5 ′ × 2.5 ′ field of view of PRIMOS at C band (∼ 5 GHz) includes both the North and Main components of Sgr B2, complicating the interpretation of the spectral profiles.Therefore, it is crucial to have interferometric observations to disentangle in detail the various possible spatial components of maser emission.
In this work, we present a rigorous imaging study of maser emission from CH 3 OH and CH 2 NH toward Sgr B2(N).We describe the interferometric observations at centimeter (cm) and millimeter (mm) wavelengths in Section 2. A description of the methodology employed to model maser emission is described in Section 3, which includes the implementation of a Bayesian inference of physical characteristics.Section 4 presents both the spectra and spatial distribution per maser species.Section 5 contains discussions on the metastable states of maser transitions (Section 5.1), pumping mechanisms inferred from the spatial correlations (Section 5.2), derived physical conditions of masing gas clumps (Section 5.3), and a caveat on modeling maser intensity (Section 5.4).

OBSERVATIONS
The interferometric data used here were acquired from two data sets involving the Karl G. Jansky Very Large Array (VLA) in C-band (nominal frequency coverage 4-8 GHz) and the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3 (nominal frequency coverage 84-116 GHz).
VLA Observations -The VLA observations were carried out on 2016 May 27 (project code: 16A-076, P.I.: B. McGuire) with the C-band receiver.The correlator was configured to deliver 24 subbands or spectral windows, with bandwidths of 4, 8, and 16 MHz, to target various spectral lines.Additionally, nine 128 MHz spectral windows were set up to help with calibration, if needed.The results presented here are based on three spectral windows that delivered channel widths between 0.485 km s −1 and 2.31 km s −1 , depending on the corresponding bandwidth and observing frequency.Detailed observational parameters of each spectral window are presented in Table 1.The phase center of the observations was targeted toward Sgr B2 (N), with the field center at α J2000 = 17 h 47 m 19 s .80,δ J2000 = −28 • 22 ′ 17 ′′ .00.The bandpass and absolute flux density scale calibrator was 3C286, while the complex gain calibrator was J1744-3116.The array was in the B configuration with 27 antennas, which covered baselines from 250 m up to 11 km and resulted in the maximum recoverable scale of 37 ′′ at 4 GHz.
The CASA software package (version 5.6.1,CASA Team et al. 2022) was used for data reduction and imaging analysis of the VLA observations.For each spectral window, the continuum emission was subtracted in the uv-plane before imaging the spectral cubes.We used the tclean task with the auto-masking algorithm for imaging with the natural weighting scheme and a cell size of 0.27 ′′ .The resulting size of the synthesized beam (FWHM) of each image cube has a median value of ∼3.8 ′′ × 1.5 ′′ .The median of the resultant typical RMS noise levels is ∼0.9 mJy beam −1 .
ALMA Observations -The ALMA observational data were acquired from the ALMA Science Archive of the survey project entitled "Exploring Molecular Complexity with ALMA (EMoCA)" (project code: 2011.0.00017.S, P.I.: A. Belloche).A detailed description of the EMoCA data has been presented in Belloche et al. (2016).In this work, two measurement sets, spectral window 0 of setup S1 and spectral window 0 of setup S4 (Xue et al. 2019), were used for imaging the 84 GHz and 95 GHz CH 3 OH transitions.Different from Xue et al. (2019), we used the tclean task with the auto-masking algorithm for the imaging analysis, consistent with the VLA imaging process.The resultant synthesized beam size is 2.23 ′′ × 1.67 ′′ at 84 GHz and 1.94 ′′ × 1.49 ′′ at 95 GHz with a median RMS value of 4 mJy beam −1 .

ANALYSIS
In order to constrain the excitation conditions quantitatively, we used the molsim program to model the molecular spectra with non-LTE radiative transfer models and to explore the physical parameter space with the Markov-Chain Monte Carlo (MCMC) approach (Lee et al. 2023).

Non-LTE Radiative Transfer Modeling
We implemented the RADEX non-LTE radiative transfer code (van der Tak et al. 2007) into molsim written in Python.The module was accelerated with Numba to provide the high performance required for MCMC analysis.
RADEX considers a gas clump consisting of a targeted molecular species and its collision partners, such as H 2 and e -.The gas clump is embedded in a local background radiation field and has a certain geometry (e.g., slab, uniform sphere, expanding sphere).With this physical picture, RADEX iteratively solves for the level population by incorporating an equilibrium with inelastic collisions and the local radiation field.Lastly, it computes the peak optical depth, excitation temperature and peak intensity of each line based on the level population and the line properties, assuming the local radiation field as the only source of background continuum emission.The model can be described with physical parameters including kinetic temperatures (T k ), densities of collision partners (n c ), background radiation field temperatures (T bg ), source velocity (V lsr ), line widths (dV ), and molecular column densities (N col ).
In our implementation, we construct a modification and extend it to support a continuum temperature (T cont ) different from the local radiation field temperature (T rad 1 ) and allow multiple components aligned perfectly along the line of sight (LOS).In the case of a population inversion, T cont serves as background radiation that stimulates the maser gain and is amplified by the foreground gas clump.We assume the continuum source and the gas components to be spatially distant from each other such that the continuum radiation is not included in the excitation analysis and the embedded radiation field of each component is only dominated by the local radiation field.Detailed justifications for this assumption will be discussed in Section 5.3.With this assumption, the level population of each component is computed with a distinct excitation condition using RADEX.The frequency-dependent optical depth and source function of the i-th component are reconstructed by summing over all lines as follows: Here, τ peak i,j and T ex,i,j are the optical depth and excitation temperature returned by RADEX for the j-th line of the i-th component, B ν is the Planck function, ϕ i is the line profile for the i-th component using the Gaussian function G, and τ i and S i are the total optical depth and source function with all lines being accounted for.Using these definitions, for an n-component model where the continuum radiation is first processed by the first component, then the second component and so on, the final intensity with the continuum subtracted can be expressed as In the case of two components, the continuumsubtracted intensity is It is worthwhile to note that the expressions reduce back to that of RADEX in the case of one component when the continuum temperature is equal to the local radiation temperature (i.e.T cont = T rad = T bg ).The calculations in molsim reproduces the results returned by RADEX down to machine precision when the results are converged.
In this work, we assume that the local radiation is contributed by the cosmic microwave background only, i.e., T rad = 2.725 K (Mather et al. 1999) and a uniform spherical geometry.We only consider H 2 as the collision partner, assuming a kinetic temperature-dependent thermal ortho-to-para ratio (OPR), which can be approximated by: We left all the rest as free parameters and allowed them to be fitted using the MCMC approach as described below.

MCMC Fitting
For line fitting, we used a Bayesian approach to compare models with different parameters and calculate their probability distributions.This approach provides an inference by conditioning the data on priors and sampling posterior distributions.In practice, we use the Markov chain Monte Carlo (MCMC) approach to explore parameter space.
We use a logarithmic likelihood distribution with the formula below: where x is the observed spectral line intensity, µ is the model line intensity, and σ is the uncertainty of the observations including an estimated 10% systematic uncertainty with the standard baseline deviation.The distribution applies to the case where the dataset has a non-uniform uncertainty, in which case we compute the uncertainties for each spectral chunk separately.
The resulting corner plots are shown in Appendix B. A corner plot shows the probability distribution for each pair of parameters and for each parameter individually.In the following, the fitted results are represented by 50 th , 16 th , and 84 th percentile as the representation value and the range of uncertainty, also known as the 68% confidence interval.This choice is equivalent to choosing the mean and +/−1σ for a normal distribution, which is the case for most of our posterior distributions.One of the advantages of the MCMC technique over a traditional least-squares fit approach is that a larger parameter space can be explored more efficiently.Rather than only returning the parameters that maximize the likelihood function (e.g.least-squares fitting2 ), it returns the actual probability of the parameters, which can be used to assess whether parameters are highly co-   variant with one another.The correlation between parameters manifests itself as non-separable distributions in the corner plots.

RESULTS
We investigate the maser activities of two transitions in this work.The detailed spectroscopic parameters for each maser transition along with the observational results of the corresponding emission features are summarized in Table 2.For CH 3 OH, we examined the molecular spectra extracted toward two distinct regions, maser spot 1 (MS1, α J2000 = 17 h 47 m 19 s .730,δ J2000 = −28 • 22 ′ 04 ′′ .580)and maser spot 2 (MS2, α J2000 = 17 h 47 m 20 s .393,δ J2000 = −28 • 22 ′ 12 ′′ .410).Toward both regions, we performed MCMC analyses with a two-chained-component model, as presented in Section 4.1.1 for MS1 and 4.1.2for MS2 respectively.For CH 2 NH, as there is no detectable signal toward MS1, we examined its spectrum extracted toward MS2 only and performed the MCMC analysis with a one-component model as presented in Section 4.2.
Empirically, a brightness temperature (T B ) of an emission line much greater than the source physical (kinetic) temperature is considered as an indication of a maser emission rather than a thermally-excited line (Voronkov et al. 2006;Gargaud et al. 2011).As elaborated in the following sections, our MCMC analyses suggest the T k for the masing clumps in Sgr B2(N) to be ∼29 − 37 K.In contrast, as summarized in Table 2, the observed T B values of the molecular emission investigated in this work (> 120 K) are significantly higher than the physical temperatures.Furthermore, we also performed a spectral analysis assuming a single excitation temperature to equal T k for comparison and found that the observed spectral profiles deviate substantially from the thermal model.Additional evidence for maser activities, including the non-Gaussian shape of line profiles, the compact emitting regions, the negative opacity values, and the metastable upper-energy levels, will be discussed more in the following sections.All combine to support the observed emission of CH 3 OH at 84 GHz toward both MS1 and MS2 and CH 2 NH at 5.29 GHz toward MS2 to consist of astronomical (weak) masers.
As a reference for the relative positions, Figure 1 shows the continuum emission from Sgr B2(N) observed at 6.8 cm (Left) and 3.5 mm (Right).The free-free continuum emission concentrations at 6.8 cm include four cores, K1 -K4, and the two bow-shaped H II regions, K5 and K6 (Corby et al. 2015, marked in green in Figure 1).At 3.5 mm, three main hot cores, N1 -N3, were revealed by their molecular emission (Bonfand et al. 2017, marked in orange in Figure 1).

Class I CH 3 OH Masers
Class I CH 3 OH masers are prevalent in massive star formation regions.In particular, the 44 GHz maser line was observed toward Sgr B2 at 3 ′′ resolution (Mehringer & Menten 1997).In contrast, there is limited information about the 5 −1 − 4 0 E Class I CH 3 OH maser at 84 GHz, a member of the same J +1 −1 −J 0 E-type transition family along with the 36 GHz line (Breen et al. 2019, and references therein).Interferometric observations of the 84 GHz line in the Milky Way are even rarer (Voronkov et al. 2006), as recent work has focused on extragalactic detections (Humire et al. 2022;McCarthy et al. 2018).
Here, we present the first interferometric ALMA observations of the 84 GHz CH 3 OH masers in a Milky Way object at a resolution of 1.5 ′′ , as well as the first detection of the 84 GHz maser toward Sgr B2.The spread of 84 GHz maser spots is manifest in the channel images shown in Figure 2 with velocities ranging from 50 to 86 km s −1 .The brightest maser spot located ∼3.6 ′′ west of the K4 core with a V lsr of 58 km s −1 (hereafter referred to as MS1) has a peak T B of 389.2 ± 0.2 K (8.454 ± 0.004 Jy beam −1 , Figure 2 (q)-(t)).In addition to the brightest spot MS1, another maser knot associated with the K6 H II region is also identified with a T B of 75.4 ± 0.4 K (1.636 ± 0.008 Jy beam −1 ) at a V lsr of 83 km s −1 (hereafter referred to as MS2), as shown in Figure 2 (b)-(e).
The molecular data file for E-type CH 3 OH was taken from the Excitation of Molecules and Atoms for Astrophysics (EMAA) database3 , with spectroscopy from the Cologne Database for Molecular Spectroscopy (CDMS) catalog4 (Müller et al. 2005) and collision rate coefficients reported by Rabli & Flower (2010) for the rotational excitation of CH 3 OH by para-and ortho-H 2 (hereafter p-H 2 and o-H 2 ).The calculations of Rabli & Flower (2010) were performed for rotational transitions within the three first torsional states of CH 3 OH, ν = 0, 1, 2, with a CH 3 OH basis including rotational states j 1 ≤ 15 for p-H 2 (j 2 = 0) and j 1 ≤ 9 for o-H 2 (j 2 = 1), where j 1 and j 2 are the rotational quantum numbers.The coupled-states approximation was used to generate rate coefficients at temperatures 10 ≤ T k ≤ 200 K.Because torsionally inelastic transitions were not considered by Rabli & Flower (2010), the EMAA dataset for E-CH 3 OH includes rotational transitions within the ground torsional state (ν = 0) only.The first excited state (ν = 1) of E-CH 3 OH opening at 204.  3.
Table 3. Summary Statistics of the Marginalized Posterior toward MS1 Notes.The quoted uncertainties represent the 16 th and 84 th percentile (1σ for a Gaussian distribution) uncertainties. (* ) nc represents the number density of the collision partners, which we take to be nH 2 as hydrogen is assumed to be the dominant partner. (†) Two velocity components have similar V lsr of ∼ 58 km s −1 along the LOS. (‡) Continuum temperature for 84-109 GHz band. (§) The MCMC analysis revealed two distinct peaks in the posterior distributions, each representing one of the two scenarios.Depending on the scenario, the background clump exhibits maser or non-inverted emission, while the foreground clump always exhibits maser emission.
the highest level of E-CH 3 OH included in the EMAA dataset is the rotational state 9 −6 (j 1 = 9, K = −6) with an energy of 203.9 cm −1 .

Masers with Dark Continuum Background -MS1
Four uncontaminated E-type CH 3 OH transitions covered in the collision data were used to constrain the physical parameters quantitatively using a non-LTE model and MCMC fitting algorithm.The observed spectra extracted toward the brightest maser spot MS1 are shown in Figure 4. We first attempted with a singleclump model but found that the observed spectra could not be reproduced with physically realistic parameters even after relaxing the constraints on T cont and T rad .In order to match the bright maser emission at 84.5 GHz, strong background radiation is needed to stimulate the maser gain in the single-clump model but it will also in-troduce significant absorption features at both 85.5 and 108.9 GHz.
Moreover, as shown in Figure 2(r), there is no significant mm-wave continuum source at MS1.While this spot also exhibits a strong CH 3 OH emission at 44 GHz (T B = 157 K, Mehringer & Menten 1997), it shows no continuum emission at cm-wavelengths either.The closest continuum source, K4, is offset by ∼ 3.6 ′′ which corresponds to ∼0.14 pc at the distance of Sgr B2 (8.2 kpc, Reid et al. 2019).In addition, as shown in Appendix A, there are no point sources associated with this maser spot.The nearest bright infrared source observed with both the Two Micron All Sky Survey (2MASS) and the Galactic Legacy Infrared Midplane Survey Extraordinaire (GLIMPSE) program is located at ∼ 5.2 ′′ to the west of the maser spot.Instead, the background radiation stimulating the maser gain may be weak and diffuse, causing it to be resolved out by interferometers, such as the free-free and dust emission in the Galactic plane.
As illustrated in Figure 3, the model we consider consists of two chained gas clumps overlapping in velocity along the LOS, where the foreground clump exhibits maser activity but the background clump can either be also an inverted maser (i.e. the double-maser scenario) or non-inverted emitting (i.e. the single-maser scenario).In both scenarios, this model can account for the strong maser emission with an absence of a bright and compact continuum source.It is analogous to the "self-amplification" model of two masing clouds that has been used to interpret water megamasers observed in circumnuclear disks in active galactic nuclei (Kartje et al. 1999).
A chained two-component model has been constructed to iteratively generate synthetic spectra, each component being described with five parameters, T k , n c , N col , dV and V lsr .We assumed the two components have coherent but not identical values of V lsr of ∼58 km s −1 although spatially separated along the LOS.As such, two V lsr parameters have the same prior probability distributions but are fit independently.Including both V lsr values, the priors for the free parameters and the MCMC analysis are discussed in detail in Appendix B. Posterior probability distributions for each parameter and their covariances were generated with 100 walkers with 50000 samples each 5 .
As revealed in the bimodal posterior distributions, the samples converge toward two distinct peaks, each corresponding to the double-maser and single-maser scenarios, respectively.The statistics of the resulting parameters for each scenario were summarized in Table 3, which were used to construct the synthetic spectra for the double-maser scenario shown with the green trace and the single-maser scenario shown with the purple trace in Figure 4.Both synthetic profiles well reproduce the observed spectra, capturing both the maser and the optically thick emission features.Our model, unfortunately, falls short in matching the low-frequency side of the emission plateau at 108.87 GHz, corresponding to the 0 0 − 1 −1 transition at a velocity range of 70 − 90 km s −1 .This divergence from the observed spectrum implies the potential presence of another opticallythick velocity component(s) along the LOS.Because neither of the two scenarios completely reproduces the observed data, it is noteworthy that the uncertainties de-  rived from the MCMC results should be considered as a lower limit of the uncertainties.See Table 4.
Notes.Similar to Table 3. ( * ) Two velocity components are identified along the LOS at a V lsr of 83 km s −1 and 64 km s −1 respectively. (†) Continuum temperature for 84-109 GHz band. (‡) Continuum temperature at 5 GHz.
In the double-maser scenario, the τ and T ex of the 5 −1 − 4 0 maser transition are found to be −2.44 and −5.8 K for the foreground clump and −0.02 and −2027.3K for the background clump.In the singlemaser scenario, the corresponding τ and T ex values are −2.52 and −6.5 K for the foreground clump and 0.13 and 245.2K for the background clump.As an T ex of 245.2 K is significantly higher than the T k value in this case, it is evident that this transition exhibits strong non-inverted emission and is on the verge of a population inversion.

Masers with Bright Continuum Background -MS2
In contrast to MS1, MS2 is aligned with the K6 H II region along the LOS, which provides background radiation to be amplified by the stimulated emission (Figure 2d).The spectra extracted toward MS2 are shown in Figure 6 with two individual velocity components being identified for each of the spectral features.We also performed the spectral analysis for MS2 with a twocomponent model, where each component has a distinct V lsr .The adopted prior and resultant posterior probability distributions for each parameter are depicted respectively in Figure 5.The priors for V lsr were chosen to be Gaussian distributions centering on 83 km s −1 and 63 km s −1 respectively, based on the difference between the observed sky frequency and the rest frequency.Meanwhile, the T cont prior also has a Gaussian distribution, based on the measured continuum level of 6.67 ± 0.52 K.The remaining parameters all have uninformative uniformly distributed priors.Following 100 walkers with 10,000 samples each, all parameters converge to primarily Gaussian-like distributions in our model.The corresponding corner plot of the posterior distributions and the parameter covariances are shown in Appendix B, with the largest covariance observed between H 2 densities and CH 3 OH column densities.This covariance aligns with the fact that H 2 densities and CH 3 OH column densities are related through the maser emission measure, ξ ∝ n c × N col , which, together with temperature, determines the overall rate coefficient for a collisionally-pumped maser transition (Leurini et al. 2016).
The posterior medians from the MCMC inference (Table 4) were used to construct the synthetic spectra shown in Figure 6.If we take the noise level measured in each passband into account, the synthetic spectral profiles fit satisfactorily with the observed spectra for both the maser and optically thick absorption features.However, the low-velocity component (64 km s −1 ) of the 6 −2 − 7 −1 transition fails to match the observed absorption feature at 85.55 GHz.This absorption feature might be attributed to potential contamination from an unidentified transition.The excitation temperature and opacity of the 5 −1 −4 0 maser transition were found to be negative, i.e., −5.6 K of T ex and −2.52 of τ for the highvelocity component (83 km s −1 ), confirming its maser activity toward MS2.

CH 2 NH Masers
The maser activity of the 1 1,0 − 1 1,1 CH 2 NH transitions at 5.29 GHz has been revealed toward both the Galactic Center and external galaxies as point source emission (Faure et al. 2018;Gorski et al. 2021).With the VLA observations, we resolve this emission spatially into three velocity components with a V lsr of 64, 73, and 83 km s −1 .The peak intensity map for each component is shown in Figure 7 respectively.The 83 km s −1 and 64 km s −1 components are both associated with the K6 H II region, while the 73-km s −1 component is concentrated toward the K5 H II region.The brightest 83-km s −1 component displays a more compact distribution than the other two components, with a T B of 121.8 ± 8.6 K (12.1 ± 0.85 mJy beam −1 ).This measured T B represents a lower limit to the peak T B since the observed emission morphology of the 83-km s −1 component may result from point sources convolved with the coarse VLA beam.Nevertheless, this component of the CH 2 NH emission is coherent with that of the 84 GHz CH 3 OH maser emission both in the velocity domain and spatial position.Therefore, the significant brightness, highly compact morphology, and coherence with other maser emission suggest maser activity in the 83-km s −1 component.
The molecular data file for CH 2 NH was taken from the EMAA database.Extended from the calculations in Faure et al. (2018), we calculated the collisional data for rotational transitions of CH 2 NH among the first 49 rotational levels, i.e. up to the state 8 2,6 at 99.5 cm −1 , and for kinetic temperatures in the range 10 -150 K.The detail calculations of the collisional rate coefficients are presented in Appendix C.
We performed the MCMC analysis with a onecomponent model using the prior probability distributions depicted in Figure B4.The VLA observations only cover two transitions of CH 2 NH while the other transition at 4.44 GHz has no signal over the noise levels.As such, the data are insufficient to simultaneously constrain the physical characteristics so the CH 2 NH MCMC analysis requires a confined choice of priors on some parameters.We opted to use E-type CH 3 OH discussed in Section 4.1.2for physical constraints, where the 83 km s −1 components between CH 3 OH and CH 2 NH are spatially coherent and are both associated with the K6 H-II region.We therefore extracted the CH 2 NH spectrum toward the same position of MS2.The posteriors on T k and n H2 obtained from CH 3 OH (Table 4; Figure 5; Figure B3) are adopted as MCMC priors for the CH 2 NH model.The priors for V lsr and dV are chosen to be Gaussian distributions, informed by the strongest emission feature (12σ) detected in the data.The Gaussian prior for T cont is based on the continuum level of 2687 ± 58 K measured by the VLA.Meanwhile, the N col of CH 2 NH has a non-informative uniformly distributed prior.
As shown in the posterior distributions in Figure B4, following 10,000 samples with each of the 100 walkers, the samplers converge at a N col of 4.76 +1.09 −0.64 ×10 14 cm −2 .The synthetic line profile is shown in Figure 8, reproducing six hyperfine components of the 1 1,0 − 1 1,1 transition within the noise level.The six hyperfine components have T ex in the range of -0.75 --0.55 K, τ 's in the range of -0.03 --0.006, and, in particular, the F = 2 − 2 transition has a T ex of -0.65 K and a τ of -0.03.

Metastable Energy Levels of Maser Transitions
The majority of the molecular lines observed toward the H II regions in Sgr B2(N) are in absorption (Corby et al. 2015, and reference therein).In contrast, specific transitions of certain species toward these regions exhibit emission features (e.g.McGuire et al. 2012;Corby et al. 2015;Remijan et al. 2022).These transitions establish population inversions between energy levels (n up g low > n low g up ), a prerequisite for maser action to occur.The mechanism of population inversion is conceptualized by spontaneous decay and collisions: spontaneous decay populates the levels determined by the Einstein coefficients, meanwhile collisions redistribute molecules due to their specific geometry and rotational energy level patterns.Figures 9 and 10 show the energylevel structure of spontaneous transitions connected to the two target maser transitions, with Einstein spontaneous emission coefficients (A ij ) shown in colors.
In the case of the CH 3 OH 5 −1 − 4 0 E transition at 84 GHz, the upper-energy level 5 −1 E is fed by the other higher levels with strong landing transitions but spon- taneously decays to the lower-energy level 4 0 E with a weak escaping transition6 .5 −1 E has a total escaping log 10 A ij /s −1 of −4.2139 (equivalent to a lifetime of ∼ 1.6 × 10 4 s), which is ∼3.0 orders of magnitude lower than its total landing A ij .As discussed in Lees (1973) and Menten (1991), CH 3 OH spontaneously decays more rapidly within the K = −1 stack as compared with decaying out of the K = −1 stack.Specifically, among the two escaping transitions, 5 −1 − 4 −1 E is stronger than 5 −1 − 4 0 E (Figures 9).In fact, the K = −1 levels are the end-points for all spontaneous decay routes (Lees 1973).For example, the lower-energy level of this maser transition (4 0 E) decays rapidly to the lower K = −1 levels.
For CH 2 NH, we sum the hyperfine components of the 14 N splittings and use the sum of A ij as the representative A ij for each transition in Figure 10.In the case of the CH 2 NH 1 1,0 − 1 1,1 transition at 5.29 GHz, the upper-energy level 1 1,0 is coupled to the higher levels 2 1,1 and 2 2,1 with strong spontaneous transitions but to the lower-energy level 1 1,1 with a weak escaping transition.Compared with the upper-energy level 1 1,0 , which has a total escaping log 10 A ij /s −1 of −4.2476, (equivalent to a lifetime of ∼ 1.8 × 10 4 s), the lower-energy level 1 1,1 has a slightly stronger escaping log 10 A ij /s −1 of −4.0303 (equivalent to a shorter lifetime of ∼ 1.1 × 10 4 s).
Meanwhile, collisions of CH 2 NH exhibit a propensity for the |∆K a | = 0 transitions compared with the inter-stack transitions |∆K a | = 1 and 2 (Faure et al. 2018).In particular, the total rate coefficient of the ∆K a = 0 collisions de-populating a K a = 1 level (∼ 10 −10 cm 3 s −1 ) is one order of magnitude higher than that of the |∆K a | = 1 and 2 collisions at a T k of 30 K.However, the inter-stack transitions from the ground state 0 0,0 complicate this multi-level system.The collisions from 0 0,0 more efficiently populate the lower level 1 1,1 than the upper level 1 1,0 , suppressing the inversion.Therefore, while these energy structures allow the temporary trapping of the population in the metastable 1 1,0 level, it only results in a weak population inversion.

Spatial Origin and Pumping Mechanism
The maser activity of several simple species was reported previously toward Sgr B2(N) with precise locations, including CH 3 OH, H 2 O, OH, and NH 3 , which are indicated by + symbols in Figures 1 and 11.The Class II CH 3 OH maser spot at 6.7 GHz is spatially associated with the N3 mm-wave core (Caswell 1996;Caswell et al. 2010).A survey performed with the Australia Telescope Compact Array (ATCA) characterized 30 H 2 O maser spots with the majority of them located within the N1-K2-K3 region (Walsh et al. 2014).Six OH maser spots have been identified with a VLA survey of OH masers with the brightest spot being associated with the K3 cmwave core (Argon et al. 2000;Qiao et al. 2014).In a recent VLA observation, nine NH 3 (9,6) maser spots were resolved in this region; most of these spots are concentrated in the N1-K2-K3 region, while the brightest spot is located toward the N3 hot molecular core (Yan et al. 2022).In contrast, the maser emission of the organic molecular species presented in this work is found to be located close to the K4 cm-wave core and the K6 H II region, being spatially offset from the spots of the radiative pumping maser species OH and Class II CH 3 OH.The 84 GHz CH 3 OH maser displays a bright but extended distribution, centered around MS1 (near the K4 core), and subtending an angular size of ∼ 15 ′′ .The emission tapers off rapidly from the center with non-Gaussian profiles, a morphology possibly attributed to maser beaming effects.
We found that MS2, superimposed on the K6 H II region, harbors maser emission of both Class I CH 3 OH at 84 GHz and CH 2 NH at 5.29 GHz (Figure 11).The 83-km s −1 component of the CH 2 NH maser at K6 remarkably resembles that of the Class I CH 3 OH maser.Their consistent morphology serves as strong evidence to support a similar collisional pumping mechanism for the CH 2 NH maser with the Class I CH 3 OH maser as well as a common physical condition for triggering these masers.By understanding the conditions that are necessary for maser activities to occur, we can also better understand how they can be used to monitor astronomical objects.For example, since these masers are pumped by intense collisions, we can expect them to be good candidates for studying regions with time-dependent gas dynamics, such as episodic outflows from massive protostars.Future surveys of these new masers in addition to such objects are needed to determine how frequently they accompany Class I CH 3 OH masers and how they might fit into proposed evolutionary timelines of masers in massive star formation (e.g., Ellingsen et al. 2012).
Alongside our observations of the Class I 84 GHz CH 3 OH maser, we also examined the CH 3 OH 8 0 −7 1 Atype transition at 95 GHz, which exhibits maser activity toward several famous star-forming regions, e.g., Orion-KL (Plambeck & Wright 1988), DR 21 (Plambeck & Menten 1990), and W51 (Pratap & Menten 1993).A recent single-dish survey of ATLASGAL clumps shows that the 84 GHz line is generally weaker than the 95 GHz The 83-km s −1 components of the 84 GHz Class I CH3OH maser and 5.29 GHz CH2NH maser are shown with contours in pink and orange respectively, and superimposed on the free-free continuum emission at 6.8 cm.Yellow shading indicates the overlap region where the CH3OH emission exceeds 40 σ and the CH2NH emission exceed 4 σ.The 84 GHz Class I CH3OH maser spots, MS1 and MS2, are marked with + symbols in megenta while the brightest maser spots of OH, H2O, NH3, and Class II CH3OH masers are marked in other colors (Argon et al. 2000;Caswell et al. 2010;Walsh et al. 2014;Yan et al. 2022).The K4 core and arc-like K6 H II region are marked with a green circle maker and dotted line respectively.line (Yang et al. 2023).In contrast, we found that the 95 GHz line exhibits weaker emission toward the identified maser spots in Sgr B2(N) compared to the 84 GHz line.The spatial distribution for the 95 GHz emission is similar to that for other thermal molecular lines, concentrating toward the N1 mm-wave source.This distribution may be attributed to either a weak maser or thermal radiation at different emission spots.It is noteworthy that Mehringer & Menten (1997) identified both maser and quasi-thermal emission components toward Sgr B2 in the 44 GHz line, a consecutive transition (J = 6) of the (J + 1) 0 − J 1 A series with the 95 GHz transition (J = 7).A more detailed analysis of other A-type CH 3 OH transitions is needed to determine the nature of the 95 GHz emission.

Physical Conditions of Masing Clumps
The high spatial resolutions of interferometric observations enable us to resolve and segregate compact molecular emissions in diverse environments.In particular, the rich thermalized emission of CH 3 OH and CH 2 NH occurs in the hot core regions.In the N1 hot core, the ALMA observations reported an N col of 2.0 × 10 19 cm −2 for CH 3 OH and 9.0 × 10 17 cm −2 for CH 2 NH at a T ex of 230 K (Motiyenko et al. 2020;Margulès et al. 2022).In contrast, the bright stimulated maser emission of these two molecules concentrates in the less dense masing clumps.The N col of these two molecules in the masing clumps (Table 3-4) are 2 − 3 orders of magnitudes less abundant than those in the hot core.
From the MCMC analyses on the CH 3 OH 84 GHz maser presented in Section 4.1, we found the T k for the masing clumps to be ∼29 − 37 K.This range of T k is consistent with the large-scale measurement of the Sgr B2 envelope from the OCS lines (30-37 K, Armijos-Abendaño et al. 2020) as well as the dust temperature of Sgr B2(N) obtained from Herschel observations (∼ 28 K, Etxaluze et al. 2013).It is also consistent with the value of T k (30 K) derived from the non-LTE modeling of the weak masers of methyl formate detected in PRIMOS (Faure et al. 2014).
We also constrain the n H2 of the masing clumps to be (5 − 8) × 10 4 cm −3 , a level significantly lower than that derived for the hot core regions.The n H2 of the Sgr B2(N) hot cores are > 10 7 cm −3 , with the N1 core reaching ∼ (1 − 2) × 10 9 cm −3 (Sánchez-Monge et al. 2017).In contrast, the n H2 of the masing clumps is comparable to the n H2 of the colder envelope of Sgr B2, which is estimated to be 10 3 − 10 4 cm −3 based on the absorption molecular lines (Huettemeister et al. 1995;Schmiedeke et al. 2016).The slight enhancement in physical densities might be attributed to local energetic events occurring in the envelope necessary for a collisionally pumped maser.
This work presents a quantitative constraint on the physical conditions characterizing a masing clump of organic molecules.It serves as a reference for the physical conditions sufficient for exciting population inversion and collisionally-pumped masers.These conditions include a T k of ∼ 30, n H2 > 5 × 10 4 cm −3 , and moderate molecular N col (> 10 16 cm −2 for CH 3 OH and > 10 14 cm −2 for CH 2 NH).The relatively low values of temperatures and densities indicate that the masing clumps are situated in the colder and low-density external envelope.The distance from Sgr B2(N) to the edge of the external envelope is suggested to be approximately 19 pc (Schmiedeke et al. 2016).The continuum emission from Sgr B2(N) has an angular size of ∼ 10 ′′ , or 0.4 pc at a distance of 8.34 kpc (Reid et al. 2014), which subtends a solid angle of 3.5 × 10 −4 sr to the envelope.This aligns with our assumption mentioned in Section 3.1, where background continuum emission sources are far from masing clumps and, therefore, do not significantly affect the local energy distribution of the molecules.Maser activities can manifest themselves in the presence of either a bright continuum background radiation or multiple chained masing clumps in the direction of a diffuse and weak background.

Caveat on Geometric Effects
Geometric effects related to the propagation of maser radiation become increasingly significant as the optical depth becomes lower than -1.Although RADEX and MOLSIM support various geometries through escape probability approximation, other effects such as maser beaming are not treated.Such effects are particularly important in the case of unresolved clumps.The models also assume homogeneous physical conditions (i.e.one-zone models), which is inappropriate for very negative optical depths.For precisely modeling maser lines with strong amplification, a non-local treatment of radiation and level population is required.
Although our current treatment does not provide a precise solution, it is in principle appropriate for | τ |≪ 1 (the case of CH 2 NH) and still offers a reasonable approximation to the first order as | τ | is of the order of unity (the case of CH 3 OH).Incorporating corrections such as non-local treatments and other geometric effects would definitely improve the solution but it would also introduce significant computational costs.This becomes especially infeasible for MCMC calculations.The current solution represents the best that is achievable with our current resources, acknowledging that this approach is subject to future expansions and refinements.

SUMMARY
We report interferometric observations of maser activity from organic molecules toward Sgr B2 using the ALMA and VLA.Enabled by the ALMA observations, we report the first detection of the Class I CH 3 OH maser at 84 GHz toward Sgr B2.Multiple 84 GHz CH 3 OH maser knots are discovered with the brightest spot MS1 located ∼3 ′′ west of the K4 free-free continuum emission core.The peak brightness temperature of MS1 reaches 389.2 ± 0.2 K at a velocity of 58 km s −1 .In addition, another CH 3 OH maser spot MS2 coexists with the K6 UC H II region at a velocity of 83 km s −1 .With the VLA observations, we resolve the 5.29 GHz CH 2 NH maser emission into three velocity components.These three components display distinct fine-scale structures associated with the K5 and K6 bow-shaped H II regions.The 83 km s −1 component is co-spatial with MS2 and its peak brightness temperature has a lower limit of 121.8 ± 8.6 K.
We present a tool for modeling and fitting maser spectra with non-LTE calculations and a Bayesian ap-proach.By leveraging the theoretical calculations on collisional rate coefficients and cross sections for both Etype CH 3 OH and CH 2 NH transitions, this fitting technique enables us to statistically constrain the physical conditions of the masing regions.In particular, our analysis suggests a two-chain-clump maser model for explaining the intense CH 3 OH maser activity observed toward a region with weak and diffuse background radiation.This model can reproduce both the maser and the optically thick emission features of the CH 3 OH spectra observed toward the brightest maser spot MS1.
Interferometric observations have not only enabled us to characterize the spatial origin and extent of maser emission but also helped to constrain the pumping mechanism through direct comparisons of their morphologies.We found a definite association between the maser emission from complex molecules and the K6 UC H II region.They are offset from the masing regions of the radiatively pumped maser species and the infrared sources.In contrast, the spatial correlation between the activities of the 84 GHz CH 3 OH Class I maser and 5.29 GHz CH 2 NH maser suggest that the CH 2 NH maser is pumped by intense collisions, analogous to the Class I CH 3 OH maser.By understanding the nature of the new maser species, which are pumped by intense collisions, we can better understand how they may serve as good candidates for studying and monitoring astronomical objects with time-dependent gas dynamics.

DATA ACCESS & CODE
The code used to perform the analysis is part of the molsim open-source package; an archival version of the code can be accessed at Lee et al. (2023).
The spectroscopic and collisional  component, corresponding to the higher velocity peak of the background component, corresponds to the double-maser scenario, while the other peak corresponds to the single-maser scenario.To better visualize the posteriors, we isolated the two peaks using V lsr = 58 km s −1 as the criterion and plotted the posteriors separately in Figure B2.The corner plots with green and purple backgrounds correspond to the double-maser and single-maser scenarios, respectively.
In each corner, the first five parameters belong to the foreground component, while the last five are those of the background component.The 2D posterior distribution among parameters for the foreground component is shown in blue, while that for the background component is shown in red.The cross terms between the two components are shown in brown.The middle parameter is the background continuum temperature and the distribution between continuum temperature and other parameters is shown in yellow.
A similar MCMC analysis to that for the spectra extracted toward MS1 was carried out for the spectra extracted toward MS2.In contrast to Figure B1, the optimal set of parameters of MS2 converges to a single peak as shown in Figure 5.The corresponding corner plot of the parameter covariances and their distribution for the MS2 MCMC fit are shown in Figure B3.Similar to Figure B2, the parameters of the low-velocity component are shown in blue while those of the high-velocity component are shown in red.We found that the MCMC analyses yield the consistent results, regardless of whether the foreground or background velocity component is assigned as the high or low velocity along the LOS.Furthermore, the brown circular distribution, the cross terms between the two components, reveals that the two components are completely independent.This independence arises because the two components fit distinct portions of the spectrum because of their disparate central velocities compared with the line widths.In contrast, the strong degeneracy of the H 2 density and CH 3 OH column density is revealed for each individual velocity component.
CH2NH -An MCMC analysis with a one-clump model was carried out for the CH 2 NH spectra.The prior and posterior distributions are shown in Figure B4, while a corner plot of the parameter covariances for the CH 2 NH MCMC fit is shown in Figure B5.

C. COLLISIONAL RATE COEFFICIENT CALCULATIONS FOR CH 2 NH
A first set of collisional data for CH 2 NH + p-H 2 (j 2 = 0) was published by Faure et al. (2018) for applications to cold gas.The calculations of Faure et al. (2018) were performed for rotational transitions among the first 15 rotational levels of CH 2 NH, i.e. up to the state 2 2,0 (j 1 = 2, K a = 2, K c = 0) with an energy of 28.3 cm −1 .The full close-coupling method was used to generate rate coefficients at temperatures 5 ≤ T k ≤ 30 K and hyperfine-resolved rate coefficients were derived from the rotational rate coefficients using the statistical approximation, see Faure et al. (2018) for full details.In this work, we have extended the calculations of Faure et al. (2018) to higher levels and temperatures and we have thus considered o-H 2 (j 2 = 1) as a second collider.Because close-coupling calculations are computationally very demanding when the temperature increases, we had recourse to the coupled-states approximation.As in Faure et al. (2018)  Green 2012), using the same propagation parameters and rotational basis for CH 2 NH (j 1 ≤ 12) and p-H 2 (j 2 = 0, 2).The basis for o-H 2 was restricted to j 2 = 1 and the MOLSCAT parameter EMAX (used to limit basis functions to CH 2 NH plus H 2 energies less than EMAX) was set to 1200 cm −1 for p-H 2 and 1318.64 cm −1 for o-H 2 .The coupled-states calculations were performed for total energies between 2 and 800 cm −1 for p-H 2 and between 120 and 918 cm −1 for o-H 2 .Rate coefficients were obtained for all rotational transitions among the first 49 levels, i.e. up to the state 8 2,6 at 99.5 cm −1 , and for kinetic temperatures in the range 10 -150 K. Hyperfine-resolved rate coefficients were derived using the statistical approximation, as in Faure et al. (2018).Finally, the collisional dataset was combined with the energy levels and radiative rates taken from the CDMS catalog (Müller et al. 2005) in order to provide an EMAA-formatted datafile.

Figure 1 .
Figure 1.Continuum emission toward the Sgr B2(N) region obtained with the VLA at 6.8 cm and ALMA at 3.5 mm.The cm-wave emission concentrations (K1 -K6, Corby et al. 2015) are indicated by green dashed lines and circle markers whereas the molecular hot cores identified in millimeter (N1 -N3, Bonfand et al. 2017) are indicated by orange triangle markers.The + symbols mark the location of Class II CH3OH masers in blue, the brightest maser spots are marked in rose for the H2O maser(Walsh et al. 2014), in purple for the OH maser(Argon et al. 2000), and in firebrick for the NH3 maser(Yan et al. 2022).

Figure 4 .
Figure 4. Observed spectra of the E-type CH3OH transitions extracted toward MS1 are shown in black with noise levels shaded in gray.Quantum numbers are given at the top right of each panel.The synthetic spectra are overlaid in colors for different scenarios, with the green trace representing the double-maser scenario and the purple trace representing the singlemaser scenario.See Table3.

Figure 5 .
Figure 5. Prior and posterior distributions for the MCMC fit of CH3OH E-type transitions toward MS2.The prior distribution of each parameter is shown in gray.The posterior distributions for the low-velocity component are shown in blue, whereas those for the high-velocity component are shown in red.The posterior of Tcont is shown in yellow.The discrepancy of Tcont in the posterior relative to the prior could be be attributed to a certain portion of continuum flux, especially from extended components, being resolved out by ALMA.

Figure 6 .
Figure 6.Observed spectra of the E-type CH3OH transitions extracted toward MS2 are shown in black with noise levels shaded in gray.Quantum numbers are given at the top right of each panel.The rest frequency for each transition are 84423.769,84521.172,85568.131,108893.945MHz from left to right.The synthetic spectra are overlaid in colors with the red trace representing the high-velocity component (83 km s −1 ) and the blue trace representing the low-velocity component (64 km s −1 ).See Table4.

Figure 7 .Figure 8 .
Figure 7. Peak intensity maps for each velocity component of the CH2NH maser at 5.29 GHz superimposed on the continuum emission observed at 6.8 cm.The contour levels start at 4 σ and continue up to the peak value with a step size of 2 σ.The values of σ and peak values for the 83-, 73-, and 64-km s −1 components are 0.80, 0.93, 0.85 mJy beam −1 , and 12.1, 8.86, 6.60 mJy beam −1 , respectively.

Figure 9 .
Figure 9. Energy level structure of the connected rotational transitions for the CH3OH 5−1 − 40 E transition.Energy levels are classified by the K quantum number and ordered by increasing energy but are not drawn to scale.The log 10 A ij s −1 for each transition is shown with diverging color schemes, in which the critical middle value is set to be the log 10 A ij s −1 of the 5−1 − 40 E transition (−5.7).The targeted transitions are shown with gray arrows.Transitions with a smaller Aij are shown as blue dotted lines while transitions with a larger Aij are shown as red solid lines.
Figure 11.The 83-km s −1 components of the 84 GHz Class I CH3OH maser and 5.29 GHz CH2NH maser are shown with contours in pink and orange respectively, and superimposed on the free-free continuum emission at 6.8 cm.Yellow shading indicates the overlap region where the CH3OH emission exceeds 40 σ and the CH2NH emission exceed 4 σ.The 84 GHz Class I CH3OH maser spots, MS1 and MS2, are marked with + symbols in megenta while the brightest maser spots of OH, H2O, NH3, and Class II CH3OH masers are marked in other colors(Argon et al. 2000;Caswell et al. 2010;Walsh et al. 2014;Yan et al. 2022).The K4 core and arc-like K6 H II region are marked with a green circle maker and dotted line respectively.

Figure B1 .
Figure B1.Prior and posterior distributions for the MCMC fit of CH3OH E-type transitions toward MS1.The prior distribution of each parameter is shown in gray.The bimodal posterior distributions for the foreground component are shown in blue, whereas those for the background component are shown in red.The posterior of Tcont is shown in yellow.

Figure B2 .
Figure B2.Parameter covariances and marginalized posterior distributions for the MCMC fit of CH3OH E-type transitions toward MS1.The 16 th , 50 th , and 84 th confidence intervals (corresponding to ±1 sigma for a Gaussian posterior distribution) are shown as vertical lines.The distribution among the parameters for the foreground component is shown in blue and that for the background component is shown in red.The distribution among the parameters between the two components is shown in brown.The distribution between the background continuum temperature and other parameters is shown in yellow.To better visualize the bimodal posteriors, we isolated the two peaks and plotted separate posteriors for each, with the double-maser and single-maser scenarios highlighted by green and purple backgrounds respectively, as per the color scheme in Figure4.

Figure B3 .
Figure B3.Similar to Figure B2, parameter covariances and marginalized posterior distributions for the MCMC fit of CH3OH E-type transitions toward MS2.The distribution among the parameters for the low-velocity (64-km s −1 ) component are shown in blue and that for the high-velocity (83-km s −1 ) component is shown in red.The distribution among the parameters between the two components are shown in brown.The distribution between background continuum temperature and other parameters is shown in yellow.

Figure B4 .Figure B5 .
Figure B4.Prior and posterior distributions for the MCMC fit of the CH2NH transitions toward MS2.The prior distribution for each parameter is shown in grey, while the priors on T k and nH 2 are from the MCMC posteriors of CH3OH in MS2 (Figure 5).The posterior distributions for the physical characteristics are shown in red while the posterior for Tcont is shown in yellow.

Table 1 .
Summary of VLA Observation Parameters

Table 2 .
Investigated Peak brightness temperature measured at the brightest emission spot for each molecule.The brightest emission spot for 84 GHz CH3OH Class I maser is located near the K4 core at a V lsr of 58 km s −1 .The brightest emission spot for 5.29 GHz CH2NH maser is located in the K6 H II region at a V lsr of 83 km s −1 .The measured TB of CH2NH represents a lower limit to the T B,peak because of the limited angular resolution of the VLA observations.

5
Each walker represents a Markov Chain consisting of multiple samples with each sample being randomly sampled within the parameter space.
, all calculations were performed with the OpenMP extension 7 of the version 14 of the MOLSCAT code(Hutson &