No Arabic abstract
We provide the most accurate estimate yet of the bright end of the infrared (IR) luminosity functions (LFs) and the abundance of hyperluminous IR galaxies (HLIRGs) with IR luminosities > 10^13 L_solar, thanks to the combination of the high sensitivity, angular resolution, and large area of the LOFAR Deep Fields, which probes an unprecedented dynamic range of luminosity and volume. We cross-match Herschel sources and LOFAR sources in Bootes (8.63 deg^2), Lockman Hole (10.28 deg^2), and ELAIS-N1 (6.74 deg^2) with rms sensitivities of around 32, 22, and 20 mJy per beam, respectively. We divide the matched samples into unique and multiple categories. For the multiple matches, we de-blend the Herschel fluxes using the LOFAR positions and the 150-MHz flux densities as priors. We perform spectral energy distribution (SED) fitting, combined with multi-wavelength counterpart identifications and photometric redshift estimates, to derive IR luminosities. The depth of the LOFAR data allows us to identify highly complete (around 92% completeness) samples of bright Herschel sources with a simple selection based on the 250 micron flux density (45, 40, and 35 mJy in Bootes, Lockman Hole, and ELAIS-N1, respectively). Most of the bright Herschel sources fall into the unique category (i.e. a single LOFAR counterpart). For the multiple matches, there is excellent correspondence between the radio emission and the far-IR emission. We find a good agreement in the IR LFs with a previous study out to z around 6 which used de-blended Herschel data. Our sample gives the strongest and cleanest indication to date that the population of HLIRGs has surface densities of around 5 to 18 / deg^2 (with variations due to a combination of the applied flux limit and cosmic variance) and an uncertainty of a factor of 2. In comparison, the GALFORM semi-analytic model significantly under-predicts the abundance of HLIRGs.
We make use of multi-wavelength data of a large hyperluminous infrared (HLIRG) sample to derive their main physical properties, e.g., stellar mass, star-formation rate (SFR), volume density, contribution to the cosmic stellar mass density and to the cosmic SFR density. We also study the black hole (BH) growth rate and its relationship with the SFR of the host galaxy. We select 526 HLIRGs in three deep fields (Bo$o$tes, Lockman-Hole, ELAIS-N1) and adopt two spectral energy distribution (SED) fitting codes, CIGALE, which assumes energy balance, and CYGNUS, which is based on radiative transfer models and does not adopt energy balance principle. We use two different active galactic nucleus (AGN) models in CIGALE and three AGN models in CYGNUS to compare the results estimated using different SED fitting codes and different AGN models. The stellar mass, total IR luminosity and AGN luminosity agree well between different models with a typical median offset of 0.1 dex. The SFR estimates show the largest dispersions (up to 0.5 dex). This dispersion has an impact on the subsequent analysis, which may suggest that previous contradictory results could partly be due to different choices of methods. HLIRGs are ultra-massive galaxies with 99% of them having stellar masses larger than $10^{11} M_{odot}$. Our results reveal a higher space density of ultra-massive galaxies than found in previous surveys or predicted by simulations. We find that HLIRGs contribute more to the cosmic SFR density as redshift increases. In terms of BH growth, the two SED fitting methods provide different results. We can see a clear trend in which SFR decreases as AGN luminosity increases when using CYGNUS estimates, possibly implying quenching by AGN, while this trend is much weaker when using CIGALE estimates. This difference is also influenced by the dispersion between SlFR estimates obtained by the two codes.
By cross-correlating AKARI infrared (IR) sources with the SDSS galaxies, we identified 2357 infrared galaxies with a spectroscopic redshift. This is not just one of the largest samples of local IR galaxies, but AKARI provides crucial FIR bands (9, 18, 65, 90, 140, and 160um) in accurately measuring galaxy SED across the peak of the dust emission at ~100um. By fitting modern IR SED models to the AKARI photometry, we measured the total infrared luminosity (L_IR) of individual galaxies more accurately. Using this L_IR, we constructed luminosity functions of infrared galaxies at a median redshift of z=0.031, with 4 times larger sample than previous work. The LF agrees well with that at z=0.0082 (RBGS), showing smooth and continuous evolution toward higher redshift LFs measured in the AKARI NEP deep field. The derived local cosmic IR luminosity density is Omega_IR=3.8x10^8 LsunMpc^-3. We separate galaxies into AGN, star-forming, and composite by using the [NII]/Ha vs [OIII]/Hb line ratios. The fraction of AGN shows a continuous increase with increasing L_IR from 25% to 90% at 9<log L_IR<12.5. The SFR_Ha and L_[OIII] show good correlations with L_IR for SFG (star-forming galaxies) and AGN, respectively. The self-absorption corrected Ha/Hb ratio shows a weak increase with L_IR with a substantial scatter. When we separate IR LFs into contributions from AGN and SFG, the AGN contribution becomes dominant at L_IR>10^11Lsun, coinciding the break of the both SFG and AGN IR LFs. At L_IR<10^11Lsun, SFG dominates IR Lfs. Only 1.1% of Omega_IR is produced by LIRG, and only 0.03% is by ULIRG in the local Universe. This work also provides the most accurate infrared luminosity density of the local Universe to date. Compared with high redshift results from the AKARI NEP deep survey, we observed a strong evolution of Omega_IR^SFG ~(1+z)^4.1+-0.4 and Omega_IR^AGN ~(1+z)^4.1+-0.5 (abridged).
We present a catalogue of 179 hyperluminous infrared galaxies (HLIRGs) from the Imperial IRAS-FSS Redshift (IIFSCz) Catalogue. Of the 92 with detections in at least two far infrared bands, 62 are dominated by an M82-like starburst, 22 by an Arp220-like starburst and 8 by an AGN dust torus. On the basis of previous gravitational lensing studies and an examination of HST archive images for a further 5 objects, we estimate the fraction of HLIRGs that are significantly lensed to be 10-30%. We show simple infrared template fits to the SEDs of 23 HLIRGs with spectroscopic redshifts and at least 5 photometric bands. Most can be fitted with a combination of two simple templates: an AGN dust torus and an M82-like starburst. In the optical, 17 of the objects are fitted with QSO templates, 6 are fitted with galaxy templates. 20 of the 23 objects (87%) show evidence of an AGN either from the optical continuum or from the signature of an AGN dust torus, but the starburst component is the dominant contribution to bolometric luminosity in 14 out of 23 objects (61%). The implied star-formation rates, even after correcting for lensing magnification, are in excess of 1000 Mo /yr. We use infrared template-fitting models to predict fluxes for all HLIRGs at submillimetre wavelengths, and show predictions at 350 and 850 mu. Most would have 850 mu fluxes brighter than 5 mJy so should be easily detectable with current submillimetre telescopes. At least 15% should be detectable in the Planck all-sky survey at 350 mu and all Planck all-sky survey sources with z < 0.9 should be IIFSCz sources. From the luminosity-volume test we find that HLIRGs show strong evolution. A simple exponential luminosity evolution applied to all HLIRGs would be consistent with the luminosity functions found in redshift bins 0.3-0.5, 0.5-1 and 1-2.
We present infrared luminosity functions and dust mass functions for the EAGLE cosmological simulation, based on synthetic multi-wavelength observations generated with the SKIRT radiative transfer code. In the local Universe, we reproduce the observed infrared luminosity and dust mass functions very well. Some minor discrepancies are encountered, mainly in the high luminosity regime, where the EAGLE-SKIRT luminosity functions mildly but systematically underestimate the observed ones. The agreement between the EAGLE-SKIRT infrared luminosity functions and the observed ones gradually worsens with increasing lookback time. Fitting modified Schechter functions to the EAGLE-SKIRT luminosity and dust mass functions at different redshifts up to $z=1$, we find that the evolution is compatible with pure luminosity/mass evolution. The evolution is relatively mild: within this redshift range, we find an evolution of $L_{star,250}propto(1+z)^{1.68}$, $L_{star,text{TIR}}propto(1+z)^{2.51}$ and $M_{star,text{dust}}propto(1+z)^{0.83}$ for the characteristic luminosity/mass. For the luminosity/mass density we find $varepsilon_{250}propto(1+z)^{1.62}$, $varepsilon_{text{TIR}}propto(1+z)^{2.35}$ and $rho_{text{dust}}propto(1+z)^{0.80}$, respectively. The mild evolution of the dust mass density is in relatively good agreement with observations, but the slow evolution of the infrared luminosity underestimates the observed luminosity evolution significantly. We argue that these differences can be attributed to increasing limitations in the radiative transfer treatment due to increasingly poorer resolution, combined with a slower than observed evolution of the SFR density in the EAGLE simulation and the lack of AGN emission in our EAGLE-SKIRT post-processing recipe.
The primary goal of the Taiwan ECDFS Near-Infrared Survey (TENIS) is to find well screened galaxy candidates at z>7 (z dropout) in the Extended Chandra Deep Field-South (ECDFS). To this end, TENIS provides relatively deep J and Ks data (~25.3 ABmag, 5-sigma) for an area of 0.5*0.5 degree. Leveraged with existing data at mid-infrared to optical wavelengths, this allows us to screen for the most luminous high-z objects, which are rare and thus require a survey over a large field to be found. We introduce new color selection criteria to select a z>7 sample with minimal contaminations from low-z galaxies and Galactic cool stars; to reduce confusion in the relatively low angular resolution IRAC images, we introduce a novel deconvolution method to measure the IRAC fluxes of individual sources. Illustrating perhaps the effectiveness at which we screen out interlopers, we find only one z>7 candidate, TENIS-ZD1. The candidate has a weighted z_phot of 7.8, and its colors and luminosity indicate a young (45M years old) starburst galaxy with a stellar mass of 3.2*10^10 M_sun. The result matches with the observational luminosity function analysis and the semi-analytic simulation result based on the Millennium Simulations, which may over predict the volume density for high-z massive galaxies. The existence of TENIS-ZD1, if confirmed spectroscopically to be at z>7, therefore poses a challenge to current theoretical models for how so much mass can accumulate in a galaxy at such a high redshift.