Extreme equivalent width-selected low-mass starbursts at z=4-9: insights into their role in cosmic reionization

We investigate the properties of extreme emission line galaxies (EELGs) at $z=4-9$ and their role in reionization. Compact, low-mass galaxies with intense optical emission lines are linked to elevated specific star formation rates (sSFRs) and recent bursts of star formation. Feedback in these systems may enable the leakage of ionizing radiation into the intergalactic medium. Using JWST/NIRSpec spectroscopy from the CAPERS, CEERS, and RUBIES surveys, we compile 160 NIRCam-selected EELGs in the EGS field. These galaxies show extreme rest-frame equivalent widths (EWs), with a median EW([O III]+H${\beta}$)=1616\r{A} and EW(H${\alpha}$)=763\r{A}. They are low-mass (median log(M$_{\star}$/M$_{\odot}$)=8.26) with high sSFRs (median 43 Gyr$^{-1}$), above the $z\sim6$ main sequence. UV slopes are diverse, with a mean $\beta=-2.0$, and only 7% have extremely blue continua ($\beta<-2.6$). Emission-line diagnostics suggest stellar populations as the primary ionizing source, although an AGN fraction of 14% is cannot be entirely ruled out. These galaxies are efficient ionizing photon producers, with mean log($\xi_{\rm ion}$ [Hz erg$^{-1}$])=25.37, exceeding typical values at similar redshifts. Escape fractions, however, are heterogeneous: 16% of EELGs at $z<7$ show escape fractions $>$5% for both Ly${\alpha}$ and LyC photons, while 82% lack detectable Ly${\alpha}$ emission. The median inferred LyC escape fraction is modest (5%) but enhanced in compact super-Eddington systems with sSFR >25 Gyr$^{-1}$. These results indicate that EELGs contribute approximately 16-40% of the total ionizing emissivity required to sustain hydrogen reionization. EELGs are extremely compact, with a median effective radius of 0.49 kpc, and exhibit a recent star-formation burst. Our analysis indicates that sSFR and star-formation rate surface density are the primary drivers of their extreme emission line strengths.