We present details of the Automated Radio Telescope Imaging Pipeline (ARTIP) and results of a sensitive blind search for HI and OH absorbers at $z<0.4$ and $z<0.7$, respectively. ARTIP is written in Python 3.6, extensively uses the Common Astronomy Software Application (CASA) tools and tasks, and is designed to enable the geographically-distributed MeerKAT Absorption Line Survey (MALS) team to collaboratively process large volumes of radio interferometric data. We apply it to the first MALS dataset obtained using the 64-dish MeerKAT radio telescope and 32K channel mode of the correlator. With merely 40 minutes on target, we present the most sensitive spectrum of PKS1830-211 ever obtained and characterize the known HI ($z=0.19$) and OH ($z=0.89$) absorbers. We further demonstrate ARTIP's capabilities to handle realistic observing scenarios by applying it to a sample of 72 bright radio sources observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) to blindly search for HI and OH absorbers. We estimate the numbers of HI and OH absorbers per unit redshift to be $n_{21}(z\sim0.18)<$0.14 and $n_{\rm OH}(z\sim0.40)<$0.12, respectively, and constrain the cold gas covering factor of galaxies at large impact parameters (50 kpc $<\rho<$ 150 kpc) to be less than 0.022. Due to the small redshift path, $\Delta z\sim$13 for HI with column density$>5.4\times10^{19}$ cm$^{-2}$, the survey has probed only the outskirts of star-forming galaxies at $\rho>30$ kpc. MALS with the expected $\Delta z\sim10^{3-4}$ will overcome this limitation and provide stringent constraints on the cold gas fraction of galaxies in diverse environments over $0<z<1.5$.

Precise measurement of the foreground synchrotron emission, which contaminates the faint polarized cosmic microwave background radiation (CMB), is a major challenge for the next-generation of CMB experiments. To address this, dedicated foreground measurement experiments are being undertaken at radio frequencies between 2 and 40 GHz. Foreground polarized synchrotron emission measurements are particularly challenging, primarily due to the complicated frequency dependence in the presence of Faraday rotation, and are best recovered through broad fractional-bandwidth polarization measurements at frequencies $\lesssim 5$ GHz. A unique opportunity for measuring the foreground polarized synchrotron emission will be provided by the 15-m SKA-MPG telescope operating in the frequency range 1.7 to 3.5~GHz (S-Band). Here, we present the scope of a Southern sky survey in S-Band at 1 degree angular resolution and explore its added advantage for application of powerful techniques, such as, Stokes $Q$, $U$ fitting and RM-synthesis. A full Southern-sky polarization survey with this telescope, when combined with other on-going efforts at slightly higher frequencies, will provide an excellent frequency coverage for modeling and extrapolating the foreground polarized synchrotron emission to CMB frequencies ($\gtrsim80$~GHz) with rms brightness temperature better than 10~nK per 1 degree$^2$. We find that this survey will be crucial for understanding the effects of Faraday depolarization, especially in low Galactic latitude regions. This will allow better foreground cleaning and thus will contribute significantly in further improving component separation analyses and increase usable sky area for cosmological analysis of the \textit{Planck} data, and the \textit{LiteBIRD} mission in the future.

We present the first data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 103 nearby galaxies from the Palomar sample, observed with the eMERLIN array, as part of the LeMMINGs survey. This sample includes galaxies which are active (LINER and Seyfert) and quiescent (HII galaxies and Absorption line galaxies, ALG), which are reclassified based upon revised emission-line diagrams. We detect radio emission $\gtrsim$ 0.2 mJy for 47/103 galaxies (22/34 for LINERS, 4/4 for Seyferts, 16/51 for HII galaxies and 5/14 for ALGs) with radio sizes typically of $\lesssim$100 pc. We identify the radio core position within the radio structures for 41 sources. Half of the sample shows jetted morphologies. The remaining half shows single radio cores or complex morphologies. LINERs show radio structures more core-brightened than Seyferts. Radio luminosities of the sample range from 10$^{32}$ to 10$^{40}$ erg s$^{-1}$: LINERs and HII galaxies show the highest and the lowest radio powers respectively, while ALGs and Seyferts have intermediate luminosities. We find that radio core luminosities correlate with black hole (BH) mass down to $\sim$10$^{7}$ M$_{\odot}$, but a break emerges at lower masses. Using [O III] line luminosity as a proxy for the accretion luminosity, active nuclei and jetted HII galaxies follow an optical fundamental plane of BH activity, suggesting a common disc-jet relationship. In conclusion, LINER nuclei are the scaled-down version of FR I radio galaxies; Seyferts show less collimated jets; HII galaxies may host weak active BHs and/or nuclear star-forming cores; and recurrent BH activity may account for ALG properties.

MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities.

Deep galaxy surveys have revealed that the global star formation rate (SFR) density in the Universe peaks at 1 < z < 2 and sharply declines towards z = 0. But a clear picture of the underlying processes, in particular the evolution of cold atomic (~100 K) and molecular gas phases, that drive such a strong evolution is yet to emerge. MALS is designed to use MeerKAT's L- and UHF-band receivers to carry out the most sensitive (N(HI)>10$^{19}$ cm$^{-2}$) dust-unbiased search of intervening HI 21-cm and OH 18-cm absorption lines at 0 < z < 2. This will provide reliable measurements of the evolution of cold atomic and molecular gas cross-sections of galaxies, and unravel the processes driving the steep evolution in the SFR density. The large sample of HI and OH absorbers obtained from the survey will (i) lead to tightest constraints on the fundamental constants of physics, and (ii) be ideally suited to probe the evolution of magnetic fields in disks of galaxies via Zeeman Splitting or Rotation Measure synthesis. The survey will also provide an unbiased census of HI and OH absorbers, i.e. cold gas associated with powerful AGNs (>10$^{24}$ W Hz$^{-1}$) at 0 < z < 2, and will simultaneously deliver a blind HI and OH emission line survey, and radio continuum survey. Here, we describe the MALS survey design, observing plan and the science issues to be addressed under various science themes.

The galaxy NGC 4418 harbours a compact ($<20$ pc) core with a very high bolometric luminosity ($\sim10^{11}$L$_\odot$). As most of the galaxy's energy output comes from this small region, it is of interest to determine what fuels this intense activity. An interaction with VV 655 has been proposed, where gas aquired by NGC 4418 could trigger intense star formation and/or black hole accretion in the centre. We aim to constrain the interaction hypothesis by studying neutral hydrogen structures around the two galaxies. We present observations at 1.4 GHz with the Very Large Array of radio continuum as well as emission and absorption from atomic hydrogen. Gaussian distributions are fitted to observed HI emission and absorption spectra. An atomic HI bridge is seen in emission, connecting NGC 4418 to VV 655. While NGC 4418 is bright in continuum emission and seen in HI absorption, VV 655 is barely detected in the continuum but show bright HI emission (M$_\mathrm{HI}\sim10^9$ M$_\odot$). We estimate SFRs from 1.4 GHz of 3.2 M$_\odot$ yr$^{-1}$ and 0.13 M$_\odot$ yr$^{-1}$ for NGC 4418 and VV 655 respectively. Systemic HI velocities of 2202$\pm$20 km s$^{-1}$ (emission) and 2105.4$\pm$10 km s$^{-1}$ (absorption) are measured for VV 655 and NGC 4418 respectively. Redshifted HI absorption is seen towards NGC 4418, suggesting gas infall. Blueshifted HI-emission is seen north-west of NGC 4418, which we interpret as a continuation of the outflow previously discussed by Sakamoto et al. (2013). The morphology and velocity structure seen in HI is consistent with an interaction scenario, where gas was transferred from VV 655 to NGC 4418, and may fuel the activity in the centre.

The formaldehyde MegaMaser emission has been mapped for the three host galaxies IC\,860. IRAS\,15107$+$0724, and Arp\,220. Elongated emission components are found at the nuclear centres of all galaxies with an extent ranging between 30 to 100 pc. These components are superposed on the peaks of the nuclear continuum. Additional isolated emission components are found superposed in the outskirts of the radio continuum structure. The brightness temperatures of the detected features ranges from 0.6 to 13.4 $\times 10^{4}$ K, which confirms their masering nature. The masering scenario is interpreted as amplification of the radio continuum by foreground molecular gas that is pumped by far-infrared radiation fields in these starburst environments of the host galaxies.

The QueSERA Study Group (QSG) have been tasked by the RadioNet Board to produce a White Paper on the future organisation and coordination of radio astronomy in Europe. This White Paper describes the options discussed by the QSG, and our conclusions on how to move forward. We propose, that as a first step, RadioNet-work, be established as an entity that persists between EC contracts, and that takes responsibility for preparing or coordinating responses to EC opportunities specific to the field of radio astronomy research infrastructures. RadioNet-work should provide a safety net that ensures that cooperation and collaboration between the various radio astronomy partners in Europe is maintained with or without EC funding.

The Square Kilometer Array (SKA) has the potential to produce galaxy redshift surveys which will be competitive with other state of the art cosmological experiments in the next decade. In this chapter we summarise what capabilities the first and the second phases of the SKA will be able to achieve in its current state of design. We summarise the different cosmological experiments which are outlined in further detail in other chapters of this Science Book. The SKA will be able to produce competitive Baryonic Oscillation (BAOs) measurements in both its phases. The first phase of the SKA will provide similar measurements as optical and IR experiments with completely different systematic effects whereas the second phase being transformational in terms of its statistical power. The SKA will produce very accurate Redshift Space Distortions (RSD) measurements, being superior to other experiments at lower redshifts, due to the large number of galaxies. Cross correlations of the galaxy redshift data from the SKA with radio continuum surveys and optical surveys will provide extremely good calibration of photometric redshifts as well as extremely good bounds on modifications of gravity. Basing on a Principle Component Analysis (PCA) approach, we find that the SKA will be able to provide competitive constraints on dark energy and modified gravity models. Due to the large area covered the SKA it will be a transformational experiment in measuring physics from the largest scales such as non-Gaussian signals from $\textrm{f}_{\textrm{nl}}$. Finally, the SKA might produce the first real time measurement of the redshift drift. The SKA will be a transformational machine for cosmology as it grows from an early Phase 1 to its full power.

In recent years cosmology has undergone a revolution, with precise measurements of the microwave background radiation, large galaxy redshift surveys, and the discovery of the recent accelerated expansion of the Universe using observations of distant supernovae. In this light, the SKA enables us to do an ultimate test in cosmology by measuring the expansion rate of the Universe in real time. This can be done by a rather simple experiment of observing the neutral hydrogen (HI) signal of galaxies at two different epochs. The signal will encounter a change in frequency imprinted as the Universe expands over time and thus monitoring the drift in frequencies will provide a real time measure of the cosmic acceleration. Over a period of 12 years one would expected a frequency shift of the order of 0.1 Hz assuming a standard Lambda-CDM cosmology. Based on the sensitivity estimates of the SKA and the number counts of the expected HI galaxies, it is shown that the number counts are sufficiently high to compensate for the observational uncertainties of the measurements and hence allow a statistical detection of the frequency shift. [abstract abridged]

Galaxies and supermassive black holes (SMBHs) are believed to evolve through a process of hierarchical merging and accretion. Through this paradigm, multiple SMBH systems are expected to be relatively common in the Universe. However, to date there are poor observational constraints on multiple SMBHs systems with separations comparable to a SMBH gravitational sphere of influence (<< 1 kpc). In this chapter, we discuss how deep continuum observations with the SKA will make leading contributions towards understanding how multiple black hole systems impact galaxy evolution. In addition, these observations will provide constraints on and an understanding of stochastic gravitational wave background detections in the pulsar timing array sensitivity band (nanoHz -microHz). We also discuss how targets for pointed gravitational wave experiments (that cannot be resolved by VLBI) could potentially be found using the large-scale radio-jet morphology, which can be modulated by the presence of a close-pair binary SMBH system. The combination of direct imaging at high angular resolution; low-surface brightness radio-jet tracers; and pulsar timing arrays will allow the SKA to trace black hole binary evolution from separations of a galaxy virial radius down to the sub-parsec level. This large dynamic range in binary SMBH separation will ensure that the SKA plays a leading role in this observational frontier.

We study intergalactic medium (IGM) heating in a sample of five fossil galaxy groups by using their radio properties at 610 MHz and 1.4 GHz. The power by radio jets introducing mechanical heating for the sampled objects is not sufficient enough to suppress the cooling flow. Therefore, we discussed shock-, vortex heating, and conduction as alternative heating processes. Further, the 1.4 GHz and 610 MHz radio luminosities of fossil groups are compared to a sample of normal galaxy groups of the same radio brightest (BGGs), stellar mass, and total group stellar mass, quantified using the $K$-band luminosity. It appears that the fossil BGGs are under luminous at 1.4 GHz and 610 MHz for a given BGG stellar mass and luminosity, in comparison to a general population of the groups. In addition, we explore how the bolometric radio luminosity of fossil sample depends on clusters and groups characteristics. Using the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS) as a control sample we found that the large-scale behaviours of fossil galaxy groups are consistent with their relaxed and virialised nature.

Galaxies are believed to evolve through merging, which should lead to multiple supermassive black holes in some. There are four known triple black hole systems, with the closest pair being 2.4 kiloparsecs apart (the third component is more distant at 3 kiloparsecs), which is far from the gravitational sphere of influence of a black hole with mass $\sim$10$^9$ M$_\odot$ (about 100 parsecs). Previous searches for compact black hole systems concluded that they were rare, with the tightest binary system having a separation of 7 parsecs. Here we report observations of a triple black hole system at redshift z=0.39, with the closest pair separated by $\sim$140 parsecs. The presence of the tight pair is imprinted onto the properties of the large-scale radio jets, as a rotationally-symmetric helical modulation, which provides a useful way to search for other tight pairs without needing extremely high resolution observations. As we found this tight pair after searching only six galaxies, we conclude that tight pairs are more common than hitherto believed, which is an important observational constraint for low-frequency gravitational wave experiments.

The galaxy NGC 4418 contains one of the most compact obscured nuclei within a luminous infrared galaxy (LIRG) in the nearby Universe. This nucleus contains a rich molecular gas environment and an unusually high ratio of infrared to radio luminosity (q-factor). The compact nucleus is powered by either a compact starburst or an active galactic nucleus (AGN). The aim of this study is to constrain the nature of the nuclear region (starburst or AGN) within NGC 4418 via very-high-resolution radio imaging. Archival data from radio observations using the EVN and MERLIN interferometers are imaged. Sizes and flux densities are obtained by fitting Gaussian intensity distributions to the image. The average spectral index of the compact radio emission is estimated from measurements at 1.4 GHz and 5.0 GHz. The nuclear structure of NGC 4418 visible with EVN and MERLIN consists of eight compact (<49 mas i.e. <8 pc) features spread within a region of 250 mas, i.e. 41 pc. We derive an inverted spectral index $\alpha\ge0.7$ ($S_\nu\propto\nu^{\alpha}$) for the compact radio emission. Brightness temperatures $>10^{4.8}$ K indicate that these compact features cannot be HII-regions. The complex morphology and inverted spectrum of the eight detected compact features is evidence against the hypothesis that an AGN alone is powering the nucleus of NGC 4418. The compact features could be super star clusters (SSCs) with intense star formation, and their associated free-free absorption could then naturally explain both their inverted radio spectrum and the low radio to IR ratio of the nucleus. The required star formation area density is extreme, however, and close to the limit of what can be observed in a well-mixed thermal/non-thermal plasma produced by star-formation, and is also close to the limit of what can be physically sustained.

We report 1.7 GHz Very Long Baseline Interferometry (VLBI) observations of IRAS F10214+4724, a lensed z=2.3 obscured quasar with prodigious star formation. We detect what we argue to be the obscured active nucleus with an effective angular resolution of < 50 pc at z = 2.3 . The S_{1.7} = 210 micro-Jy (9-\sigma) detection of this unresolved source is located within the HST rest-frame ultraviolet/optical arc, however, >~100 mas northward of the arc centre of curvature. This leads to a source plane inversion that places the European VLBI Network detection to within milli-arcseconds of the modelled cusp caustic, resulting in a very large magnification (\mu ~70), over an order of magnitude larger than the CO (1-0) derived magnification of a spatially resolved JVLA map, using the same lens model. We estimate the quasar bolometric luminosity from a number of independent techniques and with our X-ray modelling find evidence that the AGN may be close to Compton-thick, with an intrinsic bolometric luminosity log(L_{bol,QSO} / L_sun) = 11.34 +- 0.27 dex. We make the first black hole mass estimate of IRAS F10214+4724 and find log(M_{BH}/M_sun) = 8.36 +- 0.56 which suggests a low black hole accretion rate (\lambda = \dot{M} / \dot{M}_{Edd} ~ 3\pm^7_2 percent). We find evidence for a M_{BH}/M_{spheroid} ratio that is 1-2 orders of magnitude larger than that of submillimetre galaxies (SMGs) at z~2. At face value, this suggests IRAS F10214+4724 has undergone a different evolutionary path compared to SMGs at the same epoch. A primary result of this work is the demonstration that emission regions of differing size and position can undergo significantly different magnification boosts (> 1 dex) and therefore distort our view of high-redshift, gravitationally lensed galaxies.

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

We describe a 325-MHz survey, undertaken with the Giant Metrewave Radio Telescope (GMRT), which covers a large part of the three equatorial fields at 9, 12 and 14.5 h of right ascension from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) in the area also covered by the Galaxy And Mass Assembly survey (GAMA). The full dataset, after some observed pointings were removed during the data reduction process, comprises 212 GMRT pointings covering ~90 deg^2 of sky. We have imaged and catalogued the data using a pipeline that automates the process of flagging, calibration, self-calibration and source detection for each of the survey pointings. The resulting images have resolutions of between 14 and 24 arcsec and minimum rms noise (away from bright sources) of ~1 mJy/beam, and the catalogue contains 5263 sources brighter than 5 sigma. We investigate the spectral indices of GMRT sources which are also detected at 1.4 GHz and find them to agree broadly with previously published results; there is no evidence for any flattening of the radio spectral index below S_1.4=10 mJy. This work adds to the large amount of available optical and infrared data in the H-ATLAS equatorial fields and will facilitate further study of the low-frequency radio properties of star formation and AGN activity in galaxies out to z~1.

The aim of this study is to constrain the dynamics, structure and feeding of the compact nucleous of NGC4418, and to reveal the nature of the main hidden power source: starburst or AGN. We obtained high spatial resolution observations of NGC4418 at 1.4 and 5 GHz with MERLIN, and at 230 and 270 GHz with the SMA very extended configuration. We use the continuum morphology and flux density to estimate the size of the emitting region, the star formation rate and the dust temperature. Emission lines are used to study the kinematics through position-velocity diagrams. Molecular emission is studied with population diagrams and by fitting an LTE synthetic spectrum. We detect bright 1mm line emission from CO, HC3N, HNC and C34S, and 1.4 GHz absorption from HI. The CO 2-1 emission and HI absorption can be fit by two velocity components at 2090 and 2180 km s-1. We detect vibrationally excited HC3N and HNC, with Tvib 300K. Molecular excitation is consistent with a layered temperature structure, with three main components at 80, 160 and 300 K. For the hot component we estimate a source size of less than 5 pc. The nuclear molecular gas surface density of 1e4 Msun pc-2 is extremely high, and similar to that found in the ultra-luminous infrared galaxy (ULIRG) Arp220. Our observations confirm the the presence of a molecular and atomic in-flow, previously suggested by Herschel observations, which is feeding the activity in the center of NGC4418. Molecular excitation confirms the presence of a very compact, hot dusty core. If a starburst is responsible for the observed IR flux, this has to be at least as extreme as the one in Arp220, with an age of 3-10 Myr and a star formation rate >10 Msun yr-1. If an AGN is present, it must be extremely Compton-thick.

We present first results from electronic Multi-Element Remotely Linked Interferometer Network (e-MERLIN) and electronic European VLBI Network (e-EVN) observations of a small sample of ultra-steep spectrum radio sources, defined as those sources with a spectral index alpha < -1.4 between 74 MHz and 325 MHz, which are unresolved on arcsecond scales. Such sources are currently poorly understood and a number of theories as to their origin have been proposed in the literature. The new observations described here have resulted in the first detection of two of these sources at milliarcsecond scales and show that a significant fraction of ultra-steep spectrum sources may have compact structures which can only be studied at the high resolution available with very long baseline interferometry (VLBI).

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) provides an unprecedented opportunity to search for blazars at sub-mm wavelengths. We cross-matched the FIRST radio source catalogue with the 11655 sources brighter than 35 mJy at 500{\mu}m in the \sim 135 square degrees of the sky covered by the H-ATLAS equatorial fields at 9 h and 15 h, plus half of the field at 12 h. We found that 379 of the H-ATLAS sources have a FIRST counterpart within 10 arcsec, including 8 catalogued blazars (plus one known blazar that was found at the edge of one the H-ATLAS maps). To search for additional blazar candidates we have devised new diagnostic diagrams and found that known blazars occupy a region of the log(S500{\mu}m/S350{\mu}m) vs. log(S500{\mu}m/S1.4GHz) plane separated from that of the other sub-mm sources with radio counterparts. Using this diagnostic we have selected 12 further candidates that turn out to be scattered in the (r-z) vs. (u-r) plane or in the WISE colour-colour diagram proposed by Massaro et al. (2012), where known blazars are concentrated in well defined strips. This suggests that the majority of them either are not blazars or have spectral energy distributions contaminated by their host galaxies. A significant fraction of true blazars are found to be hosted by star-forming galaxies. This finding, supported by an analysis of blazars detected in Planck 545 and 857 GHz bands, is at odds with the notion that blazar hosts are passive ellipticals and indicates that the sub-mm selection is providing a novel prospect on blazar properties. Based on an inspection of the available photometric data, including the WISE all-sky survey, the unpublished VIKING survey and new radio observations, we tentatively estimate that there are 11 blazars with synchrotron flux density S500{\mu}m > 35mJy over the considered area. This result already allows us to constrain blazar evolution models.

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