The BPT Diagram in Cosmological Galaxy Formation Simulations: Understanding the Physics Driving Offsets at High-Redshift

Abstract

The Baldwin, Philips, & Terlevich diagram of [O III]/H$\beta$ vs. [N II]/H$\alpha$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim 2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the SIMBA cosmological hydrodynamic galaxy formation simulation, using the CLOUDY photoionization code to compute the nebular line luminosities from H II regions. We find that the observed shift toward higher [O III]/H$\beta$ and [N II]/H$\alpha$ values at high redshift arises from sample selection: when we consider only the most massive galaxies $M_{\ast }\sim 10^{10-11}{M}_{\odot }$, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to $z\sim 0$ observations. Even when accounting for sample selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H II region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.

25 pages, 15 figures. Accepted for publication in ApJ

References in corpus (36)

• Modules for Experiments in Stellar Astrophysics (MESA)
• 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
• Simba: Cosmological Simulations with Black Hole Growth and Feedback
• Binary Population and Spectral Synthesis Version 2.1: construction, observational verification and new results
• Strong Nebular Line Ratios in the Spectra of z~2-3 Star-forming Galaxies: First Results from KBSS-MOSFIRE
• Southern Massive Stars at High Angular Resolution: Observational Campaign and Companion Detection
• A systematic variation of the stellar initial mass function in early-type galaxies
• The High Angular Resolution Multiplicity of Massive Stars
• Understanding Galaxy Evolution through Emission Lines
• The Atomic to Molecular Transition in Galaxies. II: HI and H_2 Column Densities
• Grackle: a Chemistry and Cooling Library for Astrophysics
• Modelling the nebular emission from primeval to present-day star-forming galaxies
• New fully empirical calibrations of strong-line metallicity indicators in star forming galaxies
• Nebular Emission Line Ratios in z~2-3 Star-Forming Galaxies with KBSS-MOSFIRE: Exploring the Impact of Ionization, Excitation, and Nitrogen-to-Oxygen Ratio
• New insights into the stellar content and physical conditions of star-forming galaxies at z = 2-3 from spectral modelling
• SDSS-IV MaNGA: The Impact of Diffuse Ionized Gas on Emission-line Ratios, Interpretation of Diagnostic Diagrams, and Gas Metallicity Measurements
• The MOSDEF Survey: Excitation Properties of $z\sim 2.3$ Star-forming Galaxies
• Ultraviolet Emission Lines in Young Low Mass Galaxies at z~2: Physical Properties and Implications for Studies at z>7
• Metallicities and Physical Conditions in Star-forming Galaxies at z~1.0-1.5
• CHAOS I: Direct Chemical Abundances for HII Regions in NGC 628
• Rest-Frame Optical Spectra of Three Strongly Lensed Galaxies at z~2
• The oxygen abundance in the inner HII regions of M101. Implications for the calibration of strong-line metallicity indicators
• Gravitationally Lensed Galaxies at 2<z<3.5: Direct Abundance Measurements of Lya Emitters
• The Nature of Extreme Emission Line Galaxies at z=1-2: Kinematics and Metallicities from Near-Infrared Spectroscopy
• The Low Mass End of the Fundamental Relation for Gravitationally Lensed Star Forming Galaxies at 1<z<6
• Synthetic nebular emission from massive galaxies I: origin of the cosmic evolution of optical emission-line ratios
• Probing Cosmic Dawn with Emission Lines: Predicting Infrared and Nebular Line Emission for ALMA and JWST
• Biases in metallicity measurements from global galaxy spectra: the effects of flux-weighting and diffuse ionized gas contamination
• The binarity of Milky Way F,G,K stars as a function of effective temperature and metallicity
• The MOSDEF-LRIS Survey: The Connection Between Massive Stars and Ionized Gas in Individual Galaxies at $z\sim 2$
• First direct metallicity measurement of a lensed star-forming galaxy at z=1.7
• Star Cluster Mass and Age Distributions of Two Fields in M83 Based on HST/WFC3 Observations
• Star-forming Clumps in Local Luminous Infrared Galaxies
• The [α/Fe]-[Fe/H] relation in the E-MOSAICS simulations: its connection to the birth place of globular clusters and the fraction of globular cluster field stars in the bulge
• Massive Star Cluster Formation and Destruction in Luminous Infrared Galaxies in GOALS
• What Drives the Redshift Evolution of Strong Emission Line Ratios?

Cited by in corpus (12)

• Quenching and the UVJ diagram in the SIMBA cosmological simulation
• The resolved chemical abundance properties within the interstellar medium of star-forming galaxies at $\text{z}\approx 1.5$
• An X-ray Census of Fast Radio Burst Host Galaxies: Constraints on AGN and X-ray Counterparts
• First Light And Reionisation Epoch Simulations (FLARES) XI: [OIII] emitting galaxies at $5<z<10$
• First Light And Reionisation Epoch Simulations (FLARES) XII: The consequences of star-dust geometry on galaxies in the EoR
• Euclid preparation. Optical emission-line predictions of intermediate-z galaxy populations in GAEA for the Euclid Deep and Wide Surveys
• TODDLERS: A New UV to Millimeter Emission Library for Star-Forming Regions. II. Star Formation Rate Indicators Using Auriga Zoom Simulations
• MAMMOTH-MOSFIRE: Environmental Effects on Galaxy Interstellar Medium at $z\sim 2$
• Cloudy-Maraston: Integrating nebular continuum and line emission with the Maraston stellar population synthesis models
• CECILIA: Ultra-Deep Rest-Optical Spectra of Faint Galaxies at Cosmic Noon
• Recovering the properties of the interstellar medium through integrated spectroscopy: application to the z~0 ECO volume-limited star-forming galaxy sample
• Euclid: A blue galaxy population and a brightest cluster galaxy in the making in a $z\sim 1.74$ MaDCoWS2 galaxy cluster candidate