Characterizing the Average Interstellar Medium Conditions of Galaxies at z ∼ 5.6–9 with Ultraviolet and Optical Nebular Lines
Weida Hu, Casey Papovich, Mark Dickinson, Robert C. Kennicutt, Lu Shen, R. Amorín, Pablo Arrabal Haro, Micaela B. Bagley, Rachana Bhatawdekar, Nikko J. Cleri, Justin W. Cole, Avishai Dekel, Alexander de la Vega, Steven L. Finkelstein, Norman A. Grogin, Nimish P. Hathi, Michaela Hirschmann, Benne W. Holwerda, Taylor A. Hutchison, Intae Jung, Anton M. Koekemoer, Jeyhan S. Kartaltepe, Ray A. Lucas, Mario Llerena, Sara Mascia, Bahram Mobasher, Lorenzo Napolitano, Jeffrey A. Newman, L. Pentericci, Pablo G. Pérez‐González, Jonathan R. Trump, Stephen M. Wilkins, L. Y. Aaron Yung
Abstract
Abstract Ultraviolet (UV; rest-frame ∼1200–2000 Å) spectra provide a wealth of diagnostics to characterize fundamental galaxy properties, such as their chemical enrichment, the nature of their stellar populations, and their amount of Lyman-continuum (LyC) radiation. In this work, we leverage publicly released JWST data to construct the rest-frame UV-to-optical composite spectrum of a sample of 63 galaxies at 5.6 < z < 9, spanning a wavelength range from 1500 to 5200 Å. Based on the composite spectrum, we derive an average dust attenuation <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mi>E</mml:mi> <mml:msub> <mml:mrow> <mml:mo stretchy="false">(</mml:mo> <mml:mi>B</mml:mi> <mml:mo>−</mml:mo> <mml:mi>V</mml:mi> <mml:mo stretchy="false">)</mml:mo> </mml:mrow> <mml:mrow> <mml:mi>gas</mml:mi> </mml:mrow> </mml:msub> <mml:mspace width="0.25em"/> <mml:mo>=</mml:mo> <mml:mspace width="0.25em"/> <mml:msubsup> <mml:mrow> <mml:mn>0.10</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.11</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.10</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> from H β /H γ , an electron density <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msub> <mml:mrow> <mml:mi>n</mml:mi> </mml:mrow> <mml:mrow> <mml:mi>e</mml:mi> </mml:mrow> </mml:msub> <mml:mo>=</mml:mo> <mml:msubsup> <mml:mrow> <mml:mn>570</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>290</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>510</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> cm −3 from the [O ii ] doublet ratio, an electron temperature <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msub> <mml:mrow> <mml:mi>T</mml:mi> </mml:mrow> <mml:mrow> <mml:mi>e</mml:mi> </mml:mrow> </mml:msub> <mml:mo>=</mml:mo> <mml:msubsup> <mml:mrow> <mml:mn>16700</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>1500</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>1500</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> K from the [O iii ] λ 4363/[O iii ] λ 5007 ratio, and an ionization parameter <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mi>log</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:mi>U</mml:mi> <mml:mo stretchy="false">)</mml:mo> <mml:mo>=</mml:mo> <mml:mo>−</mml:mo> <mml:msubsup> <mml:mrow> <mml:mn>2.15</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.03</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.03</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> from the [O iii ]/[O ii ] ratio. Using a direct T e method, we calculate an oxygen abundance <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mn>12</mml:mn> <mml:mo>+</mml:mo> <mml:mi>log</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:mi mathvariant="normal">O</mml:mi> <mml:mrow> <mml:mo stretchy="true">/</mml:mo> </mml:mrow> <mml:mi mathvariant="normal">H</mml:mi> <mml:mo stretchy="false">)</mml:mo> <mml:mspace width="0.25em"/> <mml:mo>=</mml:mo> <mml:mspace width="0.25em"/> <mml:mn>7.67</mml:mn> <mml:mo>±</mml:mo> <mml:mn>0.08</mml:mn> </mml:math> and a carbon-to-oxygen (C/O) abundance ratio <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mi>log</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:mi mathvariant="normal">C</mml:mi> <mml:mrow> <mml:mo stretchy="true">/</mml:mo> </mml:mrow> <mml:mi mathvariant="normal">O</mml:mi> <mml:mo stretchy="false">)</mml:mo> <mml:mspace width="0.25em"/> <mml:mo>=</mml:mo> <mml:mspace width="0.25em"/> <mml:mo>−</mml:mo> <mml:msubsup> <mml:mrow> <mml:mn>0.86</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.10</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.13</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> . This C/O ratio is smaller than compared to z = 0 and z = 2–4 star-forming galaxies, albeit with moderate significance. This indicates the reionization-era galaxies might be undergoing a rapid buildup of stellar mass with high specific star formation rates. A UV diagnostic based on the ratios of C iii ] λ λ 1907, 1909/He ii λ 1640 versus O iii ] λ 1666/He ii λ 1640 suggests that the star formation is the dominant source of ionization, similar to the local extreme dwarf galaxies and z ∼ 2–4 He ii –detected galaxies. The [O iii ]/[O ii ] and C iv /C iii ] ratios of the composite spectrum are marginally larger than the criteria used to select galaxies as LyC leakers, suggesting that some of the galaxies in our sample are strong contributors to the reionizing radiation.