Massive young stellar objects in the Local Group spiral galaxy M33 identified using machine learning

TL;DR

This study applies machine learning to classify stellar populations in galaxy M33, identifying young stellar objects and star-forming regions, and estimating star formation rates with high accuracy using IR data.

Contribution

It introduces a supervised machine learning approach with Probabilistic Random Forest to classify diverse stellar objects and quantify star formation in M33, expanding methods used in extragalactic stellar population studies.

Findings

Classified 162,746 sources with 86% accuracy.

Identified 4,985 young stellar objects and 68 star-forming regions.

Estimated global star formation rate of 1.42 solar masses per year.

Abstract

We present a supervised machine learning classification of stellar populations in the Local Group spiral galaxy M\,33. The Probabilistic Random Forest (PRF) methodology, previously applied to populations in NGC\,6822, utilises both near and far-IR classification features. It classifies sources into nine target classes: young stellar objects (YSOs), oxygen- and carbon-rich asymptotic giant branch stars, red giant branch and red super-giant stars, active galactic nuclei, blue stars (e.g. O-, B- and A-type main sequence stars), Wolf-Rayet stars and Galactic foreground stars. Across 100 classification runs the PRF classified 162,746 sources with an average estimated accuracy of $\sim$\,86\,per\,cent, based on confusion matrices. We identified 4985 YSOs across the disk of M\,33, applying a density-based clustering analysis to identify 68 star forming regions (SFRs) primarily in the galaxy's spiral arms. SFR counterparts to known H\,{\sc ii} regions were recovered, with $\sim$\,91\,per\,cent of SFRs spatially coincident with giant molecular clouds identified in the literature. Using photometric measurements, as well as SFRs in NGC\,6822 with an established evolutionary sequence as a benchmark, we employed a novel approach combining ratios of [H$\alpha$]$/$[24$\mu$m] and [250$\mu$m]$/$[500$\mu$m] to estimate the relative evolutionary status of all M\,33 SFRs. Masses were estimated for each YSO ranging from 6\,$-$\,27\,M$_\odot$. Using these masses, we estimate star formation rates based on direct YSO counts of 0.63\,M$_\odot$\,yr$^{-1}$ in M\,33's SFRs, 0.79\,$\pm$\,0.16\,M$_\odot$\,yr$^{-1}$ in its centre and 1.42\,$\pm$\,0.16\,M$_\odot$\,yr$^{-1}$ globally.