The Maser Monitoring Parkes Project (M2P2) is an ongoing project to observe masers towards high-mass star-forming regions (HMSFRs) using the 64 m CSIRO Parkes radio telescope, Murriyang. In this paper, we outline the project and introduce Stokes-I data from the first two years of observations. For the 63 sightlines observed in this project we identify a total of 1 514 individual maser features: 14.4% of these (203) towards 27 sightlines show significant variability. Most of these (160/203) are seen in the main-line transitions of OH at 1665 and 1667 MHz, but this data set also includes a significant number of variable features in the satellite lines at 1 612 and 1 720 MHz (33 and 10, respectively), most of which (24 and 9, respectively) appear to be associated with the HMSFRs. We divide these features into 4 broad categories based on the behaviour of their intensity over time: flares (6%), periodic (11%), long-term trends (33%), and ‘other’ (50%). Variable masers provide a unique laboratory for the modelling of local environmental conditions of HMSFRs, and follow-up publications will delve into this in more detail.

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.

We present the first unbiased survey of neutral hydrogen absorption in the Small Magellanic Cloud. The survey utilises pilot neutral hydrogen observations with the Australian Square Kilometre Array Pathfinder telescope as part of the Galactic Australian Square Kilometre Array Pathfinder neutral hydrogen project whose dataset has been processed with the Galactic Australian Square Kilometre Array Pathfinder-HI absorption pipeline, also described here. This dataset provides absorption spectra towards 229 continuum sources, a 275% increase in the number of continuum sources previously published in the Small Magellanic Cloud region, as well as an improvement in the quality of absorption spectra over previous surveys of the Small Magellanic Cloud. Our unbiased view, combined with the closely matched beam size between emission and absorption, reveals a lower cold gas faction (11%) than the 2019 ATCA survey of the Small Magellanic Cloud and is more representative of the Small Magellanic Cloud as a whole. We also find that the optical depth varies greatly between the Small Magellanic Cloud’s bar and wing regions. In the bar we find that the optical depth is generally low (correction factor to the optically thin column density assumption of $\mathcal{R}_{\mathrm{HI}} \sim 1.04$ ) but increases linearly with column density. In the wing however, there is a wide scatter in optical depth despite a tighter range of column densities.

We present the initial results from a class I 44-GHz methanol maser follow-up survey, observed with the MOPRA telescope, towards 272 sources from the Methanol Multi-beam survey (MMB). Over half (∼60%) of the 6.7 GHz class II MMB maser sources are associated with a class I 44-GHz methanol maser at a greater than 5σ detection level. We find that class II MMB masers sources with an associated class I methanol maser have stronger peak fluxes compared to regions without an associated class I maser. Furthermore, as part of the MOPRA follow-up observations we simultaneously observed SiO emission which is a known tracer of shocks and outflows in massive star forming regions. The presence of SiO emission, and potentially outflows, is found to be strongly associated with the detection of class I maser emission in these regions.

Masers are becoming increasingly important probes of high-mass star formation, revealing details about the kinematics and physical conditions at the elusive, early stages of formation. Over the last decade significant investment has been made in a number of large-scale, sensitive maser surveys targeting transitions found in the vicinity of young, high-mass stars. Individually, these searches have led to valuable insights into maser populations, their associated star formation regions, and often revealed further details such as Galactic structure. In combination, they become even more powerful, especially when considered together with complementary multi-wavelength data. Another consequence of large maser surveys has been the identification of a number of especially interesting sources that have been the subject of subsequent detailed studies. I summarize the recent plethora of maser surveys, their results, and how they are contributing to our understanding of star formation. Ongoing searches will ensure a bright future of maser surveys in the decade to come.

We have detected maser emission from the 36.2 GHz (4−1 → 30E) methanol transition towards NGC 4945. This emission has been observed in two separate epochs and is approximately five orders of magnitude more luminous than typical emission from this transition within our Galaxy. NGC 4945 is only the fourth extragalactic source observed hosting class I methanol maser emission. Extragalactic class I methanol masers do not appear to be simply highly-luminous variants of their galactic counterparts and instead appear to trace large-scale regions where low-velocity shocks are present in molecular gas.

The Australia Telescope Compact Array (ATCA) participated in a number of survey programs to search for and image common class I methanol masers (at 36 and 44 GHz) with high angular resolution. In this paper, we discuss spatial and velocity distributions revealed by these surveys. In particular, the number of maser regions is found to fall off exponentially with the linear distance from the associated young stellar object traced by the 6.7-GHz maser, and the scale of this distribution is 263±15 milliparsec. Although this relationship still needs to be understood in the context of the broader field, it can be utilised to estimate the distance using methanol masers only. This new technique has been analysed to understand its limitations and future potential. It turned out, it can be very successful to resolve the ambiguity in kinematic distances, but, in the current form, is much less accurate (than the kinematic method) if used on its own.

We applied three statistical classification techniques—linear discriminant analysis (LDA), logistic regression, and random forests—to three astronomical datasets associated with searches for interstellar masers. We compared the performance of these methods in identifying whether specific mid-infrared or millimetre continuum sources are likely to have associated interstellar masers. We also discuss the interpretability of the results of each classification technique. Non-parametric methods have the potential to make accurate predictions when there are complex relationships between critical parameters. We found that for the small datasets the parametric methods logistic regression and LDA performed best, for the largest dataset the non-parametric method of random forests performed with comparable accuracy to parametric techniques, rather than any significant improvement. This suggests that at least for the specific examples investigated here accuracy of the predictions obtained is not being limited by the use of parametric models. We also found that for LDA, transformation of the data to match a normal distribution led to a significant improvement in accuracy. The different classification techniques had significant overlap in their predictions; further astronomical observations will enable the accuracy of these predictions to be tested.

Determining an evolutionary clock for high-mass star formation is an important step towards realising a unified theory of star formation, as it will enable qualitative studies of the associated high-mass stars to be executed. We have carried out detailed studies of a large number of sources suspected of undergoing high-mass star formation and have found that common maser transitions offer the best opportunity to determine an evolutionary scheme for these objects. We have investigated the relative evolutionary phases of massive star formation associated with the presence or absence of combinations of water, methanol and main-line hydroxyl masers. The locations of the different maser species have been compared with the positions of 1.2 mm dust clumps, radio continuum, GLIMPSE point sources and Extended Green Objects. Comparison between the characteristics of coincident sources has revealed strong evidence for an evolutionary sequence for the different maser species in high-mass star formation regions. We present our proposed sequence for the presence of the common maser species associated with young high-mass stars and highlight recent advances. We discuss future investigations that will be made in this area by comparing data from the Methanol Multibeam (MMB) Survey with chemical clocks from the Millimetre Astronomy Legacy Team 90 GHz (MALT90) Survey.

The Australia Telescope Compact Array has been used to observe all the 603 6.7 GHz methanol masers detected in the Methanol Multi-Beam survey between l = 310° - 20°. To date we have measured positions with arcsecond accuracy for all the observations in the l = 6° - 20°.