skip to main content

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

The Baldwin, Philips, & Terlevich diagram of [Oiii]/Hβversus [Nii]/Hα(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 atz∼ 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 thesimbacosmological hydrodynamic galaxy formation simulation, using thecloudyphotoionization code to compute the nebular line luminosities from Hiiregions. We find that the observed shift toward higher [Oiii]/Hβand [Nii]/Hαvalues at high redshift arises from sample selection: when we consider only the most massive galaxiesM*∼ 1010–11M, 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 toz∼ 0 observations. Even when accounting for samples-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, Hiiregion densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the more » 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.

« less
Authors:
; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
2009313 1909153
Publication Date:
NSF-PAR ID:
10362890
Journal Name:
The Astrophysical Journal
Volume:
926
Issue:
1
Page Range or eLocation-ID:
Article No. 80
ISSN:
0004-637X
Publisher:
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Existing star-forming vs. active galactic nucleus (AGN) classification schemes using optical emission-line diagnostics mostly fail for low-metallicity and/or highly star-forming galaxies, missing AGN in typicalz∼ 0 dwarfs. To recover AGN in dwarfs with strong emission lines (SELs), we present a classification scheme optimizing the use of existing optical diagnostics. We use Sloan Digital Sky Survey emission-line catalogs overlapping the volume- and mass-limited REsolved Spectroscopy Of a Local VolumE (RESOLVE) and Environmental COntex (ECO) surveys to determine the AGN percentage in SEL dwarfs. Our photoionization grids show that the [Oiii]/Hβversus [Sii]/Hαdiagram (Siiplot) and [Oiii]/Hβversus [Oi]/Hαdiagram (Oiplot) are less metallicity sensitive and more successful in identifying dwarf AGN than the popular [Oiii]/Hβversus [Nii]/Hαdiagnostic (Niiplot or “BPT diagram”). We identify a new category of “star-forming AGN” (SF-AGN) classified as star-forming by the Niiplot but as AGN by the Siiand/or Oiplots. Including SF-AGN, we find thez∼ 0 AGN percentage in dwarfs with SELs to be ∼3%–16%, far exceeding most previous optical estimates (∼1%). The large range in our dwarf AGN percentage reflects differences in spectral fitting methodologies between catalogs. The highly complete nature of RESOLVE and ECO allows us to normalize strong emission-line galaxy statistics to the full galaxy population, reducing the dwarfmore »AGN percentage to ∼0.6%–3.0%. The newly identified SF-AGN are mostly gas-rich dwarfs with halo mass <1011.5M, where highly efficient cosmic gas accretion is expected. Almost all SF-AGN also have low metallicities (Z≲ 0.4Z), demonstrating the advantage of our method.

    « less
  2. Abstract

    We use Hubble Space Telescope Wide Field Camera 3 G102 and G141 grism spectroscopy to measure rest-frame optical emission-line ratios of 533 galaxies atz∼ 1.5 in the CANDELS LyαEmission at Reionization survey. We compare [Oiii]/Hβversus [Sii]/(Hα+ [Nii]) as an “unVO87” diagram for 461 galaxies and [Oiii]/Hβversus [Neiii]/[Oii] as an “OHNO” diagram for 91 galaxies. The unVO87 diagram does not effectively separate active galactic nuclei (AGN) and [Nev] sources from star-forming galaxies, indicating that the unVO87 properties of star-forming galaxies evolve with redshift and overlap with AGN emission-line signatures atz> 1. The OHNO diagram does effectively separate X-ray AGN and [Nev]-emitting galaxies from the rest of the population. We find that the [Oiii]/Hβline ratios are significantly anticorrelated with stellar mass and significantly correlated withlog(LHβ), while [Sii]/(Hα+ [Nii]) is significantly anticorrelated withlog(LHβ). Comparison with MAPPINGS V photoionization models indicates that these trends are consistent with lower metallicity and higher ionization in low-mass and high-star formation rate (SFR) galaxies. We do not find evidence for redshift evolution of the emission-line ratios outside of the correlations with mass and SFR. Our results suggest that the OHNO diagram of [Oiii]/Hβversus [Neiii]/[Oii] willmore »be a useful indicator of AGN content and gas conditions in very high-redshift galaxies to be observed by the James Webb Space Telescope.

    « less
  3. ABSTRACT We analyse the rest-optical emission-line spectra of z ∼ 2.3 star-forming galaxies in the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. In investigating the origin of the well-known offset between the sequences of high-redshift and local galaxies in the [O iii]λ5008/Hβ versus [N ii]λ6585/Hα (‘[N ii] BPT’) diagram, we define two populations of z ∼ 2.3 MOSDEF galaxies. These include the high population that is offset towards higher [O iii]λ5008/Hβ and/or [N ii]λ6585/Hα with respect to the local SDSS sequence and the low population that overlaps the SDSS sequence. These two groups are also segregated within the [O  iii]λ5008/Hβ versus [S ii]λλ6718,6733/Hα and the [O iii]λλ4960,5008/[O ii ]λλ3727,3730 (O32) versus ([O  iii]λλ4960,5008+[O ii]λλ3727,3730)/Hβ (R23) diagrams, which suggests qualitatively that star-forming regions in the more offset galaxies are characterized by harder ionizing spectra at fixed nebular oxygen abundance. We also investigate many galaxy properties of the split sample and find that the high sample is on average smaller in size and less massive, but has higher specific star formation rate (SFR) and SFR surface density values and is slightly younger compared to the low population. From Cloudy+BPASS photoionization models, we estimate that the high population has a lower stellar metallicity (i.e. harder ionizing spectrum) but slightly higher nebular metallicity and higher ionizationmore »parameter compared to the low population. While the high population is more α-enhanced (i.e. higher α/Fe) than the low population, both samples are significantly more α-enhanced compared to local star-forming galaxies with similar rest-optical line ratios. These differences must be accounted for in all high-redshift star-forming galaxies – not only those ‘offset’ from local excitation sequences.« less
  4. Abstract

    We present results on the nature of extreme ejective feedback episodes and the physical conditions of a population of massive (M*∼ 1011M), compact starburst galaxies atz= 0.4–0.7. We use data from Keck/NIRSPEC, SDSS, Gemini/GMOS, MMT, and Magellan/MagE to measure rest-frame optical and near-IR spectra of 14 starburst galaxies with extremely high star formation rate surface densities (mean ΣSFR∼ 2000Myr−1kpc−2) and powerful galactic outflows (maximum speedsv98∼ 1000–3000 km s−1). Our unique data set includes an ensemble of both emission ([Oii]λλ3726,3729, Hβ, [Oiii]λλ4959,5007, Hα, [Nii]λλ6549,6585, and [Sii]λλ6716,6731) and absorption (Mgiiλλ2796,2803, and Feiiλ2586) lines that allow us to investigate the kinematics of the cool gas phase (T∼ 104K) in the outflows. Employing a suite of line ratio diagnostic diagrams, we find that the central starbursts are characterized by high electron densities (medianne∼ 530 cm−3), and high metallicity (solar or supersolar). We show that the outflows are most likely driven by stellar feedback emerging from the extreme central starburst, rather than by an AGN. We also present multiple intriguing observational signatures suggesting that these galaxies may have substantial Lyman continuum (LyC) photon leakage, including weak [Sii]nebular emission lines. Our results imply that these galaxies may be captured in a short-lived phase of extrememore »star formation and feedback where much of their gas is violently blown out by powerful outflows that open up channels for LyC photons to escape.

    « less
  5. Abstract We present spatially resolved Hubble Space Telescope grism spectroscopy of 15 galaxies at z ∼ 0.8 drawn from the DEEP2 survey. We analyze H α +[N ii ], [S ii ], and [S iii ] emission on kiloparsec scales to explore which mechanisms are powering emission lines at high redshifts, testing which processes may be responsible for the well-known offset of high-redshift galaxies from the z ∼ 0 locus in the [O iii ]/H β versus [N ii ]/H α Baldwin—Phillips—Terlevich (BPT) excitation diagram. We study spatially resolved emission-line maps to examine evidence for active galactic nuclei (AGN), shocks, diffuse ionized gas (DIG), or escaping ionizing radiation, all of which may contribute to the BPT offsets observed in our sample. We do not find significant evidence of AGN in our sample and quantify that, on average, AGN would need to contribute ∼25% of the H α flux in the central resolution element in order to cause the observed BPT offsets. We find weak (2 σ ) evidence of DIG emission at low surface brightnesses, yielding an implied total DIG emission fraction of ∼20%, which is not significant enough to be the dominant emission line driver in our sample. Inmore »general we find that the observed emission is dominated by star-forming H ii regions. We discuss trends with demographic properties and the possible role of α -enhanced abundance patterns in the emission spectra of high-redshift galaxies. Our results indicate that photoionization modeling with stellar population synthesis inputs is a valid tool to explore the specific star formation properties which may cause BPT offsets, to be explored in future work.« less