---
title: Extreme Emission Line Galaxies (EELGs)
url: https://www.emergentmind.com/topics/extreme-emission-line-galaxies-eelgs
type: topic
---

# Extreme Emission Line Galaxies (EELGs)

Extreme emission line galaxies (EELGs) are galaxies whose spectra and broad-band spectral energy distributions are dominated by nebular emission lines with unusually large equivalent widths, most commonly in \([\mathrm{O\,III}] + \mathrm{H}\beta\) and \(\mathrm{H}\alpha\). Across local, intermediate-redshift, and JWST high-redshift samples, they are consistently described as compact, low-mass, intensely star-forming, low-metallicity systems with high ionization conditions and bursty star-formation histories, and they are widely used as local or intermediate-redshift analogues of galaxies in the epoch of reionization [2311.02158] [2404.17415] [2304.11181] [2604.09516].

## 1. Operational definitions and taxonomic scope

There is no universal equivalent-width threshold for EELGs in the literature. Instead, different surveys adopt thresholds matched to their bandpasses, spectral resolution, and scientific goals [2404.17415]. In low- and intermediate-redshift spectroscopic work, a common operational choice is an extreme \([\mathrm{O\,III}]\lambda5007\) equivalent width. In the zCOSMOS-bright survey, EELGs were defined by \(\mathrm{EW}([\mathrm{O\,III}]\,\lambda5007) \ge 100\) Å over \(0.11 \le z \le 0.93\), a cut that automatically selected galaxies with very high Balmer-line equivalent widths and very young bursts [1403.3441]. In SDSS-based stellar-population work, the class was identified with \(\mathrm{EW}(\mathrm{H}\beta) > 30\) Å and \(\mathrm{EW}([\mathrm{O\,III}]\,\lambda5007) > 100\) Å [2205.07660]. In MUSE HUDF spectroscopy, the primary sample required \(\mathrm{EW}_0 \ge 300\) Å and the extended sample \(200 \le \mathrm{EW}_0 \le 300\) Å in at least one of \([\mathrm{O\,II}]\), \([\mathrm{O\,III}]\), or \(\mathrm{H}\alpha\) [2404.17415]. In J-PLUS, the local sample was defined by \(\mathrm{EW}([\mathrm{OIII}]4959+5007) > 300\) Å at \(z<0.06\) [2112.06938].

In the early-universe context, the definition is often expressed photometrically rather than by a single line. One JWST-based formulation describes EELGs as systems where nebular emission contributes \(30\)–\(40\%\) of the flux in certain photometric bands, with rest-frame \([\mathrm{O\,III}] + \mathrm{H}\beta\) equivalent widths of several hundred to \(>800\) Å [2311.02158]. In CEERS, the photometric definition was moved to the observed frame: \(\mathrm{EW}_{\rm phot} > 5000\) Å in either \(\mathrm{H}\alpha\) or \(\mathrm{H}\beta + [\mathrm{O\,III}]\) [2602.23310].

| Study | Operational definition | Context |
|---|---|---|
| zCOSMOS [1403.3441] | \(\mathrm{EW}([\mathrm{O\,III}]\,\lambda5007) \ge 100\) Å | \(0.11<z<0.93\) spectroscopic sample |
| MUSE HUDF [2404.17415] | Primary: \(\mathrm{EW}_0 \ge 300\) Å; extended: \(200 \le \mathrm{EW}_0 \le 300\) Å | \([\mathrm{O\,II}]\), \([\mathrm{O\,III}]\), or \(\mathrm{H}\alpha\) |
| J-PLUS [2112.06938] | \(\mathrm{EW}([\mathrm{OIII}]4959+5007) > 300\) Å | Local \(z<0.06\) photometric sample |
| MOSEL/JADES [2311.02158] | Nebular emission contributes \(30\)–\(40\%\) of band flux; \([\mathrm{O\,III}] + \mathrm{H}\beta\) EWs of several hundred to \(>800\) Å | Early-universe/JWST context |
| CEERS [2602.23310] | \(\mathrm{EW}_{\rm phot} > 5000\) Å observed frame | \(\mathrm{H}\alpha\) or \(\mathrm{H}\beta + [\mathrm{O\,III}]\), \(4<z<9\) |

This diversity of criteria does not imply a fragmented class. A plausible implication is that “EELG” functions as a physically coherent label for galaxies in a short-lived, line-dominated starburst phase, while the exact boundary depends on the survey and line complex used.

## 2. Survey strategies and selection methodologies

EELG samples are now assembled by a heterogeneous set of methods, each with distinct strengths. Broad-band excess selections target the line-dominated filter directly. At \(3<z<3.7\), one method identified EELGs through a \(K_s\)-band excess relative to the best-fit stellar continuum, using \(\Delta m_{Ks} > 0.3\) mag as the primary criterion; this corresponds to \(\mathrm{EW}_{\rm obs}([\mathrm{O\,III}] + \mathrm{H}\beta) \gtrsim 1000\) Å [2010.07545]. Follow-up MOIRCS spectroscopy showed that \(21/23\) targets had clear emission lines and \(19/21\) were the intended \([\mathrm{O\,III}]\) emitters at \(z>3\), validating the broad-band excess strategy [2010.07545].

Medium- and narrow-band selections improve spectral localization. J-PLUS identifies local EELGs through an excess in the \(J0515\) filter, which covers \([\mathrm{O\,III}]\) at \(z<0.06\), and its final sample contains 466 galaxies with \([\mathrm{O\,III}]\) equivalent widths above 300 Å; the method reaches \(\approx 92\%\) completeness and \(\approx 96\%\) purity against SDSS spectroscopy [2112.06938]. J-PAS extends this logic to 56 optical bands. In a fully observed \(30\,\mathrm{deg}^2\) region, a photometric method combining narrow-band equivalent widths with machine learning identified 917 EELGs up to \(z=0.8\), with \(95\%\) purity and \(96\%\) completeness for \(i_{\rm SDSS}<22.5\) mag [2512.08484]. miniJPAS demonstrates the same principle in pilot form, using a contrast criterion corresponding to rest-frame equivalent width \(\gtrsim 300\) Å and a fully automated multiwavelength pipeline for source characterization [2401.13816].

Integral-field spectroscopy provides purely spectroscopic selection and spatial information simultaneously. In the MUSE HUDF survey, the search proceeded directly from rest-frame equivalent widths measured with pyPlatefit, followed by object-by-object inspection and remeasurement of line fluxes and EWs; the final sample comprises 13 EELGs at \(0.1<z<0.9\) [2404.17415]. This removes dependence on broad-band color preselection and yields direct kinematic and chemical diagnostics from the same data cube.

JWST medium-band photometry has generalized the idea of “spectroscopy from photometry.” In CANUCS, NIRCam medium bands were used to select 118 EELGs over \(1.7 \lesssim z \lesssim 6.7\), with median \(EW(\mathrm{H}\alpha)=893\) Å and median \(EW([\mathrm{OIII}] + \mathrm{H}\beta)=1255\) Å; NIRSpec spectroscopy of 15 objects confirmed the redshifts and equivalent widths derived from the medium bands [2304.11181]. In CEERS, a simple photometric method identified 1165 EELGs at \(4<z<9\), and NIRSpec spectroscopy of 34 objects validated the photometric identification of extreme emission; the medium-band F410M filter was found to be particularly efficient for isolating these systems [2312.07799].

At the largest scale, DESI has enabled a classification-driven rather than line-threshold-driven census. The k-MENDEL sample used automatic \(k\)-means classification of DESI spectra to isolate rare, high-EW spectral classes, yielding 15,014 EELGs at \(0.01<z<0.96\) after quality cuts; the selected classes correspond almost entirely to objects with \(\mathrm{EW}_{5007} \gtrsim 100\) Å, while the largest-\(\mathrm{EW}_{5007}>500\) Å objects occur only in the outlier classes [2604.09516].

## 3. Stellar populations, chemical abundances, and ionization conditions

Across these studies, EELGs occupy the low-mass, high-sSFR, metal-poor end of the star-forming population. In the MUSE HUDF sample, stellar masses lie in the dwarf regime, with \(\log M_\star/M_\odot \sim 6.85\)–8.60 and \(\log \mathrm{SFR}\sim -2.4\) to \(-0.5\), implying very high specific star-formation rates [2404.17415]. The zCOSMOS sample spans \(10^7 \lesssim M_\star \lesssim 10^{10}\,M_\odot\) with median \(\log(M_\star/M_\odot)\approx 8.8\), and its specific star-formation rates reach up to \(\sim 10^{-7}\,\mathrm{yr}^{-1}\) [1403.3441]. The DESI k-MENDEL sample broadens the demographic baseline to \(M_\star \sim 10^{6}-10^{10}\,M_\odot\), \(\mathrm{SFR}\sim 0.1-100\,M_\odot\,\mathrm{yr}^{-1}\), and \(\mathrm{sSFR}\sim 10-100\,\mathrm{Gyr}^{-1}\), with EELGs systematically above the local star-forming main sequence [2604.09516]. In CANUCS, the median stellar mass is \(\log(M_\star/M_\odot)=8.03\), with median metallicity \(Z=0.14\,Z_\odot\), median \(A_V=0.18\) mag, and median \(\mathrm{SSFR}=1.18\times10^{-8}\,\mathrm{yr}^{-1}\) [2304.11181].

Gas-phase metallicities are consistently subsolar and often extremely low. The MUSE HUDF sample has \(7.35 \lesssim 12+\log(\mathrm{O/H}) \lesssim 8.05\), with four galaxies below \(0.1\,Z_\odot\) [2404.17415]. The zCOSMOS EELGs have median \(12+\log(\mathrm{O/H})=8.16\), including a handful of extremely metal-deficient systems below \(10\%\) solar [1403.3441]. In the DESI sample, direct-\(T_e\) metallicities span \(12+\log(\mathrm{O/H}) \sim 6.9-8.5\), with a median near 7.85, while the full strong-line scale extends to \(\sim 8.6\) [2604.09516]. Dust attenuation is typically low: the MUSE HUDF sample has \(A_V \approx 0\)–0.45 mag [2404.17415], the \(z\simeq 3.3\) broad-band-selected sample has \(E(B-V)_\mathrm{star}\lesssim 0.1\) mag [2010.07545], and the J-PLUS local sample has \(E(B-V)\sim 0.1^{+0.2}_{-0.1}\) [2112.06938].

A defining feature of EELGs is their high-ionization nebular spectrum. In the \(z\simeq 3.3\) sample, \(O_{32}\) values are typically \(\gtrsim 2\)–3, with several sources at \(O_{32}\gtrsim 5\) and lower limits as high as \(\gtrsim 7\); the authors infer \(12+\log(\mathrm{O/H})\simeq 7.5\)–8.5 and ionization parameters higher by \(\gtrsim 1.5\) dex than local star-forming galaxies [2010.07545]. The ionization parameter is conventionally written as
$$
U \equiv \frac{Q_{\rm H}}{4\pi R^2 n_{\rm H} c},
$$
and interaction-driven starbursts can plausibly raise \(Q_{\rm H}\) while compressing the gas, thereby increasing \(U\) and line equivalent widths [2311.02158].

At the most extreme end, EELGs exhibit a very-high-ionization component that exceeds the reach of standard three-zone H II region models. Detailed UV and optical spectroscopy of two nearby archetypes shows strong nebular He II, C IV, \([\mathrm{Fe\,V}]\), and \([\mathrm{Ar\,IV}]\), motivating a four-zone ionization model with an added He\(^{+2}\) zone above 54.4 eV [2105.12765]. In that framework, traditional three-zone estimates can under-estimate the average \(\log U\) by up to 0.5 dex, while the total nebular abundances remain nearly unchanged [2105.12765]. The same study identifies a model-independent abundance dichotomy in which \(\alpha/\mathrm{H}\) abundances are consistent but N/H, C/H, and Fe/H are relatively deficient, implying \(\alpha/\mathrm{Fe}\) enhancement by \(>3\) times [2105.12765]. It also concludes that current photoionization models still cannot reproduce the observed very-high-ionization lines, leaving a high-energy ionizing photon production problem unresolved [2105.12765].

## 4. Structure, kinematics, and triggering mechanisms

EELGs are usually compact, but their morphology is not uniform. In zCOSMOS, the median half-light radius is \(R_{50}\sim 1.3\) kpc, and \(\sim 80\%\) of the galaxies show non-axisymmetric morphologies, including clumpy and tadpole systems; \(\sim 29\%\) show additional low-surface-brightness features that strongly suggest recent or ongoing interactions [1403.3441]. The local J-PLUS sample similarly finds that most objects have compact morphologies, while \(\approx 19\%\) are more extended dwarfs with clumpy or disturbed structure, underscoring that “extreme [O III]” denotes a physical state rather than a single morphological class [2112.06938].

Integral-field observations show that compactness does not imply kinematic simplicity. In the MUSE HUDF EELGs, four of the seven primary-sample galaxies and five of the six extended-sample galaxies are spatially resolved; three of the resolved primary EELGs show a clear rotating-disk pattern, whereas none of the resolved extended-sample objects do [2404.17415]. The same study cautions that mergers can mimic rotation at MUSE resolution, so higher-resolution IFU data are required to separate rotating disks, unresolved mergers, and outflows [2404.17415].

Environmental work has increasingly linked EELGs to interactions. The MOSEL survey used JADES NIRCam imaging to examine spectroscopically confirmed strong \([\mathrm{O\,III}]\) emitters at \(2.5<z<4\). Compared with control galaxies, the EELGs show a median brightest-companion mass ratio \(R_{\rm max}=0.65^{+0.11}_{-0.51}\) and total companion mass ratio \(R_\Sigma=0.71^{+0.27}_{-0.34}\), versus \(R_{\rm max}=0.05^{+0.05}_{-0.02}\) and \(R_\Sigma=0.07^{+0.04}_{-0.04}\) in the full control sample, with KS-test probabilities \(p=0.007\) and \(p=0.006\) [2311.02158]. Even after matching in both stellar mass and specific SFR, EELGs retain \(\approx 3\)–4 times higher brightest-companion mass ratios and \(\approx 3\) times higher total companion mass ratios than the matched controls [2311.02158]. TNG100 tests show that with the same projected-separation and velocity cuts, \(99.8\%\) of galaxies have at least one companion sharing the same descendant by \(z=0\), supporting the interpretation that these close massive companions represent real mergers or strong interactions rather than chance projections [2311.02158]. A plausible implication is that interaction-driven gas cooling and inflow are major triggers of the EELG phase at \(z\sim 3\).

The role of active galactic nuclei is now better constrained in the JWST era. A spectroscopic study of CEERS EELGs with observed-frame EW \(>5000\) Å finds that \(\sim 10\%\) of photometrically selected EELGs have broad Balmer lines, increasing to \(35\%\) in deep spectroscopy, so AGN are not negligible within the class [2602.23310]. However, many AGN selected photometrically as EELGs have incorrectly high photometric equivalent widths, and in sources that remain true high-EW EELGs spectroscopically, the narrow \(\mathrm{H}\alpha\) component dominates over the broad component; the same study concludes that CEERS EELGs are powered primarily by star formation rather than AGN [2602.23310]. It also finds that \([\mathrm{O\,III}]\) morphology becomes more compact at higher equivalent width, consistent with increasingly centralized ionizing sources in the most extreme systems [2602.23310].

## 5. High-redshift analogues, reionization, and contamination of galaxy selection

One of the principal scientific uses of EELGs is as analogues of the galaxies that dominated the early Universe. In JWST observations at \(z>6\), \(\sim 80\%\) of galaxies have rest-frame \([\mathrm{O\,III}] + \mathrm{H}\beta\) equivalent widths above 800 Å, almost three times the equivalent width of a typical star-forming galaxy at \(z\sim 2\) [2311.02158]. The HUDF EELGs occupy the same loci as \(z\approx 6\)–8 galaxies in the mass–metallicity relation, lying well below the local \(z\approx 0\) relation and close to the \(z=2.2\) MOSDEF and JWST high-redshift sequences [2404.17415]. The DESI k-MENDEL sample similarly follows a shallower mass–metallicity relation offset by \(0.3\)–0.5 dex from local relations and closely resembling young galaxies observed with JWST at \(z>3\)–10 [2604.09516]. Local and low-redshift photometric samples make the same point from the opposite direction: J-PLUS and J-PAS explicitly frame their EELGs as nearby analogues of reionization-era galaxies, and in J-PAS most sources exceed the canonical ionizing-efficiency threshold often associated with sustaining reionization [2112.06938] [2512.08484].

Ionizing-photon production efficiency is one of the clearest bridges between low-redshift EELGs and high-redshift galaxies. In the \(z\simeq 3.3\) broad-band-selected sample, \(\log \xi_{\rm ion,0}\) is typically \(25.1\)–25.9, with a median around 25.5 and a clear positive correlation with \(\mathrm{EW}([\mathrm{O\,III}])\); galaxies with \(\mathrm{EW}([\mathrm{O\,III}])\gtrsim 1000\) Å typically have \(\log \xi_{\rm ion,0}\gtrsim 25.5\) [2010.07545]. J-PAS finds the same qualitative behavior, with
$$
\log(\xi_{\rm ion}) = (0.69 \pm 0.15)\,\log(\mathrm{EW}_0[\mathrm{O\,III}]) + (23.35 \pm 0.06),
$$
and a significant fraction of its EELGs above the often-quoted \(\log \xi_{\rm ion}\gtrsim 25.3\) threshold [2512.08484]. This suggests that the EELG phase is closely tied to high ionizing efficiency.

The same properties that make EELGs valuable analogues also make them problematic contaminants in photometric searches for the highest-redshift galaxies. CLASH showed that EELGs with rest-frame \([\mathrm{O\,III}] + \mathrm{H}\beta\) equivalent widths \(\gtrsim 600\) Å, and in some cases \(\sim 3300\)–3700 Å, can mimic dropout colors when broad-band coverage is incomplete or signal-to-noise is limited [1412.7909]. That study identified 52 candidates in cluster-lensed HST fields and concluded that the fraction of EELGs in future high-redshift galaxy selections cannot be neglected [1412.7909]. JWST has sharpened the issue rather than removed it: the CEERS census reports examples of EELGs that could be incorrectly classified at ultra-high redshift, \(z>12\), owing to extreme \(\mathrm{H}\beta + [\mathrm{O\,III}]\) emission blended across the reddest filters [2312.07799]. Medium-band imaging mitigates the problem. In CEERS, F410M is especially effective because it captures the line complex while bracketing filters better constrain the continuum [2312.07799]; in CANUCS, medium-band color selection is explicitly advantageous because it selects by equivalent width rather than requiring strong continuum emission, and therefore recovers faint-continuum or red-continuum EELGs that broad-band selections can misclassify [2304.11181].

## 6. Modeling challenges, systematics, and open problems

The strongest methodological lesson from EELG research is that nebular emission cannot be treated as a perturbation. Self-consistent spectral modeling of 414 SDSS EELGs shows that stellar mass and mean age estimates from STARLIGHT are systematically biased toward higher values, and that adequate recovery of EELG stellar properties is only possible when nebular continuum emission is included [2205.07660]. The same paper finds that the discrepancies between stellar-only and stellar-plus-nebular synthesis correlate with specific SFR and with the summed flux of the strongest emission lines, implying that the systematic error is itself a function of star-formation intensity [2205.07660]. This point also appears in broad-band SED work: the zCOSMOS analysis found that ignoring nebular line contamination in photometry would overestimate \(M_\star\) by a median of \(\approx 0.25\) dex and by factors of 3–5 in the most extreme cases [1403.3441].

Selection systematics remain substantial. In the MOSEL/JADES environment study, companion identification depends on photometric redshifts with \(\sigma_{\rm NMAD}=0.024\) and \(\sim 13\%\) catastrophic outliers, so the fiducial \(|\Delta v|<10{,}000\ \mathrm{km\,s^{-1}}\) cut is necessarily broad [2311.02158]. In the MUSE HUDF sample, continuum estimation is intrinsically difficult because EELGs have very faint continua; varying the continuum window changed at least one object from above-threshold to borderline status, showing how sensitive equivalent widths are to small continuum errors [2404.17415]. In J-PAS, the quoted \(\xi_{\rm ion}\) values depend on the assumed attenuation law and on the choice \(f_{\rm esc}=0\); the authors note that an SMC-like law would systematically increase \(\xi_{\rm ion}\) by up to \(\sim 0.3\) dex, and that the reported values are lower limits if ionizing photons escape [2512.08484]. The CEERS AGN study adds a further caveat: broad-line AGN can acquire spuriously large photometric equivalent widths when continuum slopes are misestimated, so photometric “ultra-EELG” samples require spectroscopic vetting [2602.23310].

Beyond measurement systematics, there are unresolved physical problems. In the DESI sample, the mass–metallicity relation shows large intrinsic scatter and the scatter is not reduced by projection along the fundamental metallicity relation, indicating strong departures from simple equilibrium “bathtub” models and suggesting stochastic metal-poor inflows plus strong feedback [2604.09516]. In the nearby two-object study, even after invoking \(\alpha/\mathrm{Fe}\) enhancement and very high \(\log U\), photoionization models still fail to reproduce the observed very-high-ionization lines, leaving the high-energy ionizing photon production problem open [2105.12765]. This suggests that the hardest ionizing radiation in EELGs may require ingredients beyond conventional stellar population models, or at least beyond their current calibration in the low-metallicity, high-sSFR regime.

Taken together, these results define EELGs not as a narrowly bounded phenomenological class but as a recurrent, non-equilibrium phase in galaxy evolution. The phase is marked by very large nebular equivalent widths, compact and often centrally concentrated star formation, low metallicity, hard radiation fields, and in at least some samples a strong association with interactions or rapid gas accretion. It is now observed from the nearby Universe to \(z\gtrsim 9\), and the convergence of DESI, J-PAS, MUSE, and JWST results suggests that EELGs are simultaneously a local laboratory, a high-redshift selection challenge, and a key empirical route to the physics of the galaxies that dominated the first billion years of cosmic history.

Source: https://www.emergentmind.com/topics/extreme-emission-line-galaxies-eelgs