We present the discovery with Keck of the extremely infrared (IR) luminous transient AT 2017gbl, coincident with the Northern nucleus of the luminous infrared galaxy (LIRG) IRAS 23436+5257. Our extensive multi-wavelength follow-up spans ~900 days, including photometry and spectroscopy in the optical and IR, and (very long baseline interferometry) radio and X-ray observations. Radiative transfer modelling of the host galaxy spectral energy distribution shows the presence of a hitherto undetected dust obscured active galactic nucleus (AGN). The optical and near-IR spectra show broad 2000 km/s hydrogen, He I and O I emission features that decrease in flux over time. Radio imaging shows a fast evolving compact source of synchrotron emission spatially coincident with AT 2017gbl. We infer a lower limit for the radiated energy of 7.3 x 10^50 erg from the IR photometry. An extremely energetic supernova would satisfy this budget, but is ruled out by the radio counterpart evolution. Instead, we propose AT 2017gbl is related to an accretion event by the central supermassive black hole, where the spectral signatures originate in the AGN broad line region and the IR photometry is consistent with re-radiation by polar dust. Given the fast evolution of AT 2017gbl, we deem a tidal disruption event (TDE) of a star a more plausible scenario than a dramatic change in the AGN accretion rate. This makes AT 2017gbl the third TDE candidate to be hosted by a LIRG, in contrast to the so far considered TDE population discovered at optical wavelengths and hosted preferably by post-starburst galaxies.

We present new sub-arcsecond-resolution Karl G. Jansky Very Large Array (VLA) imaging at 10 GHz of 155 ultra-luminous ($L_{\rm bol}\sim10^{11.7-14.2} L_\odot$) and heavily obscured quasars with redshifts $z \sim0.4-3$. The sample was selected to have extremely red mid-infrared (MIR)-optical color ratios based on data from Wide-Field Infrared Survey Explorer (WISE) along with a detection of bright, unresolved radio emission from the NRAO VLA Sky Survey (NVSS) or Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) Survey. Our high-resolution VLA observations have revealed that the majority of the sources in our sample (93 out of 155) are compact on angular scales $<0.2^{\prime \prime}$ ($\leq 1.7$ kpc at $z \sim2$). The radio luminosities, linear extents, and lobe pressures of our sources are similar to young radio active galactic nuclei (AGN; e.g., Gigahertz Peaked Spectrum, GPS, and Compact Steep Spectrum, CSS, sources), but their space density is considerably lower. Application of a simple adiabatic lobe expansion model suggests relatively young dynamical ages ($\sim10^{4-7}$ years), relatively high ambient ISM densities ($\sim1-10^4$ cm$^{-3}$), and modest lobe expansion speeds ($\sim30-10,000$ km s$^{-1}$). Thus, we find our sources to be consistent with a population of newly triggered, young jets caught in a unique evolutionary stage in which they still reside within the dense gas reservoirs of their hosts. Based on their radio luminosity function and dynamical ages, we estimate only $\sim20\%$ of classical large scale FRI/II radio galaxies could have evolved directly from these objects. We speculate that the WISE-NVSS sources might first become GPS or CSS sources, of which some might ultimately evolve into larger radio galaxies.

The cosmic infrared background (CIB) provides a fundamental observational constraint on the star-formation history of galaxies over cosmic history. We estimate the contribution to the CIB from catalogued galaxies in the COSMOS field by using a novel map fitting technique on the \textit{Herschel} SPIRE maps. Prior galaxy positions are obtained using detections over a large range in wavelengths in the $K_{\rm s}$--3\,GHz range. Our method simultaneously fits the galaxies, the system foreground, and the leakage of flux from galaxies located in masked areas and corrects for an "over-fitting" effect not previously accounted for in stacking methods. We explore the contribution to the CIB as a function of galaxy survey wavelength and depth. We find high contributions to the CIB with the deep $r$ ($m_{\rm AB} \le 26.5$), $K_{\rm s}$ ($m_{\rm AB} \le 24.0$) and 3.6\,$\mu$m ($m_{\rm AB} \le 25.5$) catalogues. We combine these three deep catalogues and find a total CIB contributions of 10.5 $\pm$ 1.6, 6.7 $\pm$ 1.5 and 3.1 $\pm$ 0.7\,nWm$^{-2}$sr$^{-1}$ at 250, 350 and 500\,$\mu$m, respectively. Our CIB estimates are consistent with recent phenomenological models, prior based SPIRE number counts and with (though more precise than) the diffuse total measured by FIRAS. Our results raise the interesting prospect that the CIB contribution at $\lambda \le 500\,\mu$m from known galaxies has converged. Future large-area surveys like those with the Large Synoptic Survey Telescope are therefore likely to resolve a substantial fraction of the population responsible for the CIB at 250\,$\mu$m $\leq \lambda \leq$ 500\,$\mu$m.

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

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

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

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

Dust attenuation shapes the spectral energy distribution of galaxies. It is particularly true for dusty galaxies in which stars experience a heavy attenuation. The combination of UV-to-IR photometry with the spectroscopic measurement of the H$\alpha$ recombination line helps to quantify dust attenuation of the whole stellar population and its wavelength dependence. We selected an IR complete sample of galaxies in the COSMOS 3D-HST CANDELS field detected with the Herschel satellite with a signal to noise ratio larger than five. Optical to NIR photometry is available as well as NIR spectroscopy for each source. We reduced the sample to the redshift range $0.6 < z < 1.6$ to include the H$\alpha$ line in the G141 grism spectra. We have used a new version of the CIGALE code to fit simultaneously the continuum and H$\alpha$ line emission of the 34 selected galaxies. Using flexible attenuation laws with free parameters, we are able to measure the shape of the attenuation curve for each galaxy as well as the amount of attenuation of each stellar population, the former being in general steeper than the starburst law in the UV-optical with a large variation of the slope among galaxies. The attenuation of young stars or nebular continuum is found on average about twice the attenuation affecting older stars, again with a large variation. Our model with power-laws, based on a modification of the Charlot and Fall recipe, gives results in better agreement with the radiative transfer models than the global modification of the slope of the Calzetti law.

There is significant scientific value to be gained from combining AKARI fluxes with data at other far-IR wavelengths from the IRAS and Herschel missions. To be able to do this we must ensure that there are no systematic differences between the datasets that need to be corrected before the fluxes are compatible with each other. One such systematic effect identified in the BSCv1 data is the issue of beam corrections. We determine these for the BSCv2 data by correlating ratios of appropriate IRAS and AKARI bands with the difference in 2MASS J band extended and point source magnitudes for sources cross matched between the IRAS FSC, AKARI BSCv2 and 2MASS catalogs. We find significant correlations (p<< 10^-13) indicating that beam corrections are necessary in the 65 and 90 micron bands. We then use these corrected fluxes to supplement existing data in spectral energy distribution (SED) fits for ultraluminous infrared galaxies in the HERUS survey. The addition of AKARI fluxes makes little difference to the results of simple (T,beta) fits to the SEDs of these sources, though there is a general decrease in reduced chi^2 values. The utility of the extra AKARI data, however, is in allowing physically more realistic SED models with more parameters to be fit to the data. We also extend our analysis of beam correction issues in the AKARI data by examining the Herschel Reference Sample, which have Herschel photometry from 100 to 500 microns and which are more spatially extended than the HERUS ULIRGs. Thirty-four of the HRS sources have good Herschel SEDs and matching data from AKARI. This investigation finds that our simple 2MASS-based beam correction scheme is inadequate for these larger and more complex sources. There are also indications that additional beam corrections at 140 and 160 microns are needed for these sources, extended on scales >1 arcminute.

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH.

We present observations of NGC 1068 covering the $19.7-53.0$ $\mu$m wavelength range using FORCAST and HAWC+ onboard SOFIA. Using these observations, high-angular resolution infrared (IR) and sub-mm observations, we find an observational turn-over of the torus emission in the $30-40$ $\mu$m wavelength range with a characteristic temperature of $70-100$ K. This component is clearly different from the diffuse extended emission in the narrow line and star formation regions at 10-100 $\mu$m within the central 700 pc. We compute $2.2-432$ $\mu$m 2D images using the best inferred \textsc{clumpy} torus model based on several nuclear spectral energy distribution (SED) coverages. We find that when $1-20$ $\mu$m SED is used, the inferred result gives a small torus size ($<4$ pc radius) and a steep radial dust distribution. The computed torus using the $1-432$ $\mu$m SED provides comparable torus sizes, $5.1^{+0.4}_{-0.4}$ pc radius, and morphology to the recently resolved 432 $\mu$m ALMA observations. This result indicates that the $1-20$ $\mu$m wavelength range is not able to probe the full extent of the torus. The characterization of the turn-over emission of the torus using the $30-60$ $\mu$m wavelength range is sensitive to the detection of cold dust in the torus. The morphology of the dust emission in our 2D image at 432 $\mu$m is spatially coincident with the cloud distribution, while the morphology of the emission in the $1-20$ $\mu$m wavelength range shows an elongated morphology perpendicular to the cloud distribution. We find that our 2D \textsc{clumpy} torus image at 12 $\mu$m can produce comparable results to those observed using IR interferometry.

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

We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60${\mu}$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500${\mu}$m. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-IR spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parameterised modified black body model where temperature and emissivity $\beta$ are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9+/-4.7 K compared to 21.3+/-3.4 K) but a similar $\beta$ distribution (median $\beta$ = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al., 2014) galaxies. Dust masses are found to be in the range of 10^7.5 to 10^9 Msun significantly higher than that of Herschel Reference Sample (HRS) sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850${\mu}$m. These latter sources generally have cooler dust and/or redder 100-to-250 ${\mu}$m colours than our 60${\mu}$m-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population.

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

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

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

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

Core collapse supernova (CCSN) rates suffer from large uncertainties as many CCSNe exploding in regions of bright background emission and significant dust extinction remain unobserved. Such a shortfall is particularly prominent in luminous infrared galaxies (LIRGs), which have high star formation (and thus CCSN) rates and host bright and crowded nuclear regions, where large extinctions and reduced search detection efficiency likely lead to a significant fraction of CCSNe remaining undiscovered. We present the first results of project SUNBIRD (Supernovae UNmasked By Infra-Red Detection), where we aim to uncover CCSNe that otherwise would remain hidden in the complex nuclear regions of LIRGs, and in this way improve the constraints on the fraction that is missed by optical seeing-limited surveys. We observe in the near-infrared 2.15 {\mu}m $K_s$-band, which is less affected by dust extinction compared to the optical, using the multi-conjugate adaptive optics imager GeMS/GSAOI on Gemini South, allowing us to achieve a spatial resolution that lets us probe close in to the nuclear regions. During our pilot program and subsequent first full year we have discovered three CCSNe and one candidate with projected nuclear offsets as small as 200 pc. When compared to the total sample of LIRG CCSNe discovered in the near-IR and optical, we show that our method is singularly effective in uncovering CCSNe in nuclear regions and we conclude that the majority of CCSNe exploding in LIRGs are not detected as a result of dust obscuration and poor spatial resolution.

In this work we investigate the far-UV to NIR shape of the dust attenuation curve of a sample of IR selected dust obscured (U)LIRGs at z$\sim$2. The spectral energy distributions (SEDs) are fitted with CIGALE, a physically-motivated spectral synthesis model based on energy balance. Its flexibility allows us to test a wide range of different analytical prescriptions for the dust attenuation curve, including the well-known Calzetti and Charlot & Fall curves, and modified versions of them. The attenuation curves computed under the assumption of our reference double power-law model are in very good agreement with those derived, in previous works, with radiative transfer (RT) SED fitting. We investigate the position of our galaxies in the IRX-$\beta$ diagram and find this to be consistent with grayer slopes, on average, in the UV. We also find evidence for a flattening of the attenuation curve in the NIR with respect to more classical Calzetti-like recipes. This larger NIR attenuation yields larger derived stellar masses from SED fitting, by a median factor of $\sim$ 1.4 and up to a factor $\sim$10 for the most extreme cases. The star formation rate appears instead to be more dependent on the total amount of attenuation in the galaxy. Our analysis highlights the need for a flexible attenuation curve when reproducing the physical properties of a large variety of objects.

Nuclear starbursts and AGN activity are the main heating processes in luminous infrared galaxies (LIRGs) and their relationship is fundamental to understand galaxy evolution. In this paper, we study the star-formation and AGN activity of a sample of 11 local LIRGs imaged with subarcsecond angular resolution at radio (8.4GHz) and near-infrared ($2.2\mu$m) wavelengths. This allows us to characterize the central kpc of these galaxies with a spatial resolution of $\simeq100$pc. In general, we find a good spatial correlation between the radio and the near-IR emission, although radio emission tends to be more concentrated in the nuclear regions. Additionally, we use an MCMC code to model their multi-wavelength spectral energy distribution (SED) using template libraries of starburst, AGN and spheroidal/cirrus models, determining the luminosity contribution of each component, and finding that all sources in our sample are starburst-dominated, except for NGC6926 with an AGN contribution of $\simeq64$\%. Our sources show high star formation rates ($40$ to $167M_\odot\mathrm{yr}^{-1}$), supernova rates (0.4 to $2.0\mathrm{SN}\mathrm{yr}^{-1}$), and similar starburst ages (13 to $29\mathrm{Myr}$), except for the young starburst (9Myr) in NGC6926. A comparison of our derived star-forming parameters with estimates obtained from different IR and radio tracers shows an overall consistency among the different star formation tracers. AGN tracers based on mid-IR, high-ionization line ratios also show an overall agreement with our SED model fit estimates for the AGN. Finally, we use our wide-band VLA observations to determine pixel-by-pixel radio spectral indices for all galaxies in our sample, finding a typical median value ($\alpha\simeq-0.8$) for synchrotron-powered LIRGs.

We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least in most part from AGN activity. The AGN fractional luminosities are higher than those in either sub-millimeter galaxies, or AGN selected via other mid-infrared criteria. Adopting the $G$, M$_{20}$ and $A$ morphological parameters, together with traditional classification boundaries, infers that three quarters of the sample as mergers. Our sample do not, however, show any correlation between the considered morphological parameters and either infrared luminosity or AGN fractional luminosity. Moreover, their asymmetries and effective radii are distributed identically to those of massive galaxies at $z\sim2$. We conclude that our sample is not preferentially associated with mergers, though a significant merger fraction is still plausible. Instead, we propose that our sample are examples of the massive galaxy population at $z\sim2$ that harbor a briefly luminous, "flickering" AGN, and in which the $G$ and M$_{20}$ values have been perturbed, due to either the AGN, and/or the earliest formation stages of a bulge in an inside-out manner. Furthermore, we find that the mass assembly of the central black holes in our sample leads the mass assembly of any bulge component. Finally, we speculate that our sample represent a small fraction of the immediate antecedents of compact star-forming galaxies at $z\sim2$.

IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transient radio source with luminosity typical of bright supernovae, give further evidence of the ongoing star formation activity, which dominates the energetics of the system. At sub-parsec scales, the radio nucleus has a core-jet morphology with the jet being a newly ejected component showing a subluminal proper motion of 0.6c-1c. The AGN contributes less than two per cent of the total IR luminosity of the system. The corresponding Eddington factor is ~1E-3, suggesting this is a low-accretion rate engine, as often found in LLAGNs. However, its high bolometric luminosity (~10E44erg/s) agrees better with a normal AGN. This apparent discrepancy may just be an indication of the transition nature of the nucleus from a system dominated by star-formation, to an AGN-dominated system. The nucleus has a strongly inverted spectrum and a turnover at ~4.4GHz, thus qualifying as a candidate for the least luminous (L_5.0GHz ~ 6.3E28erg/s/Hz) and one of the youngest (~3000yr) gigahertz-peaked spectrum (GPS) sources. If the GPS origin for the IC883 nucleus is confirmed, then advanced mergers in the LIRG category are potentially key environments to unveil the evolution of GPS sources into more powerful radio galaxies.

We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671\mu$m with photometry in the SPIRE bands at 250$\mu$m, 350$\mu$m and 500$\mu$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $\mu$m, [CI] 609 $\mu$m, and [NII] 205 $\mu$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}

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

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

We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47<z<2.85), luminous (11.7 < log L(bol) (Lsun) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid-infrared (MIR) colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are detected at 870 microns, and we find that the sample has large mid- to far-infrared luminosity ratios consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are 24.7 < log P3.0 GHz (W/Hz) < 27.3, and all sources are radio-intermediate or radio-loud. BH mass estimates are 7.7 < log M(BH) (Msun) < 10.2. The rest frame 1-5 um SEDs are very similar to the "Hot DOGs" (Hot Dust Obscured Galaxies), and steeper (redder) than almost any other known extragalactic sources. ISM masses estimated for the ALMA detected sources are 9.9 < log M(ISM) (Msun) < 11.75 assuming a dust temperature of 30K. The cool dust emission is consistent with star formation rates (SFRs) reaching several thousand Msun/yr, depending on the assumed dust temperature, however we cannot rule out the alternative that the AGN powers all the emission in some cases. Our best constrained source has radiative transfer solutions with ~ equal contributions from an obscured AGN and a young (10-15 Myr) compact starburst.

We assume that dust near active galactic nuclei (AGN) is distributed in a torus-like geometry, which may be described by a clumpy medium or a homogeneous disk or as a combination of the two (i.e. a 2-phase medium). The dust particles considered are fluffy and have higher submillimeter emissivities than grains in the diffuse ISM. The dust-photon interaction is treated in a fully self-consistent three dimensional radiative transfer code. We provide an AGN library of spectral energy distributions (SEDs). Its purpose is to quickly obtain estimates of the basic parameters of the AGN, such as the intrinsic luminosity of the central source, the viewing angle, the inner radius, the volume filling factor and optical depth of the clouds, and the optical depth of the disk midplane, and to predict the flux at yet unobserved wavelengths. The procedure is simple and consists of finding an element in the library that matches the observations. We discuss the general properties of the models and in particular the 10mic. silicate band. The AGN library accounts well for the observed scatter of the feature strengths and wavelengths of the peak emission. AGN extinction curves are discussed and we find that there is no direct one-to-one link between the observed extinction and the wavelength dependence of the dust cross sections. We show that objects of the library cover the observed range of mid IR colors of known AGN. The validity of the approach is demonstrated by matching the SEDs of a number of representative objects: Four Seyferts and two quasars for which we present new Herschel photometry, two radio galaxies, and one hyperluminous infrared galaxy. Strikingly, for the five luminous objects we find pure AGN models fit the SED without a need to postulate starburst activity.

We complement our previous analysis of a sample of z~1-2 luminous and ultra-luminous infrared galaxies ((U)LIRGs), by adding deep VLA radio observations at 1.4 GHz to a large data-set from the far-UV to the sub-mm, including Spitzer and Herschel data. Given the relatively small number of (U)LIRGs in our sample with high S/N radio data, and to extend our study to a different family of galaxies, we also include 6 well sampled near IR-selected BzK galaxies at z~1.5. From our analysis based on the radiative transfer spectral synthesis code GRASIL, we find that, while the IR luminosity may be a biased tracer of the star formation rate (SFR) depending on the age of stars dominating the dust heating, the inclusion of the radio flux offers significantly tighter constraints on SFR. Our predicted SFRs are in good agreement with the estimates based on rest-frame radio luminosity and the Bell (2003) calibration. The extensive spectro-photometric coverage of our sample allows us to set important constraints on the SF history of individual objects. For essentially all galaxies we find evidence for a rather continuous SFR and a peak epoch of SF preceding that of the observation by a few Gyrs. This seems to correspond to a formation redshift of z~5-6. We finally show that our physical analysis may affect the interpretation of the SFR-M* diagram, by possibly shifting, with respect to previous works, the position of the most dust obscured objects to higher M* and lower SFRs.

The Luminous Infrared Galaxy NGC1614 hosts a prominent circumnuclear ring of star formation. However, the nature of the dominant emitting mechanism in its central ~100 pc is still under debate. We present sub-arcsecond angular resolution radio, mid-infrared, Pa-alpha, optical, and X-ray observations of NGC1614, aimed at studying in detail both the circumnuclear ring and the nuclear region. The 8.4 GHz continuum emission traced by the Very Large Array (VLA) and the Gemini/T-ReCS 8.7 micron emission, as well as the Pa-alpha line emission, show remarkable morphological similarities within the star-forming ring, suggesting that the underlying emission mechanisms are tightly related. We used an HST/NICMOS Pa-alpha map of similar resolution to our radio maps to disentangle the thermal free-free and non-thermal synchrotron radio emission, from which we obtained the intrinsic synchrotron power-law for each individual region within the central kpc of NGC1614. The radio ring surrounds a relatively faint, steep-spectrum source at the very center of the galaxy, suggesting that the central source is not powered by an AGN, but rather by a compact (r < 90 pc) starburst. Chandra X-ray data also show that the central kpc region is dominated by starburst activity, without requiring the existence of an AGN. We also used publicly available infrared data to model-fit the spectral energy distribution of both the starburst ring and a putative AGN in NGC1614. In summary, we conclude that there is no need to invoke an AGN to explain the observed bolometric properties of the galaxy.

We present a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies at $z>$4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their Spectral Energy Distributions (SED) and a composite model to fit the SEDs. The model includes a stellar emission component at $\lambda_{\rm rest} <$ 3.5$ \mu$m; a hot dust component peaking at $\lambda_{rest} \sim$ 5$\,\mu$m; and cold dust component which becomes significant for $\lambda_{\rm rest} >$ 50$\,\mu$m. Six objects in the sample are detected at 250 and 350$ \mu$m. The dust temperatures for the sources in this sample are in the range of 40$-$80 K, and their $L_{\rm FIR}$ $\sim$ 10$^{13}$ L$_{\odot}$ qualifies them as Hyper$-$Luminous Infrared Galaxies (HyperLIRGs). The mean FIR-radio index for this sample is around $< q > = 2.2$ indicating no radio excess in their radio emission. Most sources in the sample have 24$ \mu$m detections corresponding to a rest-frame 4.5$ \mu$m luminosity of Log$_{10}$(L$_{4.5}$ / L$_{\odot}$) = 11 $\sim$ 11.5. Their L$_{\rm 4.5}$/$L_{\rm FIR}$ ratios are very similar to those of starburst dominated submillimeter galaxies at $z \sim$ 2. The $L_{\rm CO}-L_{\rm FIR}$ relation for this sample is consistent with that determined for local ULIRGs and SMGs at $z \sim$ 2. We conclude that submillimeter galaxies at $z >$ 4 are hotter and more luminous in the FIR, but otherwise very similar to those at $z \sim$ 2. None of these sources show any sign of the strong QSO phase being triggered.

We present a revised model for the infrared emission of the hyperluminous infrared galaxy IRAS F10214+4724 which takes into account recent photometric data from Spitzer and Herschel that sample the peak of its spectral energy distribution. We first present and discuss a grid of smooth active galactic nucleus (AGN) torus models computed with the method of Efstathiou & Rowan-Robinson and demonstrate that the combination of these models and the starburst models of Efstathiou and coworkers, while able to give an excellent fit to the average spectrum of Seyfert 2s and spectra of individual type 2 quasars measured by Spitzer, fails to match the spectral energy distribution of IRAS F10214+4724. This is mainly due to the fact that the nuSnu distribution of the galaxy falls very steeply with increasing frequency (a characteristic that is usually indicative of heavy absorption by dust) but shows a silicate feature in emission. Such emission features are not expected in sources with optical/near-infrared type 2 AGN spectral signatures. The Herschel data show that there is more power emitted in the rest-frame 20-50um wavelength range compared with the model presented by Efstathiou which assumes three components of emission: an edge-on torus, clouds (at a temperature of 610 and 200K) that are associated with the narrow-line region (NLR) and a highly obscured starburst that dominates in the submillimetre. We present a revised version of that model that assumes an additional component of emission which we associate with NLR clouds at a temperature of 100K. The 100K dust component could also be explained by a highly obscured hot starburst. The model suggests that the NLR of IRAS F10214+4724 has an unusually high covering factor (>= 17%) or more likely the magnification of the emission from the NLR clouds is significantly higher than that of the emission from the torus.

We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model the AGN contributes about 90% of the total luminosity of 1.1 x 10^13 Lo, which is about a factor of five higher than previous estimates. The large correction of the luminosity is due to the anisotropy of the emission of the best fit torus. Similar corrections may be necessary for other local and high-z analogs. This correction implies that IRAS 08572+3915 at a redshift of 0.05835 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [CII] to IR luminosity (log L_[CII]/L_{IR} < -3.8) and a [OI]63um to [CII]158um line ratio of about 1 that supports the model presented in this letter.

We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

We report on our observations of the 79 and 119um doublet transitions of OH for 24 local (z<0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH119 profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission component decreases as the silicate absorption increases. This locates the OH outflows inside the obscured nuclei. The maximum outflow velocities for our sources range from less than 100 to 2000 km/s, with 15/24 (10/24) sources showing OH absorption at velocities exceeding 700 km/s (1000 km/s). Three sources show maximum OH outflow velocities exceeding that of Mrk231. Since outflow velocities above 500-700 km/s are thought to require an active galactic nucleus (AGN) to drive them, about 2/3 of our ULIRG sample may host AGN-driven molecular outflows. This finding is supported by the correlation we find between the maximum OH outflow velocity and the IR-derived bolometric AGN luminosity. No such correlation is found with the IR-derived star formation rate. The highest outflow velocities are found among sources which are still deeply embedded. We speculate that the molecular outflows in these sources may be in an early phase of disrupting the nuclear dust veil before these sources evolve into less obscured AGN. Four of our sources show high-velocity wings in their [C II] fine-structure line profiles implying neutral gas outflow masses of at least 2-4.5 x 10^8 Msun.

We estimate the fraction of core-collapse supernovae (CCSNe) that remain undetected by optical SN searches due to obscuration by large amounts of dust in their host galaxies. This effect is especially important in luminous and ultraluminous infrared galaxies, which are locally rare but dominate the star formation at redshifts of z~1-2. We perform a detailed investigation of the SN activity in the nearby luminous infrared galaxy Arp 299 and estimate that up to 83% of the SNe in Arp 299 and in similar galaxies in the local Universe are missed by observations at optical wavelengths. For rest-frame optical surveys we find the fraction of SNe missed due to high dust extinction to increase from the average local value of ~19% to ~38% at z~1.2 and then stay roughly constant up to z~2. It is therefore crucial to take into account the effects of obscuration by dust when determining SN rates at high redshift and when predicting the number of CCSNe detectable by future high-z surveys such as LSST, JWST, and Euclid. For a sample of nearby CCSNe (distances 6-15 Mpc) detected during the last 12 yr, we find a lower limit for the local CCSN rate of 1.5 +0.4/-0.3 x 10^-4 yr^-1 Mpc^-3, consistent with that expected from the star formation rate. Even closer, at distances less than ~6 Mpc, we find a significant increase in the CCSN rate, indicating a local overdensity of star formation caused by a small number of galaxies that have each hosted multiple SNe.

The mid-infrared spectra of ultraluminous infrared galaxies (ULIRGs) contain a variety of spectral features that can be used as diagnostics to characterise the spectra. However, such diagnostics are biased by our prior prejudices on the origin of the features. Moreover, by using only part of the spectrum they do not utilise the full information content of the spectra. Blind statistical techniques such as principal component analysis (PCA) consider the whole spectrum, find correlated features and separate them out into distinct components. We further investigate the principal components (PCs) of ULIRGs derived in Wang et al.(2011). We quantitatively show that five PCs is optimal for describing the IRS spectra. These five components (PC1-PC5) and the mean spectrum provide a template basis set that reproduces spectra of all z<0.35 ULIRGs within the noise. For comparison, the spectra are also modelled with a combination of radiative transfer models of both starbursts and the dusty torus surrounding active galactic nuclei. The five PCs typically provide better fits than the models. We argue that the radiative transfer models require a colder dust component and have difficulty in modelling strong PAH features. Aided by the models we also interpret the physical processes that the principal components represent. The third principal component is shown to indicate the nature of the dominant power source, while PC1 is related to the inclination of the AGN torus. Finally, we use the 5 PCs to define a new classification scheme using 5D Gaussian mixtures modelling and trained on widely used optical classifications. The five PCs, average spectra for the four classifications and the code to classify objects are made available at: http://www.phys.susx.ac.uk/~pdh21/PCA/

The high star formation rates of luminous infrared galaxies (LIRGs) make them ideal places for core-collapse supernova (CCSN) searches. At radio frequencies, free from dust extinction, it is possible to detect compact components within the innermost LIRG nuclear regions, such as SNe and SN remnants, as well as AGN buried deep in the LIRG nuclei. We studied the LIRG IC883 aiming at: (i) investigating its (circum-)nuclear regions using the e-EVN at 5GHz, and e-MERLIN at 6.9GHz, complemented by archival VLBI data; (ii) detecting at radio frequencies the two recently reported circumnuclear SNe 2010cu and 2011hi, which were discovered by near-IR (NIR) adaptive optics observations of IC883; and (iii) further investigating the nature of SN2011hi at NIR by means of observations with Gemini-North. The circumnuclear regions traced by e-MERLIN at 6.9GHz have an extension of ~1kpc, and show a striking double-sided structure, which very likely corresponds to a warped rotating ring, in agreement with previous studies. Our e-EVN observations at 5GHz and complementary archival VLBI data at 5GHz and 8.4GHz, reveal the presence of various milliarcsec compact components in the nucleus of IC883. A single compact source, an AGN candidate, dominates the emission at both nuclear and circumnuclear scales, as imaged with the e-EVN and e-MERLIN, respectively. The other milliarcsec components are very suggestive of ongoing nuclear CCSN activity. Our e-EVN observations also resulted in upper limits to the radio luminosity of the two SNe in IC883 recently discovered at NIR. We refine the classification of SN2011hi as a Type IIP SN according to our latest Gemini-North epoch from 2012, in agreement with a low-luminosity radio SN nature. We estimate a CCSN rate lower limit of 1.1_{-0.6}^{+1.3} yr^{-1} for the entire galaxy, based on three nuclear radio SNe and the circumnuclear SNe 2010cu and 2011hi. (abridged)

We present a new model for the infrared emission of the high redshift hyperluminous infrared galaxy IRAS F10214+4724 which takes into account recent photometric data from Spitzer and Herschel that sample the peak of its spectral energy distribution. We first demonstrate that the combination of the AGN tapered disc and starburst models of Efstathiou and coworkers, while able to give an excellent fit to the average spectrum of type 2 AGN measured by Spitzer, fails to match the spectral energy distribution of IRAS F10214+4724. This is mainly due to the fact that the nuSnu distribution of the galaxy falls very steeply with increasing frequency (a characteristic of heavy absorption by dust) but shows a silicate feature in emission. We propose a model that assumes two components of emission: clouds that are associated with the narrow-line region and a highly obscured starburst. The emission from the clouds must suffer significantly stronger gravitational lensing compared to the emission from the torus to explain the observed spectral energy distribution.

We use SDSS spectra and optical to far-infrared photometry for a sample of 31 FeLoBAL QSOs to study the relationship between the AGN-driven outflows, and obscured star formation in their host galaxies. We find that FeLoBAL QSOs invariably have IR luminosities exceeding 10^12 Solar luminosities. The AGN supplies 75% of the total IR emission, on average, but with a range from 20% to 100%. We find a clear anticorrelation between the strength of the AGN-driven outflows and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN act to curtail obscured star formation in the host galaxies of reddened QSOs to less than ~25% of the total IR luminosity. This is the most direct evidence yet obtained for `quasar mode' AGN feedback.

We present optical to far-infrared photometry of 31 reddened QSOs that show evidence for radiatively driven outflows originating from AGN in their rest-frame UV spectra. We use these data to study the relationships between the AGN-driven outflows, and the AGN and starburst infrared luminosities. We find that FeLoBAL QSOs are invariably IR-luminous, with IR luminosities exceeding 10^{12} Solar luminosities in all cases. The AGN supplies 76% of the total IR emission, on average, but with a range from 20% to 100%. We find no evidence that the absolute luminosity of obscured star formation is affected by the AGN-driven outflows. Conversely, we find an anticorrelation between the strength of AGN-driven outflows, as measured from the range of outflow velocities over which absorption exceeds a minimal threshold, and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no convincing evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN can have a dramatic, negative impact on luminous star formation in their host galaxies. We find that such outflows act to curtail star formation such that star formation contributes less than ~25% of the total IR luminosity. We also propose that the degree to which termination of star formation takes place is not deducible from the IR luminosity of the AGN.

We report the discovery of two consecutive supernovae (SNe), 2010cu and 2011hi, located at 0.37" (180 pc) and 0.79" (380 pc) projected distance respectively from the centre of the K-band nucleus of the luminous infrared galaxy IC 883. The SNe were discovered in an ongoing near-infrared K-band search for core-collapse SNe in such galaxies using the ALTAIR/NIRI adaptive optics system with laser guide star at the Gemini-North Telescope. These are thus the closest SNe yet discovered to a LIRG nucleus in optical or near-infrared wavelengths. The near-infrared light curves and colours of both SNe are consistent with core-collapse events. Both SNe seem to suffer from relatively low host galaxy extinction suggesting that regardless of their low projected galactocentric distances, they are not deeply buried in the nuclear regions of the host galaxy.

We present follow-up observations of the far-infrared (FIR) sources at 90, 150 and 180 mu detected as part of the ISOPHOT EBL project, which has recently measured the absolute surface brightness of the cosmic infrared background radiation (CIRB) for the first time independently from COBE data. We have observed the fields at the North Galactic Pole region in the optical and near-IR, and complement these data with SDSS photometry, and spectroscopy where available, and present identifications of the 25 FIR sources which reach down to ~150 mJy in all three ISOPHOT bands. Identifications are done by means of full spectral energy density fitting to all sources in the FIR error circle areas. Approximately 80 per cent are identified as star-forming or star-bursting galaxies at z<0.3. We also find that more than half of the counterparts have disturbed morphologies, with signs of past or present interactions. However, only 20 per cent of all the sources are uniquely matched with a single galaxy -- 40 per cent are blends of two or more of these nearby star-forming galaxies, while another 20 per cent are likely blends of nearby and fainter galaxies. The final 20 per cent are likely to be more luminous IR galaxies at higher redshifts. The blended sources have an effect on the FIR source counts. In particular, taking into account realistic confusion or blending of sources, the differential FIR counts move down by a factor of ~1.5 and steepen in the 100 to 400 mJy range.

We explore whether our models for starbursts, quiescent star-forming galaxies and for AGN dust tori are able to model the full range of IRS spectra measured with Spitzer. The diagnostic plot of 9.7 mu silicate optical depth versus 6.2 mu PAH equivalent width, introduced by Spoon and coworkers in 2007, gives a good indication of the age and optical depth of a starburst, and of the contribution of an AGN dust torus. However there is aliasing between age and optical depth at later times in the evolution of a starburst, and between age and the presence of an AGN dust torus. Modeling the full IRS spectra and using broad-band 25-850 mu fluxes can help to resolve these aliases. The observed spectral energy distributions require starbursts of a range of ages with initial dust optical depth ranging from 50-200, optically thin dust emission ('cirrus') illuminated by a range of surface brightnesses of the interstellar radiation field, and AGN dust tori with a range of viewing angles.

We present radiative transfer models for submillimeter galaxies with spectroscopic redshifts and mid-infrared spectroscopy from Spitzer/IRS and analyze available Spitzer/MIPS 24, 70 and 160mu data. We use two types of starburst models, a cirrus model and a model for the emission of an AGN torus in order to investigate the nature of these objects. We find that for three of the objects (25%) cirrus emission alone can account for the mid-infrared spectrum and the MIPS and submillimeter data. For the remaining objects we find that we need a combination of starburst and cirrus in order to fit simultaneously the multi--wavelength data. We find that the typical submillimeter galaxy has comparable luminosity in the starburst (median L=10^12.5 Lo) and cirrus (median L=10^12.4 Lo) components. This could arise if the galaxies have been forming stars continuously for the last 250Myr with the star formation occurring in the last 5Myr being shrouded by high-optical-depth molecular cloud dust, whereas the rest of the starlight is attenuated by diffuse dust or cirrus with an Av of about 1mag.

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 report the discovery of a confirmed supernova (SN) and a supernova-candidate in near-infrared images from the ALTAIR/NIRI adaptive optics system on the Gemini-North Telescope and NICMOS on the Hubble Space Telescope. The Gemini images were obtained as part of a near-infrared K-band search for highly-obscured SNe in the nuclear regions of luminous infrared galaxies. SN 2008cs apparent in the Gemini images is the first SN discovered using laser guide star adaptive optics. It is located at 1500 pc projected distance from the nucleus of the luminous infrared galaxy IRAS 17138-1017. The SN luminosity, JHK colors and light curve are consistent with a core-collapse event suffering from a very high host galaxy extinction of 15.7 +- 0.8 magnitudes in V-band which is to our knowledge the highest yet measured for a SN. The core-collapse nature of SN 2008cs is confirmed by its radio detection at 22.4 GHz using our Very Large Array observations 28 days after the SN discovery, indicating a prominent interaction of the SN ejecta with the circumstellar medium. An unconfirmed SN apparent in the NICMOS images from 2004 is located in the same galaxy at 660 pc projected distance from the nucleus and has a lower extinction.

(abridged) We present near-infrared (NIR) adaptive optics imaging obtained with VLT/NACO and optical spectroscopy from the Southern African Large Telescope (SALT) of a luminous infrared galaxy (LIRG) IRAS 19115-2124. These data are combined with archival HST imaging and Spitzer imaging and spectroscopy, allowing us to study this disturbed interacting/merging galaxy, dubbed the Bird, in extraordinary detail. In particular, the data reveal a triple system where the LIRG phenomenon is dominated by the smallest of the components. One nucleus is a regular barred spiral with significant rotation, while another is highly disturbed with a surface brightness distribution intermediate to that of disk and bulge systems, and hints of remaining arm/bar structure. We derive dynamical masses in the range 3-7x10^10 M_solar for both. The third component appears to be a 1-2x10^10 M_solar mass irregular galaxy. The total system exhibits HII galaxy-like optical line ratios and strengths, and no evidence for AGN activity is found from optical or mid-infrared data. The star formation rate is estimated to be 190 M_solar/yr. We search for SNe, super star clusters, and detect 100-300 km/s outflowing gas from the Bird. Overall, the Bird shows kinematic, dynamical, and emission line properties typical for cool ultra luminous IR galaxies. However, the interesting features setting it apart for future studies are its triple merger nature, and the location of its star formation peak - the strongest star formation does not come from the two major K-band nuclei, but from the third irregular component. Aided by simulations, we discuss scenarios where the irregular component is on its first high-speed encounter with the more massive components.

We report a supernova discovery in Ks-band images from the NAOS CONICA adaptive optics (AO) system on the ESO Very Large Telescope (VLT). The images were obtained as part of a near-infrared search for highly-obscured supernovae in the nuclear regions of luminous and ultraluminous infrared galaxies. SN 2004ip is located within a circumnuclear starburst at 1.4 arcsec (or 500 pc) projected distance from the K-band nucleus of the luminous infrared galaxy IRAS 18293-3413. The supernova luminosity and light curve are consistent with a core-collapse event suffering from a host galaxy extinction of up to about 40 magnitudes in V-band which is as expected for a circumnuclear starburst environment. This is the first supernova to be discovered making use of AO correction and demonstrates the potential of the current 8-meter class telescopes equipped with AO in discovering supernovae from the innermost nuclear regions of luminous and ultraluminous infrared galaxies.

A model for the infrared emission of the high redshift ultraluminous infrared galaxy FSC 10214+4724 is presented. The model assumes three components of emission: a dusty torus viewed edge-on, clouds that are associated with the narrow-line region and a highly obscured starburst. It is demonstrated that the presence of clouds in the narrow-line region, with a covering factor of 17%, can explain why the mid-infrared spectrum of FSC 10214+4724 shows a silicate feature in emission despite the fact that its torus is viewed edge-on. It is also shown that the same model, but with the torus viewed face-on, predicts a spectrum with silicate emission features that is characteristic of the spectra of quasars recently observed with Spitzer.

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

Using the large multi-wavelength data set in the chandra/SWIRE Survey (0.6 square degrees in the Lockman Hole), we show evidence for the existence of highly obscured (Compton-thick) AGN, estimate a lower limit to their surface density and characterize their multi-wavelength properties. Two independent selection methods based on the X-ray and infrared spectral properties are presented. The two selected samples contain 1) 5 X-ray sources with hard X-ray spectra and column densities > 10^24 cm-2, and 2) 120 infrared sources with red and AGN-dominated infrared spectral energy distributions (SEDs). We estimate a surface density of at least 25 Compton-thick AGN per square degree detected in the infrared in the chandra/SWIRE field of which ~40% show distinct AGN signatures in their optical/near-infrared SEDs, the remainings being dominated by the host-galaxy emission. Only ~33% of all Compton-thick AGN are detected in the X-rays at our depth (F(0.3-8 keV)>10^-15 erg/cm2/s. We report the discovery of two sources in our sample of Compton-thick AGN, SWIRE_J104409.95+585224.8 (z=2.54) and SWIRE_J104406.30+583954.1 (z=2.43), which are the most luminous Compton-thick AGN at high-z currently known. The properties of these two sources are discussed in detail with an analysis of their spectra, SEDs, luminosities and black-hole masses.

Aims. Filamentary structures of early type stars are found to be a common feature of the Magellanic Clouds formed at an age of about 0.9-2*10^8 yr. As we go to younger ages these large structures appear fragmented and sooner or later form young clusters and associations. In the optical domain we have detected 56 such large structures of young objects, known as stellar complexes in the LMC for which we give coordinates and dimensions. We also investigate star formation activity and evolution of these stellar complexes and define the term "starburst region". Methods. IR properties of these regions have been investigated using IRAS data. A colour-magnitude diagram (CMD) and a two-colour diagram from IRAS data of these regions ware compared with observations of starburst galaxies and cross-matching with HII regions and SNRs was made . Radio emission maps at 8.6-GHz and the CO (1 to 0) line were also cross correlated with the map of the stellar complexes. Results. It has been found that nearly 1/3 of the stellar complexes are extremely active resembling the IR behaviour of starburst galaxies and HII regions. These stellar complexes illustrating such properties are called here "starburst regions". They host an increased number of HII regions and SNRs. The main starburst tracers are their IR luminosity (F60 well above 5.4 Jy) and the 8.6-GHz radio emission. Finally the evolution of all stellar complexes is discussed based on the CO emission.

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 ISO mid-infrared spectrophotometry and far-infrared photometry of galaxies with Hidden Broad Line Regions (HBLR). We also present radiative transfer models of their spectral energy distributions which enable us to separate the contributions from the dusty disc of the AGN and the dusty starbursts. We find that the combination of tapered discs (discs whose thickness increases with distance from the central source in the inner part but stays constant in the outer part) and starbursts provide good fits to the data. The tapered discs dominate in the mid-infrared part of the spectrum and the starbursts in the far-infrared. After correcting the AGN luminosity for anisotropic emission we find that the ratio of the AGN luminosity to the starburst luminosity, L(AGN)/L(SB), ranges from about unity for IRAS14454-4343 to about 13 for IRAS01475-0740. Our results suggest that the warm IRAS colours of HBLR are due to the relatively high L(AGN)/L(SB). Our fits are consistent with the unified model and the idea that the infrared emission of AGN is dominated by a dusty disc in the mid-infrared and starbursts in the far-infrared.

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.

Present AGN models aiming to account for the observed infrared SEDs consider a physical description of the dust and a solution of the radiative transfer problem. Mid infrared spectra obtained at different spatial scales are presented. They show that PAH bands are detected in starburst regions but significantly reduced near the centre of AGN. This may be explained by examining the heating mechanism of PAHs after hard photon interactions. On the radiative transfer side first a most economic model is presented where three parameters, luminosity, effective size and extinction of the nucleus are varied to obtain SED fits. A full grid of model spectra are made available at: "http://www.eso.org/~rsiebenm/agn_models". This model is sufficient to account for ISO broad band data of a sample of 68 radio galaxies and quasars of the 3CR catalogue. The hot dust component detected is mainly due to small grains and PAHs. In such models, type 1 AGNs are represented by a compact dust distribution with warm grains and weak PAH emission. In AGNs of type 2, the dust appears to be more extended, relatively colder and PAH bands are strong. Realistic AGN models which are consistent with the unification need to explain the overall absence of the 9.7mic. silicate emission feature. This can be done by considering various geometries (tapered discs). Models which combine AGN and starburst activity are presented for galaxies with hidden broad line region. It is found that the AGN torus dominate the mid IR continuum emission and that the starbursts dominate the PAH band as well as the far infrared and submillimeter emission.

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

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

We present a model for the UV to submillimeter emission from stars embedded in the general interstellar dust in galaxies (the 'infrared cirrus' component). Such emission is characterized by relatively low optical depths of dust and by cool (< 30 K) dust temperatures. The model incorporates the stellar population synthesis model of Bruzual & Charlot and the dust model of Siebenmorgen & Krugel which includes the effects of small grains/PAHs. We apply the model to fit the optical to submillimeter spectral energy distributions (SEDs) of nearby galaxies which are dominated by cirrus emission and find that our simple model is quite adequate to explain the observed SEDs. We also, more controversially, apply this cirrus model to the SEDs of high redshift sources detected in blank field submillimeter surveys with SCUBA. Surprisingly, an excellent fit is found for many of these sources, with typical values for the optical depth A_V and the surface brightness of the stellar radiation field psi being only a factor 2-3 higher than for nearby galaxies. This increase is not unreasonable given the expected evolution of dust optical depth in currently favoured star-formation history models. We conclude that the tendency to interpret the high-z SCUBA galaxies as very highly obscured starbursts may be premature and that these galaxies may be more closely linked to optically selected high redshift galaxies than previously assumed.

(Abridged) We examine the power source of 41 local Ultraluminous Infrared Galaxies using archival infrared and optical photometry. We fit the observed Spectral Energy Distributions (SEDs) with starburst and AGN components; each component being drawn from a family of templates. We find all of the sample require a starburst, whereas only half require an AGN. In 90% of the sample the starburst provides over half the IR emission, with a mean fractional luminosity of 82%. When combined with other galaxy samples we find that starburst and AGN luminosities correlate over 6 decades in IR luminosity suggesting that a common factor governs both luminosities, plausibly the gas masses in the nuclear regions. We find that the mid-IR 7.7 micron line-continuum ratio is no indication of the starburst luminosity, or the fractional AGN luminosity, and therefore that this ratio is not a reliable diagnostic of the power source in ULIRGs. We propose that the scatter in the radio-IR correlation in ULIRGs is due to a skewed starburst IMF and/or relic relativistic electrons from a previous starburst, rather than contamination from an obscured AGN. We show that most ULIRGs undergo multiple starbursts during their lifetime, and by inference that mergers between more than two galaxies may be common amongst ULIRGs. Our results support the evolutionary model for ULIRGs proposed by Farrah et al 2001, where they can follow many different evolutionary paths of starburst and AGN activity in transforming merging spiral galaxies into elliptical galaxies, but that most do not go through an optical QSO phase. The lower level of AGN activity in our local sample than in z~1 HLIRGs implies that the two samples are distinct populations. We postulate that different galaxy formation processes at high-z are responsible for this difference.

We present nuclear spectral energy distributions (SEDs) in the range 0.4-16micron for an expanded CfA sample of Seyfert galaxies. The spectral indices from 1-16micron range from alpha_IR 0.9 to 3.8. The shapes of the spectra are correlated with Seyfert type in the sense that steeper nuclear SEDs (nu*f_nu increasing with increasing wavelength) tend to be found in Seyfert 2s and flatter SEDs (nu *f_nu constant) in Seyfert 1-1.5s. The galaxies optically classified as Seyferts 1.8s and 1.9s display values of alpha_IR as in type 1 objects, or values intermediate between those of Seyfert 1s and Seyfert 2s. The intermediate SEDs of many Seyfert 1.8-1.9s may be consistent with the presence of a pure Seyfert 1 viewed through a moderate amount (A_V <5mag) of foreground galaxy extinction. We find, however, that between 10 and 20% of galaxies with broad optical line components have steep infrared SEDs. Torus models usually adopt high equatorial opacities to reproduce the infrared properties of Seyfert 1s and 2s, resulting in a dichotomy of infrared SEDs (flat for type 1s, and steep for type 2s). Such a dichotomy, however, is not observed in our sample. The wide range of spectral indices observed in the type 2 objects, the lack of extremely steep SEDs, and the large numbers of objects with intermediate spectral indices cannot be reconciled with predictions from existing optically thick torus models. We discuss possible modifications to improve torus models, including low optical depth tori, clumpy dusty tori, and high-optical-depth tori with an extended optically thin component.

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.

We present sub-mm photometry for 11 Hyperluminous Infrared Galaxies (HLIRGs) and use radiative transfer models for starbursts and AGN to investigate the IR emission. In all sources both a starburst and AGN are required to explain the IR emission. The mean starburst fraction is 35%, with a range spanning 80% starburst dominated to 80% AGN dominated. In all cases the starburst dominates at rest-frame wavelengths >50 microns, with star formation rates >500 solar masses per year. The trend of increasing AGN fraction with increasing IR luminosity seen in IRAS galaxies peaks in HLIRGs, and is not higher than the fraction seen in bright ULIRGs. The AGN and starburst luminosities correlate, suggesting that a common physical factor, plausibly the dust masses, governs their luminosities. Our results suggest that the HLIRG population is comprised both of ULIRG-like galaxy mergers, and of young galaxies going through their maximal star formation periods whilst harbouring an AGN. The coeval AGN and starburst activity in our sources implies that starburst and AGN activity, and the peak starburst and AGN luminosities, can be coeval in active galaxies generally. When extrapolated to high-z our sources have comparable sub-mm fluxes to sub-mm survey sources. At least some sub-mm survey sources are therefore likely to be comprised of similar galaxy populations to those found in the HLIRG population. It is also plausible from these results that high-z sub-mm sources harbour heavily obscured AGN. The differences in X-ray and sub-mm properties between HLIRGs at z~1 and sub-mm sources at z~3 implies evolution between the two epochs. Either the mean AGN obscuration level is greater at z~3 than at z~1, or the fraction of IR-luminous sources at z~3 that contain AGN is smaller than that at z~1.

We present 7-180micron photometry of a sample of hyperluminous infrared galaxies (HyLIGs) obtained with the photometer and camera mounted on the Infrared Space Observatory. We have used radiative transfer models of obscured starbursts and dusty torii to model their spectral energy distributions (SEDs). We find that F00235+1024, F14218+3845 and F15307+3252 require a combination of starburst and AGN components to explain their mid to far-infrared emission, while for TXS0052+471 a dust torus AGN model alone is sufficient. For F00235+1024 and F14218+3845 the starburst component is the predominant contributor whereas for F15307+3252 the dust torus component dominates. The implied star formation rates (SFR) for these three sources estimated from their infrared luminosities are $\dot{M}_{*,all} > 3000M_{\odot}yr^{-1}h_{50}^{-2}$ and are amongst the highest SFRs estimated to date. We also demonstrate that the well-known radio-FIR correlation extends into both higher radio and infrared power than previously investigated. The relation for HyLIGs has a mean q value of ~1.94. The results of this study imply that better sampling of the IR spectral energy distributions of HyLIgs may reveal that both AGN and starburst components are required to explain all the emission from the NIR to the sub-millimetre.

We present results from a deep mid-infrared survey of the Hubble Deep Field South (HDF-S) region performed at 7 and 15 micron with the CAM instrument on board the Infrared Space Observatory (ISO). The final map in each band was constructed by the coaddition of four independent rasters, registered using bright sources securely detected in all rasters, with the absolute astrometry being defined by a radio source detected at both 7 and 15 micron. We sought detections of bright sources in a circular region of radius 2.5 arcmin at the centre of each map, in a manner that simulations indicated would produce highly reliable and complete source catalogues using simple selection criteria. Merging source lists in the two bands yielded a catalogue of 35 distinct sources, which we calibrated photometrically using photospheric models of late-type stars detected in our data. We present extragalactic source count results in both bands, and discuss the constraints they impose on models of galaxy evolution models, given the volume of space sampled by this galaxy population.

We present results from a deep mid-IR survey of the Hubble Deep Field South (HDF-S) region performed at 7 and 15um with the CAM instrument on board ISO. We found reliable optical/near-IR associations for 32 of the 35 sources detected in this field by Oliver et al. (2002, Paper I): eight of them were identified as stars, one is definitely an AGN, a second seems likely to be an AGN, too, while the remaining 22 appear to be normal spiral or starburst galaxies. Using model spectral energy distributions (SEDs) of similar galaxies, we compare methods for estimating the star formation rates (SFRs) in these objects, finding that an estimator based on integrated (3-1000um) IR luminosity reproduces the model SFRs best. Applying this estimator to model fits to the SEDs of our 22 spiral and starburst galaxies, we find that they are forming stars at rates of ~1-100 M_sol/yr, with a median value of ~40M_sol/yr, assuming an Einstein - de Sitter universe with a Hubble constant of 50 km/s/Mpc, and star formation taking place according to a Salpeter (1955) IMF across the mass range 0.1-100M_sol. We split the redshift range 0.0<z<0.6 into two equal-volume bins to compute raw estimates of the star formation rate density contributed by these sources, assuming the same cosmology and IMF as above and computing errors based on estimated uncertainties in the SFRs of individual galaxies. We compare these results with other estimates of the SFR density made with the same assumptions, showing them to be consistent with the results of Flores et al. (1999) from their ISO survey of the CFRS 1415+52 field. However, the relatively small volume of our survey means that our SFR density estimates suffer from a large sampling variance, implying that our results, by themselves, do not place tight constraints on the global mean SFR density.

We present a cross-correlation between the ELAIS 15$\mu$m ISO survey with the {\it ROSAT} (0.1-2 keV) all-sky survey and the pointed observations WGACAT source catalogue. The resulting sample contains 15 objects. Optical spectroscopic identifications exist for 13 objects: 6 broad-line QSOs, 4 narrow line galaxies or type-2 AGN (NLG) and 3 stars. We have used both the X-ray to IR luminosity flux ratio $f_{x}/f_{IR}$ and the X-ray hardness ratios diagnostics to estimate the amount of obscuration in these objects. The X-ray spectrum of the narrow-line galaxies does not present strong evidence for obscuration; however, the low $f_{x}/f_{IR}$ ratio combined with the high X-ray luminosities suggest that at least one of the narrow line galaxies is associated with an obscured Seyfert nucleus. 4 out of 6 QSOs present high $f_{x}/f_{IR}$ ratios and steep X-ray spectra with $\Gamma>2$. One QSO (ELAISC15-J133442+375736) at a redshift of z=1.89, has an abnormally low X-ray/IR flux ratio, with its infrared luminosity approaching that of an hyperluminous galaxy ($\sim10^{12.98} h^{-2} L_\odot$). Finally, one radio-loud QSO is the hardest X-ray source in our sample, presenting strong evidence for a high absorbing column ($N_H\sim10^{22}\rm cm^{-2}$).

We present the results of a multi-wavelength study of the 19 most significant sub-mm sources detected in the SCUBA 8-mJy survey. As described in Scott et al. (2001), this survey covers ~260 arcmin^2 using the sub-millimetre camera SCUBA, to a limiting source detection limit S(850um) ~ 8 mJy. One advantage of this relatively bright flux-density limit is that accurate astrometric positions are potentially achievable for every source using existing radio and/or mm-wave interferometers. However, an associated advantage is that SED-based redshift constraints should be more powerful than in fainter sub-mm surveys. Here we therefore exploit the parallel SCUBA 450um data, in combination with existing radio and ISO data at longer and shorter wavelengths to set constraints on the redshift of each source. We also analyse new and existing optical and near-infrared imaging of our SCUBA survey fields to select potential identifications consistent with these constraints. Our derived SED-based redshift constraints, and the lack of statistically significant associations with even moderately bright galaxies allow us to conclude that all 19 sources lie at z > 1, and at least half of them apparently lie at z > 2.

The mid infrared polarisation properties of four Ultraluminous Infrared Galaxies (ULIRGs) have been investigated by broad band filter observations with the ISOCAM instrument on board the Infrared Space Observatory (ISO). The wavelength region from 5 to 18 mic. was selected where the emission from the putative torus peaks. We report detection of polarisation for all ULIRGs studied. The fractional polarisation ranges from about 3% up to 8%. The highest polarisation is recorded in Mrk231 which has a clear AGN signature, whereas the lowest is for Arp220, which is generally thought to be powered predominantly by star formation. We discuss the various mechanisms that could give rise to the polarisation and conclude that the most likely interpretation is that it is due to magnetically aligned elongated dust grains. This is the same mechanism believed to be operating in a number of galactic sources. The position angle of polarisation could give the projected magnetic field direction and therefore constrain models for the formation of the tori.

The IRAM Plateau de Bure mm interferometer and deep K-band imaging have been used to identify the brightest submm source detected in the Lockman field of the UK 8mJy SCUBA survey. The near infrared counterpart is an extended (20-30kpc), clumpy, and extremely red object. The spectral energy distribution suggests it to be a dusty star forming object at a redshift of about 3 (2-4). Its star formation rate and near-infrared properties are consistent with Lockman850.1 being a massive elliptical in formation.

We report the discovery of the first hyperluminous infrared galaxy (HyLIG) in the course of the European Large Area ISO Survey (ELAIS). This object has been detected by ISO at 6.7, 15, and 90 microns, and is found to be a broad-line, radio-quiet quasar at a redshift: z = 1.099. From a detailed multi-component model fit of the spectral energy distribution, we derive a total infrared luminosity: L_IR (1-1000 microns) ~ 1.0 x 10E13 h_65^-2 L_sun (q_0 = 0.5), and discuss the possible existence of a starburst contributing to the far-IR output. Observations to date present no evidence for lens magnification. This galaxy is one of the very few HyLIGs with an X-ray detection. On the basis of its soft X-ray properties, we suggest that this broad-line object may be the face-on analogue of narrow-line, Seyfert-like HyLIGs.

In the context of the unified model of Seyfert galaxies, we use observations from the literature and a radiative transfer model to investigate the near-IR to mm emission produced by the presumed torus in the Circinus galaxy, from 2microns to 1.3mm. From the infrared SED modelling, we find that the total luminosity (L_{IR}) in this wavelength range, consists of similar contributions from the torus and starburst with a ratio of nuclear luminosity to starburst luminosity (L_{NUC}/L_{SB}) approx. 0.8. By using a similar torus model to that of NGC1068, but without the conical dust, we find an upper limit to the outer torus radius of approx. 12pc with a best fit of approx. 2pc. The upper limit torus size estimated from the radiative transfer modelling is consistent with the 16pc torus radius estimated from near-IR imaging polarimetry of Circinus.

We present BeppoSAX observations of the Southern S1 region in the European Large Area ISO Survey (ELAIS). These observations cover an area of ~1.7 sq. deg. and reach an on-axis (~0.7 sq. deg) 2-10 keV (HX) sensitivity of ~10E-13 cgs. This is the first HX analysis of an ISOCAM survey. We detect 9 sources with a signal to noise ratio SNR_HX>3, 4 additional sources with a 1.3 to 10 keV (T) SNR_T>3 and 2 additional sources which appear to be associated with QSOs with SNR_T>2.9. The number densities of the SNR_HX>3 sources are consistent with the ASCA and BeppoSAX logN-logS functions. Six BeppoSAX sources have reliable ISOCAM 15 micron counterparts within ~60 arcsec. All these ISOCAM sources have optical counterparts of R<20 mags. Five of these sources have been previously optically classified giving 4 QSOs and 1 BALQSO at z=2.2. The remaining unclassified source has X-ray and photometric properties consistent with that of a nearby Seyfert galaxy. One further HX source has a 2.6sigma ISOCAM counterpart associated with a galaxy at z=0.325. If this ISOCAM source is real, the HX/MIR properties suggest either an unusual QSO or a cD cluster galaxy. We have constructed MIR and HX spectral energy distributions to compute the expected HX/MIR ratios for these classes of objects up to z=3.2 and assess the HX/MIR survey depth. The BALQSO has an observed X-ray softness ratio and HX/MIR flux ratio similar to QSOs but different to those found for low redshift BALQSOs. This difference can be explained in terms of absorption, and suggests that high redshift BALQSOs should be comparatively easy to detect in the HX band, allowing their true fraction in the high redshift QSO population to be determined (ABRIDGED).

JHKL'M (1-5 micron) imaging of a sample of Seyfert (Sy) 2 galaxies is presented. We confirm that the 1-2.2 micron nuclear continuum of some Sy 2s is dominated by stellar emission, whereas the continuum emission at longer wavelengths (3-5 micron) is almost entirely non-stellar in origin. The IR non-stellar spectral energy distributions (SED) (up to 15 micron) of Sy galaxies show a variety of shapes, and they are well reproduced with the tapered disk models of Efstathiou & Rowan-Robinson (1995). We have used two models, one including an optically thin cone component and a coneless model. Although our modelling of the SEDs does not allow us to favor either model to account for all the observed SEDs, we find that the viewing angle towards the central source is well constrained by both models. We have also investigated non-stellar color-color diagrams. The colors of the Sy galaxies with viewing angles theta < 30 degree are better reproduced with the cone model. These diagrams provide a good means to separate Sy 2s with moderate obscuration (A_V < 20 mag, from hard X-ray observations) from those with high obscuration. The ground-based 4.8 micron and ISO 9.6 micron luminosities are well correlated with the hard X-ray luminosities of Sy 1s and 2s. These continuum emissions appear as a good indicator of the AGN luminosity, at least in the cases of hard X-ray Compton-thin Sy galaxies. We finally stress the finding that some Compton thick galaxies show bright non-stellar emission at 5 micron. This suggests that the near-IR emission in Sys is produced in an extended component illuminated by the central source, that is more visible from all viewing angles. We discuss possible implications of mid-IR surveys for the search of counterparts of highly obscured hard X-ray sources. (Abridged)

We present multiwaveband photometric and optical spectropolarimetric observations of the R=15.9 narrow emission line galaxy R117_A which lies on the edge of the error circle of the ROSAT X-ray source R117 (from McHardy et al 1998). The overall spectral energy distribution of the galaxy is well modelled by a combination of a normal spiral galaxy and a moderate-strength burst of star formation. The far infra-red and radio emission is extended along the major axis of the galaxy, indicating an extended starburst. On positional grounds, the galaxy is a good candidate for the identification of R117 and the observed X-ray flux is very close to what would be expected from a starburst of the observed far infra-red and radio fluxes. Although an obscured high redshift QSO cannot be entirely ruled out as contributing some fraction of the X-ray flux, we find no candidates to K=20.8 within the X-ray errorbox and so conclude that R117_A is responsible for a large fraction, if not all, of the X-ray emission from R117. Searches for indicators of an obscured AGN in R117_A have so far proven negative; deep spectropolarimetric observations show no signs of broad lines to a limit of one per cent and, for the observed far infra-red and radio emission, we would expect a ten times greater X-ray flux if the overall emission were powered by an AGN. We therefore conclude that the X-ray emission from R117 is dominated by starburst emission from the galaxy R117_A.

We present the luminosity function of 90um selected galaxies from the European Large Area ISO Survey (ELAIS), extending to z=0.3. Their luminosities are in the range 10^9 < h_65^-2 L/Lsun < 10^12, i.e. non-ultraluminous. From our sample of 37 reliably detected galaxies in the ELAIS S1 region from the Efstathiou et al. (2000) S_90 >= 100mJy database, we found optical, 15um or 1.4GHz identifications for 24 (65%). We have obtained 2dF and UK Schmidt FLAIR spectroscopy of 89% of IDs to rigid multivariate flux limits. We construct a luminosity function assuming (a) our spectroscopic subset is an unbiased sparse sample, and (b) there are no galaxies which would not be represented in our spectroscopic sample at {\it any} redshift. We argue that we can be confident of both assumptions. We find the luminosity function is well-described by the local 100um luminosity function of Rowan-Robinson, Helou & Walker (1987). {\it Assuming} this local normalisation, we derive luminosity evolution of (1+z)^{2.45\pm0.85} (95% confidence). We argue that star formation dominates the bolometric luminosities of these galaxies and we derive comoving star formation rates in broad agreement with the Flores et al. (1999) and Rowan-Robinson et al. (1997) mid-IR-based estimates.

We present results and source counts at 90micron extracted from the Preliminary Analysis of the European Large Area ISO Survey (ELAIS). The survey covered about 11.6 square degrees of the sky in four main areas and was carried out with the PHOT instrument onboard the Infrared Space Observatory (ISO). The survey is at least an order of magnitude deeper than the IRAS 100micron survey and is expected to provide constraints on the formation and evolution of galaxies. The majority of the detected sources are associated with galaxies on optical images. In some cases the optical associations are interacting pairs or small groups of galaxies suggesting the sample may include a significant fraction of luminous infrared galaxies. The source counts extracted from a reliable subset of the detected sources are in agreement with strongly evolving models of the starburst galaxy population.

We use an 850-micron SCUBA map of the HDF to study the dust properties of optically-selected starburst galaxies at high redshift. Although the 850-micron map does not allow detection of the majority of individual sources, we show that the galaxies with the highest UV star-formation rates are detected statistically, with a flux density of about S_{850}=0.2 mJy for an apparent UV star-formation rate of 1 h^{-2} msun yr^{-1}. This level of submillimetre output indicates that the total star-forming activity is on average a factor of approximately 6 times larger than the rate inferred from the UV output of these galaxies. The general population of optical starbursts is then predicted to contribute at least 25 per cent of the 850-micron background. We carry out a power-spectrum analysis of the map, which yields some evidence for angular clustering of the background source population, but at a level lower than that seen in Lyman-break galaxies. Together with other lines of argument, particularly from the NICMOS HDF data, this suggests that the 850-micron background originates over an extremely wide range of redshifts -- perhaps 1 < z < 6.

We present some recent results from the European Large Area ISO Survey (ELAIS). This survey was the largest non-serendipitous ISO field survey. A preliminary reduction has recently been completed and catalogues of sources released to the community. Early results show strongly evolving source counts. A comprehensive identification programme is underway and a number of extremely luminous objects have already been discovered. This survey provides an exciting legacy from the ISO mission and (amongst many goals) will allow us to provide important constraints on the obscured star-formation history of the Universe.

We present preliminary source counts at 6.7um and 15um from the Preliminary Analysis of the European Large Area ISO survey, with limiting flux densities of \~2mJy at 15um & ~1mJy at 6.7um. We separate the stellar contribution from the extragalactic using identifications with APM sources made with the likelihood ratio technique. We quantify the completeness & reliability of our source extraction using (a) repeated observations over small areas, (b) cross-IDs with stars of known spectral type, (c) detections of the PSF wings around bright sources, (d) comparison with independent algorithms. Flux calibration at 15um was performed using stellar IDs; the calibration does not agree with the pre-flight estimates, probably due to effects of detector hysteresis and photometric aperture correction. The 6.7um extragalactic counts are broadly reproduced in the Pearson & Rowan-Robinson model, but the Franceschini et al. (1997) model underpredicts the observed source density by ~0.5-1 dex, though the photometry at 6.7um is still preliminary. At 15um the extragalactic counts are in excellent agreement with the predictions of the Pearson & Rowan-Robinson (1996), Franceschini et al. (1994), Guiderdoni et al. (1997) and the evolving models of Xu et al. (1998), over 7 orders of magnitude in 15um flux density. The counts agree with other estimates from the ISOCAM instrument at overlapping flux densities (Elbaz et al. 1999), provided a consistent flux calibration is used. Luminosity evolution at a rate of (1+z)^3, incorporating mid-IR spectral features, provides a better fit to the 15um differential counts than (1+z)^4 density evolution. No-evolution models are excluded, and implying that below around 10mJy at 15um the source counts become dominated by an evolving cosmological population of dust-shrouded starbursts and/or active galaxies.

We present ISOPHOT_S, ISOSWS and 8 to 13 micron ground based observations of Centaurus A that show prominent PAH and silicate features. These and other data are used to construct a model for the infrared continuum. We find that in a nuclear sized aperture (~4 arcsec, ~60 pc) the SED is characteristic of emission from a starburst and AGN torus; in larger apertures an additional component of cirrus emission is required. Based on our model, the torus diameter is estimated to be 3.6 pc and the best fitting inclination angle of the torus is 45 degrees. This result has implications for the detectability of tori in low power AGN and in particular for the use of the IRAS 60/25 micron flux ratio as an indicator of the torus inclination.

We present illustrative models for the UV to millimeter emission of starburst galaxies which are treated as an ensemble of optically thick giant molecular clouds (GMCs) centrally illuminated by recently formed stars. The models follow the evolution of the GMCs due to the ionization-induced expansion of the HII regions and the evolution of the stellar population within the GMC according to the Bruzual & Charlot stellar population synthesis models. The effect of transiently heated dust grains/PAHs to the radiative transfer, as well as multiple scattering, is taken into account. The expansion of the HII regions and the formation of a narrow neutral shell naturally explains why the emission from PAHs dominates over that from hot dust in the near to mid-IR, an emerging characteristic of the infrared spectra of starburst galaxies. The models allow us to relate the observed properties of a galaxy to its age and star formation history. We find that exponentially decaying 10^7-10^8 yrs old bursts can explain the IRAS colours of starburst galaxies. The models are also shown to account satisfactorily for the multiwavelength data on the prototypical starburst galaxy M82 and NGC6090, a starburst galaxy recently observed by ISO. In M82 we find evidence for two bursts separated by 10^7yrs. In NGC6090 we find that at least part of the far-IR excess may be due to the age of the burst (6.4 x 10^7yrs). We also make predictions about the evolution of the luminosity of starbursts at different wavelengths which indicate that far-IR surveys may preferentially detect older starbursts than mid-IR surveys.

The European Large Area ISO Survey (ELAIS) has surveyed 12 square degrees of the sky at 15 and 90 microns, and subsets of this area at 6.75 and 175 microns, using the Infrared Space Observatory (ISO). This project was the largest single open time programme executed by ISO, taking 375 hours of data. A preliminary catalogue of more than 1000 galaxies has been produced. In this talk we describe the goals of the project, describe the follow-up programmes that are in progress, and present some first scientific results including a provisional number count analysis at 15 and 90 microns.

The advent of sensitive sub-mm array cameras now allows a proper census of dust-enshrouded massive star-formation in very distant galaxies, previously hidden activity to which even the faintest optical images are insensitive. We present the deepest sub-mm survey of the sky to date, taken with the SCUBA camera on the James Clerk Maxwell Telescope and centred on the Hubble Deep Field. The high source density found in this image implies that the survey is confusion-limited below a flux density of 2 mJy. However, within the central 80 arcsec radius independent analyses yield 5 reproducible sources with S(850um) > 2 mJy which simulations indicate can be ascribed to individual galaxies. We give positions and flux densities for these, and furthermore show using multi-frequency photometric data that the brightest sources in our map lie at redshifts z~3. These results lead to integral source counts which are completely inconsistent with a no-evolution model, and imply that massive star-formation activity continues at redshifts > 2. The combined brightness of the 5 most secure sources in our map is sufficient to account for 30 - 50% of the previously unresolved sub-mm background, and we estimate statistically that the entire background is resolved at about the 0.3 mJy level. Finally we discuss possible optical identifications and redshift estimates for the brightest sources. One source appears to be associated with an extreme starburst galaxy at z~1, whilst the remaining four appear to lie in the redshift range 2 < z < 4. This implies a star-formation density over this redshift range that is at least five times higher than that inferred from the ultraviolet output of HDF galaxies.