We present the first spatially resolved study of molecular gas in the vicinity of a Gamma Ray Burst, using CO(2-1) emission line observations with the Atacama Large Millimetre Array (ALMA) at ~50 pc scales. The host galaxy of GRB 980425 contains a ring of high column density HI gas which is likely to have formed due to a collision between the GRB host and its companion galaxy, within which the GRB is located. We detect eleven molecular gas clumps in the galaxy, seven of which are within the gas ring. The clump closest to the GRB position is at a projected separation of ~280 pc. Although it is plausible that the GRB progenitor was ejected from clusters formed in this clump, we argue that the in situ formation of the GRB progenitor is the most likely scenario. We measure the molecular gas masses of the clumps and find them to be sufficient for forming massive star clusters. The molecular gas depletion times of the clumps show a variation of ~2 dex, comparable with the large variation in depletion times found in starburst galaxies in the nearby Universe. This demonstrates the presence of starburst modes of star formation on local scales in the galaxy, even while the galaxy as a whole cannot be categorised as a starburst based on its global properties. Our findings suggest that the progenitor of GRB 9802425 was originated in a young massive star cluster formed in the starburst mode of star formation.

We present the detection of multiple CO line transitions with ALMA in a few tens of infrared-selected galaxies on and above the main sequence at z=1.1-1.7. We reliably detected the emission of CO(5-4), CO(2-1), and CO(7-6)+[CI](2-1) in 50, 33, and 13 galaxies, respectively, and we complemented this information with available CO(4-3) and [CI](1-0) fluxes for part of the sample, and modeling of the optical-to-mm SEDs. We retrieve a quasi-linear relation between LIR and CO(5-4) or CO(7-6) for main-sequence galaxies and starbursts, corroborating the hypothesis that these transitions can be used as SFR tracers. We find the CO excitation to steadily increase as a function of the star formation efficiency, the mean intensity of the radiation field warming the dust, the surface density of SFR, and, less distinctly, with the distance from the main sequence. This adds to the tentative evidence for higher excitation of the CO+[CI] SLED of starbursts relative to that for main-sequence objects, where the dust opacities play a minor role in shaping the high-J CO transitions in our sample. However, the distinction between the average SLED of upper main-sequence and starburst galaxies is blurred, driven by a wide variety of intrinsic shapes. LVG radiative transfer modeling demonstrates the existence of a highly excited component that elevates the CO SLED of high-redshift main-sequence and starbursting galaxies above the typical values observed in the disk of the Milky Way. This excited component is dense and it encloses ~50% of the total molecular gas mass in main-sequence objects. We interpret the observed trends involving the CO excitation as mainly driven by a combination of large SFRs and compact sizes, naturally connected with enhanced dense molecular gas fractions and higher dust and gas temperatures, due to increasing UV radiation fields, cosmic ray rates, and dust/gas coupling. [Abridged]

We present rest-frame far-infrared (FIR) and optical size measurements of AGN hosts and star-forming galaxies in the COSMOS field, enabled by high-resolution ALMA/1 mm (0.1 arcsec - 0.4 arcsec) and HST/F814W imaging (~ 0.1 arcsec). Our sample includes 27 galaxies at z<2.5, classified as infrared-selected AGN (3 sources), X-ray selected AGN (4 sources), and non-AGN star-forming galaxies (20 sources), for which high-resolution Band 6/7 ALMA images are available at 1 mm from our own observing program as well as archival observations. The sizes and SFR surface densities measured from both ALMA/1 mm and HST/F814W images show that obscured AGN host galaxies are more compact than non-AGN star-forming galaxies at similar redshift and stellar mass. This result suggests that the obscured accretion phase may be related to galaxies experiencing a compaction of their gaseous component, which could be associated with enhanced central star formation before a subsequent quenching driving the formation of compact passive galaxies. Moreover, most of the detected and stacked rest-frame FIR sizes of AGNs in our sample are similar or more compact than their rest-frame optical sizes, which is consistent with recent results of ALMA detected sources. This might be explained by the fact that the dusty starbursts take place in the compact regions, and suggests that the star formation mechanisms in the compact regions of AGN hosts are similar to those observed in star-forming galaxies observed with ALMA.

We investigate the contribution of clumps and satellites to the galaxy mass assembly. We analyzed spatially-resolved Hubble Space Telescope observations (imaging and slitless spectroscopy) of 53 star-forming galaxies at z ~ 1 - 3. We created continuum and emission line maps and pinpointed residual "blobs" detected after subtracting the galaxy disk. Those were separated into compact (unresolved) and extended (resolved) components. Extended components have sizes ~ 2 kpc and comparable stellar mass and age as the galaxy disks, whereas the compact components are 1.5 dex less massive and 0.4 dex younger than the disks. Furthermore the extended blobs are typically found at larger distances from the galaxy barycenter than the compact ones. Prompted by these observations and by the comparison with simulations, we suggest that compact blobs are in-situ formed clumps, whereas the extended ones are accreting satellites. Clumps and satellites enclose respectively ~ 20% and ~< 80% of the galaxy stellar mass, ~ 30% and ~ 20% of its star formation rate. Considering the compact blobs, we statistically estimated that massive clumps (Mstar >~ 10^9 Msun) have lifetimes of ~ 650 Myr, and the less massive ones (10^8 < Mstar < 10^9 Msun) of ~ 145 Myr. This supports simulations predicting long-lived clumps (lifetime > 100 Myr). Finally, ~< 30% (13%) of our sample galaxies are undergoing single (multiple) merger(s), they have a projected separation ~< 10 kpc, and the typical mass ratio of our satellites is 1:5 (but ranges between 1:10 and 1:1), in agreement with literature results for close pair galaxies.

We investigate the fraction of close pairs and morphologically identified mergers on and above the star-forming main sequence (MS) at 0.2$\leq z\leq$2.0. The novelty of our work lies in the use of a non-parametric morphological classification performed on resolved stellar mass maps, reducing the contamination by non-interacting, high-redshift clumpy galaxies. We find that the merger fraction rapidly rises to $\geq$70% above the MS, implying that -- already at $z{\gtrsim}1$ -- starburst (SB) events ($\Delta_{\rm MS}\geq$0.6) are almost always associated with a major merger (1:1 to 1:6 mass ratio). The majority of interacting galaxies in the SB region are morphologically disturbed, late-stage mergers. Pair fractions show little dependence on MS-offset and pairs are more prevalent than late-stage mergers only in the lower half of the MS. In our sample, major mergers on the MS occur with a roughly equal frequency of $\sim$5-10% at all masses ${\gtrsim} 10^{10}M_{\odot}$. The MS major merger fraction roughly doubles between $z=0.2$ and $z=2$, with morphological mergers driving the overall increase at $z{\gtrsim}1$. The differential redshift evolution of interacting pairs and morphologically classified mergers on the MS can be reconciled by evolving observability timescales for both pairs and morphological disturbances. The observed variation of the late-stage merger fraction with $\Delta_{\rm MS}$ follows the perturbative 2-Star Formation Mode model, where any MS galaxy can experience a continuum of different SFR enhancements. This points to a starburst-merger connection not only for extreme events, but also more moderate bursts which merely scatter galaxies upward within the MS, rather than fully elevating them above it.

Unusual stellar explosions represent an opportunity to learn about both stellar and galaxy evolution. Mapping the atomic gas in host galaxies of such transients can lead to an understanding of the conditions triggering them. We provide resolved atomic gas observations of the host galaxy, CGCG137-068, of the unusual, poorly-understood transient AT2018cow searching for clues to understand its nature. We test whether it is consistent with a recent inflow of atomic gas from the intergalactic medium, as suggested for host galaxies of gamma-ray bursts (GRBs) and some supernovae (SNe). We observed the HI hyperfine structure line of the AT2018cow host with the Giant Metrewave Radio Telescope. There is no unusual atomic gas concentration near the position of AT2018cow. The gas distribution is much more regular than those of GRB/SN hosts. The AT2018cow host has an atomic gas mass lower by 0.24 dex than predicted from its star formation rate (SFR) and is at the lower edge of the galaxy main sequence. In the continuum we detected the emission of AT2018cow and of a star-forming region in the north-eastern part of the bar (away from AT2018cow). This region hosts a third of the galaxy's SFR. The absence of atomic gas concentration close to AT2018cow, along with a normal SFR and regular HI velocity field, sets CGCG137-068 apart from GRB/SN hosts studied in HI. The environment of AT2018cow therefore suggests that its progenitor may not have been a massive star. Our findings are consistent with an origin of the transient that does not require a connection between its progenitor and gas concentration or inflow: an exploding low-mass star, a tidal disruption event, a merger of white dwarfs, or a merger between a neutron star and a giant star. We interpret the recently reported atomic gas ring in CGCG137-068 as a result of internal processes connected with gravitational resonances caused by the bar.

We report the Atacama Large Millimeter/submillimeter Array (ALMA) observations of the metal rich host galaxy of superluminous supernova (SLSN) PTF10tpz, a barred spiral galaxy at z=0.03994. We find the CO(1-0) emission to be confined within the bar of the galaxy. The distribution and kinematics of molecular gas in the host galaxy resemble gas flows along two lanes running from the tips of the bar towards the galaxy center. These gas lanes end in a gaseous structure in the inner region of the galaxy, likely associated with an inner Lindblad resonance. The interaction between the large-scale gas flows in the bar and the gas in the inner region plausibly leads to the formation of massive molecular clouds and consequently massive clusters. This in turn can result in formation of massive stars, and thus the likely progenitor of the SLSN in a young, massive cluster. This picture is consistent with SLSN PTF10tpz being located near the inner structure. We find the molecular gas in the vicinity of the SLSN to have high surface densities, comparable with those in interacting galaxies or starburst regions in nearby galaxies. This lends support to high densities being favorable conditions for formation of SLSNe progenitors.

ALMA measurements for 93 $Herschel$-selected galaxies at $1.1 \leqslant z \leqslant 1.7$ in COSMOS reveal a sizable ($>29$\%) population with compact star formation (SF) sizes, lying on average $> \times 3.6$ below the optical stellar mass ($M_{\star}$)-size relation of disks. This sample widely spans the star-forming Main Sequence (MS), having $10^{8} \leqslant M_{\star} \leqslant 10^{11.5} \ M_{\odot}$ and $20 \leqslant SFR \leqslant 680 \ M_{\odot} \rm yr^{-1}$. The 32 size measurements and 61 upper limits are measured on ALMA images that combine observations of CO(5-4), CO(4-3), CO(2-1) and $\lambda_{\rm obs} \sim 1.1-1.3 \ \rm mm$ continuum, all tracing the star-forming molecular gas. These compact galaxies have instead normally extended $K_{band}$ sizes, suggesting strong specific $SFR$ gradients. Compact galaxies comprise the $50\pm18 \%$ of MS galaxies at $M_{\star} > 10^{11} M_{\odot}$. This is not expected in standard bi-modal scenarios where MS galaxies are mostly steadily-growing extended disks. We suggest that compact MS objects are early post-starburst galaxies in which the merger-driven boost of SF has subsided. They retain their compact SF size until either further gas accretion restores pre-merger galaxy-wide SF, or until becoming quenched. The fraction of merger-affected SF inside the MS seems thus larger than anticipated and might reach $\sim 50$\% at the highest $M_{\star}$. The presence of large galaxies above the MS demonstrates an overall poor correlation between galaxy SF size and specific $SFR$.

We report Giant Metrewave Radio Telescope (GMRT) , Very Large Telescope (VLT) and Spitzer Space Telescope observations of ESO 184$-$G82, the host galaxy of GRB 980425/SN 1998bw, that yield evidence of a companion dwarf galaxy at a projected distance of 13 kpc. The companion, hereafter GALJ193510-524947, is a gas-rich, star-forming galaxy with a star formation rate of $\rm0.004\,M_{\odot}\, yr^{-1}$, a gas mass of $10^{7.1\pm0.1} M_{\odot}$, and a stellar mass of $10^{7.0\pm0.3} M_{\odot}$. The interaction between ESO 184$-$G82 and GALJ193510-524947 is evident from the extended gaseous structure between the two galaxies in the GMRT HI 21 cm map. We find a ring of high column density HI gas, passing through the actively star forming regions of ESO 184$-$G82 and the GRB location. This ring lends support to the picture in which ESO 184$-$G82 is interacting with GALJ193510-524947. The massive stars in GALJ193510-524947 have similar ages to those in star-forming regions in ESO 184$-$G82, also suggesting that the interaction may have triggered star formation in both galaxies. The gas and star formation properties of ESO 184$-$G82 favour a head-on collision with GALJ193510-524947 rather than a classical tidal interaction. We perform state-of-the art simulations of dwarf--dwarf mergers and confirm that the observed properties of ESO 184$-$G82 can be reproduced by collision with a small companion galaxy. This is a very clear case of interaction in a gamma ray burst host galaxy, and of interaction-driven star formation giving rise to a gamma ray burst in a dense environment.

(Abridged) Long gamma-ray bursts (LGRB) have been suggested as promising tracers of star formation owing to their association with the core-collapse of massive stars. The goal of this work is to characterise the population of host galaxies of LGRBs at 1 < z < 2, investigate the conditions in which LGRBs form at these redshifts and assess their use as tracers of star formation. We perform a spectro-photometric analysis to determine the stellar mass, star formation rate, specific star formation rate and metallicity of the complete, unbiased host galaxy sample of the Swift/BAT6 LGRB sample at 1 < z < 2. We compare the distribution of these properties to the ones of typical star-forming galaxies from the MOSDEF and COSMOS2015 Ultra Deep surveys, within the same redshift range. We find that, similarly to z < 1, LGRBs do not directly trace star formation at 1 < z < 2, and they tend to avoid high-mass, high-metallicity host galaxies. We also find evidence for an enhanced fraction of starbursts among the LGRB host sample with respect to the star-forming population of galaxies. Nonetheless we demonstrate that the driving factor ruling the LGRB efficiency is metallicity. The LGRB host distributions can be reconciled with the ones expected from galaxy surveys by imposing a metallicity upper limit of 12+logOH ~ 8.55. Metallicity rules the LGRB production efficiency, which is stifled at Z > 0.7 Zsun. Under this hypothesis we can expect LGRBs to trace star formation at z > 3, once the bulk of the star forming galaxy population are characterised by metallicities below this limit. The moderately high metallicity threshold found is in agreement with the conditions necessary to rapidly produce a fast-rotating Wolf-Rayet star a in close binary system, and could be accommodated by single star models under chemically homogeneous mixing with very rapid rotation and weak magnetic coupling.

The bright [CI] 1-0 and [CI] 2-1 lines of atomic carbon are becoming more and more widely employed tracers of the cold neutral gas in high-redshift galaxies. Here we present observations of these lines in the 11 galaxies of the set of Planck's Dusty GEMS, the brightest gravitationally lensed galaxies on the extragalactic submillimeter sky probed by the Planck satellite. We have [CI] 1-0 measurements for seven, and [CI] 2-1 measurements for eight galaxies, including four galaxies where both lines are measured. We use our observations to constrain the gas excitation mechanism, excitation temperatures, optical depths, atomic carbon and molecular gas masses, and carbon abundances. Ratios of L_CI/L_ FIR are similar to those found in the local Universe, and suggest that the total cooling budget through atomic carbon has not strongly changed in the last 12 Gyr. Both lines are optically thin and trace 1 - 6 x 10^7 M_sun of atomic carbon. Carbon abundance ratios with H_2, X_CI, are between 2.5 and 4 x 10^-5, for a "ULIRG" CO-to-H_2 conversion factor of alpha_CO=0.8 M_sun/ [K km s^-1 pc^2]. Ratios of molecular gas masses derived from [CI] 1-0 and CO agree within the measurement uncertainties for five galaxies, and to better than a factor of 2 for another two with [CI] 1-0 measurements, after taking CO excitation carefully into account. This does not support the idea that intense, high-redshift starburst galaxies host large quantities of "CO-dark" gas. These results also support the common assumptions underlying most molecular gas mass estimates made for massive, dusty, high-redshift starburst galaxies, although the good agreement between the masses obtained with both tracers cannot be taken as an independent confirmation of either alpha_CO or X_CI.

We present an extensive CO emission-line survey of the Planck's dusty GEMS, a small set of 11 strongly lensed dusty star-forming galaxies at z = 2-4 discovered with Planck and Herschel satellites, using EMIR on the IRAM 30-m telescope. We detected a total of 45 CO rotational lines from Jup=3 to Jup=11, and up to eight transitions per source, allowing a detailed analysis of the gas excitation and interstellar medium conditions within these extremely bright, vigorous starbursts. We applied radiative transfer models using the large velocity gradient approach to infer the spatially-averaged molecular gas densities, $n_{H_2}$~10$^{2.6-4.1}$ cm$^{-3}$, and kinetic temperatures, $T_k$~30-1000 K. In five sources, we find evidence of two distinct gas phases with different properties and model their CO ladder with two excitation components. The warm (70-320 K) and dense gas reservoirs in these galaxies are highly excited, while the cooler (15-60 K) and more extended low-excitation components cover a range of gas densities. In two sources, the latter is associated with diffuse Milky Way-like gas phases, which provides evidence that a significant fraction of the total gas masses of dusty starburst galaxies can be embedded in cool, low-density reservoirs. Finally, we show that the CO line luminosity ratios are consistent with those predicted by models of photon-dominated regions and disfavor scenarios of gas clouds irradiated by intense X-ray fields from active galactic nuclei. By combining CO, [CI] and [CII] line diagnostics, we obtain average PDR gas densities significantly higher than in normal star-forming galaxies at low-redshift, and far-ultraviolet radiation fields at least 100 times more intense than in the Milky Way. These spatially-averaged conditions are consistent with those in high-redshift SMGs and in a range of low-redshift environments.

We report the discovery of a molecular wind signature from a massive intensely star-forming clump of a few $10^9$ Msun, in the strongly gravitationally lensed submillimeter galaxy "the Emerald" (PLCK_G165.7+49.0) at z=2.236. The Emerald is amongst the brightest high-redshift galaxies on the submillimeter sky, and was initially discovered with the Planck satellite. The system contains two magnificient structures with projected lengths of 28.5" and 21" formed by multiple, near-infrared arcs, falling behind a massive galaxy cluster at z=0.35, as well as an adjacent filament that has so far escaped discovery in other wavebands. We used HST/WFC3 and CFHT optical and near-infrared imaging together with IRAM and SMA interferometry of the CO(4-3) line and 850 $\mu$m dust emission to characterize the foreground lensing mass distribution, construct a lens model with Lenstool, and calculate gravitational magnification factors between 20 and 50 in most of the source. The majority of the star formation takes place within two massive star-forming clumps which are marginally gravitationally bound and embedded in a $9 \times 10^{10}$ Msun, fragmented disk with 20% gas fraction. One of the clumps shows a pronounced blue wing in the CO(4-3) line profile, which we interpret as a wind signature. The mass outflow rates are high enough for us to suspect that the clump might become unbound within a few tens of Myr, unless the outflowing gas can be replenished by gas accretion from the surrounding disk. The velocity offset of -200 km s$^{-1}$ is above the escape velocity of the clump, but not that of the galaxy overall, suggesting that much of this material might ultimately rain back onto the galaxy and contribute to fueling subsequent star formation.

There is a lack of large samples of spectroscopically confirmed clusters and protoclusters at high redshifts, $z>$1.5. Discovering and characterizing distant (proto-)clusters is important for yielding insights into the formation of large-scale structure and on the physical processes responsible for regulating star-formation in galaxies in dense environments. The Spitzer Planck Herschel Infrared Cluster (SPHerIC) survey was initiated to identify these characteristically faint and dust-reddened sources during the epoch of their early assembly. We present Spitzer IRAC observations of 82 galaxy (proto-)cluster candidates at 1.3<$z_p$<3.0 that were vetted in a two step process: (1) using Planck to select by color those sources with the highest star-formation rates, and (2) using Herschel at higher resolution to separate out the individual red sources. The addition of the Spitzer data enables efficient detection of the central and massive brightest red cluster galaxies (BRCGs). We find that BRCGs are associated with highly significant, extended and crowded regions of IRAC sources which are more overdense than the field. This result corroborates our hypothesis that BRCGs within the Planck - Herschel sources trace some of the densest and actively star-forming proto-clusters in the early Universe. On the basis of a richness-mass proxy relation, we obtain an estimate of their mean masses which suggests our sample consists of some of the most massive clusters at z$\approx$2 and are the likely progenitors of the most massive clusters observed today.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes

The potential similarity of the powering mechanisms of relativistic SNe and GRBs allowed us to make a prediction that relativistic SNe are born in environments similar to those of GRBs, that is, ones which are rich in atomic gas. Here we embark on testing this hypothesis by analysing the properties of the host galaxy NGC 3278 of the relativistic SN 2009bb. This is the first time the atomic gas properties of a relativistic SN host are provided and the first time resolved 21 cm-hydrogen-line (HI) information is provided for a host of an SN of any type in the context of the SN position. We obtained radio observations with ATCA covering the HI line, and optical integral field unit spectroscopy observations with MUSE. The atomic gas distribution of NGC 3278 is not centred on the optical galaxy centre, but instead around a third of atomic gas resides in the region close to the SN position. This galaxy has a few times lower atomic and molecular gas masses than predicted from its SFR. SN 2009bb exploded close to the region with the highest SFR density and the lowest age (~5.5 Myr, corresponding to the initial mass of the progenitor star ~36 Mo). As for GRB hosts, the gas properties of NGC 3278 are consistent with a recent inflow of gas from the intergalactic medium, which explains the concentration of atomic gas close to the SN position and the enhanced SFR. Super-solar metallicity at the position of the SN (unlike for most GRBs) may mean that relativistic explosions signal a recent inflow of gas (and subsequent star formation), and their type (GRBs or SNe) is determined either i) by the metallicity of the inflowing gas, so that metal-poor gas results in a GRB explosion and metal-rich gas results in a relativistic SN explosion without an accompanying GRB, or ii) by the efficiency of gas mixing, or iii) by the type of the galaxy.

The Infra-Red Telescope (IRT) on board the Transient High Energy Sky and Early Universe Surveyor (THESEUS) ESA M5 candidate mission will play a key role in identifying and characterizing moderate to high redshift Gamma-Ray Bursts afterglows. The IRT is the enabling instrument on board THESEUS for measuring autonomously the redshift of the several hundreds of GRBs detected per year by the Soft X-ray Imager (SXI) and the X- and Gamma-Ray Imaging Spectrometer (XGIS), and thus allowing the big ground based telescopes to be triggered on a redshift pre-selected sample, and finally fulfilling the cosmological goals of the mission. The IRT will be composed by a primary mirror of 0.7 m of diameter coupled to a single camera in a Cassegrain design. It will work in the 0.7-1.8 {\mu}m wavelength range, and will provide a 10x10 arc min imaging field of view with sub-arc second localization capabilities, and, at the same time, a 5x5 arc min field of view with moderate (R up to ~500) spectroscopic capabilities. Its sensitivity, mainly limited by the satellite jitter, is adapted to detect all the GRBs, localized by the SXI/XGIS, and to acquire spectra for the majority of them.

We report the detection of the CO(3-2) emission line from the host galaxy of Gamma Ray Burst (GRB) 080207 at $z$ = 2.086. This is the first detection of molecular gas in emission from a GRB host galaxy beyond redshift 1. We find this galaxy to be rich in molecular gas with a mass of $1.1 \times 10^{11}\,\rm M_{\odot}$ assuming $\alpha_{\rm CO}=$ 4.36 $\rm M_{\odot}(\rm K\,km\,s^{-1}\,pc^2)^{-1}$. The molecular gas mass fraction of the galaxy is $\sim$ 0.5, typical of star forming galaxies (SFGs) with similar stellar masses and redshifts. With a $\rm SFR_{FIR}$ of 260 $\rm M_{\odot}\,yr^{-1}$, we measure a molecular-gas-depletion timescale of 0.43 Gyr, near the peak of the depletion timescale distribution of SFGs at similar redshifts. Our findings are therefore in contradiction with the proposed molecular gas deficiency in GRB host galaxies. We argue that the reported molecular gas deficiency for GRB hosts could be the artifact of improper comparisons or neglecting the effect of the typical low metallicities of GRB hosts on the CO-to-molecular-gas conversion factor. We also compare the kinematics of the CO(3-2) emission line to that of the H$\alpha$ emission line from the host galaxy. We find the H$\alpha$ emission to have contributions from two separate components, a narrow and a broad one. The narrow component matches the CO emission well in velocity space. The broad component, with a FWHM of $\sim$ 1100 $\rm km\,s^{-1}$, is separated by $+390$ $\rm km\,s^{-1}$ in velocity space from the narrow component. We speculate this broad component to be associated with a powerful outflow in the host galaxy or in an interacting system.

We present a study of the connection between black hole accretion, star formation, and galaxy morphology at z<=2.5. We focus on active galactic nuclei (AGNs) selected by their mid-IR power-law emission. By fitting optical to far-IR photometry with state-of-the-art spectral energy distribution (SED) techniques, we derive stellar masses, star formation rates, dust properties, and AGN contributions in galaxies over the whole COSMOS field. We find that obscured AGNs lie within or slightly above the star-forming sequence. We confirm our previous finding about compact host galaxies of obscured AGNs at z~1, and find that galaxies with 20-50% AGN contributions tend to have smaller sizes, by ~25-50%, compared to galaxies without AGNs. Furthermore, we find that a high merger fraction of up to 0.5 is appropriate for the most luminous (log (LIR/Lsun) ~ 12.5) AGN hosts and non-AGN galaxies, but not for the whole obscured AGN sample. Moreover, merger fraction depends on the total and star-forming infrared luminosity, rather than the decomposed AGN infrared luminosity. Our results suggest that major mergers are not the main driver of AGN activity, and therefore obscured AGNs might be triggered by internal mechanisms, such as secular processes, disk instabilities, and compaction in a particular evolutionary stage. We make the SED modeling results publicly available.

We present a comprehensive study of the relations between gas kinematics, metallicity, and stellar mass in a sample of 82 GRB-selected galaxies using absorption and emission methods. We find the velocity widths of both emission and absorption profiles to be a proxy of stellar mass. We also investigate the velocity-metallicity correlation and its evolution with redshift and find the correlation derived from emission lines to have a significantly smaller scatter compared to that found using absorption lines. Using 33 GRB hosts with measured stellar mass and metallicitiy, we study the mass-metallicity relation for GRB host galaxies in a stellar mass range of $10^{8.2} M_{\odot}$ to $10^{11.1} M_{\odot}$ and a redshift range of $ z\sim 0.3-3.4$. The GRB-selected galaxies appear to track the mass-metallicity relation of star forming galaxies but with an offset of 0.15 towards lower metallicities. This offset is comparable with the average error-bar on the metallicity measurements of the GRB sample and also the scatter on the MZ relation of the general population. It is hard to decide whether this relatively small offset is due to systematic effects or the intrinsic nature of GRB hosts. We also investigate the possibility of using absorption-line metallicity measurements of GRB hosts to study the mass-metallicity relation at high redshifts. Our analysis shows that the metallicity measurements from absorption methods can significantly differ from emission metallicities and assuming identical measurements from the two methods may result in erroneous conclusions.

We present an analysis of high-resolution ALMA interferometry of CO(4-3) line emission and dust continuum in the "Ruby" (PLCK_G244.8+54.9), a bright, gravitationally lensed galaxy at z = 3.0 discovered with the Planck all-sky survey. The Ruby is the brightest of Planck's Dusty GEMS, a sample of 11 of the brightest gravitationally lensed high-redshift galaxies on the extragalactic sub-mm sky. We resolve the high-surface-brightness continuum and CO line emission of the Ruby in several extended clumps along a partial, nearly circular Einstein ring with 1.4" diameter around a massive galaxy at z = 1.5. Local star-formation intensities are up to 2000 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, amongst the highest observed at high redshift, and clearly in the range of maximal starbursts. Gas-mass surface densities are a few $\times$ 10$^4$ M$_{\odot}$ pc$^{-2}$. The Ruby lies at, and in part even above, the starburst sequence in the Schmidt-Kennicutt diagram, and at the limit expected for star formation that is self-regulated through the kinetic energy injection from radiation pressure, stellar winds, and supernovae. We show that these processes can also inject sufficient kinetic energy and momentum into the gas to explain the turbulent line widths, which are consistent with marginally gravitationally bound molecular clouds embedded in a critically Toomre-stable disk. The star-formation efficiency is in the range 1-10% per free-fall time, consistent with the notion that the pressure balance that sets the local star-formation law in the Milky Way may well be universal out to the highest star-formation intensities. AGN feedback is not necessary to regulate the star formation in the Ruby, in agreement with the absence of a bright AGN component in the infrared and radio regimes.

We present deep ALMA CO(5-4) observations of a main sequence, clumpy galaxy at z=1.5 in the HUDF. Thanks to the ~0.5" resolution of the ALMA data, we can link stellar population properties to the CO(5-4) emission on scales of a few kpc. We detect strong CO(5-4) emission from the nuclear region of the galaxy, consistent with the observed $L_{\rm IR}$-$L^{\prime}_{\rm CO(5-4)}$ correlation and indicating on-going nuclear star formation. The CO(5-4) gas component appears more concentrated than other star formation tracers or the dust distribution in this galaxy. We discuss possible implications of this difference in terms of star formation efficiency and mass build-up at the galaxy centre. Conversely, we do not detect any CO(5-4) emission from the UV-bright clumps. This might imply that clumps have a high star formation efficiency (although they do not display unusually high specific star formation rates) and are not entirely gas dominated, with gas fractions no larger than that of their host galaxy (~50%). Stellar feedback and disk instability torques funnelling gas towards the galaxy centre could contribute to the relatively low gas content. Alternatively, clumps could fall in a more standard star formation efficiency regime if their actual star-formation rates are lower than generally assumed. We find that clump star-formation rates derived with several different, plausible methods can vary by up to an order of magnitude. The lowest estimates would be compatible with a CO(5-4) non-detection even for main-sequence like values of star formation efficiency and gas content.

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]

With the increasing number of deep multi-wavelength galaxy surveys, the spectral energy distribution (SED) of galaxies has become an invaluable tool for studying the formation of their structures and their evolution. In this context, standard analysis relies on simple spectro-photometric selection criteria based on a few SED colors. If this fully supervised classification already yielded clear achievements, it is not optimal to extract relevant information from the data. In this article, we propose to employ very recent advances in machine learning, and more precisely in feature learning, to derive a data-driven diagram. We show that the proposed approach based on denoising autoencoders recovers the bi-modality in the galaxy population in an unsupervised manner, without using any prior knowledge on galaxy SED classification. This technique has been compared to principal component analysis (PCA) and to standard color/color representations. In addition, preliminary results illustrate that this enables the capturing of extra physically meaningful information, such as redshift dependence, galaxy mass evolution and variation over the specific star formation rate. PCA also results in an unsupervised representation with physical properties, such as mass and sSFR, although this representation separates out. less other characteristics (bimodality, redshift evolution) than denoising autoencoders.

We study the properties of the foreground galaxy of the Ruby, the brightest gravitationally lensed high-redshift galaxy on the sub-millimeter sky as probed by the Planck satellite, and part of our sample of Planck's Dusty GEMS. The Ruby consists of an Einstein ring of 1.4" diameter at z = 3.005 observed with ALMA at 0.1" resolution, centered on a faint, red, massive lensing galaxy seen with HST/WFC3, which itself has an exceptionally high redshift, z = 1.525 $\pm$ 0.001, as confirmed with VLT/X-Shooter spectroscopy. Here we focus on the properties of the lens and the lensing model obtained with LENSTOOL. The rest-frame optical morphology of this system is strongly dominated by the lens, while the Ruby itself is highly obscured, and contributes less than 10% to the photometry out to the K band. The foreground galaxy has a lensing mass of (3.70 $\pm$ 0.35) $\times$ 10$^{11}$ M$_{\odot}$. Magnification factors are between 7 and 38 for individual clumps forming two image families along the Einstein ring. We present a decomposition of the foreground and background sources in the WFC3 images, and stellar population synthesis modeling with a range of star-formation histories for Chabrier and Salpeter initial mass functions (IMFs). Only the stellar mass range obtained with the latter agrees well with the lensing mass. This is consistent with the bottom-heavy IMFs of massive high-redshift galaxies expected from detailed studies of the stellar masses and mass profiles of their low-redshift descendants, and from models of turbulent gas fragmentation. This may be the first direct constraint on the IMF in a lens at z = 1.5, which is not a cluster central galaxy.

We investigate of the properties of $\sim$2000 Herschel/SPIRE-selected galaxies from $0<z<4$ using a combination of extensive spectroscopy, deep imaging from CFHT, VLA, Spitzer, XMM-Newton, and Herschel, and well-calibrated SED fitting. Herschel galaxies are observed to span a range of stellar masses, colors, and absolute magnitudes equivalent to galaxies undetected in SPIRE. Though many Herschel galaxies appear to be in transition, such galaxies are largely consistent with normal star-forming galaxies when rest-frame colors are utilized. The nature of the star-forming "main sequence" is studied and we warn against adopting this framework unless the main sequence is determined precisely. Herschel galaxies at different total infrared luminosities ($L_{TIR}$) are compared. Bluer colors, larger nebular extinctions, and larger contributions from younger stellar populations are observed for galaxies with larger $L_{TIR}$, suggesting that low-$L_{TIR}$ galaxies are undergoing rejuvenated starbursts while galaxies with higher $L_{TIR}$ are forming a larger percentage of their stellar mass. A variety of methods are used to select powerful active galactic nuclei (AGN). Galaxies hosting AGN are observed to be undergoing starbursts more commonly and vigorously than a matched sample of galaxies without powerful AGN and, additionally, the fraction of galaxies with an AGN increases with increasing star formation rate at all redshifts. At all redshifts ($0<z<4$) the most prodigious star-forming galaxies are found to contain the highest fraction of powerful AGN. For redshift bins that allow a comparison ($z>0.5$), the highest $L_{TIR}$ galaxies in a given redshift bin are unobserved by SPIRE at subsequently lower redshifts, a trend linked to downsizing. In conjunction with other results, this evidence is used to argue for prevalent AGN-driven quenching in starburst galaxies across cosmic time.

We present a structural study of 182 obscured Active Galactic Nuclei (AGNs) at z<=1.5, selected in the COSMOS field from their extreme infrared to X-ray luminosity ratio and their negligible emission at optical wavelengths. We fit optical to far-infrared spectral energy distributions and analyze deep HST imaging to derive the physical and morphological properties of their host galaxies. We find that such galaxies are more compact than normal star-forming sources at similar redshift and stellar mass, and we show that it is not an observational bias related to the emission of the AGN. Based on the distribution of their UVJ colors, we also argue that this increased compactness is not due to the additional contribution of a passive bulge. We thus postulate that a vast majority of obscured AGNs reside in galaxies undergoing dynamical compaction, similar to processes recently invoked to explain the formation of compact star-forming sources at high redshift.

We present spatially resolved ALMA [CII] observations of the bright (flux density S=400 mJy at 350 microns), gravitationally lensed, starburst galaxy PLCK G045.1+61.1 at z=3.427, the "Garnet". This source is part of our set of "Planck's Dusty GEMS", discovered with the Planck's all-sky survey. Two emission-line clouds with a relative velocity offset of ~600 km/s extend towards north-east and south-west, respectively, of a small, intensely star-forming clump with a star-formation intensity of 220 Msun/yr/kpc^2, akin to maximal starbursts. [CII] is also seen in absorption, with a redshift of +350 km/s relative to the brightest CO component. [CII] absorption has previously only been found in the Milky Way along sightlines toward bright high-mass star-forming regions, and this is the first detection in another galaxy. Similar to Galactic environments, the [CII] absorption feature is associated with [CI] emission, implying that this is diffuse gas shielded from the UV radiation of the clump, and likely at large distances from the clump. Since absorption can only be seen in front of a continuum source, the gas in this structure can definitely be attributed to gas flowing towards the clump. The absorber could be part of a cosmic filament or merger debris being accreted onto the galaxy. We discuss our results also in light of the on-going debate of the origin of the [CII] deficit in dusty star-forming galaxies.

We have recently suggested that gas accretion can be studied using host galaxies of gamma-ray bursts (GRBs). We obtained the first ever far-infrared (FIR) line observations of a GRB host, namely Herschel/PACS resolved [CII] 158 um and [OI] 63 um spectroscopy, as well as APEX CO(2-1) and ALMA CO(1-0) observations of the GRB 980425 host. It has elevated [CII]/FIR and [OI]/FIR ratios and higher values of star formation rate (SFR) derived from line ([CII], [OI], Ha) than from continuum (UV, IR, radio) indicators. [CII] emission exhibits a normal morphology, peaking at the galaxy center, whereas [OI] is concentrated close to the GRB position and the nearby Wolf-Rayet region. The high [OI] flux indicates high radiation field and gas density. The [CII]/CO luminosity ratio of the GRB 980425 host is close to the highest values found for local star-forming galaxies. Its CO-derived molecular gas mass is low given its SFR and metallicity, but the [CII]-derived molecular gas mass is close to the expected value. The [OI] and HI concentrations, and the high radiation field and density are consistent with the hypothesis of a very recent (at most a few tens of Myr ago) inflow of atomic gas triggering star formation. Dust has not had time to build up (explaining high line-to-continuum ratios). Such a recent enhancement of star-formation would indeed manifest itself in high SFR_line/SFR_continuum ratios, because the line indicators are sensitive only to recent (<10 Myr) activity, whereas the continuum indicators measure the SFR averaged over much longer periods (~100 Myr). Other GRB hosts exhibit a mean SFR_line/SFR_continuum of 1.74+-0.32. This is consistent with a very recent enhancement of star formation being common among GRB hosts, so galaxies which have recently experienced inflow of gas may preferentially host stars exploding as GRBs. Hence GRB hosts may be used to investigate recent gas accretion.

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

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

Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and could thus be a potentially powerful tool to trace cosmic star formation. However, especially at low redshifts (z < 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the impact of the environment on GRB formation. We study host galaxy spectra of the Swift/BAT6 complete sample of 14 z < 1 bright LGRBs. We use the detected nebular emission lines to measure the dust extinction, star formation rate (SFR) and nebular metallicity (Z) of the hosts and supplement the data set with previously measured stellar masses M$_{\star}$. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M$_{\star}$ relations) are compared to samples of field star-forming galaxies.We find that LGRB hosts at z < 1 have on average lower SFRs than if they were direct star-formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to Z $\sim 8.4-8.5$, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are not inconsistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of increased incidence of starburst galaxies in the Swift/BAT6 z < 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevent us from reaching more quantitative conclusions.

We present an analysis of the morphological classification of 89 luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS) sample using non-parametric coefficients and compare their morphology as a function of wavelength. We rely on images obtained in the optical (B- and I-band) as well as in the infrared (H-band and 5.8$\mu$m). Our classification is based on the calculation of $Gini$ and the second order of light ($M_{20}$) non-parametric coefficients which we explore as a function of stellar mass ($M_\star$), infrared luminosity ($L_{IR}$) and star formation rate (SFR). We investigate the relation between $M_{20}$, the specific SFR (sSFR) and the dust temperature ($T_{dust}$) in our galaxy sample. We find that $M_{20}$ is a better morphological tracer than $Gini$, as it allows to distinguish systems formed by double systems from isolated and post-merger LIRGs. The multi-wavelength analysis allows us to identify a region in the $Gini$-$M_{20}$ parameter space where ongoing mergers reside, regardless of the band used to calculate the coefficients. In particular when measured in the H-band, this region can be used to identify ongoing mergers, with a minimal contamination from LIRGs in other stages. We also find that while the sSFR is positively correlated with $M_{20}$ when measured in the mid-infrared, i.e. star-bursting galaxies show more compact emission, it is anti-correlated with the B-band based $M_{20}$. We interpret this as the spatial decoupling between obscured and un-obscured star formation, whereby the ultraviolet/optical size of a LIRGs experience an intense dust enshrouded central starburst is larger than in the one in the mid-infrared since the contrast between the nuclear to the extended disk emission is smaller in the mid-infrared. This has important implications for high redshift surveys of dusty sources. [abridged]

We present new data for four candidate obscured Compton-Thick (CT) quasars at z $\sim$1-2.5 observed with SINFONI VLT spectrograph in AO mode. These sources were selected from a 24$\mu$m Spitzer MIPS survey of the COSMOS field, on the basis of red mid-infrared-to-optical and optical-to-near-infrared colours, with the intention of identifying active galactic nuclei (AGNs) in dust enshrouded environments, where most of the black hole mass is assembled in dust enshrouded environments. Near infrared spectra were analyzed in order to check for emission line features and to search for broad components in the [OIII]-H$\beta$ and H$\alpha$-[NII] regions. X-ray spectral analysis, radio and MIR diagnostics, and SED fitting have also been employed to study the nature of the sources. We successfully identified three objects for which we had only a photometric redshift estimate. Based on their emission line diagnostics and on ancillary multi-wavelength constraints, we find that all four targets harbor obscured AGNs. Broad profiles that could be attributed to the effects of outflows are revealed in only one target, MIRO20581. In particular, we clearly resolved a fast ($\sim$1600 km/s) and extended ($\sim$5 kpc) outflow in the [OIII]5007 emission line. This feature, the commonly used indicator for ionised outflowing gas, was sampled and detected only for this target; hence, we can not exclude the presence of outflows in the other sources. Overall, the constraints we obtain from our targets and from other comparative samples from the literature suggest that these optically faint luminous infrared galaxies, hosting obscured AGNs, may represent a brief evolutionary phase between the post-merger starburst and the unobscured QSO phases.

Long gamma-ray bursts (GRBs), among the most energetic events in the Universe, are explosions of massive and short-lived stars, so they pinpoint locations of recent star formation. However, several GRB host galaxies have recently been found to be deficient in molecular gas (H2), believed to be the fuel of star formation. Moreover, optical spectroscopy of GRB afterglows implies that the molecular phase constitutes only a small fraction of the gas along the GRB line-of-sight. Here we report the first ever 21 cm line observations of GRB host galaxies, using the Australia Telescope Compact Array, implying high levels of atomic hydrogen (HI), which suggests that the connection between atomic gas and star formation is stronger than previously thought, with star formation being potentially directly fuelled by atomic gas (or with very efficient HI-to-H2 conversion and rapid exhaustion of molecular gas), as has been theoretically shown to be possible. This can happen in low metallicity gas near the onset of star formation, because cooling of gas (necessary for star formation) is faster than the HI-to-H2 conversion. Indeed, large atomic gas reservoirs, together with low molecular gas masses, stellar and dust masses are consistent with GRB hosts being preferentially galaxies which have very recently started a star formation episode after accreting metal-poor gas from the intergalactic medium. This provides a natural route for forming GRBs in low-metallicity environments. The gas inflow scenario is also consistent with the existence of the companion HI object with no optical counterpart ~19 kpc from the GRB 060505 host, and with the fact that the HI centroids of the GRB 980425 and 060505 hosts do not coincide with optical centres of these galaxies, but are located close to the GRB positions.

LGRBs are associated with massive stars and are therefore linked to star formation. The conditions necessary to produce LGRBs can affect the relation between the LGRB rate and star formation. By using the power of a complete LGRB sample, our aim is to understand whether such a bias exists and, if it does, what is its origin. In this first paper, we build the SED of the z<1 host galaxies of the BAT6 LGRB sample, and determine their stellar masses from SED fitting. We compare the resulting stellar mass distribution (i) with star-forming galaxies observed in deep surveys (UltraVISTA); (ii) with semi-analitical models of the z<1 star forming galaxy population and (iii) with numerical simulations of LGRB hosts having different metallicity thresholds for the progenitor star environment. We find that at z<1 LGRBs tend to avoid massive galaxies and are powerful in selecting faint low-mass star-forming galaxies. The stellar mass distribution of the hosts is not consistent with that of the UltraVISTA star-forming galaxies weighted for their SFR. This implies that, at least at z<1, LGRBs are not unbiased tracers of star formation. To make the two distributions consistent, a much steeper faint-end of the mass function would be required, or a very shallow SFR-Mass relation for the low mass galaxy population. GRB host galaxy simulations indicates that, to reproduce the stellar mass distribution, a metallicity threshold of the order of Z_th=0.3-0.5Z_sun is necessary. Models without a metallicity threshold or with an extreme threshold of Z_th = 0.1Z_sun are excluded at z<1. The use of the BAT6 complete sample makes this result not affected by possible biases which could have influenced past results based on incomplete samples. The preference for low metallicities (Z<~0.5Z_sun) can be a consequence of the particular conditions needed for the progenitor star to produce a GRB. (Abridged)

We study the properties of satellites in the environment of massive star-forming galaxies at z~1.8 in the COSMOS field, using a sample of 215 galaxies on the main sequence of star formation with an average mass of 10^11 Msun. At z>1.5, these galaxies typically trace halos of mass >10^13 Msun. We use optical-near-infrared photometry to estimate stellar masses and star formation rates (SFR) of centrals and satellites down to ~6*10^9 Msun. We stack data around 215 central galaxies to statistically detect their satellite halos, finding an average of ~3 galaxies in excess of the background density. We fit the radial profiles of satellites with simple beta-models, and compare their integrated properties to model predictions. We find that the total stellar mass of satellites amounts to 68% of the central galaxy, while SED modeling and far-infrared photometry consistently show their total SFR to be 25-35% of the central's rate. We also see significant variation in the specific SFR of satellites within the halo with, in particular, a sharp decrease at <100 kpc. After considering different potential explanations, we conclude that this is likely an environmental signature of the hot inner halo. This effect can be explained in the first order by a simple free-fall scenario, suggesting that these low-mass environments can shut down star formation in satellites on relatively short timescales of ~0.3 Gyr.

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

At z>1, the distinction between merging and 'normal' star-forming galaxies based on single band morphology is often hampered by the presence of large clumps which result in a disturbed, merger-like appearance even in rotationally supported disks. In this paper we discuss how a classification based on canonical, non-parametric structural indices measured on resolved stellar mass maps, rather than on single-band images, reduces the misclassification of clumpy but not merging galaxies. We calibrate the mass-based selection of mergers using the MIRAGE hydrodynamical numerical simulations of isolated and merging galaxies which span a stellar mass range of $10^{9.8}$-$10^{10.6}M_{sun}$ and merger ratios between 1:1-1:6.3. These simulations are processed to reproduce the typical depth and spatial resolution of observed HUDF data. We test our approach on a sample of real z~2 galaxies with kinematic classification into disks or mergers and on ~100 galaxies in the HUDF field with photometric/spectroscopic redshift between 1.5$\leqslant z \leqslant$3 and $M>10^{9.4}M_{sun}$. We find that a combination of the asymmetry $A_{MASS}$ and $M_{20, MASS}$ indices measured on the stellar mass maps can efficiently identify real (major) mergers with $\lesssim$20% contamination from clumpy disks in the merger sample. This mass-based classification cannot be reproduced in star-forming galaxies by $H-$band measurements alone, which instead result in a contamination from clumpy galaxies that can be as high as 50%. Moreover, we find that the mass-based classification always results in a lower contamination from clumpy galaxies than an $H-$band classification, regardless of the depth of the imaging used (e.g., CANDELS vs. HUDF).

When the cosmic star formation history peaks (z ~ 2), galaxies vigorously fed by cosmic reservoirs are gas dominated and contain massive star-forming clumps, thought to form by violent gravitational instabilities in highly turbulent gas-rich disks. However, a clump formation event has not been witnessed yet, and it is debated whether clumps survive energetic feedback from young stars, thus migrating inwards to form galaxy bulges. Here we report spatially resolved spectroscopy of a bright off-nuclear emission line region in a galaxy at z = 1.987. Although this region dominates the star formation in the galaxy disk, its stellar continuum remains undetected in deep imaging, revealing an extremely young (age < 10 Myr) massive clump, forming through the gravitational collapse of > 10$^9$ M$_{\odot}$ of gas. Gas consumption in this young clump is > 10 times faster than in the host galaxy, displaying high star formation efficiency during this phase, in agreement with our hydrodynamic simulations. The frequency of older clumps with similar masses coupled with our initial estimate of their formation rate (~ 2.5 Gyr$^{-1}$) supports long lifetimes (~ 500 Myr), favouring scenarios where clumps survive feedback and grow the bulges of present-day galaxies.

[Abridged] We use the Planck all-sky submm and mm maps to search for rare sources distinguished by extreme brightness, a few hundreds of mJy, and their potential for being situated at high redshift. These "cold" Planck sources, selected using the High Frequency Instrument (HFI) directly from the maps and from the Planck Catalogue of Compact Sources (PCCS), all satisfy the criterion of having their rest-frame far-infrared peak redshifted to the frequency range 353 and 857 GHz. This colour-selection favours galaxies in the redshift range z=2-4, which we consider as cold peaks in the cosmic infrared background (CIB). We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500um sources, in comparison to reference SPIRE fields. About 94% of the SPIRE sources in the Planck fields are consistent with being overdensities of galaxies peaking at 350um. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z>2.2. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10sigma. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z~2. Under the Td=35K assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4x10^12 Lsun, yielding star formation rates of typically 700 Msun.yr^-1. If the observed overdensities are actual gravitationally-bound structures, the total total star formation rates reaches 7x10^3 Msun.yr^-1. Taken together, these sources show the signatures of high-z (z>$) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.

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.

(Abridged) The knowledge of the number and of the physical nature of low-metallicity massive galaxies is crucial for the determination and interpretation of the mass-metallicity relation (MZR). Using VLT-ISAAC near-infrared (NIR) spectroscopy of 39 zCOSMOS z~0.7 galaxies, we have measured Halpha and [NII] emission line fluxes for galaxies with [OII], Hbeta and [OIII] available from VIMOS optical spectroscopy. The NIR spectroscopy enables us to break the degeneracy of the R23 method to derive unambiguously O/H gas metallicities, and also SFRs from extinction corrected Halpha. Using, as a benchmark, the position in the D4000 vs. [OIII]/Hbeta diagram of galaxies with reliable O/Hs from NIR spectroscopy, we were able to break the lower/upper branch R23 degeneracy of additional 900 zCOSMOS z~0.7 galaxies. Additionally, the Halpha-based SFR measurements were used to find the best SFR calibration based on [OII] for the zCOSMOS z~0.7 galaxies without Halpha measurements. We find a fraction of 19% of lower mass 9.5<logM/Msun<10.3 zCOSMOS galaxies which shows a larger evolution of the MZR relation, compared to higher mass galaxies, being more metal poor at a given mass by a factor of 2-3 compared to SDSS. This indicates that the low-mass MZR slope is getting steeper at z~0.7 compared to local galaxies. The existence of these metal-poor galaxies at z~0.7 can be interpreted as the chemical version of galaxy downsizing. Moreover, the sample of zCOSMOS galaxies shows direct evidence that SFR influences the MZR at these redshifts. The comparison of the measured metallicities for the zCOSMOS sample with the values expected for a non-evolving fundamental metallicity relation (FMR) shows broadly agreement, and reveals that also galaxies with lower metallicities and typically higher (specific) SFRs, as found in our zCOSMOS sample at z~0.7, are in agreement with the predictions of a non-evolving Z(M,SFR).

Even 10 billion years ago, the cores of the first galaxy clusters are often found to host a characteristic population of massive galaxies with already suppressed star formation. Here we search for distant cluster candidates at z~2 using massive passive galaxies as tracers. With a sample of ~40 spectroscopically confirmed passive galaxies at 1.3<z<2.1, we tune photometric redshifts of several thousands passive sources in the full 2 sq.deg. COSMOS field. This allows us to map their density in redshift slices, probing the large scale structure in the COSMOS field as traced by passive sources. We report here on the three strongest passive galaxy overdensities that we identify in the redshift range 1.5<z<2.5. While the actual nature of these concentrations is still to be confirmed, we discuss their identification procedure, and the arguments supporting them as candidate galaxy clusters (likely mid-10^13 M_sun range). Although this search approach is likely biased towards more evolved structures, it has the potential to select still rare, cluster-like environments close to their epoch of first appearance, enabling new investigations of the evolution of galaxies in the context of structure growth.

[abridged] We use the latest release of CIGALE, a galaxy SED fitting model relying on energy balance, to study the influence of an AGN in estimating both the SFR and stellar mass in galaxies, as well as the contribution of the AGN to the power output of the host. Using the galaxy formation SAM GALFORM, we create mock galaxy SEDs using realistic star formation histories (SFH) and add an AGN of Type 1, Type 2, or intermediate type whose contribution to the bolometric luminosity can be variable. We perform an SED fitting of these catalogues with CIGALE assuming three different SFHs: a single- and double-exponentially-decreasing, and a delayed SFH. Constraining thecontribution of an AGN to the LIR (fracAGN) is very challenging for fracAGN<20%, with uncertainties of ~5-30% for higher fractions depending on the AGN type, while FIR and sub-mm are essential. The AGN power has an impact on the estimation of $M_*$ in Type 1 and intermediate type AGNs but has no effect for galaxies hosting Type 2 AGNs. We find that in the absence of AGN emission, the best estimates of $M_*$ are obtained using the double-exponentially-decreasing model but at the expense of realistic ages of the stellar population. The delayed SFH model provides good estimates of $M_*$ and SFR, with a maximum offset of 10% as well as better estimates of the age. Our analysis shows that the underestimation of the SFR increases with fracAGN for Type 1 systems, as well as for low contributions of an intermediate AGN type, but it is quite insensitive to the emission of Type 2 AGNs up to fracAGN~45%. Similarly the UV emission is critical in accurately retrieving the $M_*$ for Type 1 and intermediate type AGN, and the SFR of all of the three AGN types. We show that the presence of AGN emission introduces a scatter to the SFR-$M_*$ main sequence relation derived from SED fitting, which is driven by the uncertainties on $M_*$.

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 present stringent constraints on the average mid-, far-infrared and radio emissions of $\sim$14200 quiescent galaxies (QGs), identified out to $z=3$ in the COSMOS field via their rest-frame NUV$-$r and r$-$J colors, and with stellar masses $M_{\star}=10^{9.8-12.2} \,M_{\odot} $. Stacking in deep Spitzer (MIPS $24\,\mu$m), Herschel (PACS and SPIRE), and VLA (1.4 GHz) maps reveals extremely low dust-obscured star formation rates for QGs (SFR $<0.1-3\,M_{\odot}$yr$^{-1}$ at $z \leqslant 2$ and $<6-18\,M_{\odot}$yr$^{-1}$ at $z > 2$), consistent with the low unobscured SFRs ($<0.01-1.2\,M_{\odot}$yr$^{-1}$) inferred from modeling their ultraviolet-to-near-infrared photometry. The average SFRs of QGs are $>10\times$ below those of star-forming galaxies (SFGs) within the $M_{\star}$- and $z$-ranges considered. The stacked 1.4 GHz signals (S/N $> 5$) are, if attributed solely to star formation, in excess of the total (obscured plus unobscured) SFR limits, suggestive of a widespread presence of low-luminosity active galactic nuclei (AGN) among QGs. Our results reaffirm the existence of a significant population QGs out to $z = 3$, thus corroborating the need for powerful quenching mechanism(s) to terminate star formation in galaxies at earlier epochs.

[Abridged] We study the evolution of the radio spectral index and far-infrared/radio correlation (FRC) across the star-formation rate-stellar masse (i.e. SFR-M*) plane up to z 2. We start from a M*-selected sample of galaxies with reliable SFR and redshift estimates. We then grid the SFR-M* plane in several redshift ranges and measure the infrared luminosity, radio luminosity, radio spectral index, and ultimately the FRC index (i.e. qFIR) of each SFR-M*-z bin. The infrared luminosities of our SFR-M*-z bins are estimated using their stacked far-infrared flux densities inferred from observations obtained with Herschel. Their radio luminosities and radio spectral indices (i.e. alpha, where Snu nu^-alpha) are estimated using their stacked 1.4GHz and 610MHz flux densities from the VLA and GMRT, respectively. Our far-infrared and radio observations include the most widely studied blank extragalactic fields -GOODS-N/S, ECDFS, and COSMOS- covering a sky area of 2deg^2. Using this methodology, we constrain the radio spectral index and FRC index of star-forming galaxies with M*>10^10Msun and 0<z<2.3. We find that alpha^1.4GHz_610MHz does not evolve significantly with redshift or with the distance of a galaxy with respect to the main sequence (MS) of the SFR-M* plane (i.e. Delta_log(SSFR)_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). Instead, star-forming galaxies have a radio spectral index consistent with a canonical value of 0.8, which suggests that their radio spectra are dominated by non-thermal optically thin synchrotron emission. We find that qFIR displays a moderate but statistically significant redshift evolution as qFIR(z)=(2.35+/-0.08)*(1+z)^(-0.12+/-0.04), consistent with some previous literature. Finally, we find no significant correlation between qFIR and Delta_log(SSFR)_MS, though a weak positive trend, as observed in one of our redshift bins, cannot be firmly ruled out using our dataset.

Galaxies' rest-frame ultraviolet (UV) properties are often used to directly infer the degree to which dust obscuration affects the measurement of star formation rates. While much recent work has focused on calibrating dust attenuation in galaxies selected at rest-frame ultraviolet wavelengths, locally and at high-$z$, here we investigate attenuation in dusty, star-forming galaxies (DSFGs) selected at far-infrared wavelengths. By combining multiwavelength coverage across 0.15--500\,$\mu$m in the COSMOS field, in particular making use of {\it Herschel} imaging, and a rich dataset on local galaxies, we find a empirical variation in the relationship between rest-frame UV slope ($\beta$) and ratio of infrared-to-ultraviolet emission ($L_{\rm IR}/L_{\rm UV}\equiv\,IRX$) as a function of infrared luminosity, or total star formation rate, SFR. Both locally and at high-$z$, galaxies above SFR$\gt$50\,M$_\odot$\,yr$^{-1}$ deviate from the nominal $IRX-\beta$ relation towards bluer colors by a factor proportional to their increasing IR luminosity. We also estimate contamination rates of DSFGs on high-$z$ dropout searches of $\ll1$\%\ at $z\lt4-10$, providing independent verification that contamination from very dusty foreground galaxies is low in LBG searches. Overall, our results are consistent with the physical interpretation that DSFGs, e.g. galaxies with $>50$\,M$_\odot$\,yr$^{-1}$, are dominated at all epochs by short-lived, extreme burst events, producing many young O and B stars that are primarily, yet not entirely, enshrouded in thick dust cocoons. The blue rest-frame UV slopes of DSFGs are inconsistent with the suggestion that most DSFGs at $z\sim2$ exhibit steady-state star formation in secular disks.

We compare the morphologies of a sample of 20 LIRGs from the Great Observatories All-sky LIRG Survey (GOALS) in the FUV, B, I and H bands, using the Gini (G) and M20 parameters to quantitatively estimate the distribution and concentration of flux as a function of wavelength. HST images provide an average spatial resolution of ~80 pc. While our LIRGs can be reliably classified as mergers across the entire range of wavelengths studied here, there is a clear shift toward more negative M20 (more bulge-dominated) and a less significant decrease in G values at longer wavelengths. We find no correlation between the derived FUV G-M20 parameters and the global measures of the IR to FUV flux ratio, IRX. Given the fine resolution in our HST data, this suggests either that the UV morphology and IRX are correlated on very small scales, or that the regions emitting the bulk of the IR emission emit almost no FUV light. We use our multi-wavelength data to simulate how merging LIRGs would appear from z~0.5-3 in deep optical and near-infrared images such as the HUDF, and use these simulations to measure the G-M20 at these redshifts. Our simulations indicate a noticeable decrease in G, which flattens at z >= 2 by as much as 40%, resulting in mis-classifying our LIRGs as disk-like, even in the rest-frame FUV. The higher redshift values of M20 for the GOALS sources do not appear to change more than about 10% from the values at z~0. The change in G-M20 is caused by the surface brightness dimming of extended tidal features and asymmetries, and also the decreased spatial resolution which reduced the number of individual clumps identified. This effect, seen as early as z~0.5, could easily lead to an underestimate of the number of merging galaxies at high-redshift in the rest-frame FUV.

Using the first 50% of data collected for the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH) observations on the 1.8 deg$^2$ Cosmological Evolution Survey (COSMOS) we estimate the masses and star formation rates of 3398 $M_*>10^{10}M_\odot $ star-forming galaxies at $4 < z < 6$ with a substantial population up to $M_* \gtrsim 10^{11.5} M_\odot$. We find that the strong correlation between stellar mass and star formation rate seen at lower redshift (the "main sequence" of star-forming galaxies) extends to $z\sim6$. The observed relation and scatter is consistent with a continued increase in star formation rate at fixed mass in line with extrapolations from lower-redshift observations. It is difficult to explain this continued correlation, especially for the most massive systems, unless the most massive galaxies are forming stars near their Eddington-limited rate from their first collapse. Furthermore, we find no evidence for moderate quenching at higher masses, indicating quenching either has not occurred prior to $z \sim 6$ or else occurs rapidly, so that few galaxies are visible in transition between star-forming and quenched.

Gamma-ray bursts (GRBs) have been proposed as a tool for studying star formation in the Universe, so it is crucial to investigate whether their host galaxies and immediate environments are in any way special compared with other star-forming galaxies. Here we present spatially resolved maps of dust emission of the host galaxy of the closest known GRB 980425 at z=0.0085 using our new high-resolution observations from Herschel, Atacama Pathfinder Experiment (APEX), Atacama Large Millimeter Array (ALMA) and Australia Telescope Compact Array (ATCA). We modelled the spectral energy distributions of the host and of the star-forming region displaying the Wolf-Rayet signatures in the spectrum (WR region), located 800 pc from the GRB position. The host is characterised by low dust content and a high fraction of UV-visible star formation, similar to other dwarf galaxies. These galaxies are abundant in the local universe, so it is not surprising to find a GRB in one of them, assuming the correspondence between the GRB rate and star formation. The WR region contributes substantially to the host emission at the far-infrared, millimetre, and radio wavelengths and we propose that this is a consequence of its high gas density. If dense environments are also found close to the positions of other GRBs, then the ISM density should also be considered, along with metallicity, an important factor influencing whether a given stellar population can produce a GRB.

Emission line diagnostic diagrams probing the ionization sources in galaxies, such as the Baldwin-Phillips-Terlevich (BPT) diagram, have been used extensively to distinguish AGN from purely star-forming galaxies. Yet, they remain poorly understood at higher redshifts. We shed light on this issue with an empirical approach based on a z~0 reference sample built from ~300,000 SDSS galaxies, from which we mimic selection effects due to typical emission line detection limits at higher redshift. We combine this low-redshift reference sample with a simple prescription for luminosity evolution of the global galaxy population to predict the loci of high-redshift galaxies on the BPT and Mass-Excitation (MEx) diagnostic diagrams. The predicted bivariate distributions agree remarkably well with direct observations of galaxies out to z~1.5, including the observed stellar mass-metallicity (MZ) relation evolution. As a result, we infer that high-redshift star-forming galaxies are consistent with having "normal" ISM properties out to z~1.5, after accounting for selection effects and line luminosity evolution. Namely, their optical line ratios and gas-phase metallicities are comparable to that of low-redshift galaxies with equivalent emission-line luminosities. In contrast, AGN narrow-line regions may show a shift toward lower metallicities at higher redshift. While a physical evolution of the ISM conditions is not ruled out for purely star-forming galaxies, and may be more important starting at z>2, we find that reliably quantifying this evolution is hindered by selections effects. The recipes provided here may serve as a basis for future studies toward this goal. Code to predict the loci of galaxies on the BPT and MEx diagnostic diagrams, and the MZ relation as a function of emission line luminosity limits, is made publicly available.

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.

By exploiting the VLA-COSMOS and the Herschel-PEP surveys, we investigate the Far Infrared (FIR) properties of radio-selected AGN. To this purpose, from VLA-COSMOS we considered the 1537, F[1.4 GHz]>0.06 mJy sources with a reliable redshift estimate, and sub-divided them into star-forming galaxies and AGN solely on the basis of their radio luminosity. The AGN sample is complete with respect to radio selection at all z<~3.5. 832 radio sources have a counterpart in the PEP catalogue. 175 are AGN. Their redshift distribution closely resembles that of the total radio-selected AGN population, and exhibits two marked peaks at z~0.9 and z~2.5. We find that the probability for a radio-selected AGN to be detected at FIR wavelengths is both a function of radio power and redshift, whereby powerful sources are more likely to be FIR emitters at earlier epochs. This is due to two distinct effects: 1) at all radio luminosities, FIR activity monotonically increases with look-back time and 2) radio activity of AGN origin is increasingly less effective at inhibiting FIR emission. Radio-selected AGN with FIR emission are preferentially located in galaxies which are smaller than those hosting FIR-inactive sources. Furthermore, at all z<~2, there seems to be a preferential (stellar) mass scale M ~[10^{10}-10^{11}] Msun which maximizes the chances for FIR emission. We find such FIR (and MIR) emission to be due to processes indistinguishable from those which power star-forming galaxies. It follows that radio emission in at least 35% of the entire AGN population is the sum of two contributions: AGN accretion and star-forming processes within the host galaxy.

We present a Herschel far-IR and sub-mm study of a sample of 120 galaxies in 28 Hickson Compact Groups. Fitting their UV to sub-mm spectral energy distributions with the model of da Cunha et al. (2008), we accurately estimate the dust masses, luminosities and temperatures of the individual galaxies. We find that nearly half of the late-type galaxies in dynamically "old" groups, those with more than 25% of early-type members and redder UV-optical colours, have also significantly lower dust-to-stellar mass ratios compared to those of actively star-forming galaxies of the same mass found both in HCGs and the field. Examining their dust-to-gas mass ratios we conclude that dust was stripped out of these systems as a result of the gravitational and hydrodynamic interactions, experienced due to previous encounters with other group members. About 40% of the early-type galaxies (mostly lenticulars), in dynamically "old" groups, display dust properties similar to those of the UV-optical red late-type galaxies. Given their stellar masses, star formation rates and UV-optical colours, we suggest that red late-type and dusty lenticular galaxies represent transition populations between blue star-forming disk galaxies and quiescent early-type ellipticals. [...ABRIDGED...] Our deep Herschel observations also allow us to detect the presence of diffuse cold intragroup dust in 4 HCGs. We also find that the fraction of 250micron emission which is located outside of the main bodies of the red late-type galaxies as well as of the dusty lenticulars is 15-20% of their integrated emission at this band. All these findings are consistent with an evolutionary scenario in which gas dissipation, shocks and turbulence in addition to tidal interactions, shape the evolution of galaxies in compact groups.

Until recently, dust emission has been detected in very few host galaxies of gamma-ray bursts (GRBHs). With Herschel, we have now observed 17 GRBHs up to redshift z~3 and detected seven of them at infrared (IR) wavelengths. This relatively high detection rate (41%) may be due to the composition of our sample which at a median redshift of 1.1 is dominated by the hosts of dark GRBs. Although the numbers are small, statistics suggest that dark GRBs are more likely to be detected in the IR than their optically-bright counterparts. Combining our IR data with optical, near-infrared, and radio data from our own datasets and from the literature, we have constructed spectral energy distributions (SEDs) which span up to 6 orders of magnitude in wavelength. By fitting the SEDs, we have obtained stellar masses, dust masses, star-formation rate (SFR), and extinctions for our sample galaxies. We find that GRBHs are galaxies that tend to have a high specfic SFR (sSFR), and like other star-forming galaxies, their ratios of dust-to-stellar mass are well correlated with sSFR. We incorporate our Herschel sample into a larger compilation of GRBHs, and compare this combined sample to SFR-weighted median stellar masses of the widest, deepest galaxy survey to date. This is done in order to establish whether or not GRBs can be used as an unbiased tracer of cosmic comoving SFR density (SFRD) in the universe. In contrast with previous results, this comparison shows that GRBHs are medium-sized galaxies with relatively high sSFRs; stellar masses and sSFRs of GRBHs as a function of redshift are similar to what is expected for star-forming galaxy populations at similar redshifts. We conclude that there is no strong evidence that GRBs are biased tracers of SFRD; thus they should be able to reliably probe the SFRD to early epochs.

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

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

[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 study the incidence of nuclear obscuration on a complete sample of 1310 AGN selected on the basis of their rest-frame 2-10 keV X-ray flux from the XMM-COSMOS survey, in the redshift range 0.3<z<3.5. We classify the AGN as obscured or un-obscured on the basis of either the optical spectral properties and the overall SED or the shape of the X-ray spectrum. The two classifications agree in about 70% of the objects, and the remaining 30% can be further subdivided into two distinct classes: at low luminosities X-ray un-obscured AGN do not always show signs of broad lines or blue/UV continuum emission in their optical spectra, most likely due to galaxy dilution effects; at high luminosities broad line AGN may have absorbed X-ray spectra, which hints at an increased incidence of small-scale (sub-parsec) dust-free obscuration. We confirm that the fraction of obscured AGN is a decreasing function of the intrinsic X-ray luminosity, while the incidence of absorption shows significant evolution only for the most luminous AGN, which appear to be more commonly obscured at higher redshift. We find no significant difference between the mean stellar masses and star formation rates of obscured and un-obscured AGN hosts. We conclude that the physical state of the medium responsible for obscuration in AGN is complex, and mainly determined by the radiation environment (nuclear luminosity) in a small region enclosed within the gravitational sphere of influence of the central black hole, but is largely insensitive to the wider scale galactic conditions.

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.

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.

We present deep 450um and 850um observations of a large, uniformly covered 394arcmin^2 area in the COSMOS field obtained with the SCUBA-2 instrument on the James Clerk Maxwell Telescope (JCMT). We achieve root-mean-square noise values of 4.13mJy at 450um and 0.80mJy at 850um. The differential and cumulative number counts are presented and compared to similar previous works. Individual point sources are identified at >3.6sigma significance, a threshold corresponding to a 3-5% sample contamination rate. We identify 78 sources at 450um and 99 at 850um, with flux densities S450=13-37mJy and S850=2-16mJy. Only 62-76% of 450um sources are 850um detected and 61-81% of 850um sources are 450um detected. The positional uncertainties at 450um are small (1-2.5") and therefore allow a precise identification of multiwavelength counterparts without reliance on detection at 24um or radio wavelengths; we find that only 44% of 450um-selected galaxies and 60% of 850um-sources have 24um or radio counterparts. 450um-selected galaxies peak at <z>=1.95+-0.19 and 850um=selected galaxies peak at <z>=2.16+-0.11. The two samples occupy similar parameter space in redshift and luminosity, while their median SED peak wavelengths differ by ~10-50um (translating to deltaTdust =8-12K, where 450um-selected galaxies are warmer). The similarities of the 450um and 850um populations, yet lack of direct overlap between them, suggests that submillimeter surveys conducted at any single far-infrared wavelength will be significantly incomplete (~>30%) at censusing infrared-luminous star formation at high-z.

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.

This is the first paper of a series aiming at investigating galaxy formation and evolution in the giant-void class of the Lemaitre-Tolman-Bondi (LTB) models that best fits current cosmological observations. Here we investigate the Luminosity Function (LF) methodology, and how its estimates would be affected by a change on the cosmological model assumed in its computation. Are the current observational constraints on the allowed Cosmology enough to yield robust LF results? We use the far-infrared source catalogues built on the observations performed with the Herschel/PACS instrument, and selected as part of the PACS evolutionary probe (PEP) survey. Schechter profiles are obtained in redshift bins up to z approximately 4, assuming comoving volumes in both the standard model, that is, Friedmann-Lemaitre-Robertson-Walker metric with a perfect fluid energy-momentum tensor, and non-homogeneous LTB dust models, parametrized to fit the current combination of results stemming from the observations of supernovae Ia, the cosmic microwave background, and baryonic acoustic oscillations. We find that the luminosity functions computed assuming both the standard model and LTB void models show in general good agreement. However, the faint-end slope in the void models shows a significant departure from the standard model up to redshift 0.4. We demonstrate that this result is not artificially caused by the used LF estimator which turns out to be robust under the differences in matter-energy density profiles of the models. The differences found in the LF slopes at the faint end are due to variation in the luminosities of the sources, which depend on the geometrical part of the model. It follows that either the standard model is over-estimating the number density of faint sources or the void models are under-estimating it.

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 describe the far-infrared (FIR; rest-frame 8--1000\mu m) properties of a sample of 443 H\alpha-selected star-forming galaxies in the COSMOS and UDS fields detected by the HiZELS imaging survey. Sources are identified using narrow-band filters in combination with broad-band photometry to uniformly select H\alpha\ (and [OII] if available) emitters in a narrow redshift slice at z = 1.47+/-0.02. We use a stacking approach in Spitzer, Herschel (from PEP and HerMES surveys) and AzTEC images to describe their typical FIR properties. We find that HiZELS galaxies with observed H\alpha\ luminosities of ~ 10^{8.1-9.1} Lo have bolometric FIR luminosities of typical LIRGs, L_FIR ~ 10^{11.48+/-0.05} Lo. Combining the H\alpha\ and FIR luminosities, we derive median SFR = 32+/-5 Mo/yr and H\alpha\ extinctions of A(H\alpha) = 1.0+/-0.2 mag. Perhaps surprisingly, little difference is seen in typical HiZELS extinction levels compared to local star-forming galaxies. We confirm previous empirical stellar mass (M*) to A(H\alpha) relations and the little or no evolution up to z = 1.47. For HiZELS galaxies, we provide an empirical parametrisation of the SFR as a function of (u-z)_rest colours and 3.6\mu m photometry. We find that the observed H\alpha\ luminosity is a dominant SFR tracer when (u-z)_rest ~< 0.9 mag or when 3.6\mu m photometry > 22 mag (Vega) or when M* < 10^9.7 Mo. We do not find any correlation between the [OII]/H\alpha\ and FIR luminosity, suggesting that this emission line ratio does not trace the extinction of the most obscured star-forming regions. The luminosity-limited HiZELS sample tends to lie above of the so-called `main sequence' for star-forming galaxies, especially at low M*. This work suggests that obscured star formation is linked to the assembly of M*, with deeper potential wells in massive galaxies providing dense, heavily obscured environments in which stars can form rapidly.

Long Gamma Ray Bursts (LGRBs) are related to the final stages of evolution of massive stars. As such, they should follow the star formation rate (SFR) of galaxies. We can use them to probe for star-forming galaxies (SFGs) in the distant universe following this assumption. The relation between the rate of LGRBs in a given galaxy and its SFR (that we call the LGRB "bias") may be complex, as we have indications that the LGRB hosts are not perfect analogues to the general population of SFGs. In this work, we try to quantify the dependence of the LGRB bias on physical parameters of their host galaxy such as the SFR or the stellar mass. We propose an empirical method based on the comparison of stellar mass functions (and SFR distributions) of LGRB hosts and of SFGs in order to find how the bias depends on the stellar mass or the SFR. By applying this method to a sample of LGRB hosts at redshifts lower than 1.1, where the properties of SFGs are well established, and where the properties of LGRB host galaxies can be deduced from observations (for stellar masses larger than 10**9.25 Msun and SFR larger than 1.8 Msun / yr), we find that the LGRB bias depends on both the stellar mass and SFR. We find that the bias decreases with the SFR, i.e. we see no preference for highly SFGs, once taken into account the larger number of massive stars in galaxies with larger SFR. We do not find any trend with the specific star formation rate (SSFR) but the dynamical range in SSFR in our study is narrow. Although through an indirect method, we relate these trends to a possible decrease of the LGRBs rate / SFR ratio with the metallicity. The method we propose suggests trends that may be useful to constrain models of LGRB progenitors, showing a clear decrease of the LGRB bias with the metallicity. This is promising for the future as the number of LGRB hosts studied will increase.

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 a measurement of the average supermassive black hole accretion rate (BHAR) as a function of star formation rate (SFR) for galaxies in the redshift range 0.25<z<0.8. We study a sample of 1,767 far-IR selected star-forming galaxies in the 9 deg^2 Bo\"otes multiwavelength survey field. The SFR is estimated using 250 micron observations from the Herschel Space Observatory, for which the contribution from the AGN is minimal. In this sample, 121 AGNs are directly identified using X-ray or mid-IR selection criteria. We combined these detected AGNs and an X-ray stacking analysis for undetected sources to study the average BHAR for all of the star-forming galaxies in our sample. We find an almost linear relation between the average BHAR (in M_sun/year) and the SFR (in M_sun/year) for galaxies across a wide SFR range 0.85<log SFR<2.56 : log BHAR=(-3.72\pm0.52)+(1.05\pm0.33) log SFR. This global correlation between SFR and average BHAR is consistent with a simple picture in which SFR and AGN activity are tightly linked over galaxy evolution timescales.

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

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)

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

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

(abridged) Far-infrared Herschel photometry from the PEP and HerMES programs is combined with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields. Based on this rich dataset, we reproduce the restframe UV to FIR ten-colors distribution of galaxies using a superposition of multi-variate Gaussian modes. The median SED of each mode is then fitted with a modified version of the MAGPHYS code that combines stellar light, emission from dust heated by stars and a possible warm dust contribution heated by an AGN. The defined Gaussian grouping is also used to identify rare sources. The zoology of outliers includes Herschel-detected ellipticals, very blue z~1 Ly-break galaxies, quiescent spirals, and torus-dominated AGN with star formation. Out of these groups and outliers, a new template library is assembled, consisting of 32 SEDs describing the intrinsic scatter in the restframe UV-to-submm colors of infrared galaxies. This library is tested against L(IR) estimates with and without Herschel data included, and compared to eight other popular methods often adopted in the literature. When implementing Herschel photometry, these approaches produce L(IR) values consistent with each other within a median absolute deviation of 10-20%, the scatter being dominated more by fine tuning of the codes, rather than by the choice of SED templates. Finally, the library is used to classify 24 micron detected sources in PEP GOODS fields. AGN appear to be distributed in the stellar mass (M*) vs. star formation rate (SFR) space along with all other galaxies, regardless of the amount of infrared luminosity they are powering, with the tendency to lie on the high SFR side of the "main sequence". The incidence of warmer star-forming sources grows for objects with higher specific star formation rates (sSFR), and they tend to populate the "off-sequence" region of the M*-SFR-z space.

We measure the clustering of Extremely Red Objects (EROs) in ~8 deg^2 of the NOAO Deep Wide Field Survey Bo\"otes field in order to establish robust links between ERO z~1.2 and local galaxy z<0.1 populations. Three different color selection criteria from the literature are analyzed to assess the consequences of using different criteria for selecting EROs. Specifically, our samples are (R-K_s)>5.0 (28,724 galaxies), (I-K_s)>4.0 (22,451 galaxies) and (I-[3.6])>5.0 (64,370 galaxies). Magnitude-limited samples show the correlation length (r_0) to increase for more luminous EROs, implying a correlation with stellar mass. We can separate star-forming and passive ERO populations using the (K_s-[24]) and ([3.6]-[24]) colors to K_s=18.4 and [3.6]=17.5, respectively. Star-forming and passive EROs in magnitude limited samples have different clustering properties and host dark halo masses, and cannot be simply understood as a single population. Based on the clustering, we find that bright passive EROs are the likely progenitors of >4L^* elliptical galaxies. Bright EROs with ongoing star formation were found to occupy denser environments than star-forming galaxies in the local Universe, making these the likely progenitors of >L^* local ellipticals. This suggests that the progenitors of massive >4L^* local ellipticals had stopped forming stars by z>1.2, but that the progenitors of less massive ellipticals (down to L^*) can still show significant star formation at this epoch.

We present the mid-infrared spectra of seven of the most powerful radio-galaxies known to exist at 1.5 < z < 2.6. The radio emission of these sources is dominated by the AGN with 500 MHz luminosities in the range 10^27.8 - 10^29.1 W/Hz. The AGN signature is clearly evident in the mid-infrared spectra, however, we also detect polycyclic aromatic hydrocarbons emission, indicative of prodigious star formation at a rate of up to ~1000 Msun/yr. Interestingly, we observe no significant correlation between AGN power and star formation in the host galaxy. We also find most of these radio galaxies to have weak 9.7 um silicate absorption features (tau_{9.7 um} < 0.8) which implies that their mid-IR obscuration is predominantly due to the dusty torus that surrounds the central engine, rather than the host galaxy. The tori are likely to have an inhomogeneous distribution with the obscuring structure consisting of individual clouds. We estimate that these radio galaxies have already formed the bulk of their stellar mass and appear to lie at a stage in their evolution where the obscured AGN dominates the energy output of the system but star formation is also prevalent.

Many theoretical models require powerful active galactic nuclei (AGNs) to suppress star formation in distant galaxies and reproduce the observed properties of today's massive galaxies. A recent study based on Herschel-SPIRE submillimetre observations claimed to provide direct support for this picture, reporting a significant decrease in the mean star-formation rates (SFRs) of the most luminous AGNs (Lx>10^44 erg/s) at z=1-3 in the Chandra Deep Field-North (CDF-N). In this letter we extend these results using Herschel-SPIRE 250um data in the COSMOS and CDF-S fields to achieve an order of magnitude improvement in the number of sources at Lx>10^44 erg/s. On the basis of our analysis, we find no strong evidence for suppressed star formation in Lx>10^44 erg/s AGNs at z=1-3. The mean SFRs of the AGNs are constant over the broad X-ray luminosity range of Lx~10^43-10^45 erg/s (with mean SFRs consistent with typical star-forming galaxies at z~2; <SFRs>~100-200 Msol/yr). We suggest that the previous CDF-N results were likely due to low number statistics. We discuss our results in the context of current theoretical models.

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.

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.

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.

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.

Using the wide multi-band photometry available in the COSMOS field we explore the host galaxy properties of a large sample of Active Galactic Nuclei (AGN) obtained by combining X-ray and optical spectroscopic selections. Based on a careful study of their Spectral Energy Distribution (SED), which has been parametrized using a 2-component (AGN+galaxy) model fit, we derived dust-corrected rest-frame magnitudes, colors, stellar masses and star formation rates (SFRs). We find that AGN hosts span a large range of stellar masses and SFRs. No color-bimodality is seen at any redshift in the AGN hosts, which are found to be mainly massive, red galaxies. Once accounting for the color-mass degeneracy in well defined mass-matched samples, we find a residual marginal enhancement of AGN incidence in redder galaxies with lower specific star formation rates, and we argue that this result might emerge because of our ability to properly account for AGN light contamination and dust extinction. Interestingly, we find that the probability for a galaxy to host a black hole growing at any given "specific accretion rate" (i.e. the ratio of X-ray luminosity to the host stellar mass) is almost independent of the host galaxy mass, while it decreases as a power-law with Lx/M. By analyzing the normalization of such probability distribution, we show how the incidence of AGN increases with redshift as rapidly as (1+z)^4, in close resemblance with the overall evolution of the specific star formation rate of the entire galaxy population. Although AGN activity and star formation in galaxies do appear to have a common triggering mechanism, at least in a statistical sense, within the COSMOS sample we do not find strong evidence of any 'smoking gun' signaling powerful AGN influence on the star-forming properties of their hosts galaxies.

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.

We present the results of a Herschel-PACS study of a sample of 97 LINERs at redshift z\sim 0.3 selected from the zCOSMOS survey. Of these sources, 34 are detected in a least one PACS band, enabling reliable estimates of the far-infrared L(FIR) luminosities, and a comparison to the FIR luminosities of local LINERs. Many of our PACS-detected LINERs are also UV sources detected by Galex. Assuming that the FIR is produced in young dusty star-forming regions, the typical star-formation rates (SFRs) for the host galaxies in our sample is \sim 10 M_Sun yr-1, 1-2 orders of magnitude larger than in many local LINERs. Given stellar masses inferred from optical/NIR photometry of the (unobscured) evolved stellar populations, we find that the entire sample lies close to the star-forming "main sequence" for galaxies at redshift 0.3. For young star-forming regions, the H\alpha- and UV-based estimates of the SFRs are much smaller than the FIR-based estimates, by factors \sim 30, even assuming that all of the H\alpha\ emission is produced by O-star ionization rather than by the AGNs. These discrepancies may be due to large (and uncertain) extinctions towards the young stellar systems. Alternatively, the H\alpha\ and UV emissions could be tracing residual star-formation in an older less obscured population with decaying star formation. We also compare L(SF) and L(AGN) in local LINERs and in our sample and comment on the problematic use of several line diagnostic diagrams in cases similar to the sample under study.

We provide evidence for a correlation between the presence of giant clumps and the occurrence of active galactic nuclei (AGN) in disk galaxies. Giant clumps of 10^8-9 Msun arise from violent gravitational instability in gas-rich galaxies, and it has been proposed that this instability could feed supermassive black holes (BH). We use emission line diagnostics to compare a sample of 14 clumpy (unstable) disks and a sample of 13 smoother (stable) disks at redshift z~0.7. The majority of clumpy disks in our sample have a high probability of containing AGN. Their [OIII] emission line is strongly excited, inconsistent with low-metallicity star formation alone. [NeIII] excitation is also higher. Stable disks rarely have such properties. Stacking ultra sensitive Chandra observations (4 Ms) reveals an X-ray excess in clumpy galaxies, which cannot be solely due to star formation and confirms the presence of AGN. The clumpy galaxies in our intermediate-redshift sample have properties typical of gas-rich disk galaxies rather than mergers, being in particular on the Main Sequence of star formation. This suggests that our findings apply to the physically-similar and numerous gas-rich unstable disks at z>1. Using the observed [OIII] and X-ray luminosities, we conservatively estimate that AGN hosted by clumpy disks have typical bolometric luminosities of the order of a few 10^43 erg/s, BH growth rates ~10^-2 Msun/yr, and that these AGN are substantially obscured in X-rays. This moderate-luminosity mode could be sufficient to provide a large fraction of today's BH mass over a couple of Gyr given that our observations suggest a high duty cycle (>10%), accretion bursts with higher luminosities being possible over shorter phases. The observed evolution of disk instabilities with mass and redshift could explain the simultaneous downsizing of star formation and of BH growth.

Bolometric luminosities and Eddington ratios of both X-ray selected broad-line (Type-1) and narrow-line (Type-2) AGN from the XMM-Newton survey in the COSMOS field are presented. The sample is composed by 929 AGN (382 Type-1 AGN and 547 Type-2 AGN) and it covers a wide range of redshifts, X-ray luminosities and absorbing column densities. About 65% of the sources are spectroscopically identified as either Type-1 or Type-2 AGN (83% and 52% respectively), while accurate photometric redshifts are available for the rest of the sample. The study of such a large sample of X-ray selected AGN with a high quality multi-wavelength coverage from the far-infrared (now with the inclusion of Herschel data at 100 micron and 160 micron) to the optical-UV allows us to obtain accurate estimates of bolometric luminosities, bolometric corrections and Eddington ratios. The kbol-Lbol relations derived in the present work are calibrated for the first time against a sizable AGN sample, and rely on observed redshifts, X-ray luminosities and column density distributions. We find that kbol is significantly lower at high Lbol with respect to previous estimates by Marconi et al. (2004) and Hopkins et al. (2007). Black hole masses and Eddington ratios are available for 170 Type-1 AGN, while black hole masses for Type-2 AGN are computed for 481 objects using the black hole mass-stellar mass relation and the morphological information. We confirm a trend between kbol and lambda_Edd, with lower hard X-ray bolometric corrections at lower Eddington ratios for both Type-1 and Type-2 AGN. We find that, on average, Eddington ratio increases with redshift for all Types of AGN at any given Mbh, while no clear evolution with redshift is seen at any given Lbol.

We present infrared observations of the supernova remnant G21.5-0.9 with the Very Large Telescope, the Canada-France-Hawaii Telescope and the Spitzer Space Telescope. Using the VLT/ISAAC camera equipped with a narrow-band [FeII] 1.64um filter the entire pulsar wind nebula in SNR G21.5-0.9 was imaged. This led to detection of iron line-emitting material in the shape of a broken ring-like structure following the nebula's edge. The detected emission is limb-brightened. We also detect the compact nebula surrounding PSR J1833-1034, both through imaging with the CFHT/AOB-KIR instrument (K' band) and the IRAC camera (all bands) and also through polarimetric observations performed with VLT/ISAAC (Ks band). The emission from the compact nebula is highly polarised with an average value of the linear polarisation fraction $P_{L}^{avg} \simeq 0.47$, and the swing of the electric vector across the nebula can be observed. The infrared spectrum of the compact nebula can be described as a power law of index $\alpha_{IR} = 0.7 \pm 0.3$, and suggests that the spectrum flattens between the infrared and X-ray bands.

We present observations of CO J=2-1 line emission in infrared-luminous cluster galaxies at z~1 using the IRAM Plateau de Bure Interferometer. Our two primary targets are optically faint, dust-obscured galaxies (DOGs) found to lie within 2 Mpc of the centers of two massive (>10^14 Msun) galaxy clusters. CO line emission is not detected in either DOG. We calculate 3-sigma upper limits to the CO J=2-1 line luminosities, L'_CO < 6.08x10^9 and < 6.63x10^9 K km/s pc^2. Assuming a CO-to-H_2 conversion factor derived for ultraluminous infrared galaxies in the local Universe, this translates to limits on the cold molecular gas mass of M_H_2 < 4.86x10^9 Msun and M_H_2 < 5.30x10^9 Msun. Both DOGs exhibit mid-infrared continuum emission that follows a power-law, suggesting that an AGN contributes to the dust heating. As such, estimates of the star formation efficiencies in these DOGs are uncertain. A third cluster member with an infrared luminosity, L_IR < 7.4x10^11 Lsun, is serendipitously detected in CO J=2-1 line emission in the field of one of the DOGs located roughly two virial radii away from the cluster center. The optical spectrum of this object suggests that it is likely an obscured AGN, and the measured CO line luminosity is L'_CO = (1.94 +/- 0.35)x10^10 K km/s pc^2, which leads to an estimated cold molecular gas mass M_H_2 = (1.55+/-0.28)x10^10 Msun. A significant reservoir of molecular gas in a z~1 galaxy located away from the cluster center demonstrates that the fuel can exist to drive an increase in star-formation and AGN activity at the outskirts of high-redshift clusters.

Detailed analysis of the substructure of Lya nebulae can put important constraints on the physical mechanisms at work and the properties of galaxies forming within them. Using high resolution HST imaging of a Lya nebula at z~2.656, we have taken a census of the compact galaxies in the vicinity, used optical/near-infrared colors to select system members, and put constraints on the morphology of the spatially-extended emission. The system is characterized by (a) a population of compact, low luminosity (~0.1 L*) sources --- 17 primarily young, small (Re~1-2 kpc), disky galaxies including an obscured AGN --- that are all substantially offset (>20 kpc) from the line-emitting nebula; (b) the lack of a central galaxy at or near the peak of the Lya emission; and (c) several nearly coincident, spatially extended emission components --- Lya, HeII, and UV continuum --- that are extremely smooth. These morphological findings are difficult to reconcile with theoretical models that invoke outflows, cold flows, or resonant scattering, suggesting that while all of these physical phenomena may be occurring, they are not sufficient to explain the powering and large extent of Lya nebulae. In addition, although the compact galaxies within the system are irrelevant as power sources, the region is significantly overdense relative to the field galaxy population (by at least a factor of 4). These observations provide the first estimate of the luminosity function of galaxies within an individual Lya nebula system, and suggest that large Lya nebulae may be the seeds of galaxy groups or low-mass clusters.

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.

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 present a characterization of the close environment of GRB980425 based on 5-160mic spectro-imaging obtained with Spitzer. The Gamma-Ray Burst GRB980425 occurred in a nearby (z=0.0085) SBc-type dwarf galaxy, at a projected distance of 900pc from an HII region with strong signatures of Wolf-Rayet (WR) stars. While this "WR region" produces less than 5% of the B-band emission of the host, we find that it is responsible for 45+/-10% of the total infrared luminosity, with a maximum contribution reaching 75% at 25-30mic. This atypical property is rarely observed among morphologically-relaxed dwarves, suggesting a strong causal link with the GRB event. The luminosity of the WR region (L_8-1000mic=4.6x10^8 Lsol), the peak of its spectral energy distribution at <~100mic and the presence of highly-ionized emission lines (e.g., [NeIII]) also reveal extremely young (<5Myr) star-forming activity, with a typical time-scale of only 47Myr to double the stellar mass already built. Finally, the mid-IR over B-band luminosity ratio in this region is substantially higher than in star-forming galaxies with similar L_IR, but it is lower than in young dust-enshrouded stellar clusters. Considering the modest obscuration measured from the silicate features (tau_9.7mic ~ 0.015), this suggests that the WR region is dominated by one or several star clusters that have either partly escaped or cleared out their parent molecular cloud. Combined with the properties characterizing the whole population of GRB hosts, our results reinforce the idea that long GRBs mostly happen within or in the vicinity of relatively unobscured galactic regions harboring very recent star formation.

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.

In order to constrain the bolometric luminosities, dust properties and molecular gas content of giant Lyman alpha nebulae, the so-called Lyman alpha blobs, we have carried out a study of dust continuum and CO line emission in two well-studied representatives of this population at z ~ 3: a Lya blob discovered by its strong Spitzer MIPS 24um detection (LABd05; Dey et al. 2005) and the Steidel blob 1 (SSA22-LAB01; Steidel et al. 2000). We find that the spectral energy distribution of LABd05 is well described by an AGN-starburst composite template with L(FIR) = (4.0 +/- 0.5) x 10^12 Lsun, comparable to high-z sub-millimeter galaxies and ultraluminous infrared galaxies. New APEX/LABOCA 870um measurements rule out the reported SCUBA detection of the SSA22-LAB01 (S[850um] = 16.8 mJy) at the > 4sigma level. Consistent with this, ultra-deep Plateau de Bure Interferometer (PdBI) observations with ~2arcsec spatial resolution also fail to detect any 1.2mm continuum source down to ~0.45mJy per beam (3sigma). Combined with the existing (sub)mm observations in the literature, we conclude that the FIR luminosity of SSA22-LAB01 remains uncertain. No CO line is detected in either case down to integrated flux limits of (Snu dV) < 0.25--1.0 Jy km/s, indicating a modest molecular gas reservoir, M(H_2) < 1--3 x 10^10 Msun. The non-detections exclude, with high significance (12 sigma), the previous tentative detection of a CO(4-3) line in the SSA22-LAB01. The increased sensitivity afforded by ALMA will be critical in studying molecular gas and dust in these interesting systems.

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

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.

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.

We present a comprehensive study on the impact of the environment of compact galaxy groups on the evolution of their members using a multi-wavelength analysis, from the UV to the infrared, for a sample of 32 Hickson compact groups (HCGs) containing 135 galaxies. Fitting the SEDs of all galaxies with the state-of-the-art model of da Cunha (2008) we can accurately calculate their mass, SFR, and extinction, as well as estimate their infrared luminosity and dust content. We compare our findings with samples of field galaxies, early-stage interacting pairs, and cluster galaxies with similar data. We find that classifying the groups as dynamically "old" or "young", depending on whether or not at least one quarter of their members are early-type systems, is physical and consistent with past classifications of HCGs based on their atomic gas content. [...ABRIDGED...] We also examine their SF properties, UV-optical and mid-IR colors, and we conclude that all the evidence point to an evolutionary scenario in which the effects of the group environment and the properties of the galaxy members are not instantaneous. Early on, the influence of close companions to group galaxies is similar to the one of galaxy pairs in the field. However, as the time progresses, the effects of tidal torques and minor merging, shape the morphology and star formation history of the group galaxies, leading to an increase of the fraction of early-type members and a rapid built up of the stellar mass in the remaining late-type galaxies.

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.

We present Hubble Space Telescope (HST) imaging of 22 ultra-luminous infrared galaxies (ULIRGs) at z~2 with extremely red R-[24] colors (called dust-obscured galaxies, or DOGs) which have a local maximum in their spectral energy distribution (SED) at rest-frame 1.6um associated with stellar emission. These sources, which we call "bump DOGs", have star-formation rates of 400-4000 Msun/yr and have redshifts derived from mid-IR spectra which show strong polycyclic aromatic hydrocarbon emission --- a sign of vigorous on-going star-formation. Using a uniform morphological analysis, we look for quantifiable differences between bump DOGs, power-law DOGs (Spitzer-selected ULIRGs with mid-IR SEDs dominated by a power-law and spectral features that are more typical of obscured active galactic nuclei than starbursts), sub-millimeter selected galaxies (SMGs), and other less-reddened ULIRGs from the Spitzer extragalactic First Look Survey (XFLS). Bump DOGs are larger than power-law DOGs (median Petrosian radius of 8.4 +/- 2.7 kpc vs. 5.5 +/- 2.3 kpc) and exhibit more diffuse and irregular morphologies (median M_20 of -1.08 +/- 0.05 vs. -1.48 +/- 0.05). These trends are qualitatively consistent with expectations from simulations of major mergers in which merging systems during the peak star-formation rate period evolve from M_20 = -1.0 to M_20 = -1.7. Less obscured ULIRGs (i.e., non-DOGs) tend to have more regular, centrally peaked, single-object morphologies rather than diffuse and irregular morphologies. This distinction in morphologies may imply that less obscured ULIRGs sample the merger near the end of the peak star-formation rate period. Alternatively, it may indicate that the intense star-formation in these less-obscured ULIRGs is not the result of a recent major merger.

Using data from the Spitzer Space Telescope, we analyze the mid-infrared (3-70 micron) spectral energy distributions of dry merger candidates in the Bootes field of the NOAO Deep Wide-Field Survey. These candidates were selected by previous authors to be luminous, red, early-type galaxies with morphological evidence of recent tidal interactions. We find that a significant fraction of these candidates exhibit 8 and 24 micron excesses compared to expectations for old stellar populations. We estimate that a quarter of dry merger candidates have mid-infrared-derived star formation rates greater than ~1 MSun/yr. This represents a "frosting" on top of a large old stellar population, and has been seen in previous studies of elliptical galaxies. Further, the dry merger candidates include a higher fraction of starforming galaxies relative to a control sample without tidal features. We therefore conclude that the star formation in these massive ellipticals is likely triggered by merger activity. Our data suggest that the mergers responsible for the observed tidal features were not completely dry, and may be minor mergers involving a gas-rich dwarf galaxy.

We present the first measurement of the spatial clustering of mid-infrared selected obscured and unobscured quasars, using a sample in the redshift range 0.7 < z < 1.8 selected from the 9 deg^2 Bootes multiwavelength survey. Recently the Spitzer Space Telescope and X-ray observations have revealed large populations of obscured quasars that have been inferred from models of the X-ray background and supermassive black hole evolution. To date, little is known about obscured quasar clustering, which allows us to measure the masses of their host dark matter halos and explore their role in the cosmic evolution of black holes and galaxies. In this study we use a sample of 806 mid-infrared selected quasars and ~250,000 galaxies to calculate the projected quasar-galaxy cross-correlation function w_p(R). The observed clustering yields characteristic dark matter halo masses of log (M_halo [h^-1 M_sun]) = 12.7^+0.4_-0.6 and 13.3^+0.3_-0.4 for unobscured quasars (QSO-1s) and obscured quasars (Obs-QSOs), respectively. The results for QSO-1s are in excellent agreement with previous measurements for optically-selected quasars, while we conclude that the Obs-QSOs are at least as strongly clustered as the QSO-1s. We test for the effects of photometric redshift errors on the optically-faint Obs-QSOs, and find that our method yields a robust lower limit on the clustering; photo-z errors may cause us to underestimate the clustering amplitude of the Obs-QSOs by at most ~20%. We compare our results to previous studies, and speculate on physical implications of stronger clustering for obscured quasars.

By cross-correlating AKARI infrared (IR) sources with the SDSS galaxies, we identified 2357 infrared galaxies with a spectroscopic redshift. This is not just one of the largest samples of local IR galaxies, but AKARI provides crucial FIR bands (9, 18, 65, 90, 140, and 160um) in accurately measuring galaxy SED across the peak of the dust emission at ~100um. By fitting modern IR SED models to the AKARI photometry, we measured the total infrared luminosity (L_IR) of individual galaxies more accurately. Using this L_IR, we constructed luminosity functions of infrared galaxies at a median redshift of z=0.031, with 4 times larger sample than previous work. The LF agrees well with that at z=0.0082 (RBGS), showing smooth and continuous evolution toward higher redshift LFs measured in the AKARI NEP deep field. The derived local cosmic IR luminosity density is Omega_IR=3.8x10^8 LsunMpc^-3. We separate galaxies into AGN, star-forming, and composite by using the [NII]/Ha vs [OIII]/Hb line ratios. The fraction of AGN shows a continuous increase with increasing L_IR from 25% to 90% at 9<log L_IR<12.5. The SFR_Ha and L_[OIII] show good correlations with L_IR for SFG (star-forming galaxies) and AGN, respectively. The self-absorption corrected Ha/Hb ratio shows a weak increase with L_IR with a substantial scatter. When we separate IR LFs into contributions from AGN and SFG, the AGN contribution becomes dominant at L_IR>10^11Lsun, coinciding the break of the both SFG and AGN IR LFs. At L_IR<10^11Lsun, SFG dominates IR Lfs. Only 1.1% of Omega_IR is produced by LIRG, and only 0.03% is by ULIRG in the local Universe. This work also provides the most accurate infrared luminosity density of the local Universe to date. Compared with high redshift results from the AKARI NEP deep survey, we observed a strong evolution of Omega_IR^SFG ~(1+z)^4.1+-0.4 and Omega_IR^AGN ~(1+z)^4.1+-0.5 (abridged).

We 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 the J-band luminosity function of 1838 mid-infrared and X-ray selected AGNs in the redshift range 0<z<5.85. These luminosity functions are constructed by combining the deep multi-wavelength broad-band observations from the UV to the mid-IR of the NDWFS Bootes field with the X-ray observations of the XBootes survey and the spectroscopic observations of the same field by AGES. Our sample is primarily composed of IRAC-selected AGNs, targeted using modifications of the Stern et al.(2005) criteria, complemented by MIPS 24 microns and X-ray selected AGNs to alleviate the biases of IRAC mid-IR selection against z~4.5 quasars and AGNs faint with respect to their hosts. This sample provides an accurate link between low and high redshift AGN luminosity functions and does not suffer from the usual incompleteness of optical samples at z~3. We find that the space density of the brightest quasars strongly decreases from z=3 to z=0, while the space density of faint quasars is at least flat, and possibly increasing, over the same redshift range. At z>3 we observe a decrease in the space density of quasars of all brightnesses. We model the luminosity function by a double power-law and find that its evolution cannot be described by either pure luminosity or pure density evolution, but must be a combination of both. Our best-fit model has bright and faint power-law indices consistent with the low redshift measurements based on the 2QZ and 2SLAQ surveys and it generally agrees with the number of bright quasars predicted by other LFs at all redshifts. If we construct the QSO luminosity function using only the IRAC-selected AGNs, we find that the biases inherent to this selection method significantly modify the behavior of phi*(z) only for z<1 and have no significant impact upon the characteristic magnitude M*_J(z).

We report a sizable class of type 1 active galactic nuclei (AGNs) with unusually weak near-infrared (1-3{\mu}m) emission in the XMM-COSMOS type 1 AGN sample. The fraction of these "hot-dust-poor" AGNs increases with redshift from 6% at lowredshift (z < 2) to 20% at moderate high redshift (2 < z < 3.5). There is no clear trend of the fraction with other parameters: bolometric luminosity, Eddington ratio, black hole mass, and X-ray luminosity. The 3{\mu}m emission relative to the 1{\mu}m emission is a factor of 2-4 smaller than the typical Elvis et al. AGN spectral energy distribution (SED), which indicates a "torus" covering factor of 2%-29%, a factor of 3-40 smaller than required by unified models. The weak hot dust emission seems to expose an extension of the accretion disk continuum in some of the source SEDs. We estimate the outer edge of their accretion disks to lie at $(0.3-2.0) /times 10^4$ Schwarzschild radii, ~10-23 times the gravitational stability radii. Formation scenarios for these sources are discussed.

We present Spitzer 16 micron imaging of the Great Observatories Origins Deep Survey (GOODS) fields. We survey 150 square arcminutes in each of the two GOODS fields (North and South), to an average 3 sigma depth of 40 and 65 micro-Jy respectively. We detect about 1300 sources in both fields combined. We validate the photometry using the 3-24 micron spectral energy distribution of stars in the fields compared to Spitzer spectroscopic templates. Comparison with ISOCAM and AKARI observations in the same fields show reasonable agreement, though the uncertainties are large. We provide a catalog of photometry, with sources cross correlated with available Spitzer, Chandra, and HST data. Galaxy number counts show good agreement with previous results from ISOCAM and AKARI, with improved uncertainties. We examine the 16 to 24 micron flux ratio and find that for most sources it lies within the expected locus for starbursts and infrared luminous galaxies. A color cut of S_{16}/S_{24}>1.4 selects mostly sources which lie at 1.1<z<1.6, where the 24 micron passband contains both the redshifted 9.7 micron silicate absorption and the minimum between PAH emission peaks. We measure the integrated galaxy light of 16 micron sources, and find a lower limit on the galaxy contribution to the extragalactic background light at this wavelength to be 2.2\pm 0.2$ nW m^{-2} sr^{-1}.

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

We present Spitzer 7-38um spectra for a 24um flux limited sample of galaxies at z~0.7 in the COSMOS field. The detailed high-quality spectra allow us to cleanly separate star formation (SF) and active galactic nucleus (AGN) in individual galaxies. We first decompose mid-infrared Luminosity Functions (LFs). We find that the SF 8um and 15um LFs are well described by Schechter functions. AGNs dominate the space density at high luminosities, which leads to the shallow bright-end slope of the overall mid-infrared LFs. The total infrared (8-1000um) LF from 70um selected galaxies shows a shallower bright-end slope than the bolometrically corrected SF 15um LF, owing to the intrinsic dispersion in the mid-to-far-infrared spectral energy distributions. We then study the contemporary growth of galaxies and their supermassive black holes (BHs). Seven of the 31 Luminous Infrared Galaxies with Spitzer spectra host luminous AGNs, implying an AGN duty cycle of 23+/-9%. The time-averaged ratio of BH accretion rate and SF rate matches the local M_BH-M_bulge relation and the M_BH-M_host relation at z ~ 1. These results favor co-evolution scenarios in which BH growth and intense SF happen in the same event but the former spans a shorter lifetime than the latter. Finally, we compare our mid-infrared spectroscopic selection with other AGN identification methods and discuss candidate Compton-thick AGNs in the sample. While only half of the mid-infrared spectroscopically selected AGNs are detected in X-ray, ~90% of them can be identified with their near-infrared spectral indices.

Hickson Compact Groups (HCGs) are among the densest galaxy environments of the local universe. To examine the effects of the environment on the infrared properties of these systems, we present an analysis of Spitzer and ISO mid-infrared imaging as well as deep ground based near-infrared imaging of 14 HCGs containing a total of 69 galaxies. Based on mid-infrared color diagnostics we identify the galaxies which appear to host an active nucleus, while using a suite of templates, we fit the complete infrared spectral energy distribution for each group member. We compare our estimates of galaxy mass, star formation rate, total infrared luminosities, and specific star formation rates (sSFR) for our HCG sample, to samples of isolated galaxies and interacting pairs and find that overall there is no discernible difference among them. However, HCGs which can be considered as dynamically "old", host late-type galaxies with a slightly lower sSFR than the one found in dynamically "young" groups. This could be attributed to multiple past interactions among the galaxies in old groups, that have led to the build up of their stellar mass. It is also consistent with our prediction for the presence of diffuse cold dust in the intergalactic medium of 9 of the dynamically "old" groups.

Infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution. The AKARI IR space telescope performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160um) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, AKARI can much more precisely measure the total infrared luminosity (L_TIR) of individual galaxies, and thus, the total infrared luminosity density in the local Universe. By fitting IR SED models, we have re-measured L_TIR of the IRAS Revised Bright Galaxy Sample. We present mid-IR monochromatic luminosity to L_TIR conversions for Spitzer 8,24um, AKARI 9,18um, IRAS 12um, WISE 12,22um, and ISO 15um filters, with scatter ranging 13-44%. The resulting AKARI IR luminosity function (LF) agrees well with that from the IRAS. We integrate the LF weighted by L_TIR to obtain a cosmic IR luminosity density of Omega_TIR= (8.5^{+1.5}_{-2.3})x 10^7 L Mpc^-3, of which 7+-1% is produced by LIRGs, and only 0.4+-0.1% is from ULIRGs in the local Universe. Once IR contributions from AGN and star-forming galaxies (SFG) are separated, SFG IR LF shows a steep decline at the bright-end. Compared with high-redshift results from the AKARI NEP deep survey, these data show a strong evolution of Omega_TIRSF propto (1+z)^4.0+-0.5, and Omega_TIRAGN propto (1+z)^4.4+-0.4. For Omega_TIRAGN, the ULIRG contribution exceeds that from LIRG already by z~1. A rapid evolution in both Omega_TIRAGN and Omega_TIRSFG suggests the correlation between star formation and black hole accretion rate continues up to higher redshifts. We compare the evolution of Omega_TIRAGN to that of X-ray luminosity density. The Omega_TIRAGN/Omega_X-rayAGN ratio shows a possible increase at z>1, suggesting an increase of obscured AGN at z>1.

We present a search for Herschel-PACS counterparts of dust-obscured, high-redshift objects previously selected at submillimeter and millimeter wavelengths in the Great Observatories Origins Deep Survey North field. We detect 22 of 56 submillimeter galaxies (39%) with a SNR of >=3 at 100 micron down to 3.0 mJy, and/or at 160 micron down to 5.7 mJy. The fraction of SMGs seen at 160 micron is higher than that at 100 micron. About 50% of radio-identified SMGs are associated with PACS sources. We find a trend between the SCUBA/PACS flux ratio and redshift, suggesting that these flux ratios could be used as a coarse redshift indicator. PACS undetected submm/mm selected sources tend to lie at higher redshifts than the PACS detected ones. A total of 12 sources (21% of our SMG sample) remain unidentified and the fact that they are blank fields at Herschel-PACS and VLA 20 cm wavelength may imply higher redshifts for them than for the average SMG population (e.g., z>3-4). The Herschel-PACS imaging of these dust-obscured starbursts at high-redshifts suggests that their far-infrared spectral energy distributions have significantly different shapes than template libraries of local infrared galaxies.

We investigate the effects of the environment on star-formation in a sample of massive luminous and ultra-luminous infrared galaxies (LIRGs and ULIRGs) with S(24 micron)>80 uJy and i+<24 in the COSMOS field. We exploit the accurate photometric redshifts in COSMOS to characterize the galaxy environment and study the evolution of the fraction of LIRGs and ULIRGs in different environments in the redshift range z=0.3-1.2 and in bins of stellar mass. We find that the environment plays a role in the star formation processes and evolution of LIRGs and ULIRGs. We find an overall increase of the ULIRG+LIRG fraction in an optically-selected sample with increasing redshift, as expected from the evolution of the star formation rate density. We find that the ULIRG+LIRG fraction decreases with increasing density up to z~1, and that the dependence on density flattens with increasing redshift. We do not observe the reversal of the star-formation rate density relation up to z=1 in massive LIRGs and ULIRGs, suggesting that such reversal might occur at higher redshift in this infrared luminosity range.

We investigate the (large-scale) bar fraction in a mass-complete sample of M > 10^10.5 Msun disk galaxies at 0.2 < z < 0.6 in the COSMOS field. The fraction of barred disks strongly depends on mass, disk morphology, and specific star formation rate (SSFR). At intermediate stellar mass (10^10.5 < M < 10^11 Msun) the bar fraction in early-type disks is much higher, at all redshifts, by a factor ~2, than that in late-type disks. This trend is reversed at higher stellar mass (M > 10^11 Msun), where the fraction of bars in early-type disks becomes significantly lower, at all redshifts, than that in late-type disks. The bar fractions for galaxies with low and high SSFRs closely follow those of the morphologically-selected early-type and late-type populations, respectively. This indicates a close correspondence between morphology and SSFR in disk galaxies at these earlier epochs. Interestingly, the total bar fraction in 10^10.5 < M < 10^11 Msun disks is built up by a factor of ~2 over the redshift interval explored, while for M > 10^11 Msun disks it remains roughly constant. This indicates that, already by z ~ 0.6, spectral and morphological transformations in the most massive disk galaxies have largely converged to the familiar Hubble sequence that we observe in the local Universe, while for intermediate mass disks this convergence is ongoing until at least z ~ 0.2. Moreover, these results highlight the importance of employing mass-limited samples for quantifying the evolution of barred galaxies. Finally, the evolution of the barred galaxy populations investigated does not depend on the large-scale environmental density (at least, on the scales which can be probed with the available photometric redshifts).

We analyze the morphological properties of a large sample of 1503 70 micron selected galaxies in the COSMOS field spanning the redshift range 0.01<z< 3.5 with a median redshift of 0.5 and an infrared luminosity range of 10^8<L_IR<10^14L_sun with a median luminosity of 10^11.4 L_sun. In general these galaxies are massive, with a stellar mass range of 10^10-10^12 M_sun, and luminous, with -25<M_K<-20. We find a strong correlation between the fraction of major mergers and L_IR, with the fraction at the highest luminosity being up to 50%. We also find that the fraction of spirals drops dramatically with L_IR. Minor mergers likely play a role in boosting the infrared luminosity for sources with low luminosities. The precise fraction of mergers in any given L_IR bin varies by redshift due to sources at z>1 being difficult to classify and subject to the effects of band pass shifting, therefore, these numbers can only be considered lower limits. At z<1, where the morphological classifications are most robust, major mergers clearly dominate the ULIRG population (50-80%) and are important for the LIRG population (25-40%). At z>1 the fraction of major mergers is at least 30-40% for ULIRGs. Although the general morphological trends agree with what has been observed for local (U)LIRGs, the fraction of major mergers is slightly lower than seen locally. This is in part due to the difficulty of identifying merger signatures at high redshift. We argue that given the number of major gas-rich mergers observed and the relatively short timescale that they would be observable in the (U)LIRG phase that it is plausible for the observed red sequence of massive ellipticals (<10^12 M_sun) to have been formed entirely by gas-rich major mergers.

We present the 24 micron rest-frame luminosity function (LF) of star-forming galaxies in the redshift range 0.0 < z < 0.6 constructed from 4047 spectroscopic redshifts from the AGN and Galaxy Evolution Survey of 24 micron selected sources in the Bootes field of the NOAO Deep Wide-Field Survey. This sample provides the best available combination of large area (9 deg^2), depth, and statistically complete spectroscopic observations, allowing us to probe the evolution of the 24 micron LF of galaxies at low and intermediate redshifts while minimizing the effects of cosmic variance. In order to use the observed 24 micron luminosity as a tracer for star formation, active galactic nuclei (AGNs) that could contribute significantly at 24 micron are identified and excluded from our star-forming galaxy sample based on their mid-IR spectral energy distributions or the detection of X-ray emission. The evolution of the 24 micron LF of star-forming galaxies for redshifts of z < 0.65 is consistent with a pure luminosity evolution where the characteristic 24 micron luminosity evolves as (1+z)^(3.8+/-0.3). We extend our evolutionary study to encompass 0.0 < z < 1.2 by combining our data with that of the Far-Infrared Deep Extragalactic Legacy Survey. Over this entire redshift range the evolution of the characteristic 24 micron luminosity is described by a slightly shallower power law of (1+z)^(3.4+/-0.2). We find a local star formation rate density of (1.09+/-0.21) x 10^-2 Msun/yr/Mpc^-3, and that it evolves as (1+z)^(3.5+/-0.2) over 0.0 < z < 1.2. These estimates are in good agreement with the rates using optical and UV fluxes corrected for the effects of intrinsic extinction in the observed sources. This agreement confirms that star formation at z <~ 1.2 is robustly traced by 24 micron observations and that it largely occurs in obscured regions of galaxies. (Abridged)

We present ultra-deep mid-IR spectra of 48 infrared-luminous galaxies in the GOODS-South field obtained with the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. These galaxies are selected among faint infrared sources (0.14 - 0.5 mJy at 24 um) in two redshift bins (0.76-1.05 and 1.75-2.4) to sample the major contributors to the cosmic infrared background at the most active epochs. We estimate redshifts for 92% of the sample using PAH and Si absorption features. Only few of these galaxies (5% at z~1 and 12% at z~2) have their total infrared luminosity dominated by emission from AGN. The averaged mid-IR spectra of the z~1 LIRGs and of the z~2 ULIRGs are very similar to the averaged spectrum of local starbursts and HII-like ULIRGs, respectively. We find that 6.2um PAH equivalent widths reach a plateau of ~1 um for L(24 mu) < 1E11 L(sun). At higher luminosities, EW (6.2 mu) anti-correlates with L(24 um). Intriguingly, high-z ULIRGs and SMG lie above the local EW (6.2 um) - L(24 um) relationship suggesting that, at a given luminosity, high-z ULIRGs have AGN contributions to their dust emission lower than those of local counterparts. A quantitative analysis of their morphology shows that most of the luminous IR galaxies have morphologies similar to those of IR-quiet galaxies at the same redshift. All z~2 ULIRGs of our sample are IR-excess BzK galaxies and most of them have L(FIR)/L(1600A) ratios higher than those of starburst galaxies at a given UV slope. The ``IR excess'' (Daddi et al. 2007) is mostly due to strong 7.7 um PAH emission and under-estimation of UV dust extinction. On the basis of the AGN-powered L (6 um) continuum measured directly from the mid-IR spectra, we estimate an average intrinsic X-ray AGN luminosity of L(2-10 keV) = (0.1 +/- 0.6) 1E43 erg/s, a value substantially lower than the prediction by Daddi et al. (2007).

We present a study of the average properties of luminous infrared galaxies detected directly at 24 $\mu$m in the COSMOS field using a median stacking analysis at 70$\mu$m and 160 $\mu$m. Over 35000 sources spanning 0<z<3 and 0.06 mJy<S_{24}<3.0 mJy are stacked, divided into bins of both photometric redshift and 24 $\mu$m flux. We find no correlation of $S_{70}/S_{24}$ flux density ratio with $S_{24}$, but find that galaxies with higher $S_{24}$ have a lower $S_{160}/S_{24}$ flux density ratio. These observed ratios suggest that 24 $\mu$m selected galaxies have warmer SEDs at higher mid-IR fluxes, and therefore have a possible higher fraction of AGN. Comparisons of the average $S_{70}/S_{24}$ and $S_{160}/S_{24}$ colors with various empirical templates and theoretical models show that the galaxies detected at 24 $\mu$m are consistent with "normal" star-forming galaxies and warm mid-IR galaxies such as Mrk 231, but inconsistent with heavily obscured galaxies such as Arp 220. We perform a $\chi^{2}$ analysis to determine best fit galactic model SEDs and total IR luminosities for each of our bins. We compare our results to previous methods of estimating $L_{\rm{IR}}$ and find that previous methods show considerable agreement over the full redshift range, except for the brightest $S_{24}$ sources, where previous methods overpredict the bolometric IR luminosity at high redshift, most likely due to their warmer dust SED. We present a table that can be used as a more accurate and robust method for estimating bolometric infrared luminosity from 24 $\mu$m flux densities.

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

We exploit the deepest existing far-infrared (FIR) data obtained so far by Herschel at 100 and 160 um in the GOODS-N, as part of the PACS Evolutionary Probe (PEP) survey, to derive for the first time the evolution of the rest-frame 60-um, 90-um, and total IR luminosity functions (LFs) of galaxies and AGNs from z=0 to unprecedented high redshifts (z~2-3). The PEP LFs were computed using the 1/Vmax method. The FIR sources were classified by means of a detailed broad- band SED-fitting analysis and spectral characterisation. Based on the best-fit model results, k-correction and total IR (8-1000 um) luminosity were obtained for each source. LFs (monochromatic and total) were then derived for various IR populations separately in different redshift bins and compared to backward evolution model predictions. We detect strong evolution in the LF to at least z~2. Objects with SEDs similar to local spiral galaxies are the major contributors to the star formation density (SFD) at z< 0.3, then, as redshift increases, moderate SF galaxies - most likely containing a low-luminosity AGN - start dominating up to z ~= 1.5. At >1.5 the SFD is dominated by the contributions of starburst galaxies. In agreement with previous findings, the comoving IR LD derived from our data evolves approximately as (1 + z)^(3.8+/-0.3) up to z~1, there being some evidence of flattening up to z~2.

We report the final optical identifications of the medium-depth (~60 ksec), contiguous (2 deg^2) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ~800 X-ray sources down to limiting fluxes of ~5x10^{-16}, ~3x10^{-15}, and ~7x10^{-15} erg/cm2/s in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. The work is complemented by an extensive collection of multi-wavelength data from 24 micron to UV, available from the COSMOS survey, for each of the X-ray sources, including spectroscopic redshifts for ~50% of the sample, and high-quality photometric redshifts for the rest. The XMM and multiwavelength flux limits are well matched: 1760 (98%) of the X-ray sources have optical counterparts, 1711 (~95%) have IRAC counterparts, and 1394 (~78%) have MIPS 24micron detections. Thanks to the redshift completeness (almost 100%) we were able to constrain the high-luminosity tail of the X-ray luminosity function confirming that the peak of the number density of logL_X>44.5 AGN is at z~2. Spectroscopically-identified obscured and unobscured AGN, as well as normal and starforming galaxies, present well-defined optical and infrared properties. We devised a robust method to identify a sample of ~150 high redshift (z>1), obscured AGN candidates for which optical spectroscopy is not available. We were able to determine that the fraction of the obscured AGN population at the highest (L_X>10^{44} erg s^{-1}) X-ray luminosity is ~15-30% when selection effects are taken into account, providing an important observational constraint for X-ray background synthesis. We studied in detail the optical spectrum and the overall spectral energy distribution of a prototypical Type 2 QSO, caught in a stage transitioning from being starburst dominated to AGN dominated, which was possible to isolate only thanks to the combination of X-ray and infrared observations.

Previous observational studies of the infrared (IR)-radio relation out to high redshift employed any detectable star forming systems at a given redshift within the restricted area of cosmological survey fields. Consequently, the evolution inferred relies on a comparison between the average IR/radio properties of (i) very IR-luminous high-z sources and (ii) more heterogeneous low(er)-z samples that often lack the strongest IR emitters. In this report we consider populations of objects with comparable luminosities over the last 10 Gyr by taking advantage of deep IR (esp. Spitzer 24 micron) and VLA 1.4 GHz observations of the COSMOS field. Consistent with recent model predictions, both Ultra Luminous Infrared Galaxies (ULIRGs) and galaxies on the bright end of the evolving IR luminosity function do not display any change in their average IR/radio ratios out to z~2 when corrected for bias. Uncorrected data suggested ~0.3 dex of positive evolution.

We present a set of low resolution empirical SED templates for AGNs and galaxies in the wavelength range from 0.03 to 30 microns based on the multi-wavelength photometric observations of the NOAO Deep-Wide Field Survey Bootes field and the spectroscopic observations of the AGN and Galaxy Evolution Survey. Our training sample is comprised of 14448 galaxies in the redshift range 0<~z<~1 and 5347 likely AGNs in the range 0<~z<~5.58. We use our templates to determine photometric redshifts for galaxies and AGNs. While they are relatively accurate for galaxies, their accuracies for AGNs are a strong function of the luminosity ratio between the AGN and galaxy components. Somewhat surprisingly, the relative luminosities of the AGN and its host are well determined even when the photometric redshift is significantly in error. We also use our templates to study the mid-IR AGN selection criteria developed by Stern et al.(2005) and Lacy et al.(2004). We find that the Stern et al.(2005) criteria suffers from significant incompleteness when there is a strong host galaxy component and at z =~ 4.5, when the broad Halpha emission line is redshifted into the [3.6] band, but that it is little contaminated by low and intermediate redshift galaxies. The Lacy et al.(2004) criterion is not affected by incompleteness at z =~ 4.5 and is somewhat less affected by strong galaxy host components, but is heavily contaminated by low redshift star forming galaxies. Finally, we use our templates to predict the color-color distribution of sources in the upcoming WISE mission and define a color criterion to select AGNs analogous to those developed for IRAC photometry. We estimate that in between 640,000 and 1,700,000 AGNs will be identified by these criteria, but will have serious completeness problems for z >~ 3.4.

We present a detailed study of a peculiar source in the COSMOS survey at z=0.359. Source CXOCJ100043.1+020637 (CID-42) presents two compact optical sources embedded in the same galaxy. The distance between the 2, measured in the HST/ACS image, is 0.495" that, at the redshift of the source, corresponds to a projected separation of 2.46 kpc. A large (~1200 km/s) velocity offset between the narrow and broad components of Hbeta has been measured in three different optical spectra from the VLT/VIMOS and Magellan/IMACS instruments. CID-42 is also the only X-ray source having in its X-ray spectra a strong redshifted broad absorption iron line, and an iron emission line, drawing an inverted P-Cygni profile. The Chandra and XMM data show that the absorption line is variable in energy by 500 eV over 4 years and that the absorber has to be highly ionized, in order not to leave a signature in the soft X-ray spectrum. That these features occur in the same source is unlikely to be a coincidence. We envisage two possible explanations: (1) a gravitational wave recoiling black hole (BH), caught 1-10 Myr after merging, (2) a Type 1/ Type 2 system in the same galaxy where the Type 1 is recoiling due to slingshot effect produced by a triple BH system. The first possibility gives us a candidate gravitational waves recoiling BH with both spectroscopic and imaging signatures. In the second case, the X-ray absorption line can be explained as a BAL-like outflow from the foreground nucleus (a Type 2 AGN) at the rearer one (a Type 1 AGN), which illuminates the otherwise undetectable wind, giving us the first opportunity to show that fast winds are present in obscured AGN.

VLA 1.4 GHz (rms noise ~0.012 mJy) and MIPS 24 and 70 micron (rms noise ~0.02 and ~1.7 mJy, respectively) observations covering the 2 square degree COSMOS field are combined with an extensive multi-wavelength data set to study the evolution of the IR-radio relation at intermediate and high redshift. With ~4500 sources -- of which ~30% have spectroscopic redshifts -- the current sample is significantly larger than previous ones used for the same purpose. Both monochromatic IR/radio flux ratios (q24 & q70), as well as the ratio of the total IR and the 1.4 GHz luminosity (qTIR) are used as indicators for the IR/radio properties of star forming galaxies and AGN. Using a sample jointly selected at IR and radio wavelengths in order to reduce selection biases, we provide firm support for previous findings that the IR-radio relation remains unchanged out to at least z~1.4. Moreover, based on data from ~150 objects we also find that the local relation likely still holds at 2.5<z<5. At redshift z<1.4 we observe that radio-quiet AGN populate the locus of the IR-radio relation in similar numbers as star forming sources. In our analysis we employ the methods of survival analysis in order to ensure a statistically sound treatment of flux limits arising from non-detections. We determine the observed shift in average IR/radio properties of IR- and radio- selected populations and show that it can reconcile apparently discrepant measurements presented in the literature. Finally, we also investigate variations of the IR/radio ratio with IR and radio luminosity and find that it hardly varies with IR luminosity but is a decreasing function of radio luminosity.

Dust-obscured star-formation becomes much more important with increasing intensity, and increasing redshift. We aim to reveal cosmic star-formation history obscured by dust using deep infrared observation with the AKARI. We construct restframe 8um, 12um, and total infrared (TIR) luminosity functions (LFs) at 0.15<z<2.2 using 4128 infrared sources in the AKARI NEP-Deep field. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24um) by the AKARI satellite allows us to estimate restframe 8um and 12um luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. We have found that all 8um (0.38<z<2.2), 12um (0.15<z<1.16), and TIR LFs (0.2<z<1.6), show a continuous and strong evolution toward higher redshift. In terms of cosmic infrared luminosity density (Omega_IR), which was obtained by integrating analytic fits to the LFs, we found a good agreement with previous work at z<1.2, and that the Omega_IR evolves as propto (1+z)^4.4+-1.0. When we separate contributions to Omega_IR by LIRGs and ULIRGs, we found more IR luminous sources are increasingly more important at higher redshift. We found that the ULIRG (LIRG) contribution increases by a factor of 10 (1.8) from z=0.35 to z=1.4.

We present a study of the X-ray to optical properties of a sample of 545 X-ray selected Type 1 AGN, from the XMM-COSMOS survey, over a wide range of redshifts ($0.04<\z<4.25$) and X-ray luminosities ($40.6 \leq \Log \Lhard \leq 45.3$). About 60% of them are spectroscopically identified Type 1 AGN, while the others have a reliable photometric redshift and are classified as Type 1 AGN on the basis of their multi-band Spectral Energy Distributions. We discuss the relationship between UV and X-ray luminosity, as parameterized by the $\alphaox$ spectral slope, and its dependence on redshift and luminosity. We compare our findings with previous investigations of optically selected broad-line AGN (mostly from SDSS). A highly significant correlation between $\alphaox$ and $\lo$ is found, in agreement with previous investigations of optically selected samples. We calculate bolometric corrections, $\kbol$, for the whole sample using hard X-ray luminosities ($\Lhard$), and the Eddington ratios for a subsample of 150 objects for which black hole mass estimates are available. We confirm the trend of increasing bolometric correction with increasing Eddington ratio as proposed in previous works. A tight correlation is found between $\alphaox$ and $\kbol$, which can be used to estimate accurate bolometric corrections using only optical and X-ray data. We find a significant correlation between $\alphaox$ and Eddington ratio, in which $\alphaox$ increases for increasing Eddington ratios.

We follow the galaxy stellar mass assembly by morphological and spectral type in the COSMOS 2-deg^2 field. We derive the stellar mass functions and stellar mass densities from z=2 to z=0.2 using 196,000 galaxies selected at F(3.6 micron) > 1 microJy with accurate photometric redshifts (sigma_((zp-zs)/(1+zs))=0.008 at i<22.5). Using a spectral classification, we find that z~1 is an epoch of transition in the stellar mass assembly of quiescent galaxies. Their stellar mass density increases by 1.1 dex between z=1.5-2 and z=0.8-1 (Delta t ~2.5 Gyr), but only by 0.3 dex between z=0.8-1 and z~0.1 (Delta t ~ 6 Gyr). Then, we add the morphological information and find that 80-90% of the massive quiescent galaxies (log(M)~11) have an elliptical morphology at z<0.8. Therefore, a dominant mechanism links the shutdown of star formation and the acquisition of an elliptical morphology in massive galaxies. Still, a significant fraction of quiescent galaxies present a Spi/Irr morphology at low mass (40-60% at log(M)~9.5), but this fraction is smaller than predicted by semi-analytical models using a ``halo quenching'' recipe. We also analyze the evolution of star-forming galaxies and split them into ``intermediate activity'' and ``high activity'' galaxies. We find that the most massive ``high activity'' galaxies end their high star formation rate phase first. Finally, the space density of massive star-forming galaxies becomes lower than the space density of massive elliptical galaxies at z<1. As a consequence, the rate of ``wet mergers'' involved in the formation of the most massive ellipticals must decline very rapidly at z<1, which could explain the observed slow down in the assembly of these quiescent and massive sources.

We present a large robust sample of 1503 reliable and unconfused 70microm selected sources from the multiwavelength data set of the Cosmic Evolution Survey (COSMOS). Using the Spitzer IRAC and MIPS photometry, we estimate the total infrared luminosity, L_IR (8--1000 microns), by finding the best fit template from several different template libraries. The long wavelength 70 and 160 micron data allow us to obtain a reliable estimate of L_IR, accurate to within 0.2 and 0.05 dex, respectively. The 70 micron data point enables a significant improvement over the luminosity estimates possible with only a 24 micron detection. The full sample spans a wide range in L_IR, L_IR ~ 10^8-10^14 L_sun, with a median luminosity of 10^11.4 L_sun. We identify a total of 687 luminous, 303 ultraluminous, and 31 hyperluminous infrared galaxies (LIRGs, ULIRGs, and HyLIRGs) over the redshift range 0.01<z<3.5 with a median redshift of 0.5. Presented here are the full spectral energy distributions for each of the sources compiled from the extensive multiwavelength data set from the ultraviolet (UV) to the far-infrared (FIR). Using SED fits we find possible evidence for a subset of cooler ultraluminous objects than observed locally. However, until direct observations at longer wavelengths are obtained, the peak of emission and the dust temperature cannot be well constrained. We use these SEDs, along with the deep radio and X-ray coverage of the field, to identify a large sample of candidate active galactic nuclei (AGN). We find that the fraction of AGN increases strongly with L_IR, as it does in the local universe, and that nearly 70% of ULIRGs and all HyLIRGs likely host a powerful AGN.

We use ~8,600 >5e10 Msol COSMOS galaxies to study how the morphological mix of massive ellipticals, bulge-dominated disks, intermediate-bulge disks, bulge-less disks and irregular galaxies evolves from z=0.2 to z=1. The morphological evolution depends strongly on mass. At M>3e11 Msol, no evolution is detected in the morphological mix: ellipticals dominate since z=1, and the Hubble sequence has quantitatively settled down by this epoch. At the 1e11 Msol mass scale, little evolution is detected, which can be entirely explained with major mergers. Most of the morphological evolution from z=1 to z=0.2 takes place at masses 5e10 - 1e11 Msol, where: (i) The fraction of spirals substantially drops and the contribution of early-types increases. This increase is mostly produced by the growth of bulge-dominated disks, which vary their contribution from ~10% at z=1 to >30% at z=0.2 (cf. the elliptical fraction grows from ~15% to ~20%). Thus, at these masses, transformations from late- to early-types result in disk-less elliptical morphologies with a statistical frequency of only 30% - 40%. Otherwise, the processes which are responsible for the transformations either retain or produce a non-negligible disk component. (ii) The bulge-less disk galaxies, which contribute ~15% to the intermediate-mass galaxy population at z=1, virtually disappear by z=0.2. The merger rate since z=1 is too low to account for the disappearance of these massive bulge-less disks, which most likely grow a bulge via secular evolution.

We study the properties of a sample of 211 heavily-obscured Active Galactic Nucleus (AGN) candidates in the Extended Chandra Deep Field-South selecting objects with f_24/f_R>1000 and R-K>4.5. Of these, 18 were detected in X-rays and found to be obscured AGN with neutral hydrogen column densities of ~10^23 cm^-2. In the X-ray undetected sample, the following evidence suggests a large fraction of heavily-obscured (Compton Thick) AGN: (i) The stacked X-ray signal of the sample is strong, with an observed ratio of soft to hard X-ray counts consistent with a population of ~90% heavily obscured AGN combined with 10% star-forming galaxies. (ii) The X-ray to mid-IR ratios for these sources are significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller than for the general AGN population, suggesting column densities of N_H>5x10^24 cm^-2. (iii) The Spitzer near- and mid-IR colors of these sources are consistent with those of the X-ray-detected sample if the effects of dust self-absorption are considered. Spectral fitting to the rest-frame UV/optical light (dominated by the host galaxy) returns stellar masses of ~10^11 M_sun and <E(B-V)> =0.5, and reveals evidence for a significant young stellar population, indicating that these sources are experiencing considerable star-formation. This sample of heavily-obscured AGN candidates implies a space density at z~2 of ~10^-5 Mpc^-3, finding a strong evolution in the number of L_X>10^44 erg/s sources from z=1.5 to 2.5, possibly consistent with a short-lived heavily-obscured phase before an unobscured quasar is visible.

We have used the zCOSMOS-bright 10k sample to identify 3244 Spitzer/MIPS 24-micron-selected galaxies with 0.06< S(24um)< 0.50 mJy and I(AB)<22.5, over 1.5 deg^2 of the COSMOS field, and studied different spectral properties, depending on redshift. At 0.2<z<0.3, we found that different reddening laws of common use in the literature explain the dust extinction properties of around 80% of our infrared (IR) sources, within the error bars. For up to 16% of objects, instead, the Halpha/Hbeta ratios are too high for their IR/UV attenuations, which is probably a consequence of inhomogenous dust distributions. In only a few of our galaxies at 0.2<z<0.3 the IR emission could be mainly produced by dust heated by old rather than young stars. Besides, the line ratios of ~22% of our galaxies suggest that they might be star-formation/nuclear-activity composite systems. At 0.5<z<0.7, we estimated galaxy metallicities for 301 galaxies: at least 12% of them are securely below the upper-branch mass-metallicity trend, which is consistent with the local relation. Finally, we performed a combined analysis of the Hdelta equivalent-width versus Dn(4000) diagram for 1722 faint and bright 24um galaxies at 0.6<z<1.0, spanning two decades in mid-IR luminosity. We found that, while secondary bursts of star formation are necessary to explain the position of the most luminous IR galaxies in that diagram, quiescent, exponentially-declining star formation histories can well reproduce the spectral properties of ~40% of the less luminous sources. Our results suggest a transition in the possible modes of star formation at total IR luminosities L(TIR)=(3 +/-2)x10^11 Lsun.

[abridged] We present constraints on the nature of the first galaxies selected at 350 microns. The sample includes galaxies discovered in the deepest blank-field survey at 350 microns (in the Bootes Deep Field) and also later serendipitous detections in the Lockman Hole. Spectral energy distribution templates are fit to identified counterparts, and the sample is found to comprise IR-luminous galaxies at 1<z<3 predominantly powered by star formation. The first spectrum of a 350-micron-selected galaxy provides an additional confirmation, showing prominent dust grain features typically associated with star-forming galaxies. Compared to submillimeter galaxies selected at 850 and 1100 microns, galaxies selected at 350 microns have a similar range of far-infrared color temperatures. However, no 350-micron-selected sources are reliably detected at 850 or 1100 microns. Galaxies in our sample with redshifts 1<z<2 show a tight correlation between the far- and mid-infrared flux densities, but galaxies at higher redshifts show a large dispersion in their mid- to far-infrared colors. The 350-micron flux densities (15<S(350)<40 mJy) place these objects near the Herschel/SPIRE 350-micron confusion threshold, with the lower limit on the star formation rate density suggesting the bulk of the 350-micron contribution will come from less luminous infrared sources and normal galaxies.

We present the first results obtained from the identification of ~30,000 sources in the 24mic observations of the COSMOS field at S24>~80muJy. Using accurate photo-z's and extrapolations of the counts at faint fluxes we resolve with unprecedented detail the build-up of the mid-IR background across cosmic ages. We find that ~50% and ~80% of the 24mic background intensity originate from galaxies at z<~1 and z<~2 respectively, supporting the scenario where highly obscured sources at very high redshifts (z>~2) contribute only marginally to the CIB. Assuming flux-limited selections at optical wavelengths, we also find that the fraction of i-band sources with 24mic detection strongly increases up to z~2 as a consequence of the rapid evolution that star-forming galaxies have undergone with lookback time. Nonetheless this rising trend shows a clear break at z~1.3, probably due to k-correction effects implied by the complexity of SEDs in the mid-IR. Finally, we compare our results with the predictions from different models of galaxy formation. We note that semi-analytical formalisms fail to reproduce the redshift distributions observed at 24mic. Furthermore the simulated galaxies exhibit R-K colors much bluer than observed and the predicted K-band fluxes are systematically underestimated at z>0.5. Unless these discrepancies result from an incorrect treatment of extinction they may reflect an underestimate of the predicted density of high-z massive sources with strong star formation, which would point to more fundamental processes that are still not fully controlled in the simulations. The most recent backward evolution scenarios reproduce reasonably well the flux/redshift distribution of 24mic sources up to z~3, although none of them is able to exactly match our results at all redshifts. [Abridged]

We investigate the properties of 12 ultra-massive passively evolving early type galaxies (ETGs) at z_phot>1.4 in the COSMOS 2 deg^2 field. These 12 ETGs were selected as pBzKs, have accurate 1.4<= z_phot <=1.7,high Sersic index profiles typical of ellipticals, no detection at 24 micron, resulting in a complete ETG sample at M*>2.5x10^11 M_sun (Chabrier IMF). Contrary to previous claims, the half light radii estimated in very high S/N imaging data from HST+ACS are found to be large for most of the sample, consistent with local ellipticals. If the high redshift ETGs with M*<2.5x10^11 M_sun are really small in size and compact as reported in previous studies, our result may suggest a "downsizing" scenario, whereby the most massive ETGs reach their final structure earlier and faster than lower mass ones. However, simulating galaxies with morphological properties fixed to those of local ETGs with the same stellar mass show that the few compact galaxies that we still recover in our sample can be understood in term of fluctuations due to noise preventing the recovery of the extended low surface brightness halos in the light profile. Such halos, typical of Sersic profiles, extending even up to 40 kpc, are indeed seen in our sample.

We present SHARC-II 350um imaging of twelve 24um-bright (F_24um > 0.8 mJy) Dust-Obscured Galaxies (DOGs) and CARMA 1mm imaging of a subset of 2 DOGs, all selected from the Bootes field of the NOAO Deep Wide-Field Survey. Detections of 4 DOGs at 350um imply IR luminosities which are consistent within a factor of 2 of expectations based on a warm dust spectral energy distribution (SED) scaled to the observed 24um flux density. The 350um upper limits for the 8 non-detected DOGs are consistent with both Mrk231 and M82 (warm dust SEDs), but exclude cold dust (Arp220) SEDs. The two DOGs targeted at 1mm were not detected in our CARMA observations, placing strong constraints on the dust temperature: T_dust > 35-60 K. Assuming these dust properties apply to the entire sample, we find dust masses of ~3x10^8 M_sun. In comparison to other dusty z ~ 2 galaxy populations such as sub-millimeter galaxies (SMGs) and other Spitzer-selected high-redshift sources, this sample of DOGs has higher IR luminosities (2x10^13 L_sun vs. 6x10^12 L_sun for the other galaxy populations), warmer dust temperatures (>35-60 K vs. ~30 K), and lower inferred dust masses (3x10^8 M_sun vs. 3x10^9 M_sun). Herschel and SCUBA-2 surveys should be able to detect hundreds of these power-law dominated DOGs. We use HST and Spitzer/IRAC data to estimate stellar masses of these sources and find that the stellar to gas mass ratio may be higher in our 24um-bright sample of DOGs than in SMGs and other Spitzer-selected sources. Although larger sample sizes are needed to provide a definitive conclusion, the data are consistent with an evolutionary trend in which the formation of massive galaxies at z~2 involves a sub-millimeter bright, cold-dust and star-formation dominated phase followed by a 24um-bright, warm-dust and AGN-dominated phase.

We present Spitzer 70 and 160 micron observations of the COSMOS Spitzer survey (S-COSMOS). The data processing techniques are discussed for the publicly released products consisting of images and source catalogs. We present accurate 70 and 160 micron source counts of the COSMOS field and find reasonable agreement with measurements in other fields and with model predictions. The previously reported counts for GOODS-North and the extragalactic First Look Survey are updated with the latest calibration, and counts are measured based on the large area SWIRE survey to constrain the bright source counts. We measure an extragalactic confusion noise level of sigma_c = 9.4+/-3.3 mJy (q=5) for the MIPS 160-micron band based on the deep S-COSMOS data and report an updated confusion noise level of sigma_c = 0.35+/-0.15 mJy (q=5) for the MIPS 70-micron band.

We present the results from an optical and near-infrared spectroscopic study of the ultraviolet-luminous z = 2.73 galaxy, the 8 o'clock arc. Due to gravitational lensing, this galaxy is magnified by a factor of > 10, allowing in-depth measurements which are usually unfeasible at such redshifts. In the optical spectra, we measured the systemic redshift of the galaxy, z = 2.7322 +/- 0.0012, using stellar photospheric lines. This differs from the redshift of absorption lines in the interstellar medium, z = 2.7302 +/- 0.0006, implying gas outflows on the order of 160 km/s. With H and K-band near-infrared spectra, we have measured nebular emission lines of Halpha, Hbeta, Hgamma, [N II] and [O III], which have a redshift z = 2.7333 +/- 0.0001, consistent with the derived systemic redshift. From the Balmer decrement, we measured the dust extinction to be A_5500 = 1.17 +/- 0.36 mag. Correcting Halpha for dust extinction and the assumed lensing factor, we measure a star-formation rate of ~ 270 Msol/yr, which is higher than ~ 85% of star-forming galaxies at z ~ 2-3. Using combinations of all detected emission lines, we find that the 8 o'clock arc has a gas-phase metallicity of ~ 0.8 Zsol, showing that enrichment at high-redshift is not rare, even in blue, star-forming galaxies. Studying spectra from two of the arc components separately, we find that one component dominates the dust extinction and star-formation rate, although the metallicities between the two components are similar. We derive the mass via stellar population modeling, and find that the arc has a total stellar mass of ~ 4 x 10^11 Msol, which falls on the mass-metallicity relation at z ~ 2. Finally, we estimate the total gas mass, and find it to be only ~ 12% of the stellar mass, implying that the 8 o'clock arc is likely nearing the end of a starburst.

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we present low-resolution (64 < lambda / dlambda < 124), mid-infrared (20-38 micron) spectra of 23 high-redshift ULIRGs detected in the Bootes field of the NOAO Deep Wide-Field Survey. All of the sources were selected to have 1) fnu(24 micron) > 0.5 mJy; 2) R-[24] > 14 Vega mag; and 3) a prominent rest-frame 1.6 micron stellar photospheric feature redshifted into Spitzer's 3-8 micron IRAC bands. Of these, 20 show emission from polycyclic aromatic hydrocarbons (PAHs), usually interpreted as signatures of star formation. The PAH features indicate redshifts in the range 1.5 < z < 3.0, with a mean of <z>=1.96 and a dispersion of 0.30. Based on local templates, these sources have extremely large infrared luminosities, comparable to that of submillimeter galaxies. Our results confirm previous indications that the rest-frame 1.6 micron stellar bump can be efficiently used to select highly obscured starforming galaxies at z~2, and that the fraction of starburst-dominated ULIRGs increases to faint 24 micron flux densities. Using local templates, we find that the observed narrow redshift distribution is due to the fact that the 24 micron detectability of PAH-rich sources peaks sharply at z = 1.9. We can analogously explain the broader redshift distribution of Spitzer-detected AGN-dominated ULIRGs based on the shapes of their SEDs. Finally, we conclude that z~2 sources with a detectable 1.6 micron stellar opacity feature lack sufficient AGN emission to veil the 7.7 micron PAH band.

Optically bright Galactic stars (V < 13 mag) having fv(24 um) > 1 mJy are identified in Spitzer mid-infrared surveys within 8.2 square degrees for the Bootes field of the NOAO Deep Wide-Field Survey and within 5.5 square degrees for the First Look Survey (FLS). 128 stars are identified in Bootes and 140 in the FLS, and their photometry is given. (K-[24]) colors are determined using K magnitudes from the 2MASS survey for all stars in order to search for excess 24 um luminosity compared to that arising from the stellar photosphere. Of the combined sample of 268 stars, 141 are of spectral types F, G, or K, and 17 of these 141 stars have 24 um excesses with (K-[24]) > 0.2 mag. Using limits on absolute magnitude derived from proper motions, at least 8 of the FGK stars with excesses are main sequence stars, and estimates derived from the distribution of apparent magnitudes indicate that all 17 are main sequence stars. These estimates lead to the conclusion that between 9% and 17% of the main sequence FGK field stars in these samples have 24 um infrared excesses. This result is statistically similar to the fraction of stars with debris disks found among previous Spitzer targeted observations of much brighter, main sequence field stars.

We analyze a sample of galaxies chosen to have F24>0.5 mJy and satisfy a certain IRAC color criterion. IRS spectra yield redshifts, spectral types, and PAH luminosities, to which we add multi-wavelength broadband photometry. Stellar population modeling and IRS spectra together demonstrate that the double criteria used to select this sample have efficiently isolated massive star-forming galaxies at z~1.9. This is the first starburst-dominated ULIRG sample at high redshift with total infrared luminosity measured directly from FIR and millimeter photometry, and as such gives us the first accurate view of broadband SEDs for starburst galaxies at extremely high luminosity and at all wavelengths. Similar broadband data are assembled for three other galaxy samples -- local starburst galaxies, local AGN/ULIRGS, and a second 24mu-luminous z~1.9 sample dominated by AGN. L(PAH)/L(IR) for the new z~1.9 starburst sample is the highest ever seen, some three times higher than in local starbursts, whereas in AGNs this ratio is depressed below the starburst trend, often severely. Several pieces of evidence imply that AGNs exist in this starburst dominated sample, except two of which even host very strong AGN, while they still have very strong PAH emission. The ACS images show most objects have very extended morphologies in the rest-frame UV band, thus extended distribution of PAH molecules. Such an extended distribution prevents further destruction PAH molecules by central AGNs. We conclude that objects in this sample are ULIRGs powered mainly by starburst; and the total infrared luminosity density contributed by this type of objects is 0.9-2.6x 10^7 Lsol/Mpc^3.

We combine new Spitzer Space Telescope observations in the mid- and far-infrared with SCUBA 850 micron observations to improve the measurement of dust temperatures, masses and luminosities for 11 galaxies of the SCUBA Local Universe Galaxy Survey (SLUGS). By fitting dust models we measure typical dust masses of 10E7.9 M_sol and dust luminosities of ~ 10E10 L_sol, for galaxies with modest star formation rates. The data presented in this paper combined with previous observations show that cold dust is present in all types of spiral galaxies and is a major contributor to their total luminosity. Because of the lower dust temperature of the SCUBA sources measured in this paper, they have flatter Far-IR nu F_nu(160um)/nu F_nu(850um) slopes than the larger Spitzer Nearby Galaxies Survey (SINGS), the sample that provides the best measurements of the dust properties of galaxies in the nearby universe. The new data presented here added to SINGS extend the parameter space that is well covered by local galaxies, providing a comprehensive set of templates that can be used to interpret the observations of nearby and distant galaxies.

We present high spatial resolution optical and near-infrared imaging obtained using the ACS, WFPC2 and NICMOS cameras aboard the Hubble Space Telescope of 31 24um--bright z~2 Dust Obscured Galaxies (DOGs) identified in the Bootes Field of the NOAO Deep Wide-Field Survey. Although this subset of DOGs have mid-IR spectral energy distributions dominated by a power-law component suggestive of an AGN, all but one of the galaxies are spatially extended and not dominated by an unresolved component at rest-frame UV or optical wavelengths. The observed V-H and I-H colors of the extended components are 0.2-3 magnitudes redder than normal star-forming galaxies. All but 1 have axial ratios >0.3, making it unlikely that DOGs are composed of an edge-on star-forming disk. We model the spatially extended component of the surface brightness distributions of the DOGs with a Sersic profile and find effective radii of 1-6 kpc. This sample of DOGs is smaller than most sub-millimeter galaxies (SMGs), but larger than quiescent high-redshift galaxies. Non-parametric measures (Gini and M20) of DOG morphologies suggest that these galaxies are more dynamically relaxed than local ULIRGs. We estimate lower limits to the stellar masses of DOGs based on the rest-frame optical photometry and find that these range from ~10^(9-11) M_sun. If major mergers are the progenitors of DOGs, then these observations suggest that DOGs may represent a post-merger evolutionary stage.

Observations with the Spitzer Space Telescope have revealed a population of red-sequence galaxies with a significant excess in their 24-micron emission compared to what is expected from an old stellar population. We identify 900 red galaxies with 0.15<z<0.3 from the AGN and Galaxy Evolution Survey (AGES) selected from the NOAO Deep Wide-Field Survey Bootes field. Using Spitzer/MIPS, we classify 89 (~10%) with 24-micron infrared excess (f24>0.3 mJy). We determine the prevalence of AGN and star-formation activity in all the AGES galaxies using optical line diagnostics and mid-IR color-color criteria. Using the IRAC color-color diagram from the IRAC Shallow Survey, we find that 64% of the 24-micron excess red galaxies are likely to have strong PAH emission features in the 8-micron IRAC band. This fraction is significantly larger than the 5% of red galaxies with f24<0.3 mJy that are estimated to have strong PAH emission, suggesting that the infrared emission is largely due to star-formation processes. Only 15% of the 24-micron excess red galaxies have optical line diagnostics characteristic of star-formation (64% are classified as AGN and 21% are unclassifiable). The difference between the optical and infrared results suggest that both AGN and star-formation activity is occurring simultaneously in many of the 24-micron excess red galaxies. These results should serve as a warning to studies that exclusively use optical line diagnostics to determine the dominant emission mechanism in the infrared and other bands. We find that ~40% of the 24-micron excess red galaxies are edge-on spiral galaxies with high optical extinctions. The remaining sources are likely to be red galaxies whose 24-micron emission comes from a combination of obscured AGN and star-formation activity.

(abridged) We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24$\mu$m fluxes higher than 550$\mu$Jy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7--1.2 is $(3.7\pm1.1) \times10^{-5}$ Mpc$^{-3}$, corresponding to $\sim67%$ of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is $49\pm10%$, a value significantly higher than that previously estimated at similar 24$\mu$m fluxes. We discuss how our findings can constrain AGN feedback models.

We present the angular autocorrelation function of 2603 dust-obscured galaxies (DOGs) in the Bootes field of the NOAO Deep Wide-Field Survey. DOGs are red, obscured galaxies, defined as having R-[24] \ge 14 (F_24/F_R \ga 1000). Spectroscopy indicates that they are located at 1.5 \la z \la 2.5. We find strong clustering, with r_0 = 7.40^{+1.27}_{-0.84} Mpc/h for the full F_24 > 0.3 mJy sample. The clustering and space density of the DOGs are consistent with those of submillimeter galaxies, suggestive of a connection between these populations. We find evidence for luminosity-dependent clustering, with the correlation length increasing to r_0 = 12.97^{+4.26}_{-2.64} Mpc/h for brighter (F_24 > 0.6 mJy) DOGs. Bright DOGs also reside in richer environments than fainter ones, suggesting these subsamples may not be drawn from the same parent population. The clustering amplitudes imply average halo masses of log M = 12.2^{+0.3}_{-0.2} Msun for the full DOG sample, rising to log M = 13.0^{+0.4}_{-0.3} Msun for brighter DOGs. In a biased structure formation scenario, the full DOG sample will, on average, evolve into ~ 3 L* present-day galaxies, whereas the most luminous DOGs may evolve into brightest cluster galaxies.

We investigate the close environment of 203 Spitzer 24 micron-selected sources at 0.6<z<1.0 using zCOSMOS-bright redshifts and spectra of I<22.5 AB mag galaxies, over 1.5 sq. deg. of the COSMOS field. We quantify the degree of passivity of the LIRG and ULIRG environments by analysing the fraction of close neighbours with Dn(4000)>1.4. We find that LIRGs at 0.6<z<0.8 live in more passive environments than those of other optical galaxies that have the same stellar mass distribution. Instead, ULIRGs inhabit more active regions (e.g. LIRGs and ULIRGs at 0.6<z<0.8 have, respectively, (42.0 +/- 4.9)% and (24.5 +/- 5.9)% of neighbours with Dn (4000)>1.4 within 1 Mpc and +/- 500 km/s). The contrast between the activities of the close environments of LIRGs and ULIRGs appears especially enhanced in the COSMOS field density peak at z~0.67, because LIRGs on this peak have a larger fraction of passive neighbours, while ULIRGs have as active close environments as those outside the large-scale structure. The differential environmental activity is related to the differences in the distributions of stellar mass ratios between LIRGs/ULIRGs and their close neighbours, as well as in the general local density fields. At 0.8<z<1.0, instead, we find no differences in the environment densities of ULIRGs and other similarly massive galaxies, in spite of the differential activities. We discuss a possible scenario to explain these findings.

We present accurate photometric redshifts in the 2-deg2 COSMOS field. The redshifts are computed with 30 broad, intermediate, and narrow bands covering the UV (GALEX), Visible-NIR (Subaru, CFHT, UKIRT and NOAO) and mid-IR (Spitzer/IRAC). A chi2 template-fitting method (Le Phare) was used and calibrated with large spectroscopic samples from VLT-VIMOS and Keck-DEIMOS. We develop and implement a new method which accounts for the contributions from emission lines (OII, Hbeta, Halpha and Ly) to the spectral energy distributions (SEDs). The treatment of emission lines improves the photo-z accuracy by a factor of 2.5. Comparison of the derived photo-z with 4148 spectroscopic redshifts (i.e. Delta z = zs - zp) indicates a dispersion of sigma_{Delta z/(1+zs)}=0.007 at i<22.5, a factor of 2-6 times more accurate than earlier photo-z in the COSMOS, CFHTLS and COMBO-17 survey fields. At fainter magnitudes i<24 and z<1.25, the accuracy is sigma_{Delta z/(1+zs)}=0.012. The deep NIR and IRAC coverage enables the photo-z to be extended to z~2 albeit with a lower accuracy (sigma_{Delta z/(1+zs)}=0.06 at i~24). The redshift distribution of large magnitude-selected samples is derived and the median redshift is found to range from z=0.66 at 22<i<22.5 to z=1.06 at 24.5<i<25. At i<26.0, the multi-wavelength COSMOS catalog includes approximately 607,617 objects. The COSMOS-30 photo-z enable the full exploitation of this survey for studies of galaxy and large scale structure evolution at high redshift.

The collapsar model predicts that the progenitors of Gamma-ray Bursts (GRBs) are metal poor in Fe group elements. The existence of low metallicity stellar populations could manifest itself through special characteristics of the immediate environment of the GRB site in the host galaxy. We analyse the strong emission lines from the sub-luminous host galaxy of GRB 980425, which showed the first connection with a supernova explosion (SN 1998bw). The host has a sufficient size to allow detailed resolved spectroscopy of individual HII regions and to look for regions with special properties close to the the GRB site. Using integral field spectroscopy with VIMOS we cover most of the high surface brightness part of the host including the HII region where the supernova and GRB occurred. The star formation rate, reddening, equivalent width and stellar mass in the GRB region is similar to other HII regions in the host. Extreme values arise in the only region that shows emission lines from Wolf-Rayet stars; a region that lies 800 pc in projection from the GRB site. Strong emission line diagnostics of all HII regions imply oxygen abundances between 0.3 and 0.8 solar with the lowest values arising in the WR and GRB regions. Comparing the measured emission line ratios for a sample of low redshift GRB hosts to theoretical models and to observations of field galaxies, we find that GRBs arise in a range of metallicity environments while the regions are consistently very young. The similar line ratios of GRB hosts compared to the WR region can arise in spatially unresolved galaxies with bright HII regions near the GRB location. (Abridged)

We have combined data from the NOAO Deep Wide-Field Survey in Bootes and the Spitzer Space Telescope to determine basic properties for sixteen optically "invisible" MIPS 24um (OIMS) and thirty-five optically "invisible" radio (OIRS) soruces, including their spectral energy distributions (SED) and luminosities. Most OIMSs possess steep power-law SEDs over lambda(rest) = 1-10 um, indicating the presence of obscured AGN in agreement with Spitzer spectroscopy. These objects are extremely luminous at rest-frame near and mid-IR (nu L_nu(5um) ~ 10^{38}-10^{39} W), consistent with accretion near the Eddington limit and further implying that they host buried QSOs. The majority of the IRAC detected OIRSs have flat 3.6 to 24um SEDs, implying comparable emission from stellar photospheres and hot AGN illuminated dust. This may reflect relatively small amounts of dust close to the central engine or current low mass accretion rates. A small subset of OIRSs appear to be starburst dominated with photometric redshifts from 1.0 to 4.5. The OIMSs and OIRSs with significant starburst components have similar L_K and stellar masses (M* ~ 10^{11} M_solar) assuming minimal AGN contribution. Roughly half of the OIRSs are not detected by Spitzer's IRAC or MIPS. These are most likely z > 2 radio galaxies. The IRAC detected OIRSs are more likely than OIMSs to appear non point-like in the 3.6um and 4.5um images, suggesting that interactions play a role in triggering their activity. The AGN powered OIMSs may represent sub-millimeter galaxies making the transition from starburst to accretion dominance in their evolution to current epoch massive ellipticals.

We study zCOSMOS-bright optical spectra for 609 Spitzer/MIPS 24 micron-selected galaxies with S(24um)> 0.30 mJy and I<22.5 (AB mag) over 1.5 sq. deg. of the COSMOS field. From emission-line diagnostics we find that: 1) star-formation rates (SFR) derived from the observed Halpha and Hbeta lines underestimate, on average, the total SFR by factors ~5 and 10, respectively; 2) both the Calzetti et al. and the Milky Way reddening laws are suitable to describe the extinction observed in infrared (IR) sources in most cases; 3) some IR galaxies at z<0.3 have low abundances, but many others with similar IR luminosities and redshifts are chemically enriched; 4) The average [OIII]/Hbeta ratios of nuLnu(24um)>10^11 Lsun galaxies at 0.6<z<0.7 are ~0.6 dex higher than the average ratio of all zCOSMOS galaxies at similar redshifts. Massive star formation and active galactic nuclei (AGN) could simultaneously be present in those galaxies with the highest ionising fluxes; 5) ~1/3 of the galaxies with metallicity measurements at 0.5<z<0.7 lie below the general mass-metallicity relation at the corresponding redshifts. The strengths of the 4000 Angstrom break and the Hdelta EW of our galaxies show that secondary bursts of star formation are needed to explain the spectral properties of most IR sources. The LIRG and ULIRG phases occur, on average, between 10^7 and 10^8 years after the onset of a starburst on top of underlying older stellar populations. These results are valid for galaxies of different IR luminosities at 0.6<z<1.0 and seem independent of the mechanisms triggering star formation.

We investigate whether the mean star formation activity of star-forming galaxies from z=0 to z=0.7 in the GOODS-S field can be reproduced by simple evolution models of these systems. In this case, such models might be used as first order references for studies at higher z to decipher when and to what extent a secular evolution is sufficient to explain the star formation history in galaxies. We selected star-forming galaxies at z=0 and at z=0.7 in IR and in UV to have access to all the recent star formation. We focused on galaxies with a stellar mass ranging between 10^{10} and 10^{11} M_sun for which the results are not biased by the selections. We compared the data to chemical evolution models developed for spiral galaxies and originally built to reproduce the main characteristics of the Milky Way and nearby spirals without fine-tuning them for the present analysis. We find a shallow decrease in the specific star formation rate (SSFR) when the stellar mass increases. The evolution of the SSFR characterizing both UV and IR selected galaxies from z=0 to z=0.7 is consistent with the models built to reproduce the present spiral galaxies. There is no need to strongly modify of the physical conditions in galaxies to explain the average evolution of their star formation from z=0 to z=0.7. We use the models to predict the evolution of the star formation rate and the metallicity on a wider range of redshift and we compare these predictions with the results of semi-analytical models.

We present the redshift distribution of a complete, unbiased sample of 24 micron sources down to fnu(24 micron) = 300 uJy (5-sigma). The sample consists of 591 sources detected in the Bootes field of the NOAO Deep Wide-Field Survey. We have obtained optical spectroscopic redshifts for 421 sources (71%). These have a redshift distribution peaking at z~0.3, with a possible additional peak at z~0.9, and objects detected out to z=4.5. The spectra of the remaining 170 (29%) exhibit no strong emission lines from which to determine a redshift. We develop an algorithm to estimate the redshift distribution of these sources, based on the assumption that they have emission lines but that these lines are not observable due to the limited wavelength coverage of our spectroscopic observations. The redshift distribution derived from all 591 sources exhibits an additional peak of extremely luminous (L(8-1000 micron) > 3 x 10^{12} Lsun) objects at z~2, consisting primarily of sources without observable emission lines. We use optical line diagnostics and IRAC colors to estimate that 55% of the sources within this peak are AGN-dominated. We compare our results to published models of the evolution of infrared luminous galaxies. The models which best reproduce our observations predict a large population of star-formation dominated ULIRGs at z > 1.5 rather than the AGN-dominated sources we observe.

We present mid-infrared spectroscopy of a sample of 16 optically faint infrared luminous galaxies obtained with the Infrared Spectrograph (IRS) on the Spitzer Space Telescope. These sources were jointly selected from Spitzer and Chandra imaging surveys in the NDWFS Bootes field and were selected from their bright X-ray fluxes to host luminous AGN. None of the spectra show significant emission from polycyclic aromatic hydrocarbons (PAHs; 6.2um equivalent widths <0.2um), consistent with their infrared emission being dominated by AGN. Nine of the X-ray sources show 9.7um silicate absorption features. Their redshifts are in the range 0.9<z<2.6, implying infrared luminosities of log(L{IR})=12.5-13.6 solar luminosities. The average silicate absorption strength is not as strong as that of previously targeted optically faint infrared luminous galaxies with similar mid-infrared luminosities implying that the X-ray selection favors sources behind a smaller column of Si-rich dust than non-X-ray selection. Seven of the X-ray sources have featureless power-law mid-IR spectra. We argue that the featureless spectra likely result from the sources having weak or absent silicate and PAH features rather than the sources lying at higher redshifts where these features are shifted out of the IRS spectral window. We investigate whether there are any correlations between X-ray and infrared properties and find that sources with silicate absorption features tend to have fainter X-ray fluxes and harder X-ray spectra, indicating a weak relation between the amount of silicate absorption and column density of X-ray-absorbing gas.

We present the observed correlations between rest-frame 8, 24, 70 and 160 um monochromatic luminosities and measured total infrared luminosities L_IR of galaxies detected by Spitzer. Our sample consists of 372 star-forming galaxies with individual detections and flux measurements at 8, 24, 70 and 160 um. We have spectroscopic redshifts for 93% of these sources, and accurate photometric redshifts for the remainder. We also used a stacking analysis to measure the IR fluxes of fainter sources at higher redshifts. We show that the monochromatic mid and far-infrared luminosities are strongly correlated with the total infrared luminosity and our stacking analysis confirms that these correlations also hold at higher redshifts. We provide relations between monochromatic luminosities and total infrared luminosities L_IR that should be reliable up to z~2 (z~1.1) for ULIRGs (LIRGs). In particular, we can predict L_IR with accuracies of 37% and 54% from the 8 and 24 um fluxes, while the best tracer is the 70 um flux. Combining bands leads to slightly more accurate estimates. For example, combining the 8 and 24 um luminosities predicts L_IR with an accuracy of 34%. Our results are generally compatible with previous studies, and the small changes are probably due to differences in the sample selection criteria. We can rule out strong evolution in dust properties with redshift up to z~1. Finally, we show that infrared and sub-millimeter observations are complementary means of building complete samples of star-forming galaxies, with the former being more sensitive for z<~2 and the latter at higher z>~2.

We have used the Spitzer Space Telescope to study the dust properties of a sample of star-forming dwarf galaxies. The differences in the mid-infrared spectral energy distributions for these galaxies which, in general, are low metallicity systems, indicate differences in the physical properties, heating, and/or distribution of the dust. Specifically, these galaxies have more hot dust and/or very small grains and less PAH emission than either spiral or higher luminosity starburst galaxies. As has been shown in previous studies, there is a gradual decrease in PAH emission as a function of metallicity. Because much of the energy from star formation in galaxies is re-radiated in the mid-infrared, star-formation rate indicators based on both line and continuum measurements in this wavelength range are coming into more common usage. We show that the variations in the interstellar medium properties of galaxies in our sample, as measured in the mid-infrared, result in over an order of magnitude spread in the computed star-formation rates.

We present mid-infrared spectroscopy obtained with the Spitzer Space Telescope of a sample of 11 optically faint, infrared luminous galaxies selected from a Spitzer MIPS 70um imaging survey of the NDWFS Bootes field. These are the first Spitzer IRS spectra presented of distant 70um-selected sources. All the galaxies lie at redshifts 0.3<z<1.3 and have very large infrared luminosities of L_IR~ 0.1-17 x 10^12 solar luminosities. Seven of the galaxies exhibit strong emission features attributed to polycyclic aromatic hydrocarbons (PAHs). The average IRS spectrum of these sources is characteristic of classical starburst galaxies, but with much larger infrared luminosities. The PAH luminosities of L(7.7) ~ 0.4 - 7 x 10^11 solar luminosities imply star formation rates of ~ 40 - 720 solar masses per year. Four of the galaxies show deep 9.7um silicate absorption features and no significant PAH emission features (6.2um equivalent widths < 0.03um). The large infrared luminosities and low f70/f24 flux density ratios suggests that these sources have AGN as the dominant origin of their large mid-infrared luminosities, although deeply embedded but luminous starbursts cannot be ruled out. If the absorbed sources are AGN-dominated, a significant fraction of all far-infrared bright, optically faint sources may be dominated by AGN.

We present the quantitative rest-frame B morphological evolution and galaxy merger fractions at 0.2 < z < 1.2 as observed by the All-wavelength Extended Groth Strip International Survey (AEGIS). We use the Gini coefficent and M_20 to identify major mergers and classify galaxy morphology for a volume-limited sample of 3009 galaxies brighter than 0.4 L_B^*, assuming pure luminosity evolution of 1.3 M_B per unit redshift. We find that the merger fraction remains roughly constant at 10 +/- 2% for 0.2 < z < 1.2. The fraction of E/S0/Sa increases from 21+/- 3% at z ~ 1.1 to 44 +/- 9% at z ~ 0.3, while the fraction of Sb-Ir decreases from 64 +/- 6% at z ~ 1.1 to 47 +/- 9% at z ~ 0.3. The majority of z < 1.2 Spitzer MIPS 24 micron sources with L(IR) > 10^11 L_sun are disk galaxies, and only ~ 15% are classified as major merger candidates. Edge-on and dusty disk galaxies (Sb-Ir) are almost a third of the red sequence at z ~ 1.1, while E/S0/Sa makeup over 90% of the red sequence at z ~ 0.3. Approximately 2% of our full sample are red mergers. We conclude (1) the galaxy merger rate does not evolve strongly between 0.2 < z < 1.2; (2) the decrease in the volume-averaged star-formation rate density since z ~ 1 is a result of declining star-formation in disk galaxies rather than a disappearing population of major mergers; (3) the build-up of the red sequence at z < 1 can be explained by a doubling in the number of spheroidal galaxies since z ~ 1.2.

Ultraviolet (UV) galaxies have been selected from GALEX. The presence of a FUV-dropout in their spectral energy distributions proved to be a very complete (83.3%) but not very efficient (21.4%) tool for identifying Lyman Break Galaxies (LBGs) at z~1. We divide the LBG sample into two sub-classes: red LBGs (RLBGs) detected at 24 micron which are mainly Luminous IR Galaxies (LIRGs) and blue LBGs (BLBGs) undetected at 24 microns down to 83 microJy. Two of the RLBGs are also detected at 70 micron. The median SED of the RLBGs is similar (above lambda~1 micron) to the dusty starburst HR10. However, unlike local (U)LIRGs, RLBGs are UV bright objects. We suggest that these objects contain a large amount of dust but that some bare stellar populations are also directly visible. The median SED of the BLBGs is consistent with their containing the same stellar population as the RLBGs but with a lower dust content. The luminosity function of our LBG sample at z~1 is similar to the luminosity function of NUV-selected galaxies at the same redshift. The integrated luminosity densities of z~1 LBGs and NUV-selected galaxies are very consistent. We show that star formation rates (SFRs) estimated from UV measurements and corrected using the IRX-beta method provide average total SFR_TOT in agreement with SFR_UV + SFR_dust. However, IRX-beta-based SFR_TOT shows a large dispersion. Summing up the detected UV (1150A rest-frame) and IR-based star formation rates of the detected objects, we find that only one third of the total (i.e. UV + dust) LBG SFR resides in BLBGs and two thirds in RLBGs, even though most LBGs at z~1 are BLBGs. On the other hand, the total SFR of LBGs accounts for only 11% of the total SFR at z~1. Finally, we observe a regular decrease of L_TIR / L_FUV from z=0 to z~2 for UV-selected samples.

A complete flux-limited sample of 50 galaxies is presented having f_{\nu}(24um) > 10mJy, chosen from a survey with the Multiband Imaging Photometer on Spitzer (MIPS) of 8.2 deg^{2} within the NOAO Deep Wide-Field Survey region in Bootes (NDWFS). Spectra obtained with the low-resolution modules of the Infrared Spectrograph on Spitzer (IRS) are described for 36 galaxies within this sample; 25 show strong PAH emission features characteristic of starbursts, and 11 show silicate absorption or emission, emission lines, or featureless spectra characteristic of AGN. Infrared or optical spectral classifications are available for 48 of the entire sample of 50; 33 galaxies are classified as starbursts and 15 as AGN. (There are an additional 19 Galactic stars with f_{\nu}(24um) > 10mJy in the survey area.) Using a relation between 7.7um PAH luminosity and star formation rate derived from previous IRS observations of starbursts, the star formation rate per unit volume of the local universe (SFRD) is determined from the complete sample and is found to be 0.008 \mdot Mpc^{-3}. Individual sources in the sample have star formation rates from 0.14 to 160 \mdot. The derived value for the local SFRD is about half that of the local SFRD deduced from bolometric luminosities of the IRAS 60um Bright Galaxy Sample, with the deficiency being at lower luminosities and arising primarily from the small number of low luminosity sources in the 10 mJy sample. The agreement for higher luminosities confirms the validity of using the 7.7um PAH feature as a measure of SFRD in the high redshift universe, where this is often the only indicator available for faint sources.

We study the infrared (IR) properties of high-redshift galaxies using deep Spitzer 24, 70, and 160 micron data. Our primary interest is to improve the constraints on the total IR luminosities, L(IR), of these galaxies. We combine the Spitzer data in the southern Extended Chandra Deep Field with a K-band-selected galaxy sample and photometric redshifts from the Multiwavelength Survey by Yale-Chile. We used a stacking analysis to measure the average 70 and 160 micron flux densities of 1.5 < z < 2.5 galaxies as a function of 24 micron flux density, X-ray activity, and rest-frame near-IR color. Galaxies with 1.5 < z < 2.5 and S(24)=53-250 micro-Jy have L(IR) derived from their average 24-160 micron flux densities within factors of 2-3 of those derived from the 24 micron flux densities only. However, L(IR) derived from the average 24-160 micron flux densities for galaxies with S(24) > 250 micro-Jy and 1.5 < z < 2.5 are lower than those derived using only the 24 micron flux density by factors of 2-10. Galaxies with S(24) > 250 micro-Jy have S(70)/S(24) flux ratios comparable to sources with X-ray detections or red rest-frame IR colors, suggesting that warm dust possibly heated by AGN may contribute to the high 24 micron emission. Based on the average 24-160 micron flux densities, nearly all 24 micron-selected galaxies at 1.5 < z < 2.5 have L(IR) < 6 x 10^12 solar luminosities, which if attributed to star formation corresponds to < 1000 solar masses per year. This suggests that high redshift galaxies may have similar star formation efficiencies and feedback processes as local analogs. Objects with L(IR) > 6 x 10^12 solar luminosities are quite rare, with a surface density ~ 30 +/- 10 per sq. deg, corresponding to ~ 2 +/- 1 x 10^-6 Mpc^-3 over 1.5 < z < 2.5.

We analyze the infrared (IR) spectral energy distributions (SEDs) for 10micron < lambda(rest) < 100micron for ~600 galaxies at z~0.7 in the extended Chandra Deep Field South by stacking their Spitzer 24, 70 and 160micron images. We place interesting constraints on the average IR SED shape in two bins: the brightest 25% of z~0.7 galaxies detected at 24micron, and the remaining 75% of individually-detected galaxies. Galaxies without individual detections at 24micron were not well-detected at 70micron and 160micron even through stacking. We find that the average IR SEDs of z~0.7 star-forming galaxies fall within the diversity of z~0 templates. While dust obscuration Lir/Luv seems to be only a function of star formation rate (SFR; ~ Lir+Luv), not of redshift, the dust temperature of star-forming galaxies (with SFR ~ 10 solar mass per year) at a given IR luminosity was lower at z~0.7 than today. We suggest an interpretation of this phenomenology in terms of dust geometry: intensely star-forming galaxies at z~0 are typically interacting, and host dense centrally-concentrated bursts of star formation and warm dust temperatures. At z~0.7, the bulk of intensely star-forming galaxies are relatively undisturbed spirals and irregulars, and we postulate that they have large amounts of widespread lower-density star formation, yielding lower dust temperatures for a given IR luminosity. We recommend what IR SEDs are most suitable for modeling intermediate redshift galaxies with different SFRs.

The Hubble Ultra Deep Field source JD2 presented in Mobasher et al. (2005) is an unusual galaxy that is very faint at all wavelengths shortward of 1.1 micron. Photometric redshift fits to data at 0.4 to 8 microns yield a significant probability that it is an extremely massive galaxy at z~6.5. In this paper we present new photometry at 16 microns and 22 micron from Spitzer Infrared Spectrograph (IRS) peak-up imaging of the Great Observatories Origins Deep Survey (GOODS) fields. We find that the spectral energy distribution shows a factor of ~4 rise in flux density between the 16 micron and 22 micron bandpass which is most likely due to the entrance of polycyclic aromatic hydrocarbon emission features into the 22 micron and 24 micron passbands. The flux ratio between these bandpasses can be best fit by a z~1.7 luminous infrared galaxy with a bolometric luminosity of (2-6)\times10^{11} L_sun corresponding to a star-formation rate of 80 M_sun/yr. The predicted flux density values at other longer wavelengths are below the detection limits of current instrumentation but such sources could potentially be detected in lensed submillimeter surveys. Re-evaluation of the optical/near-infrared photometry continues to favor z>6 photometric redshift solutions, but we argue that the consistency of the multiwavelength parameters of this galaxy with other dusty starbursts favor the z~2 mid-infrared photometric redshift. The data presented here provide evidence that optically undetected near-infrared sources which are detected at 24 microns are most likely dusty, starburst galaxies at a redshift of z~2 with stellar masses >10^{10} M_sun.

A SST survey in the NOAO Deep-Wide Field in Bo\"otes provides a complete, 8-micron-selected sample of galaxies to a limiting (Vega) magnitude of 13.5. In the 6.88 deg$^2$ field sampled, 79% of the 4867 galaxies have spectroscopic redshifts, allowing an accurate determination of the local (z<0.3) galaxy luminosity function. Stellar and dust emission can be separated on the basis of observed galaxy colors. Dust emission (mostly PAH) accounts for 80% of the 8 micron luminosity, stellar photospheres account for 19%, and AGN emission accounts for roughly 1 %. A sub-sample of the 8 micron-selected galaxies have blue, early-type colors, but even most of these have significant PAH emission. The luminosity functions for the total 8 micron luminosity and for the dust emission alone are both well fit by Schechter functions. For the 8 micron luminosity function, the characteristic luminosity is \nu L_{\nu}^*(8.0 \micron) = 1.8 \times 10^{10}$ \Lsun while for the dust emission alone it is 1.6 x 10^{10}$ \Lsun\null. The average 8 \micron luminosity density at z<0.3 is 3.1 x 10^7 \Lsun Mpc^{-3}, and the average luminosity density from dust alone is 2.5 x 10^7 \Lsun Mpc^{-3}. This luminos ity arises predominantly from galaxies with 8 \micron luminosities ($\nu L_{\nu}$) between $2\times 10^9$ and $2 x 10^{10}$ \Lsun, i.e., normal galaxies, not LIRGs or ULIRGs.

We have imaged a $\sim$6 arcminute$^2$ region in the Bo\"otes Deep Field using the 350 $\mu$m-optimised second generation Submillimeter High Angular Resolution Camera (SHARC II), achieving a peak 1$\sigma$ sensitivity of $\sim$5 mJy. We detect three sources above 3$\sigma$, and determine a spurious source detection rate of 1.09 in our maps. In the absence of $5\sigma$ detections, we rely on deep 24 $\mu$m and 20 cm imaging to deduce which sources are most likely to be genuine, giving two real sources. From this we derive an integral source count of 0.84$^{+1.39}_{-0.61}$ sources arcmin$^{-2}$ at $S>13$ mJy, which is consistent with 350 $\mu$m source count models that have an IR-luminous galaxy population evolving with redshift. We use these constraints to consider the future for ground-based short-submillimetre surveys.

Using data from the DEEP2 Galaxy Redshift Survey and HST/ACS imaging in the Extended Groth Strip, we select nearly 100 interacting galaxy systems including kinematic close pairs and morphologically identified merging galaxies. Spitzer MIPS 24 micron fluxes of these systems reflect the current dusty star formation activity, and at a fixed stellar mass (M_{*}) the median infrared luminosity (L_{IR}) among merging galaxies and close pairs of blue galaxies is twice (1.9 +/- 0.4) that of control pairs drawn from isolated blue galaxies. Enhancement declines with galaxy separation, being strongest in close pairs and mergers and weaker in wide pairs compared to the control sample. At z ~ 0.9, 7.1% +/- 4.3% of massive interacting galaxies (M_{*} > 2*10^{10} M_{solar}) are found to be ULIRGs, compared to 2.6% +/- 0.7% in the control sample. The large spread of IR luminosity to stellar mass ratio among interacting galaxies suggests that this enhancement may depend on the merger stage as well as other as yet unidentified factors (e.g., galaxy structure, mass ratio, orbital characteristics, presence of AGN or bar). The contribution of interacting systems to the total IR luminosity density is moderate (<= 36 %).

We analyze star formation (SF) as a function of stellar mass (M*) and redshift z in the All Wavelength Extended Groth Strip International Survey (AEGIS). For 2905 field galaxies, complete to 10^10(10^10.8) Msun at z<0.7(1), with Keck spectroscopic redshifts out to z=1.1, we compile SF rates (SFR) from emission lines, GALEX, and Spitzer MIPS 24 micron photometry, optical-NIR M* measurements, and HST morphologies. Galaxies with reliable signs of SF form a distinct "main sequence (MS)", with a limited range of SFR at a given M* and z (1 sigma < +-0.3 dex), and log(SFR) approximately proportional to log(M*). The range of log(SFR) remains constant to z>1, while the MS as a whole moves to higher SFR as z increases. The range of SFR along the MS constrains the amplitude of episodic variations of SF, and the effect of mergers on SFR. Typical galaxies spend ~67(95)% of their lifetime since z=1 within a factor of <~ 2(4) of their average SFR at a given M* and z. The dominant mode of the evolution of SF since z~1 is apparently a gradual decline of the average SFR in most individual galaxies, not a decreasing frequency of starburst episodes, or a decreasing factor by which SFR are enhanced in starbursts. LIRGs at z~1 seem to mostly reflect the high SFR typical for massive galaxies at that epoch. The smooth MS may reflect that the same set of few physical processes governs star formation prior to additional quenching processes. A gradual process like gas exhaustion may play a dominant role.

We analyze star formation (SF) as a function of stellar mass (M_*) and redshift z in the All Wavelength Extended Groth Strip International Survey (AEGIS), for star-forming field galaxies with M_* >~ 10^10 M_sun out to z=1.1. The data indicate that the high specific SF rates (SFR) of many less massive galaxies do not represent late, irregular or recurrent, starbursts in evolved galaxies. They rather seem to reflect the onset (initial burst) of the dominant SF episode of galaxies, after which SF gradually declines on Gyr timescales to z=0 and forms the bulk of a galaxy's M_*. With decreasing mass, this onset of major SF shifts to decreasing z for an increasing fraction of galaxies (staged galaxy formation). This process may be an important component of the ``downsizing'' phenomenon. We find that the predominantly gradual decline of SFR (Noeske et al. 2007, this volume) can be reproduced by exponential SF histories (tau models), if less massive galaxies have systematically longer e-folding times tau, and a later onset of SF (z_f). Our model can provide a first parametrization of SFR as a function of M_* and z, and quantify mass-dependences of tau and z_f, from direct observations of M_* and SFR up to z>1. The observed evolution of SF in galaxies can plausibly reflect the dominance of gradual gas exhaustion. The data are also consistent with the history of cosmological accretion onto Dark Matter halos.

The TIR and FUV luminosity functions of galaxies and the related luminosity densities rho(TIR) and rho(FUV) are known to evolve from z=0 to z~1 but with a different rate: the galaxy populations appear brighter in the past at both wavelengths but the evolution in TIR is larger than in FUV. It leads to an increase of the ratio of TIR to FUV luminosity densities rho(TIR)/rho(FUV) which can be interpreted as a global increase of the dust attenuation from z=0 to z~1. Our aim is to understand the origin of this increase. We focus on infrared galaxies with L(TIR>10^{11} Lsun at z~0.7 observed by SPITZER/MIPS and we measure their UV emission at 2310 A from GALEX. These Luminous InfraRed Galaxies (LIRGs) represent the bulk of the TIR luminosity density at intermediate redshift. Some evolution of L(TIR)/L(FUV) and therefore of dust attenuation is found for these galaxies: LIRGs at z=0.7 span a larger range of L(TIR)/L(FUV) ratios than at z=0 and their mean dust attenuation at FUV is found ~0.5 mag lower than for their local counterparts. The decrease of dust attenuation is found lower than that reported in other studies for bright galaxies selected in UV rest-frame at z=1 and 2. From a semi-quantitative analysis we find that the slight decrease of dust attenuation for LIRGs at z=0.7 remains consistent with the increase of rho(TIR)/rho(FUV) with redshift.

We present near-infrared spectroscopic observations for a sample of ten optically faint luminous infrared galaxies (R-[24]> 14) using Keck NIRSPEC and Gemini NIRI. The sample is selected from a 24 micron Spitzer MIPS imaging survey of the NDWFS Bootes field. We measure accurate redshifts in the range 1.3<z<3.4. Based on either emission line widths or line diagnostics, we find that all ten galaxies harbor luminous AGN. Seven sources are type I AGN, exhibiting broad (>1900 km/s) Halpha or Hbeta emission lines; the remaining three are type II AGN. Given their large mid-IR luminosities and faint optical magnitudes, we might expect these sources to be heavily extincted quasars, and therefore only visible as type II AGN. The visibility of broad lines in 70% of the sources suggests that it is unlikely that these AGN are being viewed through the mid-plane of a dusty torus. For four of the sources we constrain the Halpha/Hbeta Balmer decrement and estimate the extinction to the emission line region to be large for both type I and type II AGN, with A_Halpha > 2.4-5 mag. Since the narrow-line region is also extincted and the UV continuum emission from the host galaxies is extremely faint, this suggests that much of the obscuration is contributed by dust on large (~kpc) scales within the host galaxies. These sources may be examples of "host-obscured" AGN which could have space densities comparable or greater to that of optically luminous type I AGN with similar bolometric luminosities.

We use deep GALEX images of CDFS in UV to define the first large sample of 420 Lyman Break Galaxies at z~1. We use a PSF fitting to estimate UV magnitudes on these deep crowded images. Deep Spitzer IRAC and MIPS provide the first detection of a large sample of Lyman Break Galaxies in the mid- to far-infrared range. We are therefore able to study and compare the UV and TIR emission of Lyman Break Galaxies. We find that about 15% of the LBG sample are strong emitters at 24 microns (Red LBGs). Most of them are Luminous IR Galaxies (LIRGs) while the rest (Blue LBGs) are undetected at the 83 microJy level of MIPS GTO image. We find that Blue LBGs have a Spectral Energy Distribution similar to high redshift Lyman Break Galaxies. Finally, the dust-to-FUV ratio of this sample is compared with similar ratios at z=1 and z~2. This work suggests an evolution (decrease) of the dust-to-FUV ratio with the redshift.

We present the rest-frame 8 micron luminosity function (LF) at redshifts z=1 and ~2, computed from Spitzer 24 micron-selected galaxies in the GOODS fields over an area of 291 sq. arcmin. Using classification criteria based on X-ray data and IRAC colours, we identify the AGN in our sample. The rest-frame 8 micron LF for star-forming galaxies at redshifts z=1 and ~2 have the same shape as at z~0, but with a strong positive luminosity evolution. The number density of star-forming galaxies with log_{10}(nu L_nu(8 micron))>11 increases by a factor >250 from redshift z~0 to 1, and is basically the same at z=1 and ~2. The resulting rest-frame 8 micron luminosity densities associated with star formation at z=1 and ~2 are more than four and two times larger than at z~0, respectively. We also compute the total rest-frame 8 micron LF for star-forming galaxies and AGN at z~2 and show that AGN dominate its bright end, which is well-described by a power-law. Using a new calibration based on Spitzer star-forming galaxies at 0<z<0.6 and validated at higher redshifts through stacking analysis, we compute the bolometric infrared (IR) LF for star-forming galaxies at z=1 and ~2. We find that the respective bolometric IR luminosity densities are (1.2+/-0.2) x 10^9 and (6.6^{+1.2}_{-1.0}) x 10^8 L_sun Mpc^{-3}, in agreement with previous studies within the error bars. At z~2, around 90% of the IR luminosity density associated with star formation is produced by luminous and ultraluminous IR galaxies (LIRG and ULIRG), with the two populations contributing in roughly similar amounts. Finally, we discuss the consistency of our findings with other existing observational results on galaxy evolution.

We determine the rest-frame 8 micron luminosity function of type I quasars over the redshift range 1<z<5. Our sample consists of 292 24 micron sources brighter than 1 mJy selected from 7.17 square degrees of the Spitzer Space Telescope MIPS survey of the NOAO Deep Wide-Field Survey Bootes field. The AGN and Galaxy Evolution Survey (AGES) has measured redshifts for 270 of the R<21.7 sources and we estimate that the contamination of the remaining 22 sources by stars and galaxies is low. We are able to select quasars missed by ultra-violet excess quasar surveys, including reddened type I quasars and 2.2<z<3.0 quasars with optical colors similar to main sequence stars. We find reddened type I quasars comprise 20% of the type I quasar population. Nonetheless, the shape, normalization, and evolution of the rest-frame 8 micron luminosity function is comparable to that of quasars selected from optical surveys. The 8 micron luminosity function of type I quasars is well approximated by a power-law with index -2.75(+/-0.14). We directly measure the peak of the quasar space density to be at z=2.6(+/-0.3).

As part of the "All Wavelength Extended Groth Strip International Survey" (AEGIS), we describe the panchromatic characterization of an X-ray luminous active galactic nucleus (AGN) in a merging galaxy at z=1.15. This object is detected at infrared (8mic, 24mic, 70mic, 160mic), submillimeter (850mic) and radio wavelengths, from which we derive a bolometric luminosity L_bol ~ 9x10^12 Lsol. We find that the AGN clearly dominates the hot dust emission below 40mic but its total energetic power inferred from the hard X-rays is substantially less than the bolometric output of the system. About 50% of the infrared luminosity is indeed produced by a cold dust component that probably originates from enshrouded star formation in the host galaxy. In the context of a coeval growth of stellar bulges and massive black holes, this source might represent a ``transition'' object sharing properties with both quasars and luminous starbursts. Study of such composite galaxies will help address how the star formation and disk-accretion phenomena may have regulated each other at high redshift and how this coordination may have participated to the build-up of the relationship observed locally between the masses of black holes and stellar spheroids.

Aims: We aim at detecting and determining the properties of the host galaxy of the dark GRB 050223. Methods: We use VLT optical/NIR images coupled to Swift X-ray positioning, and optical spectra of the host galaxy to measure its properties. Results: We find a single galaxy within the Swift error box of GRB 050223. It is located at z = 0.584 and its luminosity is L ~ 0.4 L*. Emission lines in the galaxy spectrum imply an intrinsic SFR > 7 Msun/yr, and a large extinction A_V > 2 mag within it. We also detect absorption lines, which reveal an underlying stellar population with an age between 40 Myr and 1.5 Gyr. Conclusions: The identification of a host galaxy with atypical properties using only the X-ray transient suggests that a bias may be present in the former sample of host galaxies. Dust obscuration together with intrinsic faintness are the most probable causes for the darkness of this burst.

We report the detection of rest--frame 6.2 and 7.7 \micron emission features arising from Polycyclic Aromatic Hydrocarbons (PAH) in the Spitzer/IRS spectrum of an infrared-luminous Lyman break galaxy at z=3.01. This is currently the highest redshift galaxy where these PAH emission features have been detected. The total infrared luminosity inferred from the MIPS 24 \micron and radio flux density is 2$\times10^{13}$ L$_{\odot}$, which qualifies this object as a so--called hyperluminous infrared galaxy (HyLIRG). However, unlike local HyLIRGs which are generally associated with QSO/AGNs and have weak or absent PAH emission features, this HyLIRG has very strong 6.2 and 7.7 \micron PAH emission. We argue that intense star formation dominates the infrared emission of this source, although we cannot rule out the presence of a deeply obscured AGN. This LBG appears to be a distorted system in the HST ACS F606W and F814W images, possibly indicating that a significant merger or interaction is driving the large IR luminosity.

We present Spitzer observations of thirteen z~6 quasars using the Infrared Array Camera (IRAC) and Multiband Imaging Photometer for Spitzer (MIPS). All the quasars except SDSS J000552.34-000655.8 (SDSS J0005-0006) were detected with high S/N in the four IRAC channels and the MIPS 24um band, while SDSS J0005-0006 was marginally detected in the IRAC 8.0um band, and not detected in the MIPS 24um band. We find that most of these quasars have prominent emission from hot dust as evidenced by the observed 24um fluxes. Their spectral energy distributions (SEDs) are similar to those of low-redshift quasars at rest-frame 0.15-3.5 um, suggesting that accretion disks and hot-dust structures for these sources already have reached maturity. However, SDSS J0005-0006 has an unusual SED that lies significantly below low-redshift SED templates at rest-frame 1 and 3.5 um, and thus shows a strong near-IR (NIR) deficit and no hot-dust emission. Type I quasars with extremely small NIR-to-optical flux ratios like SDSS J0005-0006 are not found in low-redshift quasar samples, indicating that SDSS J0005-0006 has different dust properties at high redshift. We combine the Spitzer observations with X-ray, UV/optical, mm/submm and radio observations to determine bolometric luminosities for all the quasars. We find that the four quasars with central black-hole mass measurements have Eddington ratios of order unity.

In this the first of a series of Letters, we present a description of the panchromatic data sets that have been acquired in the Extended Groth Strip region of the sky. Our survey, the All-wavelength Extended Groth strip International Survey (AEGIS), is intended to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray (0.5 - 10 keV), GALEX ultraviolet (1200 - 2500 Angstrom), CFHT/MegaCam Legacy Survey optical (3600 - 9000 Angstroms), CFHT/CFH12K optical (4500 - 9000 Angstroms), Hubble Space Telescope/ACS optical (4400 - 8500 Angstroms), Palomar/WIRC near-infrared (1.2 - 2.2 microns), Spitzer/IRAC mid-infrared (3.6 - 8.0 microns), Spitzer/MIPS far-infrared (24 - 70 microns), and VLA radio continuum (6 - 20 cm). In addition, this region of the sky has been targeted for extensive spectroscopy using the DEIMOS spectrograph on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage.

We study the distribution of Eddington luminosity ratios, L_bol/L_edd, of active galactic nuclei (AGNs) discovered in the AGN and Galaxy Evolution Survey (AGES). We combine H-beta, MgII, and CIV line widths with continuum luminosities to estimate black hole (BH) masses in 407 AGNs, covering the redshift range z~0.3-4 and the bolometric luminosity range L_bol~10^45-10^47 erg/s. The sample consists of X-ray or mid-infrared (24 micron) point sources with optical magnitude R<=21.5 mag and optical emission line spectra characteristic of AGNs. For the range of luminosity and redshift probed by AGES, the distribution of estimated Eddington ratios is well described as log-normal with a peak at L_bol/L_edd ~ 1/4 and a dispersion of 0.3 dex. Since additional sources of scatter are minimal, this dispersion must account for contributions from the scatter between estimated and true BH mass and the scatter between estimated and true bolometric luminosity. Therefore, we conclude that: (1) neither of these sources of error can contribute more than \~0.3 dex rms; and (2) the true Eddington ratios of optically luminous AGNs are even more sharply peaked. Because the mass estimation errors must be smaller than ~0.3 dex, we can also investigate the distribution of Eddington ratios at fixed BH mass. We show for the first time that the distribution of Eddington ratios at fixed BH mass is peaked, and that the dearth of AGNs at a factor ~10 below Eddington is real and not an artifact of sample selection. These results provide strong evidence that supermassive BHs gain most of their mass while radiating close to the Eddington limit, and they suggest that the fueling rates in luminous AGNs are ultimately determined by BH self-regulation of the accretion flow rather than galactic scale dynamical disturbances.

We have obtained deep spectroscopic observations of several nearby gamma-ray burst (GRB) host galaxies revealing for the first time the presence of Wolf-Rayet (WR) stars and numerous O stars located in rich and compact clusters or star forming regions. Surprisingly, high spatial resolution imaging shows that the GRBs and the associated supernovae did not occur in these regions, but several hundreds of parsec away. Considering various scenarios for GRB progenitors, we do not find any simple explanation of why they should be preferentially born in regions with low stellar densities. All the examined GRBs and associated SNe have occurred 400 to 800 pc from very high density stellar environments including large numbers of WR stars. Such distances can be travelled through at velocities of 100 km/s or larger, assuming the travel time to be the typical life time of WR stars. It leads us to suggest that GRB progenitors may be runaway massive stars ejected from compact massive star clusters. The ejection from such super star clusters may lead to a spin-up of these stars, producing the loss of the hydrogen and/or helium envelopes leading to the origin of the type Ibc supernovae associated with GRBs. If this scenario applies tocd text/Sc all GRBs, it provides a natural explanation of the very small fraction of massive stars that emit a GRB at the end of their life. An alternative to this scenario could be a binary origin for GRBs, but this still requires an explanation of why it would preferentially occur in low stellar density regions.

We present a brief discussion of the evolution of disk galaxy stellar masses, sizes, rotation velocities, and star formation rates over the last eight billion years. Recent observations have failed to detect significant evolution in the stellar mass Tully-Fisher relation, stellar mass-size relation, and the stellar mass function of disk galaxies. Yet, most z<1 star formation is in disks, and this star formation would be expected to drive a rapid growth of the total stellar mass (and therefore mass function) of disks in the last eight billion years. Such a build-up is not seen; instead, a rapid build-up in the total stellar mass in non-star-forming spheroid-dominated galaxies is observed. Large numbers of disk-dominated galaxies are systematically shutting off their star formation and building up a spheroid (or losing a disk) in the epoch 0<z<1.

We investigate the role of the luminous infrared galaxy (LIRG) and ultra-luminous infrared galaxy (ULIRG) phases in the evolution of Ks-selected galaxies and, in particular, Extremely Red Galaxies (ERGs). With this aim, we compare the properties of a sample of 2905 Ks<21.5 (Vega mag) galaxies in the GOODS/CDFS with the sub-sample of those 696 sources which are detected at 24 microns. We find that LIRGs constitute 30% of the galaxies with stellar mass M>1x10^{11} Msun assembled at redshift z=0.5. A minimum of 65% of the galaxies with M>2.5x10^{11} Msun at z~2-3 are ULIRGs at those redshifts. 60% of the ULIRGs in our sample have the characteristic colours of ERGs. Conversely, 40% of the ERGs with stellar mass M>1.3x10^{11} Msun at 1.5<z<2.0 and a minimum of 52% of those with the same mass cut at 2.0<z<3.0 are ULIRGs. The average optical/near-IR properties of the massive ERGs at similar redshifts that are identified with ULIRGs and that are not have basically no difference, suggesting that both populations contain the same kind of objects in different phases of their lives. LIRGs and ULIRGs have an important role in galaxy evolution and mass assembly, and, although they are only able to trace a fraction of the massive (M>1x10^{11} Msun) galaxies present in the Universe at a given time, this fraction becomes very significant (>50%) at redshifts z>~2.

We have studied the infrared properties of the brightest cluster galaxies (BCGs) located in the cores of X-ray-luminous clusters at 0.15 < z < 0.35. The majority of the BCGs are not particularly infrared-luminous compared with other massive early-type galaxies, suggesting that the cluster environment has little influence on the infrared luminosities of the BCGs. The exceptions, however, are the BCGs in the three X-ray-brightest clusters in the sample, A1835, Z3146, and A2390. These BCGs have a prominent far-infrared peak in their spectral energy distributions (SEDs), and two of them (those in A1835 and Z3146) can be classified as luminous infrared galaxies (LIRGs: L_{IR} > 10^{11} L_{sun}). Although radio AGNs are found to be prevalent among the BCGs, the infrared luminosities of these three BCGs, judged from the infrared SED signatures, are likely to be powered by star formation. Considering the overall trend that clusters with shorter radiative gas cooling times harbor more infrared-luminous BCGs, the enhanced star formation may be caused by the cooling cluster gas accreting onto the BCGs.

We are currently undertaking a Spitzer GTO program to image ~30 massive lensing clusters at moderate redshift with both IRAC and MIPS. By taking advantage of the gravitatinoal lensing power of these clusters, we will study the population of faint galaxies that are below the nominal Spitzer detection limits. Here, we present a few examples of our science programs.

We present a spectroscopic survey using the MMT/Hectospec fiber spectrograph of 24 micron sources selected with the Spitzer Space Telescope in the Spitzer First Look Survey. We report 1296 new redshifts for 24 micron sources, including 599 with f(24micron) > 1 mJy. Combined with 291 additional redshifts for sources from the Sloan Digital Sky Survey (SDSS), our observing program was highly efficient and is ~90% complete for i' < 21 mag and f(24micron) > 1 mJy, and is 35% complete for i' < 20.5 mag and 0.3 mJy < f(24micron) < 1 mJy. Our Hectospec survey includes 1078 and 168 objects spectroscopically classified as galaxies and QSOs, respectively. Combining the Hectospec and SDSS samples, we find 24 micron-selected galaxies to z < 0.98 and QSOs to z < 3.6, with mean redshifts of <z(gal)> = 0.27 and <z(QSO)> =1.1. As part of this publication, we include the redshift catalogs and the reduced spectra; these are also available online (http://mips.as.arizona.edu/~papovich/fls) and through the NASA/IPAC Infrared Science Archive (http://irsa.ipac.caltech.edu).

Aims. We quantify the contributions of 24um galaxies to the Far-Infrared (FIR) Background at 70 and 160um. We provide new estimates of the Cosmic Infrared Background (CIB), and compare it with the Cosmic Optical Background (COB). Methods. Using Spitzer data at 24, 70 and 160um in three deep fields, we stacked more than 19000 MIPS 24um sources with S24>60uJy at 70 and 160um, and measured the resulting FIR flux densities. Results. This method allows a gain up to one order of magnitude in depth in the FIR. We find that the Mid-Infrared (MIR) 24um selected sources contribute to more than 70% of the CIB at 70 and 160um. This is the first direct measurement of the contribution of MIR-selected galaxies to the FIR CIB. Galaxies contributing the most to the total CIB are thus z~1 luminous infrared galaxies, which have intermediate stellar masses. We estimate that the CIB will be resolved at 0.9 mJy at 70 and 3 mJy at 160um. By combining the extrapolation of the 24um source counts below 60uJy, with 160/24 and 70/24 colors as measured with the stacking analysis, we obtain lower limits of 7.1+/-1.0 and 13.4+/-1.7 nW/m2/sr for the CIB at 70 and 160um, respectively. Conclusions. The MIPS surveys have resolved more than three quarters of the MIR and FIR CIB. By carefully integrating the Extragalactic Background Light (EBL) SED, we also find that the CIB has the same brightness as the COB, around 24 nW/m2/sr. The EBL is produced on average by 115 infrared photons for one visible photon. Finally, the galaxy formation and evolution processes emitted a brightness equivalent to 5% of the primordial electromagnetic background (CMB).

We present high resolution 890 micron images of two 20 mJy submillimeter galaxies, SMMJ123711+622212 and MIPS~J142824.0+352619, obtained using the Submillimeter Array (SMA). Using submillimeter interferometric observations with an angular resolution of 2.5", the coordinates of these high redshift sources are determined with an accuracy of 0.2". The new SMA data on SMMJ123711+622212 reveal an unresolved submm source offset to the east by 0.8" from an optical galaxy found in deep HST images, suggesting either a large galaxy with a dusty central region, or an interacting galaxy system. The SMA image of hyper-luminous (LFIR = 3.2 x 10^{13} Lsun) source MIPS~J142824.0+352619 provides a firm upper limit to the source size of < 1.2. This constraint provides evidence that the foreground lens is only weakly affecting the observed high FIR luminosity.

We determine the contribution of AGN to the mid-IR emission of luminous infrared galaxies (LIRGs) at z>0.6 by measuring the mid-IR dust continuum slope of 20,039 mid-IR sources. The 24 micron sources are selected from a Spitzer/MIPS survey of the NOAO Deep Wide-Field Survey Bo\"otes field and have corresponding 8 micron data from the IRAC Shallow Survey. There is a clear bimodal distribution in the 24 micron to 8 micron flux ratio. The X-ray detected sources fall within the peak corresponding to a flat spectrum in nufnu, implying that it is populated by AGN-dominated LIRGs, whereas the peak corresponding to a higher 24 micron to 8 micron flux ratio is likely due to LIRGs whose infrared emission is powered by starbursts. The 24 micron emission is increasingly dominated by AGN at higher 24 micron flux densities (f_24): the AGN fraction of the z>0.6 sources increases from ~9% at f_24 ~ 0.35 mJy to 74+/-20% at f_24 ~ 3 mJy in good agreement with model predictions. Deep 24 micron, small area surveys, like GOODS, will be strongly dominated by starburst galaxies. AGN are responsible for ~ 3-7% of the total 24 micron background.

We report on IRAC-4.5mic, IRAC-8.0mic and MIPS-24mic deep observations of 16 Gamma-Ray Burst (GRBs) host galaxies performed with the Spitzer Space Telescope, and we investigate in the thermal infrared the presence of evolved stellar populations and dust-enshrouded star-forming activity associated with these objects. Our sample is derived from GRBs that were identified with sub-arcsec localization between 1997 and 2001, and only a very small fraction (~20%) of the targeted sources is detected down to f_4.5mic ~3.5microJy and f_24mic ~85microJy (3sigma). This likely argues against a population dominated by massive and strongly-starbursting (i.e., SFR > ~100 Msol/yr) galaxies as it has been recently suggested from submillimeter/radio and optical studies of similarly-selected GRB hosts. Furthermore we find evidence that some GRBs do not occur in the most infrared-luminous regions -- hence the most actively star-forming environments -- of their host galaxies. Should the GRB hosts be representative of all star-forming galaxies at high redshift, models of infrared galaxy evolution indicate that > ~50% of GRB hosts should have f_24mic > ~100microJy. Unless the identification of GRBs prior to 2001 was prone to strong selection effects biasing our sample against dusty galaxies, we infer in this context that the GRBs identified with the current techniques can not be directly used as unbiased probes of the global and integrated star formation history of the Universe.

We report on the first mid-infrared observations of 16 GRB host galaxies performed with the Spitzer Space Telescope, and investigate the presence of evolved stellar populations and dust-enshrouded star-forming activity associated with GRBs. Only a very small fraction of our sample is detected by Spitzer, which is not consistent with recent works suggesting the presence of a GRB host population dominated by massive and strongly-starbursting galaxies (SFR > ~100 Msol/yr). Should the GRB hosts be representative of star-forming galaxies at high redshift, models of galaxy evolution indicate that >~50% of GRB hosts would be easily detected at the depth of our mid-infrared observations. Unless our sample suffers from a strong observational bias which remains to be understood, we infer in this context that the GRBs identified with the current techniques can not be directly used as unbiased probes of the global and integrated star formation history of the Universe.

We present low-resolution (64 < R < 124) mid-infrared (8--38 micron) Spitzer/IRS spectra of two z~1.3 ultraluminous infrared galaxies (LFIR~10^13) discovered in a Spitzer/MIPS survey of the Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). MIPS J142824.0+352619 is a bright 160 micron source with a large infrared-to-optical flux density ratio and a possible lensing amplification of <~10. The 6.2, 7.7, 11.3, and 12.8 micron PAH emission bands in its IRS spectrum indicate a redshift of z~1.3. The large equivalent width of the 6.2 micron PAH feature indicates that at least 50% of the mid-infrared energy is generated in a starburst, an interpretation that is supported by a large [NeII]/[NeIII] ratio and a low upper limit on the X-ray luminosity. SST24 J142827.19+354127.71 has the brightest 24 micron flux (10.55 mJy) among optically faint (R > 20) galaxies in the NDWFS. Its mid-infrared spectrum lacks emission features, but the broad 9.7 micron silicate absorption band places this source at z~1.3. Given this redshift, SST24 J142827.19+354127.71 has among the largest rest-frame 5 micron luminosities known. The similarity of its SED to those of known AGN-dominated ULIRGs and its lack of either PAH features or large amounts of cool dust indicate that the powerful mid-infrared emission is dominated by an active nucleus rather than a starburst. Our results illustrate the power of the IRS in identifying massive galaxies in the ``redshift desert'' and in discerning their power sources. Because they are bright, MIPS J142824.0+352619 (pending future observations to constrain its lensing amplification) and SST24 J142827.19+354127.71 are useful z>1 templates of a high luminosity starburst and AGN, respectively.

We stack Spitzer 24 micron images for ~7000 galaxies with 0.1<z<1 in the Chandra Deep Field South to probe the thermal dust emission in low-luminosity galaxies over this redshift range. Through stacking, we can detect mean 24 micron fluxes that are more than an order of magnitude below the individual detection limit. We find that the correlations for low and moderate luminosity galaxies between the average L_IR/L_UV and rest-frame B-band luminosity, and between the star formation rate (SFR) and L_IR/L_UV, are similar to those in the local Universe. This verifies that oft-used assumption in deep UV/optical surveys that the dust obscuration-SFR relation for galaxies with SFR < 20 solar mass per year varies little with epoch. We have used this relation to derive the cosmic IR luminosity density from z=1 to z=0.1. The results also demonstrate directly that little of the bolometric luminosity of the galaxy population arises from the faint end of the luminosity function, indicating a relatively flat faint-end slope of the IR luminosity function with a power law index of 1.2+-0.3.

We investigate the nature of a sample of 92 Spitzer/MIPS 24 micron selected galaxies in the CDFS, showing power law-like emission in the Spitzer/IRAC 3.6-8 micron bands. The main goal is to determine whether the galaxies not detected in X-rays (47% of the sample) are part of the hypothetical population of obscured AGN not detected even in deep X-ray surveys. The majority of the IR power-law galaxies are ULIRGs at z>1, and those with LIRG-like IR luminosities are usually detected in X-rays. The optical to IR spectral energy distributions (SEDs) of the X-ray detected galaxies are almost equally divided between a BLAGN SED class (similar to an optically selected QSO) and a NLAGN SED (similar to the BLAGN SED but with an obscured UV/optical continuum). A small fraction of SEDs resemble warm ULIRG galaxies (e.g., Mrk231). Most galaxies not detected in X-rays have SEDs in the NLAGN+ULIRG class as they tend to be optically fainter, and possibly more obscured. Moreover, the IR power-law galaxies have SEDs significantly different from those of high-z (z_sp>1) IR (24 micron) selected and optically bright (VVDS I_AB<=24) star-forming galaxies whose SEDs show a very prominent stellar bump at 1.6 micron. The galaxies detected in X-rays have 2-8 keV rest-frame luminosities typical of AGN. The galaxies not detected in X-rays have global X-ray to mid-IR SED properties that make them good candidates to contain IR bright X-ray absorbed AGN. If all these sources are actually obscured AGN, we would observe a ratio of obscured to unobscured 24 micron detected AGN of 2:1, whereas models predict a ratio of up to 3:1. Additional studies using Spitzer to detect X-ray-quiet AGN are likely to find more such obscured sources.

We investigate the properties of massive galaxies at z=1-3.5 using HST observations, ground-based near-IR imaging, and Spitzer Space Telescope observations at 3-24 micron. We identify 153 distant red galaxies (DRGs) with J-K > 2.3 mag (Vega) in the southern GOODS field. This sample is approximately complete in stellar mass for passively evolving galaxies above 10^11 solar masses and z < 3. The galaxies identified by this selection are roughly split between objects whose optical and near-IR rest-frame light is dominated by evolved stars combined with ongoing star formation, and galaxies whose light is dominated by heavily reddened starbursts. Very few of the galaxies (< 10%) have no indication of current star formation. Using SFR estimates that include the reradiated IR emission, the DRGs at z=1.5-3 with stellar masses > 10^11 solar masses have specific SFRs (SFRs per unit stellar mass) ranging from 0.2 to 10 Gyr^-1, with a mean value of ~2.4 Gyr^-1. The DRGs with stellar masses > 10^11 solar masses and 1.5 < z < 3 have integrated specific SFRs greater the global value over all galaxies. In contrast, we find that galaxies at z = 0.3-0.75 with these stellar masses have integrated specific SFRs less than the global value, and more than an order of magnitude lower than that for massive DRGs at z = 1.5-3. At z < 1, lower-mass galaxies dominate the overall cosmic mass assembly. This suggests that the bulk of star formation in massive galaxies occurs at early cosmic epochs and is largely complete by z~1.5. [Abridged]

Spectra have been obtained with the Infrared Spectrograph on the Spitzer Space Telescope for 18 optically faint sources (R > 23.9,mag) having f(nu) (24um) > 1.0,mJy and having radio detections at 20 cm to a limit of 115 microJy. Sources are within the Spitzer First Look Survey. Redshifts are determined for 14 sources from strong silicate absorption features (12 sources) or strong PAH emission features (2 sources), with median redshift of 2.1. Results confirm that optically faint sources of ~1 mJy at 24um are typically at redshifts z ~ 2, verifying the high efficiency in selecting high redshift sources based on extreme infrared to optical flux ratio, and indicate that 24um sources which also have radio counterparts are not systematically different than samples chosen only by their infrared to optical flux ratios. Using the parameter q = log[f(nu)(24um)/f(nu)(20 cm)] 17 of the 18 sources observed have values of 0<q<1, in the range expected for starburst-powered sources, but only a few of these show strong PAH emission as expected from starbursts, with the remainder showing absorbed or power-law spectra consistent with an AGN luminosity source. This confirms previous indications that optically faint Spitzer sources with f(nu)(24um) > 1.0mJy are predominately AGN and represent the upper end of the luminosity function of dusty sources at z ~ 2. Based on the characteristics of the sources observed so far, we predict that the nature of sources selected at 24um will change for f(nu)(24um) < 0.5 mJy to sources dominated primarily by starbursts.

We present deep Ks<21.5 (Vega) identifications, redshifts and stellar masses for most of the sources composing the bulk of the 24 micron background in the GOODS/CDFS. Our identified sample consists of 747 Spitzer/MIPS 24 micron objects, and includes ~94% of all the 24 micron sources in the GOODS-South field which have fluxes Snu(24)>83 microJy (the 80% completeness limit of the Spitzer/GTO 24 micron catalog). 36% of our galaxies have spectroscopic redshifts (mostly at z<1.5) and the remaining ones have photometric redshifts of very good quality, with a median of |dz|=|zspec-zphot|/(1+zspec)=0.02. We find that MIPS 24 micron galaxies span the redshift range z~0-4, and that a substantial fraction (28%) lie at high redshifts z>1.5. We determine the existence of a bump in the redshift distribution at z~1.9, indicating the presence of a significant population of galaxies with PAH emission at these redshifts. Massive (M>10^11 Msun) star-forming galaxies at redshifts 2<z<3 are characterized by very high star-formation rates (SFR>500 Msun/yr), and some of them are able to construct a mass of 10^10-10^11 Msun in a single burst lifetime (~0.01-0.1 Gyr). At lower redshifts z<2, massive star-forming galaxies are also present, but appear to be building their stars on long timescales, either quiescently or in multiple modest burst-like episodes. At redshifts z~1-2, the ability of the burst-like mode to produce entire galaxies in a single event is limited to some lower (M<7x10^10 Msun) mass systems, and it is basically negligible at z<1. Our results support a scenario where star-formation activity is differential with assembled stellar mass and redshift, and where the relative importance of the burst-like mode proceeds in a down-sizing way from high to low redshifts. (abridged)

Using the 16$\mu$m peakup imager on the Infrared Spectrograph (IRS) on Spitzer, we present a serendipitous survey of 0.0392 deg$^{2}$ within the area of the NOAO Deep Wide Field Survey in Bootes. Combining our results with the available Multiband Imaging Photometer for Spitzer (MIPS) 24$\mu$m survey of this area, we produce a catalog of 150 16$\mu$m sources brighter than 0.18 mJy (3$\sigma$) for which we derive measures or limits on the 16/24$\mu$m colors. Such colors are especially useful in determining redshifts for sources whose mid infrared spectra contain strong emission or absorption features that characterize these colors as a function of redshift. We find that the 9.7$\mu$m silicate absorption feature in Ultraluminous Infrared Galaxies (ULIRGs) results in sources brighter at 16$\mu$m than at 24$\mu$m at z $\sim$ 1--1.8 by at least 20%. With a threshold flux ratio of 1.2, restricting our analysis to $>5\sigma$ detections at 16$\mu$m, and using a $3\sigma$ limit on 24$\mu$m non-detections, the number of silicate-absorbed ULIRG candidates is 36. This defines a strong upper limit of $\sim$920 sources deg$^{-2}$, on the population of silicate-absorbed ULIRGs at z $\sim$ 1--1.8. This source count is about half of the total number of sources predicted at z $\sim$ 1--2 by various phenomenological models. We note that the high 16/24$\mu$m colors measured cannot be reproduced by any of the mid-IR spectral energy distributions assumed by these models, which points to the strong limitations currently affecting our phenomenological and theoretical understanding of infrared galaxy evolution.

Using the SHARC-II camera at the Caltech Submillimeter Observatory to obtain 350 micron images of sources detected with the MIPS instrument on Spitzer, we have discovered a remarkable object at z=1.325+/-0.002 with an apparent Far-Infrared luminosity of 3.2(+/-0.7) x 10^13 Lsun. Unlike other z>1 sources of comparable luminosity selected from mid-IR surveys, MIPS J142824.0+352619 lacks any trace of AGN activity, and is likely a luminous analog of galaxies selected locally by IRAS, or at high redshift in the submillimeter. This source appears to be lensed by a foreground elliptical galaxy at z=1.034, although the amplification is likely modest (~10). We argue that the contribution to the observed optical/Near-IR emission from the foreground galaxy is small, and hence are able to present the rest-frame UV through radio Spectral Energy Distribution of this galaxy. Due to its unusually high luminosity, MIPS J142824.0+352619 presents a unique chance to study a high redshift dusty starburst galaxy in great detail.

We combine optical morphologies and photometry from HST, redshifts from Keck, and mid-infrared luminosities from Spitzer for an optically selected sample of~800 galaxies in GOODS-N to track morphology evolution of infrared luminous galaxies (LIRGs) since redshift z=1. We find a 50% decline in the number of LIRGs from z~1 to lower redshift, in agreement with previous studies. In addition, there is evidence for a morphological evolution of the populations of LIRGs. Above z=0.5, roughly half of all LIRGs are spiral, the peculiar/irregular to spiral ratio is ~0.7, and both classes span a similar range of L_{IR} and M_B. At low-z, spirals account for one-third of LIRGs, the peculiar to spiral fraction rises to 1.3, and for a given M_B spirals tend to have lower IR luminosity than peculiars. Only a few percent of LIRGs at any redshift are red early-type galaxies. For blue galaxies (U-B < 0.2), M_B is well correlated with log(L_{IR}) with an RMS scatter (about a bivariate linear fit) of ~0.25 dex in IR luminosity. Among blue galaxies that are brighter than M_B = -21, 75% are LIRGs, regardless of redshift. These results can be explained by a scenario in which at high-z, most large spirals experience an elevated star formation rate as LIRGs. Gas consumption results in a decline of LIRGs, especially in spirals, to lower redshifts.

We analyze a sample of ~2600 MIPS/Spitzer 24mic sources brighter than ~80muJy and located in the Chandra Deep Field South to characterize the evolution of the comoving infrared (IR) energy density of the Universe up to z~1. Using published ancillary optical data we first obtain a nearly complete redshift determination for the 24mic objects associated with R<24 counterparts at z<1. We find that the 24mic population at 0.5<z<1 is dominated by ``Luminous Infrared Galaxies'' (i.e., 10^11 L_sol < L_IR < 10^12 L_sol), the counterparts of which appear to be also luminous at optical wavelengths and tend to be more massive than the majority of optically-selected galaxies. We finally derive 15mic and total IR luminosity functions (LFs) up to z~1. In agreement with the previous results from ISO and SCUBA and as expected from the MIPS source number counts, we find very strong evolution of the contribution of the IR-selected population with lookback time. Pure evolution in density is firmly excluded by the data, but we find considerable degeneracy between strict evolution in luminosity and a combination of increases in both density and luminosity (L*_IR prop. to (1+z)^{3.2_{-0.2}^{+0.7}}, Phi*_IR prop. to (1+z)^{0.7_{-0.6}^{+0.2}}). Our results imply that the comoving IR energy density of the Universe evolves as (1+z)^(3.9+/-0.4) up to z~1 and that galaxies luminous in the infrared (i.e., L_IR > 10^11 L_IR) are responsible for 70+/-15% of this energy density at z~1. Taking into account the contribution of the UV luminosity evolving as (1+z)^~2.5, we infer that these IR-luminous sources dominate the star-forming activity beyond z~0.7. The uncertainties affecting these conclusions are largely dominated by the errors in the k-corrections used to convert 24mic fluxes into luminosities.

The combination of both contributions from the observed UV emission and the absorbed radiations reprocessed in the infrared represents the ideal approach to constrain the activity of massive star formation in galaxies. Using recent results from GALEX and Spitzer, we compare the evolutions of the UV and IR energy densities with redshift as well as their contributions to the star formation history at 0<z<1. We find that the comoving IR luminosity is characterized by a much faster evolution than seen in the UV. Our results also indicate that ~70% of the star-forming activity at z~1 is produced by the so-called IR-luminous sources (L_IR > 10^11 L_sol).

Two images of Cassiopeia A obtained at 24 micrometer with the Spitzer Space Telescope over a one year time interval show moving structures outside the shell of the supernova remnant to a distance of more than 20 arcmin. Individual features exhibit apparent motions of 10 to 20 arcsec per year, independently confirmed by near-infrared observations. The observed tangential velocities are at roughly the speed of light. It is likely that the moving structures are infrared echoes, in which interstellar dust is heated by the explosion and by flares from the compact object near the center of the remnant.

We present ISOCAM 15 micron and MIPS 24 micron photometry of a sample of 16 distant Luminous Infrared Galaxies (LIRGs) characterized by a median luminosity L(IR) 2x10^11 Lsol and redshift z = 0.7 (distributed from z = 0.1 to 1.2). While some sources display 24/15 micron flux ratios also consistent with a featureless continuum dominating their mid-infrared (MIR) spectral energy distributions (SEDs), the presence of prominent emission features such as the Polycyclic Aromatic Hydrocarbons is clearly required to explain the observed colors for more than half of the sample. As a result, a general good agreement is observed between the data and predictions from the local starburst-dominated SEDs that have been used so far to constrain IR galaxy evolution. This is consistent with the star-forming nature of LIRGs derived from previous works, even though our approach cannot rule out the dominance of an AGN in some cases. Our study also supports the possibility of tracing the total IR luminosity of distant galaxies (up to z ~ 1) from their MIR emission.

We present a first analysis of deep 24 micron observations with the Spitzer Space Telescope of a sample of nearly 1500 galaxies in a thin redshift slice, 0.65<z<0.75. We combine the infrared data with redshifts, rest-frame luminosities, and colors from COMBO-17, and with morphologies from Hubble Space Telescope images collected by the GEMS and GOODS projects. To characterize the decline in star-formation rate (SFR) since z~0.7, we estimate the total thermal infrared (IR) luminosities, SFRs, and stellar masses for the galaxies in this sample. At z~0.7, nearly 40% of intermediate and high-mass galaxies (with stellar masses >2x10^10 solar masses) are undergoing a period of intense star formation above their past-averaged SFR. In contrast, less than 1% of equally-massive galaxies in the local universe have similarly intense star formation activity. Morphologically-undisturbed galaxies dominate the total infrared luminosity density and SFR density: at z~0.7, more than half of the intensely star-forming galaxies have spiral morphologies, whereas less than \~30% are strongly interacting. Thus, a decline in major-merger rate is not the underlying cause of the rapid decline in cosmic SFR since z~0.7. Physical properties that do not strongly affect galaxy morphology - for example, gas consumption and weak interactions with small satellite galaxies - appear to be responsible.

We have surveyed a field covering 9.0 degrees^2 within the NOAO Deep Wide-Field Survey region in Bootes with the Multiband Imaging Photometer on the Spitzer Space Telescope (SST) to a limiting 24 um flux density of 0.3 mJy. Thirty one sources from this survey with F(24um) > 0.75 mJy which are optically very faint (R > 24.5 mag) have been observed with the low-resolution modules of the Infrared Spectrograph on SST. Redshifts derived primarily from strong silicate absorption features are reported here for 17 of these sources; 10 of these are optically invisible (R > 26 mag), with no counterpart in B_W, R, or I. The observed redshifts for 16 sources are 1.7 < z < 2.8. These represent a newly discovered population of highly obscured sources at high redshift with extreme infrared to optical ratios. Using IRS spectra of local galaxies as templates, we find that a majority of the sources have mid-infrared spectral shapes most similar to ultraluminous infrared galaxies powered primarily by AGN. Assuming the same templates also apply at longer wavelengths, bolometric luminosities exceed 10^13 L(solar).

We have combined a survey at 24 microns to 0.3 mJy with the Multiband Imaging Photometer (MIPS) on the Spitzer Space Telescope, a 20 cm A-array VLA survey covering 0.5 deg^2, and an existing 172 ks Chandra X-Ray Observatory exposure to investigate the nature of optically faint radio and X-ray sources in the NOAO Deep Wide-Field Survey (NDWFS) in Bootes. We find little overlap between the radio and infrared selected populations: only 9% of the infrared sources are detected in the radio and only 33% of the radio sources are detected in the infrared. Thirty-six (10%) of the 377 compact radio sources lack optical counterparts in the NDWFS Bw, R, & I images. We refer to these objects as optically invisible radio sources (OIRS). Only four of the thirty-one OIRSs observed with MIPS are detected at 24 microns. Comparisons of the radio and infrared properties of the OIRSs with various galaxy spectral energy distributions demonstrate that most of these sources are powered by AGN rather than starbursts. Similarly, eleven X-ray sources observed by both MIPS and the VLA are classified as optically invisible X-ray sources (OIXS). None are detected at 24 microns or 20 cm. All seven OIXSs detected in Chandra's 0.5-2 keV band have infrared to X-ray flux ratios consistent with their luminosity being dominated by an unobscured AGN. From these results we conclude that both the optically invisible radio and X-ray source populations are primarily AGN, relatively unaffected by dust and most likely at z > 1. No OIRSs are detected in X-ray emission and no OIXSs are detected at 20 cm. However, given the wide range in radio and X-ray properties of known AGN and the size of our samples, this lack of overlap does not necessarily imply AGN source populations.

Using deep observations of the Chandra Deep Field South obtained with MIPS at 24mic, we present our preliminary estimates on the evolution of the infrared (IR) luminosity density of the Universe from z=0 to z~1. We find that a pure density evolution of the IR luminosity function is clearly excluded by the data. The characteristic luminosity L_IR* evolves at least by (1+z)^3.5 with lookback time, but our monochromatic approach does not allow us to break the degeneracy between a pure evolution in luminosity or an evolution in both density and luminosity. Our results imply that IR luminous systems (L_IR > 10^11 L_sol) become the dominant population contributing to the comoving IR energy density beyond z~0.5-0.6. The uncertainties affecting our measurements are largely dominated by the poor constraints on the spectral energy distributions that are used to translate the observed 24mic flux into luminosities.

We report the detection of a z~7 galaxy strongly lensed by the massive galaxy cluster Abell 2218 (z=0.175) at 3.6 and 4.5 um using the Spitzer Observatory and at 1.1 um using the Hubble Space Telescope. The new data indicate a refined photometric redshift in the range of 6.6-6.8 depending on the presence of Ly-alpha emission. The spectral energy distribution is consistent with having a significant Balmer break, suggesting that the galaxy is in the poststarburst stage with an age of at least ~50 Myr and quite possibly a few hundred Myr. This suggests the possibility that a mature stellar population is already in place at such a high redshift. Compared with typical Lyman break galaxies at z~3-4, the stellar mass is an order of magnitude smaller (~10^{9} Msun), but the specific star formation rate (star formation rate/M_{star}) is similarly large (> 10^{-9} yr^{-1}), indicating equally vigorous star-forming activity.

New surveys with the Spitzer space telescope identify distant star-forming and active galaxies by their strong emission at far-infrared wavelengths, which provides strong constraints on these galaxies' bolometric energy. Using early results from Spitzer surveys at 24 micron, we argue that the faint sources correspond to the existence of a population of infrared-luminous galaxies at z > 1 that are not expected from predictions based on previous observations from ISO and IRAS. Combining Spitzer images with deep ground-based optical and Hubble Space Telescope imaging, we discuss the properties of galaxies selected at 24 micron in the region of the Chandra Deep Field South, including redshift and morphological distributions. Galaxies with z < 1 constitute roughly half of the faint 24 micron sources. Infrared-luminous galaxies at these redshifts span a wide variety of normal to strongly interacting/merging morphologies, which suggests that a range of mechanisms produce infrared activity. Large-area, joint surveys between Spitzer and HST are needed to understand the complex relation between galaxy morphology, structure, environment and activity level, and how this evolves with cosmic time. We briefly discuss strategies for constructing surveys to maximize the legacy of these missions.

We derive galaxy source counts at 70 and 160 microns using the Multiband Imaging Photometer for Spitzer (MIPS) to map the Chandra Deep Field South (CDFS) and other fields. At 70 microns, our observations extend upwards about 2 orders of magnitude in flux density from a threshold of 15 mJy, and at 160 microns they extend about an order of magnitude upward from 50 mJy. The counts are consistent with previous observations on the bright end. Significant evolution is detected at the faint end of the counts in both bands, by factors of 2-3 over no-evolution models. This evolution agrees well with models that indicate most ofthe faint galaxies lie at redshifts between 0.7 and 0.9. The new Spitzer data already resolve about 23% of the Cosmic Far Infrared Background at 70 microns and about 7% at 160 microns.

We present Spitzer IRAC (3.6, 4.5, 5.8, 8.0 um) and MIPS (24 um) observations of the SCUBA submillimeter sources and uJy VLA radio sources in a 5'x5' area in the Lockman Hole East region. Out of the ~40 SCUBA/VLA sources in the field, Spitzer counterparts were detected for nearly all except for the few low-weight SCUBA detections. We show that the majority (80-90%) of the detected sources are cold (i.e., starburst-like) infrared-luminous galaxies (L_{IR} > 10^{11} Lsun) at redshift 0.5<z<3.5, whose star-formation rate density (SFRD) is comparable to that of the optically-selected star-forming galaxies.

We investigate the luminous X-ray sources in the Lockman Hole (LH) and the Extended Groth Strip (EGS) detected at 24microns using MIPS and also with IRAC on board Spitzer. We assemble optical/infrared spectral energy distributions (SEDs) for 45 X-ray/24micron sources in the EGS and LH. Only about 1/4 of the hard X-ray/24micron sources show pure type 1 AGN SEDs. More than half of the X-ray/24micron sources have stellar-emission-dominated or obscured SEDs, similar to those of local type 2 AGN and spiral/starburst galaxies. One-third of the sources detected in hard X-rays do not have a 24micron counterpart. Two such sources in the LH have SEDs resembling those of S0/elliptical galaxies. The broad variety of SEDs in the optical-to-Spitzer bands of X-ray selected AGN means that AGN selected according to the behavior in the optical/infrared will have to be supplemented by other kinds of data (e.g., X-ray) to produce unbiased samples of AGN.

We present preliminary results on 24micron detections of luminous infrared galaxies at z>1 with the Multiband Imaging Photometer for Spitzer (MIPS). Observations were performed in the Lockman Hole and the Extended Groth Strip (EGS), and were supplemented by data obtained with the Infrared Array Camera (IRAC) between 3 and 9microns. The positional accuracy of ~2arcsec for most MIPS/IRAC detections provides unambiguous identifications of their optical counterparts. Using spectroscopic redshifts from the Deep Extragalactic Evolutionary Probe survey, we identify 24micron sources at z>1 in the EGS, while the combination of the MIPS/IRAC observations with $BVRIJHK$ ancillary data in the Lockman Hole also shows very clear cases of galaxies with photometric redshifts at 1<z<2.5. The observed 24micron fluxes indicate infrared luminosities greater than 10^11 L_sol, while the data at shorter wavelengths reveal rather red and probably massive (M>=M*) galaxy counterparts. It is the first time that this population of luminous objects is detected up to z~2.5 in the infrared. Our work demonstrates the ability of the MIPS instrument to probe the dusty Universe at very high redshift, and illustrates how the forthcoming Spitzer deep surveys will offer a unique opportunity to illuminate a dark side of cosmic history not explored by previous infrared experiments.

IRAC imaging of a 4'7x4'7 area in the Lockman Hole detected over 400 galaxies in the IRAC 3.6 micron and 4.5 micron bands, 120 in the 5.8 micron, and 80 in the 8 micron bandin 30 minutes of observing time. Color-color diagrams suggest that about half of these galaxies are at redshifts 0.6<z<1.3 with about a quarter at higher redshifts (z>1.3). We also detect IRAC counterparts for 6 of the 7 SCUBA sources and all 9 XMM sources in this area. The detection of the counterparts of the SCUBA sources and galaxies at z>1.3 demonstrates the ability of IRAC to probe the universe at very high redshifts.

We examine the 24 micron to X-ray color of 157 X-ray--selected AGN as a function of X-ray obscuration and optical classification in the Chandra Deep Field South. The sample consists of the Chandra hard--band detections with hard flux >10^-15 erg/s/cm^2. A deep 24 micron mosaic obtained with Spitzer provides mid--infrared fluxes for the sample. Since obscured AGN locally have higher 24 micron / 2-8 keV ratios than unobscured AGN, and since X-ray background models predict a large population of obscured AGN, we expect to find many X-ray--hard, IR--bright AGN. Instead, we find that the 24 micron to X-ray flux ratio does not depend on X-ray hardness in the full sample, nor does it differ between narrow--line and broad--line AGN. We identify 5 nearly Compton--thick AGN, and find they have similar 24 micron to X-ray flux ratios compared to the full sample. We consider AGN in the narrow redshift spikes at z~0.7; for these AGN, there is some evidence that the flux ratio increases with X-ray hardness. The redshift slice also shows an odd trend, which is also prominent in the full sample: a group of X-ray--hard AGN with very low 24 micron to X-ray flux ratios.

We present an updated phenomenological galaxy evolution model to fit the Spitzer 24, 70 and 160 microns number counts as well as all the previous mid and far infrared observations. Only a minor change of the co-moving luminosity density distribution in the previous model (Lagache, Dole, Puget 2003), combined with a slight modification of the starburst template spectra mainly between 12 and 30 microns, are required to fit all the data available. We show that the peak in the MIPS 24 micron counts is dominated by galaxies with redshift between 1 and 2, with a non negligible contribution from the z>2 galaxies (~30% at S=0.2 mJy). The very close agreement between the model and number counts at 15 and 24 microns strikingly implies that (1) the PAHs (Policyclic Aromatic Hydrocarbons) features remain prominent in the redshift band 0.5 to 2.5 and (2) the IR energy output has to be dominated by ~3 10^11 Lo to ~3 10^12 Lo galaxies from redshift 0.5 to 2.5. Combining Spitzer with the Infrared Space Observatory (ISO) deep cosmological surveys gives for the first time an unbiased view of the infrared Universe from z=0 to z=2.5.

We use the source counts measured with the Multiband Imaging Photometer for Spitzer (MIPS) at 24, 70, and 160 microns to determine the 5-sigma confusion limits due to extragalactic sources: 56 micro-Jy, 3.2 and 40 mJy at 24, 70 and 160 microns, respectively. We also make predictions for confusion limits for a number of proposed far infrared missions of larger aperture (3.5 to 10m diameter).

We present sub-millimeter statistical detections of galaxies discovered in the 5'x5' Spitzer Early Release Observations (to 4-15 microJy 5 sigma at 3.6-8 microns, 170 microJy at 24 microns) through a stacking analysis of our reanalysed SCUBA 8mJy survey maps, and a Spitzer identification of a new sub-millimeter point source in the 8mJy survey region. For sources detected at 5.8 or 8 microns (154 and 111 sources respectively), we detect positive skews in the sub-millimeter flux distributions at 99.2-99.8% confidence using Kolmogorov-Smirnov tests, at both 850 microns and 450 microns. We also marginally detect the Spitzer 24 micron galaxies at 850 microns at 97% confidence, and place limits on the mean sub-millimeter fluxes of the 3.6 and 4.5 micron sources. Integrating the sub-millimeter fluxes of the Spitzer populations, we find the 5.8 micron galaxies contribute 0.12 +/- 0.05 nW/m^2/sr to the 850 micron background, and 2.4 +/- 0.7 nW/m^2/sr to the 450 micron background; similar contributions are made by the 8 micron-selected sample. We infer that the populations dominating the 5.8 and 8 micron extragalactic background light also contribute around a quarter of the 850 micron background and the majority of the 450 micron background.

Galaxy source counts in the infrared provide strong constraints on the evolution of the bolometric energy output from distant galaxy populations. We present the results from deep 24 micron imaging from Spitzer surveys, which include approximately 50,000 sources to an 80% completeness of 60 uJy. The 24 micron counts rapidly rise at near-Euclidean rates down to 5 mJy, increase with a super-Euclidean rate between 0.4 - 4 mJy, and converge below 0.3 mJy. The 24 micron counts exceed expectations from non-evolving models by a factor >10 at 0.1 mJy. The peak in the differential number counts corresponds to a population of faint sources that is not expected from predictions based on 15 micron counts from ISO. We argue that this implies the existence of a previously undetected population of infrared-luminous galaxies at z ~ 1-3. Integrating the counts to 60 uJy, we derive a lower limit on the 24 micron background intensity of 1.9 +/- 0.6 nW m-2 sr-1 of which the majority (~ 60%) stems from sources fainter than 0.4 mJy. Extrapolating to fainter flux densities, sources below 60 uJy contribute 0.8 {+0.9/-0.4} nW m-2 sr-1 to the background, which provides an estimate of the total 24 micron background of 2.7 {+1.1/-0.7} nW m-2 sr-1.

We investigate Extremely Red Objects (EROs) using near- and mid-infrared observations in five passbands (3.6 to 24 micron) obtained from the Spitzer Space Telescope, and deep ground-based R and K imaging. The great sensitivity of the IRAC camera allows us to detect 64 EROs in only 12 minutes of IRAC exposure time, by means of an R-[3.6] color cut (analogous to the traditional red R-K cut). A pure infrared K-[3.6] red cut detects a somewhat different population and may be more effective at selecting z > 1.3 EROs. We find 17% of all galaxies detected by IRAC at 3.6 or 4.5 micron to be EROs. These percentages rise to about 40% at 5.8 micron, and about 60% at 8.0 micron. We utilize the spectral bump at 1.6 micron to divide the EROs into broad redshift slices using only near-infrared colors (2.2/3.6/4.5 micron). We conclude that two-thirds of all EROs lie at redshift z > 1.3. Detections at 24 micron imply that at least 11% of 0.6 < z < 1.3 EROs and at least 22% of z > 1.3 EROs are dusty star-forming galaxies.

The polarization curve of GRB 020813 is discussed and compared to different models for the structure, evolution and magnetisation properties of the jet and the interstellar medium onto which the fireball impacts. GRB 020813 is best suited for this kind of analysis for the smoothness of its afterglow light curve, ensuring the applicability of current models. The polarization dataset allows us to rule out the standard GRB jet, in which the energy and Lorentz factor have a well defined value inside the jet opening angle and the magnetic field is generated at the shock front. We explore alternative models finding that a structured jet or a jet with a toroidal component of the magnetic field can fit equally well the polarization curve. Stronger conclusions cannot be drawn due to the incomplete sampling of the polarization curve. A more dense sampling, especially at early times, is required to pin down the structure of the jet and the geometry of its magnetic field.

We analyze cross-correlation functions between Gamma-Ray Burst (GRB) hosts and surrounding galaxies. We have used data obtained with the Very Large Telescope at Cerro Paranal (Chile), as well as public Hubble Space Telescope data. Our results indicate that Gamma-Ray Burst host galaxies do not reside in high galaxy density environments. Moreover, the host-galaxy cross-correlations show a relatively low amplitude. Our results are in agreement with the cross-correlation function between star-forming galaxies and surrounding objects in the HDF-N.

We present evidence for the dramatically different morphology between the rest frame UV and 7micron mid-IR emission of VV114 and Arp299, two nearby (z~0) violently interacting infrared luminous galaxies (LIRGs). Nearly all LIRGs are interacting systems and it is currently accepted that they dominate the IR emission at z>1. Luminous IR galaxies located at z=1-2 could easily be detected as unresolved sources in deep optical/near-IR ground based surveys, as well as in upcoming 24micron surveys with the Space Infrared Telescope Facility. We demonstrate that the spatial resolution of these surveys will result in blending of the emission from unresolved interacting components. An increased scatter will thus be introduced in the observed optical to mid-IR colors of these galaxies, leading to a systematic underestimation of their dust content.

We present mid-infrared spectro-imaging (5 - 16 microns) observations of the infrared luminous interacting system Arp 299 (=Mrk171 =IC694+NGC3690) obtained with the ISOCAM instrument aboard ISO. Our observations show that nearly 40% of the total emission at 7 and 15 microns is diffuse, originating from the interacting disks of the galaxies. Moreover, they indicate the presence of large amounts of hot dust in the main infrared sources of the system and large extinctions toward the nuclei. While the observed spectra have an overall similar shape, mainly composed of Unidentified Infrared Bands (UIB) in the short wavelength domain, a strong continuum at ~ 13 microns and a deep silicate absorption band at 10 microns, their differences reveal the varying physical conditions of each component. For each source, the spectral energy distribution (SED) can be reproduced by a linear combination of a UIB "canonical" spectral template and a hot dust continuum due to a 230-300 K black body, after independently applying an extinction correction to both of them. We find that the UIB extinction does not vary much throughout the system (A_V ~ 5 mag) suggesting that most UIBs originate from less enshrouded regions. IC694 appears to dominate the infrared emission of the system and our observations support the interpretation of a deeply embedded nuclear starburst located behind an absorption of about 40 mag. The central region of NGC3690 displays a hard radiation field characterized by a [NeIII]/[NeII] ratio > 1.8. It also hosts a strong continuum from 5 to 16 microns which can be explained as thermal emission from a deeply embedded (A_V ~ 60 mag) compact source, consistent with the mid-infrared signature of an active galactic nucleus (AGN), and in agreement with recent X-ray findings.

Polarization measurements for the optical counterpart to GRB 021004 are presented and discussed. Our observations were performed with the TNG and the VLT-UT3 (Melipal) during the first and fourth night after the gamma-ray burst discovery. We find robust evidence of temporal evolution of the polarization, which is therefore, at least partially, intrinsic to the optical transient. We do not find convincing evidence of wavelength dependence for the intrinsic polarization of the transient, in agreement with current polarization models for optical afterglows. We discuss the role of dust, both in our galaxy and in the host, in modifying the transmitted polarization vector, showing how a sizable fraction of the observed polarized flux is due to Galactic selective extinction, while it is not possible to single out any clear contribution from dust in the host galaxy. We discuss how our data compare to those obtained by different groups showing that a two-component model is required to describe the complete dataset. This is not surprising given the complex lightcurve of GRB 021004.

We present K-band imaging observations of ten Gamma-Ray Burst (GRB) host galaxies. We compare their observed and absolute K magnitudes as well as their R-K colours with those of other distant sources detected in various optical, near-infrared, mid-infrared and submillimeter deep surveys. We find that the GRB host galaxies, most of them lying at 0.5<z<1.5, exhibit very blue colours, comparable to those of the faint blue star-forming sources at high redshift. They are sub-luminous in the K-band, suggesting a low stellar mass content. We do not find any GRB hosts harbouring R- and K-band properties similar to those characterizing the luminous infrared/submillimeter sources and the extremely red starbursts. Should GRBs be regarded as an unbiased probe of star-forming activity, this lack of luminous and/or reddened objects among the GRB host sample might reveal that the detection of GRB optical afterglows is likely biased toward unobscured galaxies. It would moreover support the idea that a large fraction of the optically-dark GRBs occur within dust-enshrouded regions of star formation. On the other hand, our result might also simply reflect intrinsic properties of GRB host galaxies experiencing a first episode of very massive star formation and characterized by a rather weak underlying stellar population. Finally, we compute the absolute B magnitudes for the whole sample of GRB host galaxies with known redshifts and detected at optical wavelengths. We find that the latter appear statistically even less luminous than the sub-luminous blue sources which mostly contributed to the B-band light emitted at high redshift. This indicates that the formation of GRBs could be favoured in particular systems with very low luminosities and, therefore, low metallicities. (Abridged)

We present VLT spectroscopic observations of Gamma-Ray Burst (GRB) 990705 host galaxy and highlight the benefits provided by the prompt-phase features of GRBs to derive the redshifts of the latter. In the host spectrum, we indeed detect an emission feature which we attribute to the [OII] 3726/3729 Angstrom doublet and derive an unambiguous redshift z=0.8424+/-0.0002 for this galaxy. This is in full agreement with the value z~0.86+/-0.17 previously derived using a transient absorption edge discovered in the X-ray spectrum of GRB990705. This burst is therefore the first GRB for which a reliable redshift was derived from the prompt phase emission itself, as opposed to redshift determinations performed using putative host galaxy emission lines or intestellar absorption lines in the GRB afterglows. Deep and high resolution images of the host of GRB990705 with the HST/STIS camera reveal that the burst occured in a nearly face-on Sc spiral galaxy typical of disk-dominated systems at 0.75<z<1. Assuming a cosmology with Ho=65km/s/Mpc, Omega_m=0.3 and Omega_lambda=0.7, we derive an absolute B magnitude M_B=-21.75 for this galaxy and a star formation rate SFR ~ 5-8~Msol/yr. Finally, we discuss the implications of using X-ray transient features to derive GRB redshifts with larger burst samples, and especially examine the case of short and dark-long GRBs.

Cosmological Gamma-Ray Bursts (GRBs) offer a unique perspective to probe the evolution of distant galaxies. We discuss one of the multiple benefits of this approach, i.e. the detection of transient edges in the GRB prompt phase emission. These absorption features can be used to directly derive the redshift of GRBs and their host galaxies without the need of any optical spectroscopic follow-up.

Polarization measurements for the optical counterpart to GRB 020405 are presented and discussed. Our observations were performed with the VLT-UT3 (Melipal) during the second and third night after the gamma-ray burst discovery. The polarization degree (and the position angle) appears to be constant between our two observations at a level around (1.5 - 2)%. The polarization can be intrinsic but it is not possible to unambiguously exclude that a substantial fraction of it is induced by dust in the host galaxy.

There is now an increasing number of evidence supporting the idea that the cosmic Gamma-ray Bursts (GRBs) originate from the collapse of massive stars in distant star-forming galaxies. Because GRBs are likely detectable up to very high redshift, and because the gamma-rays are not attenuated by intervening columns of gas and dust, these phenomena thus offer a unique perspective to probe the star formation in the early Universe independently of the biases associated with dust extinction. We report on the R-K colors of the GRB host galaxies. We find that they are very blue sources, with colors comparable to those of the faint blue star-forming galaxies at high redshift.