Principal Component Analysis and Radiative Transfer modelling of Spitzer IRS Spectra of Ultra Luminous Infrared Galaxies

Abstract

The mid-infrared spectra of ultraluminous infrared galaxies (ULIRGs) contain a variety of spectral features that can be used as diagnostics to characterise the spectra. However, such diagnostics are biased by our prior prejudices on the origin of the features. Moreover, by using only part of the spectrum they do not utilise the full information content of the spectra. Blind statistical techniques such as principal component analysis (PCA) consider the whole spectrum, find correlated features and separate them out into distinct components. We further investigate the principal components (PCs) of ULIRGs derived in Wang et al.(2011). We quantitatively show that five PCs is optimal for describing the IRS spectra. These five components (PC1-PC5) and the mean spectrum provide a template basis set that reproduces spectra of all z<0.35 ULIRGs within the noise. For comparison, the spectra are also modelled with a combination of radiative transfer models of both starbursts and the dusty torus surrounding active galactic nuclei. The five PCs typically provide better fits than the models. We argue that the radiative transfer models require a colder dust component and have difficulty in modelling strong PAH features. Aided by the models we also interpret the physical processes that the principal components represent. The third principal component is shown to indicate the nature of the dominant power source, while PC1 is related to the inclination of the AGN torus. Finally, we use the 5 PCs to define a new classification scheme using 5D Gaussian mixtures modelling and trained on widely used optical classifications. The five PCs, average spectra for the four classifications and the code to classify objects are made available at: http://www.phys.susx.ac.uk/~pdh21/PCA/

11 pages, 12 figures, accepted for publication in MNRAS

References in corpus (10)

• Mid-IR Galaxy Classification Based on Silicate Obscuration and PAH Equivalent Width
• CASSIS: The Cornell Atlas of Spitzer/IRS Sources
• High resolution mid-infrared spectroscopy of ultraluminous infrared galaxies
• Dust in starburst nuclei and ULIRGs: SED models for observers
• A Spitzer IRS Low-Resolution Spectroscopic Search for Buried AGNs in Nearby Ultraluminous Infrared Galaxies - A Constraint on Geometry between Energy Sources and Dust -
• The nature of star formation in distant ultraluminous infrared galaxies selected in a remarkably narrow redshift range
• Direct Evidence for Termination of Obscured Star Formation by Radiatively Driven Outflows in Reddened QSOs
• A model for the infrared emission of FSC 10214+4724
• 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

Cited by in corpus (7)

• CASSIS: The Cornell Atlas of Spitzer/Infrared Spectrograph Sources. II. High-resolution observations
• X-ray Emission from Strongly Asymmetric Circumstellar Material in the Remnant of Kepler's Supernova
• Analysis of Spitzer-IRS spectra of hyperluminous infrared galaxies
• X-ray Study of Spatial Structures in Tycho's Supernova Remnant Using Unsupervised Deep Learning
• A Principal Component Analysis of polycyclic aromatic hydrocarbon emission in NGC 7023
• Classification of local ultraluminous infrared galaxies and quasars with kernel principal component analysis
• Graph Theoretical Analysis of local ultraluminous infrared galaxies and quasars