Mattia Vaccari mattia@mattiavaccari mattiavaccari

1. A&A Herschel Special Issue A&A Herschel Special Issue http://hermes.sussex.ac.uk Mattia Vaccari mattia@mattiavaccari.net www.mattiavaccari.net A) Summary We present the first measurement of the Submillimeter Local Luminosity Function based on HerMES data. Herschel SPIRE observations are combined with Spitzer, 2MASS & SDSS datasets in the SWIRE Lockman Hole (LH) and Spitzer Extragalactic First Look Survey (XFLS) fields (totaling 14.7 deg2) to optimally detect and identify Herschel sources. We compute the 250, 350, 500 µm as well as the IR bolometric (8-1000 µm) local (0<z<0.2) luminosity functionand thus derive local benchmarks for models of the formation and evolution of submillimeter galaxies at a very early stage of the Herschel mission. HerMES ~70 deg2 • B) HerMES SPIRE Source Extraction and Cross-Identification • (Roseboom et al., MNRAS, in preparation) • At the depth of our SDP imaging there are ~1-10 beams/src, significantly below the traditional confusion limit of 30 beams/src (see Figure 1). Thus we need to take account of source blending in performing source photometry and cross-identification • We achieve this by making use of the strong correspondence between 24 μm imaging from Spitzer and the SPIRE bands. The SPIRE fluxes of 24 μm sources are estimated using a linear inversion method which finds the best fit set of fluxes considering the 24 μm positions and the SPIRE maps (see mat). Similar methods have been used previously on Spitzer and BLAST data (Magnelli et al. 2009, Bethermin et al. 2010, Chapin et al. 2010) • Testing on simulations (see Figures 2/3 & 4) shows that our method is both reliable and returns high completeness for even the faintest SPIRE sources • While our method requires 24 μm detections, we estimate that in our deepest fields we are missing at most 15% of the faint SPIRE population with extreme S250/S24 flux ratios • For local (0<z<0.2) sources we estimate to be complete down to the adopted flux limits C) Optical/NIR Counterparts & SED Fitting The HerMES SPIRE Submillimeter Luminosity FunctionMattia Vaccari & Lucia Marchetti & Alberto Franceschini (University of Padova)Isaac Roseboom (University of Sussex) & the HerMES Consortium - http://hermes.sussex.ac.uk The accurate 24 μm positions used for Herschel source extraction allows us to reliably determine the SDSS and 2MASS counterparts of all our 0<z<0.2 sources. SDSS, NED and SWIRE follow-up spectroscopic redshifts are used along with SDSS photometric redshifts to obtain a complete redshift information. Archival multi-wavelength spectrophotometry (including IRAC and MIPS) then allows us to determine accurate photometric redshifts and IR luminosities for all 0<z<0.2 sources. D) SPIRE 250, 350, 500 µm & IR Bolometric (8-1000 µm) LLF • We evaluate the Monochromatic & IR Bolometric LLFs using the 1/Vmax estimator • We compare our estimates with models and measurements from recent literature • (Poisson) errors are estimated in each field and a weighted mean is then computed Figure 1. GOODS-N 250 μm image (left) & R-band optical image (right). Green squares are 24 μm sources. Circle is beam at 250 μm. It is clear that even if we consider only 24 μm sources there are still ~1 beam/src. The strong correspondence between 250 μm and 24 μm can be seen in the left image. • Testing was performed on two simulated datasets by taking the mock catalogs of Fernandez-Conde et al. (2008) and producing maps which match the observed properties (i.e. noise, PRF) of our GOODS-N (deep) and LH (shallow) data • Comparisons are performed for the shallow data between our XID method and two existing techniques; the ubiquitous p-statistic coupled with Sussextractor derived source catalogues, and a simpler map-based approach based on Bethermin et al. (2010) • At 250 μm the HerMES XID method is seen to clearly outperform the others • At the longer wavelengths the advantage currently offered by our method is limited Basic Linear Method If you can describe the flux of a point source as Empty Diamonds Sanders et al. 2003 Then the data (i.e. map) can be described as • Φ250,350,500 plots : black solid line is total model and color lines are different populations of Franceschini et al. (2010), dashed black line is model of Xu et al. (2001), dot-dot-dot-dashed black line is model of Negrello et al. (2007) • ΦIR plot : dashed, dot-dashed, dot-dot-dot-dashed and dotted black lines indicate predictions by Xu et al. (2001), Lagache et al. (2004), Negrello et al. (2007) and Valiante et al. (2009), dashed red line is modified Schechter (double exponential) best fit Figures 2/3. Simulations used for XID: GOODS-N 250 μm simulated map (left) & real map (right) HerMES XID Mathematical Framework Maximum likelihood solution N is covariance matrix A • E) Conclusions and future work • We provide useful local (0<z<0.2) benchmarks for submillimeter galaxy formation and evolution studies at a very early stage of the Herschel mission and thus pave the way for wider-area data soon to be provided by Herschel surveys such as HerMES and H-ATLAS • Using 24 μm positions and combining linear inversion and model selection techniques we reliably detect Herschel sources and identify them in multi-wavelength images • We find a slightly more abundant local submillimeter population than predicted by most models in recent literature and an IR bolometric (8-1000 μm) LLD of 1.3 108 L at z~0.1 • The continuation of the Herschel mission will yield larger samples and improved SED templates, providing better IR bolometric luminosity estimates and stronger constraints on models for galaxy evolution and dust emission from the local to the distant Universe