MCSED: A flexible spectral energy distribution fitting code and its application to $z\sim 2$ emission-line galaxies

Abstract

We present MCSED, a new spectral energy distribution (SED)-fitting code, which mates flexible stellar evolution calculations with the Markov Chain Monte Carlo algorithms of the software package emcee. MCSED takes broad, intermediate, and narrow-band photometry, emission-line fluxes, and/or absorption line spectral indices, and returns probability distributions and co-variance plots for all model parameters. MCSED includes a variety of dust attenuation curves with parameters for varying the UV slopes and bump strengths, a prescription for continuum and PAH emission from dust, models for continuum and line emission from ionized gas, options for fixed and variable stellar metallicity, and a selection of star formation rate (SFR) histories. The code is well-suited for exploring parameter inter-dependencies in sets of galaxies with known redshifts, for which there is multi-band photometry and/or spectroscopy. We apply MCSED to a sample of $\sim 2000$ $1.90<z<2.35$ galaxies in the five CANDELS fields, which were selected via their strong [O III] $\lambda 5007$ emission, and explore the systematic behavior of their SEDs. We find the galaxies become redder with stellar mass, due to both increasing internal attenuation and a greater population of older stars. The slope of the UV extinction curve also changes with stellar mass, and at least some galaxies exhibit an extinction excess at 2175 Angstroms. Finally, we demonstrate that below $M\lesssim 10^9{M}_{\odot }$), the shape of the star-forming galaxy main sequence is highly dependent on the galaxies' assumed SFR history, as calculations which assume a constant SFR produce stellar masses that are $\sim 1$ dex smaller than those found using more realistic SFR histories.

27 pages, 18 figures, accepted to ApJ

References in corpus (86)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Astropy: A Community Python Package for Astronomy
• emcee: The MCMC Hammer
• The NumPy array: a structure for efficient numerical computation
• Stellar population synthesis at the resolution of 2003
• Galactic Stellar and Substellar Initial Mass Function
• On the variation of the Initial Mass Function
• Star Formation in Galaxies Along the Hubble Sequence
• The Dust Content and Opacity of Actively Star-Forming Galaxies
• Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
• Star Formation in the Milky Way and Nearby Galaxies
• Theoretical Modeling of Starburst Galaxies
• Evolutionary tracks and isochrones for low- and intermediate-mass stars: from 0.15 to 7 M_sun, and from Z=0.0004 to 0.03
• The UKIRT Infrared Deep Sky Survey (UKIDSS)
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results
• The Average Star Formation Histories of Galaxies in Dark Matter Halos from z=0-8
• Infrared Emission from Interstellar Dust. IV. The Silicate-Graphite-PAH Model in the Post-Spitzer Era
• The Cosmic Evolution Survey (COSMOS) -- Overview
• The propagation of uncertainties in stellar population synthesis modeling I: The relevance of uncertain aspects of stellar evolution and the IMF to the derived physical properties of galaxies
• Stellar mass-to-light ratios and the Tully-Fisher relation
• A Simple Model for the Absorption of Starlight by Dust in Galaxies
• The Great Observatories Origins Deep Survey: Initial Results From Optical and Near-Infrared Imaging
• A Highly Consistent Framework for the Evolution of the Star-Forming "Main Sequence" from z~0-6
• Evolution of asymptotic giant branch stars II. Optical to far-infrared isochrones with improved TP-AGB models
• A simple model to interpret the ultraviolet, optical and infrared emission from galaxies
• UV Luminosity Functions at redshifts z~4 to z~10: 10000 Galaxies from HST Legacy Fields
• The propagation of uncertainties in stellar population synthesis modeling III: model calibration, comparison, and evaluation
• Using Strong Lines to Estimate Abundances in Extragalactic HII Regions and Starburst Galaxies
• Dust Masses, PAH Abundances, and Starlight Intensities in the SINGS Galaxy Sample
• An ultra-deep near-infrared spectrum of a compact quiescent galaxy at z=2.2
• 3D-HST WFC3-selected Photometric Catalogs in the Five CANDELS/3D-HST Fields: Photometry, Photometric Redshifts and Stellar Masses
• Analysis of galaxy SEDs from far-UV to far-IR with CIGALE: Studying a SINGS test sample
• Modeling the Panchromatic Spectral Energy Distributions of Galaxies
• Constraining the Low-Mass Slope of the Star Formation Sequence at 0.5<z<2.5
• Inferring the star-formation histories of massive quiescent galaxies with BAGPIPES: Evidence for multiple quenching mechanisms
• The 3D-HST Survey: Hubble Space Telescope WFC3/G141 grism spectra, redshifts, and emission line measurements for $\sim 100,000$ galaxies
• The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years
• Strong Nebular Line Ratios in the Spectra of z~2-3 Star-forming Galaxies: First Results from KBSS-MOSFIRE
• The All-wavelength Extended Groth Strip International Survey (AEGIS) Data Sets
• PRIMUS: Constraints on Star Formation Quenching and Galaxy Merging, and the Evolution of the Stellar Mass Function From z=0-1
• Star Formation and Dust Attenuation Properties in Galaxies from a Statistical UV-to-FIR Analysis
• Understanding Galaxy Evolution through Emission Lines
• Deriving Physical Properties from Broadband Photometry with Prospector: Description of the Model and a Demonstration of its Accuracy Using 129 Galaxies in the Local Universe
• The effect of massive binaries on stellar populations and supernova progenitors
• The Origin and Evolution of the Galaxy Mass-Metallicity Relation
• UV Continuum Slope and Dust Obscuration from z~6 to z~2: The Star Formation Rate Density at High Redshift
• [OII] as a Star Formation Rate Indicator
• Modelling and interpreting spectral energy distributions of galaxies with BEAGLE
• Nebular Continuum and Line Emission in Stellar Population Synthesis Models
• The Dust Attenuation Law in Distant Galaxies: Evidence for Variation with Spectral Type
• Stellar population models at high spectral resolution
• Fitting the integrated Spectral Energy Distributions of Galaxies
• The Relation Between SFR and Stellar Mass for Galaxies at 3.5 $\le z\le$ 6.5 in CANDELS
• Optical Star-Formation Rate Indicators
• The global dust modelling framework THEMIS (The Heterogeneous dust Evolution Model for Interstellar Solids)
• The Characteristic Star Formation Histories of Galaxies at Redshifts z~2-7
• New Strong Line Abundance Diagnostics for \HII Regions: Effects of $\kappa$-Distributed Electron Energies and New Atomic Data
• The infrared side of galaxy formation. I. The local universe in the semi-analytical framework
• The MOSDEF Survey: Excitation Properties of $z\sim 2.3$ Star-forming Galaxies
• Empirical determination of the shape of dust attenuation curves in star-forming galaxies
• Great Optically Luminous Dropout Research Using Subaru HSC (GOLDRUSH). I. UV Luminosity Functions at $z\sim 4-7$ Derived with the Half-Million Dropouts on the 100 deg${}^2$ Sky
• GOODS-Herschel:dust attenuation properties of UV selected high redshift galaxies
• Firefly (Fitting IteRativEly For Likelihood analYsis): a full spectral fitting code
• The CANDELS/SHARDS multi-wavelength catalog in GOODS-N: Photometry, Photometric Redshifts, Stellar Masses, Emission line fluxes and Star Formation Rates
• How well can we really estimate the stellar masses of galaxies from broadband photometry?
• Emission-line diagnostics of nearby HII regions including interacting binary populations
• Reconstruction of Galaxy Star Formation Histories through SED Fitting: The Dense Basis Approach
• Should we believe the results of UV-mm galaxy SED modelling?
• The Nature of Extreme Emission Line Galaxies at z=1-2: Kinematics and Metallicities from Near-Infrared Spectroscopy
• Uv-to-fir analysis of spitzer/irac sources in the extended groth strip i: Multi-wavelength photometry and spectral energy distributions
• Dust attenuation in disk-dominated galaxies: evidence for the 2175A dust feature
• Characterizing Dust Attenuation in Local Star Forming Galaxies: UV and Optical Reddening
• The dependence of the galaxy mass-metallicity relation on environment and the implied metallicity of the IGM
• SED fitting with Markov Chain Monte Carlo: Methodology and Application to z=3.1 Lyman Alpha Emitting Galaxies
• The rise and fall of the star formation histories of blue galaxies at redshifts 0.2<z<1.4
• CIGALEMC: Galaxy Parameter Estimation using a Markov Chain Monte Carlo Approach with Cigale
• Simultaneous Modeling of the Stellar and Dust Emission in Distant Galaxies: Implications for Star Formation Rate Measurements
• Measurement of [OIII] Emission in Lyman Break Galaxies
• BayeSED: A General Approach to Fitting the Spectral Energy Distribution of Galaxies
• Interpreting high [O III]/Hbeta ratios with maturing starbursts
• Characterizing Dust Attenuation in Local Star Forming Galaxies: Inclination Effects and the 2175Å Feature
• The Dust Attenuation Curve versus Stellar Mass for Emission Line Galaxies at z ~ 2
• Swift UV/Optical Telescope Imaging of Star Forming Regions in M81 and Holmberg IX
• The Relation Between [OIII]/H$\beta$ and Specific Star Formation Rate in Galaxies at $z\sim 2$
• Galaxies of the $z\sim 2$ Universe. I. Grism-Selected Rest-Frame Optical Emission Line Galaxies
• A Cautionary Tale of Attenuation in Star Forming Regions

Cited by in corpus (13)

• The Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques
• Value Added Catalog of physical properties of more than 1.3 million galaxies from the DESI Survey
• Stochastic Modelling of Star Formation Histories III. Constraints from Physically-Motivated Gaussian Processes
• The HETDEX Survey: The Ly$\alpha$ Escape Fraction from 3D-HST Emission Line Galaxies at $z\sim 2$
• Lyα Halos around [O III]-Selected Galaxies in HETDEX
• HETDEX [OIII] Emitters I: A spectroscopically selected low-redshift population of low-mass, low-metallicity galaxies
• New-generation dust emission templates for star-forming galaxies
• Measuring Stellar Masses of Emission Line Galaxies at 1.2<z<1.9
• The H$\alpha$ and [O III] $\lambda 5007$ Luminosity Functions of $1.2<z<1.9$ Emission-Line Galaxies from HST Grism Spectroscopy
• Measuring Dust Attenuation Curves of SINGS/KINGFISH Galaxies Using Swift/UVOT Photometry
• Stellar, Gas, and Dust Emission of Star Forming Galaxies out to z~2
• The $z\sim 2$ $[\mathrm{O}\mathrm{III}]$ Luminosity Function of Grism-selected Emission-line Galaxies
• Revisiting the Local Star-Forming Galaxies Observed in the HETDEX Pilot Survey