We provide a coherent, uniform measurement of the evolution of the logarithmic star formation rate (SFR) - stellar mass ($M_*$) relation, called the main sequence of star-forming galaxies (MS), for galaxies out to $z\sim5$. We measure the MS using mean stacks of 3 GHz radio continuum images to derive average SFRs for $\sim$200,000 mass-selected galaxies at $z>0.3$ in the COSMOS field. We describe the MS relation adopting a new model that incorporates a linear relation at low stellar mass (log($M_*$/M$_\odot$)$<$10) and a flattening at high stellar mass that becomes more prominent at low redshift ($z<1.5$). We find that the SFR density peaks at $1.5<z<2$ and at each epoch there is a characteristic stellar mass ($M_* = 1 - 4 \times 10^{10}\mathrm{M}_\odot$) that contributes the most to the overall SFR density. This characteristic mass increases with redshift, at least to $z\sim2.5$. We find no significant evidence for variations in the MS relation for galaxies in different environments traced by the galaxy number density at $0.3<z<3$, nor for galaxies in X-ray groups at $z\sim0.75$. We confirm that massive bulge-dominated galaxies have lower SFRs than disk-dominated galaxies at a fixed stellar mass at $z<1.2$. As a consequence, the increase in bulge-dominated galaxies in the local star-forming population leads to a flattening of the MS at high stellar masses. This indicates that "mass-quenching" is linked with changes in the morphological composition of galaxies at a fixed stellar mass.

We measure the 850-$\mu$m source densities of 46 candidate protoclusters selected from the Planck High-z catalogue (PHz) and the Planck Catalogue of Compact Sources (PCCS) that were followed up with Herschel-SPIRE and SCUBA-2. This paper aims to search for overdensities of 850-$\mu$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protoclusters, 25 have significant overdensities ($>$5 times the field counts), 11 have intermediate overdensities (3--5 times the field counts) and 10 have no overdensity ($<$3 times the field counts) of 850-$\mu$m sources. We find that the enhanced number densities are unlikely to be the result of sample variance. Compared with the number counts of another sample selected from Planck's compact source catalogues, this [PHz+PCCS]-selected sample has a higher fraction of candidate protoclusters with significant overdensities, though both samples show overdensities of 850-$\mu$m sources above intermediate level. Based on the estimated star-formation rate densities (SFRDs), we suggest that both samples can efficiently select protoclusters with starbursting galaxies near the redshift at which the global field SFRD peaks ($2 < z < 3$). Based on the confirmation of overdensities found here, future follow-up observations on other PHz targets may greatly increase the number of genuine DSFG-rich clusters/protoclusters.

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 VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey, in synergy with the SUDARE survey, is a deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). The observations for the CDFS field comprise about 4.38 deg$^2$ down to $r\sim26$ mag. The total on-sky time spans over four years in this field, distributed over four adjacent sub-fields. In this paper, we use the multi-epoch $r$-band imaging data to measure the variability of the detected objects and search for transients. We perform careful astrometric and photometric calibrations and point spread function (PSF) modeling. A new method, referring to as differential running-average photometry, is proposed to measure the light curves of the detected objects. With the method, the difference of PSFs between different epochs can be reduced, and the background fluctuations are also suppressed. Detailed uncertainty analysis and detrending corrections on the light curves are performed. We visually inspect the light curves to select variable objects, and present some objects with interesting light curves. Further investigation of these objects in combination with multi-band data will be presented in our forthcoming paper.

Variability has proven to be a powerful tool to detect active galactic nuclei (AGN) in multi-epoch surveys. The new-generation facilities expected to become operational in the next few years will mark a new era in time-domain astronomy and their wide-field multi-epoch campaigns will favor extensive variability studies. We present our analysis of AGN variability in the second half of the VST survey of the Wide Chandra Deep Field South (W-CDF-S), performed in the r band and covering a 2 sq. deg. area. The analysis complements a previous work, in which the first half of the area was investigated. We provide a reliable catalog of variable AGN candidates, which will be critical targets in future variability studies. We selected a sample of optically variable sources and made use of infrared data from the Spitzer mission to validate their nature by means of color-based diagnostics. We obtain a sample of 782 AGN candidates among which 12 are classified as supernovae, 54 as stars, and 232 as AGN. We estimate a contamination $\lesssim 20\%$ and a completeness $\sim 38\%$ with respect to mid-infrared selected samples.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

We use photometric redshifts derived from new $u$-band through 4.5$\mu$m Spitzer IRAC photometry in the 4.8\,deg$^2$ of the XMM-LSS field to construct surface density maps in the redshift range 0.1-1.5. Our density maps show evidence for large-scale structure in the form of filaments spanning several tens of Mpc. Using these maps, we identify 339 overdensities that our simulated lightcone analysis suggests are likely associated with dark matter haloes with masses, $M_{\rm halo}$, log($M_{\rm halo}/M_{\odot})>$13.7. From this list of overdensities we recover 43 of 70 known X-ray detected and spectroscopically confirmed clusters. The missing X-ray clusters are largely at lower redshifts and lower masses than our target log($M_{\rm halo}/M_{\odot})>$13.7. The bulk of the overdensities are compact, but a quarter show extended morphologies which include likely projection effects, clusters embedded in apparent filaments as well as at least one potential cluster merger (at $z\sim1.28$). The strongest overdensity in our highest redshift slice (at $z\sim1.5$) shows a compact red galaxy core potentially implying a massive evolved cluster.

We study the properties of star-forming galaxies selected at 610 MHz with the GMRT in a survey covering $\sim$1.86 deg$^2$ down to a noise of $\sim$7.1\,$\mu$Jy / beam. These were identified by combining multiple classification diagnostics: optical, X-ray, infrared and radio data. Of the 1685 SFGs from the GMRT sample, 496 have spectroscopic redshifts whereas 1189 have photometric redshifts. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4 GHz radio luminosity ratio $\rm{q_{IR}}$, decreases with increasing redshift: $\rm{q_{IR}\,=\,2.86\pm0.04(1\,+\,z)^{-0.20\pm0.02}}$ out to $z \sim 1.8$. We use the $\rm{V/V_{max}}$ statistic to quantify the evolution of the co-moving space density of the SFG sample. Averaged over luminosity our results indicate $\rm{\langle V/V_{max} \rangle}$ to be $\rm{0.51\,\pm\, 0.06}$, which is consistent with no evolution in overall space density. However we find $\rm V/V_{max}$ to be a function of radio luminosity, indicating strong luminosity evolution with redshift. We explore the evolution of the SFGs radio luminosity function by separating the source into five redshift bins and comparing to theoretical model predictions. We find a strong redshift trend that can be fitted with a pure luminosity evolution of the form $\rm{L_{610\,MHz}\,\propto\,(\,1+\,z)^{(2.95\pm0.19)-(0.50\pm0.15)z}}$. We calculate the cosmic SFR density since $\rm{z \sim 1.5}$ by integrating the parametric fits of the evolved 610\,MHz luminosity function. Our sample reproduces the expected steep decline in the star formation rate density since $\rm{z\,\sim\,1}$.

The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future.

We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170 million objects is constructed by a positional cross match of 51 public surveys. This high resolution optical, near-infrared, and mid-infrared catalogue is designed for photometric redshift estimation, extraction of fluxes in lower resolution far-infrared maps, and spectral energy distribution modelling. It collates, standardises, and provides value added derived quantities including corrected aperture magnitudes and astrometry correction over the Herschel extragalactic wide fields for the first time. $grizy$ fluxes are available on all fields with $g$ band data reaching $5\sigma$ point-source depths in a 2 arcsec aperture of 23.5, 24.4, and 24.6 (AB) mag at the 25th, 50th, and 75th percentiles, by area covered, across all HELP fields. It has $K$ or $K_s$ coverage over 1146 deg.$^2$ with depth percentiles of 20.2, 20.4, and 21.0 mag respectively. The IRAC Ch 1 band is available over 273 deg.$^2$ with depth percentiles of 17.7, 21.4, and 22.2 mag respectively. This paper defines the "masterlist" objects for the first data release (DR1) of HELP. This large sample of standardised total and corrected aperture fluxes, uniform quality flags, and completeness measures provides large well understood statistical samples over the full Herschel extragalactic sky.

This is the first of a series of papers based on sensitive 610 MHz observations of the ELAIS N1 field using the Giant Metrewave Radio Telescope. We describe the observations, processing and source catalogue extraction from a deep image with area of 1.86\,deg$^{2}$ and minimum noise of $\sim$7.1\,$\mu$Jy/beam. We compile a catalogue of 4290 sources with flux densities in the range 28.9\,$\mu$Jy - 0.503\,Jy, and derive the Euclidean-normalized differential source counts for sources with flux densities brighter than $\rm{35.5\,\mu\,Jy}$. Our counts show a flattening at 610 MHz flux densities below 1 mJy. Below the break the counts are higher than previous observations at this frequency, but generally consistent with recent models of the low-frequency source population. The radio catalogue is cross-matched against multi-wavelength data leading to identifications for 92\% and reliable redshifts for 72\% of our sample, with 19\% of the redshifts based on spectroscopy. For the sources with redshifts we use radio and X-ray luminosity, optical spectroscopy and mid-infrared colours to search for evidence of the presence of an Active Galactic Nucleus (AGN). We compare our identifications to predictions of the flux density distributions of star forming galaxies (SFGs) and AGN, and find a good agreement assuming the majority of the sources without redshifts are SFGs. We derive spectral index distributions for a sub-sample. The majority of the sources are steep spectra, with a median spectral index that steepens with frequency; $\mathrm{\alpha^{325}_{610}\,=\,-0.80\,\pm\,0.29}$, $\mathrm{\alpha^{610}_{1400}\,=\,-0.83\,\pm\,0.31}$ and $\mathrm{\alpha^{610}_{5000}\,=\,-1.12\,\pm\,0.15}$.

We present SCUBA-2 850-$\mu$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$\mu$m sources in 9/13 of these candidate protoclusters show significant overdensities compared to the field, with the probability $<$10$^{-2}$ assuming the sources are randomly distributed in the sky. Using the 250-, 350-, 500- and 850-$\mu$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution (SED) templates with an MCMC method. The photometric redshift distribution, peaking at $2<z<3$, is consistent with that of known $z>2$ protoclusters and the peak of the cosmic star-formation rate density (SFRD). We find that the 850-$\mu$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim$10$^{12}L_{\odot}$ and star-formation rates of SFR=(500-1,500)$M_{\odot}$yr$^{-1}$. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorised into four groups: six of them being high-redshift starbursting protoclusters, one being a lower-redshift cluster/protocluster, three being protoclusters that contain lensed DSFG(s) or are rich in 850-$\mu$m sources, and three regions without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850-$\mu$m sources.

We present an enhanced version of the multiwavelength spectral modeling code MAGPHYS that allows the estimation of galaxy photometric redshift and physical properties (e.g., stellar mass, star formation rate, dust attenuation) simultaneously, together with robust characterization of their uncertainties. The self-consistent modeling over ultraviolet to radio wavelengths in MAGPHYS+photo-z is unique compared to standard photometric redshift codes. The broader wavelength consideration is particularly useful for breaking certain degeneracies in color vs. redshift for dusty galaxies with limited observer-frame ultraviolet and optical data (or upper limits). We demonstrate the success of the code in estimating redshifts and physical properties for over 4,000 infrared-detected galaxies at 0.4<z<6.0 in the COSMOS field with robust spectroscopic redshifts. We achieve high photo-z precision ($\sigma_{\Delta z/(1+z_{spec})}\lesssim0.04$), high accuracy (i.e., minimal offset biases; median$(\Delta z/(1+z_{spec}))\lesssim0.02$), and low catastrophic failure rates ($\eta\simeq4\%$) over all redshifts. Interestingly, we find that a weak 2175A absorption feature in the attenuation curve models is required to remove a subtle systematic photo-z offset ($z_{phot}-z_{spec}\simeq-0.03$) that occurs when this feature is not included. As expected, the accuracy of derived physical properties in MAGPHYS+photo-z decreases strongly as redshift uncertainty increases. The all-in-one treatment of uncertainties afforded with this code is beneficial for accurately interpreting physical properties of galaxies in large photometric datasets. Finally, we emphasize that MAGPHYS+photo-z is not intended to replace existing photo-z codes, but rather offer flexibility to robustly interpret physical properties when spectroscopic redshifts are unavailable. The MAGPHYS+photo-z code is publicly available online.

A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

The analysis of the variability of active galactic nuclei (AGNs) at different wavelengths and the study of possible correlations among different spectral windows are nowadays a major field of inquiry. Optical variability has been largely used to identify AGNs in multivisit surveys. The strength of a selection based on optical variability lies in the chance to analyze data from surveys of large sky areas by ground-based telescopes. However the effectiveness of optical variability selection, with respect to other multiwavelength techniques, has been poorly studied down to the depth expected from next generation surveys. Here we present the results of our r-band analysis of a sample of 299 optically variable AGN candidates in the VST survey of the COSMOS field, counting 54 visits spread over three observing seasons spanning > 3 yr. This dataset is > 3 times larger in size than the one presented in our previous analysis (De Cicco et al. 2015), and the observing baseline is ~8 times longer. We push towards deeper magnitudes (r(AB) ~23.5 mag) compared to past studies; we make wide use of ancillary multiwavelength catalogs in order to confirm the nature of our AGN candidates, and constrain the accuracy of the method based on spectroscopic and photometric diagnostics. We also perform tests aimed at assessing the relevance of dense sampling in view of future wide-field surveys. We demonstrate that the method allows the selection of high-purity (> 86%) samples. We take advantage of the longer observing baseline to achieve great improvement in the completeness of our sample with respect to X-ray and spectroscopically confirmed samples of AGNs (59%, vs. ~15% in our previous work), as well as in the completeness of unobscured and obscured AGNs. The effectiveness of the method confirms the importance to develop future, more refined techniques for the automated analysis of larger datasets.

Future radio surveys will generate catalogues of tens of millions of radio sources, for which redshift estimates will be essential to achieve many of the science goals. However, spectroscopic data will be available for only a small fraction of these sources, and in most cases even the optical and infrared photometry will be of limited quality. Furthermore, radio sources tend to be at higher redshift than most optical sources and so a significant fraction of radio sources hosts differ from those for which most photometric redshift templates are designed. We therefore need to develop new techniques for estimating the redshifts of radio sources. As a starting point in this process, we evaluate a number of machine-learning techniques for estimating redshift, together with a conventional template-fitting technique. We pay special attention to how the performance is affected by the incompleteness of the training sample and by sparseness of the parameter space or by limited availability of ancillary multi-wavelength data. As expected, we find that the quality of the photometric-redshift degrades as the quality of the photometry decreases, but that even with the limited quality of photometry available for all sky-surveys, useful redshift information is available for the majority of sources, particularly at low redshift. We find that a template-fitting technique performs best with high-quality and almost complete multi-band photometry, especially if radio sources that are also X-ray emitting are treated separately. When we reduced the quality of photometry to match that available for the EMU all-sky radio survey, the quality of the template-fitting degraded and became comparable to some of the machine learning methods. Machine learning techniques currently perform better at low redshift than at high redshift, because of incompleteness of the currently available training data at high redshifts.

U-band observations with the LSST have yet to be fully optimized in cadence. The straw man survey design is a simple coverage of the medium-deep-fast survey. Here we argue that deep coverage of the four deep drilling fields (XMM-LSS, ECDFS, ELAIS-S1 and COSMOS) has a much higher scientific return, given that these are also the target of the Southern Hemisphere's Square Kilometer Array Pathfinder, the MeerKAT specifically, deep radio observations.

Quadruply lensed quasars are extremely rare objects, but incredibly powerful cosmological tools. Only few dozen are known in the whole sky. Here we present the spectroscopic confirmation of two new quadruplets WG0214-2105 and WG2100-4452 discovered by Agnello & Spiniello (2018) within the Dark Energy Survey (DES) public footprints. We have conducted spectroscopic follow-up of these systems with the Southern African Large Telescope as part of a program that aims at confirming the largest possible number of optically selected strong gravitational lensing systems in the Equatorial and Southern Hemisphere. For both systems, we present the spectra for the sources and deflectors that allowed us to estimate the source redshifts and unambiguously confirm their lensing nature. For the brighter deflector (WG2100-4452), we measure the stellar velocity dispersion from the spectrum. We also obtain photometry for both lenses, directly from DES multi-band images, isolating the lens galaxies from the quasar images. One of the quadruplets, WG0214-2105, was also observed by Pan-STARRS, allowing us to estimate the apparent brightness of each quasar image at two different epochs, and thus to find evidence for flux variability. This result could suggest a microlensing event for the faintest components, although intrinsic variability cannot be excluded with only two epochs. Finally, we present simple lens models for both quadruplets, obtaining Einstein radii, SIE velocity dispersions, ellipticities, and position angles of the lens systems, as well as time delay predictions assuming a concordance cosmological model.

This paper is part of a series discussing the results obtained in the framework of a wide international collaboration - the Lockman Hole Project - aimed at improving the extensive multiband coverage available in the Lockman Hole region, through novel deep, wide-area, multifrequency (60, 150, 350 MHz, and 1.4 GHz) radio surveys. This multifrequency, multi-band information will be exploited to get a comprehensive view of star formation and active galactic nucleus activities in the high-redshift Universe from a radio perspective. In this paper, we present novel 1.4 GHz mosaic observations obtained with the Westerbork Synthesis Radio Telescope. With an area coverage of 6.6 deg2, this is the largest survey reaching an rms noise of 11 uJy/beam. In this paper, we present the source catalogue (~6000 sources with flux densities S>55 uJy (5sigma), and we discuss the 1.4 GHz source counts derived from it. Our source counts provide very robust statistics in the flux range 0.1<S<1 mJy, and are in excellent agreement with other robust determinations obtained at lower and higher flux densities. A clear excess is found with respect to the counts predicted by the semi-empirical radio sky simulations developed in the framework of the Square Kilometre Array Simulated Skies project. A preliminary analysis of the identified (and classified) sources suggests this excess is to be ascribed to star-forming galaxies, which seem to show a steeper evolution than predicted.

We use two catalogues, a Herschel catalogue selected at 500 mu (HerMES) and an IRAS catalogue selected at 60 mu (RIFSCz), to contrast the sky at these two wavelengths. Both surveys demonstrate the existence of extreme starbursts, with star-formation rates (SFRs) > 5000 Msun/yr. The maximum intrinsic star-formation rate appears to be ~30,000 Msun/yr. The sources with apparent SFR estimates higher than this are in all cases either lensed systems, blazars, or erroneous photometric redshifts. At redshifts of 3 to 5, the time-scale for the Herschel galaxies to make their current mass of stars at their present rate of formation ~ 10^8 yrs, so these galaxies are making a significant fraction of their stars in the current star-formation episode. Using dust mass as a proxy for gas mass, the Herschel galaxies at redshift 3 to 5 have gas masses comparable to their mass in stars. Of the 38 extreme starbursts in our Herschel survey for which we have more complete SED information, over 50% show evidence for QSO-like optical emission, or exhibit AGN dust tori in the mid-infrared SEDs. In all cases however the infrared luminosity is dominated by a starburst component. We derive a mean covering factor for AGN dust as a function of redshift and derive black hole masses and black hole accretion rates. There is a universal ratio of black-hole mass to stellar mass, ~ 10^{-3}, driven by the strong period of star-formation and black-hole growth at z = 1-5.

The HELP project focuses on the data from ESA's Herschel mission, which covered over 1300$deg^2$ and is preparing to publish a multi-wavelength catalogue of millions of objects. Our main goal is to find the best approach to simultaneously fitting SEDs of millions of galaxies across a wide redshift range to obtain homogeneous estimates of the main physical parameters of detected IR galaxies. We perform SED fitting on the UV/NIR to FIR emission of 42 047 galaxies from the pilot HELP field: ELAIS N1. We use the latest release of CIGALE, a galaxy SED fitting code relying on energy balance, to deliver the main physical parameters such as $M_{star}$, SFR, and $L_{dust}$. We implement additional quality criteria to the fits by calculating $\chi^2$ values for the stellar and dust part of the spectra independently. These criteria allow us to identify the best fits and to identify peculiar galaxies. We perform the SED fitting by assuming three different dust attenuation laws separately allowing us to test the impact of the assumed law on estimated physical parameters. We implemented two additional quality value checks for the SED fitting method based on $M_{star}$ estimation and energy budget. This method allows us to identify possible objects with incorrect matching in the catalogue and peculiar galaxies; we found 351 possible candidates of lensed galaxies using two complementary $\chi^2$s criteria (stellar and IR) and $z_{phot}$ calculated for the IR part of the spectrum only. We find that the attenuation law has an important impact on the $M_{star}$ estimate (on average leading to disparities of a factor of two). We derive the relation between $M_{star}$ estimates obtained by different attenuation laws and we find the best recipe for our sample. We also make independent estimates of the total $L_{dust}$ parameter from stellar emission by fitting the galaxies with and without IR data separately.

Upcoming surveys with new radio observatories such as the Square Kilometer Array will generate a wealth of imaging data containing large numbers of radio galaxies. Different classes of radio galaxies can be used as tracers of the cosmic environment, including the dark matter density field, to address key cosmological questions. Classifying these galaxies based on morphology is thus an important step toward achieving the science goals of next generation radio surveys. Radio galaxies have been traditionally been classified as Fanaroff-Riley (FR) I and II, although some exhibit more complex 'bent' morphologies arising from environmental factors or intrinsic properties. In this work we present the FIRST Classifier, an on-line system for automated classification of Compact and Extended radio sources. We developed the FIRST Classifier based on a trained Deep Convolutional Neural Network Model to automate the morphological classification of compact and extended radio sources observed in the FIRST radio survey. Our model achieved an overall accuracy of 97% and a recall of 98%, 100%, 98% and 93% for Compact, BENT, FRI and FRII galaxies respectively. The current version of the FIRST classifier is able to predict the morphological class for a single source or for a list of sources as Compact or Extended (FRI, FRII and BENT).

We present images taken using the Gemini South Adaptive Optics Imager (GSAOI) with the Gemini Multiconjugate Adaptive Optics System (GeMS) in three 2 arcmin$^2$ fields in the Spitzer Extragalactic Representative Volume Survey. These GeMS/GSAOI observations are among the first $\approx 0.1^{''}$ resolution data in the near-infrared spanning extragalactic fields exceeding $1.5^{\prime}$ in size. We use these data to estimate galaxy sizes, obtaining results similar to those from studies with the Hubble Space Telescope, though we find a higher fraction of compact star forming galaxies at $z>2$. To disentangle the star-forming galaxies from active galactic nuclei (AGN), we use multiwavelength data from surveys in the optical and infrared, including far-infrared data from Herschel, as well as new radio continuum data from the Australia Telescope Compact Array and Very Large Array. We identify ultraluminous infrared galaxies (ULIRGs) at $z \sim 1-3$, which consist of a combination of pure starburst galaxies and Active Galactic Nuclei (AGN)/starburst composites. The ULIRGs show signs of recent merger activity, such as highly disturbed morphologies and include a rare candidate triple AGN. We find that AGN tend to reside in hosts with smaller scale sizes than purely star-forming galaxies of similar infrared luminosity. Our observations demonstrate the potential for MCAO to complement the deeper galaxy surveys to be made with the James Webb Space Telescope.

We have combined spectrosopic and photometric data from the Sloan Digital Sky Survey (SDSS) with $1.4$ GHz radio observations, conducted as part of the Stripe 82 $1-2$ GHz Snapshot Survey using the Karl G. Jansky Very Large Array (VLA), which covers $\sim100$ sq degrees, to a flux limit of 88 $\mu$Jy rms. Cross-matching the $11\,768$ radio source components with optical data via visual inspection results in a final sample of $4\,795$ cross-matched objects, of which $1\,996$ have spectroscopic redshifts and $2\,799$ objects have photometric redshifts. Three previously undiscovered Giant Radio Galaxies (GRGs) were found during the cross-matching process, which would have been missed using automated techniques. For the objects with spectroscopy we separate radio-loud Active Galactic Nuclei (AGN) and star-forming galaxies (SFGs) using three diagnostics and then further divide our radio-loud AGN into the HERG and LERG populations. A control matched sample of HERGs and LERGs, matched on stellar mass, redshift and radio luminosity, reveals that the host galaxies of LERGs are redder and more concentrated than HERGs. By combining with near-infrared data, we demonstrate that LERGs also follow a tight $K-z$ relationship. These results imply the LERG population are hosted by population of massive, passively evolving early-type galaxies. We go on to show that HERGs, LERGs, QSOs and star-forming galaxies in our sample all reside in different regions of a WISE colour-colour diagram. This cross-matched sample bridges the gap between previous `wide but shallow' and `deep but narrow' samples and will be useful for a number of future investigations.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is a Guaranteed Time program carried out with the ESO/VST telescope to provide deep optical imaging over two 4 deg$^2$ patches of the sky centred on the CDFS and ES1 pointings. We present the cosmic shear measurement over the 4 deg$^2$ covering the CDFS region in the $r$-band using LensFit. Each of the four tiles of 1 deg$^2$ has more than one hundred exposures, of which more than 50 exposures passed a series of image quality selection criteria for weak lensing study. The $5\sigma$ limiting magnitude in $r$- band is 26.1 for point sources, which is $\sim$1 mag deeper than other weak lensing survey in the literature (e.g. the Kilo Degree Survey, KiDS, at VST). The photometric redshifts are estimated using the VOICE $u,g,r,i$ together with near-infrared VIDEO data $Y,J,H,K_s$. The mean redshift of the shear catalogue is 0.87, considering the shear weight. The effective galaxy number density is 16.35 gal/arcmin$^2$, which is nearly twice the one of KiDS. The performance of LensFit on such a deep dataset was calibrated using VOICE-like mock image simulations. Furthermore, we have analyzed the reliability of the shear catalogue by calculating the star-galaxy cross-correlations, the tomographic shear correlations of two redshift bins and the contaminations of the blended galaxies. As a further sanity check, we have constrained cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. For a flat $\Lambda$CDM model we have obtained $\Sigma_8$ = $\sigma_8(\Omega_m/0.3)^{0.5}$ = $0.68^{+0.11}_{-0.15}$.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is proposed to obtain deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). At present, the observations for the CDFS field have been completed, and comprise in total about 4.9 deg$^2$ down to $r_\mathrm{AB}$$\sim$26 mag. In the companion paper by Fu et al. (2018), we present the weak lensing shear measurements for $r$-band images with seeing $\le$ 0.9 arcsec. In this paper, we perform image simulations to calibrate possible biases of the measured shear signals. Statistically, the properties of the simulated point spread function (PSF) and galaxies show good agreements with those of observations. The multiplicative bias is calibrated to reach an accuracy of $\sim$3.0%. We study the bias sensitivities to the undetected faint galaxies and to the neighboring galaxies. We find that undetected galaxies contribute to the multiplicative bias at the level of $\sim$0.3%. Further analysis shows that galaxies with lower signal-to-noise ratio (SNR) are impacted more significantly because the undetected galaxies skew the background noise distribution. For the neighboring galaxies, we find that although most have been rejected in the shape measurement procedure, about one third of them still remain in the final shear sample. They show a larger ellipticity dispersion and contribute to $\sim$0.2% of the multiplicative bias. Such a bias can be removed by further eliminating these neighboring galaxies. But the effective number density of the galaxies can be reduced considerably. Therefore efficient methods should be developed for future weak lensing deep surveys.

We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling $2.7$ Ms of flare-filtered exposure, with a $46$ ks median PN exposure time. We provide an X-ray catalog of 5242 sources detected in the soft (0.5-2 keV), hard (2-10 keV), and/or full (0.5-10 keV) bands with a 1% expected spurious fraction determined from simulations. A total of 2381 new X-ray sources are detected compared to previous source catalogs in the same area. Our survey has flux limits of $1.7\times10^{-15}$, $1.3\times10^{-14}$, and $6.5\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the soft, hard, and full bands, respectively, which is comparable to those of the XMM-COSMOS survey. We identify multiwavelength counterpart candidates for 99.9% of the X-ray sources, of which 93% are considered as reliable based on their matching likelihood ratios. The reliabilities of these high-likelihood-ratio counterparts are further confirmed to be $\approx 97$% reliable based on deep Chandra coverage over $\approx 5$% of the XMM-LSS region. Results of multiwavelength identifications are also included in the source catalog, along with basic optical-to-infrared photometry and spectroscopic redshifts from publicly available surveys. We compute photometric redshifts for X-ray sources in 4.5 deg$^2$ of our field where forced-aperture multi-band photometry is available; $>70$% of the X-ray sources in this subfield have either spectroscopic or high-quality photometric redshifts.

For a sample of star forming galaxies in the redshift interval 0.15$<$z$<$0.3, we study how both the relative strength of the AGN infra-red emission, compared to that due to the star formation (SF), and the numerical fraction of AGNs, change as a function of the total stellar mass of the hosting galaxy group (M$^{*}_{\mathrm{group}}$), between $10^{10.25}$ and $10^{11.9}$M$_{\odot}$. Using a multi-component SED fitting analysis, we separate the contribution of stars, AGN torus and star formation to the total emission at different wavelengths. This technique is applied to a new multi-wavelength data-set in the SIMES field (23 not redundant photometric bands), spanning the wavelength range from the UV (GALEX) to the far-IR (Herschel) and including crucial AKARI and WISE mid-IR observations (4.5 \mu m$<\lambda<$24 \mu m), where the BH thermal emission is stronger. This new photometric catalog, that includes our best photo-z estimates, is released through the NASA/IPAC Infrared Science Archive (IRSA). Groups are identified through a friends of friends algorithm ($\sim$62% purity, $\sim$51% completeness). We identified a total of 45 galaxies requiring an AGN emission component, 35 of which in groups and 10 in the field. We find BHAR$\propto ($M$^{*}_{\mathrm{group}})^{1.21\pm0.27}$ and (BHAR/SFR)$\propto ($M$^{*}_{\mathrm{group}})^{1.04\pm0.24}$ while, in the same range of M$^{*}_{\mathrm{group}}$, we do not observe any sensible change in the numerical fraction of AGNs. Our results indicate that the nuclear activity (i.e. the BHAR and the BHAR/SFR ratio) is enhanced when galaxies are located in more massive and richer groups.

With the advent of wide-area submillimeter surveys, a large number of high-redshift gravitationally lensed dusty star-forming galaxies (DSFGs) has been revealed. Due to the simplicity of the selection criteria for candidate lensed sources in such surveys, identified as those with $S_{500\mu m} > 100$ mJy, uncertainties associated with the modelling of the selection function are expunged. The combination of these attributes makes submillimeter surveys ideal for the study of strong lens statistics. We carried out a pilot study of the lensing statistics of submillimetre-selected sources by making observations with the Atacama Large Millimetre Array (ALMA) of a sample of strongly-lensed sources selected from surveys carried out with the Herschel Space Observatory. We attempted to reproduce the distribution of image separations for the lensed sources using a halo mass function taken from a numerical simulation which contains both dark matter and baryons. We used three different density distributions, one based on analytical fits to the halos formed in the EAGLE simulation and two density distributions (Singular Isothermal Sphere (SIS) and SISSA) that have been used before in lensing studies. We found that we could reproduce the observed distribution with all three density distributions, as long as we imposed an upper mass transition of $\sim$$10^{13} M_{\odot}$ for the SIS and SISSA models, above which we assumed that the density distribution could be represented by an NFW profile. We show that we would need a sample of $\sim$500 lensed sources to distinguish between the density distributions, which is practical given the predicted number of lensed sources in the Herschel surveys.

High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$).

The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic surveys of galaxies that could be performed by SPICA. The surveys are designed to address the question how star formation and black hole accretion activities evolved over cosmic time through spectral diagnostics of the physical conditions of the interstellar/circumnuclear media in galaxies. On the basis of results obtained with Herschel far-IR photometric surveys of distant galaxies and Spitzer and AKARI near- to mid-IR spectroscopic observations of nearby galaxies, we estimate the numbers of the galaxies at redshift z > 0.5, which are expected to be detected in the PAH features or dust continuum by a wide (10 deg^2) or deep (1 deg^2) blind survey, both for a given observation time of 600 hours. As by-products of the wide blind survey, we also expect to detect debris disks, through the mid-IR excess above the photospheric emission of nearby main-sequence stars, and we estimate their number. We demonstrate that the SPICA mid-IR surveys will efficiently provide us with unprecedentedly large spectral samples, which can be studied further in the far-IR with SPICA.

We describe the observing strategy, data reduction tools and early results of a supernova (SN) search project, named SUDARE, conducted with the ESO VST telescope aimed at measuring the rate of the different types of SNe in the redshift range 0.2<z<0.8. The search was performed in two of the best-studied extragalactic fields, CDFS and COSMOS, for which a wealth of ancillary data are available in the literature or public archives. (abridged) We obtained a final sample of 117 SNe, most of which are SNIa (57%) and the remaining core collapse events of which 44% type II, 22% type IIn and 34% type Ib/c. In order to link the transients, we built a catalog of ~1.3x10^5 galaxies in the redshift range 0<z<1 with a limiting magnitude K_AB=23.5 mag. We measured the SN rate per unit volume for SN Ia and core collapse SNe in different bin of redshifts. The values are consistent with other measurements from the literature. The dispersion of the rate measurements for SNe Ia is comparable with the scatter of the theoretical tracks for single (SD) and double degenerate (DD) binary systems models, therefore the data do not allow to disentangle among the two different progenitor scenarios. However, we may notice that among the three tested models, SD and two flavours of DD, either with a steep (DDC) or a wide (DDW) delay time distribution, the SD gives a better fit across the whole redshift range whereas the DDC better matches the steep rise up to redshift ~1.2. The DDW appears instead less favoured. The core collapse SN rate is fully consistent, unlike recent claims, with the prediction based on recent estimates of the star formation history, and standard progenitor mass range.

Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation.

Cross-identifying complex radio sources with optical or infra red (IR) counterparts in surveys such as the Australia Telescope Large Area Survey (ATLAS) has traditionally been performed manually. However, with new surveys from the Australian Square Kilometre Array Pathfinder (ASKAP) detecting many tens of million of radio sources such an approach is no longer feasible. This paper presents new software (LRPY - Likelihood Ratio in PYthon) to automate the process of cross-identifying radio sources with catalogues at other wavelengths. LRPY implements the Likelihood Ratio (LR) technique with a modification to account for two galaxies contributing to a sole measured radio component. We demonstrate LRPY by applying it to ATLAS DR3 and a {\it Spitzer}-based multi-wavelength fusion catalogue, identifying 3,848 matched sources via our LR-based selection criteria. A subset of 1987 sources have flux density values for all IRAC bands which allow us to use criteria to distinguish between active galactic nuclei (AGN) and star-forming galaxies (SFG). We find that 936 radio sources ($\approx47\,\%$) meet both of the Lacy and Stern AGN selection criteria. Of the matched sources, 295 have spectroscopic redshifts and we examine the radio to IR flux ratio vs redshift, proposing an AGN selection criterion below the Elvis radio-loud (RL) AGN limit for this dataset. Taking the union of all three AGN selection criteria we identify 956 as AGN ($\approx48\,\%$). From this dataset, we find a decreasing fraction of AGN with lower radio flux densities consistent with other results in the literature.

We discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

We model the $z \geq 6.6$ Ly$\alpha$ luminosity function to estimate the number of lensed high$-z$ Ly$\alpha$ emitters that may be detected by the Euclid Deep Survey. To span the whole range of possible predictions we exploit two Ly$\alpha$ luminosity function models and two strong gravitational lensing models from the literature. We show that the planned Euclid Deep Survey observing 40 deg$^2$ over the 920-1850 nm wavelength range down to a flux limit of $F_{lim}=5\times10^{-17}\,$erg s$^{-1}\,$cm$^{-2}$ will enable us to find between $\sim 0.85$ and $\sim 1.82$ deg$^{-2}$ lensed Ly$\alpha$ emitters at $z \geq 6.6$ depending on the adopted Ly$\alpha$ luminosity function and strong gravitational lensing model. The obvious [OII], [OIII] and H$\beta$ contaminants of the Ly$\alpha$ lensed population will be identified with the help of Euclid's spectral resolving power, while the SKA will enable the identification of the interloper population of H$\alpha$ emitters. By combining Euclid and the SKA, we will thus be able to identify, for the first time, a sample of $\sim 34$ to $\sim 73$ lensed Ly$\alpha$ emitters at $z \geq 6.6$.

We apply The Tractor image modeling code to improve upon existing multi-band photometry for the Spitzer Extragalactic Representative Volume Survey (SERVS). SERVS consists of post-cryogenic Spitzer observations at 3.6 and 4.5 micron over five well-studied deep fields spanning 18 square degrees. In concert with data from ground-based near-infrared (NIR) and optical surveys, SERVS aims to provide a census of the properties of massive galaxies out to z ~ 5. To accomplish this, we are using The Tractor to perform "forced photometry." This technique employs prior measurements of source positions and surface brightness profiles from a high-resolution fiducial band from the VISTA Deep Extragalactic Observations (VIDEO) survey to model and fit the fluxes at lower-resolution bands. We discuss our implementation of The Tractor over a square degree test region within the XMM-LSS field with deep imaging in 12 NIR/optical bands. Our new multi-band source catalogs offer a number of advantages over traditional position-matched catalogs, including 1) consistent source cross-identification between bands, 2) de-blending of sources that are clearly resolved in the fiducial band but blended in the lower-resolution SERVS data, 3) a higher source detection fraction in each band, 4) a larger number of candidate galaxies in the redshift range 5 < z < 6, and 5) a statistically significant improvement in the photometric redshift accuracy as evidenced by the significant decrease in the fraction of outliers compared to spectroscopic redshifts. Thus, forced photometry using The Tractor offers a means of improving the accuracy of multi-band extragalactic surveys designed for galaxy evolution studies. We will extend our application of this technique to the full SERVS footprint in the future.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

We present the VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey, a VST INAF Guaranteed Time program designed to provide optical coverage of two 4 deg$^2$ cosmic windows in the Southern hemisphere. VOICE provides the first, multi-band deep optical imaging of these sky regions, thus complementing and enhancing the rich legacy of longer-wavelength surveys with VISTA, Spitzer, Herschel and ATCA available in these areas and paving the way for upcoming observations with facilities such as the LSST, MeerKAT and the SKA. VOICE exploits VST's OmegaCAM optical imaging capabilities and completes the reduction of WFI data available within the ES1 fields as part of the ESO-Spitzer Imaging Extragalactic Survey (ESIS) program providing $ugri$ and $uBVR$ coverage of 4 and 4 deg$^2$ areas within the CDFS and ES1 field respectively. We present the survey's science rationale and observing strategy, the data reduction and multi-wavelength data fusion pipeline. Survey data products and their future updates will be released at http://www.mattiavaccari.net/voice/ and on CDS/VizieR.

The SCUBA-2 Ambitious Sky Survey (SASSy) is composed of shallow 850-$\umu$m imaging using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) on the James Clerk Maxwell Telescope. Here we describe the extraction of a catalogue of beam-sized sources from a roughly $120\,{\rm deg}^2$ region of the Galactic plane mapped uniformly (to an rms level of about 40\,mJy), covering longitude 120\degr\,$<$\,\textit{l}\,$<$\,140\degr\ and latitude $\abs{\textit{b}}$\,$<$\,2.9\degr. We used a matched-filtering approach to increase the signal-to-noise (S/N) ratio in these noisy maps and tested the efficiency of our extraction procedure through estimates of the false discovery rate, as well as by adding artificial sources to the real images. The primary catalogue contains a total of 189 sources at 850\,$\umu$m, down to a S/N threshold of approximately 4.6. Additionally, we list 136 sources detected down to ${\rm S/N}=4.3$, but recognise that as we go lower in S/N, the reliability of the catalogue rapidly diminishes. We perform follow-up observations of some of our lower significance sources through small targeted SCUBA-2 images, and list 265 sources detected in these maps down to ${\rm S/N}=5$. This illustrates the real power of SASSy: inspecting the shallow maps for regions of 850-$\umu$m emission and then using deeper targeted images to efficiently find fainter sources. We also perform a comparison of the SASSy sources with the Planck Catalogue of Compact Sources and the \textit{IRAS} Point Source Catalogue, to determine which sources discovered in this field might be new, and hence potentially cold regions at an early stage of star formation.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We present a multi-wavelength study into the nature of faint radio sources in a deep radio image with the Giant Meterwave Radio Telescope at 612\,MHz covering 1.2 deg$^2$ of the ELAIS N1 region. We detect 2800 sources above 50~$\mu$Jy\,beam$^{-1}$. By matching to multi-wavelength data we obtain a redshift estimate for 63\%, with 29\% based on spectroscopy. For 1526 of the sources with redshifts we use radio and X-ray luminosity, optical spectroscopy, mid-infrared colors, and 24$\mu$m and IR to radio flux ratios to search for the presence of an AGN. The analysis reveals a rapid change in the population as flux density decreases from $\sim$500\,$\mu$Jy to $\sim$100\,$\mu$Jy. We find that 80.3\% of the objects show no evidence of AGN and have multi-wavelength properties consistent with radio emission from star forming galaxies (SFG). We classify 11.4\% as Radio Quiet (RQ) AGN and the remaining 8.3\% as Radio Loud (RL) AGN. The redshift of all populations extends to $z > 3$ with a median of $\sim$1. The median radio and far-IR luminosity increases systematically from SFG, to RQ AGN and RL AGN. The median $q_{\rm 24 \mu m}$ for SFG, $0.89\pm0.01$ is slightly below that for RQ AGN, $1.05\pm0.03$, and both differ substantially from the value for RL AGN of $-0.06\pm0.07$. However, SFG and RQ AGN show no significant difference in far-IR/radio ratios and have statistically indistinguishable star formation rates inferred from radio and far-IR luminosities. We conclude that radio emission from host galaxies of RQ AGN in this flux density regime result primarily from star formation activity.

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

This is the second paper of a series in which we present measurements of the Supernova (SN) rates from the SUDARE survey. In this paper, we study the trend of the SN rates with the intrinsic colours, the star formation activity and the mass of the parent galaxies. We have considered a sample of about 130000 galaxies and a SN sample of about 50 events. We found that the SN Ia rate per unit mass is higher by a factor of six in the star-forming galaxies with respect to the passive galaxies. The SN Ia rate per unit mass is also higher in the less massive galaxies that are also younger. These results suggest a distribution of the delay times (DTD) less populated at long delay times than at short delays. The CC SN rate per unit mass is proportional to both the sSFR and the galaxy mass. The trends of the Type Ia and CC SN rates as a function of the sSFR and the galaxy mass that we observed from SUDARE data are in agreement with literature results at different redshifts. The expected number of SNe Ia is in agreement with the observed one for all four DTD models considered both in passive and star-forming galaxies so we can not discriminate between different progenitor scenarios. The expected number of CC SNe is higher than the observed one, suggesting a higher limit for the minimum progenitor mass. We also compare the expected and observed trends of the SN Ia rate with the intrinsic U - J colour of the parent galaxy, assumed as a tracer of the age distribution. While the slope of the relation between the SN Ia rate and the U - J color in star-forming galaxies can be reproduced well by all four DTD models considered, only the steepest of them is able to account for the rates and colour in star-forming and passive galaxies with the same value of the SN Ia production efficiency.

We present the Spitzer Data Fusion, a database incorporating far-ultraviolet to far-infrared flux measurements as well as photometric and spectroscopic redshifts for 4.4 million IRAC-selected sources detected over 8 extragalactic fields covering 65 deg$^2$ observed by Spitzer in all IRAC and MIPS bands during its cryogenic mission. Deeper Spitzer observations carried out during its warm mission over 5 sub-fields as part of the SERVS project are also presented and analysed in a similar fashion, detecting 2.8 million IRAC-selected sources over 18 deg$^2$ and merging them with multi-wavelength catalogues within the SERVS Data Fusion. When combined with Herschel SPIRE surveys and radio continuum observations over the same fields, the Spitzer Data Fusion and the SERVS Data Fusion provide an invaluable resource for multi-wavelength galaxy formation and evolution studies at infrared/millimetre/radio wavelengths. The catalogues and their future updates will be released at http://www.mattiavaccari.net/df/ and on CDS/VizieR.

We have developed a new prior-based source extraction tool, XID+, to carry out photometry in the Herschel SPIRE maps at the positions of known sources. XID+ is developed using a probabilistic Bayesian framework which provides a natural framework in which to include prior information, and uses the Bayesian inference tool, Stan, to obtain the full posterior probability distribution on flux estimates. In this paper, we discuss the details of XID+ and demonstrate the basic capabilities and performance by running it on simulated SPIRE maps resembling the COSMOS field, and comparing to the current prior-based source extraction tool DESPHOT. We show that not only does XID+ perform better on metrics such as flux accuracy and flux uncertainty accuracy, we illustrate how obtaining the posterior probability distribution can help overcome some of the issues inherent with maximum likelihood based source extraction routines. We run XID+ on the COSMOS SPIRE maps from HerMES, using a 24 $\mathrm{\mu m}$ catalogue as a prior and show the marginalised SPIRE colour-colour plot and marginalised contribution to the cosmic infrared background at the SPIRE wavelengths. XID+ is a core tool arising from the Herschel Extragalactic Legacy Project (HELP) and we discuss how additional work within HELP providing prior information on fluxes can and will be utilised. The software is available at https://github.com/H-E-L-P/XID_plus. We also provide the data product for COSMOS. We believe this is the first time that the full posterior probability of galaxy photometry has been provided as a data product.

The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust obscured star formation over 1300 deg$^2$. With this unprecedented combination of data sets, it is possible to investigate how the star formation vs stellar mass relation (main-sequence) of star-forming galaxies depends on environment. In this pilot study we explore this question between 0.1 < z < 3.2 using data in the COSMOS field. We estimate the local environment from a smoothed galaxy density field using the full photometric redshift probability distribution. We estimate star formation rates by stacking the SPIRE data from the Herschel Multi-tiered Extragalactic Survey (HerMES). Our analysis rules out the hypothesis that the main-sequence for star-forming systems is independent of environment at 1.5 < z < 2, while a simple model in which the mean specific star formation rate declines with increasing environmental density gives a better description. However, we cannot exclude a simple hypothesis in which the main-sequence for star-forming systems is independent of environment at z < 1.5 and z > 2. We also estimate the evolution of the star formation rate density in the COSMOS field and our results are consistent with previous measurements at z < 1.5 and z > 2 but we find a $1.4^{+0.3}_{-0.2}$ times higher peak value of the star formation rate density at $z \sim 1.9$.

We use 3035 Herschel-SPIRE 500$\mu$m sources from 20.3 sq deg of sky in the HerMES Lockman, ES1 and XMM-LSS areas to estimate the star-formation rate density at z = 1-6. 500 mu sources are associated first with 350 and 250 mu sources, and then with Spitzer 24 mu sources from the SWIRE photometric redshift catalogue. The infrared and submillimetre data are fitted with a set of radiative-transfer templates corresponding to cirrus (quiescent) and starburst galaxies. Lensing candidates are removed via a set of colour-colour and colour-redshift constraints. Star-formation rates are found to extend from < 1 to 20,000 Mo/yr. Such high values were also seen in the all-sky IRAS Faint Source Survey. Star-formation rate functions are derived in a series of redshift bins from 0-6, combined with earlier far-infrared estimates, where available, and fitted with a Saunders et al (1990) functional form. The star-formation-rate density as a function of redshift is derived and compared with other estimates. There is reasonable agreement with both infrared and ultraviolet estimates for z < 3, but we find higher star-formation-rate densities than ultraviolet estimates at z = 3-6. Given the considerable uncertainties in the submillimetre estimates, we can not rule out the possibility that the ultraviolet estimates are correct. But the possibility that the ultraviolet estimates have seriously underestimated the contribution of dust-shrouded star-formation can also not be excluded.

How did galaxies form and evolve? This is one of the most challenging questions in astronomy to- day. Answering it requires a careful combination of observational and theoretical work to reliably determine the observed properties of cosmic bodies over large portions of the distant Universe on the one hand, and accurately model the physical processes driving their evolution on the other. Most importantly, it requires bringing together disparate multi-wavelength and multi-resolution spectro-photometric datasets in an homogeneous and well-characterized manner so that they are suitable for a rigorous statistical analysis. The Herschel Extragalactic Legacy Project (HELP) funded by the EC FP7 SPACE program aims to achieve this goal by combining the expertise of optical, infrared and radio astronomers to provide a multi-wavelength database for the dis- tant Universe as an accessible value-added resource for the astronomical community. It will do so by bringing together multi-wavelength datasets covering the 1000 deg2 mapped by Herschel extragalactic surveys and thus creating a joint lasting legacy from several ambitious sky surveys.

We present the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $\chi^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$\sigma$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $\sigma_{\Delta z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $\eta=0.5$%. At $3<z<6$, using the unique database of spectroscopic redshifts in COSMOS, we find $\sigma_{\Delta z/(1+z_s)}$ = 0.021 and $\eta=13.2\% $. The deepest regions reach a 90\% completeness limit of 10$^{10}M_\odot$ to $z=4$. Detailed comparisons of the color distributions, number counts, and clustering show excellent agreement with the literature in the same mass ranges. COSMOS2015 represents a unique, publicly available, valuable resource with which to investigate the evolution of galaxies within their environment back to the earliest stages of the history of the Universe. The COSMOS2015 catalog is distributed via anonymous ftp (ftp://ftp.iap.fr/pub/from_users/hjmcc/COSMOS2015/) and through the usual astronomical archive systems (CDS, ESO Phase 3, IRSA).

We present a series of results from a clustering analysis of the first data release of the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey. VIDEO is the only survey currently capable of probing the bulk of stellar mass in galaxies at redshifts corresponding to the peak of star formation on degree scales. Galaxy clustering is measured with the two-point correlation function, which is calculated using a non parametric kernel based density estimator. We use our measurements to investigate the connection between the galaxies and the host dark matter halo using a halo occupation distribution methodology, deriving bias, satellite fractions, and typical host halo masses for stellar masses between $10^{9.35}M_{\odot}$ and $10^{10.85}M_{\odot}$, at redshifts $0.5<z<1.7$. Our results show typical halo mass increasing with stellar mass (with moderate scatter) and bias increasing with stellar mass and redshift consistent with previous studies. We find the satellite fraction increased towards low redshifts, increasing from $\sim 5\%$ at $z\sim 1.5$, to $\sim 20\%$ at $z\sim 0.6$, also increasing for lower mass galaxies. We combine our results to derive the stellar mass to halo mass ratio for both satellites and centrals over a range of halo masses and find the peak corresponding to the halo mass with maximum star formation efficiency to be $ \sim 2 \times10^{12} M_{\odot}$ over cosmic time, finding no evidence for evolution.

We present a list of candidate gravitationally lensed dusty star-forming galaxies (DSFGs) from the HerMES Large Mode Survey (HeLMS) and the Herschel Stripe 82 Survey (HerS). Together, these partially overlapping surveys cover 372 deg$^{2}$ on the sky. After removing local spiral galaxies and known radio-loud blazars, our candidate list of lensed DSFGs is composed of 77 sources with 500 $\mu$m flux densities ($S_{500}$) greater than 100 mJy. Such sources are dusty starburst galaxies similar to the first bright Sub Millimeter Galaxies (SMGs) discovered with SCUBA. We expect a large fraction of this list to be strongly lensed, with a small fraction made up of bright SMG-SMG mergers that appear as Hyper-Luminous Infrared Galaxies (HyLIRGs; $\rm L_{IR}>10^{13} L_{\odot}$). Thirteen of the 77 candidates have spectroscopic redshifts from CO spectroscopy with ground-based interferometers, putting them at $z>1$ and well above the redshift of the foreground lensing galaxies. The surface density of our sample is 0.21 $\pm$ 0.03 deg$^{-2}$. We present follow-up imaging of a few of the candidates confirming their lensing nature. The sample presented here is an ideal tool for higher resolution imaging and spectroscopic observations to understand detailed properties of starburst phenomena in distant galaxies.

We present the Spitzer-IRAC/MIPS Extragalactic survey (SIMES) in the South Ecliptic Pole (SEP) field. The large area covered (7.7 deg$^2$), together with one of the lowest Galactic cirrus emissions in the entire sky and a very extensive coverage by Spitzer, Herschel, Akari, and GALEX, make the SIMES field ideal for extragalactic studies. The elongated geometry of the SIMES area ($\approx$4:1), allowing for a significant cosmic variance reduction, further improves the quality of statistical studies in this field. Here we present the reduction and photometric measurements of the Spitzer/IRAC data. The survey reaches a depth of 1.93 and 1.75 $\mu$Jy (1$\sigma$) at 3.6 and 4.5 $\mu$m, respectively. We discuss the multiwavelength IRAC--based catalog, completed with optical, mid-- and far--IR observations. We detect 341,000 sources with F$_{3.6\mu m} \geq 3\sigma$. Of these, 10% have an associated 24 $\mu$m counterpart, while 2.7% have an associated SPIRE source. We release the catalog through the NASA/IPAC Infrared Science Archive (IRSA). Two scientific applications of these IRAC data are presented in this paper: first we compute integral number counts at 3.6 $\mu$m. Second, we use the [3.6]--[4.5] color index to identify galaxy clusters at z$>$1.3. We select 27 clusters in the full area, a result consistent with previous studies at similar depth.

We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

A complete, flux density limited sample of 96 faint ($> 0.5$ mJy) radio sources is selected from the 10C survey at 15.7 GHz in the Lockman Hole. We have matched this sample to a range of multi-wavelength catalogues, including SERVS, SWIRE, UKIDSS and optical data; multi-wavelength counterparts are found for 80 of the 96 sources and spectroscopic redshifts are available for 24 sources. Photometric reshifts are estimated for the sources with multi-wavelength data available; the median redshift of the sample is 0.91 with an interquartile range of 0.84. Radio-to-optical ratios show that at least 94 per cent of the sample are radio loud, indicating that the 10C sample is dominated by radio galaxies. This is in contrast to samples selected at lower frequencies, where radio-quiet AGN and starforming galaxies are present in significant numbers at these flux density levels. All six radio-quiet sources have rising radio spectra, suggesting that they are dominated by AGN emission. These results confirm the conclusions of Paper I that the faint, flat-spectrum sources which are found to dominate the 10C sample below $\sim 1$ mJy are the cores of radio galaxies. The properties of the 10C sample are compared to the SKADS Simulated Skies; a population of low-redshift starforming galaxies predicted by the simulation is not found in the observed sample.

How did galaxies form and evolve? This is one of the most challenging questions in astronomy today. Answering it requires a careful combination of observational and theoretical work to reliably determine the observed properties of cosmic bodies over large portions of the distant Universe on the one hand, and accurately model the physical processes driving their evolution on the other. Most importantly, it requires bringing together disparate multi-wavelength and multi-resolution spectro-photometric datasets in an homogeneous and well-characterized manner so that they are suitable for a rigorous statistical analysis. The Herschel Extragalactic Legacy Project (HELP) funded by the EC FP7 SPACE program aims to achieve this goal by combining the expertise of optical, infrared and radio astronomers to provide a multi-wavelength database for the distant Universe as an accessible value-added resource for the astronomical community. It will do so by bringing together multi-wavelength datasets covering the 1000 deg$^2$ mapped by Herschel extragalactic surveys in an homogeneous and well-characterized manner, creating a joint lasting legacy from several ambitious sky surveys.

We investigate the star-formation rate (SFR) and stellar mass ($M_*$) relation of a star-forming (SF) galaxy sample in the XMM-LSS field to $z\sim 3.0$ using the near-infrared data from the VISTA Deep Extragalactic Observations (VIDEO) survey. Combining VIDEO with broad-band photometry, we use the SED fitting algorithm CIGALE to derive SFRs and $M_*$ and have adapted it to account for the full photometric redshift PDF uncertainty. Applying a SF selection using the D4000 index, we find evidence for strong evolution in the normalisation of the SFR-$M_*$ relation out to $z\sim 3$ and a roughly constant slope of (SFR $\propto M_*^{\alpha}$) $\alpha=0.69\pm0.02$ to $z\sim 1.7$. We find this increases close to unity toward $z\sim2.65$. Alternatively, if we apply a colour selection, we find a distinct turnover in the SFR-$M_*$ relation between $0.7\lesssim z\lesssim2.0$ at the high mass end, and suggest that this is due to an increased contamination from passive galaxies. We find evolution of the specific SFR $\propto(1+z)^{2.60}$ at $\log(M_*)\sim$10.5, out to $z\lesssim2.4$ with an observed flattening beyond $z\sim$ 2 with increased stellar mass. Comparing to a range of simulations we find the analytical scaling relation approaches, that invoke an equilibrium model, a good fit to our data, suggesting that a continual smooth accretion regulated by continual outflows may be a key driver in the overall growth of SFGs.

We separate the extragalactic radio source population above ~50 uJy into active galactic nuclei (AGN) and star-forming sources. The primary method of our approach is to fit the infrared spectral energy distributions (SEDs), constructed using Spitzer/IRAC and MIPS and Herschel/SPIRE photometry, of 380 radio sources in the Extended Chandra Deep Field-South. From the fitted SEDs, we determine the relative AGN and star-forming contributions to their infrared emission. With the inclusion of other AGN diagnostics such as X-ray luminosity, Spitzer/IRAC colours, radio spectral index and the ratio of star-forming total infrared flux to k-corrected 1.4 GHz flux density, qIR, we determine whether the radio emission in these sources is powered by star formation or by an AGN. The majority of these radio sources (60 per cent) show the signature of an AGN at some wavelength. Of the sources with AGN signatures, 58 per cent are hybrid systems for which the radio emission is being powered by star formation. This implies that radio sources which have likely been selected on their star formation have a high AGN fraction. Below a 1.4 GHz flux density of 1 mJy, along with finding a strong contribution to the source counts from pure star-forming sources, we find that hybrid sources constitute 20-65 per cent of the sources. This result suggests that hybrid sources have a significant contribution, along with sources that do not host a detectable AGN, to the observed flattening of the source counts at ~1mJy for the extragalactic radio source population.

This work makes use of the VST observations to select variable sources. We use also the IR photometry, SED fitting and X-ray information where available to confirm the nature of the AGN candidates. The IR data, available over the full survey area, allow to confirm the consistency of the variability selection with the IR color selection method, while the detection of variability may prove useful to detect the presence of an AGN in IR selected starburst galaxies.

One of the most peculiar characteristics of Active Galactic Nuclei (AGN) is their variability over all wavelengths. This property has been used in the past to select AGN samples and is foreseen to be one of the detection techniques applied in future multi-epoch surveys, complementing photometric and spectroscopic methods. In this paper, we aim to construct and characterise an AGN sample using a multi-epoch dataset in the r band from the SUDARE-VOICE survey. Our work makes use of the VST monitoring program of an area surrounding the Chandra Deep Field South to select variable sources. We use data spanning a six month period over an area of 2 square degrees, to identify AGN based on their photometric variability. The selected sample includes 175 AGN candidates with magnitude r < 23 mag. We distinguish different classes of variable sources through their lightcurves, as well as X-ray, spectroscopic, SED, optical and IR information overlapping with our survey. We find that 12% of the sample (21/175) is represented by SN. Of the remaining sources, 4% (6/154) are stars, while 66% (102/154) are likely AGNs based on the available diagnostics. We estimate an upper limit to the contamination of the variability selected AGN sample of about 34%, but we point out that restricting the analysis to the sources with available multi-wavelength ancillary information, the purity of our sample is close to 80% (102 AGN out of 128 non-SN sources with multi-wavelength diagnostics). Our work thus confirms the efficiency of the variability selection method in agreement with our previous work on the COSMOS field; in addition we show that the variability approach is roughly consistent with the infrared selection.

We combine a deep 0.5~deg$^2$, 1.4~GHz deep radio survey in the Lockman Hole with infrared and optical data in the same field, including the SERVS and UKIDSS near-infrared surveys, to make the largest study to date of the host galaxies of radio sources with typical radio flux densities $\sim 50 \;\mu$Jy. 87% (1274/1467) of radio sources have identifications in SERVS to $AB\approx 23.1$ at 3.6 or 4.5$\mu$m, and 9% are blended with bright objects (mostly stars), leaving only 4% (59 objects) which are too faint to confidently identify in the near-infrared. We are able to estimate photometric redshifts for 68% of the radio sources. We use mid-infrared diagnostics to show that the source population consists of a mixture of star forming galaxies, rapidly accreting (cold mode) AGN and low accretion rate, hot mode AGN, with neither AGN nor starforming galaxies clearly dominating. We see the breakdown in the $K-z$ relation in faint radio source samples, and show that it is due to radio source populations becoming dominated by sources with radio luminosities $\sim 10^{23}\;{\rm WHz^{-1}}$. At these luminosities, both the star forming galaxies and the cold mode AGN have hosts with stellar luminosities about a factor of two lower than those of hot mode AGN, which continue to reside in only the most massive hosts. We show that out to at least $z\sim 2$, galaxies with stellar masses $>10^{11.5}\, M_{\odot}$ have a radio-loud fraction up to $\sim 30$%. This is consistent with there being a sufficient number of radio sources that radio-mode feedback could play a role in galaxy evolution.

We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs.

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

In this Chapter we present the motivation for undertaking both a wide and deep survey with the SKA in the context of studying AGN activity across cosmic time. With an rms down to 1 $\mu$Jy/beam at 1 GHz over 1,000 - 5,000 deg$^2$ in 1 year (wide tier band 1/2) and an rms down to 200 nJy/beam over 10 - 30 deg$^2$ in 2000 hours (deep tier band 1/2), these surveys will directly detect faint radio-loud and radio-quiet AGN (down to a 1 GHz radio luminosity of about $2\times10^{23}$ W/Hz at $z=6$). For the first time, this will enable us to conduct detailed studies of the cosmic evolution of radio AGN activity to the cosmic dawn ($z\gtrsim6$), covering all environmental densities.

It has become apparent that active galactic nuclei (AGN) may have a significant impact on the growth and evolution of their host galaxies and vice versa but a detailed understanding of the interplay between these processes remains elusive. Deep radio surveys provide a powerful, obscuration-independent tool for measuring both star formation and AGN activity in high-redshift galaxies. Multiwavelength studies of deep radio fields show a composite population of star-forming galaxies and AGN, with the former dominating at the lowest flux densities (S$_{1.4\mathrm{GHz}}<$100~$\mu$Jy). The sensitivity and resolution of the SKA will allow us to identify, and separately trace, the total star formation in the bulges of individual high-redshift galaxies, the related nuclear activity and any star formation occurring on larger scales within a disc. We will therefore gain a detailed picture of the apparently simultaneous development of stellar populations and black holes in the redshift range where both star-formation and AGN activity peak (1$\leq$z$\leq$4). In this chapter we discuss the role of the SKA in studying the connection between AGN activity and galaxy evolution, and the most critical technical requirements for such of studies

Radio wavelengths offer the unique possibility of tracing the total star-formation rate in galaxies, both obscured and unobscured. As such, they may provide the most robust measurement of the star-formation history of the Universe. In this chapter we highlight the constraints that the SKA can place on the evolution of the star-formation history of the Universe, the survey area required to overcome sample variance, the spatial resolution requirements, along with the multi-wavelength ancillary data that will play a major role in maximising the scientific promise of the SKA. The required combination of depth and resolution means that a survey to trace the star formation in the Universe should be carried out with a facility that has a resolution of at least ~0.5arcsec, with high sensitivity at < 1 GHz. We also suggest a strategy that will enable new parameter space to be explored as the SKA expands over the coming decade.

Optical variability has proven to be an effective way of detecting AGNs in imaging surveys, lasting from weeks to years. In the present work we test its use as a tool to identify AGNs in the VST multi-epoch survey of the COSMOS field, originally tailored to detect supernova events. We make use of the multi-wavelength data provided by other COSMOS surveys to discuss the reliability of the method and the nature of our AGN candidates. Our selection returns a sample of 83 AGN candidates; based on a number of diagnostics, we conclude that 67 of them are confirmed AGNs (81% purity), 12 are classified as supernovae, while the nature of the remaining 4 is unknown. For the subsample of AGNs with some spectroscopic classification, we find that Type 1 are prevalent (89%) compared to Type 2 AGNs (11%). Overall, our approach is able to retrieve on average 15% of all AGNs in the field identified by means of spectroscopic or X-ray classification, with a strong dependence on the source apparent magnitude. In particular, the completeness for Type 1 AGNs is 25%, while it drops to 6% for Type 2 AGNs. The rest of the X-ray selected AGN population presents on average a larger r.m.s. variability than the bulk of non variable sources, indicating that variability detection for at least some of these objects is prevented only by the photometric accuracy of the data. We show how a longer observing baseline would return a larger sample of AGN candidates. Our results allow us to assess the usefulness of this AGN selection technique in view of future wide-field surveys.

We present a multi-wavelength analysis of star-forming galaxies in the massive cluster MS0451.6-0305 at z $\sim$ 0.54 to shed new light on the evolution of the far-infrared-radio relationship in distant rich clusters. We have derived total infrared luminosities for a spectroscopically confirmed sample of cluster and field galaxies through an empirical relation based on $Spitzer$ MIPS 24 $\mu$m photometry. The radio flux densities were measured from deep Very Large Array 1.4 GHz radio continuum observations. We find the ratio of far-infrared to radio luminosity for galaxies in an intermediate redshift cluster to be $q_{\rm FIR}$ = 1.80$\pm$0.15 with a dispersion of 0.53. Due to the large intrinsic dispersion, we do not find any observable change in this value with either redshift or environment. However, a higher percentage of galaxies in this cluster show an excess in their radio fluxes when compared to low redshift clusters ($27^{+23}_{-13}\%$ to $11\%$), suggestive of a cluster enhancement of radio-excess sources at this earlier epoch. In addition, the far-infrared-radio relationship for blue galaxies, where $q_{\rm FIR}$ = 2.01$\pm$0.14 with a dispersion of 0.35, is consistent with the predicted value from the field relationship, although these results are based on a sample from a single cluster.

We present Herschel far-IR photometry and spectroscopy as well as ground based CO observations of an intermediate redshift (0.21 < z < 0.88) sample of Herschel-selected (ultra)-luminous infrared galaxies (L_IR > 10^11.5L_sun). With these measurements we trace the dust continuum, far-IR atomic line emission, in particular [CII]\,157.7microns, as well as the molecular gas of z~0.3 (U)LIRGs and perform a detailed investigation of the interstellar medium of the population. We find that the majority of Herschel-selected intermediate redshift (U)LIRGs have L_CII/L_FIR ratios that are a factor of about 10 higher than that of local ULIRGs and comparable to that of local normal and high-$z$ star forming galaxies. Using our sample to bridge local and high-z [CII] observations, we find that the majority of galaxies at all redshifts and all luminosities follow a L_CII-L_FIR relation with a slope of unity, from which local ULIRGs and high-z AGN dominated sources are clear outliers. We also confirm that the strong anti-correlation between the L_CII/L_FIR ratio and the far-IR color L_60/L_100 observed in the local Universe holds over a broad range of redshifts and luminosities, in the sense that warmer sources exhibit lower L_CII/L_FIR at any epoch. Intermediate redshift ULIRGs are also characterised by large molecular gas reservoirs and by lower star formation efficiencies compared to that of local ULIRGs. The high L_CII/L_FIR ratios, the moderate star formation efficiencies (L_LIR/L_CO or L_IR/M_gas) and the relatively low dust temperatures of our sample (which are also common characteristics of high-z star forming galaxies with ULIRG-like luminosities) indicate that the evolution of the physical properties of (U)LIRGs between the present day and z > 1 is already significant by z ~ 0.3.

We have studied in detail a sample of 967 SPIRE sources with 5-sigma detections at 350 and 500 micron and associations with Spitzer-SWIRE 24 micron galaxies in the HerMES-Lockman survey area, fitting their mid- and far-infrared, and submillimetre, SEDs in an automatic search with a set of six infrared templates. For almost 300 galaxies we have modelled their SEDs individually to ensure the physicality of the fits. We confirm the need for the new cool and cold cirrus templates, and also of the young starburst template, introduced in earlier work. We also identify 109 lensing candidates via their anomalous SEDs and provide a set of colour-redshift constraints which allow lensing candidates to be identified from combined Herschel and Spitzer data. The picture that emerges of the submillimetre galaxy population is complex, comprising ultraluminous and hyperluminous starbursts, lower luminosity galaxies dominated by interstellar dust emission, lensed galaxies and galaxies with surprisingly cold (10-13K) dust. 11 % of 500$\mu$m selected sources are lensing candidates. 70 % of the unlensed sources are ultraluminous infrared galaxies and 26 % are hyperluminous. 34 % are dominated by optically thin interstellar dust ('cirrus') emission, but most of these are due to cooler dust than is characteristic of our Galaxy. At the highest infrared luminosities we see SEDs dominated by M82, Arp220 and young starburst types, in roughly equal proportions.

We use 10,387 galaxies from the Herschel Astrophysical TeraHertz Large Area Survey (H-ATLAS) to probe the far-infrared radio correlation (FIRC) of star forming galaxies as a function of redshift, wavelength, and effective dust temperature. All of the sources in our 250 {\mu}m-selected sample have spectroscopic redshifts, as well as 1.4 GHz flux density estimates measured from the Faint Images of the Radio Sky at Twenty centimetres (FIRST) survey. This enables us to study not only individual sources, but also the average properties of the 250 {\mu}m selected population using median stacking techniques. We find that individual sources detected at $\geq 5\sigma$ in both the H-ATLAS and FIRST data have logarithmic flux ratios (i.e. FIRC $q_\lambda$ parameters) consistent with previous studies of the FIRC. In contrast, the stacked values show larger $q_\lambda$, suggesting excess far-IR flux density/luminosity in 250{\mu}m selected sources above what has been seen in previous analyses. In addition, we find evidence that 250 {\mu}m sources with warm dust SEDs have a larger 1.4 GHz luminosity than the cooler sources in our sample. Though we find no evidence for redshift evolution of the monochromatic FIRC, our analysis reveals significant temperature dependence. Whilst the FIRC is reasonably constant with temperature at 100 {\mu}m, we find increasing inverse correlation with temperature as we probe longer PACS and SPIRE wavelengths. These results may have important implications for the use of monochromatic dust luminosity as a star formation rate indicator in star-forming galaxies, and in the future, for using radio data to determine galaxy star formation rates.

We present the most complete study to date of the X-ray emission from star-formation in high redshift (median z=0.7; z<1.5), IR-luminous (L_IR=10^10-10^13 L_sun) galaxies detected by Herschel's PACS and SPIRE instruments. For our purpose we take advantage of the deepest X-ray data to date, the Chandra deep fields (North and South). Sources which host AGN are removed from our analysis by means of multiple AGN indicators. We find an AGN fraction of 18+/-2 per cent amongst our sample and note that AGN entirely dominate at values of log[L_X/L_IR]>-3 in both hard and soft X-ray bands. From the sources which are star-formation dominated, only a small fraction are individually X-ray detected and for the bulk of the sample we calculate average X-ray luminosities through stacking. We find an average soft X-ray to infrared ratio of log[L_SX/L_IR]=-4.3 and an average hard X-ray to infrared ratio of log[L_HX/L_IR]=-3.8. We report that the X-ray/IR correlation is approximately linear through the entire range of L_IR and z probed and, although broadly consistent with the local (z<0.1) one, it does display some discrepancies. We suggest that these discrepancies are unlikely to be physical, i.e. due to an intrinsic change in the X-ray properties of star-forming galaxies with cosmic time, as there is no significant evidence for evolution of the L_X/L_IR ratio with redshift. Instead they are possibly due to selection effects and remaining AGN contamination. We also examine whether dust obscuration in the galaxy plays a role in attenuating X-rays from star-formation, by investigating changes in the L_X/L_IR ratio as a function of the average dust temperature. We conclude that X-rays do not suffer any measurable attenuation in the host galaxy.

We present Herschel/PACS 100 and 160 micron integrated photometry for the 323 galaxies in the Herschel Reference Survey (HRS), a K-band-, volume-limited sample of galaxies in the local Universe. Once combined with the Herschel/SPIRE observations already available, these data make the HRS the largest representative sample of nearby galaxies with homogeneous coverage across the 100-500 micron wavelength range. In this paper, we take advantage of this unique dataset to investigate the properties and shape of the far-infrared/sub-millimeter spectral energy distribution in nearby galaxies. We show that, in the stellar mass range covered by the HRS (8<log(M*/Msun)<12), the far-infrared/sub-millimeter colours are inconsistent with a single modified black-body having the same dust emissivity index beta for all galaxies. In particular, either beta decreases, or multiple temperature components are needed, when moving from metal-rich/gas-poor to metal-poor/gas-rich galaxies. We thus investigate how the dust temperature and mass obtained from a single modified black-body depend on the assumptions made on beta. We show that, while the correlations between dust temperature, galaxy structure and star formation rate are strongly model dependent, the dust mass scaling relations are much more reliable, and variations of beta only change the strength of the observed trends.

We gather infrared (IR) photometric data from 8 to 500 microns (Spitzer, WISE, IRAS and Herschel) for all of the HRS galaxies. Draine & Li (2007) models are fit to the data from which the stellar contribution has been carefully removed. We find that our photometric coverage is sufficient to constrain all of the models parameters and that a strong constraint on the 20-60 microns range is mandatory to estimate the relative contribution of the photo-dissociation regions to the IR SED. The SED models tend to systematically under-estimate the observed 500 microns flux densities, especially for low mass systems. We provide the output parameters for all of the galaxies: the minimum intensity of the interstellar radiation field (ISRF), the fraction of PAH, the relative contribution of PDR and evolved stellar population to the dust heating, the $M_{dust}$ and the $L_{IR}$. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the integrated properties of galaxies, such as $M_*$, SFR, metallicity, H$\alpha$ and H-band surface brightness, and the FUV attenuation. A good correlation between the fraction of PAH and the metallicity is found implying a weakening of the PAH emission in galaxies with low metallicities. The intensity of the IRSF and the H-band and H$\alpha$ surface brightnesses are correlated, suggesting that the diffuse dust component is heated by both the young stars in star forming regions and the diffuse evolved population. We use these results to provide a new set of IR templates calibrated with Herschel observations on nearby galaxies and a mean SED template to provide the z=0 reference for cosmological studies. For the same purpose, we put our sample on the SFR-$M_*$ diagram. The templates are compared to the most popular IR SED libraries, enlightening a large discrepancy between all of them in the 20-100 microns range.

We study a sample of Herschel-PACS selected galaxies within the GOODS-South and the COSMOS fields in the framework of the PACS Evolutionary Probe (PEP) project. Starting from the rich multi-wavelength photometric data-sets available in both fields, we perform a broad-band Spectral Energy Distribution (SED) decomposition to disentangle the possible active galactic nucleus (AGN) contribution from that related to the host galaxy. We find that 37 per cent of the Herschel-selected sample shows signatures of nuclear activity at the 99 per cent confidence level. The probability to reveal AGN activity increases for bright ($L_{\rm 1-1000} > 10^{11} \rm L_{\odot}$) star-forming galaxies at $z>0.3$, becoming about 80 per cent for the brightest ($L_{\rm 1-1000} > 10^{12} \rm L_{\odot}$) infrared (IR) galaxies at $z \geq 1$. Finally, we reconstruct the AGN bolometric luminosity function and the super-massive black hole growth rate across cosmic time up to $z \sim 3$ from a Far-Infrared (FIR) perspective. This work shows general agreement with most of the panchromatic estimates from the literature, with the global black hole growth peaking at $z \sim 2$ and reproducing the observed local black hole mass density with consistent values of the radiative efficiency $\epsilon_{\rm rad}$ ($\sim$0.07).

We report the first results from a spectroscopic survey of the [CII] 158um line from a sample of intermediate redshift (0.2<z<0.8) (ultra)-luminous infrared galaxies, (U)LIRGs (LIR>10^11.5 Lsun), using the SPIRE-Fourier Transform Spectrometer (FTS) on board the Herschel Space Observatory. This is the first survey of [CII] emission, an important tracer of star-formation, at a redshift range where the star-formation rate density of the Universe increases rapidly. We detect strong [CII] 158um line emission from over 80% of the sample. We find that the [CII] line is luminous, in the range (0.8-4)x10^(-3) of the far-infrared continuum luminosity of our sources, and appears to arise from photodissociation regions on the surface of molecular clouds. The L[CII]/LIR ratio in our intermediate redshift (U)LIRGs is on average ~10 times larger than that of local ULIRGs. Furthermore, we find that the L[CII]/LIR and L[CII]/LCO(1-0) ratios in our sample are similar to those of local normal galaxies and high-z star-forming galaxies. ULIRGs at z~0.5 show many similarities to the properties of local normal and high-z star forming galaxies. Our findings strongly suggest that rapid evolution in the properties of the star forming regions of luminous infrared galaxies is likely to have occurred in the last 5 billion years.

The Herschel Multi-tiered Extragalactic Survey (HerMES) is the largest Guaranteed Time Key Programme on the Herschel Space Observatory. With a wedding cake survey strategy, it consists of nested fields with varying depth and area totalling ~380 deg^2. In this paper, we present deep point source catalogues extracted from Herschel-SPIRE observations of all HerMES fields, except for the later addition of the 270 deg^2 HeLMS field. These catalogues constitute the second Data Release (DR2) made in October 2013. A subset of these catalogues, which consists of bright sources extracted from Herschel-SPIRE observations completed by May 1, 2010 (covering ~ 74 deg^2) were released earlier in the first extensive Data Release (DR1) in March 2012. Two different methods are used to generate the point source catalogues, the SUSSEXtractor (SXT) point source extractor used in two earlier data releases (EDR and EDR2) and a new source detection and photometry method. The latter combines an iterative source detection algorithm, StarFinder (SF), and a De-blended SPIRE Photometry (DESPHOT) algorithm. We use end-to-end Herschel-SPIRE simulations with realistic number counts and clustering properties to characterise basic properties of the point source catalogues, such as the completeness, reliability, photometric and positional accuracy. Over 500, 000 catalogue entries in HerMES fields (except HeLMS) are released to the public through the HeDAM website (http://hedam.oamp.fr/herMES).

The potential for Planck to detect clusters of dusty, star-forming galaxies at z greater than 1 is tested by examining the Herschel-SPIRE images of Planck Early Release Compact Source Catalog (ERCSC) sources lying in fields observed by the HerMES survey. Of the 16 Planck sources that lie in the roughly 90 sq. deg. examined, we find that twelve are associated with single bright Herschel sources. The remaining four are associated with overdensities of Herschel sources, making them candidate clusters of dusty, starforming galaxies. We use complementary optical and NIR data for these clumps to test this idea, and find evidence for the presence of galaxy clusters in all four cases. We use photometric redshifts and red sequence galaxies to estimate the redshifts of these clusters, finding that they range from 0.8 to 2.3. These redshifts imply that the Herschel sources in these clusters, which contribute to the detected Planck flux, are forming stars very rapidly, with typical total cluster star formation rates greater than 1000Msun per yr. The high redshift clusters discovered in these observations are used to constrain the epoch of cluster galaxy formation, finding that the galaxies in our clusters are 1 to 1.5 Gy old at z about 1 to 2. Prospects for the discovery of further clusters of dusty galaxies are discussed, using not only all sky Planck surveys, but also deeper, smaller area, Herschel surveys.

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

We quantify the fraction of the cosmic infrared background (CIB) that originates from galaxies identified in the UV/optical/near-infrared by stacking 81,250 (~35.7 arcmin^2) K-selected sources (K_AB < 24.0), split according to their rest-frame U - V vs. V - J colors into 72,216 star-forming and 9,034 quiescent galaxies, on maps from Spitzer/MIPS (24um), Herschel/SPIRE (250, 350, 500um), Herschel/PACS (100, 160um), and AzTEC (1100um). The fraction of the CIB resolved by our catalog is (69 $\pm$ 15)% at 24um, (78 $\pm$ 17)% at 70um, (58 $\pm$ 13)% at 100um, (78 $\pm$ 18)% at 160um, (80 $\pm$ 17)% at 250um, (69 $\pm$ 14)% at 350um, (65 $\pm$ 12)% at 500um, and (45 $\pm$ 8)% at 1100um. Of that total, about 95% originates from star-forming galaxies, while the remaining 5% is from apparently quiescent galaxies. The CIB at $\lambda$ < 200um is sourced predominantly from galaxies at z < 1, while at $\lambda$ > 200um the bulk originates from 1 < z < 2. Galaxies with stellar masses log(M/ M_sun)=9.5-11 are responsible for the majority of the CIB, with those in the log(M/ M_sun)=9.5-10 contributing mostly at $\lambda$ < 250um, and those in the log(M/ M_sun)=10.5-11 bin dominating at $\lambda$ > 350um. The contribution from galaxies in the log(M/ M_sun)=9.0-9.5 and log(M/ M_sun)=11.0-12.0 stellar mass bins contribute the least, both of order 5%, although the highest stellar-mass bin is a significant contributor to the luminosity density at z > 2. The luminosities of the galaxies responsible for the CIB shifts from a combination of "normal" and luminous infrared galaxies (LIRGs) at $\lambda$ < 160um, to LIRGs at 160um < $\lambda$ < 500um, to finally LIRGs and ultra-luminous infrared galaxies (ULIRGs) at $\lambda$ > 500um. Stacking analyses were performed with SIMSTACK (available at http://www.astro.caltech.edu/~viero/viero_homepage/toolbox.html) which accounts for possible biases due to clustering.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

We present new photometric data from our Herschel Key Programme, the Dwarf Galaxy Survey (DGS), dedicated to the observation of the gas and dust in 48 low-metallicity environments. They were observed with PACS and SPIRE onboard Herschel at 70,100,160,250,350, and 500 microns. We focus on a systematic comparison of the derived FIR properties (FIR luminosity, dust mass, dust temperature and emissivity index) with more metal-rich galaxies and investigate the detection of a potential submm excess. The data reduction method is adapted for each galaxy to derive the most reliable photometry from the final maps. PACS flux densities are compared with the MIPS 70 and 160 microns bands. We use colour-colour diagrams and modified blackbody fitting procedures to determine the dust properties of the DGS galaxies. We also include galaxies from the Herschel KINGFISH sample, containing more metal-rich environments, totalling 109 galaxies. The location of the DGS galaxies on Herschel colour-colour diagrams highlights the differences in global environments of low-metallicity galaxies. The dust in DGS galaxies is generally warmer than in KINGFISH galaxies (T_DGS~32 K, T_KINGFISH~23 K). The emissivity index, beta, is ~1.7 in the DGS, but metallicity does not make a strong effect on beta. The dust-to-stellar mass ratio is lower in low-metallicity galaxies: M_dust/M_star~0.02% for the DGS vs 0.1% for KINGFISH. Per unit dust mass, dwarf galaxies emit ~6 times more in the FIR than higher metallicity galaxies. Out of the 22 DGS galaxies detected at 500 micron, 41% present an excess in the submm not explained by our dust SED model. The excess mainly appears in lower metallicity galaxies (12+log(O/H) < 8.3), and the strongest excesses are detected in the most metal-poor galaxies. We stress the need for observations longwards of the Herschel wavelengths to detect any submm excess appearing beyond 500 micron.

We propose a co-ordinated multi-observatory survey at the North Ecliptic Pole. This field is the natural extragalactic deep field location for most space observatories (e.g. containing the deepest Planck, WISE and eROSITA data), is in the continuous viewing zones for e.g. Herschel, HST, JWST, and is a natural high-visibility field for the L2 halo orbit of SPICA with deep and wide-field legacy surveys already planned. The field is also a likely deep survey location for the forthcoming Euclid mission. It is already a multi-wavelength legacy field in its own right (e.g. AKARI, LOFAR, SCUBA-2): the outstanding and unparalleled continuous mid-IR photometric coverage in this field and nowhere else enables a wide range of galaxy evolution diagnostics unachievable in any other survey field, by spanning the wavelengths of redshifted PAH and silicate features and the peak energy output of AGN hot dust. We argue from the science needs of Euclid and JWST, and from the comparative multiwavelength depths, that the logical approach is (1) a deep (H-UDF) UV/optical tile in the NEP over ~10 square arcminutes, and (2) an overlapping wide-field UV/optical HST survey tier covering >100 square arcminutes, with co-ordinated submm SPIRE mapping up to or beyond the submm point source confusion limit over a wider area and PACS data over the shallower HST tier.

Dust-obscured galaxies (DOGs) are a UV-faint, IR-bright galaxy population that reside at z~2 and are believed to be in a phase of dusty star-forming and AGN activity. We present far-IR observations of a complete sample of DOGs in the 2 deg2 of COSMOS. The 3077 DOGs have <z>=1.9+/-0.3 and are selected from 24um and r+ observations using a color cut of r+ - [24]>=7.5 (AB mag) and S24>=100uJy. Based on the near-IR SEDs, 47% are star-formation dominated and 10% are AGN-dominated. We use SPIRE far-IR photometry from HerMES to calculate the IR luminosity and characteristic dust temperature for the 1572 (51%) DOGs that are detected at 250um (>=3{\sigma}). For the remaining 1505 (49%) that are undetected, we perform a median stacking analysis to probe fainter luminosities. Detected and undetected DOGs have average IR luminosities of (2.8+/-0.4) x 1012 LSun and (0.77+-0.08)x10^12LSun, and dust temperatures of 34+/-7 K and 37+/-3 K, respectively. The IR luminosity function of DOGs with S24>=100uJy is calculated, using far-IR observations and stacking. DOGs contribute 10-30% to the total star formation rate density of the Universe at z=1.5-2.5, dominated by 250um detected and bump DOGs. For comparison, DOGs contribute 30% to the star-formation rate density of the Universe for all galaxies at z=1.5-2.5 with S24>=100uJy. DOGs have a large scatter about the star-formation main sequence and their specific star-formation rates show that the observed phase of star-formation could be responsible for their observed stellar mass at z~2.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

In this work we explore the impact of the presence of an active galactic nucleus (AGN) on the mid- and far-infrared (IR) properties of galaxies as well as the effects of simultaneous AGN and starburst activity in these same galaxies. To do this we apply a multi-component, multi-band spectral synthesis technique to a sample of 250 micron selected galaxies of the Herschel Multi-tiered Extragalactic Survey, with IRS spectra available for all galaxies. Our results confirm that the inclusion of the IRS spectra plays a crucial role in the spectral analysis of galaxies with an AGN component improving the selection of the best-fit hot dust model (torus). We find a correlation between the obscured star formation rate (SFR) derived from the IR luminosity of the starburst component, SFR_IR and SFR_PAH, derived from the luminosity of the PAH features, L_PAH, with SFR_FIR taking higher values than SFR_PAH. The correlation is different for AGN- and starburst-dominated objects. The ratio of L_PAH to that of the starburst component, L_PAH/L_SB, is almost constant for AGN-dominated objects but decreases with increasing L_SB for starburst-dominated objects. SFR_FIR increases with the accretion luminosity, L_acc, with the increase less prominent for the very brightest, unobscured AGN-dominated sources. We find no correlation between the masses of the hot and cold dust components. We interpret this as a non-constant fraction of gas driven by the gravitational effects to the AGN while the starburst is ongoing. We also find no evidence of the AGN affecting the temperature of the cold dust component, though this conclusion is mostly based on objects with a non-dominant AGN component. We conclude that our findings do not provide evidence that the presence of AGN affects the star formation process in the host galaxy, but rather that the two phenomena occur simultaneously over a wide range of luminosities.

The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.

Using new homogeneous LFs in the FUV and in the FIR Herschel/PEP and Herschel/HerMES, we study the evolution of the dust attenuation with redshift. With this information in hand, we are able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density SFRD_TOT. By integrating SFRD_TOT, we follow the mass building and analyze the redshift evolution of the stellar mass density (SMD). This letter aims at providing a complete view of star formation from the local universe to z = 4 and, using assumptions on earlier star formation history, compares this evolution to what was known before in an attempt to draw a homogeneous picture of the global evolution of star formation in galaxies. The main conclusions of this letter are: 1) the dust attenuation A_FUV is found to increase from z = 0 to z \sim 1.2 and then starts to decrease up to our last data point at z = 3.6; 2) the estimated SFRD confirms published results up to z = 2. At z > 2, we observe either a plateau or a small increase up to z = 3 and then a likely decrease up to z = 3.6; 3) the peak of A_FUV is delayed with respect to the plateau of SFRD_TOT and a likely origin might be found in the evolution of the bright ends of the FUV and FIR LFs; 4) using assumptions (namely exponential rise and linear rise with time) for the evolution of the star formation density from z = 3.6 to z_form = 10, we integrate SFRD_TOT and find a good agreement with the published SMDs.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $\mu$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $\Delta$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $\Delta$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

We have revised the SWIRE Photometric Redshift Catalogue to take account of new optical photometry in several of the SWIRE areas, and incorporating 2MASS and UKIDSS near infrared data. Aperture matching is an important issue for combining near infrared and optical data, and we have explored a number of methods of doing this. The increased number of photometric bands available for the redshift solution results in improvements both in the rms error and, especially, in the outlier rate. We have also found that incorporating the dust torus emission into the QSO templates improves the performance for QSO redshift estimation. Our revised redshift catalogue contains over 1 million extragalactic objects, of which 26288 are QSOs.

The infrared (IR) emission of M_* galaxies (10^{10.4} < M_{star} < 10^{11.0} M_\sun) in galaxy pairs, derived using data obtained in Herschel (PEP/HerMES) and Spitzer (S-COSMOS) surveys, is compared to that of single disk galaxies in well matched control samples to study the cosmic evolution of the star-formation enhancement induced by galaxy-galaxy interaction. Both the mean IR SED and mean IR luminosity of star-forming galaxies (SFGs) in SFG+SFG (S+S) pairs in the redshift bin of 0.6 < z < 1 are consistent with no star-formation enhancement. SFGs in S+S pairs in a lower redshift bin of 0.2 < z < 0.6 show marginal evidence for a weak star-formation enhancement. Together with the significant and strong sSFR enhancement shown by SFGs in a local sample of S+S pairs (obtained using previously published Spitzer observations), our results reveal a trend for the star-formation enhancement in S+S pairs to decrease with increasing redshift. Between z=0 and z=1, this decline of interaction-induced star-formation enhancement occurs in parallel with the dramatic increase (by a factor of ~10) of the sSFR of single SFGs, both can be explained by the higher gas fraction in higher z disks. SFGs in mixed pairs (S+E pairs) do not show any significant star-formation enhancement at any redshift. The difference between SFGs in S+S pairs and in S+E pairs suggests a modulation of the sSFR by the inter-galactic medium IGM in the dark matter halos (DMH) hosting these pairs.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

We present new observational determination of the evolution of the rest-frame 70 and 160 micron and total infrared (TIR) galaxy luminosity functions (LFs) using 70 micron data from the Spitzer Wide-area Infrared Extragalactic Legacy Survey (SWIRE). The LFs were constructed for sources with spectroscopic redshifts only in the XMM-LSS and Lockman Hole fields from the SWIRE photometric redshift catalogue. The 70 micron and TIR LFs were constructed in the redshift range 0<z<1.2 and the 160 micron LF was constructed in the redshift range 0<z<0.5 using a parametric Bayesian and the vmax methods. We assume in our models, that the faint-end power-law index of the LF does not evolve with redshifts. We find the the double power-law model is a better representation of the IR LF than the more commonly used power-law and Gaussian model. We model the evolution of the FIR LFs as a function of redshift where where the characteristic luminosity, $L^\ast$ evolve as $\propto(1+z)^{\alpha_\textsc{l}}$. The rest-frame 70 micron LF shows a strong luminosity evolution out to z=1.2 with alpha_l=3.41^{+0.18}_{-0.25}. The rest-frame 160 micron LF also showed rapid luminosity evolution with alpha_l=5.53^{+0.28}_{-0.23} out to z=0.5. The rate of evolution in luminosity is consistent with values estimated from previous studies using data from IRAS, ISO and Spitzer. The TIR LF evolves in luminosity with alpha_l=3.82^{+0.28}_{-0.16} which is in agreement with previous results from Spitzer 24 micron which find strong luminosity evolution. By integrating the LF we calculated the co-moving IR luminosity density out to z=1.2, which confirm the rapid evolution in number density of LIRGs and ULIRGs which contribute ~68^{+10}_{-07} % to the co-moving star formation rate density at z=1.2. Our results based on 70 micron data confirms that the bulk of the star formation at z=1 takes place in dust obscured objects.

We discuss the clustering properties of galaxies with signs of ongoing star formation detected by the Spitzer Space Telescope at 24mum band in the SWIRE Lockman Hole field. The sample of mid-IR-selected galaxies includes ~20,000 objects detected above a flux threshold of S24mum=310muJy. We adopt optical/near-IR color selection criteria to split the sample into the lower-redshift and higher-redshift galaxy populations. We measure the angular correlation function on scales of theta=0.01-3.5 deg, from which, using the Limber inversion along with the redshift distribution established for similarly selected source populations in the GOODS fields (Rodighiero et al. 2010), we obtain comoving correlation lengths of r0=4.98+-0.28 h^-1 Mpc and r0 =8.04+-0.69 h^-1 Mpc for the low-z (<z>=0.7) and high-z (<z>=1.7) subsamples, respectively. Comparing these measurements with the correlation functions of dark matter halos identified in the Bolshoi cosmological simulation (Klypin et al. 2011}, we find that the high-redshift objects reside in progressively more massive halos reaching Mtot>3e12 h^-1 Msun, compared to Mtot>7e11 h^-1 Msun for the low-redshift population. Approximate estimates of the IR luminosities based on the catalogs of 24mum sources in the GOODS fields show that our high-z subsample represents a population of "distant ULIRGs" with LIR>10^12Lsun, while the low-z subsample mainly consists of "LIRGs", LIR~10^11Lsun. The comparison of number density of the 24mum selected galaxies and of dark matter halos with derived minimum mass Mtot shows that only 20% of such halos may host star-forming galaxies.

The wide-area imaging surveys with the {\it Herschel} Space Observatory at sub-mm wavelengths have now resulted in catalogs of order one hundred thousand dusty, star-burst galaxies. We make a statistical estimate of $N(z)$ using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 $\mu$m from the Herschel Multi-tiered Extragalactic Survey (HerMES) centered on the Bo\"{o}tes field. We cross-correlate {\it Herschel} galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey (SDSS), the NOAO Deep Wide Field Survey (NDWFS) and the Spitzer Deep Wide Field Survey (SDWFS). We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs of the Cosmological Evolution Survey (COSMOS), as the COSMOS field has superior spectroscopy coverage. We model-fit the clustering auto and cross-correlations of {\it Herschel} and optical/IR galaxy samples to estimate $N(z)$ and clustering bias factors. The $S_{350} > 20$ mJy galaxies have a bias factor varying with redshift as $b(z)=1.0^{+1.0}_{-0.5}(1+z)^{1.2^{+0.3}_{-0.7}}$. This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 10$^{12}$ to 10$^{13}$ M$_{\sun}$. We find that the redshift distribution peaks around $z \sim 0.5$ to 1 for galaxies selected at 250 $\mu$m with an average redshift of $< z > = 1.8 \pm 0.2$. For 350 and 500 $\mu$m-selected SPIRE samples the peak shifts to higher redshift, but the average redshift remains the same with a value of $1.9 \pm 0.2$.

The Herschel Reference Survey (HRS) is a guaranteed time Herschel key project aimed at studying the physical properties of the interstellar medium in galaxies of the nearby universe. This volume limited, K-band selected sample is composed of galaxies spanning the whole range of morphological types (from ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster). We present flux density measurements of the whole sample of 323 galaxies of the HRS in the three bands of the Spectral and Photometric Imaging Receiver (SPIRE), at 250, 350 and 500 microns. Aperture photometry is performed on extended galaxies and point spread function (PSF) fitting on timeline data for unresolved objects; we carefully estimate errors and upper limits. The flux densities are found to be in good agreement with those of the HeViCS and KINGFISH key projects in all SPIRE bands, and of the Planck consortium at 350 and 550 microns, for the galaxies in common. This submillimetre catalogue of nearby galaxies is a benchmark for the study of the dust properties in the local universe, giving the zero redshift reference for any cosmological survey.

We use Herschel-PACS far-infrared data, combined with previous multi-band information and mid-IR spectra, to properly account for the presence of an active nucleus and constrain its energetic contribution in luminous infrared (IR) sources at z\sim2. The sample is composed of 24 sources in the GOODS-South field, with typical IR luminosity of 10^{12} Lo. Data from the 4 Ms Chandra X-ray imaging in this field are also used to identify and characterize AGN emission. We reproduce the observed spectral energy distribution (SED), decomposed into a host-galaxy and an AGN component. A smooth-torus model for circum-nuclear dust is used to account for the direct and re-processed contribution from the AGN. We confirm that galaxies with typical L_{8-1000um}\sim10^{12}Lo at z\sim2 are powered predominantly by star-formation. An AGN component is present in nine objects (\sim35% of the sample) at the 3sigma confidence level, but its contribution to the 8-1000 um emission accounts for only \sim5% of the energy budget. The AGN contribution rises to \sim23% over the 5-30 um range (in agreement with Spitzer IRS results) and to \sim60% over the narrow 2-6 um range. The presence of an AGN is confirmed by X-ray data for 3 (out of nine) sources, with X-ray spectral analysis indicating the presence of significant absorption, i.e. NH\sim10^{23}-10^{24} cm^{-2}. An additional source shows indications of obscured AGN emission from X-ray data. The comparison between the mid-IR--derived X-ray luminosities and those obtained from X-ray data suggests that obscuration is likely present also in the remaining six sources that harbour an AGN according to the SED-fitting analysis.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We use the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) to explore the specific star-formation activity of galaxies and their evolution near the peak of the cosmic far-infrared (FIR) background at 70 and 160um. We use a stacking analysis to determine the mean FIR properties of well defined subsets of galaxies at flux levels well below the FIR catalogue detection limits of SWIRE and other Spitzer surveys. We tabulate the contribution of different subsets of galaxies to the FIR background at 70um and 160um. These long wavelengths provide a good constraint on the bolometric, obscured emission. The large area provides good constraints at low z and in finer redshift bins than previous work. At all redshifts we find that the specific FIR Luminosity (sLFIR) decreases with increasing mass, following a trend L_FIR/M* propto M_* ^beta with beta =-0.38\pm0.14. This is a more continuous change than expected from the {Delucia2007} semi-analytic model suggesting modifications to the feedback prescriptions. We see an increase in the sLFIR by about a factor of ~100 from 0<z<2 and find that the sLFIR evolves as (1+z)^alpha with alpha=4.4\pm0.3 for galaxies with 10.5 < log M*/Msun < 12. This is considerably steeper than the {Delucia2007} semi-analytic model (alpha \sim 2.5). When separating galaxies into early and late types on the basis of the optical/IR spectral energy distributions we find that the decrease in sLFIR with stellar mass is stronger in early type galaxies (beta ~ -0.46), while late type galaxies exhibit a flatter trend (beta \sim -0.15). The evolution is strong for both classes but stronger for the early type galaxies. The early types show a trend of decreasing strength of evolution as we move from lower to higher masses while the evolution of the late type galaxies has little dependence on stellar mass. We suggest that in late-type galaxies we are seeing a consistently declining sSFR..

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We present broad band photometry and photometric redshifts for 187611 sources located in ~0.5deg^2 in the Lockman Hole area. The catalog includes 389 X-ray detected sources identified with the very deep XMM-Newton observations available for an area of 0.2 deg^2. The source detection was performed on the Rc, z' and B band images and the available photometry is spanning from the far ultraviolet to the mid infrared, reaching in the best case scenario 21 bands. Astrometry corrections and photometric cross-calibrations over the entire dataset allowed the computation of accurate photometric redshifts. Special treatment is undertaken for the X-ray sources, the majority of which is active galactic nuclei. Comparing the photometric redshifts to the available spectroscopic redshifts we achieve for normal galaxies an accuracy of \sigma_{\Delta z/(1+z)}=0.036, with 12.7% outliers, while for the X-ray detected sources the accuracy is \sigma_{\Delta z/(1+z)}=0.069, with 18.3% outliers, where the outliers are defined as sources with |z_{phot}-z_{spec}|>0.15 (1+z_{spec})}. These results are a significant improvement over the previously available photometric redshifts for normal galaxies in the Lockman Hole, while it is the first time that photometric redshifts are computed and made public for AGN for this field.

We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

Ultra Steep Spectrum (USS) radio sources have been successfully used to select powerful radio sources at high redshifts (z>~2). Typically restricted to large-sky surveys and relatively bright radio flux densities, it has gradually become possible to extend the USS search to sub-mJy levels, thanks to the recent appearance of sensitive low-frequency radio facilities. Here a first detailed analysis of the nature of the faintest USS sources is presented. By using Giant Metrewave Radio Telescope and Very Large Array radio observations of the Lockman Hole at 610 MHz and 1.4 GHz, a sample of 58 USS sources, with 610 MHz integrated fluxes above 100 microJy, is assembled. Deep infrared data at 3.6 and 4.5 micron from the Spitzer Extragalactic Representative Volume Survey (SERVS) is used to reliably identify counterparts for 48 (83%) of these sources, showing an average total magnitude of [3.6](AB)=19.8 mag. Spectroscopic redshifts for 14 USS sources, together with photometric redshift estimates, improved by the use of the deep SERVS data, for a further 19 objects, show redshifts ranging from z=0.1 to z=2.8, peaking at z~0.6 and tailing off at high redshifts. The remaining 25 USS sources, with no redshift estimate, include the faintest [3.6] magnitudes, with 10 sources undetected at 3.6 and 4.5 micron (typically [3.6]>22-23 mag, from local measurements), which suggests the likely existence of higher redshifts among the sub-mJy USS population. The comparison with the Square Kilometre Array Design Studies Simulated Skies models indicate that Fanaroff-Riley type I radio sources and radio-quiet Active Galactic Nuclei may constitute the bulk of the faintest USS population, and raises the possibility that the high efficiency of the USS technique for the selection of high redshift sources remains even at the sub-mJy level.

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

This paper presents a study of the environments of SDSS Quasi-Stellar-Objects (QSOs) in the Spitzer Extragalactic Representative Volume Survey (SERVS). We concentrate on the high-redshift QSOs as these have not been studied in large numbers with data of this depth before. We use the IRAC 3.6-4.5{\mu}m colour of objects and ancillary r-band data to filter out as much foreground contamination as possible. This technique allows us to find a significant (> 4-{\sigma}) over-density of galaxies around QSOs in a redshift bin centred on z ~ 2.0 and a (> 2-{\sigma}) over-density of galaxies around QSOs in a redshift bin centred on z ~ 3.3. We compare our findings to the predictions of a semi-analytic galaxy formation model, based on the {\Lambda}CDM millennium simulation, and find for both redshift bins that the model predictions match well the source-density we have measured from the SERVS data.

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

We present spectroscopic follow-up observations of 70{\mu}m selected galaxies from the SWIRE XMMLSS and Lockman Hole fields. We have measured spectroscopic redshifts for 293 new sources down to a 70{\mu}m flux limit of 9mJy and R < 22. The redshift distribution peaks at z ~ 0.3 and has a high redshift tail out to z = 3.5. We perform emission line diagnostics for 91 sources where [OIII], H{\beta}, [NII], H{\alpha} and [SII] emission lines are available to determine their power source. We find in our sample 13 QSOs, 1 Seyfert II galaxy, 33 star forming galaxies, 30 composite galaxies, 5 LINERs and 21 ambiguous galaxies. We fit single temperature dust spectral energy distributions (SEDs) to 81 70{\mu}m sources with 160{\mu}m photometry to estimate dust temperatures and masses. Assuming the dust emissivity factor ({\beta}) as 1.5, we determine temperatures in the range ~ 20-60K and dust masses with a range of 10{^6}-10{^9} M{_\odot}. Plotting these objects in the luminosity-temperature diagram suggests that these objects have lower dust temperatures than local IR luminous galaxies. The Herschel Space Observatory will be crucial in understanding the nature of these sources and to accurately determining the shape of the Rayleigh-Jeans tail of the dust SED. We then model SEDs from optical to far-IR for each source using a set of galaxy and quasar templates in the optical and near-IR (NIR) and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst and AGN dust torus) in the mid-IR (MIR) to far-IR (FIR). The number of objects fit with each dust template are: 57 Arp 220, 127 M82, 9 cirrus, 1 AGN dust torus, 70 M82 and cirrus, 26 M82 and AGN dust torus and 3 Arp 220 and AGN dust torus. We determine the total IR luminosity (LIR) in range 10{^8}-10{^{15}} L{_\odot} by integrating the SED models from 8 to 1000{\mu}m.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

We study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $\beta$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = \^a H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &\mu&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

An analysis of large-area CO J=3-2 maps from the James Clerk Maxwell Telescope for 12 nearby spiral galaxies reveals low velocity dispersions in the molecular component of the interstellar medium. The three lowest luminosity galaxies show a relatively flat velocity dispersion as a function of radius while the remaining nine galaxies show a central peak with a radial fall-off within 0.2-0.4 r(25). Correcting for the average contribution due to the internal velocitydispersions of a population of giant molecular clouds, the average cloud-cloud velocity dispersion across the galactic disks is 6.1 +/- 1.0 km/s (standard deviation 2.9 km/s), in reasonable agreement with previous measurements for the Galaxy andM33. The cloud-cloud velocity dispersion derived from the CO data is on average two times smaller than the HI velocity dispersion measured in the same galaxies. The low cloud-cloudvelocity dispersion implies that the molecular gas is the critical component determining the stability of the galactic disk against gravitational collapse, especially in those regions of the disk which are H2 dominated. The cloud-cloud velocity dispersion shows a significant positivecorrelation with both the far-infrared luminosity, which traces the star formation activity, and the K-band absolute magnitude, which traces the total stellar mass. For three galaxies in the Virgo cluster, smoothing the data to a resolution of 4.5 kpc (to match the typical resolution of high redshift CO observations) increases the measured velocity dispersion by roughly a factor of two, comparable to the dispersion measured recently in a normal galaxy at z=1. This comparison suggests that the mass and star formation rate surface densities may be similar in galaxies from z=0-1 and that the high star formation rates seen at z=1 may be partly due to the presence of physically larger molecular gas disks.

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We describe the reduction of data taken with the PACS instrument on board the Herschel Space Observatory in the Science Demonstration Phase of the Herschel-ATLAS (H-ATLAS) survey, specifically data obtained for a 4x4-deg^2 region using Herschel's fast-scan (60 arcsec/s) parallel mode. We describe in detail a pipeline for data reduction using customised procedures within HIPE from data retrieval to the production of science-quality images. We found that the standard procedure for removing Cosmic-Ray glitches also removed parts of bright sources and so implemented an effective two-stage process to minimise these problems. The pronounced 1/f noise is removed from the timelines using 3.4- and 2.5-arcmin boxcar high-pass filters at 100 and 160-um. Empirical measurements of the point-spread function (PSF) are used to determine the encircled energy fraction as a function of aperture size. For the 100- and 160-um bands, the effective PSFs are ~9 and ~13 arcsec (FWHM), and the 90-per-cent encircled energy radii are 13 and 18 arcsec. Astrometric accuracy is good to ~<2 arcsec. The noise in the final maps is correlated between neighbouring pixels and rather higher than advertised prior to launch. For a pair of cross-scans, the 5-sigma point-source sensitivities are 125-165 mJy for 9-13-arcsec radius apertures at 100-um and 150-240 mJy for 13-18-arcsec radius apertures at 160-um.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

By combining Herschel-SPIRE data with archival Spitzer, HI, and CO maps, we investigate the spatial distribution of gas and dust in the two famous grand-design spirals M99 and M100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SPIRE photometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices f350/f500 and f250/f350 decrease radially consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of both galaxies.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

We present Herschel-SPIRE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, supporting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investigations. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of environmental effects on galaxy evolution.