Excitation and Ionization Properties of Star-forming Galaxies at z = 2.0-9.3 with JWST/NIRSpec

We utilize medium-resolution JWST/NIRSpec observations of 164 galaxies at z = 2.0-9.3 from the Cosmic Evolution Early Release Science (CEERS) survey to investigate the evolution of the excitation and ionization properties of galaxies at high redshifts. Our results represent the first statistical constraints on the evolution of the [O III]/Hβ versus [N II]/Hα, [S II]/Hα, and [O I]/Hα "BPT" diagrams at z > 2.7, and the first analysis of the O₃₂ versus R₂₃ diagram at z > 4 with a large sample. We divide the sample into five redshift bins containing 30-40 galaxies each. The subsamples at z ~ 2.3, z ~ 3.3, and z ~ 4.5 are representative of the main-sequence star-forming galaxy population at these redshifts, while the z ~ 5.6 and z ~ 7.5 samples are likely biased toward high specific star formation rate, due to selection effects. Using composite spectra, we find that each subsample at z = 2.0-6.5 falls on the same excitation sequence in the [N II] and [S II] BPT diagrams and the O₃₂-R₂₃ diagram on average, and is offset from the sequences followed by z = 0 H II regions in the same diagrams. The direction of these offsets are consistent with high-redshift star-forming galaxies uniformly having harder ionizing spectra than typical local galaxies at fixed nebular metallicity. The similarity of the average line ratios suggests that the ionization conditions of the interstellar medium do not strongly evolve between z ~ 2 and z ~ 6. Overall, the rest-optical line ratios suggest the z = 2.7-9.3 CEERS/NIRSpec galaxies at log(M _*/M _⊙) ~ 7.5-10 have high degrees of ionization and moderately low oxygen abundances (~0.1-0.3 Z _⊙), but are not extremely metal-poor (<0.1 Z _⊙) even at z > 6.5.