Learning the Fundamental MIR Spectral Components of Galaxies with Non-Negative Matrix Factorisation

Abstract

The mid-infrared (MIR) spectra observed with the \textit{Spitzer} Infrared Spectrograph (IRS) provide a valuable dataset for untangling the physical processes and conditions within galaxies. This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorisation (NMF). NMF is a recently developed multivariate technique shown to be successful in blind source separation problems. Unlike the more popular multivariate analysis technique, principal component analysis, NMF imposes the condition that weights and spectral components are non-negative. This more closely resembles the physical process of emission in the mid-infrared, resulting in physically intuitive components. By applying NMF to galaxy spectra in the Cornell Atlas of Spitzer/IRS sources (CASSIS), we find similar components amongst different NMF sets. These similar components include two for AGN emission and one for star formation. [... ABBREVIATED...] We show an NMF set with seven components can reconstruct the general spectral shape of a wide variety of objects, though struggle to fit the varying strength of emission lines. We also show that the seven components can be used to separate out different types of objects. We model this separation with Gaussian Mixtures modelling and use the result to provide a classification tool. We also show the NMF components can be used to separate out the emission from AGN and star formation regions and define a new star formation/AGN diagnostic which is consistent with all mid-infrared diagnostics already in use but has the advantage that it can be applied to mid-infrared spectra with low signal to noise or with limited spectral range. The 7 NMF components and code for classification are made public on arxiv and are available at: \url{https://github.com/pdh21/NMF_software/}

23 pages, 19 figures, accepted for publication in MNRAS

References in corpus (14)

• K-corrections and filter transformations in the ultraviolet, optical, and near infrared
• Mid-IR Galaxy Classification Based on Silicate Obscuration and PAH Equivalent Width
• Mid-Infrared Spectral Diagnosis of Submillimeter Galaxies
• High resolution mid-infrared spectroscopy of ultraluminous infrared galaxies
• The nature of star formation in distant ultraluminous infrared galaxies selected in a remarkably narrow redshift range
• Spitzer Mid-Infrared Spectroscopy of Infrared Luminous Galaxies at z~2 II: Diagnostics
• Mid-Infrared Spectral Diagnostics of Nuclear and Extra-Nuclear Regions in Nearby Galaxies
• GOODS-Herschel: Ultra-deep XMM-Newton observations reveal AGN/star-formation connection
• Extended silicate dust emission in PG QSOs
• Mid-Infrared Diagnostics of LINERs
• The Origin of the Silicate Emission Features in the Seyfert 2 Galaxy, NGC 2110
• An Evolutionary Paradigm for Dusty Active Galaxies at Low Redshift
• Parametric modelling of the 3.6um to 8um colour distributions of galaxies in the SWIRE Survey
• Single parameter galaxy classification: The Principal Curve through the multi-dimensional space of galaxy properties

Cited by in corpus (8)

• How Many Components? Quantifying the Complexity of the Metallicity Distribution in the Milky Way Bulge with APOGEE
• SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry
• Multivariate Approaches to Classification in Extragalactic Astronomy
• A statistical study of the optical spectral variability in gamma-ray blazars
• Classification of local ultraluminous infrared galaxies and quasars with kernel principal component analysis
• Sengi: a small, fast, interactive viewer for spectral outputs from stellar population synthesis models
• Accurate spectroscopic redshift estimation using non-negative matrix factorization: application to MUSE spectra
• SpectralUnmix: A Torch-Based Regularized Non-negative Matrix Factorization