We present a multi-wavelength study (from X-ray to millimetre) of the nearby low-luminosity active galactic nucleus (LLAGN) NGC 7213. We combine the information from the different bands to characterize the source in terms of contribution from the AGN and the host-galaxy interstellar medium (ISM). This approach allows us to provide a coherent picture of the role of the AGN and its impact, if any, on the star formation and molecular gas properties of the host galaxy. We focused our study on archival ALMA Cycle 1 observations, where the CO(2-1) emission line has been used as a tracer of the molecular gas. Using the 3DBarolo code on ALMA data, we performed the modelling of the molecular gas kinematics traced by the CO(2-1) emission, finding a rotationally dominated pattern. The host-galaxy molecular gas mass was estimated from the integrated CO(2-1) emission line obtained with APEX data, assuming an $\alpha_{CO}$ conversion factor. By using the ALMA data, we would have underestimated the gas masses by a factor $\sim$3, given the filtering out of the large scale emission in interferometric observations. We also performed a complete X-ray spectral analysis on archival observations, revealing a relatively faint and unobscured AGN. The AGN results to be too faint to significantly affect the properties of the host-galaxy, such as star formation activity and molecular gas kinematics and distribution.

A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 \mu m}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy.

The sub-mm sky is a unique window for probing the architecture of the Universe and structures within it. From the discovery of dusty sub-mm galaxies, to the ringed nature of protostellar disks, our understanding of the formation, destruction, and evolution of objects in the Universe requires a comprehensive view of the sub-mm sky. The current generation single-dish sub-mm facilities have shown of the potential for discovery, while interferometers have presented a high resolution view into the finer details. However, our understanding of large-scale structure and our full use of these interferometers is now hampered by the limited sensitivity of our sub-mm view of the universe at larger scales. Thus, now is the time to start planning the next generation of sub-mm single dish facilities, to build on these revolutions in our understanding of the sub-mm sky. Here we present the case for the Atacama Large Aperture Submillimeter Telescope (AtLAST), a concept for a 50m class single dish telescope. We envision AtLAST as a facility operating as an international partnership with a suite of instruments to deliver the transformative science described in many Astro2020 science white papers. A 50m telescope with a high throughput and 1$^\circ$ FoV with a full complement of advanced instrumentation, including highly multiplexed high-resolution spectrometers, continuum cameras and Integral Field Units, AtLAST will have mapping speeds thousands of times greater than any current or planned facility. It will reach confusion limits below $L_*$ in the distant universe and resolve low-mass protostellar cores at the distance of the Galactic Center, providing synergies with upcoming facilities across the spectrum. Located on the Atacama plateau, to observe frequencies un-obtainable by other observatories, AtLAST will enable a fundamentally new understanding of the sub-mm universe at unprecedented depths.

We present a detailed study of a high-redshift iron low-ionization broad absorption line (FeLoBAL) quasar (SDSS1214 at $z = 1.046$), including new interferometric $^{12}$CO $J$=2-1 observations, optical through far-infrared photometry, and mid-infrared spectroscopy. The CO line is well-fit by a single Gaussian centered 40 kms$^{-1}$ away from the systemic velocity and implies a total molecular gas mass of $M_\textrm{gas} = 7.3 \times 10^{10} \textrm{M}_\odot$. The infrared SED requires three components: an active galactic nucleus (AGN) torus, an AGN polar dust component, and a starburst. The starburst dominates the infrared emission with a luminosity of log($L_\textrm{SB}[\textrm{L}_\odot]) = 12.91^{+0.02}_{-0.02}$, implying a star formation rate of about 2000 $\textrm{M}_{\odot}$yr$^{-1}$, the highest known among FeLoBAL quasars. The AGN torus and polar dust components are less luminous, at log($L_\textrm{AGN}[\textrm{L}_\odot]) = 12.36^{+0.14}_{-0.15}$ and log($L_\textrm{dust}[\textrm{L}_\odot]) = 11.75^{+0.26}_{-0.46}$, respectively. If all of the molecular gas is used to fuel the ongoing star formation, then the lower limit on the subsequent duration of the starburst is 40 Myr. We do not find conclusive evidence that the AGN is affecting the CO gas reservoir. The properties of SDSS1214 are consistent with it representing the endpoint of an obscured starburst transitioning through a LoBAL phase to that of a classical quasar.

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

We present ALMA Atacama Compact Array (ACA) 870 micron continuum maps of 28 infrared-bright SDSS quasars with Herschel/SPIRE detections at redshifts 2-4, the largest such sample ever observed with ALMA. The ACA detections are centred on the SDSS coordinates to within 1 arcsec for about 80 per cent of the sample. Larger offsets indicate that the far-infrared (FIR) emission detected by Herschel might come from a companion source. The majority of the objects (about 70 per cent) have unique ACA counterparts within the SPIRE beam down to 3-4 arcsec resolution. Only 30 per cent of the sample shows clear evidence for multiple sources with secondary counterparts contributing to the total 870 micron flux within the SPIRE beam to at least 25 per cent. We discuss the limitations of the data based on simulated pairs of point-like sources at the resolution of the ACA and present an extensive comparison of our findings with recent works on the multiplicities of sub-millimetre galaxies. We conclude that, despite the coarse resolution of the ACA, our data support the idea that, for a large fraction of FIR-bright quasars, the sub-mm emission comes from single sources. Our results suggest that, on average, optically bright quasars with strong FIR emission are not triggered by early-stage mergers but are, instead, together with their associated star formation rates, the outcome of either late-stage mergers or secular processes.

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

We use a semi-empirical model to reproduce the 0.1-10um spectral energy distribution (SED) of a sample of 85 luminous quasars. In the model, the continuum emission from the accretion disk as well as the nebular lines are represented by a single empirical template (disk), where differences in the optical spectral index are reproduced by varying the amount of extinction. The near- and mid-infrared emission of the AGN-heated dust is modelled as the combination of two black-bodies (dust). The model fitting shows that the disk and dust components are remarkably uniform among individual quasars, with differences in the observed SED largely accounted for by varying levels of obscuration in the disk as well as differences in the relative luminosity of the disk and dust components. By combining the disk-subtracted SEDs of the 85 quasars, we generate a template for the 1-10um emission of the AGN-heated dust. Additionally, we use a sample of local Seyfert 1 galaxies with full spectroscopic coverage in the 0.37um to 39um range to demonstrate a method for stitching together spectral segments obtained with different PSF and extraction apertures. We show that the disk and dust templates obtained from luminous quasars also reproduce the optical-to-mid-infrared spectra of local Seyfert 1s when the contribution from the host galaxy is properly subtracted.

A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths.

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

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

Active Galactic Nuclei (AGN) are energetic astrophysical sources powered by accretion onto supermassive black holes in galaxies, and present unique observational signatures that cover the full electromagnetic spectrum over more than twenty orders of magnitude in frequency. The rich phenomenology of AGN has resulted in a large number of different "flavours" in the literature that now comprise a complex and confusing AGN "zoo". It is increasingly clear that these classifications are only partially related to intrinsic differences between AGN, and primarily reflect variations in a relatively small number of astrophysical parameters as well the method by which each class of AGN is selected. Taken together, observations in different electromagnetic bands as well as variations over time provide complementary windows on the physics of different sub-structures in the AGN. In this review, we present an overview of AGN multi-wavelength properties with the aim of painting their "big picture" through observations in each electromagnetic band from radio to gamma-rays as well as AGN variability. We address what we can learn from each observational method, the impact of selection effects, the physics behind the emission at each wavelength, and the potential for future studies. To conclude we use these observations to piece together the basic architecture of AGN, discuss our current understanding of unification models, and highlight some open questions that present opportunities for future observational and theoretical progress.

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map the CO(3-2) and the underlying continuum emissions around the type 1 low-luminosity active galactic nucleus (LLAGN; bolometric luminosity $\lesssim 10^{42}$ erg~s$^{-1}$) of NGC 1097 at $\sim 10$ pc resolution. These observations revealed a detailed cold gas distribution within a $\sim 100$ pc of this LLAGN. In contrast to the luminous Seyfert galaxy NGC 1068, where a $\sim 7$ pc cold molecular torus was recently revealed, a distinctively dense and compact torus is missing in our CO(3-2) integrated intensity map of NGC 1097. Based on the CO(3-2) flux, the gas mass of the torus of NGC 1097 would be a factor of $\gtrsim 2-3$ less than that found for NGC 1068 by using the same CO-to-H$_2$ conversion factor, which implies less active nuclear star formation and/or inflows in NGC 1097. Our dynamical modeling of the CO(3-2) velocity field implies that the cold molecular gas is concentrated in a thin layer as compared to the hot gas traced by the 2.12 $\mu$m H$_2$ emission in and around the torus. Furthermore, we suggest that NGC 1097 hosts a geometrically thinner torus than NGC 1068. Although the physical origin of the torus thickness remains unclear, our observations support a theoretical prediction that geometrically thick tori with high opacity will become deficient as AGNs evolve from luminous Seyferts to LLAGNs.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

We analyse a sample of 85 luminous (log(nuLnu(3\um)/erg s-1)>45.5) quasars with restframe ~2-11 \um spectroscopy from AKARI and Spitzer. Their high luminosity allows a direct determination of the near-infrared quasar spectrum free from host galaxy emission. A semi-empirical model consisting of a single template for the accretion disk and two blackbodies for the dust emission successfully reproduces the 0.1-10 \um spectral energy distributions (SEDs). Excess emission at 1-2 \um over the best-fitting model suggests that hotter dust is necessary in addition to the ~1200 K blackbody and the disk to reproduce the entire near-infrared spectrum. Variation in the extinction affecting the disk and in the relative strength of the disk and dust components accounts for the diversity of individual SEDs. Quasars with higher dust-to-disk luminosity ratios show slightly redder infrared continua and less prominent silicate emission. We find no luminosity dependence in the shape of the average infrared quasar spectrum. We generate a new quasar template that covers the restframe range 0.1-11 \um, and separate templates for the disk and dust components. Comparison with other infrared quasar composites suggests that previous ones are less reliable in the 2-4 \um range. Our template is the first one to provide a detailed view of the infrared emission on both sides of the 4 \um bump.

We explore the relationship between active galactic nuclei and star formation in a sample of 513 optically luminous type 1 quasars up to redshifts of $\sim$4 hosting extremely high star formation rates (SFRs). The quasars are selected to be individually detected by the \textit{Herschel} SPIRE instrument at $> $3$\sigma$ at 250 $\mu$m, leading to typical SFRs of order of 1000 M$_{\odot}$yr$^{-1}$. We find the average SFRs to increase by almost a factor 10 from $z\sim0.5$ to $z\sim3$, mirroring the rise in the comoving SFR density over the same epoch. However, we find that the SFRs remain approximately constant with increasing accretion luminosity for accretion luminosities above 10$^{12}$ L$_{\odot}$. We also find that the SFRs do not correlate with black hole mass. Both of these results are most plausibly explained by the existence of a self-regulation process by the starburst at high SFRs, which controls SFRs on time-scales comparable to or shorter than the AGN or starburst duty cycles. We additionally find that SFRs do not depend on Eddington ratio at any redshift, consistent with no relation between SFR and black hole growth rate per unit black hole mass. Finally, we find that high-ionisation broad absorption line (HiBAL) quasars have indistinguishable far-infrared properties to those of classical quasars, consistent with HiBAL quasars being normal quasars observed along a particular line of sight, with the outflows in HiBAL quasars not having any measurable effect on the star formation in their hosts.

We present global VLBI observations of the 21-cm transition of atomic hydrogen seen in absorption against the radio source J0855+5751. The foreground absorber (SDSS~J085519.05+575140.7) is a dwarf galaxy at $z$ = 0.026. As the background source is heavily resolved by VLBI, the data allow us to map the properties of the foreground HI gas with a spatial resolution of 2pc. The absorbing gas corresponds to a single coherent structure with an extent $>$35pc, but we also detect significant and coherent variations, including a change in the HI optical depth by a factor of five across a distance of $\leq$6pc. The large size of the structure provides support for the Heiles & Troland model of the ISM, as well as its applicability to external galaxies. The large variations in HI optical depth also suggest that caution should be applied when interpreting $T_S$ measurements from radio-detected DLAs. In addition, the distorted appearance of the background radio source is indicative of a strong jet-cloud interaction in its host galaxy. We have measured its redshift ($z$ = 0.54186) using optical spectroscopy on the William Herschel Telescope and this confirms that J0855+5751 is a FRII radio source with a physical extent of $<$1kpc and supports the previous identification of this source as a Compact Symmetric Object. These sources often show absorption associated with the host galaxy and we suggest that both HI and OH should be searched for in J0855+5751.

We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects.

In recent years there has been a paradigm shift from centralised to geographically distributed resources. Individual entities are no longer able to host or afford the necessary expertise in-house, and, as a consequence, society increasingly relies on widespread collaborations. Although such collaborations are now the norm for scientific projects, more technical structures providing support to a distributed scientific community without direct financial or other material benefits are scarce. The network of European ALMA Regional Centre (ARC) nodes is an example of such an internationally distributed user support network. It is an organised effort to provide the European ALMA user community with uniform expert support to enable optimal usage and scientific output of the ALMA facility. The network model for the European ARC nodes is described in terms of its organisation, communication strategies and user support.

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

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

We present an unprecedented measurement of the disc stability and local instability scales in the luminous infrared Seyfert 1 host, NGC7469, based on ALMA observations of dense gas tracers and with a synthesized beam of 165 x 132 pc. While we confirm that non-circular motions are not significant in redistributing the dense interstellar gas in this galaxy, we find compelling evidence that the dense gas is a suitable tracer for studying the origin of its intensely high-mass star forming ring-like structure. Our derived disc stability parameter accounts for a thick disc structure and its value falls below unity at the radii in which intense star formation is found. Furthermore, we derive the characteristic instability scale and find a striking agreement between our measured scale of ~ 180 pc, and the typical sizes of individual complexes of young and massive star clusters seen in high-resolution images.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

We present results on the spectral decomposition of 118 Spitzer Infrared Spectrograph (IRS) spectra from local active galactic nuclei (AGN) using a large set of Spitzer/IRS spectra as templates. The templates are themselves IRS spectra from extreme cases where a single physical component (stellar, interstellar, or AGN) completely dominates the integrated mid-infrared emission. We show that a linear combination of one template for each physical component reproduces the observed IRS spectra of AGN hosts with unprecedented fidelity for a template fitting method, with no need to model extinction separately. We use full probability distribution functions to estimate expectation values and uncertainties for observables, and find that the decomposition results are robust against degeneracies. Furthermore, we compare the AGN spectra derived from the spectral decomposition with sub-arcsecond resolution nuclear photometry and spectroscopy from ground-based observations. We find that the AGN component derived from the decomposition closely matches the nuclear spectrum, with a 1-sigma dispersion of 0.12 dex in luminosity and typical uncertainties of ~0.19 in the spectral index and ~0.1 in the silicate strength. We conclude that the emission from the host galaxy can be reliably removed from the IRS spectra of AGN. This allows for unbiased studies of the AGN emission in intermediate and high redshift galaxies -currently inaccesible to ground-based observations- with archival Spitzer/IRS data and in the future with the Mid-InfraRed Instrument of the James Webb Space Telescope. The decomposition code and templates are available at http://www.denebola.org/ahc/deblendIRS.

We present a comprehensive study of the silicate features at 9.7 and 18 micron of a sample of almost 800 active galactic nuclei (AGN) with available spectra from the Spitzer InfraRed Spectrograph (IRS). We measure the strength of the silicate feature at 9.7 micron, S9.7, before and after subtracting the host galaxy emission from the IRS spectra. The numbers of type 1 and 2 AGN with the feature in emission increase by 20 and 50%, respectively, once the host galaxy is removed, while 35% of objects with the feature originally in absorption exhibit it in even deeper absorption. The peak of S9.7, lambda_peak, has a bimodal distribution when the feature is in emission, with about 65% of the cases showing lambda_peak > 10.2 micron. Silicates can appear in emission in objects with mid-infrared (MIR) luminosity spanning over six orders of magnitude. The derived distributions of the strength of the silicate features at 9.7 and 18 micron provide a solid test bed for modeling the dust distribution in AGN. Clumpiness is needed in order to produce absorption features in unobscured AGN and can also cause the silicates to be in absorption at 9.7 micron and in emission at 18 micron in type 1 sources. We find the `cosmic' silicates of Ossenkopf et al. to be more consistent with the observations than Draine's `astronomical' silicates. Finally, we discuss the possibility of a foreground absorber to explain the deep silicate absorption features in the MIR spectra of some type 2 AGN.

We present an Analytic Model of Intergalactic-medium and GAlaxy evolution since the dark ages. AMIGA is in the spirit of the popular semi-analytic models of galaxy formation, although it does not use halo merger trees but interpolates halo properties in grids that are progressively built. This strategy is less memory-demanding and allows one to start the modeling at redshifts high enough and halo masses low enough to have trivial boundary conditions. The number of free parameters is minimized by making the causal connection between physical processes usually treated as independent from each other, which leads to more reliable predictions. But the strongest points of AMIGA are: i) the inclusion of molecular cooling and metal-poor, population III (Pop III) stars, with the most dramatic feedback, and ii) the accurate follow-up of the temperature and volume filling factor of neutral, singly, and doubly ionized regions, taking into account the distinct halo mass functions in those environments. We find the following general results. Massive Pop III stars determine the IGM metallicity and temperature, and the growth of spheroids and disks is self-regulated by that of massive black holes developed from the remnants of those stars. Yet, the properties of normal galaxies and active galactic nuclei appear to be quite insensitive to Pop III star properties owing to the much higher yield of ordinary stars compared to Pop III stars and the dramatic growth of MBHs when normal galaxies begin to develop, which cause the memory loss of the initial conditions.

The Photometry Data Model (PhotDM) standard describes photometry filters, photometric systems, magnitude systems, zero points and its interrelation with the other IVOA data models through a simple data model. Particular attention is given necessarily to optical photometry where specifications of magnitude systems and photometric zero points are required to convert photometric measurements into physical flux density units.

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

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

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

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

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

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

We present a kinematic analysis of the dense molecular gas in the central 200 parsecs of the nearby galaxy NGC1097, based on Cycle 0 observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA). We use the HCN(4-3) line to trace the densest interstellar molecular gas, and quantify its kinematics, and estimate an inflow rate for the molecular gas. We find a striking similarity between the ALMA kinematic data and the analytic spiral inflow model that we have previously constructed based on ionized gas velocity fields on larger scales. We are able to follow dense gas streaming down to 40 pc distance from the supermassive black hole in this Seyfert 1 galaxy. In order to fulfill marginal stability, we deduce that the dense gas is confined to a very thin disc, and we derive a dense gas inflow rate of 0.09 Msun/yr at 40 pc radius. Combined with previous values from the Ha and CO gas, we calculate a combined molecular and ionized gas inflow rate of 0.2 Msun/yr at 40 pc distance from the central supermassive black hole of NGC1097.

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.

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

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

We present the Seyfert and star formation Activity in the Far-InfraRed (SAFIR) project, a small (15.1h) Herschel guaranteed time proposal performing PACS and SPIRE imaging of a small sample of nearby Seyfert galaxies. This project is aimed at studying the physical nature of the nuclear IR emission by means of multi-component spectral energy distribution (SED) fitting and the star formation properties of AGN hosts, as traced by cold dust. We summarize the results achieved so far and outline the on-going work.

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

The geometry of the dust distribution within the inner regions of Active Galactic Nuclei (AGN) is still a debated issue and relates directly with the AGN unified scheme. Traditionally, models discussed in the literature assume one of two distinct dust distributions in what is believed to be a toroidal region around the Supermassive Black Holes: a continuous distribution, customarily referred to as smooth, and a concentration of dust in clumps or clouds, referred to as clumpy. In this paper we perform a thorough comparison between two of the most popular models in the literature, namely the smooth models by Fritz. et al. 2006 and the clumpy models by Nenkova et al. 2008a, in their common parameters space. Particular attention is paid to the silicate features at ~9.7 and ~18 micron, the width of the infrared bump, the near-infrared index and the luminosity at 12.3 micron, all previously reported as possible diagnostic tools to distinguish between the two dust distributions. We find that, due to the different dust chemical compositions used in the two models, the behaviour of the silicate features at 9.7 and 18 micron is quite distinct between the two models. The width of the infrared bump and the peak of the infrared emission can take comparable values, their distributions do, however, vary. The near-infrared index is also quite different, due partly to the primary sources adopted by the two models. Models with matched parameters do not produce similar SEDs and virtually no random parameter combinations can result in seemingly identical SEDs.

We present new far-infrared (70-500micron) Herschel PACS and SPIRE imaging observations as well as new mid-IR Gemini/T-ReCS imaging (8.7 and 18.3micron) and spectroscopy of the inner Lindblad resonance (ILR) region (R<2.5kpc) of the spiral galaxy NGC1365. We complemented these observations with archival Spitzer imaging and spectral mapping observations. The ILR region of NGC1365 contains a Seyfert 1.5 nucleus and a ring of star formation with an approximate diameter of 2kpc. The strong star formation activity in the ring is resolved by the Herschel/PACS imaging data, as well as by the Spitzer 24micron continuum emission, [NeII]12.81micron line emission, and 6.2 and 11.3micron PAH emission. The AGN is the brightest source in the central regions up to lambda~24micron, but it becomes increasingly fainter in the far-infrared when compared to the emission originating in the infrared clusters (or groups of them) located in the ring. We modeled the AGN unresolved infrared emission with the CLUMPY torus models and estimated that the AGN contributes only to a small fraction (~5%) of the infrared emission produced in the inner ~5kpc. We fitted the non-AGN 24-500micron spectral energy distribution of the ILR region and found that the dust temperatures and mass are similar to those of other nuclear and circumnuclear starburst regions. Finally we showed that within the ILR region of NGC1365 most of the on-going star formation activity is taking place in dusty regions as probed by the 24micron emission.

We report Herschel SPIRE (250, 350, and 500 micron) detections of 32 quasars with redshifts 0.5 < z < 3.6 from the Herschel Multi-tiered Extragalactic Survey. These sources are from a MIPS 24 micron flux-limited sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field permit construction of the rest-frame Spectral Energy Distributions (SEDs)from ultraviolet to the mid-infrared for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show dust temperatures in the range of 25K to 60K, with a mean of 34K. The FIR luminosities range from 10^{11.3} to 10^{13.5} Lsun, qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars observed at ~ 1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3--20 micron, rest-frame), and the bolometric luminosities derived using the 5100 A index may be underestimated for these FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities for quasars with higher near-infrared luminosities.

We compute the exponential disc scalelength for 686 disc galaxies with spectroscopic redshifts out to redshift 5.8 based on Hubble Space Telescope archival data. We compare the results with our previous measurements based on 30000 nearby galaxies from the Sloan Digital Sky Survey. Our results confirm the presence of a dominating exponential component in galaxies out to this redshift. At the highest redshifts, the disc scalelength for the brightest galaxies with absolute magnitude between -24 and -22 is up to a factor 8 smaller compared to that in the local Universe. This observed scalelength decrease is significantly greater than the value predicted by a cosmological picture in which baryonic disc scalelength scales with the virial radius of the dark matter halo.

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.

Mrk938 is a luminous infrared galaxy in the local Universe believed to be the remnant of a galaxy merger. It shows a Seyfert 2 nucleus and intense star formation according to optical spectroscopic observations. We have studied this galaxy using new Herschel far-IR imaging data in addition to archival X-ray, UV, optical, near-IR and mid-IR data. Mid- and far-IR data are crucial to characterise the starburst contribution, allowing us to shed new light on its nature and to study the coexistence of AGN and starburst activity in the local Universe. The decomposition of the mid-IR Spitzer spectrum shows that the AGN bolometric contribution to the mid-IR and total infrared luminosity is small (Lbol(AGN)/LIR~0.02), which agrees with previous estimations. We have characterised the physical nature of its strong infrared emission and constrained it to a relatively compact emitting region of <2kpc. It is in this obscured region where most of the current star formation activity is taking place as expected for LIRGs. We have used Herschel imaging data for the first time to constrain the cold dust emission with unprecedented accuracy. We have fitted the integrated far-IR spectral energy distribution and derived the properties of the dust, obtaining a dust mass of 3x10^7Msun. The far-IR is dominated by emission at 35K, consistent with dust heated by the on-going star formation activity.

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

We report far-infrared (FIR) imaging of the Seyfert 2 galaxy NGC 3081 in the range 70-500 micron, obtained with an unprecedented angular resolution, using the Herschel Space Observatory instruments PACS and SPIRE. The 11 kpc (~70 arcsec) diameter star-forming ring of the galaxy appears resolved up to 250 micron. We extracted infrared (1.6-500 micron) nuclear fluxes, that is active nucleus-dominated fluxes, and fitted them with clumpy torus models, which successfully reproduce the FIR emission with small torus sizes. Adding the FIR data to the near- and mid-infrared spectral energy distribution (SED) results in a torus radial extent of Ro=4(+2/-1) pc, as well as in a flat radial distribution of the clouds (i.e. the q parameter). At wavelengths beyond 200 micron, cold dust emission at T=28+/-1 K from the circumnuclear star-forming ring of 2.3 kpc (~15 arcsec) in diameter starts making a contribution to the nuclear emission. The dust in the outer parts of the galaxy is heated by the interstellar radiation field (19+/-3 K).

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

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^{\theta} (1 + z)^{\zeta}, we find {\theta} = 0.22 +/- 0.08 and {\zeta} = 1.6 +/- 0.4, although there is substantial additional uncertainty in {\zeta} of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

We present a panoramic atlas of Spitzer/IRS spectra of extragalactic sources collected from the recent literature, with value added measurements of their spectral features obtained in a homogeneous and concise manner. The atlas covers the full spectrum of the extragalactic universe and includes star forming galaxies, obscured and unobscured active galaxies, luminous and ultra-luminous infrared galaxies, and hybrid objects. Measured features such as the polycyclic aromatic hydrocarbons, the strength of the silicates in emission or absorption around 9.7 micron, rest-frame monochromatic luminosities or colours, combined with measurements derived from spectral decomposition, are used to establish diagnostics that allow for classification of sources, based on their infrared properties alone. Average templates of the various classes are also derived. The full atlas with the value added measurements and ancillary archival data are publicly available at http://www.denebola.org/atlas, with full references to the original data.

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

We use new mid-infrared (mid-IR) photometry from the Spitzer Space Telescope to study the relations between low-frequency radio luminosity density L_151MHz, mid-IR (12um rest-frame) luminosity L_12um, and optical-emission-line ([OII]) luminosity L_[OII], for a complete sample of z~1 radio galaxies from the 3CRR, 6CE, 6C*, 7CRS and TOOT00 surveys. The narrow redshift span of our sample (0.9<z<1.1) means that it is unbiased to evolutionary effects. We find evidence that these three quantities are positively correlated. The scaling between L_12um and L_[OII] is similar to that seen in other AGN samples, consistent with both L_12um and L_[OII] tracing accretion rate. We show that the positive correlation between L_12um and L_151MHz implies that there is a genuine lack of objects with low values of L_12um at high values of L_151MHz. Given that L_12um traces accretion rate, while L_151MHz traces jet power, this can be understood in terms of a minimum accretion rate being necessary to produce a given jet power. This implies that there is a maximum efficiency with which accreted energy can be chanelled into jet power and that this efficiency is of order unity.

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

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

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

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

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

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

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

We have derived disk scale lengths for 30374 non-interacting disk galaxies in all five SDSS bands. Virtual Observatory methods and tools were used to define, retrieve, and analyse the images for this unprecedentedly large sample classified as disk/spiral galaxies in the LEDA catalogue. Cross correlation of the SDSS sample with the LEDA catalogue allowed us to investigate the variation of the scale lengths for different types of disk/spiral galaxies. We further investigat asymmetry, concentration, and central velocity dispersion as indicators of morphological type, and are able to assess how the scale length varies with respect to galaxy type. We note however, that the concentration and asymmetry parameters have to be used with caution when investigating type dependence of structural parameters in galaxies. Here, we present the scale length derivation method and numerous tests that we have carried out to investigate the reliability of our results. The average r-band disk scale length is 3.79 kpc, with an RMS dispersion of 2.05 kpc, and this is a typical value irrespective of passband and galaxy morphology, concentration, and asymmetry. The derived scale lengths presented here are representative for a typical galaxy mass of $10^{10.8\pm 0.54} \rm{M}_\odot$, and the RMS dispersion is larger for more massive galaxies. Distributions and typical trends of scale lengths have also been derived in all the other SDSS bands with linear relations that indicate the relation that connect scale lengths in one passband to another. Such transformations could be used to test the results of forthcoming cosmological simulations of galaxy formation and evolution of the Hubble sequence.

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

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

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

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

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

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

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

We present an analysis of a large sample of AGN environments at z~1 using stacked Spitzer data at 3.6um. The sample contains type-1 and type-2 AGN in the form of quasars and radio galaxies, and spans a large range in both optical and radio luminosity. We find, on average, that 2 to 3 massive galaxies containing a substantial evolved stellar population lie within a 200-300 kpc radius of the AGN, constituting a >8-sigma excess relative to the field. Secondly, we find evidence for the environmental source density to increase with the radio luminosity of AGN, but not with black-hole mass. This is shown first by dividing the AGN into their classical AGN types, where we see more significant over-densities in the fields of the radio-loud AGN. If instead we dispense with the classical AGN definitions, we find that the source over-density as a function of radio luminosity for all our AGN exhibits a positive correlation. One interpretation of this result is that the Mpc-scale environment is in some way influencing the radio emission that we observe from AGN. This could be explained by the confinement of radio jets in dense environments leading to enhanced radio emission or, alternatively, may be linked to more rapid black-hole spin brought on by galaxy mergers.

The virtual observatory (VO) is a collection of interoperable data archives, tools and applications that together form an environment in which original astronomical research can be carried out. The VO is opening up new ways of exploiting the huge amount of data provided by the ever-growing number of ground-based and space facilities, as well as by computer simulations. This presentation summarises a variety of scientific results spanning various fields of astronomy, obtained thanks to the VO, after highlighting its structure, infrastructure and various capabilities.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

As a part of a Euro-VO research initiative, we have undertaken a programme aimed at studying the scale length of 54909 Sa-Sd spiral galaxies from the SDSS DR6 catalogue. We have retrieved u,g,r,i,z-band images for all galaxies in order to derive the light profiles. We also calculate asymmetry parameters to select non-disturbed disks for which we will derive exponential disk scale lengths. As images in different bands probe different optical depths and stellar populations, it is likely that a derived scale length value should depend on waveband, and our goal is to use the scale length variations with band pass, inclination, galaxy type, redshift, and surface brightness, in order to better understand the nature of spiral galaxies.

(abridged) This paper is the second part of a work investigating the properties of dusty tori in AGN by means of multi-component SED fitting. It focuses on low luminosity, low redshift (z < 0.25) AGN selected among emission line galaxies as well as X-ray, radio and mid-infrared selected type 2 AGN samples from the literature. The available multi-band photometry covers the spectral range from the u-band up to 160 um. The observed SED of each object is fit to a set of multi-component models comprising a stellar component, a high optical depth torus and cold emission from a starburst (SB). The contribution of the various components (stars, torus, SB) is reflected in the position of the objects on the IRAC colour diagram. The comparison of type 1 (as derived from Hatziminaoglou et al. 2008) and type 2 AGN properties is broadly consistent with the Unified Scheme. The estimated ratio between type 2 and type 1 objects is about 2-2.5:1. The AGN accretion-to-infrared luminosity ratio is an indicator of the obscuration of the AGN since it scales down with the covering factor. We find evidence supporting the receding torus paradigm, with the estimated fraction of obscured AGN decreasing with increasing optical luminosity over four orders of magnitude. The average star formation rates are of ~ 10, 40 and 115 Mo/yr for the low-z, type 2 and quasar samples, respectively; but this might simply reflect observational biases, as our quasars are one to two orders of magnitude more luminous than the type 2 AGN. For the majority of objects with 70 and/or 160 um detections an SB component was needed in order to reproduce the data points, implying that the far-infrared emission in AGN arises mostly from star formation; moreover, the SB-to-AGN luminosity ratio shows a slight trend with increasing luminosity.

We present submillimetre observations of the z=3.653 quasar SDSS160705+533558 together with data in the optical and infrared. The object is unusually bright in the far-IR and submm with an IR luminosity of ~10^14 L_sun. We ascribe this luminosity to a combination of AGN and starburst emission, with the starburst forming stars at a rate of a few thousand solar masses per year. Submillimetre Array (SMA) imaging observations with a resolution ~1" show that the submm (850 micron) emission is extended on scales of 10--35kpc and is offset from the optical position by ~10 kpc. This morphology is dissimilar to that found in submm galaxies, which are generally un- or marginally resolved on arcsecond scales, or submm-luminous AGN where the AGN lies at the peak of the submm or molecular emission. The simplest explanation is that the object is in the early stages of a merger between a gas rich galaxy, which hosts the starburst, and a gas-poor AGN-host galaxy, which is responsible for the quasar emission. It is also possible that jet induced star formation might contribute to the unusual morphology.

The radio-quiet neutron star 1E1207.4-5209 has been the target of several XMM-Newton observations, with a total exposure of ~350 ks. The source is located at intermediate galactic latitude (b~10 degrees), i.e. in a sky region with an extremely interesting mix of both galactic and extra-galactic X-ray sources. The aim of our work is to investigate the properties of both the intermediate-latitude galactic and extra-galactic X-ray source populations in the 1E1207.4-5209 field. We performed a coherent analysis of the whole XMM-Newton observation data set to build a catalogue of serendipitous X-ray sources detected with high confidence and to derive information on the source flux, spectra, and time variability. In addition, we performed a complete multi-band (UBVRI) optical coverage of the field with the Wide Field Imager (WFI) of the ESO/MPG 2.2m telescope (La Silla) to search for candidate optical counterparts to the X-ray sources, down to a V-band limiting magnitude of ~24.5. We detected a total of 144 serendipitous X-ray sources. Thanks to the refined X-ray positions and to the WFI observations, we found candidate optical counterparts for most of the X-ray sources in our compilation. For most of the brightest ones we proposed a likely classification based on both the X-ray spectra and the optical colours. Our results indicate that at intermediate galactic latitude the X-ray source population is dominated by the extra-galactic component, but with a significant contribution from the galactic component in the soft energy band, below 2 keV.

AIMS: We present a sample of candidate quasars selected using the KX-technique. The data cover 0.68 deg^2 of the X-ray Multi-Mirror (XMM) Large-Scale Structure (LSS) survey area where overlapping multi-wavelength imaging data permits an investigation of the physical nature of selected sources. METHODS: The KX method identifies quasars on the basis of their optical (R and z') to near-infrared (Ks) photometry and point-like morphology. We combine these data with optical (u*,g'r',i',z') and mid-infrared (3.6-24 micron) wavebands to reconstruct the spectral energy distributions (SEDs) of candidate quasars. RESULTS: Of 93 sources selected as candidate quasars by the KX method, 25 are classified as quasars by the subsequent SED analysis. Spectroscopic observations are available for 12/25 of these sources and confirm the quasar hypothesis in each case. Even more, 90% of the SED-classified quasars show X-ray emission, a property not shared by any of the false candidates in the KX-selected sample. Applying a photometric redshift analysis to the sources without spectroscopy indicates that the 25 sources classified as quasars occupy the interval 0.7 < z < 2.5. The remaining 68/93 sources are classified as stars and unresolved galaxies.

We present results on low-resolution mid-infrared (MIR) spectra of 70 infrared-luminous galaxies obtained with the Infrared Spectrograph (IRS) onboard Spitzer. We selected sources from the European Large Area Infrared Survey (ELAIS) with S15 > 0.8 mJy and photometric or spectroscopic z > 1. About half of the sample are QSOs in the optical, while the remaining sources are galaxies, comprising both obscured AGN and starbursts. We classify the spectra using well-known infrared diagnostics, as well as a new one that we propose, into three types of source: those dominated by an unobscured AGN (QSOs), obscured AGN, and starburst-dominated sources. Starbursts concentrate at z ~ 0.6-1.0 favored by the shift of the 7.7-micron PAH band into the selection 15 micron band, while AGN spread over the 0.5 < z < 3.1 range. Star formation rates (SFR) are estimated for individual sources from the luminosity of the PAH features. An estimate of the average PAH luminosity in QSOs and obscured AGN is obtained from the composite spectrum of all sources with reliable redshifts. The estimated mean SFR in the QSOs is 50-100 Mo yr^-1, but the implied FIR luminosity is 3-10 times lower than that obtained from stacking analysis of the FIR photometry, suggesting destruction of the PAH carriers by energetic photons from the AGN. The SFR estimated in obscured AGN is 2-3 times higher than in QSOs of similar MIR luminosity. This discrepancy might not be due to luminosity effects or selection bias alone, but could instead indicate a connection between obscuration and star formation. However, the observed correlation between silicate absorption and the slope of the near- to mid-infrared spectrum is compatible with the obscuration of the AGN emission in these sources being produced in a dust torus.

We have used far-infrared data from IRAS, ISO, SWIRE, SCUBA and MAMBO to constrain statistically the mean far-infrared luminosities of quasars. Our quasar compilation at redshifts 0<z<6.5 and I-band luminosities -20<I(AB)<-32 is the first to distinguish evolution from quasar luminosity dependence in such a study. We carefully cross-calibrate IRAS against Spitzer and ISO, finding evidence that IRAS 100um fluxes at <1Jy are overestimated by ~30%. We find evidence for a correlation between star formation in quasar hosts and the quasar optical luminosities, varying as SFR proportional to L_opt^(0.44+/-0.07) at any fixed redshift below z=2. We also find evidence for evolution of the mean star formation rate in quasar host galaxies, scaling as (1+z)^(1.6+/-0.3) at z<2 for any fixed quasar I-band absolute magnitude fainter than -28. We find no evidence for any correlation between star formation rate and black hole mass at 0.5<z<4. Our data are consistent with feedback from black hole accretion regulating stellar mass assembly at all redshifts.

We derive the properties of dusty tori in Active Galactic Nuclei (AGN) from the comparison of observed Spectral Energy Distributions (SEDs) of SDSS quasars and a precomputed grid of torus models. The observed SEDs comprise SDSS photometry, 2MASS J, H, and K data, whenever available and mid-Infrared (MIR) data from the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. The adopted model is that of Fritz et al., 2006. The fit is performed by standard chi^2 minimisation, the model however can be multi-component comprising a stellar and a starburst components, whenever necessary. Models with low equatorial optical depth, tau_9.7, were allowed as well as ``traditional'' models with tau_9.7 > 1.0, corresponding to A_V > 22 and the results were compared. Fits using high optical depth tori models only produced dust more compactly distributed than in the configuration where all tau_9.7 models were permitted. Tori with decreasing dust density with the distance from the centre were favoured while there was no clear preference for models with or without angular variation of the dust density. The computed outer radii of the tori are of some tens of parsecs large but can reach, in a few cases, a few hundreds of parsecs. The mass of dust, M_Dust, and infrared luminosity, L_IR, integrated in the wavelength range between 1 and 1000 micron, do not show significant variations with redshift, once the observational biases are taken into account. Objects with 70 micron detections, representing 25% of the sample, are studied separately and the starburst contribution (whenever present) to the IR luminosity can reach, in the most extreme but very few cases, 80%.

We present the SWIRE Photometric Redshift Catalogue, 1025119 redshifts of unprecedented reliability and accuracy. Our method is based on fixed galaxy and QSO templates applied to data at 0.36-4.5 mu, and on a set of 4 infrared emission templates fitted to infrared excess data at 3.6-170 mu. The code involves two passes through the data, to try to optimize recognition of AGN dust tori. A few carefully justified priors are used and are the key to supression of outliers. Extinction, A_V, is allowed as a free parameter. We use a set of 5982 spectroscopic redshifts, taken from the literature and from our own spectroscopic surveys, to analyze the performance of our method as a function of the number of photometric bands used in the solution and the reduced chi^2. For 7 photometric bands the rms value of (z_{phot}-z_{spec})/(1+z_{spec}) is 3.5%, and the percentage of catastrophic outliers is ~1%. We discuss the redshift distributions at 3.6 and 24 mu. In individual fields, structure in the redshift distribution corresponds to clusters which can be seen in the spectroscopic redshift distribution. 10% of sources in the SWIRE photometric redshift catalogue have z >2, and 4% have z>3, so this catalogue is a huge resource for high redshift galaxies. A key parameter for understanding the evolutionary status of infrared galaxies is L_{ir}/L_{opt}, which can be interpreted as the specific star-formation rate for starbursts. For dust tori around Type 1 AGN, L_{tor}/L_{opt} is a measure of the torus covering factor and we deduce a mean covering factor of 40%.

In a previous paper (Gavignaud et al. 2006), we presented the type-1 Active Galactic Nuclei (AGN) sample obtained from the first epoch data of the VIMOS-VLT Deep Survey (VVDS). The sample consists of 130 faint, broad-line AGN with redshift up to z=5 and 17.5< I <24.0, selected on the basis of their spectra. In this paper we present the measurement of the Optical Luminosity Function up to z=3.6 derived from this sample, we compare our results with previous results from brighter samples both at low and at high redshift. Our data, more than one magnitude fainter than previous optical surveys, allow us to constrain the faint part of the luminosity function up to high redshift. By combining our faint VVDS sample with the large sample of bright AGN extracted from the SDSS DR3 (Richards et al., 2006b) and testing a number of different evolutionary models, we find that the model which better represents the combined luminosity functions, over a wide range of redshift and luminosity, is a luminosity dependent density evolution (LDDE) model, similar to those derived from the major X-surveys. Such a parameterization allows the redshift of the AGN space density peak to change as a function of luminosity and explains the excess of faint AGN that we find at 1.0< z <1.5. On the basis of this model we find, for the first time from the analysis of optically selected samples, that the peak of the AGN space density shifts significantly towards lower redshift going to lower luminosity objects. This result, already found in a number of X-ray selected samples of AGN, is consistent with a scenario of "AGN cosmic downsizing", in which the density of more luminous AGN, possibly associated to more massive black holes, peaks earlier in the history of the Universe, than that of low luminosity ones.

We report results from a Spitzer GO-1 program of IRS spectroscopy of a large sample of Luminous Infrared Galaxies and quasars selected from the European Large Area ISO Survey (ELAIS) which have a wide multiwavelength coverage, including ISOCAM, ISOPHOT, IRAC and MIPS (from SWIRE), and optical photometry. We present the sample selection and results from the IRS spectroscopy.

This work focuses on the properties of dusty tori in active galactic nuclei (AGN) derived from the comparison of SDSS type 1 quasars with mid-Infrared (MIR) counterparts and a new, detailed torus model. The infrared data were taken by the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. Basic model parameters are constraint, such as the density law of the graphite and silicate grains, the torus size and its opening angle. A whole variety of optical depths is supported. The favoured models are those with decreasing density with distance from the centre, while there is no clear tendency as to the covering factor, i.e. small, medium and large covering factors are almost equally distributed. Based on the models that better describe the observed SEDs, properties such as the accretion luminosity, the mass of dust, the inner to outer radius ratio and the hydrogen column density are computed. The properties of the tori, as derived fitting the observed SEDs, are independent of the redshift, once observational biases are taken into account.

The limited angular resolution of gamma-ray telescopes prevents a direct identification of the majority of sources detected so far. This is particularly true for the low latitude, probably galactic, ones only 10 % of which has been identified. Most counterparts of the identified low-latitude gamma-ray sources are Isolated Neutron Stars (INS), both radio-loud and radio-quiet (Geminga-like) objects, which are characterised by an extremely high value of the X-ray-to-optical flux ratio f_X/f_opt. Therefore, the systematic X-ray and optical coverage of low-latitude unidentified gamma-ray sources aiming at high f_X/f_opt sources seems one of the most promising ways to spot INS candidate counterparts. Since low latitude sources are heavily affected by the interstellar absorption at both X-ray and optical wavelengths, we have focussed on two middle-latitude, probably galactic, GRO/EGRET sources: 3EG J0616-3310 and 3EG J1249-8330. These two sources, which could belong to a local galactic population, have been selected owing to their relatively good positional accuracy, spectral shape and lack of candidate extragalactic radio counterparts. Here we report on X-ray observations of the two gamma-ray error boxes performed with XMM-Newton and on their optical follow-up carried on with the Wide Field Imager at the ESO/MPG 2.2m telescope. Less than half of the ~300 sources detected by the X-ray coverage have no optical counterparts. Among those, we have selected few interesting sources with f_X/f_opt > 100, which we consider promising INS candidates.

We present the rest-frame optical and infrared colours of a complete sample of 1114 z<0.3 galaxies from the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) and the Sloan Digital Sky Survey (SDSS). We discuss the optical and infrared colours of our sample and analyse in detail the contribution of dusty star-forming galaxies and AGN to optically selected red sequence galaxies. We propose that the optical (g-r) colour and infrared log(L_{24}/L_{3.6}) colour of galaxies in our sample are determined primarily by a bulge-to-disk ratio. The (g-r) colour is found to be sensitive to the bulge-to-disk ratio for disk-dominated galaxies, whereas the log(L_{24}/L_{3.6}) colour is more sensitive for bulge-dominated systems. We identify ~18% (195 sources) of our sample as having red optical colours and infrared excess. Typically, the infrared luminosities of these galaxies are found to be at the high end of star-forming galaxies with blue optical colours. Using emission line diagnostic diagrams, 78 are found to have an AGN contribution, and 117 are identified as star-forming systems. The red (g-r) colour of the star-forming galaxies could be explained by extinction. However, their high optical luminosities cannot. We conclude that they have a significant bulge component. The number densities of optically red star-forming galaxies are found to correspond to ~13% of the total number density of our sample. In addition, these systems contribute ~13% of the total optical luminosity density, and 28% of the total infrared luminosity density of our SWIRE/SDSS sample. These objects may reduce the need for "dry-mergers".

We combine bright XMM data with the Chandra Deep Field South observations in order to explore the behavior of the intrinsic AGN absorption, as a function of redshift and luminosity.Our sample consists of 359 sources selected in the hard 2-8 keV band, spanning the flux range 6\times10^{-16}-$3\times10^{-13} erg s^-1 cm^-2 with a high rate of spectroscopic or photometric redshift completeness (100 and 85 per cent respectively for the Chandra and XMM data. We derive the column density values using X-ray spectral fits. We find that the fraction of obscured AGN falls with increasing luminosity in agreement with previous findings. The fraction of obscured AGN shows an apparent increase at high redshifts (z>2). Simulations show that this effect can be most probably attributed to the fact that at high redshifts the column densities are overestimated.

We construct rest-frame luminosity functions at 3.6, 4.5, 5.8, 8 and 24 microns over the redshift range 0<z<2 for galaxies and 0<z<4 for optical QSOs, using optical and infrared data from the Spitzer Wide-area InfraRed Extragalactic survey. The 3.6 and 4.5 micron galaxy LFs show evidence for moderate positive luminosity evolution up to z~1.5, consistent with the passive ageing of evolved stellar populations. Their comoving luminosity density was found to evolve passively, gradually increasing out to z~0.5-1 but flattening, or even declining, at higher redshift. Conversely, the 24 micron galaxy LF, which is more sensitive to obscured star formation and/or AGN activity, undergoes strong positive evolution, with the derived IR energy density and SFR density proportional to (1+z)^gamma with gamma=4.5^{+0.7}_{-0.6} and the majority of this evolution occurring since z~1. Optical QSOs, however, show positive luminosity evolution in all bands, out to the highest redshifts (3<z<4). Modelling as L* proportional to (1+z)^gamma gave gamma=1.3^{+0.1}_{-0.1} at 3.6 micron, gamma=1.0^{+0.1}_{-0.1} at 4.5 micron and stronger evolution at the longer wavelengths (5.8, 8 and 24 micron), of gamma~3. Comparison of the galaxy LFs to predictions from a semi-analytic model based on CDM indicate that an IMF skewed towards higher mass star formation in bursts compared to locally is preferred. As a result the currently inferred massive star formation rates in distant sub-mm sources may require substantial downwards revision.

We discuss spectral energy distributions, photometric redshifts, redshift distributions, luminosity functions, source-counts and the far infrared to optical luminosity ratio for sources in the SWIRE Legacy Survey. The spectral energy distributions of selected SWIRE sources are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 5%. While further spectroscopic data are needed to confirm this result, the prospect of determining good photometric redshifts for the 2 million extragalactic objects in SWIRE is excellent. The distribution of the different infrared sed types in the L{ir}/L{opt} versus L{ir} plane, where L{ir} and L{opt} are the infrared and optical bolometric luminosities, is discussed. Source-counts at 24, 70 and 160 mu are discussed, and luminosity functions at 3.6 and 24 mu are presented.

This work focuses on the properties of dusty tori in active galactic nuclei (AGN) derived from the comparison of SDSS type 1 quasars with mid-Infrared (MIR) counterparts and a new, detailed torus model. The infrared data were taken by the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. Basic model parameters are constraint, such as the density law of the graphite and silicate grains, the torus size and its opening angle. A whole variety of optical depths is supported. The favoured models are those with decreasing density with distance from the centre, while there is no clear tendency as to the covering factor, ie small, medium and large covering factors are almost equally distributed. Based on the models that better describe the observed SEDs, properties such as the accretion luminosity, the mass of dust, the inner to outer radius ratio and the hydrogen column density are computed.

We constrain the evolution of the galaxy mass and luminosity functions from the analysis of (public) multi-wavelength data in the Chandra Deep Field South (CDFS) area, obtained from the GOODS and other projects, and including very deep high-resolution imaging by HST/ACS. Our reference catalogue of faint high-redshift galaxies, which we have thoroughly tested for completeness and reliability, comes from a deep (S(3.6micron)>1 microJy) image by IRAC on the Spitzer Observatory. These imaging data in the field are complemented with extensive optical spectroscopy by the ESO VLT/FORS2 and VIMOS spectrographs, while deep K-band VLT/ISAAC imaging is also used to derive further complementary statistical constraints and to assist the source identification and SED analysis. We have selected a highly reliable IRAC 3.6micron sub-sample of 1478 galaxies with S(3.6)>10microJy, 47% of which have spectroscopic redshift, while for the remaining objects both COMBO-17 and Hyperz are used to estimate the photometric redshift. This very extensive dataset is exploited to assess evolutionary effects in the galaxy luminosity and stellar mass functions, while luminosity/density evolution is further constrained with the number counts and redshift distributions. The deep ACS imaging allows us to differentiate these evolutionary paths by morphological type, which our simulations show to be reliable at least up to z=1.5 for the two main early- (E/S0) and late-type (Sp/Irr) classes. These data, as well as our direct estimate of the stellar mass function above M=10^(10)M_sun for the spheroidal subclass, consistently evidence a progressive dearth of such objects to occur starting at z=0.7, paralleled by an increase in luminosity. (abridged)

We describe improved modelling of the emission by dust in a toroidal--like structure heated by a central illuminating source within Active Galactic Nuclei (AGN). We chose a simple but realistic torus geometry, a flared disc, and a dust grain distribution function including a full range of grain sizes. The optical depth within the torus is computed in detail taking into account the different sublimation temperatures of the silicate and graphite grains, which solves previously reported inconsistencies in the silicate emission feature in type-1 AGN. We exploit this model to study the spectral energy distributions (SEDs) of 58 extragalactic (both type-1 and type-2) sources using archival optical and infrared (IR) data. We find that both AGN and starburst contributions are often required to reproduce the observed SEDs, although in a few cases they are very well fitted by a pure AGN component. The AGN contribution to the far-IR luminosity is found to be higher in type-1 sources, with all the type-2 requiring a substantial contribution from a circum-nuclear starburst. Our results appear in agreement with the AGN Unified Scheme, since the distributions of key parameters of the torus models turn out to be compatible for type-1 and type-2 AGN. Further support to the unification concept comes from comparison with medium-resolution IR spectra of type-1 AGN by the Spitzer observatory, showing evidence for a moderate silicate emission around 10 \mums which our code reproduces. From our analysis we infer accretion flows in the inner nucleus of local AGN characterized by high equatorial optical depths ($A_V\simeq 100$), moderate sizes ($R_{max}<100 pc$) and very high covering factors ($f\simeq 80$ per cent) on average.

We present ground-based optical and Spitzer infrared imaging observations of the interacting galaxy UGC 10214, the "Tadpole Galaxy" (z = 0.0310), focusing on the star formation activity in the nuclear, disk, spiral arms and tidal tail regions. The major findings of this study are that the Tadpole is actively forming stars in the main disk outside of the nucleus and in the tidal plume, with an estimated mean star formation rate of ~2 to 4 M_sun/yr. The most prominent sites of mid-infrared emission define a "ring" morphology that, combined with the overall morphology of the system, suggest the interaction may belong to the rare class of off-center collisional ring systems that form both shock-induced rings of star formation and tidal plumes. The nuclear emission is solely powered by older stars, with little evidence for ongoing star formation at the center of the Tadpole. Extra-nuclear star formation accounts for >50% of the total star formation in the disk and spiral arms, featuring infrared-bright 'hot spots' that exhibit strong PAH emission, whose band strength is comparable to that of late-type star-forming disk galaxies. The tidal tail, which extends 2 arcmin (~75 kpc) into the intergalactic medium, is populated by super massive star clusters likely triggered by the galaxy-galaxy interaction that has distorted UGC 10214 into its current "tadpole" shape.

We present results of a morphological analysis of a small subset of the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) galaxy population. The analysis is based on public ACS data taken inside the SWIRE N1 field, which are the deepest optical high-resolution imaging available within the SWIRE fields as of today. Our reference sample includes 156 galaxies detected by both ACS and SWIRE. Among the various galaxy morphologies, we disentangle two main classes, spheroids (or bulge-dominated galaxies) and disk-dominated ones, for which we compute the number counts as a function of flux. We then limit our sample to objects with IRAC fluxes brighter than 10 microJy, estimated ~90% completeness limit of the SWIRE catalogues, and compare the observed counts to model predictions. We find that the observed counts of the spheroidal population agree with the expectations of a hierarchical model while a monolithic scenario predicts steeper counts. Both scenaria, however, under-predict the number of late-type galaxies. These observations show that the large majority (close to 80 per cent) of the 3.6 and 4.5 micron galaxy population, even at these moderately faint fluxes, is dominated by spiral and irregular galaxies or mergers.

We describe TRILEGAL, a new populations synthesis code for simulating the stellar photometry of any Galaxy field. The code attempts to improve upon several technical aspects of star count models, by: dealing with very complete input libraries of evolutionary tracks; using a stellar spectral library to simulate the photometry in any broad-band system; being very versatile allowing easy changes in the input libraries and in the description of all of its ingredients -- like the SFR, AMR, IMF, and geometry of Galaxy components. In a previous paper (Groenewegen et al. 2002), the code was first applied to describe the very deep star counts of the CDFS stellar catalogue. Here, we briefly describe its initial calibration using EIS-deep and DMS star counts, which are adequate samples to probe both the halo and the disc components of largest scale heights (oldest ages). We then present the changes in the calibration that were necessary to cope with some improvements in the model input data, and the use of more extensive photometry datasets: the relatively shallower 2MASS catalogue, which probes mostly the disc at intermediate ages, and the immediate solar neighbourhood as sampled by Hipparcos, which contains a somewhat larger fraction of younger stars than deeper surveys. Remarkably, the same model calibration can reproduce well the star counts in all the above-mentioned data sets, that span from the very deep magnitudes of CDFS (16<R<23) to the very shallow ones of Hipparcos (V<8). Significant deviations are found just for fields close to the Galactic Center and Plane, and for a single set of South Galactic Pole data. The TRILEGAL code is ready to use for the variety of wide-angle surveys in the optical/infrared that will become available in the coming years.

We exploit deep combined observations with Spitzer and Chandra of the SWIRE survey in the ELAIS-N1 region, to investigate the nature of the faint X-ray and IR sources in common, to identify AGN/starburst diagnostics, and to study the sources of the X-ray and IR cosmic backgrounds. In the 17'x17' area of the Chandra ACIS-I image there are 3400 SWIRE near-IR sources with 4 sigma detections in at least 2 IRAC bands and 988 sources detected at 24micron with MIPS brighter than 0.1 mJy. Of these, 102 IRAC and 59 MIPS sources have Chandra counterparts, out of a total of 122 X-ray sources present in the area with S(0.5-8 kev)>10^(-15) erg/cm^2/s. We have constructed SEDs for each source using data from the 4 IRAC wavebands, Chandra fluxes, and optical follow-up data in the wavebands U, g', r', i', Z, and H. We fit a number of spectral templates to the SEDs at optical and infrared wavelengths to determine photometric redshifts and spectral categories, and also make use of diagnostics based on the X-ray luminosities, hardness ratios, X-ray to infrared spectral slopes and optical morphologies. Although we have spectroscopic redshifts for only a minority of the Chandra sources, the available SEDs constrain the redshifts for most of the sample sources, which turn out to be typically at 0.5<z<2. We find that 39% of the Chandra sources are dominated by type-1 AGN emission, 23% display optical/IR spectra typical of type-2 AGNs, while the remaining 38% fraction show starburst-like or even normal galaxy spectra. Since we prove that all these galaxies are dominated by AGN emission in X-rays this brings the fraction of type-1 AGNs to be 80% of the type-2: even assuming that all the Chandra sources undetected by Spitzer are type-2 AGNs, the type-1 fraction would exceed 1/3 of the total population (abridged).

(Abridged) In this paper we present photometric redshift estimates for a sample of X-ray selected sources detected in the wide field (~2 deg^2), bright [f_{X} (0.5-8 keV)~10^{-14} cgs] XMM-Newton/2dF survey. Unlike deeper X-ray samples comprising a large fraction of sources with colours dominated by the host galaxy, our bright survey primarily probes the QSO X-ray population. Therefore photometric redshift methods employing both galaxy and QSO templates need to be used. We employ the photometric redshift technique of Hatziminaoglou, Mathez & Pello (2000) using 5-band photometry from the SDSS. We separate our X-ray sources according to their optical profile to point-like and extended. We apply QSO and galaxy templates to the point-like and extended sources respectively. X-ray sources associated with Galactic stars are identified and discarded from our point-like sample on the basis of their low X-ray--to--optical flux ratio and their broad band colours that are best fit by stellar templates. Comparison of our results with spectroscopic redshifts available, allows calibration of our method and estimation of the photometric redshift accuracy. For ~70 per cent of the point-like sources photometric redshifts are correct within dz <= 0.3 (or ~75 per cent have dz/(1+z) <= 0.2), and the rms scatter is estimated to be sigma_z = 0.30. For the optically extended objects the photometric redshifts work only in the case of red (g - r > 0.5 mag) sources yielding dz <= 0.15 and dz/(1+z) <= 0.2 for 73 and 93 per cent respectively. However, we find that the above photometric redshift technique does not work in the case of extended sources with blue colours (g - r < 0.5): such sources cannot be fit successfully by QSO or galaxy templates, or any linear combination of the two.

We present a mid-infrared analysis of 35 quasars with spectroscopic redshifts selected from the Spitzer Wide-area InfraRed Extragalactic Survey (SWIRE). We discuss their optical and mid-infrared (MIR) colors, and show that these quasars occupy well defined regions in MIR color-color space. We examine the issue of type-I AGN candidate selection in detail and propose new selection methods based on mid-IR colors. The available multi-band data allows us to construct two new, well-sampled quasar templates, covering wavelengths from the ultraviolet to the MIR.

Most counterparts of the identified low-latitude gamma-ray sources are Isolated Neutron Stars (INSs). They are characterized by an extremely high value of the X/optical flux ratio, therefore a systematic X-ray/optical coverage of the fields of unidentified low-latitude gamma-ray sources is the best way to unveil possible INS counterparts. Since the low-latitude sources are heavily affected by the interstellar absorption in both the X-ray and optical bands, we decided to apply the above strategy to two middle-latitude EGRET sources, which could belong to a local galactic population: 3EG J0616-3310 and 3EG J1249-8330. Here we report on the global X-ray characterisation of about 300 objects, on their candidate optical counterparts and on the preliminary results of their identification.

This paper discusses properties of the European Large Area ISO Survey 15 micron quasars and tries to establish a robust method of quasar selection for future use within the Spitzer Wide-Area Infrared Extragalactic Survey (SWIRE) framework. The importance of good quality ground-based optical data is stressed, both for the candidates selection and for the photometric redshifts estimates. Colour-colour plots and template fitting are used for these purposes. The properties of the 15 micron quasars sample are studied, including variability and black hole masses and compared to the properties of other quasars that lie within the same fields but have no mid-infrared counterparts. The two subsamples do not present substantial differences and are believed to come from the same parent population.

The limited angular resolution of gamma-ray telescopes prevents the straight identification of the majority of the sources detected so far. This is particularly true for the low latitude, probably galactic ones, only 10 % of which has been identified. Most of the counterparts are Isolated Neutron Stars, both radio-loud and radio-quiet, characterised by an extremely high value of X-ray-to-optical flux ratio. The best way to search for INSs in the error boxes of unidentified EGRET sources is to perform the X-ray coverage of the gamma-ray field, followed by the optical characterization of each X-ray source. We applied this procedure to two EGRET sources, which could belong to a local galactic population: 3EG J0616-3310 and 3EG J1249-8330. Here we report on the analysis of about 300 X-ray objects, as well as on their optical study.

The technique of estimating redshifts using photometric rather than spectroscopic observations has recently received great attention due to its simplicity and the accuracy of the results obtained. In this work, we estimate photometric redshifts for an X-ray selected QSO sample. This is the first time this technique is applied on such a sample. We first calculate the accuracy of the results obtained by comparing photometric to spectroscopic redshifts for a sub-sample of our QSO sample: for the majority (~67%) of the objects in this sub-sample, photometric redshift estimates are correct within Dz<0.3. We then derive the photometric redshift distribution for the whole QSO sample. In the future, we expect to use the photometric redshift distribution in order to derive the distributions of properties such as the Hardness Ratio and hence the hydrogen column density, the luminosity function etc. As an example, we estimate here the dependence of the Hardness Ratio of the QSO sample on photometric redshift.

We construct a physically motivated analytical model for the quasar luminosity function, based on the joint star formation and feeding of massive black holes suggested by the observed correlation between the black hole mass and the stellar mass of the hosting spheroids. The parallel growth of massive black holes and host galaxies is assumed to be triggered by major mergers of haloes. The halo major merger rate is computed in the frame of the extended Press-Schechter model. The evolution of black holes on cosmological timescales is achieved by the integration of the governing set of differential equations, established from a few reasonable assumptions that account for the distinct (Eddington-limited or supply-limited) accretion regimes. Finally, the typical lightcurves of the reactivated quasars are obtained under the assumption that, in such accretion episodes, the fall of matter onto the black hole is achieved in a self-regulated stationary way. The predicted quasar luminosity function is compared to the luminosity functions of the 2dF QSO sample and other, higher redshift data. We find good agreement in all cases, except for z<1 where the basic assumption of our model is likely to break down.

(abridged) Stellar catalogues in five passbands (UBVRI) over an area of approximately 0.3 deg^2, comprising about 1200 objects, and in seven passbands (UBVRIJK) over approximately 0.1 deg^2, comprising about 400 objects, in the direction of the Chandra Deep Field South are presented. The 90% completeness level of the number counts is reached at approximately U = 23.8, B = 24.0, V = 23.5, R = 23.0, I = 21.0, J = 20.5, K = 19.0. A scheme is presented to select point sources from these catalogues, by combining the SExtractor parameter CLASS_STAR from all available passbands. Probable QSOs and unresolved galaxies are identified by using the previously developed \chi^2-technique (Hatziminaoglou et al 2002), that fits the overall spectral energy distributions to template spectra and determines the best fitting template. The observed number counts, colour-magnitude diagrams, colour-colour diagrams and colour distributions are presented and, to judge the quality of the data, compared to simulations based on the predictions of a Galactic Model convolved with the estimated completeness functions and the error model used to describe the photometric errors of the data. The resulting stellar catalogues and the objects identified as likely QSOs and unresolved galaxies with coordinates, observed magnitudes with errors and assigned spectral types by the $\chi^2$-technique are presented and are publicly available.

The limited angular resolution of gamma-ray telescopes prevents the straight identification of the majority of the sources detected so far. While the only certified galactic gamma-ray emitters are associated with pulsars, more than 90% of the nearly 100 low latitude sources detected by EGRET on the GRO lack an identification. The best identification strategy devised over the years relies on the use of X-ray and optical data to single out possible counterparts, rating them on the basis of their extreme F_x/F_opt. Here we describe our multiwavelength programme based on the EPIC mapping of selected EGRET sources complemented by the optical coverage of the same fields obtained through the Wide Field Imager (WFI) of the ESO 2.2 m telescope. The field of view of the WFI, comparable to that of EPIC, will allow to directly cross-correlate X and optical data speeding up significantly the selection of interesting candidates worth follow-up optical (and X-ray) studies.

This paper describes the methodology currently being implemented in the EIS pipeline for analysing optical/infrared multi-colour data. The aim is to identify different classes of objects as well as possible undesirable features associated with the construction of colour catalogues. The classification method used is based on the xi^2-fitting of template spectra to the observed SEDs, as measured through broad-band filters. Its main advantage is the simultaneous use of all colours, properly weighted by the photometric errors. In addition, it provides basic information on the properties of the classified objects (eg redshift, effective temperature). These characteristics make the xi^2-technique ideal for handling large multi-band datasets. The results are compared to the more traditional colour-colour selection and, whenever possible, to model predictions. In order to identify objects with odd colours, either associated with rare populations or to possible problems in the catalogue, outliers are searched for in the multi-dimensional colour space using a nearest-neighbour criterion. Outliers with large xi^2-values are individually inspected to further investigate their nature. The tools developed are used for a preliminary analysis of the multi-colour point source catalogue constructed from the optical/infrared imaging data obtained for the Chandra Deep Field South (CDF-S). These data are publicly available, representing the first installment of the ongoing EIS Deep Public Survey.

This paper presents multi-passband optical data obtained from observations of the Chandra Deep Field South (CDF-S), located at alpha ~ 3h 32m, delta ~ -27d 48m. The observations were conducted at the ESO/MPG 2.2m telescope at La Silla using the 8kx8k Wide-Field Imager (WFI). This data set, taken over a period of one year, represents the first field to be completed by the ongoing Deep Public Survey (DPS) being carried out by the ESO Imaging Survey (EIS) project. This paper describes the optical observations, the techniques employed for un-supervised pipeline processing and the general characteristics of the final data set. The paper includes data taken in six different filters U'UBVRI. The data cover an area of about 0.25 square degrees reaching 5 sigma limiting magnitudes of U'_AB=26.0, U_AB=25.7, B_AB=26.4$, V_AB=25.4, R_AB=25.5 and I_AB= 24.7 mag, as measured within a 2xFWHM aperture. The optical data covers the area of ~ 0.1

This paper presents the first set of fully calibrated images and associated stellar catalogs of the Pre-FLAMES survey being carried out by the ESO Imaging Survey (EIS) project. The primary goal of this survey is to provide the ESO community with data sets from which suitable target lists can be extracted for follow-up observations with the new VLT facility FLAMES. For this purpose 160 stellar fields have been selected for observations in B, V and I using the 8kx8k Wide Field Imager (WFI) at the MPG/ESO 2.2 m telescope at La Silla. These data have been used to assess the observing strategy adopted and to define suitable reduction techniques and procedures for the preparation of input catalogs for FLAMES. The images and catalogs presented here are publicly available and can be requested from the URL address ``http://www.eso.org/eis''.

This paper presents new J and Ks near-infrared data obtained from observations of the Chandra Deep Field South (CDF-S) conducted at the ESO 3.5m New Technology Telescope (NTT). These data were taken as part of the ongoing Deep Public Survey (DPS) being carried out by the ESO Imaging Survey (EIS) program, extending the EIS-DEEP survey. Combined these surveys now provide a contiguous coverage over an area of 400 square arcmin in the near-infrared, nearly matching that covered by the deep X-ray observations of Chandra, four times the area of the original EIS-DEEP survey. The paper briefly describes the observations and the new techniques being employed for pipeline processing jittered infrared observations, which include unbiased de-fringing and sky-background subtraction, pixel-based astrometry and stacking and pixel registration based on a multi-resolution decomposition of the images. The astrometric solution is based on a pre-release of the GSC-II catalog and has an accuracy of < 0.15 arcsec. The final images for 12 pointings presented here reach median 5 sigma limiting magnitudes of J_AB~23.4 and K_AB~22.6 as measured within an aperture 2xFWHM. The frame to frame variation of the photometric zero-point is estimated to be <0.09 mag. The data are publicly available in the form of fully calibrated J and Ks pixel maps and source lists extracted for each pointing. These data can be requested through the URL ``http://www.eso.org/eis''.

The growth of Super Massive Black Holes and the parallel development of activity in galactic nuclei are implemented in an analytic code of hierarchical clustering. The evolution of the luminosity function of quasars and AGN will be computed with special attention paid to the connection between quasars and Seyfert galaxies. One of the major interests of the model is the parallel study of quasar formation and evolution and the History of Star Formation.

We have developed a consistent analytical model to describe the observed evolution of the quasar luminosity function. Our model combines black hole mass distributions based on the Press - Schechter theory of the structure formation in the Universe with quasar luminosity functions resulting from a physics-based emission model that takes into account the time-dependent phenomena occurring in the accretion disks. Quasar evolution and CDM models are mutually constraining, therefore our model gives an estimation of the exponent, n, of the power spectrum, P(k), which is found to be -1.8 < n < -1.6. We were able to reject a generally assumed hypothesis of a constant ratio between Dark Matter Halo and the Black Hole mass, since the observed data could not be fitted under this assumption. We found that the relation between the Dark Matter Halos and Black Hole masses is better described by M_{BH}=M_{DMH}^{0.668}. This model provides a reasonable fit to the observed quasar luminosity function at redshifts higher than ~2.0. We suggest that the disagreement at lower redshift is due to mergers. Based on the agreement at high redshift, we estimated the merger rate at lower redshift, and argue that this rate should depend on the redshift, like (1+z)^3.

We present a quasar candidate identification technique based on multicolor photometry. The traditional multi-dimensional method (2 $\times$ N dimensions, where N is the number of the color-color diagrams) is reduced to a one-dimensional technique, which consists in a standard fitting procedure, where the observed spectral energy distributions are compared to quasar simulated spectra and stellar templates. This new multicolor approach is firstly applied to simulated catalogues and its {\it efficiency} is examined in various redshift ranges, as a function of the filter combination and the available observing time for spectroscopy. We conclude that this method is better suited than the usual multicolor selection techniques to quasar identification, especially for high-redshift quasars. The application of the method to real quasar samples found in the literature results in an {\it efficiency} comparable to the one obtained from the use of color-color diagrams. The major advantage of the new method is the estimation of the {\it photometric} redshift of quasar candidates, enabling, in almost all cases, spectroscopy to be targeted to best suited wavelength ranges.

Large, homogeneous quasar samples are necessary tools for the study of QSO statistics, cosmological tests, large scale structure and AGN evolution. These samples must be complete within well defined flux limits at all redshifts. An observational strategy without previous photometric selection of quasar candidates is described, based essentially on the VIRMOS Survey spectroscopy.

The evolution of optically selected quasars is usually supposed to be well described by a single constant parameter, either k_L or k_D, depending whether we refer to luminosity or density evolution. In this paper we present a study of the variations of the evolution parameters with redshift, for different cosmological models, in order to probe the differential evolution with redshift. Two different quasar samples have been analyzed, the AAT Boyle's et al. and the LBQS catalogues. Basically, these samples are divided in redshift intervals and in each of them k_L and k_D are estimated by forcing that <V/Vmax>=0.5. The dependence with respect to the cosmological parameters is small. Both AAT and LBQS show roughly the same tendencies. LBQS, however, shows strong fluctuations, whose origin is not statistical but rather due to the selection criteria. A discussion on selection techniques, biases and binning effects explains the differences between these results. We finally conclude that the evolution parameter is almost constant in the redshift range 0.7 < z < 1.7, at least within 2 sigma, while it decreases slightly afterwards. Results depend on the binning chosen (but not in a very significant way). The method has been tested with Monte-Carlo simulated catalogues in order to give a better understanding of the results coming from the real catalogues. A correlation between k_L (k_D) and <V/Vmax> is also derived and is used for the calculation of the error bars on the evolution parameter.