We present (sub)millimeter imaging at 0.″5 resolution of the massive star-forming region G358.93−0.03 acquired in multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz CH₃OH maser emission. Using the Submillimeter Array and Atacama Large Millimeter/submillimeter Array, we have discovered 14 new Class II CH₃OH maser lines ranging in frequency from 199 to 361 GHz, which originate mostly from = 1 torsionally excited transitions and include one = 2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the = 2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure 1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral profiles evolved significantly over one month, and the intensities dropped by factors of 3.0–7.2, with the = 2 line showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1 has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core (MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of the sources during the one month interval, and the total protocluster emission remains comparable to prior single-dish measurements.

We report the detection of new 12.178, 12.229, 20.347, and 23.121 GHz methanol masers in the massive star-forming region G358.93−0.03, which are flaring on similarly short time-scales (days) as the 6.668 GHz methanol masers also associated with this source. The brightest 12.178 GHz channel increased by a factor of over 700 in just 50 d. The masers found in the 12.229 and 20.347 GHz methanol transitions are the first ever reported and this is only the fourth object to exhibit associated 23.121 GHz methanol masers. The 12.178 GHz methanol maser emission appears to have a higher flux density than that of the 6.668 GHz emission, which is unusual. No associated near-infrared flare counterpart was found, suggesting that the energy source of the flare is deeply embedded.

We report the first detection of isotopic methanol (¹³CH₃OH) maser emission in interstellar space. The emission was detected toward the high-mass young stellar object G358.93-0.03 during monitoring of a flare in the 6.7 GHz methanol (CH₃OH) maser emission in this source. We find that the spectral and spatial distribution of the ¹³CH₃OH masers differs from the CH₃OH masers imaged at the same epoch, contrary to expectations from similarity of their pumping. This conclusively demonstrates that isotopic methanol masers are bright under different physical conditions and suggests that they can provide additional, complementary information to the CH₃OH masers from the same source. We detect a rapid decay of the ¹³CH₃OH maser lines suggesting that they are transient phenomena (masing for only a few months), likely associated with rapid changes in radiation field due to an accretion burst induced by massive disk fragmentation. Changes in the line flux density are faster than required to achieve equilibrium in the energy level populations, indicating that the pumping of these masers is likely variable.

In 2019 September, a sudden flare of the 6.7 GHz methanol maser was observed toward the high-mass young stellar object (HMYSO) G24.33+0.14. This may represent the fourth detection of a transient mass accretion event in an HMYSO after S255IR NIRS3, NGC 6334I-MM1, and G358.93−0.03-MM1. G24.33+0.14 is unique among these sources as it clearly shows a repeating flare with an 8 yr interval. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the millimeter continuum and molecular lines toward G24.33+0.14 in the pre-flare phase in 2016 August (ALMA Cycle 3) and the mid-flare phase in 2019 September (ALMA Cycle 6). We identified three continuum sources in G24.33+0.14, and the brightest source, C1, which is closely associated with the 6.7 GHz maser emission, shows only a marginal increase in flux density with a flux ratio (Cycle 6|$/$|Cycle 3) of 1.16 ± 0.01, considering an additional absolute flux calibration uncertainty of |$10\%$|⁠. We identified 26 transitions from 13 molecular species other than methanol, and they exhibit similar levels of flux differences with an average flux ratio of 1.12 ± 0.15. In contrast, eight methanol lines observed in Cycle 6 are brighter than those in Cycle 3 with an average flux ratio of 1.23 ± 0.13, and the higher excitation lines tend to show a larger flux increase. If this systematic increasing trend is real, it would suggest radiative heating close to the central HMYSO due to an accretion event which could expand the size of the emission region and/or change the excitation conditions. Given the low brightness temperatures and small flux changes, most of the methanol emission is likely to be predominantly thermal, except for the 229.759 GHz (8₋₁–7₀ E) line known as a class I methanol maser. The flux change in the millimeter continuum of G24.33+0.14 is smaller than in S255IR NIRS3 and NGC 6334I-MM1 but is comparable with that in G358.93−0.03-MM1, suggesting different amounts of accreted mass in these events.