2022ApJS..262...38L

Based on Sloan Digital Sky Survey Data Release 16, we have detected the large-scale structure of Lyα emission in the universe at redshifts z = 2-3.5 by cross-correlating quasar positions and Lyα emission imprinted in the residual spectra of luminous red galaxies. We apply an analytical model to fit the corresponding Lyα surface brightness profile and multipoles of the redshift-space quasar-Lyα emission cross-correlation function. The model suggests an average cosmic Lyα luminosity density of ${6.6}_{-3.1}^{+3.3}\times {10}^{40}\,\mathrm{erg}\,{{\rm{s}}}^{-1}\,{\mathrm{cMpc}}^{-3}$, a ∼2σ detection with a median value about 8-9 times those estimated from deep narrowband surveys of Lyα emitters at similar redshifts. Although the low signal-to-noise ratio prevents us from a significant detection of the Lyα forest-Lyα emission cross-correlation, the measurement is consistent with the prediction of our best-fit model from quasar-Lyα emission cross-correlation within current uncertainties. We rule out the scenario where the Lyα photons mainly originate from quasars. We find that Lyα emission from star-forming galaxies, including contributions from that concentrated around the galaxy centers and that in diffuse Lyα-emitting halos, is able to explain the bulk of the Lyα luminosity density inferred from our measurements. Ongoing and future surveys can further improve the measurements and advance our understanding of the cosmic Lyα emission field.