Ultraluminous infrared galaxies: mergers of sub-L* galaxies?

A sample of 27 low-redshift, mostly cool, ultraluminous infrared galaxies (ULIRGs) has been imaged at 1.6 μm with the Hubble Space Telescope (HST) Near-Infrared Camera and Multi-Object Spectrometer (NICMOS). The majority (67%) of the sample's galaxies are multiple-nucleus galaxies with projected separations of up to 17 kpc, and the rest of the sample (33%) are single-nucleus galaxies, as determined by the NICMOS angular resolution limit. The average observed, integrated (host+nucleus) H magnitude of our HST H sample ULIRGs is -24.3, slightly above that of an L* galaxy (MH = -24.2), and 52% of the sample's galaxies have sub-L* luminosities. The ULIRGs in the HST H sample are not generated as a result of the merging of two luminous (i.e., ≥L*) spiral galaxies. Instead, the interactions and mergers occur in general between two, or in some cases more, less massive sub-L* (0.3-0.5L*) galaxies. Only one out of the 49 nuclei identified in the entire HST H sample has the properties of a bright quasar-like nucleus. On average, the brightest nuclei in the HST H sample galaxies (i.e., cool ULIRGs) are 1.2 mag fainter than warm ULIRGs and low-luminosity Bright Quasar Survey quasars (BQS QSOs) and 2.6 mag fainter than high-luminosity BQS QSOs. Since the progenitor galaxies involved in the merger are sub-L* galaxies, the mass of the central black hole in these ULIRGs would be only about (1-2) × 107 M☉, if the bulge-to-black hole mass ratio of nearby galaxies holds for ULIRGs. The estimated mass of the central black hole is similar to that of nearby Seyfert 2 galaxies but at least 1 order of magnitude lower than the massive black holes thought to be located at the center of high-luminosity QSOs. Massive nuclear starbursts with constant star formation rates of 10-40 M☉ yr-1 could contribute significantly to the nuclear H-band flux and are consistent with the observed nuclear H-band magnitudes of the ULIRGs in the HST H sample. An evolutionary merging scenario is proposed for the generation of the different types of ULIRGs and QSOs on the basis of the masses of the progenitors involved in the merging process. According to this scenario, cool ULIRGs would be the end product of the merging of two or more low-mass (0.3L*-0.5L*) disk galaxies. Warm ULIRGs and low-luminosity QSOs would be generated by a merger involving intermediate-mass (0.5 L*) disk galaxies. Under this scenario, warm ULIRGs could still be the dust-enshrouded phases of UV-bright low-luminosity QSOs, but cool ULIRGs, which are most ULIRGs, would not evolve into QSOs

Timeline analysis and wavelet multiscale analysis of the AKARI All-Sky Survey at 90 micron

We present a careful analysis of the point source detection limit of the AKARI All-Sky Survey in the WIDE-S 90 $\mu$m band near the North Ecliptic Pole (NEP). Timeline Analysis is used to detect IRAS sources and then a conversion factor is derived to transform the peak timeline signal to the interpolated 90 $\mu$m flux of a source. Combined with a robust noise measurement, the point source flux detection limit at S/N $>5$ for a single detector row is $1.1\pm0.1$ Jy which corresponds to a point source detection limit of the survey of $\sim$0.4 Jy. Wavelet transform offers a multiscale representation of the Time Series Data (TSD). We calculate the continuous wavelet transform of the TSD and then search for significant wavelet coefficients considered as potential source detections. To discriminate real sources from spurious or moving objects, only sources with confirmation are selected. In our multiscale analysis, IRAS sources selected above $4\sigma$ can be identified as the only real sources at the Point Source Scales. We also investigate the correlation between the non-IRAS sources detected in Timeline Analysis and cirrus emission using wavelet transform and contour plots of wavelet power spectrum. It is shown that the non-IRAS sources are most likely to be caused by excessive noise over a large range of spatial scales rather than real extended structures such as cirrus clouds.Comment: 16 pages, 19 figures, 5 tables, accepted for publication in MNRA

A Population of Dust-rich Quasars at z ~ 1.5

We report Herschel SPIRE (250, 350, and 500 μm) detections of 32 quasars with redshifts 0.5 ≤z < 3.6 from the Herschel Multi-tiered Extragalactic Survey (HerMES). These sources are from a MIPS 24 μm flux-limited sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field permit construction of the rest-frame spectral energy distributions (SEDs) from ultraviolet to the mid-infrared for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show dust temperatures in the range of 25-60 K, with a mean of 34 K. The FIR luminosities range from 10^(11.3) to 10^(13.5) L_☉, qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars observed at ~1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3-20 μm, rest frame), and the bolometric luminosities derived using the 5100 Å index may be underestimated for these FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities for quasars with higher near-infrared luminosities

Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo.Comment: 20 pages, 11 figures, 6 tables; to appear in the Astrophysical Journal (accepted November 29, 2020

Galaxy Clustering in Far-Infrared SWIRE Fields

We present measurements of galaxy clustering detected in the six SWIRE fields in all MIPS channels at 24, 70, and 160 microns. The measurements include the low-order 2-point angular correlation functions, and high-order probes including probability distribution functions and Rényi information

The extragalactic sub-mm population: predictions for the SCUBA Half-Degree Extragalactic Survey (SHADES)

We present predictions for the angular correlation function and redshift distribution for SHADES, the SCUBA HAlf-Degree Extragalactic Survey, which will yield a sample of around 300 sub-mm sources in the 850 micron waveband in two separate fields. Complete and unbiased photometric redshift information on these sub-mm sources will be derived by combining the SCUBA data with i) deep radio imaging already obtained with the VLA, ii) guaranteed-time Spitzer data at mid-infrared wavelengths, and iii) far-infrared maps to be produced by BLAST, the Balloon-borne Large-Aperture Sub-millimeter Telescope. Predictions for the redshift distribution and clustering properties of the final anticipated SHADES sample have been computed for a wide variety of models, each constrained to fit the observed number counts. Since we are dealing with around 150 sources per field, we use the sky-averaged angular correlation function to produce a more robust fit of a power-law shape w(theta)=(theta/A)^{-delta} to the model data. Comparing the predicted distributions of redshift and of the clustering amplitude A and slope delta, we find that models can be constrained from the combined SHADES data with the expected photometric redshift information.Comment: updated and improved version, accepted for publication in the MNRA

Overcoming Confusion Noise with Hyperspectral Imaging from PRIMAger

The PRobe far-Infrared Mission for Astrophysics (PRIMA) concept aims to perform mapping with spectral coverage and sensitivities inaccessible to previous FIR space telescopes. PRIMA's imaging instrument, PRIMAger, provides unique hyperspectral imaging simultaneously covering 25-235 $\mu$m. We synthesise images representing a deep, 1500 hr deg$^{-2}$ PRIMAger survey, with realistic instrumental and confusion noise. We demonstrate that we can construct catalogues of galaxies with a high purity ($>95$ per cent) at a source density of 42k deg$^{-2}$ using PRIMAger data alone. Using the XID+ deblending tool we show that we measure fluxes with an accuracy better than 20 per cent to flux levels of 0.16, 0.80, 9.7 and 15 mJy at 47.4, 79.7, 172, 235 $\mu$m respectively. These are a factor of $\sim$2 and $\sim$3 fainter than the classical confusion limits for 72-96 $\mu$m and 126-235 $\mu$m, respectively. At $1.5 \leq z \leq 2$, we detect and accurately measure fluxes in 8-10 of the 10 channels covering 47-235 $\mu$m for sources with $2 \leq$ log(SFR) $\leq 2.5$, a 0.5 dex improvement on what might be expected from the classical confusion limit. Recognising that PRIMager will operate in a context where high quality data will be available at other wavelengths, we investigate the benefits of introducing additional prior information. We show that by introducing even weak prior flux information when employing a higher source density catalogue (more than one source per beam) we can obtain accurate fluxes an order of magnitude below the classical confusion limit for 96-235 $\mu$m.Comment: 14 pages, 11 figure

The Far-Infrared Properties of Spatially Resolved AKARI Observations

We present the spatially resolved observations of IRAS sources from the Japanese infrared astronomy satellite AKARI All-Sky Survey during the performance verification (PV) phase of the mission. We extracted reliable point sources matched with IRAS point source catalogue. By comparing IRAS and AKARI fluxes, we found that the flux measurements of some IRAS sources could have been over or underestimated and affected by the local background rather than the global background. We also found possible candidates for new AKARI sources and confirmed that AKARI observations resolved IRAS sources into multiple sources. All-Sky Survey observations are expected to verify the accuracies of IRAS flux measurements and to find new extragalactic point sources.Comment: 11 pages, 7 figures, accepted publication in PASJ AKARI special issu

High-Redshift QSOs in the SWIRE Survey and the z~3 QSO Luminosity Function

We use a simple optical/infrared (IR) photometric selection of high-redshift QSOs that identifies a Lyman Break in the optical photometry and requires a red IR color to distinguish QSOs from common interlopers. The search yields 100 z~3 (U-dropout) QSO candidates with 19<r'<22 over 11.7 deg^2 in the ELAIS-N1 (EN1) and ELAIS-N2 (EN2) fields of the Spitzer Wide-area Infrared Extragalactic (SWIRE) Legacy Survey. The z~3 selection is reliable, with spectroscopic follow-up of 10 candidates confirming they are all QSOs at 2.83<z<3.44. We find that our z~4$ (g'-dropout) sample suffers from both unreliability and incompleteness but present 7 previously unidentified QSOs at 3.50<z<3.89. Detailed simulations show our z~3 completeness to be ~80-90% from 3.0<z<3.5, significantly better than the ~30-80% completeness of the SDSS at these redshifts. The resulting luminosity function extends two magnitudes fainter than SDSS and has a faint end slope of beta=-1.42 +- 0.15, consistent with values measured at lower redshift. Therefore, we see no evidence for evolution of the faint end slope of the QSO luminosity function. Including the SDSS QSO sample, we have now directly measured the space density of QSOs responsible for ~70% of the QSO UV luminosity density at z~3. We derive a maximum rate of HI photoionization from QSOs at z~3.2, Gamma = 4.8x10^-13 s^-1, about half of the total rate inferred through studies of the Ly-alpha forest. Therefore, star-forming galaxies and QSOs must contribute comparably to the photoionization of HI in the intergalactic medium at z~3.Comment: Accepted for publication in ApJ. emulateapj format. 23 pages, 17 figure