The Physical Origin of Extreme Emission Line Galaxies at High redshifts: Strong {sc [Oiii]} Emission Lines Produced by Obscured AGNs

We present deep Subaru/FOCAS spectra for two extreme emission line galaxies (EELGs) at $z\sim 1$ with strong {\sc[Oiii]}$\lambda$5007 emission lines, exhibiting equivalent widths (EWs) of $2905^{+946}_{-578}$ \AA\ and $2000^{+188}_{-159}$ \AA, comparable to those of EELGs at high redshifts that are now routinely identified with JWST spectroscopy. Adding a similarly large {\sc [Oiii]} EW ($2508^{+1487}_{-689}$ \AA) EELG found at $z\sim 2$ in the JWST CEERS survey to our sample, we explore for the physical origins of the large {\sc [Oiii]} EWs of these three galaxies with the Subaru spectra and various public data including JWST/NIRSpec, NIRCam, and MIRI data. While there are no clear signatures of AGN identified by the optical line diagnostics, we find that two out of two galaxies covered by the MIRI data show strong near-infrared excess in the spectral energy distributions (SEDs) indicating obscured AGN. Because none of the three galaxies show clear broad H$\beta$ lines, the upper limits on the flux ratios of broad-H$\beta$ to {\sc [Oiii]} lines are small, $\lesssim 0.15$ that are comparable with Seyfert $1.8-2.0$ galaxies. We conduct \texttt{Cloudy} modeling with the stellar and AGN incident spectra, allowing a wide range of parameters including metallicities and ionization parameters. We find that the large {\sc [Oiii]} EWs are not self-consistently reproduced by the spectra of stars or unobscured AGN, but obscured AGN that efficiently produces O$^{++}$ ionizing photons with weak nuclear and stellar continua that are consistent with the SED shapes.