SMART: spectral energy distributions Markov chain analysis with radiative transfer models

Charalambia Varnava,Andreas Efstathiou
DOI: https://doi.org/10.1093/mnras/stae1141
2024-05-28
Abstract:In this paper we present the publicly available open-source spectral energy distribution (SED) fitting code SMART (Spectral energy distributions Markov chain Analysis with Radiative Transfer models). Implementing a Bayesian Markov chain Monte Carlo (MCMC) method, SMART fits the ultraviolet to millimetre SEDs of galaxies exclusively with radiative transfer models that currently constitute four types of pre-computed libraries, which describe the starburst, active galactic nucleus (AGN) torus, host galaxy and polar dust components. An important novelty of SMART is that, although it fits SEDs exclusively with radiative transfer models, it takes comparable time to popular energy balance methods to run. Here we describe the key features of SMART and test it by fitting the multi-wavelength SEDs of the 42 local ultraluminous infrared galaxies (ULIRGs) that constitute the HERschel Ultraluminous Infrared Galaxy Survey (HERUS) sample. The Spitzer spectroscopy data of the HERUS ULIRGs are included in the fitting at a spectral resolution, which is matched to that of the radiative transfer models. We also present other results that highlight the performance and versatility of SMART. SMART promises to be a useful tool for studying galaxy evolution in the JWST era. SMART is developed in PYTHON and is available at <a class="link-external link-https" href="https://github.com/ch-var/SMART.git" rel="external noopener nofollow">this https URL</a>.
Astrophysics of Galaxies,Cosmology and Nongalactic Astrophysics
What problem does this paper attempt to address?
This paper presents an open-source code called SMART (Spectral Energy Distributions Markov chain Analysis with Radiative Transfer models) for analyzing the spectral energy distribution of celestial objects. SMART utilizes Bayesian Markov chain Monte Carlo (MCMC) methods specifically designed for modeling galaxy spectra spanning from ultraviolet to millimeter wavelengths. It includes four pre-computed model libraries for starbursts, active galactic nuclei (AGN) dust tori, host galaxies, and polar dust. The paper highlights that a key innovation of SMART is that despite being based entirely on radiative transfer models, its computational time is comparable to popular energy balance methods. SMART is applicable for studying multi-wavelength observations of local galaxies with Spitzer and JWST data as well as distant galaxies. It considers the crucial role of infrared spectral photometry in constraining the energy sources, such as starbursts and AGN activities. The paper demonstrates the performance and versatility of SMART by testing multi-wavelength spectra of 42 local ultra-luminous infrared galaxies (ULIRGs). SMART can be used for investigating galaxy evolution, particularly in the era of JWST. Furthermore, the paper discusses the comparison of SMART with other existing tools like CIGALE and mentions the future development directions of SMART. The uniqueness of SMART lies in: 1. It uses the full range of model parameters, especially for key parameters like AGN dust torus inclination. 2. It is parallelized for faster computation compared to other MCMC codes. 3. It provides four different types of AGN model libraries. 4. Optionally, it can include polar dust components. 5. It can fit a specific part of a galaxy with a time-dependent model similar to spheroidal or disk-shaped host galaxy models. The structure of the paper includes methods introduction, testing SMART, case studies, and conclusions. SMART is developed in PYTHON and is publicly available on GitHub.