Monte Carlo Method for Calculating Oxygen Abundances and Their Uncertainties from Strong-Line Flux Measurements

Abstract

We present the open-source Python code pyMCZ that determines oxygen abundance and its distribution from strong emission lines in the standard metallicity calibrators, based on the original IDL code of Kewley & Dopita (2002) with updates from Kewley & Ellison (2008), and expanded to include more recently developed calibrators. The standard strong-line diagnostics have been used to estimate the oxygen abundance in the interstellar medium through various emission line ratios in many areas of astrophysics, including galaxy evolution and supernova host galaxy studies. We introduce a Python implementation of these methods that, through Monte Carlo sampling, better characterizes the statistical oxygen abundance confidence region including the effect due to the propagation of observational uncertainties. These uncertainties are likely to dominate the error budget in the case of distant galaxies, hosts of cosmic explosions. Given line flux measurements and their uncertainties, our code produces synthetic distributions for the oxygen abundance in up to 15 metallicity calibrators simultaneously, as well as for E(B-V), and estimates their median values and their 68% confidence regions. We test our code on emission line measurements from a sample of nearby supernova host galaxies (z < 0.15) and compare our metallicity results with those from previous methods. Our metallicity estimates are consistent with previous methods but yield smaller statistical uncertainties. Systematic uncertainties are not taken into account. We offer visualization tools to assess the spread of the oxygen abundance in the different calibrators, as well as the shape of the estimated oxygen abundance distribution in each calibrator, and develop robust metrics for determining the appropriate Monte Carlo sample size. The code is open access and open source and can be found at https://github.com/nyusngroup/pyMCZ (Abridged)

Accepted for publication in Astronomy& Computing, 2016. The code is open-access, open-source, and suggestions and improvements are welcome! The python module is available at https://github.com/nyusngroup/pyMCZ

References in corpus (17)

• The Host Galaxies and Classification of Active Galactic Nuclei
• Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory
• The MAPPINGS III Library of Fast Radiative Shock Models
• Measured Metallicities at the Sites of Nearby Broad-Lined Type Ic Supernovae and Implications for the SN-GRB Connection
• Hydrogen-Poor Superluminous Supernovae and Long-Duration Gamma-Ray Bursts Have Similar Host Galaxies
• Spectroscopy of superluminous supernova host galaxies. A preference of hydrogen-poor events for extreme emission line galaxies
• Deriving model-based T${}_e$-consistent chemical abundances in ionised gaseous nebulae
• Resolving the electron temperature discrepancies in HII Regions and Planetary Nebulae: kappa-distributed electrons
• Progenitor Diagnostics for Stripped Core-Collapse Supernovae: Measured Metallicities at Explosion Sites
• IZI: Inferring the Gas Phase Metallicity (Z) and Ionization Parameter (q) of Ionized Nebulae using Bayesian Statistics
• The Host Galaxies of Type Ia Supernovae Discovered by the Palomar Transient Factory
• CHAOS II: Gas-Phase Abundances in NGC 5194
• Imprints of galaxy evolution on H ii regions Memory of the past uncovered by the CALIFA survey
• The properties of SN Ib/c locations
• Error estimation in astronomy: A guide
• Testing the existence of non-Maxwellian electron distributions in H II regions after assessing atomic data accuracy
• Data analysis recipes: Choosing the binning for a histogram

Cited by in corpus (30)

• Host-Galaxy Properties of 32 Low-Redshift Superluminous Supernovae from the Palomar Transient Factory
• Energetic eruptions leading to a peculiar hydrogen-rich explosion of a massive star
• The Koala: A Fast Blue Optical Transient with Luminous Radio Emission from a Starburst Dwarf Galaxy at $z=0.27$
• A hot and fast ultra-stripped supernova that likely formed a compact neutron star binary
• Complexity in the light curves and spectra of slow-evolving superluminous supernovae
• Discovery of Strongly Inverted Metallicity Gradients in Dwarf Galaxies at $z$$\sim$2
• The Grism Lens-Amplified Survey from Space (GLASS) X. Sub-kpc resolution gas-phase metallicity maps at cosmic noon behind the Hubble Frontier Fields cluster MACS1149.6+2223
• iPTF16asu: A Luminous, Rapidly-Evolving, and High-Velocity Supernova
• Luminous Type II Short-Plateau Supernovae 2006Y, 2006ai, and 2016egz: A Transitional Class from Stripped Massive Red Supergiants
• The evolution of superluminous supernova LSQ14mo and its interacting host galaxy system
• Interrogating Seyferts with NebulaBayes: Spatially probing the narrow-line region radiation fields and chemical abundances
• Host Galaxies of Type Ic and Broad-lined Type Ic Supernovae from the Palomar Transient Factory: Implication for Jet Production
• A Recalibration of Strong Line Oxygen Abundance Diagnostics via the Direct Method and Implications for the High Redshift Universe
• Metallicity, ionization parameter, and pressure variations of HII regions in the TYPHOON spiral galaxies
• Spatially resolved MaNGA observations of the host galaxy of superluminous supernova 2017egm
• Ionised gas structure of 100 kpc in an over-dense region of the galaxy group COSMOS-Gr30 at z ~ 0.7
• iPTF 16hgs: A double-peaked Ca-rich gap transient in a metal poor, star forming dwarf galaxy
• A low-energy explosion yields the underluminous Type IIP SN 2020cxd
• A Revised Host Galaxy Association for GRB 020819B: A High-Redshift Dusty Starburst, Not a Low-Redshift Gas-Poor Spiral
• The luminous host galaxy, faint supernova and rapid afterglow rebrightening of GRB 100418A
• Swift UVOT Grism Observations of Nearby Type Ia Supernovae -- II. Probing the Progenitor Metallicity of SNe Ia with Ultraviolet Spectra
• The Bimodal Absorption System Imaging Campaign (BASIC) I. A Dual Population of Low-metallicity Absorbers at z $<1$
• KSP-SN-2016kf: a long-rising H-rich Type II Supernova with unusually high ${}^{56}$Ni mass discovered in the KMTNet Supernova Program
• The Exotic Type Ic Broad-Lined Supernova SN 2018gep: Blurring the Line Between Supernovae and Fast Optical Transients
• Extraplanar H II Regions in Spiral Galaxies. I. Low-Metallicity Gas Accreting through the Disk-Halo Interface of NGC 4013
• The host of the Type I SLSN 2017egm: A young, sub-solar metallicity environment in a massive spiral galaxy
• Close, bright and boxy: the superluminous SN 2018hti
• SN 2020zbf: A fast-rising hydrogen-poor superluminous supernova with strong carbon lines
• Surface density: a new parameter in the fundamental metallicity relation of star-forming galaxies
• SN 2019hcc: A Type II Supernova Displaying Early O II Lines