In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

The main limiting factor of cosmological analyses based on thermal Sunyaev-Zel'dovich (SZ) cluster statistics comes from the bias and systematic uncertainties that affect the estimates of the mass of galaxy clusters. High-angular resolution SZ observations at high redshift are needed to study a potential redshift or morphology dependence of both the mean pressure profile and of the mass-observable scaling relation used in SZ cosmological analyses. The NIKA2 camera is a new generation continuum instrument installed at the IRAM 30-m telescope. With a large field of view, a high angular resolution and a high-sensitivity, the NIKA2 camera has unique SZ mapping capabilities. In this paper, we present the NIKA2 SZ large program, aiming at observing a large sample of clusters at redshifts between 0.5 and 0.9, and the characterization of the first cluster oberved with NIKA2.

NIKA2 is a dual-band millimetric continuum camera of 2900 Kinetic Inductance Detectors (KID), operating at $150$ and $260\,\rm{GHz}$, installed at the IRAM 30-meter telescope. We present the performance assessment of NIKA2 after one year of observation using a dedicated point-source calibration method, referred to as the \emph{baseline} method. Using a large data set acquired between January 2017 and February 2018 that span the whole range of observing elevations and atmospheric conditions encountered at the IRAM 30-m telescope, we test the stability of the performance parameters. We report an instantaneous field of view (FOV) of 6.5' in diameter, filled with an average fraction of $84\%$ and $90\%$ of valid detectors at $150$ and $260\,\rm{GHz}$, respectively. The beam pattern is characterized by a FWHM of $17.6'' \pm 0.1''$ and $11.1''\pm 0.2''$, and a beam efficiency of $77\% \pm 2\%$ and $55\% \pm 3\%$ at $150$ and $260\,\rm{GHz}$, respectively. The rms calibration uncertainties are about $3\%$ at $150\,\rm{GHz}$ and $6\%$ at $260\,\rm{GHz}$. The absolute calibration uncertainties are of $5\%$ and the systematic calibration uncertainties evaluated at the IRAM 30-m reference Winter observing conditions are below $1\%$ in both channels. The noise equivalent flux density (NEFD) at $150$ and $260\,\rm{GHz}$ are of $9 \pm 1\, \rm{mJy}\cdot s^{1/2}$ and $30 \pm 3\, \rm{mJy}\cdot s^{1/2}$. This state-of-the-art performance confers NIKA2 with mapping speeds of $1388 \pm 174$ and $111 \pm 11 \,\rm{arcmin}^2\cdot \rm{mJy}^{-2}\cdot \rm{h}^{-1}$ at $150$ and $260\,\rm{GHz}$. With these unique capabilities of fast dual-band mapping at high (better that 18'') angular resolution, NIKA2 is providing an unprecedented view of the millimetre Universe.

The NIKA2 polarization channel at 260 GHz (1.15 mm) has been proposed primarily to observe galactic star-forming regions and probe the critical scales between 0.01-0.05 pc at which magnetic field lines may channel the matter of interstellar filaments into growing dense cores. The NIKA2 polarimeter consists of a room temperature continuously rotating multi-mesh HWP and a cold polarizer that separates the two orthogonal polarizations onto two 260 GHz KIDs arrays. We describe in this paper the preliminary results obtained during the most recent commissioning campaign performed in December 2018. We concentrate here on the analysis of the extended sources, while the observation of compact sources is presented in a companion paper [12]. We present preliminary NIKA2 polarization maps of the Crab nebula. We find that the integrated polarization intensity flux measured by NIKA2 is consistent with expectations.In terms of polarization angle, we are still limited by systematic uncertainties that will be further investigated in the forthcoming commissioning campaigns.

We present here the polarization system of the NIKA camera and give a summary of the main results obtained and performed studies on Orion and the Crab nebula. The polarization system was equipped with a room temperature continuously rotating multi-mesh half wave plate and a grid polarizer facing the NIKA cryostat window. NIKA even though less sensitive than NIKA2 had polarization capability in both 1 and 2 millimiter bands. NIKA polarization observations demonstrated the ability of such a technology in detecting the polarization of different targets, compact and extended sources like the Crab nebula and Orion Molecular Cloud region OMC-1. These measurements together with the developed techniques to deal with systematics, opened the way to the current observations of NIKA2 in polarization that will provide important advances in the studies of galactic and extra-galactic emission and magnetic fields.

The new NIKA2 camera at the IRAM 30m radiotelescope was used to observe three known debris disks in order to constrain the SED of their dust emission in the millimeter wavelength domain. We have found that the spectral index between the two NIKA2 bands (1mm and 2mm) is consistent with the Rayleigh-Jeans regime (lambda^{-2}), unlike the steeper spectra (lambda^{-3}) measured in the submillimeter-wavelength domain for two of the three disks $-$ around the stars Vega and HD107146. We provide a succesful proof of concept to model this spectral inversion in using two populations of dust grains, those smaller and those larger than a grain radius a0 of 0.5mm. This is obtained in breaking the slope of the size distribution and the functional form of the absorption coefficient of the standard model at a0. The third disk - around the star HR8799 - does not exhibit this spectral inversion but is also the youngest.

The steps taken to tailor to NIKA2 observations the Scanamorphos algorithm (initially developed to subtract low-frequency noise from Herschel on-the-fly observations) are described, focussing on the consequences of the different instrument architecture and observation strategy. The method, making the most extensive use of the redundancy built in the multi-scan coverage with large arrays of a given region of the sky, is applicable to extended sources, while the pipeline is so far optimized for compact sources. An example of application is given. A related tool to build consistent broadband SEDs from 60 microns to 2 mm, combining Herschel and NIKA2 data, has also been developed. Its main task is to process the data least affected by low-frequency noise and coverage limitations (i.e. the Herschel data) through the same transfer function as the NIKA2 data, simulating the same scan geometry and applying the same noise and atmospheric signal as extracted from the 1 mm and 2 mm data.

Clusters of galaxies, the largest bound objects in the Universe, constitute a cosmological probe of choice, which is sensitive to both dark matter and dark energy. Within this framework, the Sunyaev-Zel'dovich (SZ) effect has opened a new window for the detection of clusters of galaxies and for the characterization of their physical properties such as mass, pressure and temperature. NIKA, a KID-based dual band camera installed at the IRAM 30-m telescope, was particularly well adapted in terms of frequency, angular resolution, field-of-view and sensitivity, for the mapping of the thermal and kinetic SZ effect in high-redshift clusters. In this paper, we present the NIKA cluster sample and a review of the main results obtained via the measurement of the SZ effect on those clusters: reconstruction of the cluster radial pressure profile, mass, temperature and velocity.

Understanding where and when the mass of stars is determined is one of the fundamental, mostly unsolved, questions in astronomy. Here, we present the first results of GASTON, the Galactic Star Formation with NIKA2 large programme on the IRAM 30m telescope, that aims to identify new populations of low-brightness sources to tackle the question of stellar mass determination across all masses. In this paper, we focus on the high-mass star formation part of the project, for which we map a $\sim2$ deg$^2$ region of the Galactic plane around $l=24^\circ$ in both 1.2 mm and 2.0 mm continuum. Half-way through the project, we reach a sensitivity of 3.7 mJy/beam at 1.2mm. Even though larger than our target sensitivity of 2 mJy, the current sensitivity already allows the identification of a new population of cold, compact sources that remained undetected in any (sub-)mm Galactic plane survey so far. In fact, about 25% of the $\sim 1600$ compact sources identified in the 1.2 mm GASTON image are new detections. We present a preliminary analysis of the physical properties of the GASTON sources as a function of their evolutionary stage, arguing for a potential evolution of the mass distribution of these sources with time.

We present the results of the analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ based on high angular resolution NIKA2 Sunyaev-Zel'dovich (SZ) and $Chandra$ X-ray data. This multi-wavelength analysis enables us to estimate the shape of the temperature profile with unprecedented precision at this redshift and to obtain a map of the gas entropy distribution averaged along the line of sight. The comparison between the cluster morphological properties observed in the NIKA2 and $Chandra$ maps together with the analysis of the entropy map allows us to conclude that MOO J1142$+$1527 is an on-going merger hosting a cool-core at the position of the X-ray peak. This work demonstrates how the addition of spatially-resolved SZ observations to low signal-to-noise X-ray data can bring valuable insights on the intracluster medium thermodynamic properties at $z>1$.

High-resolution mapping of the hot gas in galaxy clusters is a key tool for cluster-based cosmological analyses. Taking advantage of the NIKA2 millimeter camera operated at the IRAM 30-m telescope, the NIKA2 SZ Large Program seeks to get a high-resolution follow-up of 45 galaxy clusters covering a wide mass range at high redshift in order to re-calibrate some of the tools needed for the cosmological exploitation of SZ surveys. We present the second cluster analysis of this program, targeting one of the faintest sources of the sample in order to tackle the difficulties in data reduction for such faint, low-SNR clusters. In this study, the main challenge is the precise estimation of the contamination by sub-millimetric point sources, which greatly affects the tSZ map of the cluster. We account for this contamination by performing a joint fit of the SZ signal and of the flux density of the compact sources. A prior knowledge of these fluxes is given by the adjustment of the SED of each source using data from both NIKA2 and the \textit{Herschel} satellite. The first results are very promising and demonstrate the possibility to estimate thermodynamic properties with NIKA2, even in a compact cluster heavily contaminated by point sources.

The characterization of the Intra-Cluster Medium (ICM) properties of high-redshift galaxy clusters is fundamental to our understanding of large-scale structure formation processes. We present the results of a multi-wavelength analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ discovered as part of the Massive and Distant Clusters of WISE Survey (MaDCoWS). This analysis is based on high angular resolution $Chandra$ X-ray and NIKA2 Sunyaev-Zel'dovich (SZ) data. Although the X-ray data have only about 1700 counts, we are able to determine the ICM thermodynamic radial profiles, namely temperature, entropy, and hydrostatic mass. These have been obtained with unprecedented precision at this redshift and up to $0.7R_{500}$, thanks to the combination of high-resolution X-ray and SZ data. The comparison between the galaxy distribution mapped in infrared by $Spitzer$ and the morphological properties of the ICM derived from the combined analysis of the $Chandra$ and NIKA2 data leads us to the conclusion that the cluster is an on-going merger. We measure the hydrostatic mass profile of the cluster in four angular sectors centered on the large-scale X-ray centroid. This allows us to estimate a systematic uncertainty on the cluster total mass that characterizes both the impact of the observed deviations from spherical symmetry and of the core dynamics on the mass profile. We further combine the X-ray and SZ data at the pixel level to obtain maps of the temperature and entropy distributions averaged along the line of sight. We find a relatively low entropy core at the position of the X-ray peak and high temperature regions located on its south and west sides. The increase in ICM temperature at the location of the SZ peak is expected given the merger dynamics. (abridged)

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

We present new ALMA dust continuum observations of 101 $\log(\mathrm{M}_* / \mathrm{M}_\odot)$ > 9.5 galaxies in the COSMOS field to study the effect of environment on the interstellar medium at z ~ 0.7. At this redshift, our targets span a wide range of environments allowing for a diverse sample of galaxies with densities, $\Sigma$ = 0.16-10.5 Mpc$^{-2}$ (per $\Delta$ z = 0.024). Using the ALMA observations, we calculate the total ISM mass and look for depletion as a function of galaxy density in order to understand the quenching or triggering of star formation in galaxies in different environments. ISM mass is found to have a small dependence on environment, while the depletion timescale remains constant (~200 Myrs) across all environments. We find elevated ISM mass values at intermediate densities and lower values at high densities compared to low (field) densities. Our observed evolution in gas fraction with density in this single redshift slice is equivalent to the observed evolution with cosmic time over 2-3 Gyr. To explain the change in gas mass fraction seen in galaxies in intermediate and high densities, these results suggest environmental processes such as mergers and ram pressure stripping are likely playing a role in dense filamentary-cluster environments.

NIKA2 is a dual-band millimetric camera of thousands of Kinetic Inductance Detectors (KID) installed at the IRAM 30-meter telescope in the Spanish Sierra Nevada. The instrument commissioning was completed in September 2017, and NIKA2 is now open to the scientific community and will operate for the next decade. NIKA2 has well-adapted instrumental design and performance to produce high-resolution maps of the thermal Sunyaev-Zel'dovich (SZ) effect toward intermediate and high redshift galaxy clusters. Moreover, it benefits from a guaranteed time large program dedicated to mapping a representative sample of galaxy clusters via SZ and that includes X-ray follow-ups. The main expected outputs of the SZ large program are the constraints on the redshift evolution of the pressure profile and the mass-observable relation. The first SZ mapping of a galaxy cluster with NIKA2 was produced, as part of the SZ large program. We found a sizable impact of the intracluster medium dynamics on the integrated SZ observables. This shows NIKA2 capabilities for the precise characterisation of the mass-observable relation that is required for accurate cosmology with galaxy clusters.

We study the abundance of star-forming and quiescent galaxies in a sample of 21 massive clusters at 0.5<z<0.7, detected with the Planck satellite. We measure the cluster galaxy stellar mass function (SMF), which is a fundamental observable to study and constrain the formation and evolution of galaxies. Our measurements are based on homogeneous and deep multi-band photometry spanning u- to the Ks-band for each cluster and are supported by spectroscopic data from different programs. The galaxy population is separated between quiescent and star-forming galaxies based on their rest-frame U-V and V-J colours. The SMF is compared to that of field galaxies at the same redshifts, using data from the COSMOS/UltraVISTA survey. We find that the shape of the SMF of star-forming galaxies does not depend on environment, while the SMF of quiescent galaxies has a significantly steeper low-mass slope in the clusters compared to the field. We estimate the environmental quenching efficiency (f_EQ), i.e. the probability for a galaxy that would normally be star forming in the field, to be quenched due to its environment. The f_EQ shows no stellar-mass dependence in any environment, but it increases from 40% in the cluster outskirts to ~90% in the cluster centres. The radial signature of f_EQ provides constraints on where the dominant quenching mechanism operates in these clusters and on what timescale. Exploring these using a simple model based on galaxy orbits obtained from an N-body simulation, we find a clear degeneracy between both parameters. For example, the quenching process may either be triggered on a long (~3 Gyr) time scale at large radii (r~8R_500), or happen well within 1 Gyr at r<R_500. The radius where quenching is triggered is at least r_quench> 0.67R_500 (95%CL). The ICM density at this location suggests that ram-pressure stripping of the cold gas is a likely cause of quenching. [Abridged]

The complete characterization of the pressure profile of high-redshift galaxy clusters, from their core to their outskirts, is a major issue for the study of the formation of large-scale structures. It is essential to constrain a potential redshift evolution of both the slope and scatter of the mass-observable scaling relations used in cosmology studies based on cluster statistics. In this paper, we present the first thermal Sunyaev-Zel'dovich (tSZ) mapping of a cluster from the sample of the NIKA2 SZ large program that aims at constraining the redshift evolution of cluster pressure profiles and the tSZ-mass scaling relation. We have observed the galaxy cluster PSZ2 G144.83+25.11 at redshift $z=0.58$ with the NIKA2 camera, a dual-band (150 and 260 GHz) instrument operated at the IRAM 30-meter telescope. We identify a thermal pressure excess in the south-west region of PSZ2 G144.83+25.11 and a high redshift sub-millimeter point source that affect the intracluster medium (ICM) morphology of the cluster. The NIKA2 data are used jointly with tSZ data acquired by the MUSTANG, Bolocam and $Planck$ experiments in order to non-parametrically set the best constraints on the electronic pressure distribution from the cluster core ($\rm{R} \sim 0.02 \rm{R_{500}}$) to its outskirts ($\rm{R} \sim 3 \rm{R_{500}} $). We investigate the impact of the over-pressure region on the shape of the pressure profile and on the constraints on the integrated Compton parameter $\rm{Y_{500}}$. A hydrostatic mass analysis is also performed by combining the tSZ-constrained pressure profile with the deprojected electronic density profile from XMM-$Newton$. This allows us to conclude that the estimates of $\rm{Y_{500}}$ and $\rm{M_{500}}$ obtained from the analysis with and without masking the disturbed ICM region differ by 65 and 79% respectively. (abridged)

The thermal emission of dust grains is a powerful tool for probing cold, dense regions of molecular gas in the ISM, and so constraining dust properties is key to obtaining accurate measurements of dust mass and temperature. By placing constraints on the dust emissivity spectral index, beta, towards two star-forming infrared dark clouds, SDC18.888 and SDC24.489, we evaluate the role of mass concentration in the associated star-formation activity. We exploit the simultaneous 1.2mm and 2.0mm imaging capability of NIKA on the IRAM 30m telescope to construct maps of beta for both clouds, and by incorporating Herschel observations, we create H2 column density maps with 13" resolution. While we find no significant systematic radial variations around the most massive clumps in either cloud on >0.1 pc scales, their mean beta values are significantly different, with beta = 2.07 +/- 0.09 (rand) +/- 0.25 (syst) for SDC18.888 and beta = 1.71 +/- 0.09 (rand) +/- 0.25 (syst) for SDC24.489. These differences could be a consequence of the very different environments in which both clouds lie, and we suggest that the proximity of SDC18.888 to the W39 HII region may raise beta on scales of 1 pc. We also find that the mass in SDC24.489 is more centrally concentrated and circularly symmetric than in SDC18.888, and is consistent with a scenario in which spherical globally-collapsing clouds concentrate a higher fraction of their mass into a single core than elongated clouds that will more easily fragment, distributing their mass into many cores. We demonstrate that beta variations towards interstellar clouds can be robustly constrained with high-SNR NIKA observations, providing more accurate estimates of their masses. The methods presented here will be applied to the Galactic Star Formation with NIKA2 (GASTON) large programme, extending our analysis to a statistically significant sample of star-forming clouds.

The New IRAM KID Arrays 2 (NIKA2) consortium has just finished installing and commissioning a millimetre camera on the IRAM 30 m telescope. It is a dual-band camera operating with three frequency multiplexed kilo-pixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK, designed to observe the intensity and polarisation of the sky at 260 and 150 GHz (1.15 and 2 mm). NIKA2 is today an IRAM resident instrument for millimetre astronomy, such as Intra Cluster Medium from intermediate to distant clusters and so for the follow-up of Planck satellite detected clusters, high redshift sources and quasars, early stages of star formation and nearby galaxies emission. We present an overview of the instrument performance as it has been evaluated at the end of the commissioning phase.

Millimeter-wave continuum astronomy is today an indispensable tool for both general Astrophysics studies and Cosmology. General purpose, large field-of-view instruments are needed to map the sky at intermediate angular scales not accessible by the high-resolution interferometers and by the coarse angular resolution space-borne or ground-based surveys. These instruments have to be installed at the focal plane of the largest single-dish telescopes. In this context, we have constructed and deployed a multi-thousands pixels dual-band (150 and 260 GHz, respectively 2mm and 1.15mm wavelengths) camera to image an instantaneous field-of-view of 6.5arc-min and configurable to map the linear polarization at 260GHz. We are providing a detailed description of this instrument, named NIKA2 (New IRAM KID Arrays 2), in particular focusing on the cryogenics, the optics, the focal plane arrays based on Kinetic Inductance Detectors (KID) and the readout electronics. We are presenting the performance measured on the sky during the commissioning runs that took place between October 2015 and April 2017 at the 30-meter IRAM (Institut of Millimetric Radio Astronomy) telescope at Pico Veleta. NIKA2 has been successfully deployed and commissioned, performing in-line with the ambitious expectations. In particular, NIKA2 exhibits FWHM angular resolutions of around 11 and 17.5 arc-seconds at respectively 260 and 150GHz. The NEFD (Noise Equivalent Flux Densities) demonstrated on the maps are, at these two respective frequencies, 33 and 8 mJy*sqrt(s). A first successful science verification run has been achieved in April 2017. The instrument is currently offered to the astronomical community during the coming winter and will remain available for at least the next ten years.

Galaxy clusters constitute a major cosmological probe. However, Planck 2015 results have shown a weak tension between CMB-derived and cluster-derived cosmological parameters. This tension might be due to poor knowledge of the cluster mass and observable relationship. As for now, arcmin resolution Sunyaev-Zeldovich (SZ) observations ({\it e.g.} SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters ($z<0.2$). For high redshift clusters ($z>0.5$) high resolution and high sensitivity SZ observations are needed. With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain) is particularly well adapted for these observations. The NIKA2 SZ observation program will map a large sample of clusters (50) at redshifts between 0.5 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.}

We carried out targeted ALMA observations of 129 fields in the COSMOS region at 1.25 mm, detecting 152 galaxies at S/N$\geq$5 with an average continuum RMS of 150 $\mu$Jy. These fields represent a S/N-limited sample of AzTEC / ASTE sources with 1.1 mm S/N$\geq$4 over an area of 0.72 square degrees. Given ALMA's fine resolution and the exceptional spectroscopic and multiwavelength photometric data available in COSMOS, this survey allows us unprecedented power in identifying submillimeter galaxy counterparts and determining their redshifts through spectroscopic or photometric means. In addition to 30 sources with prior spectroscopic redshifts, we identified redshifts for 113 galaxies through photometric methods and an additional nine sources with lower limits, which allowed a statistically robust determination of the redshift distribution. We have resolved 33 AzTEC sources into multi-component systems and our redshifts suggest that nine are likely to be physically associated. Our overall redshift distribution peaks at $z\sim$2.0 with a high redshift tail skewing the median redshift to $\tilde{z}$=2.48$\pm$0.05. We find that brighter millimeter sources are preferentially found at higher redshifts. Our faintest sources, with S$_{1.25 \rm mm}$<1.25 mJy, have a median redshift of $\tilde{z}$=2.18$\pm$0.09, while the brightest sources, S$_{1.25 \rm mm}$>1.8 mJy, have a median redshift of $\tilde{z}$=3.08$\pm$0.17. After accounting for spectral energy distribution shape and selection effects these results are consistent with several previous submillimeter galaxy surveys, and moreover, support the conclusion that the submillimeter galaxy redshift distribution is sensitive to survey depth.

The development of precision cosmology with clusters of galaxies requires high-angular resolution Sunyaev-Zel'dovich (SZ) observations. As for now, arcmin resolution SZ observations (e.g. SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters (z<0.2). With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain), will bring valuable information in the field of SZ imaging of clusters of galaxies. The NIKA2 SZ observation program will allow us to observe a large sample of clusters (50) at redshifts between 0.4 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.

Context: In the past decade, sensitive, resolved Sunyaev-Zel'dovich (SZ) studies of galaxy clusters have become common. Whereas many previous SZ studies have parameterized the pressure profiles of galaxy clusters, non-parametric reconstructions will provide insights into the thermodynamic state of the intracluster medium (ICM). Aims: We seek to recover the non-parametric pressure profiles of the high redshift ($z=0.89$) galaxy cluster CLJ 1226.9+3332 as inferred from SZ data from the MUSTANG, NIKA, Bolocam, and Planck instruments, which all probe different angular scales. Methods: Our non-parametric algorithm makes use of logarithmic interpolation, which under the assumption of ellipsoidal symmetry is analytically integrable. For MUSTANG, NIKA, and Bolocam we derive a non-parametric pressure profile independently and find good agreement among the instruments. In particular, we find that the non-parametric profiles are consistent with a fitted gNFW profile. Given the ability of Planck to constrain the total signal, we include a prior on the integrated Compton Y parameter as determined by Planck. Results: For a given instrument, constraints on the pressure profile diminish rapidly beyond the field of view. The overlap in spatial scales probed by these four datasets is therefore critical in checking for consistency between instruments. By using multiple instruments, our analysis of CLJ 1226.9+3332 covers a large radial range, from the central regions to the cluster outskirts: $0.05 R_{500} < r < 1.1 R_{500}$. This is a wider range of spatial scales than is typical recovered by SZ instruments. Similar analyses will be possible with the new generation of SZ instruments such as NIKA2 and MUSTANG2.

We determine the physical properties of a sample of SMGs in the COSMOS field that were pre-selected at the observed wavelength of $\lambda_{\rm obs}=1.1$ mm, and followed up at $\lambda_{\rm obs}=1.3$ mm with ALMA. We used MAGPHYS to fit the panchromatic (ultraviolet to radio) SEDs of 124 of the target SMGs, 19.4% of which are spectroscopically confirmed. The SED analysis was complemented by estimating the gas masses of the SMGs by using the $\lambda_{\rm obs}=1.3$ mm emission as a tracer of the molecular gas. The sample median and 16th-84th percentile ranges of the stellar masses, SFRs, dust temperatures, and dust and gas masses were derived to be $\log(M_{\star}/{\rm M}_{\odot})=11.09^{+0.41}_{-0.53}$, ${\rm SFR}=402^{+661}_{-233}$ ${\rm M}_{\odot}~{\rm yr}^{-1}$, $T_{\rm dust}=39.7^{+9.7}_{-7.4}$ K, $\log(M_{\rm dust}/{\rm M}_{\odot})=9.01^{+0.20}_{-0.31}$, and $\log(M_{\rm gas}/{\rm M}_{\odot})=11.34^{+0.20}_{-0.23}$, respectively. The median gas-to-dust ratio and gas fraction were found to be $120^{+73}_{-30}$ and $0.62^{+0.27}_{-0.23}$, respectively. We found that 57.3% of our SMGs populate the main sequence (MS) of star-forming galaxies, while 41.9% of the sources lie above the MS by a factor of >3 (one source lies below the MS). The largest 3 GHz radio sizes are found among the MS sources. Those SMGs that appear irregular in the rest-frame UV are predominantly starbursts, while the MS SMGs are mostly disk-like. The larger radio-emitting sizes of the MS SMGs compared to starbursts is a likely indication of their more widespread, less intense star formation. The irregular UV morphologies of the starburst SMGs are likely to echo their merger nature. Our results suggest that the transition from high-$z$ SMGs to local ellipticals via compact, quiescent galaxies (cQGs) at $z \sim 2$ might not be universal, and the latter population might also descend from the so-called blue nuggets.

We measure the stellar mass function (SMF) of galaxies in the COSMOS field up to $z\sim6$. We select them in the near-IR bands of the COSMOS2015 catalogue, which includes ultra-deep photometry from UltraVISTA-DR2, SPLASH, and Subaru/Hyper-SuprimeCam. At $z>2.5$ we use new precise photometric redshifts with error $\sigma_z=0.03(1+z)$ and an outlier fraction of $12\%$, estimated by means of the unique spectroscopic sample of COSMOS. The increased exposure time in the DR2, along with our panchromatic detection strategy, allow us to improve the stellar mass completeness at high $z$ with respect to previous UltraVISTA catalogues. We also identify passive galaxies through a robust colour-colour selection, extending their SMF estimate up to $z=4$. Our work provides a comprehensive view of galaxy stellar mass assembly between $z=0.1$ and 6, for the first time using consistent estimates across the entire redshift range. We fit these measurements with a Schechter function, correcting for Eddington bias. We compare the SMF fit with the halo mass function predicted from $\Lambda$CDM simulations. We find that at $z>3$ both functions decline with a similar slope in the high-mass end. This feature could be explained assuming that the mechanisms that quench star formation in massive haloes become less effective at high redshift; however further work needs to be done to confirm this scenario. Concerning the SMF low-mass end, it shows a progressive steepening as moving towards higher redshifts, with $\alpha$ decreasing from $-1.47_{-0.02}^{+0.02}$ at $z\simeq0.1$ to $-2.11_{-0.13}^{+0.30}$ at $z\simeq5$. This slope depends on the characterisation of the observational uncertainties, which is crucial to properly remove the Eddington bias. We show that there is currently no consensus on the method to quantify such errors: different error models result in different best-fit Schechter parameters. [Abridged]

ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.

We determine the radio size distribution of a large sample of 152 SMGs in COSMOS that were detected with ALMA at 1.3 mm. For this purpose, we used the observations taken by the VLA-COSMOS 3 GHz Large Project. One hundred and fifteen of the 152 target SMGs were found to have a 3 GHz counterpart. The median value of the major axis FWHM at 3 GHz is derived to be $4.6\pm0.4$ kpc. The radio sizes show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be $\alpha=-0.67$ and $T_{\rm B}=12.6\pm2$ K. Three of the target SMGs, which are also detected with the VLBA, show clearly higher brightness temperatures than the typical values. Although the observed radio emission appears to be predominantly powered by star formation and supernova activity, our results provide a strong indication of the presence of an AGN in the VLBA and X-ray-detected SMG AzTEC/C61. The median radio-emitting size we have derived is 1.5-3 times larger than the typical FIR dust-emitting sizes of SMGs, but similar to that of the SMGs' molecular gas component traced through mid-$J$ line emission of CO. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known IR-radio correlation.

The New IRAM KID Array 2 (NIKA 2) is a dual-band camera operating with three frequency-multiplexed kilopixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK. NIKA 2 is designed to observe the intensity and polarisation of the sky at 1.15 and 2.0 mm wavelength from the IRAM 30 m telescope. The NIKA 2 instrument represents a huge step in performance as compared to the NIKA pathfinder instrument, which has already shown state-of-the-art detector and photometric performance. After the commissioning planned to be accomplished at the end of 2016, NIKA 2 will be an IRAM resident instrument for the next ten years or more. NIKA 2 should allow the astrophysical community to tackle a large number of open questions reaching from the role of the Galactic magnetic field in star formation to the discrepancy between cluster-based and CMB-based cosmology possibly induced by the unknown cluster physics. We present an overview of the commissioning phase together with some first results.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

ALMA Cycle 2 observations of the long wavelength dust emission in 145 star-forming galaxies are used to probe the evolution of star-forming ISM. We also develop the physical basis and empirical calibration (with 72 low-z and z ~ 2 galaxies) for using the dust continuum as a quantitative probe of interstellar medium (ISM) masses. The galaxies with highest star formation rates (SFRs) at <z> = 2.2 and 4.4 have gas masses up to 100 times that of the Milky Way and gas mass fractions reaching 50 to 80%, i.e. gas masses 1 - 4 times their stellar masses. We find a single high-z star formation law: SFR = 35 M_ mol^0.89 x (1+z)_{z=2}^0.95 x (sSFR)_{MS}^0.23 \msun yr^-1 -- an approximately linear dependence on the ISM mass and an increased star formation efficiency per unit gas mass at higher redshift. Galaxies above the Main Sequence (MS) have larger gas masses but are converting their ISM into stars on a timescale only slightly shorter than those on the MS -- thus these 'starbursts' are largely the result of having greatly increased gas masses rather than and increased efficiency for converting gas to stars. At z $> 1$, the entire population of star-forming galaxies has $\sim$ 2 - 5 times shorter gas depletion times than low-z galaxies. These shorter depletion times indicate a different mode of star formation in the early universe -- most likely dynamically driven by compressive, high-dispersion gas motions -- a natural consequence of the high gas accretion rates.

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).

ALMA Cycle 2 observations of the long wavelength dust emission in 180 star-forming (SF) galaxies are used to investigate the evolution of ISM masses at z = 1 to 6.4. The ISM masses exhibit strong increases from z = 0 to $\rm <z>$ = 1.15 and further to $\rm <z>$ = 2.2 and 4.8, particularly amongst galaxies above the SF galaxy main sequence (MS). The galaxies with highest SFRs at $\rm <z>$ = 2.2 and 4.8 have gas masses 100 times that of the Milky Way and gas mass fractions reaching 50 to 80\%, i.e. gas masses 1 - 4$\times$ their stellar masses. For the full sample of galaxies, we find a single, very simple SF law: $\rm SFR \propto M_{\rm ISM}^{0.9}$, i.e. a `linear' dependence on the ISM mass -- on and above the MS. Thus, the galaxies above the MS are converting their larger ISM masses into stars on a timescale similar to those on the MS. At z $> 1$, the entire population of star-forming galaxies has $\sim$5 - 10$\times$ shorter gas depletion times ($\sim0.2$ Gyr) than galaxies at low redshift. These {\bf shorter depletion times are due to a different, dominant mode of SF in the early universe} -- dynamically driven by compressive, high dispersion gas motions and/or galaxy interactions. The dispersive gas motions are a natural consequence of the extraordinarily high gas accretion rates which must occur to maintain the prodigious SF.

The Planck catalogue of SZ sources limits itself to a significance threshold of 4.5 to ensure a low contamination rate by false cluster candidates. This means that only the most massive clusters at redshift z>0.5, and in particular z>0.7, are expected to enter into the catalogue, with a large number of systems in that redshift regime being expected around and just below that threshold. In this paper, we follow-up a sample of SZ sources from the Planck SZ catalogues from 2013 and 2015. In the latter maps, we consider detections around and at lower significance than the threshold adopted by the Planck Collaboration. To keep the contamination rate low, our 28 candidates are chosen to have significant WISE detections, in combination with non-detections in SDSS/DSS, which effectively selects galaxy cluster candidates at redshifts $z\gtrsim0.5$. By taking r- and z-band imaging with MegaCam@CFHT, we bridge the 4000A rest-frame break over a significant redshift range, thus allowing accurate redshift estimates of red-sequence cluster galaxies up to z~0.8. After discussing the possibility that an overdensity of galaxies coincides -by chance- with a Planck SZ detection, we confirm that 16 of the candidates have likely optical counterparts to their SZ signals, 13 (6) of which have an estimated redshift z>0.5 (z>0.7). The richnesses of these systems are generally lower than expected given the halo masses estimated from the Planck maps. However, when we follow a simplistic model to correct for Eddington bias in the SZ halo mass proxy, the richnesses are consistent with a reference mass-richness relation established for clusters detected at higher significance. This illustrates the benefit of an optical follow-up, not only to obtain redshift estimates, but also to provide an independent mass proxy that is not based on the same data the clusters are detected with, and thus not subject to Eddington bias.

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$.

Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}\sim$10$^{12.5-14.2}$$M_{\odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $\sim10^{10.4}-10^{10.6}$$M_{\odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $z\sim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.

We present both the observations and the data reduction procedures of the Subaru COSMOS 20 project that is an optical imaging survey of the HST COSMOS field, carried out by using Suprime-Cam on the Subaru Telescope with the following 20 optical filters: 6 broad-band (B, g', V, r', i', and z'), 2 narrow-band (NB711 and NB816), and 12 intermediate-band filters (IA427, IA464, IA484, IA505, IA527, IA574, IA624, IA679, IA709, IA738, IA767, and IA827). A part of this project is described in Taniguchi et al. (2007) and Capak et al. (2007) for the six broad-band and one narrow-band (NB816) filter data. In this paper, we present details of the observations and data reduction for remaining 13 filters (the 12 IA filters and NB711). In particular, we describe the accuracy of both photometry and astrometry in all the filter bands. We also present optical properties of the Suprime-Cam IA filter system in Appendix.

Dust-Obscured galaxies (DOGs) are bright 24 um-selected sources with extreme obscuration at optical wavelengths. They are typically characterized by a rising power-law continuum of hot dust (T_D ~ 200-1000K) in the near-IR indicating that their mid-IR luminosity is dominated by an an active galactic nucleus (AGN). DOGs with a fainter 24 um flux display a stellar bump in the near-IR and their mid-IR luminosity appears to be mainly powered by dusty star formation. Alternatively, it may be that the mid-IR emission arising from AGN activity is dominant but the torus is sufficiently opaque to make the near-IR emission from the AGN negligible with respect to the emission from the host component. In an effort to characterize the astrophysical nature of the processes responsible for the IR emission in DOGs, this paper exploits Herschel data (PACS + SPIRE) on a sample of 95 DOGs within the COSMOS field. We derive a wealth of far-IR properties (e.g., total IR luminosities; mid-to-far IR colors; dust temperatures and masses) based on SED fitting. Of particular interest are the 24 um-bright DOGs (F_24um > 1mJy). They present bluer far-IR/mid-IR colors than the rest of the sample, unveiling the potential presence of an AGN. The AGN contribution to the total 8-1000 um flux increases as a function of the rest-frame 8 um-luminosity irrespective of the redshift. This confirms that faint DOGs (L_8 um< 10^12 L_sun) are dominated by star-formation while brighter DOGs show a larger contribution from an AGN.

We present observations and analysis of a sample of 123 galaxy clusters from the 2013 Planck catalogue of Sunyaev-Zel'dovich sources with the Arcminute Microkelvin Imager (AMI), a ground-based radio interferometer. AMI provides an independent measurement with higher angular resolution, 3 arcmin compared to the Planck beams of 5-10 arcmin. The AMI observations thus provide validation of the cluster detections, improved positional estimates, and a consistency check on the fitted 'size' ($\theta_{s}$) and 'flux' ($Y_{\rm tot}$) parameters in the Generalised Navarro, Frenk and White (GNFW) model. We detect 99 of the clusters. We use the AMI positional estimates to check the positional estimates and error-bars produced by the Planck algorithms PowellSnakes and MMF3. We find that $Y_{\rm tot}$ values as measured by AMI are biased downwards with respect to the Planck constraints, especially for high Planck-SNR clusters. We perform simulations to show that this can be explained by deviation from the 'universal' pressure profile shape used to model the clusters. We show that AMI data can constrain the $\alpha$ and $\beta$ parameters describing the shape of the profile in the GNFW model for individual clusters provided careful attention is paid to the degeneracies between parameters, but one requires information on a wider range of angular scales than are present in AMI data alone to correctly constrain all parameters simultaneously.

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

The relation between the stellar mass and the star formation rate characterizes how the instantaneous star formation is determined by the galaxy past star formation history and by the growth of the dark matter structures. We deconstruct the M-SFR plane by measuring the specific SFR functions in several stellar mass bins from z=0.2 out to z=1.4. Our analysis is primary based on a MIPS 24$\mu m$ selected catalogue combining the COSMOS and GOODS surveys. We estimate the SFR by combining mid- and far-infrared data for 20500 galaxies. The sSFR functions are derived in four stellar mass bins within the range 9.5<log(M/Msun)<11.5. First, we demonstrate the importance of taking into account selection effects when studying the M-SFR relation. Secondly, we find a mass-dependent evolution of the median sSFR with redshift varying as $sSFR \propto (1+z)^{b}$, with $b$ increasing from $b=2.88$ to $b=3.78$ between $M=10^{9.75}Msun$ and $M=10^{11.1}Msun$, respectively. At low masses, this evolution is consistent with the cosmological accretion rate and predictions from semi-analytical models (SAM). This agreement breaks down for more massive galaxies showing the need for a more comprehensive description of the star-formation history in massive galaxies. Third, we obtain that the shape of the sSFR function is invariant with time at z<1.4 but depends on the mass. We observe a broadening of the sSFR function ranging from 0.28 dex at $M=10^{9.75}Msun$ to 0.46 dex at $M=10^{11.1}Msun$. Such increase in the scatter of the M-SFR relation suggests an increasing diversity of SFHs as the stellar mass increases. Finally, we find a gradual decline of the sSFR with mass as $log(sSFR) \propto -0.17M$. We discuss the numerous physical processes, as gas exhaustion in hot gas halos or secular evolution, which can gradually reduce the sSFR and increase the SFH diversity.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr.

The relationship between galaxy star formation rates (SFR) and stellar masses ($M_\ast$) is re-examined using a mass-selected sample of $\sim$62,000 star-forming galaxies at $z \le 1.3$ in the COSMOS 2-deg$^2$ field. Using new far-infrared photometry from $Herschel$-PACS and SPIRE and $Spitzer$-MIPS 24 $\mu$m, along with derived infrared luminosities from the NRK method based on galaxies' locations in the restframe color-color diagram $(NUV - r)$ vs. $(r - K)$, we are able to more accurately determine total SFRs for our complete sample. At all redshifts, the relationship between median $SFR$ and $M_\ast$ follows a power-law at low stellar masses, and flattens to nearly constant SFR at high stellar masses. We describe a new parameterization that provides the best fit to the main sequence and characterizes the low mass power-law slope, turnover mass, and overall scaling. The turnover in the main sequence occurs at a characteristic mass of about $M_{0} \sim 10^{10} M_{\odot}$ at all redshifts. The low mass power-law slope ranges from 0.9-1.3 and the overall scaling rises in SFR as a function of $(1+z)^{4.12 \pm 0.10}$. A broken power-law fit below and above the turnover mass gives relationships of $SFR \propto M_{*}^{0.88 \pm 0.06}$ below the turnover mass and $SFR \propto M_{*}^{0.27 \pm 0.04}$ above the turnover mass. Galaxies more massive than $M_\ast \gtrsim 10^{10}\ M_{\rm \odot}$ have on average, a much lower specific star formation rate (sSFR) than would be expected by simply extrapolating the traditional linear fit to the main sequence found for less massive galaxies.

We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes.

We use percent-level precision photometric redshifts in the UltraVISTA-DR1 near-infrared survey to investigate the changing relationship be- tween galaxy stellar mass and the dark matter haloes hosting them to $z\sim2$. We achieve this by measuring the clustering properties and abundances of a series of volume-limited galaxy samples selected by stellar mass and star-formation activity. We interpret these results in the framework of a phenomenological halo model and numerical simulations. Our measurements span a uniquely large range in stellar mass and redshift and reach below the characteristic stellar mass to $z\sim2$. Our results are: 1. At fixed redshift and scale, clustering amplitude depends monotonically on sample stellar mass threshold; 2. At fixed angular scale, the projected clustering amplitude decreases with redshift but the co-moving correlation length remains constant; 3. Characteristic halo masses and galaxy bias increase with increasing median stellar mass of the sample; 4. The slope of these relationships is modified in lower mass haloes; 5. Concerning the passive galaxy population, characteristic halo masses are consistent with a simply less-abundant version of the full galaxy sample, but at lower redshifts the fraction of satellite galaxies in the passive population is very different from the full galaxy sample; 6. Finally we find that the ratio between the characteristic halo mass and median stellar mass at each redshift bin reaches a peak at ${\log(M_h/M_\odot)}\sim12.2$ and the position of this peak remains constant out to $z\sim2$. The behaviour of the full and passively evolving galaxy samples can be understood qualitatively by considering the slow evolution of the characteristic stellar mass in the redshift range probed by our survey.

The most striking feature of the Cosmic Star Formation History (CSFH) of the Universe is a dramatic drop of the star formation (SF) activity, since z~1. In this work we investigate if the very same process of assembly and growth of structures is one of the major drivers of the observed decline. We study the contribution to the CSFH of galaxies in halos of different masses. This is done by studying the total SFR-halo mass-redshift plane from redshift 0 to redshift z~1.6 in a sample of 57 groups and clusters by using the deepest available mid- and far-infrared surveys conducted with Spitzer MIPS and Herschel PACS and SPIRE. Our results show that low mass groups provide a 60-80% contribution to the CSFH at z~1. Such contribution declines faster than the CSFH in the last 8 billion years to less than 10% at z<0.3, where the overall SF activity is sustained by lower mass halos. More massive systems provide only a marginal contribution (<10%) at any epoch. A simplified abundance matching method shows that the large contribution of low mass groups at z~1 is due to a large fraction (>50%) of very massive, highly star forming Main Sequence galaxies. Below z~1 a quenching process must take place in massive halos to cause the observed faster suppression of their SF activity. Such process must be a slow one though, as most of the models implementing a rapid quenching of the SF activity in accreting satellites significantly underpredicts the observed SF level in massive halos at any redshift. Starvation or the transition from cold to hot accretion would provide a quenching timescale of 1 Gyrs more consistent with the observations. Our results suggest a scenario in which, due to the structure formation process, more and more galaxies experience the group environment and, thus, the associated quenching process. This leads to the progressive suppression of their SF activity shaping the CSFH below z~1.

There is now a large consensus that the current epoch of the Cosmic Star Formation History (CSFH) is dominated by low mass galaxies while the most active phase at 1<z<2 is dominated by more massive galaxies, which undergo a faster evolution. Massive galaxies tend to inhabit very massive halos such as galaxy groups and clusters. We aim to understand whether the observed "galaxy downsizing" could be interpreted as a "halo downsizing", whereas the most massive halos, and their galaxy populations, evolve more rapidly than the halos of lower mass. Thus, we study the contribution to the CSFH of galaxies inhabiting group-sized halos. This is done through the study of the evolution of the Infra-Red (IR) luminosity function of group galaxies from redshift 0 to ~1.6. We use a sample of 39 X-ray selected groups in the Extended Chandra Deep Field South (ECDFS), the Chandra Deep Field North (CDFN), and the COSMOS field, where the deepest available mid- and far-IR surveys have been conducted with Spitzer MIPS and Hersche PACS. Groups at low redshift lack the brightest, rarest, and most star forming IR-emitting galaxies observed in the field. Their IR-emitting galaxies contribute <10% of the comoving volume density of the whole IR galaxy population in the local Universe. At redshift >~1, the most IR-luminous galaxies (LIRGs and ULIRGs) are preferentially located in groups, and this is consistent with a reversal of the star-formation rate vs .density anti-correlation observed in the nearby Universe. At these redshifts, group galaxies contribute 60-80% of the CSFH, i.e. much more than at lower redshifts. Below z~1, the comoving number and SFR densities of IR-emitting galaxies in groups decline significantly faster than those of all IR-emitting galaxies. Our results are consistent with a "halo downsizing" scenario and highlight the significant role of "environment" quenching in shaping the CSFH.

We study the evolution of the total star formation (SF) activity, total stellar mass and halo occupation distribution in massive halos by using one of the largest X-ray selected sample of galaxy groups with secure spectroscopic identification in the major blank field surveys (ECDFS, CDFN, COSMOS, AEGIS). We provide an accurate measurement of SFR for the bulk of the star-forming galaxies using very deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations. For undetected IR sources, we provide a well-calibrated SFR from SED fitting. We observe a clear evolution in the level of SF activity in galaxy groups. The total SF activity in the high redshift groups (0.5<z<1.1) is higher with respect to the low redshift (0.15<z<0.5) sample at any mass by 0.8+/-0.12 dex. A milder difference (0.35+/-0.1 dex) is observed between the low redshift bin and the groups at z~0. We show that the level of SF activity is declining more rapidly in the more massive halos than in the more common lower mass halos. We do not observe any evolution in the halo occupation distribution and total stellar mass- halo mass relations in groups. The picture emerging from our findings suggests that the galaxy population in the most massive systems is evolving faster than galaxies in lower mass halos, consistently with a "halo downsizing" scenario.

Star-forming galaxies (SFGs) display a continuous distribution of specific star formation rates (sSFR) which can be approximated by the superposition of two log-normal distributions. The 1st of these encompasses the main sequence (MS) of SFGs, the 2nd one a rarer population of starbursts (SB). We show that the sSFR-distribution of SBs can be regarded as the result of a physical process (plausibly merging) taking the mathematical form of a log-normal boosting kernel and enhancing star formation activity. We explore the utility of splitting the star-forming population into MS and SB galaxies - an approach we term "2-Star Formation Mode" (2-SFM) framework - for understanding their molecular gas properties. Variations of star formation efficiency (SFE) and gas fractions among SFGs take a simple, redshift-independent form, once these quantities are normalized to the value of an average MS galaxy. The change in SFE for galaxies undergoing a starburst event scales supra-linearly with the SFR-increase, as expected for mergers. This implies a continuous distribution of galaxies in the Schmidt-Kennicutt plane that separates more clearly into loci for SBs and normal galaxies than observed in SFR vs. M* space. SBs with the largest deviations (>10-fold) from the MS, like many local ULIRGs, are not average SBs, but even rarer events with progenitors having larger gas fractions than typical MS galaxies. We statistically infer that the gas fractions of typical SBs decrease by a factor of 2 to 3 with respect to their direct MS progenitors, as expected to occur in short-term SFR-boosts during which internal gas reservoirs are drained more quickly than gas is accreted from the cosmic web. We predict variations of the CO-to-H2 conversion factor in the SFR-M* plane and provide evidence that the higher sSFR of distant galaxies is a direct consequence of larger gas fraction in these systems.

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

We present the activities carried out to calibrate and characterise the performance of the elements of attitude control and measurement on board the Herschel spacecraft. The main calibration parameters and the evolution of the indicators of the pointing performance are described, from the initial values derived from the observations carried out in the performance verification phase to those attained in the last year and half of mission, an absolute pointing error around or even below 1 arcsec, a spatial relative pointing error of some 1 arcsec and a pointing stability below 0.2 arsec. The actions carried out at the ground segment to improve the spacecraft pointing measurements are outlined. On-going and future developments towards a final refinement of the Herschel astrometry are also summarised. A brief description of the different components of the attitude control and measurement system (both in the space and in the ground segments) is also given for reference. We stress the importance of the cooperation between the different actors (scientists, flight dynamics and systems engineers, attitude control and measurement hardware designers, star-tracker manufacturers, etc.) to attain the final level of performance.

Using a sample of 25683 star-forming and 2821 passive galaxies at $z\sim2$, selected in the COSMOS field following the BzK color criterion, we study the hosting halo mass and environment of galaxies as a function of their physical properties. Spitzer and Herschel provide accurate SFR estimates for starburst galaxies. We measure the auto- and cross-correlation functions of various galaxy sub-samples and infer the properties of their hosting halos using both an HOD model and the linear bias at large scale. We find that passive and star-forming galaxies obey a similarly rising relation between the halo and stellar mass. The mean host halo mass of star forming galaxies increases with the star formation rate between 30 and 200 M$_\odot$.yr$^{-1}$, but flattens for higher values, except if we select only main-sequence galaxies. This reflects the expected transition from a regime of secular co-evolution of the halos and the galaxies to a regime of episodic starburst. We find similar large scale biases for main-sequence, passive, and starburst galaxies at equal stellar mass, suggesting that these populations live in halos of the same mass. We detect an excess of clustering on small scales for passive galaxies and showed, by measuring the large-scale bias of close pairs, that this excess is caused by a small fraction ($\sim16%$) of passive galaxies being hosted by massive halos ($\sim 3 \times 10^{13}$ M$_\odot$) as satellites. Finally, extrapolating the growth of halos hosting the z$\sim$2 population, we show that M$_\star \sim 10^{10}$ M$_\odot$ galaxies at z$\sim$2 will evolve, on average, into massive (M$_\star \sim 10^{11}$ M$_\odot$), field galaxies in the local Universe and M$_\star \sim 10^{11}$ M$_\odot$ galaxies at z=2 into local, massive, group galaxies. The most massive main-sequence galaxies and close pairs of massive, passive galaxies end up in today's clusters.

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

The use of submm dust continuum emission to probe the mass of interstellar dust and gas in galaxies is empirically calibrated using samples of local star forming galaxies, Planck observations of the Milky Way and high redshift submm galaxies (SMGs). All of these objects suggest a similar calibration, strongly supporting the view that the Rayleigh-Jeans (RJ) tail of the dust emission can be used as an accurate and very fast probe of the ISM in galaxies. We present ALMA Cycle 0 observations of the Band 7 (350 GHz) dust emission in 107 galaxies from z = 0.2 to 2.5. Three samples of galaxies with a total of 101 galaxies were stellar mass-selected from COSMOS to have $M_* \simeq10^{11}$\msun: 37 at z$\sim0.4$, 33 at z$\sim0.9$ and 31 at z$=2$. A fourth sample with 6 IR luminous galaxies at z = 2 was observed for comparison with the purely mass-selected samples. From the fluxes detected in the stacked images for each sample, we find that the ISM content has decreased a factor $\sim 6$ from $1 - 2 \times 10^{10}$\msun at both z = 2 and 0.9 down to $\sim 2 \times 10^9$\msun at z = 0.4. The IR luminous sample at z = 2 shows a further $\sim 4$ times increase in M$_{ISM}$ compared to the equivalent non-IR bright sample at the same redshift. The gas mass fractions are $\sim 2\pm0.5, 12\pm3, 14\pm2 ~\rm{and} ~53\pm3$ $%$ for the four subsamples (z = 0.4, 0.9, 2 and IR bright galaxies).

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.

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.

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 investigate the evolution of the star formation rate (SFR)-density relation in the Extended Chandra Deep Field South (ECDFS) and the Great Observatories Origin Deep Survey (GOODS) fields up to z~1.6. In addition to the "traditional method", in which the environment is defined according to a statistical measurement of the local galaxy density, we use a "dynamical" approach, where galaxies are classified according to three different environment regimes: group, "filament-like", and field. Both methods show no evidence of a SFR-density reversal. Moreover, group galaxies show a mean SFR lower than other environments up to z~1, while at earlier epochs group and field galaxies exhibit consistent levels of star formation (SF) activity. We find that processes related to a massive dark matter halo must be dominant in the suppression of the SF below z~1, with respect to purely density-related processes. We confirm this finding by studying the distribution of galaxies in different environments with respect to the so-called Main Sequence (MS) of star-forming galaxies. Galaxies in both group and "filament-like" environments preferentially lie below the MS up to z~1, with group galaxies exhibiting lower levels of star-forming activity at a given mass. At z>1, the star-forming galaxies in groups reside on the MS. Groups exhibit the highest fraction of quiescent galaxies up to z~1, after which group, "filament-like", and field environments have a similar mix of galaxy types. We conclude that groups are the most efficient locus for star-formation quenching. Thus, a fundamental difference exists between bound and unbound objects, or between dark matter haloes of different masses.

We combine Herschel PACS and SPIRE maps of the full 2 deg^2 COSMOS field with existing multi-wavelength data to obtain template and model-independent optical-to-far-infrared spectral energy distributions (SEDs) for 4,218 Herschel-selected sources with log(L_IR/L_sun) = 9.4-13.6 and z = 0.02-3.54. Median SEDs are created by binning the optical to far-infrared (FIR) bands available in COSMOS as a function of infrared luminosity. Herschel probes rest-frame wavelengths where the bulk of the infrared radiation is emitted, allowing us to more accurately determine fundamental dust properties of our sample of infrared luminous galaxies. We find that the SED peak wavelength (lambda_peak) decreases and the dust mass (M_dust) increases with increasing total infrared luminosity (L_IR). In the lowest infrared luminosity galaxies (log(L_IR/L_sun) = 10.0-11.5), we see evidence of Polycyclic Aromatic Hydrocarbons (PAH) features (lambda=7-9 um), while in the highest infrared luminosity galaxies (L_IR > 10^12 L_sun) we see an increasing contribution of hot dust and/or power-law emission, consistent with the presence of heating from an active galactic nucleus (AGN). We study the relationship between stellar mass and star formation rate of our sample of infrared luminous galaxies and find no evidence that Herschel-selected galaxies follow the SFR/M_* "main sequence" as previously determined from studies of optically selected, star-forming galaxies. Finally, we compare the mid-infrared (MIR) to FIR properties of our infrared luminous galaxies using the previously defined diagnostic, IR8 = L_IR / L_8, and find that galaxies with L_IR > 10^11.3 L_sun tend to systematically lie above (~3-5) the IR8 "infrared main sequence", suggesting either suppressed PAH emission or an increasing contribution from AGN heating.

We present an empirical method of assessing the star formation rate (SFR) of star-forming galaxies based on their locations in the rest-frame color-color diagram (NUV-r) vs (r-K). By using the Spitzer 24 micron sample in the COSMOS field (~16400 galaxies with 0.2 < z < 1.3) and a local GALEX-SDSS-SWIRE sample (~700 galaxies with z < 0.2), we show that the mean infrared excess <IRX>= < L_IR / L_UV > can be described by a single vector, NRK, that combines the two colors. The calibration between <IRX> and NRK allows us to recover the IR luminosity, L_IR, with an accuracy of ~0.21 dex for the COSMOS sample and ~0.27 dex for the local one. The SFRs derived with this method agree with the ones based on the observed (UV+IR) luminosities and on the spectral energy distribution fitting for the vast majority (~85 %) of the star-forming population. Thanks to a library of model galaxy SEDs with realistic prescriptions for the star formation history, we show that we need to include a two-component dust model (i.e., birth clouds and diffuse ISM) and a full distribution of galaxy inclinations in order to reproduce the behavior of the <IRX> stripes in the NUVrK diagram. In conclusion, the NRK method, based only on rest-frame UV and optical colors available in most of the extragalactic fields, offers a simple alternative of assessing the SFR of star-forming galaxies in the absence of far-IR or spectral diagnostic observations.

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.

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$\beta$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

We explore the role of environment in the evolution of galaxies over 0.1<z<0.7 using the final zCOSMOS-bright data set. Using the red fraction of galaxies as a proxy for the quenched population, we find that the fraction of red galaxies increases with the environmental overdensity and with the stellar mass, consistent with previous works. As at lower redshift, the red fraction appears to be separable in mass and environment, suggesting the action of two processes: mass and environmental quenching. The parameters describing these appear to be essentially the same at z~0.7 as locally. We explore the relation between red fraction, mass and environment also for the central and satellite galaxies separately, paying close attention to the effects of impurities in the central-satellite classification and using carefully constructed samples matched in stellar mass. There is little evidence for a dependence of the red fraction of centrals on overdensity. Satellites are consistently redder at all overdensities, and the satellite quenching efficiency increases with overdensity at 0.1<z<0.4. This is less marked at higher redshift, but both are nevertheless consistent with the equivalent local measurements. At a given stellar mass, the fraction of galaxies that are satellites also increases with the overdensity. At a given overdensity and mass, the obtained relation between the environmental quenching and the satellite fraction agrees well with the satellite quenching efficiency, demonstrating that the environmental quenching in the overall population is consistent with being entirely produced through the satellite quenching process at least up to z=0.7. However, despite the unprecedented size of our high redshift samples, the associated statistical uncertainties are still significant and our statements should be understood as approximations to physical reality, rather than physically exact formulae.

In the local Universe, galaxy properties show a strong dependence on environment. In cluster cores, early type galaxies dominate, whereas star-forming galaxies are more and more common in the outskirts. At higher redshifts and in somewhat less dense environments (e.g. galaxy groups), the situation is less clear. One open issue is that of whether and how the star formation rate (SFR) of galaxies in groups depends on the distance from the centre of mass. To shed light on this topic, we have built a sample of X-ray selected galaxy groups at 0<z<1.6 in various blank fields (ECDFS, COSMOS, GOODS). We use a sample of spectroscopically confirmed group members with stellar mass M >10^10.3 M_sun in order to have a high spectroscopic completeness. As we use only spectroscopic redshifts, our results are not affected by uncertainties due to projection effects. We use several SFR indicators to link the star formation (SF) activity to the galaxy environment. Taking advantage of the extremely deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations, we have an accurate, broad-band measure of the SFR for the bulk of the star-forming galaxies. We use multi-wavelength SED fitting techniques to estimate the stellar masses of all objects and the SFR of the MIPS and PACS undetected galaxies. We analyse the dependence of the SF activity, stellar mass and specific SFR on the group-centric distance, up to z~1.6, for the first time. We do not find any correlation between the mean SFR and group-centric distance at any redshift. We do not observe any strong mass segregation either, in agreement with predictions from simulations. Our results suggest that either groups have a much smaller spread in accretion times with respect to the clusters and that the relaxation time is longer than the group crossing time.

[Abridged] Recent studies revealed a strong correlation between the star formation rate (SFR) and stellar mass of star-forming galaxies, the so-called star-forming main sequence. An empirical modeling approach (2-SFM) which distinguishes between the main sequence and rarer starburst galaxies is capable of reproducing most statistical properties of infrared galaxies. In this paper, we extend this approach by establishing a connection between stellar mass and halo mass with the technique of abundance matching. Based on a few, simple assumptions and a physically motivated formalism, our model successfully predicts the (cross-)power spectra of the cosmic infrared background (CIB), the cross-correlation between CIB and cosmic microwave background (CMB) lensing, and the correlation functions of bright, resolved infrared galaxies measured by Herschel, Planck, ACT and SPT. We use this model to infer the redshift distribution these observables, as well as the level of correlation between CIB-anisotropies at different wavelengths. We also predict that more than 90% of cosmic star formation activity occurs in halos with masses between 10^11.5 and 10^13.5 Msun. Taking into account subsequent mass growth of halos, this implies that the majority of stars were initially (at z>3) formed in the progenitors of clusters, then in groups at 0.5<z<3 and finally in Milky-Way-like halos at z<0.5. At all redshifts, the dominant contribution to the star formation rate density stems from halos of mass ~10^12 Msun, in which the instantaneous star formation efficiency is maximal (~70%). The strong redshift-evolution of SFR in the galaxies that dominate the CIB is thus plausibly driven by increased accretion from the cosmic web onto halos of this characteristic mass scale. Material available at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

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.

By making use of Herschel-PEP observations of the COSMOS and Extended Groth Strip fields, we have estimated the dependence of the clustering properties of FIR-selected sources on their 100um fluxes. Our analysis shows a tendency for the clustering strength to decrease with limiting fluxes: r0(S100um >8 mJy)~4.3 Mpc and r0(S100um >5 mJy)~5.8 Mpc. These values convert into minimum halo masses Mmin~10^{11.6} Msun for sources brighter than 8 mJy and Mmin~10^{12.4} Msun for S100um > 5 mJy galaxies. We show such an increase of the clustering strength to be due to an intervening population of z~2 sources, which are very strongly clustered and whose relative contribution, equal to about 10% of the total counts at S100um > 2 mJy, rapidly decreases for brighter flux cuts. By removing such a contribution, we find that z <~ 1 FIR galaxies have approximately the same clustering properties, irrespective of their flux level. The above results were then used to investigate the intrinsic dependence on cosmic epoch of the clustering strength of dusty star-forming galaxies between z~0 and z~2.5. In order to remove any bias in the selection process, the adopted sample only includes galaxies observed at the same rest-frame wavelength, lambda~60 um, which have comparable luminosities and therefore star-formation rates (SFR>~100 Msun/yr). Our analysis shows that the same amount of (intense) star forming activity takes place in extremely different environments at the different cosmological epochs. For z<~1 the hosts of such star forming systems are small, Mmin~10^{11} Msun, isolated galaxies. High (z~2) redshift star formation instead seems to uniquely take place in extremely massive/cluster-like halos, Mmin~10^{13.5} Msun, which are associated with the highest peaks of the density fluctuation field at those epochs. (abridged)

We present results from the deepest Herschel-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extragalactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and the construction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3-sigma depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 um, respectively, and resolve ~75% of the cosmic infrared background at 100um and 160um into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths and to determine the infrared luminosity function of galaxies down to LIR=10^11 Lsun at z~1 and LIR=10^12 Lsun at z~2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3-sigma) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

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.

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

Galactic black hole binaries produce powerful outflows with emit over almost the entire electromagnetic spectrum. Here, we report the first detection with the Herschel observatory of a variable far-infrared source associated with the compact jets of the black hole transient GX 339-4 during the decay of its recent 2010-2011 outburst, after the transition to the hard state. We also outline the results of very sensitive radio observations conducted with the Australia Telescope Compact Array, along with a series of near-infrared, optical (OIR) and X-ray observations, allowing for the first time the re-ignition of the compact jets to be observed over a wide range of wavelengths. The compact jets first turn on at radio frequencies with an optically thin spectrum that later evolves to optically thick synchrotron emission. An OIR reflare is observed about ten days after the onset of radio and hard X-ray emission, likely reflecting the necessary time to build up enough density, as well as to have acceleration (e.g. through shocks) along an extended region in the jets. The Herschel measurements are consistent with an extrapolation of the radio inverted power-law spectrum, but they highlight a more complex radio to OIR spectral energy distribution for the jets.

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.

(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.

We have compiled a large sample of 151 high redshift (z=0.5-4) galaxies selected at 24 microns (S24>100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-infrared spectrum into contributions from star formation and activity in the galactic nuclei. In addition, we have a wealth of photometric data from Spitzer IRAC/MIPS and Herschel PACS/SPIRE. We explore how effective different infrared color combinations are at separating our mid-IR spectroscopically determined active galactic nuclei from our star forming galaxies. We look in depth at existing IRAC color diagnostics, and we explore new color-color diagnostics combining mid-IR, far-IR, and near-IR photometry, since these combinations provide the most detail about the shape of a source's IR spectrum. An added benefit of using a color that combines far-IR and mid-IR photometry is that it is indicative of the power source driving the IR luminosity. For our data set, the optimal color selections are S250/S24 vs. S8.0/S3.6 and S100/S24 vs. S8.0/S3.6; both diagnostics have ~10% contamination rate in the regions occupied primarily by star forming galaxies and active galactic nuclei, respectively. Based on the low contamination rate, these two new IR color-color diagnostics are ideal for estimating both the mid-IR power source of a galaxy when spectroscopy is unavailable and the dominant power source contributing to the IR luminosity. In the absence of far-IR data, we present color diagnostics using the WISE mid-IR bands which can efficiently select out high z (z~2) star forming galaxies.

We study the far-infrared (IR) and sub-millimeter properties of a sample of ultraviolet (UV) selected galaxies at z\sim1.5. Using stacking at 250, 350 and 500 um from Herschel Space Observatory SPIRE imaging of the COSMOS field obtained within the HerMES key program, we derive the mean IR luminosity as a function of both UV luminosity and slope of the UV continuum beta. The IR to UV luminosity ratio is roughly constant over most of the UV luminosity range we explore. We also find that the IR to UV luminosity ratio is correlated with beta. We observe a correlation that underestimates the correlation derived from low-redshift starburst galaxies, but is in good agreement with the correlation derived from local normal star-forming galaxies. Using these results we reconstruct the IR luminosity function of our UV-selected sample. This luminosity function recovers the IR luminosity functions measured from IR selected samples at the faintest luminosities (Lir ~ 10^{11} L_sun), but might underestimate them at the bright-end (Lir > 5.10^{11} L_sun). For galaxies with 10^{11}<Lir/L_sun<10^{13}, the IR luminosity function of a UV selection recovers (given the differences in IR-based estimates) 52-65 to 89-112 per cent of the star-formation rate density derived from an IR selection. The cosmic star-formation rate density derived from this IR luminosity function is 61-76 to 100-133 per cent of the density derived from IR selections at the same epoch. Assuming the latest Herschel results and conservative stacking measurements, we use a toy model to fully reproduce the far IR luminosity function from our UV selection at z\sim 1.5. This suggests that a sample around 4 magnitudes deeper (i.e. reaching u \sim 30 mag) and a large dispersion of the IR to UV luminosity ratio are required.

We investigate the effect of the "high-pass filter" data reduction technique on the Herschel PACS PSF and noise of the PACS maps at the 70, 100 and 160 um bands and in medium and fast scan speeds. This branch of the PACS Photometer pipeline is the most used for cosmological observations and for point-source observations.The calibration of the flux loss due to the median removal applied by the PACS pipeline (high-pass filter) is done via dedicated simulations obtained by "polluting" real PACS timelines with fake sources at different flux levels. The effect of the data reduction parameter settings on the final map noise is done by using selected observations of blank fields with high data redundancy. We show that the running median removal can cause significant flux losses at any flux level. We analyse the advantages and disadvantages of several masking strategies and suggest that a mask based on putting circular patches on prior positions is the best solution to reduce the amount of flux loss. We provide a calibration of the point-source flux loss for several masking strategies in a large range of data reduction parameters, and as a function of the source flux. We also show that, for stacking analysis, the impact of the high-pass filtering effect is to reduce significantly the clustering effect. The analysis of the global noise and noise components of the PACS maps shows that the dominant parameter in determining the final noise is the high-pass filter width. We also provide simple fitting functions to build the error map from the coverage map and to estimate the cross-correlation correction factor in a representative portion of the data reduction parameter space.

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.

In this paper we measure the merger fraction and rate, both minor and major, of massive early-type galaxies (M_star >= 10^11 M_Sun) in the COSMOS field, and study their role in mass and size evolution. We use the 30-band photometric catalogue in COSMOS, complemented with the spectroscopy of the zCOSMOS survey, to define close pairs with a separation 10h^-1 kpc <= r_p <= 30h-1 kpc and a relative velocity Delta v <= 500 km s^-1. We measure both major (stellar mass ratio mu = M_star,2/M_star,1 >= 1/4) and minor (1/10 <= mu < 1/4) merger fractions of massive galaxies, and study their dependence on redshift and on morphology. The merger fraction and rate of massive galaxies evolves as a power-law (1+z)^n, with major mergers increasing with redshift, n_MM = 1.4, and minor mergers showing little evolution, n_mm ~ 0. When split by their morphology, the minor merger fraction for early types is higher by a factor of three than that for spirals, and both are nearly constant with redshift. Our results show that massive early-type galaxies have undergone 0.89 mergers (0.43 major and 0.46 minor) since z ~ 1, leading to a mass growth of ~30%. We find that mu >= 1/10 mergers can explain ~55% of the observed size evolution of these galaxies since z ~ 1. Another ~20% is due to the progenitor bias (younger galaxies are more extended) and we estimate that very minor mergers (mu < 1/10) could contribute with an extra ~20%. The remaining ~5% should come from other processes (e.g., adiabatic expansion or observational effects). This picture also reproduces the mass growth and velocity dispersion evolution of these galaxies. We conclude from these results that merging is the main contributor to the size evolution of massive ETGs at z <= 1, accounting for ~50-75% of that evolution in the last 8 Gyr. Nearly half of the evolution due to mergers is related to minor (mu < 1/4) events.

We present deep CO observations of NGC6240 performed with the IRAM Plateau de Bure Interferometer (PdBI). NGC6240 is the prototypical example of a major galaxy merger in progress, caught at an early stage, with an extended, strongly-disturbed butterfly-like morphology and the presence of a heavily obscured active nucleus in the core of each progenitor galaxy. The CO line shows a skewed profile with very broad and asymmetric wings detected out to velocities of -600 km/s and +800 km/s with respect to the systemic velocity. The PdBI maps reveal the existence of two prominent structures of blueshifted CO emission. One extends eastward, i.e. approximately perpendicular to the line connecting the galactic nuclei, over scales of ~7 kpc and shows velocities up to -400 km/s. The other extends southwestward out to ~7 kpc from the nuclear region, and has a velocity of -100 km/s with respect to the systemic one. Interestingly, redshifted emission with velocities 400 to 800 km/s is detected around the two nuclei, extending in the east-west direction, and partly overlapping with the eastern blue-shifted structure, although tracing a more compact region of size ~1.7 kpc. The overlap between the southwestern CO blob and the dust lanes seen in HST images, which are interpreted as tidal tails, indicates that the molecular gas is deeply affected by galaxy interactions. The eastern blueshifted CO emission is co-spatial with an Halpha filament that is associated with strong H2 and soft X-ray emission. The analysis of Chandra X-ray data provides strong evidence for shocked gas at the position of the Halpha emission. Its association with outflowing molecular gas supports a scenario where the molecular gas is compressed into a shock wave that propagates eastward from the nuclei. If this is an outflow, the AGN are likely the driving force.

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 explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.

Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared (8um) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS-Herschel survey to show that most DOGs with 10^12 L_Sun < L_IR < 10^13 L_Sun are not abnormally bright in the mid-IR when compared to other dusty galaxies with similar IR (8--1000um) luminosities. We observe a relation between the median IR to UV luminosity ratios and the median UV continuum power-law indices for these galaxies, and we find that only 24% have specific star formation rates which indicate the dominance of compact star-forming regions. This circumstantial evidence supports the idea that the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We attribute differences in UV obscuration to either: 1) differences in the degree of alignment between the spatial distributions of dust and massive stars, or 2) differences in the total dust content.

The "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8\mu m$-- 4000\AA), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.

We present observations of the CO[3-2] emission towards two massive and infrared luminous Lyman Break Galaxies at z = 3.21 and z = 2.92, using the IRAM Plateau de Bure Interferometer, placing first constraints on the molecular gas masses (Mgas) of non-lensed LBGs. Their overall properties are consistent with those of typical (Main-Sequence) galaxies at their redshifts, with specific star formation rates ~1.6 and ~2.2 Gyr^(-1), despite their large infrared luminosities L_IR ~2-3 x 10^12 Lsun derived from Herschel. With one plausible CO detection (spurious detection probability of 10^(-3)) and one upper limit, we investigate the evolution of the molecular gas-to-stellar mass ratio (Mgas/M*) with redshift. Our data suggest that the steep evolution of Mgas/M* of normal galaxies up to z~2 is followed by a flattening at higher redshifts, providing supporting evidence for the existence of a plateau in the evolution of the specific star formation rate at z > 2.5.

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi-wavelength information of the CDFS to assess this issue. We select the AGNs from the 3Ms XMM-Newton survey and measure the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 um. Star-formation rates are obtained from spectral energy distribution fits, identifying and subtracting an AGN component. We divide the star-formation rates by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding Lx=10^43 erg/s and redshifts z~1. We do not find evidence for such a correlation for lower luminosity systems or those at lower redshifts. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming galaxies at the same redshift we find that the AGNs reside mostly in main-sequence and starburst hosts, reflecting the AGN - sSFR connection. Limiting our analysis to the highest X-ray luminosity AGNs (X-ray QSOs with Lx>10^44 erg/s), we find that the highest-redshift QSOs (with z>2) reside predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of QSOs at z~1.5 with specific star-formation rates compatible with the main-sequence, or even in the "quiescent" region. (abridged)

We report results from a deep Jansky Very Large Array (JVLA) search for CO 1-0 line emission from galaxies in a candidate galaxy cluster at z~1.55 in the COSMOS field. We target 4 galaxies with optical spectroscopic redshifts in the range z=1.47-1.59. Two of these 4 galaxies, ID51613 and ID51813, are nominally detected in CO line emission at the 3-4 sigma level. We find CO luminosities of 2.4x10^10 K km/s pc^2 and 1.3x10^10 K km/s pc^2, respectively. Taking advantage from the clustering and 2-GHz bandwidth of the JVLA, we perform a search for emission lines in the proximity of optical sources within the field of view of our observations. We limit our search to galaxies with K<23.5 (AB) and z_phot=1.2-1.8. We find 2 bright optical galaxies to be associated with significant emission line peaks (>4 sigma) in the data cube, which we identify with the CO line emission. To test the reliability of the line peaks found, we performed a parallel search for line peaks using a Bayesian inference method. Monte Carlo simulations show that such associations are statistically significant, with probabilities of chance association of 3.5% and 10.7% for ID 51207 and ID 51380, respectively. Modeling of their optical/IR SEDs indicates that the CO detected galaxies and candidates have stellar masses and SFRs in the range (0.3-1.1)x10^11 M_sun and 60-160 M_sun/yr, with SFEs comparable to that found in other star-forming galaxies at similar redshifts. By comparing the space density of CO emitters derived from our observations with the space density derived from previous CO detections at z~1.5, and with semi-analytic predictions for the CO luminosity function, we suggest that the latter tend to underestimate the number of CO galaxies detected at high-redshift. Finally, we argue about the benefits of future blind CO searches in clustered fields with upcoming submm/radio facilities.

Massive AGN-driven outflows are invoked by AGN-galaxy co-evolutionary models to suppress both star formation and black hole accretion. Massive molecular outflows have recently been revealed in some AGN hosts. However, the physical properties and structure of these AGN-driven molecular outflows are still poorly constrained. Here we present new IRAM PdBI observations of Mrk231, the closest quasar known, targeting both the CO(1-0) and CO(2-1) transitions. We detect broad wings in both transitions, tracing a massive molecular outflow with velocities up to 800 km/s. The wings are spatially resolved at high significance level (5-11 sigma), indicating that the molecular outflow extends on the kpc scale. The CO(2-1)/CO(1-0) ratio of the red broad wings is consistent with the ratio observed in the narrow core, while the blue broad wing is less excited than the core. The latter result suggests that quasar driven outflow models invoking shocks (which would predict higher gas excitation) are not appropriate to describe the bulk of the outflow in Mrk231. However, we note that within the central 700 pc the CO(2-1)/CO(1-0) ratio of the red wing is slightly, but significantly, higher than in the line core, suggesting that shocks may play a role in the central region. We also find that the average size of the outflow anticorrelates with the critical density of the transition used as a wind tracer. This indicates that, although diffuse and dense clumps coexist in the outflowing gas, dense outflowing clouds have shorter lifetime and that they evaporate into the diffuse component along the outflow or, more simply, that diffuse clouds are more efficiently accelerated to larger distances by radiation pressure.

The analysis of the physical properties of low-redshift Ly$\alpha$ emitters (LAEs) can provide clues in the study of their high-redshift analogues. At $z \sim 0.3$, LAEs are bright enough to be detected over almost the entire electromagnetic spectrum and it is possible to carry out a more precise and complete study than at higher redshifts. In this study, we examine the UV and IR emission, dust attenuation, SFR and morphology of a sample of 23 GALEX-discovered star-forming (SF) LAEs at $z \sim 0.3$ with direct UV (GALEX), optical (ACS) and FIR (PACS and MIPS) data. Using the same UV and IR limiting luminosities, we find that LAEs at $z\sim 0.3$ tend to be less dusty, have slightly higher total SFRs, have bluer UV continuum slopes, and are much smaller than other galaxies that do not exhibit Ly$\alpha$ emission in their spectrum (non-LAEs). These results suggest that at $z \sim 0.3$ Ly$\alpha$ photons tend to escape from small galaxies with low dust attenuation. Regarding their morphology, LAEs belong to Irr/merger classes, unlike non-LAEs. Size and morphology represent the most noticeable difference between LAEs and non-LAEs at $z \sim 0.3$. Furthermore, the comparison of our results with those obtained at higher redshifts indicates that either the Ly$\alpha$ technique picks up different kind of galaxies at different redshifts or that the physical properties of LAEs are evolving with redshift.

We compare the average star formation (SF) activity in X-ray selected AGN hosts with mass-matched control inactive galaxies,including star forming and quiescent sources, at 0.5<z<2.5. Recent observations carried out by PACS, the 60-210um Herschel photometric camera, in GOODS-S, GOODS-N and COSMOS allow us to unbiasedly estimate the far-IR luminosity, and hence the SF properties, of the two samples. Accurate AGN host stellar masses are measured by decomposing their total emission into the stellar and nuclear components. We find a higher average SF activity in AGN hosts with respect to non-AGNs. The level of SF enhancement is modest (~0.26dex at ~3sigma) at low X-ray luminosities (Lx<~10^43.5erg/s) and more pronounced (0.56dex at >10sigma) for bright AGNs. However, when comparing to star forming galaxies only, AGN hosts are broadly consistent with the locus of their `main sequence'. We investigate the relative far-IR luminosity distributions of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host, while luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and SF, possibly through major mergers. While an increased SF with respect to non-AGNs of similar mass is expected in the latter, we interpret the modest SF offsets measured in low-Lx AGN hosts as either a) generated by non-synchronous accretion and SF histories in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale (non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario.

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.

In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Ly$\alpha$ emitting galaxies at 2.0$\lesssim$z$\lesssim$3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160$\mu$m observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard \emph{Herschel Space Observatory}. We define LAEs as those Ly$\alpha$ emitting galaxies whose rest-frame Ly$\alpha$ equivalent widths (Ly$\alpha$ EW$_{rest-frame}$) are above 20\AA, the typical threshold in narrow-band searches. Ly$\alpha$ emitting galaxies with Ly$\alpha$ EW$_{rest-frame}$ are called non-LAEs. As a result of an individual SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below 100Myr and a wide variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0$\lesssim$z$\lesssim$3.5 are very diverse, and do not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, with A$_{1200}$$\gtrsim$4mag. This is an indication from direct FIR measurements that dust and Ly$\alpha$ emission are not mutually exclusive. This population of red and dusty LAEs is not seen at z$\sim$0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via their Ly$\alpha$ emission.

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

We combine Herschel-PACS data from the PEP program with Spitzer 24 um and 16 um photometry and ultra deep IRS mid-infrared spectra, to measure the mid- to far-infrared spectral energy distribution (SED) of 0.7<z<2.5 normal star forming galaxies around the main sequence (the redshift-dependent relation of star formation rate and stellar mass). Our deep data confirm from individual far-infrared detections that z~2 star formation rates are overestimated if based on 24 um fluxes and SED templates that are calibrated via local trends with luminosity. Galaxies with similar ratios of rest-frame nuLnu(8) to 8-1000 um infrared luminosity (LIR) tend to lie along lines of constant offset from the main sequence. We explore the relation between SED shape and offset in specific star formation rate (SSFR) from the redshift-dependent main sequence. Main sequence galaxies tend to have a similar nuLnu(8)/LIR regardless of LIR and redshift, up to z~2.5, and nuLnu(8)/LIR decreases with increasing offset above the main sequence in a consistent way at the studied redshifts. We provide a redshift-independent calibration of SED templates in the range of 8--60 um as a function of log(SSFR) offset from the main sequence. Redshift dependency enters only through the evolution of the main sequence with time. Ultra deep IRS spectra match these SED trends well and verify that they are mostly due to a change in ratio of PAH to LIR rather than continua of hidden AGN. Alternatively, we discuss the dependence of nuLnu(8)/LIR on LIR. Same nuLnu(8)/LIR is reached at increasingly higher LIR at higher redshift, with shifts relative to local by 0.5 and 0.8 dex in log(LIR) at redshifts z~1 and z~2. Corresponding SED template calibrations are provided for use if no stellar masses are in hand. For most of those z~2 star forming galaxies that also host an AGN, the mid-infrared is dominated by the star forming component.

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.

Star formation in massive galaxies is quenched at some point during hierarchical mass assembly. To understand where and when the quenching processes takes place, we study the evolution of the total star formation rate per unit total halo mass (\Sigma(SFR/M)) in three different mass scales: low mass halos (field galaxies), groups, and clusters, up to a redshift ~1.6. We use deep far-infrared PACS data at 100 and 160 um to accurately estimate the total star formation rate of the Luminous Infrared Galaxy population of 9 clusters with mass ~10^{15} M_{\odot}, and 9 groups/poor clusters with mass ~ 5 x 10^{13} M_{\odot}. Estimates of the field \Sigma(SFR/M) are derived from the literature, by dividing the star formation rate density by the mean comoving matter density of the universe. The field \Sigma(SFR/M) increases with redshift up to z~1 and it is constant thereafter. The evolution of the \Sigma(SFR/M)-z relation in galaxy systems is much faster than in the field. Up to redshift z~0.2, the field has a higher \Sigma(SFR/M) than galaxy groups and galaxy clusters. At higher redshifts, galaxy groups and the field have similar \Sigma(SFR/M), while massive clusters have significantly lower \Sigma(SFR/M) than both groups and the field. There is a hint of a reversal of the SFR activity vs. environment at z~1.6, where the group \Sigma(SFR/M) lies above the field \Sigma(SFR/M)-z relation. We discuss possible interpretations of our results in terms of the processes of downsizing, and star-formation quenching.

We analyze the dependence of galaxy structure (size and Sersic index) and mode of star formation (\Sigma_SFR and SFR_IR/SFR_UV) on the position of galaxies in the SFR versus Mass diagram. Our sample comprises roughly 640000 galaxies at z~0.1, 130000 galaxies at z~1, and 36000 galaxies at z~2. Structural measurements for all but the z~0.1 galaxies were based on HST imaging, and SFRs are derived using a Herschel-calibrated ladder of SFR indicators. We find that a correlation between the structure and stellar population of galaxies (i.e., a 'Hubble sequence') is already in place since at least z~2.5. At all epochs, typical star-forming galaxies on the main sequence are well approximated by exponential disks, while the profiles of quiescent galaxies are better described by de Vaucouleurs profiles. In the upper envelope of the main sequence, the relation between the SFR and Sersic index reverses, suggesting a rapid build-up of the central mass concentration in these starbursting outliers. We observe quiescent, moderately and highly star-forming systems to co-exist over an order of magnitude or more in stellar mass. At each mass and redshift, galaxies on the main sequence have the largest size. The rate of size growth correlates with specific SFR, and so does \Sigma_SFR at each redshift. A simple model using an empirically determined SF law and metallicity scaling, in combination with an assumed geometry for dust and stars is able to relate the observed \Sigma_SFR and SFR_IR/SFR_UV, provided a more patchy dust geometry is assumed for high-redshift galaxies.

Understanding the mechanisms that lead dense environments to host galaxies with redder colors, more spheroidal morphologies, and lower star formation rates than field populations remains an important problem. As most candidate processes ultimately depend on host halo mass, accurate characterizations of the local environment, ideally tied to halo mass estimates and spanning a range in halo mass and redshift are needed. In this work, we present and test a rigorous, probabalistic method for assigning galaxies to groups based on precise photometric redshifts and X-ray selected groups drawn from the COSMOS field. The groups have masses in the range 10^13 < M_200c/M_sun < 10^14 and span redshifts 0<z<1. We characterize our selection algorithm via tests on spectroscopic subsamples, including new data obtained at the VLT, and by applying our method to detailed mock catalogs. We find that our group member galaxy sample has a purity of 84% and completeness of 92% within 0.5 R200c. We measure the impact of uncertainties in redshifts and group centering on the quality of the member selection with simulations based on current data as well as future imaging and spectroscopic surveys. As a first application of our new group member catalog which will be made publicly available, we show that member galaxies exhibit a higher quenched fraction compared to the field at fixed stellar mass out to z~1, indicating a significant relationship between star formation and environment at group scales. We also address the suggestion that dusty star forming galaxies in such groups may impact the high-l power spectrum of the cosmic microwave background and find that such a population cannot explain the low power seen in recent SZ measurements.

We explore the connection between black hole growth at the center of obscured quasars selected from the XMM-COSMOS survey and the physical properties of their host galaxies. We study a bolometric regime (<Lbol > 8 x 10^45 erg/s) where several theoretical models invoke major galaxy mergers as the main fueling channel for black hole accretion. We confirm that obscured quasars mainly reside in massive galaxies (Mstar>10^10 Msun) and that the fraction of galaxies hosting such powerful quasars monotonically increases with the stellar mass. We stress the limitation of the use of rest-frame color-magnitude diagrams as a diagnostic tool for studying galaxy evolution and inferring the influence that AGN activity can have on such a process. We instead use the correlation between star-formation rate and stellar mass found for star-forming galaxies to discuss the physical properties of the hosts. We find that at z ~1, ~62% of Type-2 QSOs hosts are actively forming stars and that their rates are comparable to those measured for normal star-forming galaxies. The fraction of star-forming hosts increases with redshift: ~71% at z ~2, and 100% at z ~3. We also find that the the evolution from z ~1 to z ~3 of the specific star-formation rate of the Type-2 QSO hosts is in excellent agreement with that measured for star-forming galaxies. From the morphological analysis, we conclude that most of the objects are bulge-dominated galaxies, and that only a few of them exhibit signs of recent mergers or disks. Finally, bulge-dominated galaxies tend to host Type-2 QSOs with low Eddington ratios (lambda<0.1), while disk-dominated or merging galaxies have at their centers BHs accreting at high Eddington ratios (lambda > 0.1).

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

We study the impact of galaxy-galaxy interactions on the FIR properties of galaxies and its evolution at 0<z<1.2. Using the high-z galaxies in the fields of GOODS observed by Herschel in the framework of the GOODS-Herschel key program and the local IRAS or AKARI-selected galaxies in the field of SDSS DR7, we investigate the dependence of galaxy properties on the morphology of and the distance to the nearest neighbor galaxy. We find that the SFRs and the SSFRs of galaxies, on average, depend on the morphology of and the distance to the nearest neighbor galaxy in this redshift range. When a late-type galaxy has a close neighbor galaxy, the SFR and the SSFR increase as it approaches a late-type neighbor, which is supported by K-S and MC tests with a significance level of >99%. However, the SFR and the SSFR decrease or do not change much as it approaches an early-type neighbor. The bifurcations of SFRs and SSFRs depending on the neighbor's morphology seem to occur at R_n~0.5r_vir,nei, which is supported by K-S and MC tests with a significance level of >98%. For all redshift bins, the SSFRs of late-type galaxies interacting with late-type neighbors are increased by factors of about 1.8+/-0.7 and 4.0+/-1.2 compared to those of non-interacting galaxies when the pair separation is smaller than 0.5r_vir,nei and 0.1r_vir,nei, respectively. The Tdust of both local and high-z late-type galaxies strongly interacting with late-type neighbors (i.e. R_n < 0.1r_vir,nei) appears to be higher than that of non-interacting galaxies with a significance level of 96-99%. However, Tdust of local late-type galaxies strongly interacting with early-type neighbors seems to be lower than or similar to that of non-interacting galaxies. Our results suggest that galaxy-galaxy interactions and mergers have been strongly affecting the SFR and the dust properties of star-forming galaxies over at least 8 billion years.

We present a study of the infrared properties of X-ray selected, moderate luminosity (Lx=10^{42}-10^{44}ergs/s) active galactic nuclei (AGNs) up to z~3, to explore the links between star formation in galaxies and accretion onto their central black holes. We use 100um and 160um fluxes from GOODS-Herschel -the deepest survey yet undertaken by the Herschel telescope- and show that in >94 per cent of cases these fluxes are dominated by the host. We find no evidence of any correlation between the X-ray and infrared luminosities of moderate AGNs at any redshift, suggesting that star-formation is decoupled from nuclear (AGN) activity. The star formation rates of AGN hosts increase strongly with redshift; by a factor of 43 from z<0.1 to z=2-3 for AGNs with the same X-ray luminosities. This increase is consistent with the factor of 25-50 increase in the specific star formation rates (SSFRs) of normal, star-forming (main-sequence) galaxies. Indeed, the average SSFRs of AGN hosts are only marginally (20 per cent) lower than those of main-sequence galaxies, with this small deficit being due to a fraction of AGNs residing in quiescent (low-SSFR) galaxies. We estimate 79+/-10 per cent of moderate AGNs are hosted in main-sequence galaxies, 15+/-7 per cent in quiescent galaxies and <10 per cent in strongly starbursting galaxies. The fractions of all main sequence galaxies at z<2 experiencing a period of moderate nuclear activity is strongly dependent on galaxy stellar mass (Mstars); rising from a few per cent at Mstars~10^{10}Msun to >20 per cent at Mstars>10^{11}Msun. Our results indicate that it is galaxy stellar mass that is most important in dictating whether a galaxy hosts a moderate luminosity AGN. We argue that the majority of moderate nuclear activity is fuelled by internal mechanisms rather than violent mergers, suggesting that disk instabilities could be an important AGN feeding mechanism.

We explore the gas-to-dust mass ratio (G/D) and the CO luminosity-to-Mgas conversion factor (a_co) of two well studied galaxies in the GOODS-N field, that are expected to have different star forming modes, the starburst GN20 at z=4.05 and the normal star-forming galaxy BzK-21000 at z=1.52. Detailed sampling is available for their Rayleigh-Jeans emission via ground based mm interferometry (1.1-6.6mm) along with Herschel, PACS and SPIRE data that probe the peak of their infrared emission. Using the physically motivated Draine & Li (2007) models, as well as a modified black body function, we measure the dust mass (Md) of the sources and find 2.0^{+0.7}_{-0.6} x 10^{9} Msun for GN20 and 8.6^{+0.6}_{-0.9} x 10^{8} Msun for BzK-21000. The addition of mm data reduces the uncertainties of the derived Md by a factor of ~2, allowing the use of the local G\D vs metallicity relation to place constraints on the a_co values of the two sources. For GN20 we derive a conversion factor of a_co < 1.0 Msun pc^{-2}(K km s^{-1})^{-1}, consistent with that of local ULIRGs, while for BzK-21000 we find a considerably higher value, a_co ~4.0 Msun pc^{-2}(K km s^{-1})^{-1}, in agreement with an independent kinematic derivation reported previously. The implied star formation efficiency is ~25 Lsun/Msun for BzK-21000, a factor of ~5-10 lower than that of GN20. The findings for these two sources support the existence of different disk-like and starburst star-formation modes in distant galaxies, although a larger sample is required to draw statistically robust results

With this paper, we release accurate photometric redshifts for 1692 counterparts to Chandra sources in the central square degree of the COSMOS field. The availability of a large training set of spectroscopic redshifts that extends to faint magnitudes enabled photometric redshifts comparable to the highest quality results presently available for normal galaxies. We demonstrate that morphologically extended, faint X-ray sources without optical variability are more accurately described by a library of normal galaxies (corrected for emission lines) than by AGN-dominated templates, even if these sources have AGN-like X-ray luminosities. Preselecting the library on the bases of the source properties allowed us to reach an accuracy sigma_(Delta z/(1+z_spec)) \sim0.015 with a fraction of outliers of 5.8% for the entire Chandra-COSMOS sample. In addition, we release revised photometric redshifts for the 1735 optical counterparts of the XMM-detected sources over the entire 2 sq. deg.of COSMOS. For 248 sources, our updated photometric redshift differs from the previous release by Delta z>0.2. These changes are predominantly due to the inclusion of newly available deep H-band photometry H_AB=24 mag. We illustrate once again the importance of a spectroscopic training sample and how an assumption about the nature of a source together with the number and the depth of the available bands influence the accuracy of the photometric redshifts determined for AGN. These considerations should be kept in mind when defining the observational strategies of upcoming large surveys targeting AGN, such as eROSITA at X-ray energies and ASKAP/EMU in the radio band.

One remaining open question regarding the physical properties of Ly$\alpha$ emitters (LAEs) is their dust content and its evolution with redshift. The variety of results is large and with those reported by now is difficult to establish clear relations between dust, other fundamental parameters of galaxies (star-formation rate, metallicity or age) and redshift. In this Letter, we report \emph{Herschel} PACS-100$\mu$m, PACS-160$\mu$m and \emph{Spitzer} MIPS-24$\mu$m detections of a sample of spectroscopically GALEX-selected LAEs at z$\sim$0.3 and $\sim$1.0. Five out of ten and one out of two LAEs are detected in, at least, one PACS band at z$\sim$0.3 and $\sim$1.0, respectively. These measurements have a great importance given that they allow us to quantify, for the first time, the dust content of LAEs from direct FIR observations. MIPS-24$\mu$m detections allow us to determine IR properties of the PACS-undetected LAEs. We obtain that mid-IR/FIR detected star-forming (SF) LAEs at z$\sim$0.3 have dust content within 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.0, with a median value of A$_{1200\textrm{\AA}}$$\sim$1.1. This range broadens out to 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.5 when considering those LAEs at z$\sim$1.0. Only one SF LAE is undetected both in MIPS-24$\mu$m and PACS, with $A_{1200\AA}$ $\lesssim$0.75. These results seem to be larger than those reported for high-redshift LAEs and, therefore, although an evolutionary trend is not clearly seen, it could point out that low-redshift LAEs are dustier than high-redshift ones. However, the diverse methods used could introduce a systematic offset in the results.

Based on broad/narrow-band photometry and Keck DEIMOS spectroscopy we report a redshift of z=4.64-0.08+0.06 for AzTEC/COSMOS 1, the brightest sub-mm galaxy in the AzTEC/COSMOS field. In addition to the COSMOS-survey X-ray to radio data, we report observations of the source with Herschel/PACS (100, 160 micron), CSO/SHARC II (350 micron), CARMA and PdBI (3 mm). We do not detect CO(5-4) line emission in the covered redshift ranges, 4.56-4.76 (PdBI/CARMA) and 4.94-5.02 (CARMA). If the line is within this bandwidth, this sets 3sigma upper limits on the gas mass to <~8x10^9 M_Sol and <~5x10^10 M_Sol, respectively (assuming similar conditions as observed in z~2 SMGs). This could be explained by a low CO-excitation in the source. Our analysis of the UV-IR spectral energy distribution of AzTEC 1 shows that it is an extremely young (<~50 Myr), massive (M*~10^11 M_Sol), but compact (<~2 kpc) galaxy forming stars at a rate of ~1300 M_Sol/yr. Our results imply that AzTEC 1 is forming stars in a 'gravitationally bound' regime in which gravity prohibits the formation of a superwind, leading to matter accumulation within the galaxy and further generations of star formation.

Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.

Two main modes of star formation are know to control the growth of galaxies: a relatively steady one in disk-like galaxies, defining a tight star formation rate (SFR)-stellar mass sequence, and a starburst mode in outliers to such a sequence which is generally interpreted as driven by merging. Such starburst galaxies are rare but have much higher SFRs, and it is of interest to establish the relative importance of these two modes. PACS/Herschel observations over the whole COSMOS and GOODS-South fields, in conjunction with previous optical/near-IR data, have allowed us to accurately quantify for the first time the relative contribution of the two modes to the global SFR density in the redshift interval 1.5<z<2.5, i.e., at the cosmic peak of the star formation activity. The logarithmic distributions of galaxy SFRs at fixed stellar mass are well described by Gaussians, with starburst galaxies representing only a relatively minor deviation that becomes apparent for SFRs more than 4 times higher than on the main sequence. Such starburst galaxies represent only 2% of mass-selected star forming galaxies and account for only 10% of the cosmic SFR density at z~2. Only when limited to SFR>1000M(sun)/yr, off-sequence sources significantly contribute to the SFR density (46+/-20%). We conclude that merger-driven starbursts play a relatively minor role for the formation of stars in galaxies, whereas they may represent a critical phase towards the quenching of star formation and morphological transformation in galaxies.

We take advantage of the sensitivity and resolution of Herschel at 100 and 160 micron to directly image the thermal dust emission and investigate the infrared luminosities, L(IR), and dust obscuration of typical star-forming (L*) galaxies at high redshift. Our sample consists of 146 UV-selected galaxies with spectroscopic redshifts 1.5<z<2.6 in the GOODS-North field. Supplemented with deep Very Large Array (VLA) and Spitzer imaging, we construct median stacks at the positions of these galaxies at 24, 100, and 160 micron, and 1.4 GHz. The comparison between these stacked fluxes and a variety of dust templates and calibrations implies that typical star-forming galaxies with UV luminosities L(UV)>1e10 Lsun at z~2 are luminous infrared galaxies (LIRGs) with a median L(IR)=(2.2+/-0.3)e11 Lsun. Typical galaxies at 1.5<z<2.6 have a median dust obscuration L(IR)/L(UV) = 7.1+/-1.1, which corresponds to a dust correction factor, required to recover the bolometric star formation rate (SFR) from the unobscured UV SFR, of 5.2+/-0.6. This result is similar to that inferred from previous investigations of the UV, H-alpha, 24 micron, radio, and X-ray properties of the same galaxies studied here. Stacking in bins of UV slope implies that L* galaxies with redder spectral slopes are also dustier, and that the correlation between UV slope and dustiness is similar to that found for local starburst galaxies. Hence, the rest-frame 30 and 50 micron fluxes validate on average the use of the local UV attenuation curve to recover the dust attenuation of typical star-forming galaxies at high redshift. In the simplest interpretation, the agreement between the local and high redshift UV attenuation curves suggests a similarity in the dust production and stellar and dust geometries of starburst galaxies over the last 10 billion years.

Dust attenuation curves in external galaxies are useful to study their dust properties as well as to interpret their intrinsic spectral energy distributions. In particular the presence or absence of a UV bump at 2175 A remains an open issue which has consequences on the interpretation of broad band colours of distant galaxies. We study the dust attenuation curve in the UV range at z >1. In particular we search for the presence of a UV bump. We use deep photometric data of the CDFS obtained with intermediate and broad band filters by the MUSYC project to sample the UV rest-frame of galaxies with 1<z <2. Herschel/PACS and Spitzer/MIPS data are used to measure the dust emission. 30 galaxies were selected with high S/N in all bands. Their SEDs from the UV to the far-IR are fitted using the CIGALE code and the characteristics of the dust attenuation curve are obtained. The mean dust attenuation curve we derive exhibits a UV bump at 2175A whose amplitude corresponds to 35 % (76%) that of the MW (LMC2 supershell) extinction curve. An analytical expression of the average attenuation curve is given, it is found slightly steeper than the Calzetti et al. one, although at a 1 sigma level. Our galaxy sample is used to study the derivation of the slopes of the UV continuum from broad band colours, including the GALEX FUV-NUV colour. Systematic errors induced by the presence of the bump are quantified. We compare dust attenuation factors measured with CIGALE to the slope of the UV continuum and find that there is a large scatter around the relation valid for local starbursts (0.7 mag). The uncertainties on the determination of the UV slope lead to an extra systematic error of the order of 0.3 to 0.7 mag on dust attenuation when a filter overlaps the UV bump.

We compare multi-wavelength SFR indicators out to z~3 in GOODS-South. Our analysis uniquely combines U-to-8um photometry from FIREWORKS, MIPS 24um and PACS 70, 100, and 160um photometry from the PEP survey, and Ha spectroscopy from the SINS survey. We describe a set of conversions that lead to a continuity across SFR indicators. A luminosity-independent conversion from 24um to total infrared luminosity yields estimates of LIR that are in the median consistent with the LIR derived from PACS photometry, albeit with significant scatter. Dust correction methods perform well at low to intermediate levels of star formation. They fail to recover the total amount of star formation in systems with large SFR_IR/SFR_UV ratios, typically occuring at the highest SFRs (SFR_UV+IR \gtrsim 100 Msun/yr) and redshifts (z \gtrsim 2.5) probed. Finally, we confirm that Ha-based SFRs at 1.5<z<2.6 are consistent with SFR_SED and SFR_UV+IR provided extra attenuation towards HII regions is taken into account (Av,neb = Av,continuum / 0.44). With the cross-calibrated SFR indicators in hand, we perform a consistency check on the star formation histories inferred from SED modeling. We compare the observed SFR-M relations and mass functions at a range of redshifts to equivalents that are computed by evolving lower redshift galaxies backwards in time. We find evidence for underestimated stellar ages when no stringent constraints on formation epoch are applied. We demonstrate how resolved SED modeling, or alternatively deep UV data, may help to overcome this bias. The age bias is most severe for galaxies with young stellar populations, and reduces towards older systems. Finally, our analysis suggests that SFHs typically vary on timescales that are long (at least several 100 Myr) compared to the galaxies' dynamical time.

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program and its role in the complement of Herschel surveys, and the field selection which includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, EGS. We provide an account of the observing strategies and data reduction methods used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high redshift galaxy populations, thus testing and superseeding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.

The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 micron, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ~1 mJy to few hundreds mJy. Stacking of 24 micron sources and P(D) statistics extend the analysis down to ~0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts and CIB are studied up to z=5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. Employing the P(D) analysis, these fractions increase to ~65% and ~89%. More than half of the resolved CIB was emitted at redshift z<=1. The 50%-light redshifts lie at z=0.58, 0.67 and 0.73 at the three PACS wavelengths. The distribution moves towards earlier epochs at longer wavelengths: while the 70 micron CIB is mainly produced by z<=1.0 objects, the contribution of z>1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun).

A substantial fraction of the stellar mass growth across cosmic time occurred within dust-enshrouded environments. Yet, the exact amount of star-forming activity that took place in high-redshift dusty galaxies currently missed by optical surveys has been barely explored. Using the Spitzer observations of COSMOS we determined the fraction of luminous star-forming galaxies at 1.5<z<3 escaping the traditional color selection techniques because of dust extinction, as well as their contribution to the cosmic star formation density at high redshift. We find that the BzK criterion offers an almost complete (~90%) identification of the 24mic sources at 1.4<z<2.5, while the BM/BX criterion miss 50% of the MIPS population. Similarly the criterion based on the presence of a stellar bump in massive sources (so-called "IRAC peakers") miss up to 40% of the IR luminosity density and only 25% of the IR energy density at z~2 is produced by Optically-Faint IR-bright galaxies selected based on their extreme mid-IR to optical flux ratios. We conclude that color selections of distant star-forming galaxies must be used with lots of care given the substantial bias they can suffer. In particular, the effect of dust extinction strongly impacts the completeness of identifications at the bright end of the bolometric luminosity function, which implies large and uncertain extrapolations to account for the contribution of dusty galaxies missed by these selections. In the context of forthcoming facilities that will operate at long wavelengths (e.g., $JWST$, ALMA, SAFARI, EVLA, SKA), this emphasizes the importance of minimizing the extinction biases when probing the activity of star formation in the early Universe.

ABRIDGED-This paper presents the first direct estimate of the 3D clustering properties of far-infrared sources up to z~3. This has been possible thanks to the Pacs Evolutionary Probe (PEP) survey of the GOODS South field performed with the PACS instrument onboard the Herschel Satellite. An analysis of the two-point correlation function over the whole redshift range spanned by the data reports for the correlation length, r_0~6.3 Mpc and r_0~6.7 Mpc, respectively at 100um and 160um, corresponding to dark matter halo masses M>~10^{12.4} M_sun. Objects at z~2 instead seem to be more strongly clustered, with r_0~19 Mpc and r_0~17 Mpc in the two considered PACS channels. This dramatic increase of the correlation length between z~1 and z~2 is connected with the presence of a wide, M>~10^{14} M_sun, filamentary structure which includes more than 50% of the sources detected at z~2. An investigation of the properties of such sources indicates the possibility for boosted star-forming activity in those which reside within the overdense environment with respect of more isolated galaxies found in the same redshift range. Lastly, we also present our results on the evolution of the relationship between luminous and dark matter in star-forming galaxies between z~1 and z~2. We find that the increase of (average) stellar mass in galaxies <M*> between z~1 and z~2 is about a factor 10 lower than that of the dark matter haloes hosting such objects (<M*>[z~1]/<M*>[z~2] ~ 0.4 vs M_{halo}[z~1]/M_{halo}[z~2] ~ 0.04). Our findings agree with the evolutionary picture of downsizing whereby massive galaxies at z~2 were more actively forming stars than their z~1 counterparts, while at the same time contained a lower fraction of their mass in the form of luminous matter.

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 have identified a complete, flux-limited, (S_160>120 mJy), sample of 160 micron-selected sources from Spitzer observations of the 1-deg^2 ISO Deep Field region in the Lockman Hole. Ground-based UV, optical and near-infrared (NIR) photometry and optical spectroscopy have been used to determine colors, redshifts and masses for the complete sample of 40 galaxies. Spitzer-IRAC+MIPS photometry, supplemented by ISOPHOT data at 90 micron and 170 micron, has been used to calculate accurate total infrared luminosities, LIR(8-1000 micron), and to determine the IR luminosity function (LF) of luminous infrared galaxies (LIRGs). The maximum observed redshift is z~0.80 and the maximum total infrared luminosity is log(L_IR/L_Sun)=12.74. Over the luminosity range log(L_IR/L_Sun)=10-12, the LF for LIRGs in the Lockman Hole Deep Field is similar to that found previously for local sources at similar infrared luminosities. The mean host galaxy mass, log(M/M_Sun)=10.7, and dominance of HII-region spectral types, is also similar to what has been found for local LIRGs, suggesting that intense starbursts likely power the bulk of the infrared luminosity for sources in this range of LIR. However for the most luminous sources, log(L_IR/L_Sun)>12.0, we find evidence for strong evolution in the LF \propto (1+z)6\pm1, assuming pure number density evolution. These ultraluminous infrared galaxies (ULIRGs) have a larger mean host mass, log(M/M_Sun)=11.0, and exhibit disturbed morphologies consistent with strong-interactions/mergers, and they are also more likely to be characterized by starburst-AGN composite or AGN spectral types.

We present three bright z+ dropout candidates selected from deep Near-Infrared (NIR) imaging of the COSMOS 2 square degree field. All three objects match the 0.8-8um colors of other published z>7 candidates but are three magnitudes brighter, facilitating further study. Deep spectroscopy of two of the candidates covering 0.64-1.02um with Keck-DEIMOS and all three covering 0.94-1.10um and 1.52-1.80um with Keck-NIRSPEC detects weak spectral features tentatively identified as Ly-alpha at z=6.95 and z=7.69 in two of the objects. The third object is placed at z~1.6 based on a 24um and weak optical detection. A comparison with the spectral energy distributions of known z<7 galaxies, including objects with strong spectral lines, large extinction, and large systematic uncertainties in the photometry yields no objects with similar colors. However, the lambda>1um properties of all three objects can be matched to optically detected sources with photometric redshifts at z~1.8, so the non-detection in the i+ and z+ bands are the primary factors which favors a z>7 solution. If any of these objects are at z~7 the bright end of the luminosity function is significantly higher at z>7 than suggested by previous studies, but consistent within the statistical uncertainty and the dark matter halo distribution. If these objects are at low redshift, the Lyman-Break selection must be contaminated by a previously unknown population of low redshift objects with very strong breaks in their broad band spectral energy distributions and blue NIR colors. The implications of this result on luminosity function evolution at high redshift is discussed. We show that the primary limitation of z>7 galaxy searches with broad filters is the depth of the available optical data.

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.

XMM and Chandra opened a new area for the study of clusters of galaxies. Not only for cluster physics but also, for the detection of faint and distant clusters that were inaccessible with previous missions. This article presents 66 spectroscopically confirmed clusters (0.05<z<1.5) within an area of 6 deg2 enclosed in the XMM-LSS survey. Almost two thirds have been confirmed with dedicated spectroscopy only and 10% have been confirmed with dedicated spectroscopy supplemented by literature redshifts. Sub-samples, or classes, of extended-sources are defined in a two-dimensional X-ray parameter space allowing for various degrees of completeness and contamination. We describe the procedure developed to assess the reality of these cluster candidates using the CFHTLS photometric data and spectroscopic information from our own follow-up campaigns. Most of these objects are low mass clusters, hence constituting a still poorly studied population. In a second step, we quantify correlations between the optical properties such as richness or velocity dispersion and the cluster X-ray luminosities. We examine the relation of the clusters to the cosmic web. Finally, we review peculiar structures in the surveyed area like very distant clusters and fossil groups.

We study the impact of the environment on the evolution of galaxies in the zCOSMOS 10k sample in the redshift range 0.1<z<1.0 over an area of ~1.5 deg2. The considered sample of secure spectroscopic redshifts contains about 8500 galaxies, with their stellar masses estimated by SED fitting of the multiwavelength optical to NIR photometry. The evolution of the galaxy stellar mass function (GSMF) in high and low density regions provides a tool to study the mass assembly evolution in different environments; moreover, the contributions to the GSMF from different galaxy types, as defined by their SEDs and their morphologies, can be quantified. At redshift z~1, the GSMF is only slightly dependent on environment, but at lower redshifts the shapes of the GSMFs in high- and low-density environments become extremely different, with high density regions exhibiting a marked bimodality. As a result, we infer that galaxy evolution depends on both the stellar mass and the environment, the latter setting the probability of a galaxy to have a given mass: all the galaxy properties related to the stellar mass show a dependence on environment, reflecting the difference observed in the mass functions. The shapes of the GSMFs of early- and late-type galaxies are almost identical for the extremes of the density contrast we consider. The evolution toward z=0 of the mass at which the early- and late-type GSMFs match each other is more rapid in high density environments. The comparison of the observed GSMFs to the same quantities derived from a set of mock catalogues shows that blue galaxies in sparse environments are overproduced in the semi-analytical models at intermediate and high masses, because of a deficit of star formation suppression, while at z<0.5 an excess of red galaxies is present in dense environments at intermediate and low masses, because of the overquenching of satellites. ABRIDGED

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.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

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 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.

The Cosmic Far-Infrared Background (CIB) at wavelengths around 160 {\mu}m corresponds to the peak intensity of the whole Extragalactic Background Light, which is being measured with increasing accuracy. However, the build up of the CIB emission as a function of redshift, is still not well known. Our goal is to measure the CIB history at 70 {\mu}m and 160 {\mu}m at different redshifts, and provide constraints for infrared galaxy evolution models. We use complete deep Spitzer 24 {\mu}m catalogs down to about 80 {\mu}Jy, with spectroscopic and photometric redshifts identifications, from the GOODS and COSMOS deep infrared surveys covering 2 square degrees total. After cleaning the Spitzer/MIPS 70 {\mu}m and 160 {\mu}m maps from detected sources, we stacked the far-IR images at the positions of the 24 {\mu}m sources in different redshift bins. We measured the contribution of each stacked source to the total 70 and 160 {\mu}m light, and compare with model predictions and recent far-IR measurements made with Herschel/PACS on smaller fields. We have detected components of the 70 and 160 {\mu}m backgrounds in different redshift bins up to z ~ 2. The contribution to the CIB is maximum at 0.3 <= z <= 0.9 at 160{\mu}m (and z <= 0.5 at 70 {\mu}m). A total of 81% (74%) of the 70 (160) {\mu}m background was emitted at z < 1. We estimate that the AGN relative contribution to the far-IR CIB is less than about 10% at z < 1.5. We provide a comprehensive view of the CIB buildup at 24, 70, 100, 160 {\mu}m. IR galaxy models predicting a major contribution to the CIB at z < 1 are in agreement with our measurements, while our results discard other models that predict a peak of the background at higher redshifts. Our results are available online http://www.ias.u-psud.fr/irgalaxies/ .