Future surveys focusing on understanding the nature of dark energy (e.g., Euclid and WFIRST) will cover large fractions of the extragalactic sky in near-IR slitless spectroscopy. These surveys will detect a large number of galaxies that will have only one emission line in the covered spectral range. In order to maximize the scientific return of these missions, it is imperative that single emission lines are correctly identified. Using a supervised machine-learning approach, we classified a sample of single emission lines extracted from the WFC3 IR Spectroscopic Parallel survey (WISP), one of the closest existing analogs to future slitless surveys. Our automatic software integrates a SED fitting strategy with additional independent sources of information. We calibrated it and tested it on a "gold" sample of securely identified objects with multiple lines detected. The algorithm correctly classifies real emission lines with an accuracy of 82.6%, whereas the accuracy of the SED fitting technique alone is low (~50%) due to the limited amount of photometric data available (<=6 bands). While not specifically designed for the Euclid and WFIRST surveys, the algorithm represents an important precursor of similar algorithms to be used in these future missions.

Context:Blazars are the rarest and most powerful active galactic nuclei, playing a crucial and growing role in today multi-frequency and multi-messenger astrophysics. Current blazar catalogs, however, are incomplete and particularly depleted at low Galactic latitudes. Aims: We aim at augmenting the current blazar census to build a catalog of blazar candidates with homogeneous sky coverage that can provide candidate counterparts to unassociated gamma-ray sources, sources of high-energy neutrino emission, and ultra-high energy cosmic rays. Methods: Starting from the ALMA Calibrator Catalog we built a catalog of 1580 blazar candidates (ALMA Blazar Candidates, ABC) for which we collect multi-wavelength information. We also compared ABC sources with existing blazar catalogs. Results: The ABC catalogue fills the lack of low Galactic latitude sources in current blazar catalogues. ABC sources are significantly dimmer than known blazars in Gaia g band, and they appear bluer in SDSS and WISE colors. The majority of ABC sources (~ 90%) have optical spectra that classify them as QSO, while the remaining sources resulted galactic objects. ABC sources are similar in X-rays to known blazar, while in gamma-rays they are on average dimmer and softer, indicating a significant contribution of FSRQ sources. Making use of WISE colours, we classified 715 ABC sources as candidate gamma-ray blazar of different classes. Conclusions: We built a new catalogue of 1580 candidate blazars with a rich multi-wavelength data-set, filling the lack of low Galactic latitude sources in current blazar catalogues. This will be particularly important to identify the source population of high energy neutrinos or ultra-high energy cosmic rays. The data collected by the upcoming LSST surveys will provide a key tool to investigate the possible blazar nature of these sources.

We illustrate the extraordinary discovery potential for extragalactic astrophysics of a far-IR/submm all-sky spectroscopic survey with a 3m-class space telescope. Spectroscopy provides both a 3D view of the Universe and allows us to take full advantage of the sensitivity of present-day instrumentation, overcoming the spatial confusion that affects broadband far-IR/submm surveys. Emission lines powered by star formation will be detected in galaxies out to $z \simeq 8$. It will provide measurements of spectroscopic redshifts, SFRs, dust masses, and metal content for millions of galaxies at the peak epoch of cosmic star formation and of hundreds of them at the epoch of reionization. Many of these galaxies will be strongly lensed; the brightness amplification and stretching of their sizes will make it possible to investigate (by means of follow-up with high-resolution instruments) their internal structure and dynamics on the scales of giant molecular clouds. This will provide direct information on the physics driving the evolution. Furthermore, the arc-min resolution of the telescope at submm wavelengths is ideal for detecting the cores of galaxy proto-clusters, out to the epoch of reionization. Tens of millions of these galaxy-clusters-in-formation will be detected at $z \simeq 2$-3, with a tail out to $z \simeq 7$, and thousands of detections at 6 < z < 7. Their study will allow us to track the growth of the most massive halos well beyond what is possible with classical cluster surveys (mostly limited to $z < 1.5$-2), tracing the history of star formation in dense environments and teaching us how star formation and galaxy-cluster formation are related across all epochs. Such a survey will overcome the current lack of spectroscopic redshifts of dusty star-forming galaxies and galaxy proto-clusters, representing a quantum leap in far-IR/submm extragalactic astrophysics.

Recent estimates point to abundances of z > 4 sub-millimeter (sub-mm) galaxies far above model predictions. The matter is still debated. According to some analyses the excess may be substantially lower than initially thought and perhaps accounted for by flux boosting and source blending. However, there is no general agreement on this conclusion. An excess of z > 6 dusty galaxies has also been reported albeit with poor statistics. On the other hand, evidence of a top-heavy initial mass function (IMF) in high-z starburst galaxies has been reported in the past decades. This would translate into a higher sub-mm luminosity of dusty galaxies at fixed star formation rate, i.e., into a higher abundance of bright high-z sub-mm galaxies than expected for a universal Chabrier IMF. Exploiting our physical model for high-z proto-spheroidal galaxies, we find that part of the excess can be understood in terms of an IMF somewhat top-heavier than Chabrier. Such IMF is consistent with that recently proposed to account for the low 13C/18O abundance ratio in four dusty starburst galaxies at z = 2-3. However, extreme top-heavy IMFs are inconsistent with the sub-mm counts at z > 4.

Aims. We aim to study the far-infrared radio correlation (FIRC) at 150 MHz in the local Universe (at a median redshift z~0:05) and improve the use of the rest-frame 150-MHz luminosity, L150, as a star-formation rate (SFR) tracer, which is unaffected by dust extinction. Methods. We cross-match the 60-um selected Revised IRAS Faint Source Survey Redshift (RIFSCz) catalogue and the 150-MHz selected LOFAR value-added source catalogue in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring Field. We estimate L150 for the cross-matched sources and compare it with the total infrared (IR) luminosity, LIR, and various SFR tracers. Results. We find a tight linear correlation between log L150 and log LIR for star-forming galaxies, with a slope of 1.37. The median qIR value (defined as the logarithm of the LIR to L150 ratio) and its rms scatter of our main sample are 2.14 and 0.34, respectively. We also find that log L150 correlates tightly with the logarithm of SFR derived from three different tracers, i.e., SFR_Halpha based on the Halpha line luminosity, SFR_60 based on the rest-frame 60-um luminosity and SFR_IR based on LIR, with a scatter of 0.3 dex. Our best-fit relations between L150 and these SFR tracers are, log L150 (Lsun) = 1.35(0.06) x log SFR_Halpha (Msun/yr) + 3.20(0.06), log L150 (Lsun) = 1.31(0.05) x log SFR_60 (Msun/yr) + 3.14(0.06), and log L150 (Lsun) = 1.37(0.05) x log SFR_IR (Msun/yr) + 3.09(0.05), which show excellent agreement with each other.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

Planck, SPT and ACT surveys have clearly demonstrated that Cosmic Microwave Background (CMB) experiments, while optimised for cosmological measurements, have made important contributions to the field of extragalactic astrophysics in the last decade. Future CMB experiments have the potential to make even greater contributions. One example is the detection of high-z galaxies with extreme gravitational amplifications. The combination of flux boosting and of stretching of the images has allowed the investigation of the structure of galaxies at z ~3 with the astounding spatial resolution of about 60 pc. Another example is the detection of proto-clusters of dusty galaxies at high z when they may not yet possess the hot intergalactic medium allowing their detection in X-rays or via the Sunyaev-Zeldovich effect. Next generation CMB experiments, like PICO, CORE, CMB-Bharat from space and Simons Observatory and CMB-S4 from the ground, will discover several thousands of strongly lensed galaxies out to z~6 or more and of galaxy proto-clusters caught in the phase when their member galaxies where forming the bulk of their {stars. They will also detect tens of thousands of local dusty galaxies and thousands of radio sources at least up to z~5. Moreover they will measure the polarized emission of thousands of radio sources and of dusty galaxies at mm/sub-mm wavelengths.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

We present Atacama Large Millimeter/submillimiter Array (ALMA) high sensitivity ($\sigma_P \simeq 0.4\,$mJy) polarimetric observations at $97.5\,$GHz (Band 3) of a complete sample of $32$ extragalactic radio sources drawn from the faint Planck-ATCA Co-eval Observations (PACO) sample ($b<-75^\circ$, compact sources brighter than $200\,$mJy at $20\,$GHz). We achieved a detection rate of $~97\%$ at $3\,\sigma$ (only $1$ non-detection). We complement these observations with new Australia Telescope Compact Array (ATCA) data between $2.1$ and $35\,$GHz obtained within a few months and with data published in earlier papers from our collaboration. Adding the co-eval GaLactic and Extragalactic All-sky Murchison widefield array (GLEAM) survey detections between $70\,$ and $230\,$MHz for our sources, we present spectra over more than $3$ decades in frequency in total intensity and over about $1.7$ decades in polarization. The spectra of our sources are smooth over the whole frequency range, with no sign of dust emission from the host galaxy at mm wavelengths nor of a sharp high frequency decline due, for example, to electron ageing. We do however find indications of multiple emitting components and present a classification based on the number of detected components. We analyze the polarization fraction behaviour and distributions up to $97\,$GHz for different source classes. Source counts in polarization are presented at $95\,$GHz.

CMB surveys provide, for free, blindly selected samples of extragalactic radio sources at much higher frequencies than traditional radio surveys. Next-generation, ground-based CMB experiments with arcmin resolution at mm wavelengths will provide samples of thousands radio sources allowing the investigation of the evolutionary properties of blazar populations, the study of the earliest and latest stages of radio activity, the discovery of rare phenomena and of new transient sources and events. Space-borne experiments will extend to sub-mm wavelengths the determinations of the SEDs of many hundreds of blazars, in temperature and in polarization, allowing us to investigate the flow and the structure of relativistic jets close to their base, and the electron acceleration mechanisms. A real breakthrough will be achieved in the caracterization of the polarization properties. The first direct counts in polarization will be obtained, enabling a solid assessment of the extra-galactic source contamination of CMB maps and allowing us to understand structure and intensity of magnetic fields, particle densities and structures of emitting regions close to the base of the jet.

Next generation CMB experiments with arcmin resolution will, for free, lay the foundations for a real breakthrough on the study of the early evolution of galaxies and galaxy clusters, thanks to the detection of large samples of strongly gravitationally lensed galaxies and of proto-clusters of dusty galaxies up to high redshifts. This has an enormous legacy value. High resolution follow-up of strongly lensed galaxies will allow the direct investigation of their structure and kinematics up to z~6, providing direct information on physical processes driving their evolution. Follow-up of proto-clusters will allow an observational validation of the formation history of the most massive dark matter halos up to z~4, well beyond the redshift range accessible via X-ray or SZ measurements. These experiments will also allow a giant leap forward in the determination of polarization properties of extragalactic sources, and will provide a complete census of cold dust available for star formation in the local universe.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

We illustrate the extraordinary potential of the (far-IR) Origins Survey Spectrometer (OSS) on board the Origins Space Telescope (OST) to address a variety of open issues on the co-evolution of galaxies and AGNs. We present predictions for blind surveys, each of 1000 h, with different mapped areas (a shallow survey covering an area of 10 deg$^{2}$ and a deep survey of 1 deg$^{2}$) and two different concepts of the OST/OSS: with a 5.9 m telescope (Concept 2, our reference configuration) and with a 9.1 m telescope (Concept 1, previous configuration). In 1000 h, surveys with the reference concept will detect from $\sim 1.9 \times 10^{6}$ to $\sim 8.7 \times 10^{6}$ lines from $\sim 4.8 \times 10^{5}$-$2.7 \times 10^{6}$ star-forming galaxies and from $\sim 1.4 \times 10^{4}$ to $\sim 3.8 \times 10^{4}$ lines from $\sim 1.3 \times 10^{4}$-$3.5 \times 10^{4}$ AGNs. The shallow survey will detect substantially more sources than the deep one; the advantage of the latter in pushing detections to lower luminosities/higher redshifts turns out to be quite limited. The OST/OSS will reach, in the same observing time, line fluxes more than one order of magnitude fainter than the SPICA/SMI and will cover a much broader redshift range. In particular it will detect tens of thousands of galaxies at $z \geq 5$, beyond the reach of that instrument. The polycyclic aromatic hydrocarbons lines are potentially bright enough to allow the detection of hundreds of thousands of star-forming galaxies up to $z \sim 8.5$, i.e. all the way through the re-ionization epoch. The proposed surveys will allow us to explore the galaxy-AGN co-evolution up to $z\sim 5.5-6$ with very good statistics. OST Concept 1 does not offer significant advantages for the scientific goals presented here.

We present a new catalogue of ALMA observations of 3,364 bright, compact radio sources, mostly blazars, used as calibrators. These sources were observed between May 2011 and July 2018, for a total of 47,115 pointings in different bands and epochs. We have exploited the ALMA data to validate the photometry given in the new Planck Multi-frequency Catalogue of Non-thermal sources (PCNT), for which an external validation was not possible so far. We have also assessed the positional accuracy of Planck catalogues and the PCNT completeness limits, finding them to be consistent with those of the Second Planck Catalogue of Compact Sources. The ALMA continuum spectra have allowed us to extrapolate the observed radio source counts at 100 GHz to the effective frequencies of ALMA bands 4, 6, 7, 8 and 9 (145, 233, 285, 467 and 673 GHz, respectively), where direct measurements are scanty, especially at the 3 highest frequencies. The results agree with the predictions of the Tucci et al. model C2Ex, while the model C2Co is disfavoured.

ALMA products are stored in the Science Archive in the form of FITS images. It is a common idea that the FITS image headers should collect in their keywords all the information that an archive User might want to search for in order to quickly select, compare, or discard datasets. With this perspective in mind, we first present a short description of the current status of the ALMA FITS archive and images. We realized that at the moment most of the parameters that could be useful for a general User are still missing in the archived data. We then provide a CASA task generating the image header keywords that we suggest to be relevant for the scientific exploitation of the ALMA archival data. The proposed tool could be also applied to several types of interferometer data and part of it is implemented in a web interface. An example of the scientific application of the keywords is also discussed.

We present the Tiered Radio Extragalactic Continuum Simulation (T-RECS): a new simulation of the radio sky in continuum, over the 150 MHz-20 GHz range. T-RECS models two main populations of radio galaxies: Active Galactic Nuclei (AGNs) and Star-Forming Galaxies (SFGs), and corresponding sub-populations. Our model also includes polarized emission over the full frequency range, which has been characterised statistically for each population using the available information. We model the clustering properties in terms of probability distributions of hosting halo masses, and use lightcones extracted from a high-resolution cosmological simulation to determine the positions of haloes. This limits the sky area for the simulations including clustering to a 25deg2 field of view. We compare luminosity functions, number counts in total intensity and polarization, and clustering properties of our outputs to up-to-date compilations of data and find a very good agreement. We deliver a set of simulated catalogues, as well as the code to produce them, which can be used for simulating observations and predicting results from deep radio surveys with existing and forthcoming radio facilities, such as the Square Kilometre Array (SKA).

The Probe of Inflation and Cosmic Origins (PICO) is a NASA-funded study of a Probe-class mission concept. The top-level science objectives are to probe the physics of the Big Bang by measuring or constraining the energy scale of inflation, probe fundamental physics by measuring the number of light particles in the Universe and the sum of neutrino masses, to measure the reionization history of the Universe, and to understand the mechanisms driving the cosmic star formation history, and the physics of the galactic magnetic field. PICO would have multiple frequency bands between 21 and 799 GHz, and would survey the entire sky, producing maps of the polarization of the cosmic microwave background radiation, of galactic dust, of synchrotron radiation, and of various populations of point sources. Several instrument configurations, optical systems, cooling architectures, and detector and readout technologies have been and continue to be considered in the development of the mission concept. We will present a snapshot of the baseline mission concept currently under development.

We present a catalogue of ALMA flux density measurements of 754 calibrators observed between August 2012 and September 2017, for a total of 16,263 observations in different bands and epochs. The flux densities were measured reprocessing the ALMA images generated in the framework of the ALMACAL project, with a new code developed by the Italian node of the European ALMA Regional Centre. A search in the online databases yielded redshift measurements for 589 sources ($\sim$78 per cent of the total). Almost all sources are flat-spectrum, based on their low-frequency spectral index, and have properties consistent with being blazars of different types. To illustrate the properties of the sample we show the redshift and flux density distributions as well as the distributions of the number of observations of individual sources and of time spans in the source frame for sources observed in bands 3 (84$-$116 GHz) and 6 (211$-$275 GHz). As examples of the scientific investigations allowed by the catalogue we briefly discuss the variability properties of our sources in ALMA bands 3 and 6 and the frequency spectra between the effective frequencies of these bands. We find that the median variability index steadily increases with the source-frame time lag increasing from 100 to 800 days, and that the frequency spectra of BL Lacs are significantly flatter than those of flat-spectrum radio quasars. We also show the global spectral energy distributions of our sources over 17 orders of magnitude in frequency.

This lecture briefly reviews the major recent advances in radio astronomy made possible by ultra-deep surveys, reaching microJansky flux density levels. A giant step forward in many fields, including the study of the evolution of the cosmic star formation history is expected with the advent of the Square Kilometer Array (SKA).

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

We present high sensitivity ($\sigma_P \simeq 0.6\,$mJy) polarimetric observations in seven bands, from $2.1$ to $38\,$GHz, of a complete sample of $104$ compact extragalactic radio sources brighter than $200\,$mJy at $20\,$GHz. Polarization measurements in six bands, in the range $5.5-38\,$GHz, for $53$ of these objects were reported by \citet{Galluzzi2017}. We have added new measurements in the same six bands for another 51 sources and measurements at $2.1\,$GHz for the full sample of $104$ sources. Also, the previous measurements at $18$, $24$, $33$ and $38\,$GHz were re-calibrated using the updated model for the flux density absolute calibrator, PKS1934-638, not available for the earlier analysis. The observations, carried out with the Australia Telescope Compact Array (ATCA), achieved a $90\%$ detection rate (at $5\sigma$) in polarization. $89$ of our sources have a counterpart in the $72$ to $231\,$MHz GLEAM survey \citep{HurleyWalker2017}, providing an unparalleled spectral coverage of $2.7$ decades of frequency for these sources. While the total intensity data from $5.5$ to $38\,$GHz could be interpreted in terms of single component emission, a joint analysis of more extended total intensity spectra presented here, and of the polarization spectra, reveals that over $90\%$ of our sources show clear indications of at least two emission components. We interpret this as an evidence of recurrent activity. Our high sensitivity polarimetry has allowed a $5\,\sigma$ detection of the weak circular polarization for $\sim 38\%$ of the dataset, and a deeper estimate of $20\,$GHz polarization source counts than has been possible so far.

We present stacked average far-infrared spectra of a sample of 197 dusty, star-forming galaxies (DSFGs) at $0.005 < z < 4$ using close to 90% of the SPIRE Fourier Transform Spectrometer (FTS) extragalactic data archive from the Herschel Space Observatory based on 3.5 years of science operations. These spectra explore an observed-frame $\rm 447\,GHz-1568\,GHz$ ($\rm 191\,\mu m-671\,\mu m$) frequency (wavelength) range allowing us to observe the main atomic and molecular lines emitted by gas in the interstellar medium. The sample is sub-divided into five redshift bins at $0.005 < z < 0.05$, $0.05 < z < 0.2$, $0.2 < z < 0.5$, $0.8 < z <2$, and $2 < z < 4$. To study the dependence of observed spectral lines on total infrared luminosity, the sources in a subset of the redshift bins are stacked in luminosity bins. These stacked spectra are used to determine the average properties of the interstellar medium and dense molecular gas properties of DSFGs, in particular, the fine-structure line ([CII] 158 $\mu$m and [OI] 63 $\mu$m) luminosity ratios, and the line to far-IR luminosity ratios are used to model the gas density and radiation field strength in the photodissociation regions (PDRs). For the low-redshift sample, we additionally present the average spectral line energy distributions (SLED) of CO and $\rm{H_2O}$ rotational transitions and also consider PDR conditions based on a combination of [CI] 370 $\mu$m and 609 $\mu$m and $\rm CO (7-6)$ lines. For the high-z ($0.8 < z < 4$) sample PDR models suggest a molecular gas distribution in the presence of a radiation field that is at least a factor of 10$^3$ larger than the Milky-Way and with a neutral gas density of roughly 10$^3$ to 10$^5$ cm$^{-3}$. The corresponding PDR models for the low-z sample suggest a UV radiation field and gas density comparable to those at high-z.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]

The bulk of the stellar growth over cosmic time is dominated by IR luminous galaxies at cosmic noon (z=1-2), many of which harbor a hidden active galactic nucleus (AGN). We use state of the art infrared color diagnostics, combining Spitzer and Herschel observations, to separate dust-obscured AGN from dusty star forming galaxies (SFGs) in the CANDELS and COSMOS surveys. We calculate 24 micron counts of SFGs, AGN/star forming "Composites", and AGN. AGN and Composites dominate the counts above 0.8 mJy at 24 micron, and Composites form at least 25% of an IR sample even to faint detection limits. We develop methods to use the Mid-Infrared Instrument (MIRI) on JWST to identify dust-obscured AGN and Composite galaxies from z~1-2. With the sensitivity and spacing of MIRI filters, we will detect >4 times as many AGN hosts than with Spitzer/IRAC criteria. Any star formation rates based on the 7.7 micron PAH feature (likely to be applied to MIRI photometry) must be corrected for the contribution of the AGN, or the SFR will be overestimated by ~35% for cases where the AGN provides half the IR luminosity and ~50% when the AGN accounts for 90% of the luminosity. Finally, we demonstrate that our MIRI color technique can select AGN with an Eddington ratio of $\lambda_{\rm Edd}\sim0.01$ and will identify AGN hosts with a higher sSFR than X-ray techniques alone. JWST/MIRI will enable critical steps forward in identifying and understanding dust-obscured AGN and the link to their host galaxies.

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]

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

Observational investigations of the abundance of massive precursors of local galaxy clusters ("proto-clusters") allow us to test the growth of density perturbations, to constrain cosmological parameters that control it, to test the theory of non-linear collapse and how the galaxy formation takes place in dense environments. The Planck collaboration has recently published a catalogue of >~ 2000 cold extra-galactic sub-millimeter sources, i.e. with colours indicative of z >~ 2, almost all of which appear to be over-densities of star-forming galaxies. They are thus considered as proto-cluster candidates. Their number densities (or their flux densities) are far in excess of expectations from the standard scenario for the evolution of large-scale structure. Simulations based on a physically motivated galaxy evolution model show that essentially all cold peaks brighter than S_{545GHz} = 500 mJy found in Planck maps after having removed the Galactic dust emission can be interpreted as positive Poisson fluctuations of the number of high-z dusty proto-clusters within the same Planck beam, rather then being individual clumps of physically bound galaxies. This conclusion does not change if an empirical fit to the luminosity function of dusty galaxies is used instead of the physical model. The simulations accurately reproduce the statistic of the Planck detections and yield distributions of sizes and ellipticities in qualitative agreement with observations. The redshift distribution of the brightest proto-clusters contributing to the cold peaks has a broad maximum at 1.5 <~ z <~ 3. Therefore follow-up of Planck proto-cluster candidates will provide key information on the high-z evolution of large scale structure.

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin.

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

The assessment of the relationship between radio continuum luminosity and star formation rate (SFR) is of crucial importance to make reliable predictions for the forthcoming ultra-deep radio surveys and to allow a full exploitation of their results to measure the cosmic star formation history. We have addressed this issue by matching recent accurate determinations of the SFR function up to high redshifts with literature estimates of the 1.4 GHz luminosity functions of star forming galaxies (SFGs). This was done considering two options, proposed in the literature, for the relationship between the synchrotron emission ($L_{\rm synch}$), that dominates at 1.4 GHz, and the SFR: a linear relation with a decline of the $L_{\rm synch}$/SFR ratio at low luminosities or a mildly non-linear relation at all luminosities. In both cases we get good agreement with the observed radio luminosity functions but, in the non-linear case, the deviation from linearity must be small. The luminosity function data are consistent with a moderate increase of the $L_{\rm synch}$/SFR ratio with increasing redshift, indicated by other data sets, although a constant ratio cannot be ruled out. A stronger indication of such increase is provided by recent deep 1.4 GHz counts, down to $\mu$Jy levels. This is in contradiction with models predicting a decrease of that ratio due to inverse Compton cooling of relativistic electrons at high redshifts. Synchrotron losses appear to dominate up to $z\simeq 5$. We have also updated the Massardi et al. (2010) evolutionary model for radio loud AGNs.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

The James Webb Space Telescope's Medium Resolution Spectrometer (MRS), will offer nearly 2 orders of magnitude improvement in sensitivity and >3X improvement in spectral resolution over our previous space-based mid-IR spectrometer, the Spitzer IRS. In this paper, we make predictions for spectroscopic pointed observations and serendipitous detections with the MRS. Specifically, pointed observations of Herschel sources require only a few minutes on source integration for detections of several star-forming and active galactic nucleus lines, out to z$=$3 and beyond. But the same data will also include tens of serendipitous 0$\lesssim$z$\lesssim$4 galaxies per field with infrared luminosities ranging $\sim10^6-10^{13}$L$_{\odot}$. In particular, for the first time and for free we will be able to explore the $L_{IR}<10^{9}L_{\odot}$ regime out to $z\sim3$. We estimate that with $\sim$100 such fields, statistics of these detections will be sufficient to constrain the evolution of the low-$L$ end of the infrared luminosity function, and hence the star formation rate function. The above conclusions hold for a wide range in potential low-$L$ end of the IR luminosity function, and accounting for the PAH deficit in low-$L$, low-metallicity galaxies.

We present new estimates of AGN accretion and star-formation luminosity in galaxies obtained for the local 12-$\mu$m sample of Seyfert galaxies (12MGS), by performing a detailed broad-band spectral energy distribution (SED) decomposition including the emission of stars, dust heated by star formation and a possible AGN dusty torus. Thanks to the availability of data from the X-rays to the sub-millimetre, we constrain and test the contribution of the stellar, AGN and star-formation components to the SEDs. The availability of Spitzer-IRS low resolution mid-infrared (mid-IR) spectra is crucial to constrain the dusty torus component at its peak wavelengths. The results of SED-fitting are also tested against the available information in other bands: the reconstructed AGN bolometric luminosity is compared to those derived from X-rays and from the high excitation IR lines tracing AGN activity like [Ne V] and [O IV]. The IR luminosity due to star-formation (SF) and the intrinsic AGN bolometric luminosity are shown to be strongly related to the IR line luminosity. Variations of these relations with different AGN fractions are investigated, showing that the relation dispersions are mainly due to different AGN relative contribution within the galaxy. Extrapolating these local relations between line and SF or AGN luminosities to higher redshifts, by means of recent Herschel galaxy evolution results, we then obtain mid- and far-IR line luminosity functions useful to estimate how many star-forming galaxies and AGN we expect to detect in the different lines at different redshifts and luminosities with future IR facilities (e.g., JWST, SPICA).

We review aspects of Cosmic Microwave Background (CMB) spectral distortions which do not appear to have been fully explored in the literature. In particular, implications of recent evidences of heating of the intergalactic medium (IGM) by feedback from active galactic nuclei are investigated. Taking also into account the IGM heating associated to structure formation, we argue that values of the y parameter of several*10^(-6), i.e. a factor of a few below the COBE/FIRAS upper limit, are to be expected. The Compton scattering by the re-ionized plasma also re-processes primordial distortions, adding a y-type contribution. Hence no pure Bose-Einstein-like distortions are to be expected. An assessment of Galactic and extragalactic foregrounds, taking into account the latest results from the Planck satellite as well as the contributions from the strong CII and CO lines from star-forming galaxies, demonstrates that a foreground subtraction accurate enough to fully exploit the PIXIE sensitivity will be extremely challenging. Motivated by this fact we also discuss methods to detect spectral distortions not requiring absolute measurements and show that accurate determinations of the frequency spectrum of the CMB dipole amplitude may substantially improve over COBE/FIRAS limits on distortion parameters. Such improvements may be at reach of next generation CMB anisotropy experiments. The estimated amplitude of the Cosmic Infrared Background (CIB) dipole might be detectable by careful analyses of Planck maps at the highest frequencies. Thus Planck might provide interesting constraints on the CIB intensity, currently known with a ~30% uncertainty.

We have worked out predictions for the radio counts of star-forming galaxies down to nJy levels, along with redshift distributions down to the detection limits of the phase 1 Square Kilometer Array MID telescope (SKA1-MID) and of its precursors. Such predictions were obtained by coupling epoch dependent star formation rate (SFR) functions with relations between SFR and radio (synchrotron and free-free) emission. The SFR functions were derived taking into account both the dust obscured and the unobscured star-formation, by combining far-infrared (FIR), ultra-violet (UV) and H_alpha luminosity functions up to high redshifts. We have also revisited the South Pole Telescope (SPT) counts of dusty galaxies at 95\,GHz performing a detailed analysis of the Spectral Energy Distributions (SEDs). Our results show that the deepest SKA1-MID surveys will detect high-z galaxies with SFRs two orders of magnitude lower compared to Herschel surveys. The highest redshift tails of the distributions at the detection limits of planned SKA1-MID surveys comprise a substantial fraction of strongly lensed galaxies. We predict that a survey down to 0.25 microJy at 1.4 GHz will detect about 1200 strongly lensed galaxies per square degree, at redshifts of up to 10. For about 30% of them the SKA1-MID will detect at least 2 images. The SKA1-MID will thus provide a comprehensive view of the star formation history throughout the re-ionization epoch, unaffected by dust extinction. We have also provided specific predictions for the EMU/ASKAP and MIGHTEE/MeerKAT surveys.

We present predictions for number counts and redshift distributions of galaxies detectable in continuum and in emission lines with the Mid-infrared (MIR) Instrument (SMI) proposed for the Space Infrared Telescope for Cosmology and Astrophysics (SPICA). We have considered 24 MIR fine-structure lines, four Polycyclic Aromatic Hydrocarbon (PAH) bands (at 6.2, 7.7, 8.6 and 11.3$\mu$m) and two silicate bands (in emission and in absorption) at 9.7$\mu$m and 18.0$\mu$m. Six of these lines are primarily associated with Active Galactic Nuclei (AGNs), the others with star formation. A survey with the SMI spectrometers of 1 hour integration per field-of-view (FoV) over an area of $1\,\hbox{deg}^2$ will yield $5\,\sigma$ detections of $\simeq 140$ AGN lines and of $\simeq 5.2\times10^{4}$ star-forming galaxies, $\simeq 1.6\times10^{4}$ of which will be detected in at least two lines. The combination of a shallow ($20.0\,\hbox{deg}^{2}$, $1.4\times10^{-1}$ h integration per FoV) and a deep survey ($6.9\times10^{-3}\,\hbox{deg}^{2}$, $635$ h integration time), with the SMI camera, for a total of $\sim$1000 h, will accurately determine the MIR number counts of galaxies and of AGNs over five orders of magnitude in flux density, reaching values more than one order of magnitude fainter than the deepest Spitzer $24\,\mu$m surveys. This will allow us to determine the cosmic star formation rate (SFR) function down to SFRs more than 100 times fainter than reached by the Herschel Observatory.

We present new estimates of redshift-dependent luminosity functions of IR lines detectable by SPICA/SAFARI and excited both by star formation and by AGN activity. The new estimates improve over previous work by using updated evolutionary models and dealing in a self consistent way with emission of galaxies as a whole, including both the starburst and the AGN component. New relationships between line and AGN bolometric luminosity have been derived and those between line and IR luminosities of the starburst component have been updated. These ingredients were used to work out predictions for the source counts in 11 mid/far-IR emission lines partially or entirely excited by AGN activity. We find that the statistics of the emission line detection of galaxies as a whole is mainly determined by the star formation rate, because of the rarity of bright AGNs. We also find that the slope of the line integral number counts is flatter than 2 implying that the number of detections at fixed observing time increases more by extending the survey area than by going deeper. We thus propose a wide spectroscopic survey of 1 hour integration per field-of-view over an area of 5 deg$^{2}$ to detect (at 5$\sigma$) $\sim$760 AGNs in [OIV]25.89$\mu$m $-$ the brightest AGN mid-infrared line $-$ out to $z\sim2$. Pointed observations of strongly lensed or hyper-luminous galaxies previously detected by large area surveys such as those by Herschel and by the SPT can provide key information on the galaxy-AGN co-evolution out to higher redshifts.

While continuum imaging data at far-infrared to sub-millimeter wavelengths have provided tight constraints on the population properties of dusty star forming galaxies up to high redshifts, future space missions like the Space Infra-Red Telescope for Cosmology and Astrophysics (SPICA) and ground based facilities like the Cerro Chajnantor Atacama Telescope (CCAT) will allow detailed investigations of their physical properties via their mid-/far-infrared line emission. We present updated predictions for the number counts and the redshift distributions of star forming galaxies spectroscopically detectable by these future missions. These predictions exploit a recent upgrade of evolutionary models, that include the effect of strong gravitational lensing, in the light of the most recent Herschel and South Pole Telescope data. Moreover the relations between line and continuum infrared luminosity are re-assessed, considering also differences among source populations, with the support of extensive simulations that take into account dust obscuration. The derived line luminosity functions are found to be highly sensitive to the spread of the line to continuum luminosity ratios. Estimates of the expected numbers of detections per spectral line by SPICA/SAFARI and by CCAT surveys for different integration times per field of view at fixed total observing time are presented. Comparing with the earlier estimates by Spinoglio et al. (2012) we find, in the case of SPICA/SAFARI, differences within a factor of two in most cases, but occasionally much larger. More substantial differences are found for CCAT.