Using the EMIR instrument on the IRAM 30m telescope, we conducted a spectroscopic redshift search of seven z$_{\rm phot}$ $\sim$ 4 sub-millimetre bright galaxies selected from the Herschel Bright Sources (HerBS) sample with fluxes at 500 $\mu$m greater than 80 mJy. For four sources, we obtained spectroscopic redshifts between 3.4 < z < 4.1 through the detection of multiple CO-spectral lines with J $\leq$ 3. Later, we detected low-J transitions for two of these sources with the GBT including the CO(1-0) transition. For the remaining three sources, more data are needed to determine the spectroscopic redshift unambiguously. The measured CO luminosities and line widths suggest that all these sources are gravitationally lensed. These observations demonstrate that the 2 mm window is indispensable to confirm robust spectroscopic redshifts for z < 4 sources. Finally, we present an efficient graphical method to correctly identify spectroscopic redshifts.

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08<z<4.05$ with a median value of $z=2.9\pm$0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ~140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700$\pm$300 km/s and a median of 800 km/s. We analyse the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous $L_{\rm FIR} > 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies.

GravityCam is a new concept of ground-based imaging instrument capable of delivering significantly sharper images from the ground than is normally possible without adaptive optics. Advances in optical and near infrared imaging technologies allow images to be acquired at high speed without significant noise penalty. Aligning these images before they are combined can yield a 3-5 fold improvement in image resolution. By using arrays of such detectors, survey fields may be as wide as the telescope optics allows. We describe the instrument and detail its application to accelerate greatly the rate of detection of Earth size planets by gravitational microlensing. GravityCam will improve substantially the quality of weak shear studies of dark matter distribution in distant clusters of galaxies. An extensive microlensing survey will also provide a vast dataset for asteroseismology studies, and GravityCam promises to generate a unique data set on the population of the Kuiper belt and possibly the Oort cloud.

How far can we use multi-wavelength cross-identifications to deconvolve far-infrared images? In this short research note I explore a test case of CLEAN deconvolutions of simulated confused 850 micron SCUBA-2 data, and explore the possible scientific applications of combining this data with ostensibly deeper TolTEC Large Scale Structure (LSS) survey 1.1mm-2mm data. I show that the SCUBA-2 can be reconstructed to the 1.1mm LMT resolution and achieve an 850 micron deconvolved sensitivity of 0.7 mJy RMS, an improvement of at least ~1:5x over naive point source filtered images. The TolTEC/SCUBA-2 combination can constrain cold (<10K) observed-frame colour temperatures, where TolTEC alone cannot.

(Abridged summary) Responding to ESA's Voyage 2050 call to define the long-term plan for the future space missions that will address the astrophysics science questions during the 2035-2050 cycle, we propose a dedicated, ultra-deep spectroscopic survey in the near infrared (NIR), that will target a mass-limited sample of galaxies during two of the most fundamental epochs of cosmic evolution: the formation of the first galaxies (at z>6; cosmic dawn), and at the peak of galaxy formation activity (between redshift z=1 and 3; cosmic noon). By way of NIR observations, it is possible to study the Lyman-alpha region in the former, and the optical rest-frame in the latter, allowing us to extract fundamental observables such as gas and stellar kinematics, chemical abundances, and ages, providing a unique legacy database covering these two crucial stages of cosmic evolution. A dedicated, space-based facility will overcome the challenges faced by ground-based telescopes, no matter how large the aperture, or the reduced field of view and low multiplex factor of the best space-based instrument in the near future, namely NIRSpec at the JWST. Our project (codename Chronos) aims to produce about 1 million high quality spectra, with a high S/N in the continuum, where information about the underlying stellar populations is encoded.

The largest Herschel extragalactic surveys, H-ATLAS and HerMES, have selected a sample of "ultrared" dusty, star-forming galaxies (DSFGs) with rising SPIRE flux densities ($S_{500} > S_{350} > S_{250}$; so-called "500 $\mu$m-risers") as an efficient way for identifying DSFGs at higher redshift ($z > 4$). In this paper, we present a large Spitzer follow-up program of 300 Herschel ultrared DSFGs. We have obtained high-resolution ALMA, NOEMA, and SMA data for 63 of them, which allow us to securely identify the Spitzer/IRAC counterparts and classify them as gravitationally lensed or unlensed. Within the 63 ultrared sources with high-resolution data, $\sim$65% appear to be unlensed, and $\sim$27% are resolved into multiple components. We focus on analyzing the unlensed sample by directly performing multi-wavelength spectral energy distribution (SED) modeling to derive their physical properties and compare with the more numerous $z \sim 2$ DSFG population. The ultrared sample has a median redshift of 3.3, stellar mass of 3.7 $\times$ 10$^{11}$ $M_{\odot}$, star formation rate (SFR) of 730 $M_{\odot}$yr$^{-1}$, total dust luminosity of 9.0 $\times$ 10$^{12}$ $L_{\odot}$, dust mass of 2.8 $\times$ 10$^9$ $M_{\odot}$, and V-band extinction of 4.0, which are all higher than those of the ALESS DSFGs. Based on the space density, SFR density, and stellar mass density estimates, we conclude that our ultrared sample cannot account for the majority of the star-forming progenitors of the massive, quiescent galaxies found in infrared surveys. Our sample contains the rarer, intrinsically most dusty, luminous and massive galaxies in the early universe that will help us understand the physical drivers of extreme star formation.

The Euclid telescope, due for launch in 2021, will perform an imaging and slitless spectroscopy survey over half the sky, to map baryon wiggles and weak lensing. During the survey Euclid is expected to resolve 100,000 strong gravitational lens systems. This is ideal to find rare lens configurations, provided they can be identified reliably and on a reasonable timescale. For this reason we have developed a Convolutional Neural Network (CNN) that can be used to identify images containing lensing systems. CNNs have already been used for image and digit classification as well as being used in astronomy for star-galaxy classification. Here our CNN is trained and tested on Euclid-like and KiDS-like simulations from the Euclid Strong Lensing Group, successfully classifying 77% of lenses, with an area under the ROC curve of up to 0.96. Our CNN also attempts to classify the lenses in COSMOS HST F814W-band images. After convolution to the Euclid resolution, we find we can recover most systems that are identifiable by eye. The Python code is available on Github.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

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.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

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

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

GravityCam is a new concept of ground-based imaging instrument capable of delivering significantly sharper images from the ground than is normally possible without adaptive optics. Advances in optical and near infrared imaging technologies allow images to be acquired at high speed without significant noise penalty. Aligning these images before they are combined can yield a 2.5 to 3 fold improvement in image resolution. By using arrays of such detectors, survey fields may be as wide as the telescope optics allows. Consequently, GravityCam enables both wide-field high-resolution imaging and high-speed photometry. We describe the instrument and detail its application to provide demographics of planets and satellites down to Lunar mass (or even below) across the Milky Way. GravityCam is also suited to improve the quality of weak shear studies of dark matter distribution in distant clusters of galaxies and multiwavelength follow-ups of background sources that are strongly lensed by galaxy clusters. The photometric data arising from an extensive microlensing survey will also be useful for asteroseismology studies, while GravityCam can be used to monitor fast multiwavelength flaring in accreting compact objects, and promises to generate a unique data set on the population of the Kuiper belt and possibly the Oort cloud.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

A detailed analysis of Herschel-PACS observations at the North Ecliptic Pole is presented. High quality maps, covering an area of 0.44 square degrees, are produced and then used to derive potential candidate source lists. A rigorous quality control pipeline has been used to create final legacy catalogues in the PACS Green 100 micron and Red 160 micron bands, containing 1384 and 630 sources respectively. These catalogues reach to more than twice the depth of the current archival Herschel/PACS Point Source Catalogue, detecting 400 and 270 more sources in the short and long wavelength bands respectively. Galaxy source counts are constructed that extend down to flux densities of 6mJy and 19mJy (50% completeness) in the Green 100 micron and Red 160 micron bands respectively. These source counts are consistent with previously published PACS number counts in other fields across the sky. The source counts are then compared with a galaxy evolution model identifying a population of luminous infrared galaxies as responsible for the bulk of the galaxy evolution over the flux range (5-100mJy) spanned by the observed counts, contributing approximate fractions of 50% and 60% to the cosmic infrared background (CIRB) at 100 microns and 160 microns respectively.

We present structural parameters and morphological properties of faint 450-um selected submillimeter galaxies (SMGs) from the JCMT Large Program, STUDIES, in the COSMOS-CANDELS region. Their properties are compared to an 850um selected and a matched star-forming samples. We investigate stellar structures of 169 faint 450-um sources (S450=2.8-29.6mJy; S/N>4) at z<3 using HST near-infrared observations. Based on our spectral energy distribution fitting, half of such faint SMGs (LIR=10^11.65+-0.98Lsun) lie above the star-formation rate (SFR)/stellar mass plane. The size-mass relation shows that these SMGs are generally similar to less-luminous star-forming galaxies selected by NUV-r vs. r-J colors. Because of the intrinsic luminosity of the sample, their rest-frame optical emission is less extended than the 850um sources (S850>2mJy), and more extended than the star-forming galaxies in the same redshift range. For the stellar mass and SFR matched sample at z~=1 and z~=2, the size differences are marginal between faint SMGs and the matched galaxies. Moreover, faint SMGs have similar Sersic indices and projected axis ratios as star-forming galaxies with the same stellar mass and SFR. Both SMGs and the matched galaxies show high fractions (~70%) of disturbed features at z~=2, and the fractions depend on the SFRs. These suggest that their star formation activity is related to galaxy merging, and the stellar structures of SMGs are similar to those of star-forming galaxies. We show that the depths of submillimeter surveys are approaching the lower luminosity end of star-forming galaxies, allowing us to detect galaxies on the main sequence.

Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of their star formation rate (SFR) and specific SFR up to z~=0.5, by paying particular attention to their PAH properties. We conducted a low-resolution (R~=50) slitless spectroscopic survey at 5-13um of 9um flux-selected sources (>0.3 mJy) around the North Ecliptic Pole with the Infrared Camera (IRC) onboard AKARI. After removing 11 AGN candidates by using the IRC photometry, we identified 48 PAH galaxies with PAH 6.2, 7.7, and 8.6um features at z<0.5. The rest-frame optical-MIR spectral energy distributions (SEDs) based on CFHT and AKARI/IRC imaging covering 0.37-18um were produced, and analysed in conjunction with the PAH spectroscopy. We defined the PAH enhancement by using the luminosity ratio of the 7.7um PAH feature over the 3.5um stellar component of the SEDs. The rest-frame SEDs of all PAH galaxies have a universal shape with stellar and 7.7um bumps, except that the PAH enhancement significantly varies as a function of the PAH luminosities. We identified a PAH-enhanced population at z~>0.35, whose SEDs and luminosities are typical of luminous infrared galaxies. They show particularly larger PAH enhancement at high luminosity, implying that they are vigorous star-forming galaxies with elevated specific SFR. Our composite starburst model that combines a very young and optically very thick starburst with a very old population can successfully reproduce most of their SED characteristics, although we could not confirm this optically thick component from our spectral analysis.

Since their discovery, submillimetre-selected galaxies (SMGs) have revolutionized the field of galaxy formation and evolution. From the hundreds of square degrees mapped at submillimetre wavelengths, only a handful of sources have been confirmed to lie at z>5 and only two at z>6. All of these SMGs are rare examples of extreme starburst galaxies with star formation rates (SFRs) of >1000 M_sun/yr and therefore are not representative of the general population of dusty star-forming galaxies. Consequently, our understanding of the nature of these sources, at the earliest epochs, is still incomplete. Here we report the spectroscopic identification of a gravitationally amplified (mu = 9.3 +/- 1.0) dusty star-forming galaxy at z=6.027. After correcting for gravitational lensing, we derive an intrinsic less-extreme SFR of 380 +/- 50 M_sun/yr for this source and find that its gas and dust properties are similar to those measured for local Ultra Luminous Infrared Galaxies (ULIRGs), extending the local trends to a poorly explored territory in the early Universe. The star-formation efficiency of this galaxy is similar to those measured in its local analogues, despite a ~12 Gyr difference in cosmic time.

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

We have modelled high resolution ALMA imaging of six strong gravitationally lensed galaxies detected by the Herschel Space Observatory. Our modelling recovers mass properties of the lensing galaxies and, by determining magnification factors, intrinsic properties of the lensed sub-millimetre sources. We find that the lensed galaxies all have high ratios of star formation rate to dust mass, consistent with or higher than the mean ratio for high redshift sub-millimetre galaxies and low redshift ultra-luminous infra-red galaxies. Source reconstruction reveals that most galaxies exhibit disturbed morphologies. Both the cleaned image plane data and the directly observed interferometric visibilities have been modelled, enabling comparison of both approaches. In the majority of cases, the recovered lens models are consistent between methods, all six having mass density profiles that are close to isothermal. However, one system with poor signal to noise shows mildly significant differences.

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

The AKARI North Ecliptic Pole (NEP) survey consists of two survey projects: NEP-Deep (0.5 sq.deg) and NEP-Wide (5.4 sq.deg), providing with tens of thousands of galaxies. A continuous filter coverage in the mid-infrared wavelengths (7, 9, 11, 15, 18 and 24 $\mu$m) is unique to diagnose the contributions from dusty star-formation activity and AGNs. Here we present current status focused on the newly obtained optical images and near-future prospects with a new X-ray telescope. Hyper Suprime-Cam on Subaru telescope is a gigantic optical camera with huge Field of View (FoV). Thanks to the wide FoV, we successfully obtained deep optical images at g, r, i, z and Y-bands covering most of the NEP-Wide field. Using the deep optical images, we identified over 5000 optical counterparts of the mid-IR sources, presumably deeply obscured galaxies in NEP-Wide field. We also investigated properties of these infrared sources with SED-fitting. eROSITA, to be launched early 2018, is a new all-sky X-ray survey telescope, and expected to conduct ultra deep 2-10 keV imaging toward NEP. We expect unprecedentedly numerous Compton-thick AGN candidates when combined with the multi-wavelength data in NEP region.

We present a photometric catalog for Spitzer Space Telescope warm mission observations of the North Ecliptic Pole (NEP; centered at $\rm R.A.=18^h00^m00^s$, $\rm Decl.=66^d33^m38^s.552$). The observations are conducted with IRAC in 3.6 $\mu$m and 4.5 $\mu$m bands over an area of 7.04 deg$^2$ reaching 1$\sigma$ depths of 1.29 $\mu$Jy and 0.79 $\mu$Jy in the 3.6 $\mu$m and 4.5 $\mu$m bands respectively. The photometric catalog contains 380,858 sources with 3.6 $\mu$m and 4.5 $\mu$m band photometry over the full-depth NEP mosaic. Point source completeness simulations show that the catalog is 80% complete down to 19.7 AB. The accompanying catalog can be utilized in constraining the physical properties of extra-galactic objects, studying the AGN population, measuring the infrared colors of stellar objects, and studying the extra-galactic infrared background light.

We present the first detailed analysis of three extragalactic fields (IRAC Dark Field, ELAIS-N1, ADF-S) observed by the infrared satellite, AKARI, using an optimised data analysis toolkit specifically for the processing of extragalactic point sources. The InfraRed Camera (IRC) on AKARI complements the Spitzer space telescope via its comprehensive coverage between 8-24 microns filling the gap between the Spitzer IRAC and MIPS instruments. Source counts in the AKARI bands at 3.2, 4.1, 7, 11, 15 and 18 microns are presented. At near-infrared wavelengths, our source counts are consistent with counts made in other AKARI fields and in general with Spitzer/IRAC (except at 3.2 microns where our counts lie above). In the mid-infrared (11 - 18 microns) we find our counts are consistent with both previous surveys by AKARI and the Spitzer peak-up imaging survey with the InfraRed Spectrograph (IRS). Using our counts to constrain contemporary evolutionary models we find that although the models and counts are in agreement at mid-infrared wavelengths there are inconsistencies at wavelengths shortward of 7 microns, suggesting either a problem with stellar subtraction or indicating the need for refinement of the stellar population models. We have also investigated the AKARI/IRC filters, and find an AGN selection criteria out to $z<2$ on the basis of AKARI 4.1, 11, 15 and 18 microns colours.

We present images obtained with LABOCA on the APEX telescope of a sample of 22 galaxies selected via their red Herschel SPIRE 250-, 350- and $500\textrm{-}\mu\textrm{m}$ colors. We aim to see if these luminous, rare and distant galaxies are signposting dense regions in the early Universe. Our $870\textrm{-}\mu\textrm{m}$ survey covers an area of $\approx0.8\,\textrm{deg}^2$ down to an average r.m.s. of $3.9\,\textrm{mJy beam}^{-1}$, with our five deepest maps going $\approx2\times$ deeper still. We catalog 86 DSFGs around our 'signposts', detected above a significance of $3.5\sigma$. This implies a $100\pm30\%$ over-density of $S_{870}>8.5\,\textrm{mJy}$ DSFGs, excluding our signposts, when comparing our number counts to those in 'blank fields'. Thus, we are $99.93\%$ confident that our signposts are pinpointing over-dense regions in the Universe, and $\approx95\%$ confident that these regions are over-dense by a factor of at least $\ge1.5\times$. Using template SEDs and SPIRE/LABOCA photometry we derive a median photometric redshift of $z=3.2\pm0.2$ for our signposts, with an interquartile range of $z=2.8\textrm{-}3.6$. We constrain the DSFGs likely responsible for this over-density to within $|\Delta z|\le0.65$ of their respective signposts. These 'associated' DSFGs are radially distributed within $1.6\pm0.5\,\textrm{Mpc}$ of their signposts, have median SFRs of $\approx(1.0\pm0.2)\times10^3\,M_{\odot}\,\textrm{yr}^{-1}$ (for a Salpeter stellar IMF) and median gas reservoirs of $\sim1.7\times10^{11}\,M_{\odot}$. These candidate proto-clusters have average total SFRs of at least $\approx (2.3\pm0.5)\times10^3\,M_{\odot}\,\textrm{yr}^{-1}$ and space densities of $\sim9\times10^{-7}\,\textrm{Mpc}^{-3}$, consistent with the idea that their constituents may evolve to become massive ETGs in the centers of the rich galaxy clusters we see today.

The recent updates of the North Ecliptic Pole deep (0.5~deg$^2$, NEP-Deep) multi-wavelength survey covering from X-ray to radio-wave is presented. The NEP-Deep provides us with several thousands of 15~$\mu$m or 18~$\mu$m selected sample of galaxies, which is the largest sample ever made at this wavelengths. A continuous filter coverage in the mid-infrared wavelength (7, 9, 11, 15, 18, and 24~$\mu$m) is unique and vital to diagnose the contributions from starbursts and AGNs in the galaxies out to $z$=2.The new goal of the project is to resolve the nature of the cosmic star formation history at the violent epoch (e.g. $z$=1--2), and to find a clue to understand its decline from $z$=1 to present universe by utilizing the unique power of the multi-wavelength survey. The progress in this context is briefly mentioned.

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.

We built a 8um selected sample of galaxies in the NEP-AKARI field by defining 4 redshift bins with the four AKARI bands at 11, 15, 18 and 24 microns (0.15<z<0.49, 0.75<z<1.34, 1.34<z<1.7 and 1.7<z<2.05) . Our sample contains 4079 sources, 599 are securely detected with Herschel/PACS. Also adding ultraviolet (UV) data from GALEX, we fit the spectral energy distributions using the physically motivated code CIGALE to extract the star formation rate, stellar mass, dust attenuation and the AGN contribution to the total infrared luminosity (L_{IR}). We discuss the impact of the adopted attenuation curve and that of the wavelength coverage to estimate these physical parameters. We focus on galaxies with a luminosity close the characteristic L_{IR}^* in the different redshift bins to study the evolution with redshift of the dust attenuation in these galaxies.

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

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

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

The large infrastructure projects for the next decade will allow a new quantum leap in terms of new possible science. ESFRI, the European Strategy Forum on Research Infrastructures, a strategic initiative to develop the scientific integration of Europe, has identified four facilities (SKA, CTA, KM3Net and E-ELT) deserving priority in support. The ASTERICS project aims to address the cross-cutting synergies and common challenges shared by the various Astronomy ESFRI and other world-class facilities. The project (22 partners across Europe) is funded by the EU Horizon 2020 programme with 15 MEuro in 4 years. It brings together for the first time the astronomy, astrophysics and particle astrophysics communities, in addition to other related research infrastructures.

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

We present the first results of a new data analysis pipeline for processing extragalactic AKARI/IRC images. The main improvements of the pipeline over the standard analysis are the removal of Earth shine and image distortion correction. We present the differential number counts of the AKARI/IRC S11 filter IRAC validation field. The differential number counts are consistent with S11 AKARI NEP deep and 12 microns WISE NEP number counts, and with a phenomenological backward evolution galaxy model, at brighter fluxes densities. There is a detection of deeper galaxies in the IRAC validation field.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

We investigate the region around the Planck-detected z=3.26 gravitationally lensed galaxy HATLAS J114637.9-001132 (hereinafter HATLAS12-00) using both archival Herschel data from the H-ATLAS survey and using submm data obtained with both LABOCA and SCUBA2. The lensed source is found to be surrounded by a strong overdensity of both Herschel-SPIRE sources and submm sources. We detect 17 bright (S_870 >~7 mJy) sources at >4sigma closer than 5 arcmin to the lensed object at 850/870 microns. Ten of these sources have good cross-identifications with objects detected by Herschel-SPIRE which have redder colours than other sources in the field, with 350 micron flux > 250 micron flux, suggesting that they lie at high redshift. Submillimeter Array (SMA) observations localise one of these companions to ~1 arcsecond, allowing unambiguous cross identification with a 3.6 and 4.5 micron Spitzer source. The optical/near-IR spectral energy distribution (SED) of this source is measured by further observations and found to be consistent with z>2, but incompatible with lower redshifts. We conclude that this system may be a galaxy cluster/protocluster or larger scale structure that contains a number of galaxies undergoing starbursts at the same time.

Gravitational lensing has seen a surge of interest in the past few years. The handful of strong lensing systems known in the year 2000 has now been replaced with hundreds, thanks to innovative multi-wavelength selection, and there is an imminent prospect of thousands of lenses from Herschel and other sub-millimetre surveys. Euclid and the Square Kilometre Array promise tens or even hundreds of thousands. Gravitational lensing is one of the very few probes capable of mapping dark matter halo distributions. Lensing also provides independent cosmological parameter estimates and enables the study of galaxy populations that are otherwise too faint for detailed study. SALT is extremely well placed to have an enormous impact with follow-up observations of foreground lenses and background sources from e.g. Herschel, the South Pole Telescope, the Atacama Cosmology Telescope, Euclid and the Square Kilometre Array. This paper reviews the prospects for high-impact SALT science and the many constraints of galaxy evolution that can result.

Submillimetre and millimetre-wave surveys with Herschel and the South Pole Telescope have revolutionised the discovery of strong gravitational lenses. Their follow-ups have been greatly facilitated by the multi-wavelength supplementary data in the survey fields. The forthcoming Euclid optical/near-infrared space telescope will also detect strong gravitational lenses in large numbers, and orbital constraints are likely to require placing its deep survey at the North Ecliptic Pole (the natural deep field for a wide class of ground-based and space-based observatories including AKARI, JWST and SPICA). In this paper I review the current status of the multi-wavelength survey coverage in the NEP, and discuss the prospects for the detection of strong gravitational lenses in forthcoming or proposed facilities such as Euclid, FIRSPEX and SPICA.

We present new VLA 22-GHz and e-MERLIN 5-GHz observations of CLASS B1030+074, a two-image strong gravitational lens system whose background source is a compact flat-spectrum radio quasar. In such systems we expect a third image of the background source to form close to the centre of the lensing galaxy. The existence and brightness of such images is important for investigation of the central mass distributions of lensing galaxies, but only one secure detection has been made so far in a galaxy-scale lens system. The noise levels achieved in our new B1030+074 images reach 3 microJy/beam and represent an improvement in central image constraints of nearly an order of magnitude over previous work, with correspondingly better resulting limits on the shape of the central mass profile of the lensing galaxy. Simple models with an isothermal outer power law slope now require either the influence of a central supermassive black hole, or an inner power law slope very close to isothermal, in order to suppress the central image below our detection limit. Using the central mass profiles inferred from light distributions in Virgo galaxies, moved to z=0.5, and matching to the observed Einstein radius, we now find that 45% of such mass profiles should give observable central images, 10% should give central images with a flux density still below our limit, and the remaining systems have extreme demagnification produced by the central SMBH. Further observations of similar objects will therefore allow proper statistical constraints to be placed on the central properties of elliptical galaxies at high redshift.

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

Exploiting the sensitivity and spatial resolution of the Atacama Large Millimeter/submillimeter Array (ALMA), we have studied the morphology and the physical scale of the interstellar medium - both gas and dust - in SGP38326, an unlensed pair of interacting starbursts at $z= 4.425$. SGP38326 is the most luminous star bursting system known at $z > 4$ with an IR-derived ${\rm SFR \sim 4300 \,} M_\odot \, {\rm yr}^{-1}$. SGP38326 also contains a molecular gas reservoir among the most massive ever found in the early Universe, and it is the likely progenitor of a massive, red-and-dead elliptical galaxy at $z \sim 3$. Probing scales of $\sim 0.1"$ or $\sim 800 \, {\rm pc}$ we find that the smooth distribution of the continuum emission from cool dust grains contrasts with the more irregular morphology of the gas, as traced by the [CII] fine structure emission. The gas is also extended over larger physical scales than the dust. The velocity information provided by the resolved [CII] emission reveals that the dynamics of the two components of SGP38326 are compatible with disk-like, ordered rotation, but also reveals an ISM which is turbulent and unstable. Our observations support a scenario where at least a subset of the most distant extreme starbursts are highly dissipative mergers of gas-rich galaxies.

We present a study of the molecular gas properties in a sample of 98 HI - flux selected spiral galaxies within $\sim25$ Mpc, using the CO $J=3-2$ line observed with the James Clerk Maxwell Telescope. We use the technique of survival analysis to incorporate galaxies with CO upper limits into our results. Comparing the group and Virgo samples, we find a larger mean H$_{2}$ mass in the Virgo galaxies, despite their lower mean HI mass. This leads to a significantly higher H$_{2}$ to HI ratio for Virgo galaxies. Combining our data with complementary H$\alpha$ star formation rate measurements, Virgo galaxies have longer molecular gas depletion times compared to group galaxies, due to their higher H$_{2}$ masses and lower star formation rates. We suggest that the longer depletion times may be a result of heating processes in the cluster environment or differences in the turbulent pressure. From the full sample, we find that the molecular gas depletion time has a positive correlation with the stellar mass, indicative of differences in the star formation process between low and high mass galaxies, and a negative correlation between the molecular gas depletion time and the specific star formation rate.

We present infrared galaxy luminosity functions (LFs) in the AKARI North Ecliptic Pole (NEP) deep field using recently-obtained, wider CFHT optical/near-IR images. AKARI has obtained deep images in the mid-infrared (IR), covering 0.6 deg$^2$ of the NEP deep field. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of optical coverage of the full AKARI NEP survey. To rectify the situation, we recently obtained CFHT optical and near-IR images over the entire AKARI NEP deep field. These new CFHT images are used to derive accurate photometric redshifts, allowing us to fully exploit the whole AKARI NEP deep field. AKARI's deep, continuous filter coverage in the mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on spectral energy distribution (SED) fitting, which was the largest uncertainty in previous studies. Total infrared luminosity (TIR) is also obtained more reliably due to the superior filter coverage. The resulting restframe 8$\mu$m, 12$\mu$m, and TIR LFs at $0.15<z<2.2$ are consistent with previous works, but with reduced uncertainties, especially at the high luminosity-end, thanks to the wide field coverage. In terms of cosmic infrared luminosity density ($\Omega_{\mathrm{IR}}$), we found that the $\Omega_{\mathrm{IR}}$ evolves as $\propto (1+z)^{4.2\pm 0.4}$.

Using results from the Herschel Astrophysical Terrahertz Large-Area Survey and the Galaxy and Mass Assembly project, we show that, for galaxy masses above approximately 1.0e8 solar masses, 51% of the stellar mass-density in the local Universe is in early-type galaxies (ETGs: Sersic n > 2.5) while 89% of the rate of production of stellar mass-density is occurring in late-type galaxies (LTGs: Sersic n < 2.5). From this zero-redshift benchmark, we have used a calorimetric technique to quantify the importance of the morphological transformation of galaxies over the history of the Universe. The extragalactic background radiation contains all the energy generated by nuclear fusion in stars since the Big Bang. By resolving this background radiation into individual galaxies using the deepest far-infrared survey with the Herschel Space Observatory and a deep near-infrared/optical survey with the Hubble Space Telescope (HST), and using measurements of the Sersic index of these galaxies derived from the HST images, we estimate that approximately 83% of the stellar mass-density formed over the history of the Universe occurred in LTGs. The difference between this and the fraction of the stellar mass-density that is in LTGs today implies there must have been a major transformation of LTGs into ETGs after the formation of most of the stars.

(Abridged) We aim to study the evolution of dust attenuation in galaxies selected in the IR in the redshift range in which they are known to dominate the star formation activity in the universe. The comparison with other measurements of dust attenuation in samples selected using different criteria will give us a global picture of the attenuation at work in star-forming galaxies and its evolution with redshift. Using multiple filters of IRC instrument, we selected more than 4000 galaxies from their rest-frame emission at 8 microns, from z~0.2 to 2$. We built SEDs from the rest-frame UV to the far-IR by adding data in the optical-NIR and from GALEX and Herschel surveys. We fit SEDs with the physically-motivated code CIGALE. We test different templates for AGNs and recipes for dust attenuation and estimate stellar masses, SFRs, amount of dust attenuation, and AGN contribution to the total IR luminosity. The AGN contribution to the total IR luminosity is found to be on average approximately 10% with a slight increase with redshift. Dust attenuation in galaxies dominating the IR luminosity function is found to increase from z=0 to z=1 and to remain almost constant from z=1 to z=1.5. Conversely, when galaxies are selected at a fixed IR luminosity, their dust attenuation slightly decreases as redshift increases but with a large dispersion. The attenuation in our mid-IR selected sample is found ~ 2 mag higher than that found globally in the universe or in UV and Halpha line selections in the same redshift range. This difference is well explained by an increase of dust attenuation with the stellar mass, in global agreement with other recent studies. Starbursting galaxies do not systematically exhibit a high attenuation

We present the rest-frame optical spectral properties of 155 luminous quasars at 3.3<z<6.4 taken with the AKARI space telescope, including the first detection of H$\alpha$ emission line as far out as z~6. We extend the scaling relation between the rest-frame optical continuum and line luminosity of active galactic nuclei (AGNs) to the high luminosity, high redshift regime that has rarely been probed before. Remarkably, we find that a single log-linear relation can be applied to the 5100${\rm \AA}$ and H$\alpha$ AGN luminosities over a wide range of luminosity (10$^{42}$<$L_{5100}$<10$^{47}$ergs/s) or redshift (0<z<6), suggesting that the physical mechanism governing this relation is unchanged from z=0 to 6, over five decades in luminosity. Similar scaling relations are found between the optical and the UV continuum luminosities or line widths. Applying the scaling relations to the H$\beta$ black hole mass ($M_{\rm BH}$) estimator of local AGNs, we derive the $M_{\rm BH}$ estimators based on H$\alpha$, MgII, and CIV lines, finding that the UV-line based masses are overall consistent with the Balmer-line based, but with a large intrinsic scatter of 0.40dex for the CIV estimates. Our 43 $M_{\rm BH}$ estimates from H$\alpha$ confirm the existence of BHs as massive as ~10$^{10}M_{\odot}$ out to z~5, and provide a secure footing for previous MgII-line based studies that a rapid $M_{\rm BH}$ growth has occurred in the early universe.

Strong gravitational lenses provide an important tool to measure masses in the distant Universe, thus testing models for galaxy formation and dark matter; to investigate structure at the Epoch of Reionization; and to measure the Hubble constant and possibly w as a function of redshift. However, the limiting factor in all of these studies has been the currently small samples of known gravitational lenses (~10^2). The era of the SKA will transform our understanding of the Universe with gravitational lensing, particularly at radio wavelengths where the number of known gravitational lenses will increase to ~10^5. Here we discuss the technical requirements, expected outcomes and main scientific goals of a survey for strong gravitational lensing with the SKA. We find that an all-sky (3pi sr) survey carried out with the SKA1-MID array at an angular resolution of 0.25-0.5 arcsec and to a depth of 3 microJy / beam is required for studies of galaxy formation and cosmology with gravitational lensing. In addition, the capability to carryout VLBI with the SKA1 is required for tests of dark matter and studies of supermassive black holes at high redshift to be made using gravitational lensing.

We discuss the potential of a next generation space-borne CMB experiment for studies of extragalactic sources with reference to COrE+, a project submitted to ESA in response to the M4 call. We consider three possible options for the telescope size: 1m, 1.5m and 2m (although the last option is probably impractical, given the M4 boundary conditions). The proposed instrument will be far more sensitive than Planck and will have a diffraction-limited angular resolution. These properties imply that even the 1m telescope option will perform substantially better than Planck for studies of extragalactic sources. The source detection limits as a function of frequency have been estimated by means of realistic simulations. The most significant improvements over Planck results are presented for each option. COrE+ will provide much larger samples of truly local star-forming galaxies, making possible analyses of the properties of galaxies (luminosity functions, dust mass functions, star formation rate functions, dust temperature distributions, etc.) across the Hubble sequence. Even more interestingly, COrE+ will detect, at |b|> 30 deg, thousands of strongly gravitationally lensed galaxies. Such large samples are of extraordinary astrophysical and cosmological value in many fields. Moreover, COrE+ high frequency maps will be optimally suited to pick up proto-clusters of dusty galaxies, i.e. to investigate the evolution of large scale structure at larger redshifts than can be reached by other means. Thanks to its high sensitivity COrE+ will also yield a spectacular advance in the blind detection of extragalactic sources in polarization. This will open a new window for studies of radio source polarization and of the global properties of magnetic fields in star forming galaxies and of their relationships with SFRs.

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar.

We present rest-frame optical and near-infrared detections of one of the highest redshift submm-selected galaxies to date, HDF850.1. We do not detect the previously proposed counterpart HDF850.1K in new deep J and H-band HST WFC3 data, placing a strong limit of H-K>3.8, concluding that the K-band source is spurious. However, we detect 5.8um and 8um emission co-located with the submm in deblended images. After modelling and removing the flux contributions from another foreground galaxy, we constrain the stellar mass of HDF850.1 to be (2.5+/-1)x10^{11} Msun/mu for a lensing magnification mu=1.9+/-0.3, with a specific star formation rate of 8.5 Gyr^{-1}, faster than the 1-4 Gyr^{-1} observed for UV-selected galaxies at this epoch.

The Euclid space telescope will observe ~10^5 strong galaxy-galaxy gravitational lens events in its wide field imaging survey over around half the sky, but identifying the gravitational lenses from their observed morphologies requires solving the difficult problem of reliably separating the lensed sources from contaminant populations, such as tidal tails, as well as presenting challenges for spectroscopic follow-up redshift campaigns. Here I present alternative selection techniques for strong gravitational lenses in both Euclid and the Square Kilometer Array, exploiting the strong magnification bias present in the steep end of the Halpha luminosity function and HI mass function. Around 10^3 strong lensing events are detectable with this method in the Euclid wide survey. While only ~1% of the total haul of Euclid lenses, this sample has ~100% reliability, known source redshifts, high signal-to-noise and a magnification-based selection independent of assumptions of lens morphology. With the proposed Square Kilometer Array dark energy survey, the numbers of reliable strong gravitational lenses with source redshifts can reach 10^5.

We present the first galaxy counts at 18 microns using the Japanese AKARI satellite's survey at the North Ecliptic Pole (NEP), produced from the images from the NEP-Deep and NEP-Wide surveys covering 0.6 and 5.8 square degrees respectively. We describe a procedure using a point source filtering algorithm to remove background structure and a minimum variance method for our source extraction and photometry that delivers the optimum signal to noise for our extracted sources, confirming this by comparison with standard photometry methods. The final source counts are complete and reliable over three orders of magnitude in flux density, resulting in sensitivities (80 percent completeness) of 0.15mJy and 0.3mJy for the NEP-Deep and NEP-Wide surveys respectively, a factor of 1.3 deeper than previous catalogues constructed from this field. The differential source counts exhibit a characteristic upturn from Euclidean expectations at around a milliJansky and a corresponding evolutionary bump between 0.2-0.4 mJy consistent with previous mid-infrared surveys with ISO and Spitzer at 15 and 24 microns. We compare our results with galaxy evolution models confirming the striking divergence from the non-evolving scenario. The models and observations are in broad agreement implying that the source counts are consistent with a strongly evolving population of luminous infrared galaxies at redshifts higher than unity. Integrating our source counts down to the limit of the NEP survey at the 150 microJy level we calculate that AKARI has resolved approximately 55 percent of the 18 micron cosmic infrared background relative to the predictions of contemporary source count models.

We study the behaviour of polycyclic aromatic hydrocarbon emission in galaxies at z=0.3-1.4 using 1868 samples from the revised catalogue of AKARI North Ecliptic Pole Deep survey. The continuous filter coverage at 2-24um makes it possible to measure 8um luminosity, which is dominated by polycyclic aromatic hydrocarbon emission for galaxies at up to z=2. We compare the IR8 (= LIR/L(8)) and 8um to 4.5um luminosity ratio (L(8)/L(4.5)) with the starburstiness, Rsb, defined as excess of specific star -formation rate over that of main-sequence galaxy. All AGN candidates were excluded from our sample using an SED fitting. We found L(8)/L(4.5) increases with starburstiness at log Rsb < 0.5 and stays constant at higher starburstiness. On the other hand, IR8 is constant at log Rsb < 0, while it increases with starburstiness at log Rsb > 0. This behaviour is seen in all redshift range of our study. These results indicate that starburst galaxies have deficient polycyclic aromatic hydrocarbon emission compared with main-sequence galaxies. We also find that galaxies with extremely high L(8)/L(4.5) ratio have only moderate starburstiness. These results suggest that starburst galaxies have compact star-forming regions with intense radiation, which destroys PAHs and/or have dusty HII regions resulting in a lack of ionising photons.

Aims. We present an 8-band (u*, g', r', i', z', Y, J, Ks) optical to NIR deep photometric catalog based on the observations made with MegaCam and WIRCam at CFHT, and compute photometric redshifts, zp in the North Ecliptic Pole (NEP) region. Our catalog provides us to identify the counterparts, and zp for AKARI NIR/MIR sources. Results. The estimated 4sigma detection limits within an 1" aperture radius are 26.7, 25.9, 25.1, and 24.1 mag [AB] for g', r', i', and z'-bands and 23.4, 23.0, and 22.7 mag for Y, J, and Ks-bands, respectively. There are a total of 85797 sources in the band-merged catalog. An astrometric accuracy of this catalog determined by examining coordinate offsets with regard to 2MASS is 0.013" with a root mean square offset of 0.32". We distinguish 5441 secure stars from extended sources using u*-J vs. g'-Ks colors, combined with the SExtractor stellarity index of the images. Comparing with galaxy spectroscopic redshifts, we find a photometric redshift dispersion, sigma_(dz/(1+z)), of 0.032 and catastrophic failure rate, dz/(1+z)>0.15, of 5.8% at z<1, while a dispersion of 0.117 and a catastrophic failures rate of 16.6% at z>1. We extend estimate of the zp uncertainty over the full magnitude/redshift space with a redshift probability distribution function and find that our redshift are highly accurate with z'<22 at zp<2.5 and for fainter sources with z'<24 at z<1. From the investigation of photometric properties of AKARI infrared sources (23354 sources) using the g'z'Ks diagram, <5% of AKARI sources with optical counterparts are classified as high-z (1.4<z<2.5) star-forming galaxies. Among the high-z star-forming galaxies, AKARI MIR detected sources seem to be affected by stronger dust extinction compared with sources with non-detections in the AKARI MIR bands. The full, electronic version of our catalog with zp will be available at the CDS.

We report on deep near-infrared observations obtained with the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) of the first five confirmed gravitational lensing events discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We succeed in disentangling the background galaxy from the lens to gain separate photometry of the two components. The HST data allow us to significantly improve on previous constraints of the mass in stars of the lensed galaxy and to perform accurate lens modelling of these systems, as described in the accompanying paper by Dye et al. We fit the spectral energy distributions of the background sources from near-IR to millimetre wavelengths and use the magnification factors estimated by Dye et al. to derive the intrinsic properties of the lensed galaxies. We find these galaxies to have star-formation rates of approximately 400 to 2000 M_sol/yr, with approximately (6-25)x10^10 M_sol of their baryonic mass already turned into stars. At these rates of star formation, all remaining molecular gas will be exhausted in less than 100 Myr, reaching a final mass in stars of a few 10^11 M_sol. These galaxies are thus proto-ellipticals caught during their major episode of star formation, and observed at the peak epoch z=1.5-3 of the cosmic star formation history of the Universe.

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S_500>100 mJy, 21 are strongly lensed (multiply imaged), 4 are moderately lensed (singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r_half) and far-infrared luminosities (L_FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z_lens>0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500um flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L_FIR (median L_FIR=7.9x10^12 L_sun) and two decades in FIR luminosity surface density (median Sigma_FIR=6.0x10^11 L_sun kpc^-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift.

Several lines of argument support the existence of a link between activity at the nuclei of galaxies, in the form of an accreting supermassive black hole, and star-formation activity in these galaxies. The exact nature of this link is still under debate. Radio jets have long been argued to be an ideal mechanism that allows AGN to interact with their host galaxy and regulate star-formation. In this context, we are using a sample of radio sources in the North Ecliptic Pole (NEP) field to study the nature of the putative link between AGN activity and star-formation. This is done by means of spectral energy distribution (SED) fitting. We use the excellent spectral coverage of the AKARI infrared space telescope together with the rich ancillary data available in the NEP to build SEDs extending from UV to far-IR wavelengths. Through SED fitting we constrain both the AGN and host galaxy components. We find a significant AGN component in our sample of relatively faint radio-sources ($<$mJy), that increases in power with increasing radio-luminosity. At the highest radio-luminosities, the presence of powerful jets dominates the radio emission of these sources. A positive correlation is found between the luminosity of the AGN component and that of star-formation in the host galaxy, independent of the radio luminosity. By contrast, for a given redshift and AGN luminosity, we find that increasing radio-luminosity leads to a decrease in the specific star-formation rate. The most radio-loud AGN are found to lie on the main sequence of star-formation for their respective redshifts. For the first time, such a two-sided feedback process is seen in the same sample. We conclude that radio jets do suppress star-formation in their host galaxies but appear not to totally quench it. Our results therefore support the maintenance nature of "radio-mode" feedback from radio-AGN jets.

We use galaxies from the Herschel-ATLAS survey, and a suite of ancillary simulations based on an isothermal dust model, to study our ability to determine the effective dust temperature, luminosity and emissivity index of 250um selected galaxies in the local Universe (z < 0.5). As well as simple far-infrared SED fitting of individual galaxies based on chi^2 minimisation, we attempt to derive the best global isothermal properties of 13,826 galaxies with reliable optical counterparts and spectroscopic redshifts. Using our simulations, we highlight the fact that applying traditional SED fitting techniques to noisy observational data in the Herschel Space Observatory bands introduces artificial anti-correlation between derived values of dust temperature and emissivity index. This is true even for galaxies with the most robust detections in our sample, making the results hard to interpret. We apply a method to determine the best-fit global values of isothermal effective temperature and emissivity index for z < 0.5 galaxies in H-ATLAS, deriving T = 22.3 +/- 0.1K and Beta = 1.98 +/- 0.02 (or T = 23.5 +/- 0.1K and Beta = 1.82 +/- 0.02 if we attempt to correct for bias by assuming that T and Beta are independent and normally distributed). We use our technique to test for an evolving emissivity index, finding only weak evidence. The median dust luminosity of our sample is log(Ldust/Lsolar) = 10.72 +/- 0.05 which (unlike T) shows little dependence on the choice of Beta used in our analysis, including whether it is variable or fixed. We use a further suite of simulations to emphasise the importance of the H-ATLAS PACS data for deriving dust temperatures at these redshifts, even though they are less sensitive than the SPIRE data. The majority of galaxies detected by H-ATLAS are normal star-forming galaxies, though a substantial minority (~31%) fall in the Luminous Infrared Galaxy category. (Abridged)

We present tentative H-alpha emission line detections of four submillimetre-detected galaxies at z>3.5: the radio galaxies 8C1909+722 and 4C60.07 at signal-to-noise ratios (SNRs) of 3.1 and 2.5, and two submillimetre-selected galaxies (SMGs) near the first of these at SNRs of 10.0 and 2.4, made with the AKARI Space Telescope as part of the FUHYU mission program. These are the highest-redshift H-alpha detections in such galaxies, made possible by AKARI's unique near-infrared spectroscopic capability. The two radio galaxies had known redshifts and surrounding structure, and we have detected broad H-alpha components indicating the presence of dust-shrouded quasars. We conclude that powerful AGNs at z>3.5 occur in peaks of the star-formation density fields, supporting a close connection between stellar mass build-up and black hole mass assembly at this redshift. We also show that 4C60.07 is a binary AGN. The H-alpha detections of the two SMGs are the first redshift determinations for these sources, confirming their physical association around their companion radio galaxy. The H-alpha-derived star formation rates (SFRs) for the SMGs are lower than their far-infrared derived SFRs by a factor of ~10, suggesting a level of dust obscuration similar to that found in studies at ~1<z<2.7.

We report the redshift of HATLAS J132427.0+284452 (hereafter HATLAS J132427), a gravitationally lensed starburst galaxy, the first determined 'blind' by the Herschel Space Observatory. This is achieved via the detection of [C II] consistent with z = 1.68 in a far-infrared spectrum taken with the SPIRE Fourier Transform Spectrometer. We demonstrate that the [C II] redshift is secure via detections of CO J = 2 - 1 and 3 - 2 using the Combined Array for Research in Millimeter-wave Astronomy and the Institut de Radioastronomie Millimetrique's Plateau de Bure Interferometer. The intrinsic properties appear typical of high-redshift starbursts despite the high lensing-amplified fluxes, proving the ability of the FTS to probe this population with the aid of lensing. The blind detection of [C II] demonstrates the potential of the SAFARI imaging spectrometer, proposed for the much more sensitive SPICA mission, to determine redshifts of multiple dusty galaxies simultaneously without the benefit of lensing.

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

This is the first paper in a series that present a multi-wavelength analysis of the archetype Ultra-Luminous InfraRed Galaxy (ULIRG) IRAS FSC10214+4724, a gravitationally lensed, starburst/AGN at z=2.3. Here we present a new lens model and spatially-resolved radio data, as well as a deep HST F160W map. The lens modelling employs a Bayesian Markov Chain Monte Carlo algorithm with extended-source, forward ray-tracing. Using these high resolution HST, MERLIN and VLA maps, the algorithm allows us to constrain the level of distortion to the continuum spectral energy distribution resulting from emission components with differing magnification factors, due to their size and proximity to the caustic. Our lens model finds the narrow line region (NLR), and by proxy the active nucleus, is preferentially magnified. This supports previous claims that preferential magnification could mask the expected polycyclic aromatic hydrocarbon spectral features in the Spitzer mid-infrared spectrum which roughly trace the star-forming regions. Furthermore, we show the arc-to-counter-image flux ratio is not a good estimate of the magnification in this system, despite its common use in the IRAS FSC10214+4724 literature. Our lens modelling suggests magnifications of \mu ~ 15-20+-2 for the HST F814W, MERLIN 1.7 GHz and VLA 8 GHz maps, significantly lower than the canonical values of \mu = 50-100 often used for this system. Systematic errors such as the dark matter density slope and co-location of stellar and dark matter centroids dominate the uncertainties in the lens model at the 40 percent level.

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.

We combine Planck HFI data at 857, 545, 353 & 217GHz with data from WISE, Spitzer, IRAS & Herschel to investigate the properties of a flux limited sample of local star-forming galaxies. A 545GHz flux density limit was chosen so that the sample is 80% complete at this frequency, giving a sample of 234 local galaxies. We investigate the dust emission and star formation properties of the sample via various models & calculate the local dust mass function. Although 1-component modified black bodies fit the dust emission longward of 80um very well (median beta=1.83) the degeneracy between dust temp & beta also means that the SEDs are very well described by a dust emissivity index fixed at beta=2 and 10<T<25 K. Although a second, warmer dust component is required to fit shorter wavelength data, & contributes ~1/3 of the total infrared emission, its mass is negligible. No evidence is found for a very cold (6-10 K) dust component. The temp of the cold dust component is strongly influenced by the ratio of the star formation rate to the total dust mass. This implies, contrary to what is often assumed, that a significant fraction of even the emission from ~20 K dust is powered by ongoing star formation, whether or not the dust itself is associated with star forming clouds or `cirrus'. There is statistical evidence of a free-free contribution to the 217GHz flux densities of <20%. We find a median dust-to-stellar mass ratio of 0.0046; & that this ratio is anti-correlated with galaxy mass. There is good correlation between dust mass & atomic gas mass (median M_d/M_HI = 0.022), suggesting that galaxies that have more dust have more interstellar medium in general. Our derived dust mass function implies a mean dust mass density of the local Universe (for dust within galaxies), of 7.0+-1.4 x 10^5 M_solar/Mpc, significantly greater than that found in the most recent estimate using Herschel data.

Using IRAM PdBI we report the detection of H2O in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel H-ATLAS survey. The sources are detected either in the 2_{02}-1_{11} or 2_{11}-2_{02} H_2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy.km/s. The corresponding apparent luminosities are mu x L_H2O ~ 3-12 x 10^8 Lo, where mu is the lensing magnification factor (3 < mu < 12). These results confirm that H2O lines are among the strongest molecular lines in such galaxies, with intensities almost comparable to those of the high-J CO lines, and same profiles and line widths (200-900 km/s) as the latter. With the current sensitivity of PdBI, H2O can therefore easily be detected in high-z lensed galaxies (with F(500um) > 100 mJy) discovered in the Herschel surveys. Correcting the luminosities for lensing amplification, L_H2O is found to have a strong dependence on the IR luminosity, varying as ~L_IR^{1.2}. This relation which needs to be confirmed with better statistics, may indicate a role of radiative (IR) excitation of the H2O lines, and implies that high-z galaxies with L_IR >~ 10^13 Lo tend to be very strong emitters in H2O, that have no equivalent in the local universe.

We present a measurement of the angular power spectrum of the cosmic far-infrared background (CFIRB) anisotropies in one of the extragalactic fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) at 250, 350 and 500 \mu m bands. Consistent with recent measurements of the CFIRB power spectrum in Herschel-SPIRE maps, we confirm the existence of a clear one-halo term of galaxy clustering on arcminute angular scales with large-scale two-halo term of clustering at 30 arcminutes to angular scales of a few degrees. The power spectrum at the largest angular scales, especially at 250 \mu m, is contaminated by the Galactic cirrus. The angular power spectrum is modeled using a conditional luminosity function approach to describe the spatial distribution of unresolved galaxies that make up the bulk of the CFIRB. Integrating over the dusty galaxy population responsible for the background anisotropies, we find that the cosmic abundance of dust, relative to the critical density, to be between \Omega_dust=10^{-6} and 8 x 10^{-6} in the redshift range z ~ 0-3. This dust abundance is consistent with estimates of the dust content in the Universe using quasar reddening and magnification measurements in the SDSS.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

The Planck Early Release Compact Source Catalog (ERCSC) has offered the first opportunity to accurately determine the luminosity function of dusty galaxies in the very local Universe (i.e. distances <~ 100 Mpc), at several (sub-)millimetre wavelengths, using blindly selected samples of low redshift sources, unaffected by cosmological evolution. This project, however, requires careful consideration of a variety of issues including the choice of the appropriate flux density measurement, the separation of dusty galaxies from radio sources and from Galactic sources, the correction for the CO emission, the effect of density inhomogeneities, and more. We present estimates of the local luminosity functions at 857 GHz (350 microns), 545 GHz (550 microns) and 353 GHz (850 microns) extending across the characteristic luminosity L_star, and a preliminary estimate over a limited luminosity range at 217 GHz (1382 microns). At 850 microns and for luminosities L >~ L_star our results agree with previous estimates, derived from the SCUBA Local Universe Galaxy Survey (SLUGS), but are higher than the latter at L <~ L_star. We also find good agreement with estimates at 350 and 500 microns based on preliminary Herschel survey data.

We have measured the clustering properties of low-redshift (z < 0.3) sub-mm galaxies detected at 250 micron in the Herschel-ATLAS Science Demonstration Phase (SDP) field. We selected a sample for which we have high-quality spectroscopic redshifts, obtained from reliably matching the 250-micron sources to a complete (for r < 19.4) sample of galaxies from the GAMA database. Both the angular and spatial clustering strength are measured for all z < 0.3 sources as well as for five redshift slices with thickness delta z=0.05 in the range 0.05 < z < 0.3. Our measured spatial clustering length r_0 is comparable to that of optically-selected, moderately star-forming (blue) galaxies: we find values around 5 Mpc. One of the redshift bins contains an interesting structure, at z = 0.164.

We present a photometric catalog of infrared (IR) sources based on the North Ecliptic PoleWide field (NEP-Wide) survey of AKARI, which is an infrared space telescope launched by Japan. The NEP-Wide survey covered 5.4 deg2 area, a nearly circular shape centered on the North Ecliptic Pole, using nine photometric filter-bands from 2 - 25 {\mu}m of the Infrared Camera (IRC). Extensive efforts were made to reduce possible false objects due to cosmic ray hits, multiplexer bleeding phenomena around bright sources, and other artifacts. The number of detected sources varied depending on the filter band: with about 109,000 sources being cataloged in the near-IR bands at 2 - 5 {\mu}m, about 20,000 sources in the shorter parts of the mid-IR bands between 7 - 11 {\mu}m, and about 16,000 sources in the longer parts of the mid-IR bands, with \sim 4,000 sources at 24 {\mu}m. The estimated 5? detection limits are approximately 21 magnitude (mag) in the 2 - 5 {\mu}m bands, 19.5 - 19 mag in the 7 - 11 {\mu}m, and 18.8 - 18.5 mag in the 15 - 24 {\mu}m bands in the AB magnitude scale. The completenesses for those bands were evaluated as a function of magnitude: the 50% completeness limits are about 19.8 mag at 3 {\mu}m, 18.6 mag at 9 {\mu}m, and 18 mag at 18 {\mu}m band, respectively. To construct a reliable source catalog, all of the detected sources were examined by matching them with those in other wavelength data, including optical and ground-based near-IR bands. The final band-merged catalog contains about 114,800 sources detected in the IRC filter bands. The properties of the sources are presented in terms of the distributions in various color-color diagrams.

AKARI's all-sky survey resolves the far-infrared emission in many thousands of nearby galaxies, providing essential local benchmarks against which the evolution of high-redshift populations can be measured. This review presents some recent results in the resolved galaxy populations, covering some well-known nearby targets, as well as samples from major legacy surveys such as the Herschel Reference Survey and the JCMT Nearby Galaxies Survey. This review also discusses the prospects for higher redshifts surveys, including strong gravitational lens clusters and the AKARI NEP field.

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

The results of a deep 20 cm radio survey at 20 cm are reported of the AKARI Deep Field South (ADF-S) near the South Ecliptic Pole (SEP), using the Australia Telescope Compact Array telescope, ATCA. The survey has 1 sigma detection limits ranging from 18.7--50 microJy per beam over an area of ~1.1 sq degrees, and ~2.5 sq degrees to lower sensitivity. The observations, data reduction and source count analysis are presented, along with a description of the overall scientific objectives, and a catalogue containing 530 radio sources detected with a resolution of 6.2" x 4.9". The derived differential source counts show a pronounced excess of sources fainter than ~1 mJy, consistent with an emerging population of star forming galaxies. Cross-correlating the radio with AKARI sources and archival data we find 95 cross matches, with most galaxies having optical R-magnitudes in the range 18-24 mag, and 52 components lying within 1" of a radio position in at least one further catalogue (either IR or optical). We have reported redshifts for a sub-sample of our catalogue finding that they vary between galaxies in the local universe to those having redshifts of up to 0.825. Associating the radio sources with the Spitzer catalogue at 24 microns, we find 173 matches within one Spitzer pixel, of which a small sample of the identifications are clearly radio loud compared to the bulk of the galaxies. The radio luminosity plot and a colour-colour analysis suggest that the majority of the radio sources are in fact luminous star forming galaxies, rather than radio-loud AGN. There are additionally five cross matches between ASTE or BLAST submillimetre galaxies and radio sources from this survey, two of which are also detected at 90 microns, and 41 cross-matches with submillimetre sources detected in the Herschel HerMES survey Public Data release.

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

Bright submm-selected galaxies have been found to be a rich source of strong gravitational lenses. However, strong gravitational lensing of extended sources leads inevitably to differential magnification. In this paper I quantify the effect of differential magnification on simulated far-infrared and submm surveys of strong gravitational lenses, using a foreground population of Navarro-Frenk-White plus de Vaucouleurs' density profiles, with a model source resembling the Cosmic Eyelash and QSO J1148+5251. Some emission line diagnostics are surprisingly unaffected by differential magnification effects: for example, the bolometric fractions of [C II] 158um and CO(J=1-0), often used to infer densities and ionisation parameters, have typical differential magnification effects that are smaller than the measurement errors. However, the CO ladder itself is significantly affected. Far-infrared lensed galaxy surveys (e.g. at 60um) strongly select for high-redshift galaxies with caustics close to AGN, boosting the apparent bolometric contribution of AGN. The lens configuration of IRAS F10214+4724 is naturally explained in this context. Conversely, submm/mm-wave surveys (e.g. 500-1400um) strongly select for caustics close to knots of star formation boosting the latter's bolometric fraction. In general, estimates of bolometric fractions from spectral energy distributions of strongly lensed infrared galaxies are so unreliable as to be useless, unless a lens mass model is available to correct for differential magnification.

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope.

We present the first far-infrared luminosity function in the AKARI Deep Field South, a premier deep field of the AKARI Space Telescope, using spectroscopic redshifts obtained with AAOmega. To date, we have found spectroscopic redshifts for 389 galaxies in this field and have measured the local (z < 0.25) 90 micron luminosity function using about one-third of these redshifts. The results are in reasonable agreement with recent theoretical predictions.

While the selection of strongly lensed galaxies with 500{\mu}m flux density S_500>100 mJy has proven to be rather straightforward (Negrello et al. 2010), for many applications it is important to analyze samples larger than the ones obtained when confining ourselves to such a bright limit. Moreover, only by probing to fainter flux densities is possible to exploit strong lensing to investigate the bulk of the high-z star-forming galaxy population. We describe HALOS (the Herschel-ATLAS Lensed Objects Selection), a method for efficiently selecting fainter candidate strongly lensed galaxies, reaching a surface density of ~1.5-2 deg^-2, i.e. a factor of about 4 to 6 higher than that at the 100 mJy flux limit. HALOS will allow the selection of up to ~1000 candidate strongly lensed galaxies (with amplifications \mu>2) over the full H-ATLAS survey area. Applying HALOS to the H-ATLAS Science Demonstration Phase field (~14.4 deg^2) we find 31 candidate strongly lensed galaxies, whose candidate lenses are identified in the VIKING near-infrared catalog. Using the available information on candidate sources and candidate lenses we tentatively estimate a ~72% purity of the sample. The redshift distribution of the candidate lensed sources is close to that reported for most previous surveys for lensed galaxies, while that of candidate lenses extends to substantially higher redshifts than found in the other surveys. The counts of candidate strongly lensed galaxies are also in good agreement with model predictions (Lapi et al. 2011). Even though a key ingredient of the method is the deep near-infrared VIKING photometry, we show that H-ATLAS data alone allow the selection of a similarly deep sample of candidate strongly lensed galaxies with an efficiency close to 50%; a slightly lower surface density (~1.45 deg^-2) can be reached with a ~70% efficiency.

We present a new catalogue of mid-IR sources using the AKARI NEP-Deep survey. The InfraRed Camera (IRC) onboard AKARI has a comprehensive mid-IR wavelength coverage with 9 photometric bands at 2 - 24 micron. We utilized all of these bands to cover a nearly circular area adjacent to the North Ecliptic Pole (NEP). We designed the catalogue to include most of sources detected in 7, 9, 11, 15 and 18 micron bands, and found 7284 sources in a 0.67 deg^2 area. From our simulations, we estimate that the catalogue is ~80 per cent complete to 200 micro Jy at 15 - 18 micron, and ~10 per cent of sources are missed, owing to source blending. Star-galaxy separation is conducted using only AKARI photometry, as a result of which 10 per cent of catalogued sources are found to be stars. The number counts at 11, 15, 18, and 24 micron are presented for both stars and galaxies. A drastic increase in the source density is found in between 11 and 15 micron at the flux level of ~300 micro Jy. This is likely due to the redshifted PAH emission at 8 micron, given our rough estimate of redshifts from an AKARI colour-colour plot. Along with the mid-IR source catalogue, we present optical-NIR photometry for sources falling inside a Subaru/Sprime-cam image covering part of the AKARI NEP-Deep field, which is deep enough to detect most of AKARI mid-IR sources, and useful to study optical characteristics of a complete mid-IR source sample.

We have studied the star-formation and AGN activity of massive galaxies in the redshift range $z=0.4-2$, which are detected in a deep survey field using the AKARI InfraRed (IR) astronomical satellite and {\em Subaru} telescope toward the North Ecliptic Pole (NEP). The AKARI/IRC Mid-InfraRed (MIR) multiband photometry is used to trace their star-forming activities with the Polycyclic-Aromatic Hydrocarbon (PAH) emissions, which is also used to distinguish star-forming populations from AGN dominated ones and to estimate the Star Formation Rate (SFR) derived from their total emitting IR (TIR) luminosities. In combination with analyses of their stellar components, we have studied the MIR SED features of star-forming and AGN-harboring galaxies.

I summarise recent results on multi-wavelength properties of distant lensed galaxies, with a particular focus on Herschel. Submm surveys have already resulted in a breakthrough discovery of an extremely efficient selection technique for strong gravitational lenses. Benefitting from the gravitational magnification boost, blind mm-wave redshifts have been demonstrated on IRAM, SMA and GBT, and follow-up emission line detections have been made of water, [OIII], [CII] and other species, revealing the PDR/XDR/CRDR conditions. I also discuss HST imaging of submm lenses, lensed galaxy reconstruction, the prospects for ALMA and e-Merlin and the effects of differential magnification. Many emission line diagnostics are relatively unaffected by differential magnification, but SED-based estimates of bolometric fractions in lensed infrared galaxies are so unreliable as to be useless, unless a lens mass model is available to correct for differential amplification.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

Exploiting the Herschel-ATLAS Science Demonstration Phase (SDP) survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 micron and at several redshifts z>1, for bright sub-mm galaxies with star formation rates (SFR) >100 M_sun/yr. We find that the evolution of the comoving LF is strong up to z~2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M_H>3*10^12 M_sun) amounts to ~7*10^8 yr. Given the SFRs, which are in the range 10^2-10^3 M_sun/yr, this timescale implies final stellar masses of order of 10^11-10^12 M_sun. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z>1. The comparison of the statistics for sub-mm and UV selected galaxies suggests that the dust-free, UV bright phase, is >10^2 times shorter than the sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference Spectral Energy Distribution (SED; the one of the z~2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also the counts at wavelengths ranging from 250 micron to ~1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs of z>1 galaxies around the reference one is rather small.

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

We present the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel-ATLAS. The sample consists of galaxies selected at 250{\mu}m which have reliable counterparts from SDSS at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single temperature grey-body model for the spectral energy distribution and also using a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z = 0.4 - 0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust-to-stellar mass ratio was about 3-4 times larger, and the optical depth derived from fitting the UV-sub-mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption timescale together with either a more top-heavy IMF, efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall ISM mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

We present a far-infrared galaxy luminosity function for the local universe. We have obtained 389 spectroscopic redshifts for galaxies observed at 90 microns in the AKARI Deep Field South, using the AAOmega fibre spectrograph via optical identifications in the digitized sky survey and 4m-class optical imaging. For the luminosity function presented in this paper, we have used those galaxies which have redshifts 0<z<0.25, have optical magnitudes and are not part of a newly discovered cluster of galaxies (giving a total of 130 sources). Infrared and optical completeness functions were estimated using earlier Spitzer data and APM B-band optical data respectively, and the luminosity function has been prepared using the 1/Vmax method. We also separate the luminosity function between galaxies which show evidence of predominantly star-forming activity and predominantly active galactic nucleus (AGN) activity in their optical spectra. Our luminosity function is in good agreement with the previous 90 micron luminosity function from the European Large Area ISO Survey, and we also present a luminosity function with combined AKARI and ISO data. The result is in reasonable agreement with predictions based on the earlier IRAS 60 micron PSC-z catalogue, and with a recent backward evolution model.

We report the detection of a significant excess in the surface density of far-infrared sources from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) within ~1 Mpc of the centres of 66 optically-selected clusters of galaxies in the SDSS with <z>~0.25. From the analysis of the multiwavelength properties of their counterparts we conclude that the far-infrared emission is associated with dust-obscured star formation and/or active galactic nuclei within galaxies in the clusters themselves. The excess reaches a maximum at a radius of ~0.8 Mpc, where we find 1.0\pm0.3 S_250um>34 mJy sources on average per cluster above what would be expected for random field locations. If the far-infrared emission is dominated by star formation (as opposed to AGN) then this corresponds to an average star formation rate of ~7 M_sun/yr per cluster in sources with L_IR>5d10 L_sun. Although lensed sources make a negligible contribution to the excess signal, a fraction of the sources around the clusters could be gravitationally lensed, and we have identified a sample of potential cases of cluster-lensed Herschel sources that could be targeted in follow-up studies.

The Herschel-ATLAS is a survey of 550 square degrees with the Herschel Space Observatory in five far--infrared and submillimetre bands. The first data for the survey, observations of a field 4x4 sq. degrees in size, were taken during the Science Demonstration Phase, and reach a 5 sigma noise level of 33 mJy/beam at 250 microns. This paper describes the source extraction methods used to create the corresponding Science Demonstration Phase catalogue, which contains 6876 sources, selected at 250 microns, within ~14 sq. degrees. SPIRE sources are extracted using a new method specifically developed for Herschel data; PACS counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 microns have undergone flux density enhancements of up to a factor of ~2, due mainly to source confusion. Correction factors are determined for these effects. The SDP dataset and corresponding catalogue will be available from http://www.h-atlas.org/.

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 use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated.

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

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

We present a first complete 12 CO J=3-2 map of M81, observed as part of the Nearby Galaxies Legacy Survey being carried out at the James Clerk Maxwell Telescope. We detect 9 regions of significant CO emission located at different positions within the spiral arms, and confirm that the global CO emission in the galaxy is low. We combine these data with a new H-alpha map obtained using the Isaac Newton Telescope and archival HI, 24 microns and FUV images to uncover a correlation between the molecular gas and star forming regions in M81. For the nine regions detected in CO J=3-2, we combine our CO J=3-2 data with existing CO J=1-0 data to calculate line ratios. We find that the ratio J=(3-2)/(1-0) is in agreement with the range of typical values found in the literature (0.2-0.8). Making reasonable assumptions, this allows us to constrain the hydrogen density to the range (10^3-10^4) cm^{-3}. We also estimated the amount of hydrogen produced in photo-dissociation regions near the locations where CO J=3-2 was detected.

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

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) has uncovered a population of strongly-lensed submillimeter galaxies (SMGs). The Zpectrometer instrument on the Green Bank Telescope (GBT) was used to measure the redshifts and constrain the masses of the cold molecular gas reservoirs for two candidate high-redshift lensed sources. We derive CO(1-0) redshifts of z=3.042+/-0.001 and z=2.625+/-0.001, and measure molecular gas masses of (1--3)x10^{10}Msun, corrected for lens amplification and assuming a conversion factor of alhpa=0.8 Msun(K km/s pc^2)^{-1}. We find typical L(IR)/L'(CO) ratios of 120+/-40 and 140+/-50 Lsun (K km/s pc^2)^{-1}, which are consistent with those found for local ULIRGs and other high-redshift SMGs. From analysis of published data, we find no evidence for enhanced L(IR)/L'(CO(1-0)) ratios for the SMG population in comparison to local ULIRGs. The GBT results highlight the power of using the CO lines to derive blind redshifts, which is challenging for the SMG population at optical wavelengths given their high obscuration.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

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

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

We have made the first map of CO(J=3-2) emission covering the disk of the edge-on galaxy, NGC~4631, which is known for its spectacular gaseous halo. The strongest emission, which we model with a Gaussian ring,occurs within a radius of 5 kpc. Weaker disk emission is detected out to radii of 12 kpc, the most extensive molecular component yet seen in this galaxy. From comparisons with infrared data, we find that CO(J=3-2) emission more closely follows the hot dust component, rather than the cold dust,consistent with it being a good tracer of star formation. The first maps of $R_{3-2/1-0}$, H$_2$ mass surface density and SFE have been made for the inner 2.4 kpc radius region. Only 20\% of the SF occurs in this region and excitation conditions are typical of galaxy disks, rather than of central starbursts. The SFE suggests long gas consumption timescales ($>$ $10^9$ yr). The velocity field is dominated by a steeply rising rotation curve in the region of the central molecular ring followed by a flatter curve in the disk. A very steep gradient in the rotation curve is observed at the nucleus, providing the first evidence for a central concentration of mass: M$_{dyn}\,=\,5\,\times\,10^7$ M$_\odot$ within a radius of 282 pc. The velocity field shows anomalous features indicating the presence of molecular outflows; one of them is associated with a previously observed CO(J=1-0) expanding shell. Consistent with these outflows is the presence of a thick ($z$ up to $1.4$ kpc) CO(J=3-2) disk. We suggest that the interaction between NGC~4631 and its companion(s) has agitated the disk and also initiated star formation which was likely higher in the past than it is now. These may be necessary conditions for seeing prominent halos.

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

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.

We present extragalactic number counts and a lower limit estimate for the cosmic infrared background at 15 um from AKARI ultra deep mapping of the gravitational lensing cluster Abell 2218. This data is the deepest taken by any facility at this wavelength, and uniquely samples the normal galaxy population. We have de-blended our sources, to resolve photometric confusion, and de-lensed our photometry to probe beyond AKARI's blank-field sensitivity. We estimate a de-blended 5 sigma sensitivity of 28.7 uJy. The resulting 15 um galaxy number counts are a factor of three fainter than previous results, extending to a depth of ~ 0.01 mJy and providing a stronger lower limit constraint on the cosmic infrared background at 15 um of 1.9 +/- 0.5 nW m^-2 sr^-1.

The Westerbork Radio Synthesis Telescope, WSRT, has been used to make a deep radio survey of an ~ 1.7 sq degree field coinciding with the AKARI North Ecliptic Pole Deep Field. The observations, data reduction and source count analysis are presented, along with a description of the overall scientific objectives. The survey consisted of 10 pointings, mosaiced with enough overlap to maintain a similar sensitivity across the central region that reached as low as 21 microJy per beam at 1.4 GHz. A catalogue containing 462 sources detected with a resolution of 17"x15" is presented. The differential source counts calculated from the WSRT data have been compared with those from the shallow VLA-NEP survey of Kollgaard et al 1994, and show a pronounced excess for sources fainter than ~ 1 mJy, consistent with the presence of a population of star forming galaxies at sub-mJy flux levels. The AKARI North Ecliptic Pole Deep field is the focus of a major observing campaign conducted across the entire spectral region. The combination of these data sets, along with the deep nature of the radio observations will allow unique studies of a large range of topics including the redshift evolution of the luminosity function of radio sources, the clustering environment of radio galaxies, the nature of obscured radio-loud active galactic nuclei, and the radio/far-infrared correlation for distant galaxies. This catalogue provides the basic data set for a future series of paper dealing with source identifications, morphologies, and the associated properties of the identified radio sources.

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

Aims: We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods: We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results: We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

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 colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

We present measurements of the angular correlation function of galaxies selected from the first field of the H-ATLAS survey. Careful removal of the background from galactic cirrus is essential, and currently dominates the uncertainty in our measurements. For our 250 micron-selected sample we detect no significant clustering, consistent with the expectation that the 250 micron-selected sources are mostly normal galaxies at z<~ 1. For our 350 micron and 500 micron-selected samples we detect relatively strong clustering with correlation amplitudes A of 0.2 and 1.2 at 1', but with relatively large uncertainties. For samples which preferentially select high redshift galaxies at z~2-3 we detect significant strong clustering, leading to an estimate of r_0 ~ 7-11 h^{-1} Mpc. The slope of our clustering measurements is very steep, delta~2. The measurements are consistent with the idea that sub-mm sources consist of a low redshift population of normal galaxies and a high redshift population of highly clustered star-bursting galaxies.

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 use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

The AKARI FU-HYU mission program carried out mid-infrared imaging of several well studied Spitzer fields preferentially selecting fields already rich in multi-wavelength data from radio to X-ray wavelengths filling in the wavelength desert between the Spitzer IRAC and MIPS bands.We present the initial results for the FU-HYU survey in the GOODS-N field.We utilize the supreme multiwavelength coverage in the GOODS-N field to produce a multiwavelength catalogue from infrared to ultraviolet wavelengths, containing more than 4393 sources, including photometric redshifts. Using the FU-HYU catalogue we present colour-colour diagrams that map the passage of PAH features through our observation bands. We find that the longer mid-infrared bands from AKARI (IRC-L18W 18 micron band) and Spitzer (MIPS24 24 micron band) provide an accurate measure of the total MIR emission of the sources and therefore their probable total mid-infrared luminosity. We also find that colours incorporating the AKARI IRC-S11 11 micron band produce a bimodal distribution where an excess at 11 microns preferentially selects moderate redshift star-forming galaxies. These powerful colour-colour diagnostics are further used as tools to extract anomalous colour populations, in particular a population of Silicate Break galaxies from the GOODS-N field showing that dusty starbursts can be selected of specific redshift ranges (z=1.2 - 1.6) by mid-infrared drop-out techniques. The FU-HYU catalogue will be made publically available to the astronomical community.

We present galaxy counts at 15 microns using the Japanese AKARI satellite's NEP-deep and NEP-wide legacy surveys at the North Ecliptic Pole. The total number of sources detected are approximately 6700 and 10,700 down to limiting fluxes of 117 and 250 microJy (5 sigma) for the NEP-deep and NEP-wide survey respectively. We construct the Euclidean normalized differential source counts for both data sets (assuming 80 percent completeness levels of 200 and 270 microJy respectively) to produce the widest and deepest contiguous survey at 15 microns to date covering the entire flux range from the deepest to shallowest surveys made with the Infrared Space Observatory (ISO) over areas sufficiently significant to overcome cosmic variance, detecting six times as many sources as the largest survey carried out with ISO.We compare the results from AKARI with the previous surveys with ISO at the same wavelength and the Spitzer observations at 16 microns using the peek-up camera on its IRS instrument. The AKARI source counts are consistent with other results to date reproducing the steep evolutionary rise at fluxes less than a milliJansky and super-Euclidean slopes. We find the the AKARI source counts show a slight excess at fluxes fainter than 200 microJanskys which is not predicted by previous source count models at 15 microns. However, we caution that at this level we may be suffering from the effects of source confusion in our data. At brighter fluxes greater than a milliJansky, the NEP-wide survey source counts agree with the Northern ISO-ELAIS field results, resolving the discrepancy of the bright end calibration in the ISO 15 micron source counts.

We present the results of large-area CO J=3-2 emission mapping of three nearby field galaxies, NGC 628, NGC 3521, and NGC 3627, completed at the James Clerk Maxwell Telescope as part of the Nearby Galaxies Legacy Survey. These galaxies all have moderate to strong CO J=3-2 detections over large areas of the fields observed by the survey, showing resolved structure and dynamics in their warm/dense molecular gas disks. All three galaxies were part of the Spitzer Infrared Nearby Galaxies Survey sample, and as such have excellent published multi-wavelength ancillary data. These data sets allow us to examine the star formation properties, gas content, and dynamics of these galaxies on sub-kiloparsec scales. We find that the global gas depletion times for dense/warm molecular gas in these galaxies is consistent with other results for nearby spiral galaxies, indicating this may be independent of galaxy properties such as structures, gas compositions, and environments. Similar to the results from the THINGS HI survey, we do not see a correlation of the star formation efficiency with the gas surface density consistent with the Schmidt-Kennicutt law. Finally, we find that the star formation efficiency of the dense molecular gas traced by CO J=3-2 is potentially flat or slightly declining as a function of molecular gas density, the CO J=3-2/J=1-0 ratio (in contrast to the correlation found in a previous study into the starburst galaxy M83), and the fraction of total gas in molecular form.

We present early results from our multi-wavelength follow-up campaigns of the AKARI Deep Fields at the North and South Ecliptic Poles. We summarize our campaigns in this poster paper, and present three early outcomes. (a) Our AAOmega optical spectroscopy of the Deep Field South at the AAT has observed over 550 different targets, and our preliminary local luminosity function at 90 microns from the first four hours of data is in good agreement with the predictions from Serjeant & Harrison 2005. (b) Our GMRT 610 MHz imaging in the Deep Field North has reached ~30 microJy RMS, making this among the deepest images at this frequency. Our 610 MHz source counts at >200 microJy are the deepest ever derived at this frequency. (c) Comparing our GMRT data with our 1.4 GHz WSRT data, we have found two examples of radio-loud AGN that may have more than one epoch of activity.

We report the detection and measurement of the absolute brightness and spatial fluctuations of the cosmic infrared background (CIB) with the AKARI satellite. We have carried out observations at 65, 90, 140 and 160 um as a cosmological survey in AKARI Deep Field South (ADF-S), which is one of the lowest cirrus regions with contiguous area on the sky. After removing bright galaxies and subtracting zodiacal and Galactic foregrounds from the measured sky brightness, we have successfully measured the CIB brightness and its fluctuations across a wide range of angular scales from arcminutes to degrees. The measured CIB brightness is consistent with previous results reported from COBE data but significantly higher than the lower limits at 70 and 160 um obtained with the Spitzer satellite from the stacking analysis of 24-um selected sources. The discrepancy with the Spitzer result is possibly due to a new galaxy population at high redshift obscured by hot dust. From power spectrum analysis at 90 um, three components are identified: shot noise due to individual galaxies; Galactic cirrus emission dominating at the largest angular scales of a few degrees; and an additional component at an intermediate angular scale of 10-30 arcminutes, possibly due to galaxy clustering. The spectral shape of the clustering component at 90 um is very similar to that at longer wavelengths as observed by Spitzer and BLAST. Moreover, the color of the fluctuations indicates that the clustering component is as red as Ultra-luminous infrared galaxies (ULIRGs) at high redshift, These galaxies are not likely to be the majority of the CIB emission at 90 um, but responsible for the clustering component. Our results provide new constraints on the evolution and clustering properties of distant infrared galaxies.

Using an AKARI multi-wavelength mid-infrared (IR) survey, we identify luminous starburst galaxies at z> 0.5 based on the PAH luminosity, and investigate the nature of these PAH-selected starbursts. An extragalactic survey with AKARI towards the north ecliptic pole (NEP), the NEP-Deep survey, is unique in terms of a comprehensive wavelength coverage from 2 to 24um using all 9 photometric bands of the InfraRed Camera (IRC). This survey allows us to photometrically identify galaxies whose mid-IR emission is clearly dominated by PAHs. We propose a single colour selection method to identify such galaxies, using two mid-IR flux ratios at 11-to-7um and 15-to-9um (PAH-to-continuum flux ratio in the rest-frame), which are useful to identify starburst galaxies at z~0.5 and 1, respectively. We perform a fitting of the spectral energy distributions (SEDs) from optical to mid-IR wavelengths, using an evolutionary starburst model with a proper treatment of radiative transfer (SBURT), in order to investigate their nature. The SBURT model reproduces observed optical-to-mid-IR SEDs of more than a half of PAH-selected galaxies. Based on the 8um luminosity, we find ultra luminous infrared galaxies (ULIRGs) among PAH-selected galaxies. Their PAH luminosity is higher than local ULIRGs with a similar luminosity, and the PAH-to-total IR luminosity ratio is consistent with that of less luminous starburst galaxies. They are a unique galaxy population at high redshifts and we call these PAH-selected ULIRGs "PAH-luminous" galaxies. Although they are not as massive as submillimetre galaxies at z~2, they have the stellar mass of >3x10^{10} Msun and therefore moderately massive.

We present a comparison of the SCUBA Half Degree Extragalactic Survey (SHADES) at 450, 850 and 1100 microns with deep guaranteed time 15 microns AKARI FU-HYU survey data and Spitzer guaranteed time data at 3.6-24 microns in the Lockman Hole East. The AKARI data was analysed using bespoke software based in part on the drizzling and minimum-variance matched filtering developed for SHADES, and was cross calibrated against Infrared Space Observatory (ISO) fluxes. Our stacking analyses find AKARI 15um galaxies with >~200 microJy contribute >10% of the 450 micron background, but only <4% of the 1100 micron background, suggesting that different populations contribute at mm-wavelengths. We confirm our earlier result that the ultra-deep 450 micron SCUBA-2 Cosmology Survey will be dominated by populations already detected by AKARI and Spitzer mid-infrared surveys. The superb mid-infrared wavelength coverage afforded by combining Spitzer and AKARI photometry is an excellent diagnostic of AGN contributions, and we find that (23-52)% of submm-selected galaxies have AGN bolometric fractions f_AGN>0.3.

In this brief, partial, incomplete and egregiously self-citing review I will summarise some of the key results in the past few years in surveys for dusty star-forming galaxies and some exciting prospects for forthcoming surveys.

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

We aim to reveal environmental dependence of infrared luminosity functions (IR LFs) of galaxies at z~0.8 using the AKARI satellite. We construct restframe 8um IR LFs in the cluster region RXJ1716.4+6708 at z=0.81, and compare them with a blank field using the AKARI North Ecliptic Pole deep field data at the same redshift. AKARI's wide field of view (10'x10') is suitable to investigate wide range of galaxy environments. AKARI's 15um filter is advantageous here since it directly probes restframe 8um at z~0.8, without relying on a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. We have found that cluster IR LFs at restframe 8um have a factor of 2.4 smaller L^* and a steeper faint-end slope than that of the field. Confirming this trend, we also found that faint-end slopes of the cluster LFs becomes flatter and flatter with decreasing local galaxy density. These changes in LFs cannot be explained by a simple infall of field galaxy population into a cluster. Physics that can preferentially suppress IR luminous galaxies in high density regions is required to explain the observed results.

We present photometry, photometric redshifts and extra galactic number counts for ultra deep 15 micron mapping of the gravitational lensing cluster Abell 2218 (A2218), which is the deepest image taken by any facility at this wavelength. This data resolves the cosmic infrared background (CIRB) beyond the 80% that blank field AKARI surveys aim to achieve. To gain an understanding of galaxy formation and evolution over the age of the Universe a necessary step is to fully resolve the CIRB, which represents the dust-shrouded cosmic star formation history. Observing through A2218 gives magnifications of up to a factor of 10, thus allowing the sampling of a more representative spread of high redshift galaxies, which comprise the bulk of the CIRB. 19 pointed observations were taken by AKARIs IRC MIR-L channel, and a final combined image with an area of 122.3 square arcminutes and effective integration time of 8460 seconds was achieved. The 5 sigma sensitivity limit is estimated at 41.7 uJy. An initial 5 sigma catalogue of 565 sources was extracted giving 39 beams per source, which shows the image is confusion limited. Our 15 micron number counts show strong evolution consistent with galaxy evolution models that incorporate downsizing in star formation.

We used 3.6, 8.0, 70, 160 micron Spitzer Space Telescope data, James Clerk Maxwell Telescope HARP-B CO J=(3-2) data, National Radio Astronomy Observatory 12 meter telescope CO J=(1-0) data, and Very Large Array HI data to investigate the relations among PAHs, cold (~20 K) dust, molecular gas, and atomic gas within NGC 2403, an SABcd galaxy at a distance of 3.13 Mpc. The dust surface density is mainly a function of the total (atomic and molecular) gas surface density and galactocentric radius. The gas-to-dust ratio monotonically increases with radius, varying from ~100 in the nucleus to ~400 at 5.5 kpc. The slope of the gas-to-dust ratio is close to that of the oxygen abundance, suggesting that metallicity strongly affects the gas-to-dust ratio within this galaxy. The exponential scale length of the radial profile for the CO J=(3-2) emission is statistically identical to the scale length for the stellar continuum-subtracted 8 micron (PAH 8 micron) emission. However, CO J=(3-2) and PAH 8 micron surface brightnesses appear uncorrelated when examining sub-kpc sized regions.

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.

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

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.

We present the {\it AKARI} InfraRed Camera (IRC) imaging observation of early-type galaxies in A2218 at z $\simeq$ 0.175. Mid-infrared (MIR) emission from early-type galaxies traces circumstellar dust emission from AGB stars or/and residual star formation. Including the unique imaging capability at 11 and 15 $\mu$m, our {\it AKARI} data provide an effective way to investigate MIR properties of early-type galaxies in the cluster environment. Among our flux-limited sample of 22 red sequence early-type galaxies with precise dynamical and line strength measurements ($<$ 18 mag at 3 $\mu m$), we find that at least 41% have MIR-excess emission. The $N3-S11$ versus $N3$ (3 and 11 $\mu$m) color-magnitude relation shows the expected blue sequence, but the MIR-excess galaxies add a red wing to the relation especially at the fainter end. A SED analysis reveals that the dust emission from AGB stars is the most likely cause for the MIR-excess, with low level of star formation being the next possible explanation. The MIR-excess galaxies show a wide spread of $N3-S11$ colors, implying a significant spread (2--11 Gyr) in the estimated mean ages of stellar populations. We study the environmental dependence of MIR-excess early-type galaxies over an area out to a half virial radius ($\sim$1 Mpc). We find that the MIR-excess early-type galaxies are preferentially located in the outer region. From these evidences, we suggest that the fainter, MIR-excess early-type galaxies have just joined the red sequence, possibly due to the infall and subsequent morphological/spectral transformation induced by the cluster environment.

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 present large-area maps of the CO J=3-2 emission obtained at the James Clerk Maxwell Telescope for four spiral galaxies in the Virgo Cluster. We combine these data with published CO J=1-0, 24 micron, and Halpha images to measure the CO line ratios, molecular gas masses, and instantaneous gas depletion times. For three galaxies in our sample (NGC 4254, NGC4321, and NGC 4569), we obtain molecular gas masses of 7E8-3E9 Msun and disk-averaged instantaneous gas depletion times of 1.1-1.7 Gyr. We argue that the CO J=3-2 line is a better tracer of the dense star forming molecular gas than the CO J=1-0 line, as it shows a better correlation with the star formation rate surface density both within and between galaxies. NGC 4254 appears to have a larger star formation efficiency(smaller gas depletion time), perhaps because it is on its first passage through the Virgo Cluster. NGC 4569 shows a large-scale gradient in the gas properties traced by the CO J=3-2/J=1-0 line ratio, which suggests that its interaction with the intracluster medium is affecting the dense star-forming portion of the interstellar medium directly. The fourth galaxy in our sample, NGC 4579, has weak CO J=3-2 emission despite having bright 24 micron emission; however, much of the central luminosity in this galaxy may be due to the presence of a central AGN.

We investigate the photometric redshift accuracy achievable with the AKARI infrared data in deep multi-band surveys, such as in the North Ecliptic Pole field. We demonstrate that the passage of redshifted policyclic aromatic hydrocarbons and silicate features into the mid-infrared wavelength window covered by AKARI is a valuable means to recover the redshifts of starburst galaxies. To this end we have collected a sample of ~60 galaxies drawn from the GOODS-North Field with spectroscopic redshift 0.5<~z_spec<~1.5 and photometry from 3.6 to 24 micron, provided by the Spitzer, ISO and AKARI satellites. The infrared spectra are fitted using synthetic galaxy Spectral Energy Distributions which account for starburst and active nuclei emission. For ~90% of the sources in our sample the redshift is recovered with an accuracy |z_phot-z_spec|/(1+z_spec)<~10%. A similar analysis performed on different sets of simulated spectra shows that the AKARI infrared data alone can provide photometric redshifts accurate to |z_phot-z_spec|/(1+z_spec)<~10% (1-sigma) at z<~2. At higher redshifts the PAH features are shifted outside the wavelength range covered by AKARI and the photo-z estimates rely on the less prominent 1.6 micron stellar bump; the accuracy achievable in this case on (1+z) is ~10-15%, provided that the AGN contribution to the infrared emission is subdominant. Our technique is no more prone to redshift aliasing than optical-uv photo-z, and it may be possible to reduce this aliasing further with the addition of submillimetre and/or radio data.

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 a careful analysis of the point source detection limit of the AKARI All-Sky Survey in the WIDE-S 90 $\mu$m band near the North Ecliptic Pole (NEP). Timeline Analysis is used to detect IRAS sources and then a conversion factor is derived to transform the peak timeline signal to the interpolated 90 $\mu$m flux of a source. Combined with a robust noise measurement, the point source flux detection limit at S/N $>5$ for a single detector row is $1.1\pm0.1$ Jy which corresponds to a point source detection limit of the survey of $\sim$0.4 Jy. Wavelet transform offers a multiscale representation of the Time Series Data (TSD). We calculate the continuous wavelet transform of the TSD and then search for significant wavelet coefficients considered as potential source detections. To discriminate real sources from spurious or moving objects, only sources with confirmation are selected. In our multiscale analysis, IRAS sources selected above $4\sigma$ can be identified as the only real sources at the Point Source Scales. We also investigate the correlation between the non-IRAS sources detected in Timeline Analysis and cirrus emission using wavelet transform and contour plots of wavelet power spectrum. It is shown that the non-IRAS sources are most likely to be caused by excessive noise over a large range of spatial scales rather than real extended structures such as cirrus clouds.

We present a comparison between the SCUBA Half Degree Extragalactic Survey (SHADES) at 450 and 850 microns in the Lockman Hole East with a deep Spitzer Space Telescope survey at 3.6-24 microns conducted in Guaranteed Time. Using stacking analyses we demonstrate a striking correspondence between the galaxies contributing the submm extragalactic background light, with those likely to dominate the backgrounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photometric redshifts, we show that at least a third of the Spitzer-identified submm galaxies at 1<z<1.5 appear to reside in overdensities when the density field is smoothed at 0.5-2 Mpc comoving diameters, supporting the high-redshift reversal of the local star formation - galaxy density relation. We derive the dust-shrouded cosmic star formation history of galaxies as a function of assembled stellar masses. For model stellar masses <10^11 Msun, this peaks at lower redshifts than the ostensible z~2.2 maximum for submm point sources, adding to the growing consensus for ``downsizing'' in star formation. Our surveys are also consistent with ``downsizing'' in mass assembly. Both the mean star formation rates <dM/dt> and specific star formation rates <(1/M)dM/dt> are in striking disagreement with some semi-analytic predictions from the Millenium simulation. The discrepancy could either be resolved with a top-heavy initial mass function, or a significant component of the submm flux heated by the interstellar radiation field.

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

We present estimates of the photometric redshifts, stellar masses and star formation histories of sources in the SCUBA HAlf Degree Extragalactic Survey (SHADES). This paper describes the 60 SCUBA sources detected in the Lockman Hole covering an area of ~320 square arcmin. Using photometry spanning the B band to 8um, we find that the average SCUBA source forms a significant fraction of its stars in an early period of star formation and that most of the remainder forms in a shorter more intense burst around the redshift it is observed. This trend does not vary significantly with source redshift but the exact ratio of early to late mass is quite sensitive to the way extinction is treated in the modelling. However, the sources show a clear increase in stellar mass with redshift, consistent with downsizing. In terms of SED types, only two out of the 51 sources we have obtained photometric redshifts for are best fit by a quasar-like spectral energy distribution, with approximately 80 per cent of the sources being best fit with late-type spectra (Sc, Im and starburst). By including photometry at 850um, we conclude that the average SCUBA source is forming stars at a rate somewhere between 6 and 30 times the rate implied from the rest-frame optical in a dust obscured burst and that this burst creates 15-65 per cent of the total stellar mass. Using a simplistic calculation, we estimate from the average star formation history that between one in five and one in 15 bright (L_* +2 mag < L_optical < L_* -1 mag) galaxies in the field over the interval 0 < z < 3 will at some point in their lifetime experience a similar energetic dusty burst of star formation. Finally, we compute the evolution of the star formation rate density and find it peaks around z=2.

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 microns has been carried out to map the regions around several 850 micron-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared luminosities and hence star-formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350 micron flux densities for 24 850 micron sources, seven of which are detected at >2.5-sigma within a 10 arcsec search radius of the 850 micron positions. When results from the literature are included the total number of 350 micron flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-infrared (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich (Mdust~9x10^8 Msun), luminous (Lfir~2x10^12 Lsun), star-forming galaxies with intrinsic dust temperatures of ~35 K and star-formation rates of ~400 Msun/yr. We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350 micron source counts, N350(>25mJy)~200-500 deg^-2.

We present the results of optical identifications for 257 mid-infrared sources detected with a deep 15um survey over approximately 80 arcmin^2 area in the AKARI performance verification field near the North Ecliptic Pole. The 15um fluxes of the sources range from 1 mJy down to 40 uJy, approximately a half of which are below 100 uJy. Optical counterparts were searched for within a 2-3 arcsec radius in both the BVRi'z' catalog generated by using the deep Subaru/Suprime-cam field which covers one-third of the performance verification field, and the g'r'i'z' catalog based on observations made with MegaCam at CFHT. We found B-R and R-z' colours of sources with successful optical identifications are systematically redder than that of the entire optical sample in the same field. Moreover, approximately 40% of the 15um sources show colours R-L15>5, which cannot be explained by the spectral energy distribution (SED) of normal quiescent spiral galaxies, but are consistent with SEDs of redshifted (z>1) starburst or ultraluminous infrared galaxies. This result indicates that the fraction of the ultraluminous infrared galaxies in our faint 15um sample is much larger than that in our brighter 15um sources, which is consistent with the evolving mid-infrared luminosity function derived by recent studies based on the Spitzer 24um deep surveys. Based on an SED fitting technique, the nature of the faint 15um sources is further discussed for a selected number of sources with available K_s-band data.

The Far-Infrared Surveyor (FIS) is one of two focal plane instruments on the AKARI satellite. FIS has four photometric bands at 65, 90, 140, and 160 um, and uses two kinds of array detectors. The FIS arrays and optics are designed to sweep the sky with high spatial resolution and redundancy. The actual scan width is more than eight arcmin, and the pixel pitch is matches the diffraction limit of the telescope. Derived point spread functions (PSFs) from observations of asteroids are similar to the optical model. Significant excesses, however, are clearly seen around tails of the PSFs, whose contributions are about 30% of the total power. All FIS functions are operating well in orbit, and its performance meets the laboratory characterizations, except for the two longer wavelength bands, which are not performing as well as characterized. Furthermore, the FIS has a spectroscopic capability using a Fourier transform spectrometer (FTS). Because the FTS takes advantage of the optics and detectors of the photometer, it can simultaneously make a spectral map. This paper summarizes the in-flight technical and operational performance of the FIS.

We present the observations of the reflection nebulae IC4954 and IC4955 region with the Infrared Camera (IRC) and the Far-Infrared Surveyor (FIS) on board the infrared astronomical satellite AKARI during its performance verification phase. We obtained 7 band images from 7 to 160um with higher spatial resolution and higher sensitivities than previous observations. The mid-infrared color of the S9W (9um) and L18W (18um) bands shows a systematic variation around the exciting sources. The spatial variation in the mid-infrared color suggests that the star-formation in IC4954/4955 is progressing from south-west to north-east. The FIS data also clearly resolve two nebulae for the first time in the far-infrared. The FIS 4-band data from 65um to 160um allow us to correctly estimate the total infrared luminosity from the region, which is about one sixth of the energy emitted from the existing stellar sources. Five candidates for young stellar objects have been detected as point sources for the first time in the 11um image. They are located in the red S9W to L18W color regions, suggesting that current star-formation has been triggered by previous star-formation activities. A wide area map of the size of about 1 x 1 (deg^2) around the IC4954/4955 region was created from the AKARI mid-infrared all-sky survey data. Together with the HI 21cm data, it suggests a large hollow structure of a degree scale, on whose edge the IC4954/4955 region has been created, indicating star formation over three generations in largely different spatial scales.

AKARI, the first Japanese satellite dedicated to infrared astronomy, was launched on 2006 February 21, and started observations in May of the same year. AKARI has a 68.5 cm cooled telescope, together with two focal-plane instruments, which survey the sky in six wavelength bands from the mid- to far-infrared. The instruments also have the capability for imaging and spectroscopy in the wavelength range 2 - 180 micron in the pointed observation mode, occasionally inserted into the continuous survey operation. The in-orbit cryogen lifetime is expected to be one and a half years. The All-Sky Survey will cover more than 90 percent of the whole sky with higher spatial resolution and wider wavelength coverage than that of the previous IRAS all-sky survey. Point source catalogues of the All-Sky Survey will be released to the astronomical community. The pointed observations will be used for deep surveys of selected sky areas and systematic observations of important astronomical targets. These will become an additional future heritage of this mission.

We have studied the submillimetre (submm) properties of the following classes of near-infrared (NIR)-selected massive galaxies at high redshifts: BzK-selected star-forming galaxies (BzKs); distant red galaxies (DRGs); and extremely red objects (EROs). We used the SCUBA HAlf Degree Extragalactic Survey (SHADES), the largest uniform submm survey to date. Partial overlap of SIRIUS/NIR images and SHADES in SXDF has allowed us to identify 4 submm-bright NIR-selected galaxies, which are detected in the mid-infrared, 24 micron, and the radio, 1.4 GHz. We find that all of our submm-bright NIR-selected galaxies satisfy the BzK selection criteria, except for one galaxy whose B-z and z-K colours are however close to the BzK colour boundary. Two of the submm-bright NIR-selected galaxies satisfy all of the selection criteria we considered, i.e. they belong to the BzK-DRG-ERO overlapping population, or `extremely red' BzKs. Although these extremely red BzKs are rare (0.25 arcmin^{-2}), up to 20 % of this population could be submm galaxies. This fraction is significantly higher than that found for other galaxy populations studied here. Via a stacking analysis, we have detected the 850 micron flux of submm-faint BzKs and EROs in our SCUBA maps. While the contribution of z~2 BzKs to the submm background is about 10--15 % and similar to that from EROs typically at z~1, BzKs have a higher fraction (~30 %) of submm flux in resolved sources compared with EROs and submm sources as a whole. From the SED fitting analysis for both submm-bright and submm-faint BzKs, we found no clear signature that submm-bright BzKs are experiencing a specifically luminous evolutionary phase, compared with submm-faint BzKs. An alternative explanation might be that submm-bright BzKs are more massive than submm-faint ones.

We present the spatially resolved observations of IRAS sources from the Japanese infrared astronomy satellite AKARI All-Sky Survey during the performance verification (PV) phase of the mission. We extracted reliable point sources matched with IRAS point source catalogue. By comparing IRAS and AKARI fluxes, we found that the flux measurements of some IRAS sources could have been over or underestimated and affected by the local background rather than the global background. We also found possible candidates for new AKARI sources and confirmed that AKARI observations resolved IRAS sources into multiple sources. All-Sky Survey observations are expected to verify the accuracies of IRAS flux measurements and to find new extragalactic point sources.

Far-infrared (IR) images of the nearby Sb galaxy NGC2841 and the Sc galaxy NGC2976 at wavelengths of 65, 90, 140, and 160 um have been obtained with the Far-Infrared Surveyor (FIS) onboard AKARI. Both galaxies reveal similar morphologies of dust rings. They are, however, significantly different in the dust temperature: a cold (21 K) ring for NGC2841 and a warm (30 K) ring for NGC2976, which presumably reflects the difference in the origin of the ring structure for the different Hubble type of the galaxy. In addition to the dust ring structure, a warm dust component is detected from the central region of NGC2841, which may be attributed to the heating by its Low-Ionization Nuclear Emission-line Region nucleus. As for NGC2976, an extended dust component is observed along the minor axis, which shows a distribution somewhat asymmetrical to the galactic disk; this might be associated with the HI bridge in the M81/M82 group that NGC2976 belongs to. By taking advantage of a wealth of the far-IR bands of the FIS, it is demonstrated that the spectral energy distribution of NGC2841 is spatially more variable than that of NGC2976.

We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.

Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.

We present an initial analysis of AKARI 18um-selected galaxies using all 9 photometric bands at 2-24um available in the InfraRed Camera (IRC), in order to demonstrate new capabilities of AKARI cosmological surveys. We detected 72 sources at 18um in an area of 50.2 arcmin^2 in the AKARI/IRC monitor field towards the North Ecliptic Pole (NEP). From this sample, 25 galaxies with probable redshifts z>~ 0.5 are selected with a single colour cut (N2-N3>0.1) for a detailed SED analysis with ground-based BVRi'z'JK data. Using an SED radiative transfer model of starbursts covering the wavelength range UV -- submm, we derive photometric redshifts from the optical-MIR SEDs of 18um-selected galaxies. From the best-fit SED models, we show that the IRC all-band photometry is capable of tracing the steep rise in flux at the blue side of the PAH 6.2um emission feature. This indicates that the IRC all-band photometry is useful to constrain the redshift of infrared galaxies, specifically for dusty galaxies with a less prominent 4000A break. Also, we find that the flux dip between the PAH 7.7 and 11.2um emission feature is recognizable in the observed SEDs of galaxies at z~1. By using such a colour anomaly due to the PAH and silicate absorption features, unique samples of ULIRGs at z~1, `silicate-break' galaxies, can be constructed from large cosmological surveys of AKARI towards the NEP, i.e. the NEP-Deep and NEP-Wide survey. This pilot study suggests the possibility of detecting many interesting galaxy properties in the NEP-Deep and Wide surveys, such as a systematic difference in SEDs between high- and low-z ULIRGs, and a large variation of the PAH inter-band strength ratio in galaxies at high redshifts. [abridged]

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

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

We present a result of cross-correlating the Infrared Astronomical Satellite Faint Source Catalogue (IRAS FSC) with the spectroscopic catalogues of galaxies in the Fourth Data Release of Sloan Digital Sky Survey (SDSS), the Final Data Release of 2dF Galaxy Redshift Survey (2dFGRS) and the Second Data Release of 6dF Galaxy Survey (6dFGS). We have identified 324 ultraluminous infrared galaxies (ULIRGs) including 190 newly discovered ULIRGs, and 2 hyperluminous infrared galaxies (HLIRGs). Adding these new ULIRGs, we increase the number of known ULIRGs by about 30 per cent. The reliability of the cross-correlation is estimated using the likelihood ratio method. The incompleteness of our sample introduced by the identification procedure in this study is estimated to be about 5 per cent. Our sample covers the redshift range of z=0.037$-$0.517 with a median redshift of $\bar{z}$=0.223, which is larger than that ($\bar{z}$=0.184) of the sample of previously known ULIRGs.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

Taken together, the blank-field sub-mm surveys undertaken with SCUBA prior to the SHADES survey cover a total area of 460 sq. arcmin to a range of depths. However combining the results from these surveys has hitherto been complicated by the fact that different survey groups have used different methods of data reduction and source extraction. Here we present the results of re-reducing and analysing all of these blank field SCUBA in an almost identical manner to that employed in the SCUBA 8-mJy Survey. Comparative source catalogues are given which include a number of new significant source detections as well as failing to confirm some of those objects previously published. These new source catalogues have been combined to produce the most accurate number counts to date from 2 to 12.5 mJy. The cumulative number counts appear to steepen beyond S(850)=8mJy, which could indicate an intrinsic turn-over in the underlying luminosity function placing an upper limit on the luminosity (and hence mass) of a high redshift galaxy. We have also investigated the clustering properties of the bright SCUBA population by means of 2-point angular correlation functions. We find evidence for strong clustering broadly consistent with that measured for Extremely Red Objects (EROs). Nearest-neighbour analyses further support strong clustering on arcmin scales, rejecting the null hypothesis that the distribution of the bright submm sources is random at the 99% confidence level.

AKARI (formerly ASTRO-F) is an infrared space telescope designed for an all-sky survey at 10-180 (mu)m, and deep pointed surveys of selected areas at 2-180 (mu)m. The deep pointed surveys with AKARI will significantly advance our understanding of galaxy evolution, the structure formation of the Universe, the nature of the buried AGNs, and the cosmic infrared background. Here we describe the important characteristics of the AKARI mission: the orbit, and the attitude control system, and investigate the optimum survey area based on the updated pre-flight sensitivities of AKARI, taking into account the cirrus confusion noise as well as the surface density of bright stars. The North Ecliptic Pole (NEP) is concluded to be the best area for 2-26 (mu)m deep surveys, while the low-cirrus noise regions around the South Ecliptic Pole (SEP) are worth considering for 50-180 (mu)m pointed surveys to high sensitivities limited by the galaxy confusion noise. Current observational plans of these pointed surveys are described in detail. Comparing these surveys with the deep surveys with the Spitzer Space Telescope, the AKARI deep surveys are particularly unique in respect of their continuous wavelength coverage over the 2-26 (mu)m range in broad-band deep imaging, and their slitless spectroscopy mode over the same wavelength range.

The modest significance of most sources detected in current submm surveys can potentially compromise some analyses due to the inclusion of spurious sources in catalogues typically selected at >3.0-3.5sigma. Here, we develop and apply a dual-survey extraction technique to SCUBA and MAMBO images of the Lockman Hole. Cut above 5sigma, our catalogue of SMGs is more robust than previous samples, with a reduced likelihood of real, but faint SMGs (beneath and around the confusion limit) entering via superposition with noise. The effective flux limit of the survey is well matched to our deep 1.4-GHz image. The former is sensitive to luminous, dusty galaxies at extreme redshifts whilst the latter probes the z<3 regime. A high fraction of our robust SMGs (~80%) have radio counterparts which, given the ~10% contamination by spurious sources, suggests that very distant SMGs (z>>3) are unlikely to make up more than ~10% of the bright SMG population. This implies that almost all of the bright SMG population is amenable to study via the deepest current radio imaging. We use these radio counterparts to provide an empirical calibration of the positional uncertainty in SMG catalogues. We then go on to outline the acquisition of redshifts for radio-identified SMGs, from sample selection in the submm, to counterpart selection in the radio and optical/IR, to slit placement on spectrograph masks. We determine a median of z=2.05 +/- 0.41 from a sample of six secure redshifts for unambigious radio-identified submm sources and z=2.14 +/- 0.27 when we include submm sources with multiple radio counterparts and/or less reliable redshifts. These figures are consistent with previous estimates, suggesting that our knowledge of the median redshift of bright SMGs population has not been biased by the low significance of the source catalogues employed.

The SCUBA HAlf Degree Extragalactic Survey (SHADES) is a major blank-field extragalactic sub-mm survey underway at the James Clerk Maxwell telescope.SHADES aims to cover half a square degree at 450+850microns to a 4sigma depth of ~8mJy at 850microns. Two fields are being observed, the Subaru/XMM--Newton Deep Field (SXDF) and the Lockman Hole East. The survey has 3 main aims to investigate: i) the cosmic history of massive dust enshrouded star-formation activity, ii) the clustering properties of sub-mm-selected galaxies in order to determine whether these objects could be progenitors of present-day massive ellipticals, and iii) the fraction of sub-mm-selected sources that harbour AGN. To achieve these aims requires the sub-mm data be combined with co-spatial information spanning the radio-to-X-ray frequency range, the resulting extensive multi-wavelength dataset providing complete photometric redshift information accurate to delta z <~0.5 and detailed SEDs for the vast majority of the sub-mm-selected sources. In this paper, the first of a series on SHADES, we present an overview of the motivation for the survey, describe the survey strategy and provide a detailed description of the primary data analysis pipeline. As of Feb. 2004, 720arcmin^2 had been mapped with SCUBA (about 40% of the anticipated final total) to a median 1sigma depth of 2.2mJy per beam at 850microns (25mJy per beam at 450microns), and a source density of 650+/-50 sources/deg^2 >3sigma at 850microns, uncorrected for Eddington bias. A refined re-analysis of the original 8-mJy survey Lockman hole data was carried out in order to evaluate the new data reduction pipeline. Of the 17 most secure sources in the original sample, 12 have been re-confirmed, including 10 of the 11 for which radio IDs were previously secured.(Abridged)

We report on a search for the optical counterparts of 175 um - selected sources from the Far-Infrared Background (FIRBACK) survey in the European Large Area ISO Survey (ELAIS) N2 field. Applying a likelihood ratio technique to optical catalogues from the Isaac Newton Telescope - Wide Field Survey (INT--WFS), we found optical identifications for 33 out of 55 FIRBACK sources in this field. These were then reassessed in the light of associations with the ELAIS final catalogue for the the N2 field, to yield a final set of 31 associations. We have investigated the nature of this population through a comparison of their observed spectral energy distributions with predictions from radiative transfer models which simulate the emission from both cirrus and starburst components. We find the far-infrared sources to be 80 per cent star bursting galaxies with their starburst component at a high optical depth. The resulting SEDs were used to estimate far-infrared luminosities, star formation rates, dust temperatures and dust masses. The N2 FIRBACK population is found to consist of four suspected ULIRGs, a number of LIRGs and a population of low redshift quiescently star forming galaxies. We also discuss the implications of these results for current evolutionary models.

We present the Final Analysis of the European Large Area ISO Survey (ELAIS) 15 $\mu$m observations, carried out with the ISOCAM instrument on board the Infrared Space Observatory (ISO). The data reduction method, known as LARI method, is based on a mathematical model of the detector's behaviour and was specifically designed for the detection of faint sources in ISO-CAM/PHOT data. The method is fully interactive and leads to very reliable and complete source lists. The resulting catalogue includes 1923 sources detected with $S/N > 5$ in the {0.5 -- 100 mJy} flux range and over an area of 10.85 {deg$^2$} split into four fields, making it the largest non-serendipitous extragalactic source catalogue obtained to date from ISO data. This paper presents the concepts underlying the data reduction method together with its latest enhancements. The data reduction process, the production and basic properties of the resulting catalogue are discussed. The catalogue quality is assessed by means of detailed simulations, optical identifications and comparison with previous analyses.

We discuss optical associations, spectral energy distributions and photometric redshifts for SWIRE sources in the ELAIS-N1 area and the Lockman Validation Field. The band-merged IRAC (3.6, 4.5, 5.8 and 8.0 mu) and MIPS (24, 70, 160 mu) data have been associated with optical UgriZ data from the INT Wide Field Survey in ELAIS-N1, and with our own optical Ugri data in Lockman-VF. The spectral energy distributions of selected ELAIS sources in N1 detected by SWIRE, most with spectroscopic redshifts, 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 7%. The photometric redshifts are used to derive the 3.6 and 24 mu redshift distribution and to compare this with the predictions of models. For those sources with a clear mid infrared excess, relative to the galaxy starlight model used for the optical and near infrared, the mid and far infrared data are modelled in terms of the same dust emission templates. The proportions found of each template type are: cirrus 31%, M82 29%, Arp 220 10%, AGN dust tori 29%. 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.

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.

We present a re--analysis of the European Large Area ISO Survey (ELAIS) 90mum observations carried out with ISOPHOT, an instrument on board the ESA's Infrared Space Observatory (ISO). With more than 12 sq. deg., the ELAIS survey is the largest area covered by ISO in a single program and is about one order of magnitude deeper than the IRAS 100mum survey. The data analysis is presented and was mainly performed with the Phot Interactive Analysis software but using the pairwise method of Stickel et al. (2003) for signal processing from ERD (Edited Raw Data) to SCP (Signal per Chopper Plateau). The ELAIS 90mum catalogue contains 229 reliable sources with fluxes larger than 70 mJy and is available at www.blackwell-synergy.com. Number counts are presented and show an excess above the no-evolution model prediction. This confirms the strong evolution detected at shorter(15mum) and longer (170mum) wavelengths in other ISO surveys. The ELAIS counts are in agreement with previous works at 90mum and in particular with the deeper counts extracted from the Lockman hole observations. Comparison with recent evolutionary models show that the models of Franceschini et al. and Guiderdoni et al. which includes a heavily-extinguished population of galaxies give the best fit to the data. Deeper observations are nevertheless required to better discriminate between the model predictions in the far-infrared and are scheduled with the Spitzer Space Telescope which already started operating and will also be performed by ASTRO-F.

We present a cross-correlation of the European Large Area ISO survey (ELAIS) with the ELAIS Deep X-ray Survey of the N1 and N2 fields. There are 7 Chandra point sources with matches in the ELAIS Final Analysis 15um catalogue, out of a total of 28 extragalactic ISO sources present in the Chandra fields. Five of these are consistent with AGN giving an AGN fraction of ~19 per cent in the 15um flux range 0.8-6 mJy. We have co-added the hard X-ray fluxes of the individually-undetected ISO sources and find a low significance detection consistent with star formation in the remaining population. We combine our point source cross-correlation fraction with the XMM-Newton observations of the Lockman Hole and Chandra observations of the Hubble Deep Field North to constrain source count models of the mid-infrared galaxy population. The low dust-enshrouded AGN fraction in ELAIS implied by the number of cross-identifications between the ELAIS mid-infrared sample and the Chandra point sources is encouraging for the use of mid-infrared surveys to constrain the cosmic star formation history, provided there are not further large undetected populations of Compton-thick AGN.

We present new determinations of the local sub-mm luminosity functions, solving the ``sub-mm Olbers' Paradox.'' We also present predictions of source counts and luminosity functions in current and future far-infrared to sub-mm surveys. Using the sub-mm colour temperature relations from the SCUBA Local Universe Galaxy Survey, and the discovery of excess 450 micron excess emission in these galaxies, we interpolate and extrapolate the IRAS detections to make predictions of the SEDs of all 15411 PSC-z galaxies from 50-1300 microns. Despite the long extrapolations we find excellent agreement with (a) the 90 micron luminosity function of Serjeant et al. (2001), (b) the 850 micron luminosity function of Dunne et al. (2000), (c) the mm-wave photometry of Andreani & Franceschini (1996); (d) the asymptotic differential and integral source count predictions at 50-1300 microns by Rowan-Robinson (2001). We find the local 850 micron sub-mm luminosity density converges to (7.3+/-0.2)x10^{19} h_{65} W/Hz/Mpc^3. Remarkably, the local spectral luminosity density and the extragalactic background light together strongly constrain the cosmic star formation history for a wide class of evolutionary assumptions. We find that the extragalactic background light, the 850 micron 8mJy source counts, and the Omega_* constraints all independently point to a decline in the comoving star formation rate at z>1. In order to reconcile this with direct determinations, we suggest either there is a top-heavy initial mass function at high redshifts, and/or there is stronger evolution in the more luminous far-infrared galaxies than seen in the population as a whole.

We present the 90 micron luminosity function of the Final Analysis of the European Large Area ISO Survey (ELAIS), extending the sample size of our previous analysis (paper IV) by about a factor of 4. Our sample extends to z=1.1, around 50 times the comoving volume of paper IV, and 10^{7.7} < h^{-2}L/Lsun < 10^{12.5}. From our optical spectroscopy campaigns of the northern ELAIS 90 mircon survey (7.4 deg^2 in total, to S(90um)>70mJy), we obtained redshifts for 61% of the sample (151 redshifts) to B<21 identified at 7 microns, 15 microns, 20cm or with bright (B<18.5) optical identifications. The selection function is well-defined, permitting the construction of the 90 micron luminosity function of the Final Analysis catalogue in the ELAIS northern fields, which is in excellent agreement with our Preliminary Analysis luminosity function in the ELAIS S1 field from paper IV. The luminosity function is also in good agreement with the IRAS-based prediction of Serjeant & Harrison (2004), which if correct requires luminosity evolution of (1+z)^{3.4 +/- 1.0} for consistency with the source counts. This implies an evolution in comoving volume averaged star formation rate at z<~1 consistent with that derived from rest-frame optical and ultraviolet surveys.

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.

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

We present Spitzer observations in five wavebands between 3.6 and 24um of an unbiased sample of 9 luminous, dusty galaxies selected at 1200um by the MAMBO camera on the IRAM 30-m telescope, a population akin to the well-known submm or `SCUBA' galaxies (hereafter SMGs). Owing to the coarse resolution of submm/mm instrumentation, SMGs have traditionally been difficult to identify at other wavelengths. We compare our multi-wavelength catalogs to show that the overlap between 24 and 1200um must be close to complete at these flux levels. We find that all (4/4) of the most secure >=4sigma SMGs have robust >=4sigma counterparts at 1.4GHz, while the fraction drops to 7/9 using all >=3sigma SMGs. We show that combining mid-IR and marginal (>=3sigma) radio detections provides plausible identifications in the remaining cases, enabling us to identify the complete sample. Accretion onto an obscured central engine is betrayed by the shape of the mid-IR continuum emission for several sources, confirming Spitzer's potential to weed out active galaxies. We demonstrate the power of a S(24um)/S(8um) vs S(8um)/S(4.5um) color-color plot as a diagnostic for this purpose. However, we conclude that the majority (~75%) of SMGs have rest-frame mid-/far-IR SEDs commensurate with obscured starbursts. Sensitive 24-um observations are clearly a useful route to identify and characterize reliable counterparts to high-redshift far-IR-bright galaxies, complementing what is possible via deep radio imaging.

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

We present the final band-merged ELAIS catalogue at 6.7, 15, 90, and 175 $\mu$m, and the associated data at u,g,r,i,z,J,H,K, and 20cm. The origin of the survey, infrared and radio observations, data-reduction and optical identifications are briefly reviewed, and a summary of the area covered, and completeness limit for each infrared band is given. A detailed discussion of the band-merging and optical association strategy is given. The total catalogues consists of 2860 sources. For extragalactic sources observed in 3 or more infrared bands, colour-colour diagrams are presented and discussed in terms of the contributing infrared populations. Spectral energy distributions are shown for selected sources and compared with cirrus, M82 and Arp220 starburst, and AGN dust torus models. Spectroscopic redshifts are tabulated, where available. For the N1 and N2 areas, the INT ugriz Wide Field Survey, permits photometric redshifts to be estimated for galaxies and quasars. These agree well with the spectroscopic redshifts, within the uncertainty of the photometric method. The redshift distribution is given for selected ELAIS bands and colour-redshift diagrams are discussed. There is a high proportion of ultraluminous infrared galaxies in the ELAIS Catalogue ($> 10 %$ of 15 $\mu$m sources), many with Arp220-like colours. 10 hyperluminous infrared galaxies and 10 EROs are found in the survey. The large numbers of ultraluminous galaxies imply very strong evolution in the star-formation rate between z = 0 and 1.

We present a multi-wavelength catalog (15 um, R, K-band, 1.4 GHz flux) plus spectroscopic identifications for 406 15 um sources detected in the ELAIS region S1, over the flux density range 0.5<S<150 mJy. 332 (~82%) sources are optically identified down to R~23.0. Spectra or bona fide stellar identifications are obtained for 290 objects. The areal coverage, MIR and optical completeness of the sample are discussed in order to allow statistical and evolutionary analyses. Two main spectroscopic classes have been found to dominate the MIR extragalactic population: z<0.5 star-forming galaxies (from absorbed to extreme starbursts: nuL_nu~10^8-10^11 L_odot), which account for \~75% of the sources, and AGN (both type 1 and 2), which account for ~25% of the sources. About 20% of the extragalactic sources are dust-enshrouded starburst galaxies [e(a) spectra], and all the starburst galaxies appear more dust extincted in the optical than nearby normal galaxies. We also identified 91 stellar objects (~22% of the MIR sources). The counts for starburst galaxies and AGN down to 0.6 mJy have been derived. A general trend is found in the optical-MIR SED of the galaxies, where the MIR-luminous objects have larger MIR to optical luminosity ratios.

Hard X-ray surveys have uncovered a large population of heavily obscured AGN. They also reveal a population of quasars with moderate obscuration at both visible and X-ray wavelengths. We use Chandra selected samples of quasars from the ELAIS Deep X-ray Survey (EDXS) and CDF-N to investigate the obscuration towards the nuclei of moderately obscured AGN. We find an inverse correlation between the optical to X-ray flux ratio and the X-ray hardness ratio which can be interpreted as due to obscuration at visible and X-ray wavelengths. We present detailed optical and near-infrared data for a sample of optically-faint (R>23) quasars from the EDXS. These are used to constrain the amount of rest-frame UV/optical reddening towards these quasars. It is found that optically-faint quasars are mostly faint due to obscuration, not because they are intrinsically weak. After correcting for reddening, the optical magnitudes of most of these quasars are similar to the brighter quasars at these X-ray fluxes. Combining with gas column densities inferred from the X-ray observations we consider the gas-to-dust ratios of the obscuring matter. We find that the quasars generally have higher gas-to-dust absorption than that seen in the Milky Way - similar to what has been found for nearby Seyfert galaxies. We consider the possible existence of a large population of X-ray sources which have optical properties of Type 1 (unobscured) quasars, but X-ray properties of Type 2 (obscured) quasars. We show that the observed distribution of optical-to-X-ray flux ratios of quasars at z>1 is skewed to low values compared to the intrinsic distribution...

We present the optical identification of mid-IR and radio sources detected in the European Large Area ISO Survey (ELAIS) areas N1 and N2. Using the r' band optical data from the Wide Field Survey we apply a likelihood ratio method to search for the counterparts of the 1056 and 691 sources detected at 15 micron and 1.4 GHz respectively, down to flux limits of S_{15}=0.5 mJy and S_{1.4 GHz}=0.135 mJy. We find that ~92% of the 15 micron ELAIS sources have an optical counterpart down to the magnitude limit of the optical data, r'=24. All mid-IR sources with fluxes S_{15} >= 3 mJy have an optical counterpart. The magnitude distribution of the sources shows a well defined peak at relatively bright magnitudes r'~18. About 20% of the identified sources show a point-like morphology; its magnitude distribution has a peak at fainter magnitudes than those of galaxies. The mid-IR-to-optical and radio-to-optical flux diagrams are presented and discussed in terms of actual galaxy models. Objects with mid-IR-to-optical fluxes larger than 1000 are found that can only be explained as highly obscured star forming galaxies or AGNs. Blank fields being 8% of the 15 micron sample have even larger ratios suggesting that they may be associated with higher redshift and higher obscured objects.

We present an analysis of the two-point angular correlation function of the ELAIS S1 survey. The survey covers 4 deg$^2$ and contains 462 sources detected at 15$\mu$m to a 5$\sigma$ flux limit of 0.45 mJy. Using the 329 extragalactic sources not repeated in different observations, we detect a significant clustering signal; the resulting angular correlation function can be fitted by a exponential law $w(\theta) = A \theta^{1-\gamma}$ with $A = 0.014\pm0.005$ and $\gamma = 2.04\pm0.18$. Assuming a redshift distribution of the objects, we invert Limber's equation and deduce a spatial correlation length $r_0=4.3^{+0.4}_{-0.7} h^{-1}$ Mpc. This is smaller than that obtained from optical surveys but it is in agreement with results from IRAS. This extends to higher redshift the observational evidence that infrared selected surveys show smaller correlation lengths (i.e. reduced clustering amplitudes) than optical surveys.

We present deep wide-field K-band observations of six Hyperluminous Infrared Galaxies (HLIRGs) spanning a redshift range 0.44<z<1.55. The sample resides in a wide variety of environments, from the field to Abell 2 clusters, with a mean galaxy-HLIRG clustering amplitude of B_gh=190+/-45. The range in environments, and the mean clustering level, are both greater than those seen in local IR-luminous galaxies, from which we infer that the range of galaxy evolution processes driving IR-luminous galaxy evolution at z>0.5 is greater than locally, and includes mergers between gas-rich spiral galaxies in the field, but also includes encounters in clusters and hierarchical buildup. The similarity in the range of environments and mean clustering amplitude between our sample and QSOs over a similar redshift range is consistent with the interpretation where evolutionary connections between IR-luminous galaxies and QSOs are stronger at z>0.5 than locally, and that, at these redshifts, the processes that drive QSO evolution are similar to those that drive IR-luminous galaxy evolution. From comparison of the HLIRG and QSO host galaxies we further postulate that a larger fraction of IR-luminous galaxies pass through an optical QSO stage at z>0.5 than locally.

We present the K-band Hubble diagrams (K-z relations) of sub-mm-selected galaxies and hyperluminous galaxies (HLIRGs). We report the discovery of a remarkably tight K-z relation of HLIRGs, indistinguishable from that of the most luminous radiogalaxies. Like radiogalaxies, the HLIRG K-z relation at z<~3 is consistent with a passively evolving ~3L* instantaneous starburst starting from a redshift of z~10. In contrast, many sub-mm selected galaxies are >~2 magnitudes fainter, and the population has a much larger dispersion. We argue that dust obscuration and/or a larger mass range may be responsible for this scatter. The galaxies so far proved to be hyperluminous may have been biased towards higher AGN bolometric contributions than sub-mm-selected galaxies due to the 60 micron selection of some, so the location on the K-z relation may be related to the presence of the most massive AGN. Alternatively, a particular host galaxy mass range may be responsible for both extreme star formation and the most massive active nuclei.

We present new determinations of the local sub-mm luminosity functions. We find the local sub-mm luminosity density converging to 7.3+/-0.2 x 10^19 W/Hz/Mpc^3 /h_65 at 850um solving the ``sub-mm Olbers' Paradox.'' Using the sub-mm colour temperature relations from the SCUBA Local Universe Galaxy Survey, and the discovery of excess 450um excess emission in these galaxies, we interpolate and extrapolate the IRAS detections to make predictions of the SEDs of all 15411 PSC-z galaxies from 50-3000um. Despite the long extrapolations we find excellent agreement with (a) the 90um luminosity function of Serjeant et al. (2001), (b) the 850um luminosity function of Dunne et al. (2000), (c) the mm-wave photometry of Andreani & Franceschini (1996); (d) the asymptotic differential and integral source count predictions at 50-3000um by Rowan-Robinson (2001). Remarkably, the local luminosity density and the extragalactic background light together strongly constrain the cosmic star formation history for a wide class of evolutionary assumptions. We find that the extragalactic background light, the 850um 8mJy source counts, and the Omega_* constraints all independently point to a decline in the comoving star formation rate at z>1.

We present an extended analysis of the SCUBA observations of the Hubble Deep Field (HDF), expanding the areal coverage of the Hughes et al. 1998 study by a factor of ~1.8 and containing at least three further sources in addition to the five in that study. We also announce the public release of the reduced data products. The map is the deepest ever made in the sub-millimetre, obtained in excellent conditions (median 850um optical depth of 0.16). Two independent reductions were made, one with SURF and the other with a wholly algorithmic IDL analysis which we present in detail here. Of the three new sources, all appear to be at z>~0.9 and one is provisionally associated with an Extremely Red Object (I-K>5). There appears to be no significant cross-correlation signal between the 850um fluctuations and sources detected by ISOCAM, the VLA or Chandra, nor with Very Red Objects (I-K>4), nor quasars and quasar candidates in the HDF (notwithstanding a small number of individual weak candidate detections). This is consistent with interpretations where the 850um-selected galaxies are at higher redshifts than those currently probed by ISOCAM/VLA, and predominantly not Compton-thin AGN. There are only one or two compelling cases for the radio source being the sub-mm source. Nevertheless, most SCUBA-HDF point sources have a nearby radio source apparently well-separated from the sub-mm centroid.

The VLTI/GENIE configuration will operate using at least 4 of the VLTI telescopes (and possibly with one or more of the AT telescopes in the future if adaptive optics become available on them). GENIE effectively can be thought of as a 'smart' coronagraph, enabling high dynamic range imaging to be achieved at moderate spatial resolution, with high rejection of the emission of a central bright point source. However, but this bright source rejection may only provide a rather moderate image quality (due to the few baselines and transfer function on the sky). Operated in this way, only limited image reconstruction is possible since classical radio and millimeter wavelength interferometry techniques are not directly applicable to the outputs of optical interferometers because the absolute phases are generally not measured. However, measurements of visibility and closure phase could lead to situations where image reconstruction becomes possible. This paper addresses the issue of whether there are areas outside of the exoplanet search where it might be able to make a useful impact on astronomy.

We present an analysis of two deep (75 ks) Chandra observations of the European Large Area ISO Survey (ELAIS) fields N1 and N2 as the first results from the ELAIS deep X-ray survey. This survey is being conducted in well studied regions with extensive multi-wavelength coverage. Here we present the Chandra source catalogues along with an analysis of source counts, hardness ratios and optical classifications. A total of 233 X-ray point sources are detected in addition to 2 soft extended sources, which are found to be associated with galaxy clusters. An over-density of sources is found in N1 with 30% more sources than N2, which we attribute to large-scale structure. A similar variance is seen between other deep Chandra surveys. The source count statistics reveal an increasing fraction of hard sources at fainter fluxes. The number of galaxy-like counterparts also increases dramatically towards fainter fluxes, consistent with the emergence of a large population of obscured sources.

We present maps, source lists, and number counts from the largest, unbiassed, extragalactic sub-mm survey so far undertaken with the SCUBA camera on the JCMT. Our maps cover 260 sq. arcmin, to a noise level S(850)=2.5 mJy/beam. We have reduced the data using both SURF, and our own pipeline which produces zero-footprint maps and noise images. The uncorrelated noise maps produced by the latter approach have allowed application of a maximum-likelihood method to measure the statistical significance of each peak, leading to properly quantified flux-density errors for all potential sources. We detect 19 sources with S/N > 4, 38 with S/N > 3.5, and 72 with S/N > 3. To assess completeness and the impact of source confusion we have applied our source extraction algorithm to a series of simulated images. The result is a new estimate of the sub-mm source counts in the flux-density range S(850)=5-15mJy, which we compare with other estimates, and with model predictions. Our estimate of the cumulative source count at S(850) > 8 mJy is 320 (+80,-100) per square degree. Assuming that the majority of sources have z > 1.5, the co-moving number density of high-z galaxies forming stars at a rate >1000 solar masses per year is 10^-5 per Mpc^3, with only a weak dependence on the precise redshift distribution. This number density corresponds to that of massive ellipticals with L > 3-4 L* at low redshift and is also the same as the co-moving number density of comparably massive, passively-evolving objects in the redshift band 1<z<2 inferred from recent surveys of extremely red objects. Thus the bright sub-mm sources can plausibly account for the formation of all present-day massive ellipticals. Improved z constraints, and a proper measurement of sub-mm clustering can refine or refute this picture.

The sources discovered in deep hard X-ray surveys with 2-8 keV fluxes of 10^-14 erg cm^-2 s^-1 make up the bulk of the X-ray background at these energies. We present here detailed multi-wavelength observations of three such sources from the ELAIS Deep X-ray Survey. The observations include sensitive near-infrared spectroscopy with the Subaru Telescope and X-ray spectral information from the Chandra X-ray Observatory. The sources observed all have optical-to-near-IR colours redder than an unobscured quasar and comprise a reddened quasar, a radio galaxy and an optically-obscured AGN. The reddened quasar is at a redshift z=2.61 and shows a very large X-ray absorbing column of N_H approx 3.10^23 cm^-2. This contrasts with the relatively small amount of dust reddening, implying a gas-to-dust ratio along the line-of-sight a hundred times greater than that of the Milky Way. The radio galaxy at z=1.57 shows only narrow emission lines, but has a surprisingly soft X-ray spectrum. The softness of this spectrum either indicates an unusually low gas-to-dust ratio for the absorbing medium or X-ray emission related to the young radio source. The host galaxy is extremely red (R-K=6.4) and its optical/near-IR spectrum is best fit by a strongly reddened (A_V~2) starburst. The third X-ray source discussed is also extremely red (R-K=6.1) and lies in a close grouping of three other R-K>6 galaxies. No emission or absorption lines were detected from this object, but its redshift (and that of one of the nearby galaxies) are constrained by SED-fitting to be just greater than z=1. The extremely red colours of these two galaxies can be accounted for by old stellar populations. These observations illustrate the diverse properties of hard X-ray selected AGN.

We explore the relationship between the hard X-ray and sub-mm populations using deep Chandra observations of a large, contiguous SCUBA survey. In agreement with other recent findings, we confirm that the direct overlap is small. Of the 17 sub-mm sources detected in this field at 850 microns, only one is coincident with a Chandra source. The resulting limits imply that the majority of SCUBA sources are not powered by AGN, unless the central engine is obscured by Compton-thick material with a low (<1 per cent) scattered component. Furthermore, since Chandra detects only ~5 per cent of SCUBA sources, the typical obscuration would need to be almost isotropic. The X-ray upper limits are so strong that in most cases we can also rule out a starburst SED at low redshift, suggesting that the majority of SCUBA sources lie at z>1 even if they are purely starburst galaxies. Despite the low detection rate, we find evidence for strong angular clustering between the X-ray and sub-mm populations. The implication is that AGN and SCUBA sources trace the same large-scale structure but do not generally coincide. If bright sub-mm sources represent massive elliptical galaxies in formation, we suggest that (for a given galaxy) the major episode of star-formation must be distinct from the period of observable quasar activity.

We present J and K-band near-infrared photometry of a sample of mid-infrared sources detected by the Infrared Space Observatory (ISO) as part of the European Large Area ISO-Survey (ELAIS) and study their classification and star-forming properties. We have used the Preliminary ELAIS Catalogue for the 6.7 micron (LW2) and 15 micron (LW3) fluxes. All of the high-reliability LW2 sources and 80 per cent of the LW3 sources are identified in the near-IR survey reaching K = 17.5 mag. The near- to mid-IR flux ratios can effectively be used to separate stars from galaxies in mid-IR surveys. The stars detected in our survey region are used to derive a new accurate calibration for the ELAIS ISOCAM data in both the LW2 and LW3 filters. We show that near to mid-IR colour-colour diagrams can be used to further classify galaxies, as well as study star-formation. The ISOCAM ELAIS survey is found to mostly detect strongly star-forming late-type galaxies, possibly starburst powered galaxies, and it also picks out obscured AGN. The ELAIS galaxies yield an average mid-IR flux ratio LW2/LW3 = 0.67 +/- 0.27. We discuss this [6.7/15] ratio as a star formation tracer using ISO and IRAS data of a local comparison sample. We find that the [K/15] ratio is also agood indicator of activity level in galaxies and conclude that the drop in the [6.7/15] ratio seen in strongly star-forming galaxies is a result of both an increase of 15 mic emission and an apparent depletion of 6.7 mic emission. Near-IR data together with the mid-IR give the possibility to estimate the relative amount of interstellar matter in the galaxies.

The SCUBA 8-mJy survey is the largest submm extragalactic mapping survey undertaken to date, centred on the Lockman Hole and ELAIS N2 regions. Here, we present new 1.4-GHz imaging of the depth and resolution necessary to reliably identify radio counterparts for 18 of 30 submm sources. Armed with this greatly improved positional information, we present and analyse new optical, near-IR and XMM-Newton X-ray imaging to identify optical/IR host galaxies to half of the submm-selected sources in those fields. As many as 15% of the submm sources detected at 1.4 GHz are resolved by the 1.4'' beam and a further 25% have more than one radio counterpart, suggesting that radio and submm emission arise from extended starbursts and that interactions are common. We note that less than a quarter of the submm-selected sample would have been recovered by targeting optically faint radio sources, underlining the selective nature of such surveys. At least 60% of the radio-confirmed optical/IR host galaxies appear to be morphologically distorted; many are composite systems - red galaxies with relatively blue companions; just over half are found to be very red or extremely red; contrary to popular belief, most are sufficiently bright to be tackled with spectrographs on 8-m telescopes. We estimate the median redshift of the S(850um) ~ 8 mJy submm galaxy population: if the radio/far-IR correlation holds at high redshift, and our sample is unbiased, we derive a conservative median z >= 2.0, or >= 2.4 using spectral templates more representative of known submm galaxies [abridged].

Despite extensive observational efforts, the brightest sub--mm source in the Hubble Deep Field, HDF850.1, has failed to yield a convincing optical/infrared identification almost 4 years after its discovery. This failure is all the more notable given the availability of supporting multi-frequency data of unparalleled depth, and sub-arcsec positional accuracy for the sub-mm/mm source. Consequently, HDF850.1 has become a test case of the possibility that the most violently star-forming objects in the universe are too red and/or distant to be seen in the deepest optical images. Here we report the discovery of the host galaxy of HDF850.1. This object has been revealed by careful analysis of a new, deep K-prime image of the HDF obtained with the Subaru 8.2-m telescope. Its reality is confirmed by a similar analysis of the HST NICMOS F160W image of the same region. This object is extremely faint (K=23.5), clumpy (on sub-arcsec scales) and very red (I-K > 5.2; H-K = 1.4 +- 0.35). The likelihood that it is the correct identification is strongly reinforced by a reanalysis of the combined MERLIN+VLA 1.4-GHz map of the field which yields a new radio detection of HDF850.1 only 0.1 arcsec from the new near-ir counterpart, and with sufficient positional accuracy to exclude all previously considered alternative optical candidates. We have calculated new confidence limits on the estimated redshift of HDF850.1 and find z = 4.1 +- 0.5. We also calculate that the flux density of HDF850.1 has been boosted by a factor of ~3 through lensing by the intervening elliptical 3-586.0, consistent with predictions that a small but significant fraction of blank-field sub-mm sources are lensed by foreground galaxies. We discuss the wider implications of these results for the sub-mm population and cosmic star-formation history.