We discuss the role that coherence phenomena can have on the intensity variability of spectral lines associated with maser radiation. We do so by introducing the fundamental cooperative radiation phenomenon of (Dicke’s) superradiance and discuss its complementary nature to the maser action, as well as its role in the flaring behaviour of some maser sources. We will consider examples of observational diagnostics that can help discriminate between the two, and identify superradiance as the source of the latter. More precisely, we show how superradiance readily accounts for the different time-scales observed in the multi-wavelength monitoring of the periodic flaring in G9.62+0.20E.

The Maser Monitoring Organisation is a collection of researchers exploring the use of time-variable maser emission in the investigation of astrophysical phenomena. The forward directed aspects of research primarily involve using maser emission as a tool to investigate star formation. Simultaneously, these activities have deepened knowledge of maser emission itself in addition to uncovering previously unknown maser transitions. Thus a feedback loop is created where both the knowledge of astrophysical phenomena and the utilised tools of investigation themselves are iteratively sharpened. The project goals are open-ended and constantly evolving, however, the reliance on radio observatory maser monitoring campaigns persists as the fundamental enabler of research activities within the group.

Recently, remarkable progress has been made in understanding the formation of high mass stars. Observations provided direct evidence that massive young stellar objects (MYSOs), analogously to low-mass ones, form via disk-mediated accretion accompanied by episodic accretion bursts, possibly caused by disk fragmentation. In the case of MYSOs, the mechanism theoretically provides a means to overcome radiation pressure, but in practice it is poorly studied - only three accretion bursts in MYSOs have been caught in action to date. A significant contribution to the development of the theory has been made with the study of masers, which have proven to be a powerful tool for locating “bursting” MYSOs. This overview focuses on the exceptional role that masers play in the search and study of accretion bursts in massive protostars.

Long-term monitoring observations of the 6.7 GHz methanol masers by Hitachi 32-m operated by Ibaraki University, which are named as “the Ibaraki 6.7 GHz Methanol Maser Monitor (iMet)”, have revealed that the periods of the flux variability of 6.7 GHz methanol masers in the five high-mass star-forming regions G05.900-0.430, G06.795-0.257, G10.472+0.027, G12.209-0.102 and G13.657-0.599 are over 1000 days. These periods are approximately twice the longest known period of 6.7 GHz methanol masers of 668 days for G196.45-1.68. The facts that the flux variation patterns show symmetric sine curves and that the luminosity of the central protostar and periods of maser flux variation are consistent with the expected period-luminosity (PL) relation suggest that the mechanism of maser flux variability of G05.900-0.430, G10.472+0.027 and G12.209-0.102 can be explained by protostellar pulsation instability. From the period-luminosity relation, central stars of these three sources are expected to be very high-mass protostars with a mass of : ∼40 M⊚ and a mass accretion rate of : ∼2×10-2 M⊚yr-1. On the other hand, G06.795-0.257 and G13.657-0.599 have luminosities that are an order of magnitude smaller than that expected from PL relation, and the variation patterns are intermittent, suggesting a variation mechanism of these sources originated from binary system.

We report VLBI monitoring observations of the 22 GHz H2O masers toward the Mira variable BX Cam. Data from 37 epochs spanning ∼3 stellar pulsation periods were obtained between May 2018 and June 2021 with a time interval of 3–4 weeks. In particular, the VERA dual-beam system was used to measure the kinematics and parallaxes of the H2O maser features. The obtained parallax, 1.79±0.08 mas, is consistent with Gaia EDR3 and previous VLBI measurements. The position of the central star was estimated relied on Gaia EDR3 data and the center position of the 43 GHz SiO maser ring imaged with KVN. Analysis of the 3D maser kinematics revealed an expanding circumstellar envelope with a velocity of 13±4 km s−1 and significant spatial and velocity asymmetries. The H2O maser animation achieved by our dense monitoring program manifests the propagation of shock waves in the circumstellar envelope of BX Cam.

We started high-cadence monitoring observations of 6.7 GHz methanol masers from Dec. 2012 using Hitachi 32-m radio telescope (Yonekura et al.(2016). Observations have been conducted basically every day. On average, 13 hours of observations have been made per day, amounting to 4000–5000 hours per year. The cadence varies by sources: one observation in 1–50 days. In addition to already known 29 sources (Tanabe et al. 2023 and references therein), we have newly identified ∼20 sources with periodic flux variability. We have also detected 5 sources with sudden flux rises in 2019–2022, including G358.93–0.03 which was confirmed to be associated with the accretion burst.

We report on a direct comparison of VLBI maser data and ALMA thermal-emission data for the high-mass protostar G353.273+0.641. We detected a gravitationally-unstable disk by dust and a high-velocity jet traced by a thermal CO line by ALMA long-baselines (LB). 6.7 GHz CH3OH masers trace infalling streamlines inside the disk. The innermost maser ring indicates another compact accretion disk of 30 au. Such a nested system could be caused by angular momentum transfer by the spiral arms. 22 GHz H2O masers trace the jet-accelerating region, which are directly connecting the CO jet and the protostar. The recurrent maser flares imply episodic jet ejections per 1–2 yr, while typical separation of CO knots indicates a variation of outflow rate per 100 yr. Our study demonstrates that VLBI maser observations are still a powerful tool to explore detailed structures nearby high-mass protostars by combining ALMA LB.

We report on interferometric observations of a face-on accretion system around the high mass young stellar object, G353.273+0.641. The innermost accretion system of 100-au radius was resolved in a 45-GHz continuum image taken with the Jansky Very Large Array. Our SED analysis indicated that the continuum could be explained by optically-thick dust emission. 6.7 GHz CH3OH masers associated with the same system were also observed with the Australia Telescope Compact Array. The masers showed a spiral-like, non-axisymmetric distribution with a systematic velocity gradient. The line-of-sight velocity field is explained by an infall motion along a parabolic streamline that falls onto the equatorial plane of the face-on system. The streamline is quasi-radial and reaches the equatorial plane at a radius of 16 au. The physical origin of such a streamline is still an open question and will be constrained by the higher-resolution thermal continuum and line observations with ALMA long baselines.

For the high-sensitivity 6.7 GHz methanol maser observations, we developed a new technology for coherently combining the two signals from the Hitachi 32 m radio telescope and the Takahagi 32 m radio telescope of the Japanese VLBI Network. Furthermore, we compared the SNRs of the 6.7 GHz maser spectra for two methods. One is a VLBI method and the other is the newly developed digital position switching, which is a similar technology to that used in noise-cancelling headphones. We report the phase-up technique and the observation.

We performed 12CO(J = 2 − 1) and 13CO(J = 2 − 1) line surveys of the galactic plane using the AMANOGAWA telescope (Nakajima et al. 2007) which covers 5° ≤ l ≤ 180° on a 7'.5 grid. The telescope beamsize and velocity resolution are 8.7′ ± 0.4′ and 1.3 km s−1. The resultant rms noise level is typically 0.12 K in TA*. We found a linear correlation between 12CO(J = 2 − 1) and 12CO(J = 1 − 0) and a curved correlation between 12CO(J = 2 − 1) and 13CO (J = 2 − 1), although the intensity ratios of these three lines have intrinsic variation. These correlations can be reproduced with simple radiative transfer equations suggesting some restrictions on the physical quantities of molecular gas on a galactic scale.

In 2011 February, a burst of the 22 GHz H2O maser in Orion KL was reported. In order to identify the bursting maser features, we have been carrying out observations of the 22 GHz H2O maser in Orion KL with VERA, a Japanese VLBI network dedicated for astrometry. The bursting maser turns out to consist of two spatially different features at 7.58 and 6.95 km s−1. We determine their absolute positions and find that they are coincident with the shocked molecular gas called the Orion Compact Ridge. We tentatively detect the absolute proper motions of the bursting features toward the southwest direction, perpendicular to the elongation of the maser features. It is most likely that the outflow from the radio source I or another young stellar object interacting with Compact Ridge is a possible origin of the H2O maser burst. We will also carry out observations with ALMA in the cycle 0 period to monitor the submillimeter H2O maser lines in the Orion Compact Ridge region. These follow-up observations will provide novel information on the physical and chemical properties of the mastering region.