The 6.7 GHz methanol maser emission, a tracer of forming massive stars, sometimes shows enigmatic periodic flux variations over several 10-100 days. In this Letter, we propose that these periodic variations could be explained by the pulsation of massive protostars growing under rapid mass accretion with rates ofṀ * 10 −3 M yr −1 . Our stellar evolution calculations predict that the massive protostars have very large radii exceeding 100 R at maximum, and here we study the pulsational stability of such bloated protostars by way of the linear stability analysis. We show that the protostar becomes pulsationally unstable with various periods of several 10-100 days depending on different accretion rates. With the fact that the stellar luminosity when the star is pulsationally unstable also depends on the accretion rate, we derive the period-luminosity relation log(L/ L ) = 4.62 + 0.98 log(P /100 days), which is testable with future observations. Our models further show that the radius and mass of the pulsating massive protostar should also depend on the period. It would be possible to infer such protostellar properties and the accretion rate with the observed period. Measuring the maser periods enables a direct diagnosis of the structure of accreting massive protostars, which are deeply embedded in dense gas and are inaccessible with other observations.

No abstract

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 3 OH maser emission. Using SMA and ALMA, we have discovered 14 new Class II CH 3 OH maser lines ranging in frequency from 199 to 361 GHz, which originate mostly from v t =1 torsionally-excited transitions and include one v t =2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the v t =2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a con-

Context. Class II methanol masers are signposts of massive young stellar objects (MYSOs). Recent evidence shows that flares of these masers are driven by MYSO accretion bursts. Thus, maser monitoring can be used to identify such bursts which are hard to discover otherwise. Infrared observations reveal burst-induced changes in the spectral energy distribution (first and foremost a luminosity increase), which provide valuable information on a very intense phase of high-mass star formation. Aims. In mid-January 2019, flaring of the 6.7 GHz CH3OH maser (hereafter maser) of the MYSO G358.93-0.03 (hereafter G358) was reported. The international maser community initiated an extensive observational campaign which revealed extraordinary maser activity and yielded the detection of numerous new masering transitions. Interferometric imaging with the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array resolved the maser emitting core of the star forming region and proved the association of the masers with the brightest continuum source (MM1), which hosts a hot molecular core. These observations, however, failed to detect a significant rise in the (sub)millimeter dust continuum emission. Therefore, we performed near-infrared (NIR) and far-infrared (FIR) observations to prove or disprove whether the CH3OH flare was driven by an accretion burst. Methods. NIR imaging with the Gamma-Ray Burst Optical/Near-infrared Detector has been acquired and integral-field spectroscopy with the Field-Imaging Far-Infrared Line Spectrometer (FIFI-LS) aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA) was carried out on two occasions to detect possible counterparts to the (sub)millimeter sources and compare their photometry to archival measurements. The comparison of pre-burst and burst spectral energy distributions is of crucial importance to judge whether a substantial luminosity increase, caused by an accretion burst, is present and if it triggered the maser flare. Radiative transfer modeling of the spectral energy distribution (SED) of the dust continuum emission at multiple epochs provides valuable information on the bursting MYSO. Results. The FIR fluxes of MM1 measured with FIFI-LS exceed those from Herschel significantly, which clearly confirms the presence of an accretion burst. The second epoch data, taken about 16 months later, still show increased fluxes. Our radiative transfer modeling yielded major burst parameters and suggests that the MYSO features a circumstellar disk which might be transient. From the pre-burst, burst, and post-burst SEDs, conclusions on heating and cooling time-scales could be drawn. Circumstances of the burst-induced maser relocation have been explored. Conclusions. The verification of the accretion burst from G358 is another confirmation that Class II methanol maser flares represent an alert for such events. Thus, monitoring of these masers greatly enhances the chances of identifying MYSOs during periods of intense growth. The few events known to date already indicate that there is a broad range in burst strength and duration as well as environmental characteristics. The G358 event is the shortest and least luminous accretion burst known to date. According to models, bursts of this kind occur most often.