skip to main content


Title: The effects of binary stars on galaxies and metal-enriched gas during reionization
ABSTRACT Binary stars are abundant in nearby galaxies, but are typically unaccounted for in simulations of the high-redshift Universe. Stellar population synthesis models that include the effects of binary evolution result in greater relative abundances of ionizing photons that could significantly affect the ambient ionizing background during the epoch of hydrogen reionization, additionally leading to differences in galaxy gas content and star formation. We use hydrodynamic cosmological simulations including in situ multifrequency radiative transfer to evaluate the effects of a high binary fraction in reionization-era galaxies on traits of the early intergalactic medium and the abundance of H i and He ii ionizing photons. We further extend this to analyse the traits of enriched gas. In comparing metrics generated using a fiducial simulation assuming single stars with one incorporating a high binary fraction, we find that binary stars cause H i reionization to complete earlier and at an accelerated pace, while also increasing the abundances of high-ionization metals (C iv and Si iv) in simulated absorption spectra while reducing the abundance of low-ionization states (O i, Si ii, and C ii). However, through increased photoheating of galactic and circumgalactic gas, they simultaneously reduce the rate of star formation in low-mass galaxies, slowing the ongoing process of enrichment and suppressing their own ionizing background. This potentially contributes to a slower He ii reionization process at $z\ge 5$, and further indicates that self-regulation of galaxies could be underestimated when neglecting binary stellar evolution.  more » « less
Award ID(s):
2006550
NSF-PAR ID:
10340071
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
505
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2207 to 2223
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We present an updated model of the cosmic ionizing background from the UV to the X-rays. Relative to our previous model, the new model provides a better match to a large number of up-to-date empirical constraints, including: (1) new galaxy and AGN luminosity functions; (2) stellar spectra including binary stars; (3) obscured and unobscured AGN; (4) a measurement of the non-ionizing UV background; (5) measurements of the intergalactic H i and He ii photoionization rates at z ∼ 0−6; (6) the local X-ray background; and (7) improved measurements of the intergalactic opacity. In this model, AGN dominate the H i ionizing background at z ≲ 3 and star-forming galaxies dominate it at higher redshifts. Combined with the steeply declining AGN luminosity function beyond z ∼ 2, the slow evolution of the H i ionization rate inferred from the high-redshift H i Ly α forest requires an escape fraction from star-forming galaxies that increases with redshift (a population-averaged escape fraction of $\approx 1{{\ \rm per\ cent}}$ suffices to ionize the intergalactic medium at z = 3 when including the contribution from AGN). We provide effective photoionization and photoheating rates calibrated to match the Planck 2018 reionization optical depth and recent constraints from the He ii Ly α forest in hydrodynamic simulations. 
    more » « less
  2. ABSTRACT

    The ionizing photon escape fraction [Lyman continuum (LyC) fesc] of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshifts. Here we circumvent this uncertainty by inferring fesc using resolved (R > 4000) Lyman α (Lyα) profiles from the X-SHOOTER Lyα survey at z = 2 (XLS-z2). With empirically motivated criteria, we use Lyα profiles to select leakers ($f_{\mathrm{ esc}} > 20{{\ \rm per\ cent}}$) and non-leakers ($f_{\mathrm{ esc}} < 5{{\ \rm per\ cent}}$) from a representative sample of >0.2L* Lyman α emitters (LAEs). We use median stacked spectra of these subsets over λrest ≈ 1000–8000 Å to investigate the conditions for LyC fesc. Our stacks show similar mass, metallicity, MUV, and βUV. We find the following differences between leakers versus non-leakers: (i) strong nebular C iv and He ii emission versus non-detections; (ii) [O iii]/[O ii] ≈ 8.5 versus ≈3; (iii) Hα/Hβ indicating no dust versus E(B − V) ≈ 0.3; (iv) Mg ii emission close to the systemic velocity versus redshifted, optically thick Mg ii; and (v) Lyα fesc of ${\approx} 50{{\ \rm per\ cent}}$ versus ${\approx} 10{{\ \rm per\ cent}}$. The extreme equivalent widths (EWs) in leakers ([O iii]+$\mathrm{ H}\beta \approx 1100$ Å rest frame) constrain the characteristic time-scale of LyC escape to ≈3–10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers – extremely ionizing stellar populations, low column densities, a dust-free, high-ionization state interstellar medium (ISM) – occur simultaneously in the $f_{\rm esc} > 20{{\ \rm per\ cent}}$ stack, suggesting they are causally connected, and motivating why indicators like [O iii]/[O ii] may suffice to constrain fesc at z > 6 with the James Webb Space Telescope (JWST). The leakers comprise half of our sample, have a median LyC$f_{\rm esc} \approx 50{{\ \rm per\ cent}}$ (conservative range: $20\!-\!55{{\ \rm per\ cent}}$), and an ionizing production efficiency $\log ({\xi _{\rm {ion}}/\rm {Hz\ erg^{-1}}})\approx 25.9$ (conservative range: 25.7–25.9). These results show LAEs – the type of galaxies rare at z ≈ 2, but that become the norm at higher redshift – are highly efficient ionizers, with extreme ξion and prolific fesc occurring in sync.

     
    more » « less
  3. Abstract Observations of high-redshift galaxies ( z > 5) have shown that these galaxies have extreme emission lines with equivalent widths much larger than their local star-forming counterparts. Extreme emission line galaxies (EELGs) in the nearby universe are likely analogs to galaxies during the Epoch of Reionization and provide nearby laboratories to understand the physical processes important to the early universe. We use Hubble Space Telescope/Cosmic Origins Spectrograph and Large Binocular Telescope/Multi-Object Double Spectrographs spectra to study two nearby EELGs, J104457 and J141851. The far-UV spectra indicate that these two galaxies contain stellar populations with ages ≲10 Myr and metallicities ≤0.15 Z ⊙ . We use photoionization modeling to compare emission lines from models of single-age bursts of star formation to observed emission lines and find that the single-age bursts do not reproduce high-ionization lines including [O iii ] or very-high-ionization lines like He ii or O iv ]. Photoionization modeling using the stellar populations fit from the UV continuum similarly is not capable of reproducing the very-high-energy emission lines. We add a blackbody to the stellar populations fit from the UV continuum to model the necessary high-energy photons to reproduce the very-high-ionization lines of He ii and O iv ]. We find that we need a blackbody of 80,000 K and ∼45%–55% of the luminosity from the blackbody and young stellar population to reproduce the very-high-ionization lines while simultaneously reproducing the low-, intermediate-, and high-ionization emission lines. Our self-consistent model of the ionizing spectra of two nearby EELGs indicates the presence of a previously unaccounted-for source of hard ionizing photons in reionization analogs. 
    more » « less
  4. Abstract

    High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (fesc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is ≳70%. Moreover, simulations so far have yielded values offescfalling only on the lower end of the ∼10%–20% range, the amount presumed necessary to reionize the universe. Analyzing a high-resolution (4 pc) cosmological radiation-hydrodynamic simulation, we study howfescchanges when we include two different products of binary stellar evolution—stars stripped of their hydrogen envelopes and massive blue stragglers. Both produce significant amounts of ionizing photons 10–200 Myr after each starburst. We find the relative importance of these photons to be amplified with respect to escaped ionizing photons, because peaks in star formation rates (SFRs) andfescare often out of phase by this 10–200 Myr. Additionally, low-mass, bursty galaxies emit Lyman continuum radiation primarily from binary products when SFRs are low. Observations of these galaxies by the James Webb Space Telescope could provide crucial information on the evolution of binary stars as a function of redshift. Overall, including stripped stars and massive blue stragglers increases our photon-weighted mean escape fraction () by ∼13% and ∼10%, respectively, resulting in. Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization.

     
    more » « less
  5. ABSTRACT

    Nearly a decade ago, we began to see indications that reionization-era galaxies power hard radiation fields rarely seen at lower redshift. Most striking were detections of nebular C iv emission in what appeared to be typical low-mass galaxies, requiring an ample supply of 48 eV photons to triply ionize carbon. We have obtained deep JWST/NIRSpec R = 1000 spectroscopy of the two z > 6 C iv-emitting galaxies known prior to JWST. Here, we present a rest-UV to optical spectrum of one of these two systems, the multiply-imaged z = 6.1 lensed galaxy RXCJ2248-ID. NIRCam imaging reveals two compact (<22 pc) clumps separated by 220 pc, with one comprising a dense concentration of massive stars (>10 400 M⊙ yr−1 kpc−2) formed in a recent burst. We stack spectra of 3 images of the galaxy (J = 24.8–25.9), yielding a very deep spectrum providing a high-S/N template of strong emission line sources at z > 6. The spectrum reveals narrow high-ionization lines (He ii, C iv, N iv]) with line ratios consistent with powering by massive stars. The rest-optical spectrum is dominated by very strong emission lines ([O iii] EW = 2800 Å), albeit with weak emission from low-ionization transitions ([O iii]/[O ii] = 184). The electron density is found to be very high (6.4–31.0 × 104 cm−3) based on three UV transitions. The ionized gas is metal poor ($12+\log (\rm O/H)=7.43^{+0.17}_{-0.09}$), yet highly enriched in nitrogen ($\log (\rm N/O)=-0.39^{+0.11}_{-0.10}$). The spectrum appears broadly similar to that of GNz11 at z = 10.6, without showing the same AGN signatures. We suggest that the hard radiation field and rapid nitrogen enrichment may be a short-lived phase that many z > 6 galaxies go through as they undergo strong bursts of star formation. We comment on the potential link of such spectra to globular cluster formation.

     
    more » « less