Water vapor in nearby infrared galaxies as probed by Herschel

Abstract: We report the first systematic study of the submillimeter water vapor rotational emission lines in infrared (IR) galaxies based on the Fourier Transform Spectrometer (FTS) data of {\it Herschel} SPIRE. Among the 176 galaxies with publicly available FTS data, 45 have at least one H$2$O emission line detected. The H$_2$O line luminosities range from $\sim 1 \times 10^5$ $L{\odot}$ to $\sim 5 \times 10^7 L_{\odot}$ while the total IR luminosities ($L_\mathrm{IR}$) have a similar spread ($\sim 1-300 \times 10^{10} L_{\odot}$). In addition, emission lines of H$2$O$^+$ and H$_2^{18}$O are also detected. H$_2$O is found, for most galaxies, to be the strongest molecular emitter after CO in FTS spectra. The luminosity of the five most important H$_2$O lines is near-linearly correlated with $L\mathrm{IR}$, no matter whether strong active galactic nucleus signature is present or not. However, the luminosity of H$2$O($2{11}-2_{02}$) and H$2$O($2{20}-2_{11}$) appears to increase slightly faster than linear with $L_\mathrm{IR}$. Although the slope turns out to be slightly steeper when $z\sim 2-4$ ULIRGs are included, the correlation is still closely linear. We find that $L_\mathrm{H_2O}/L_\mathrm{IR}$ decreases with increasing $f_{25}/f_{60}$, but see no dependence on $f_{60}/f_{100}$, possibly indicating that very warm dust contributes little to the excitation of the submillimeter H$2$O lines. The average spectral line energy distribution (SLED) of the entire sample is consistent with individual SLEDs and the IR pumping plus collisional excitation model, showing that the strongest lines are H$_2$O($2{02}-1_{11}$) and H$2$O($3{21}-3_{12}$).