Modelling the anisotropies in the CIB on all the scales is a challenging task due to the complex nature of the galaxy evolution and thus often requires too many parameters in order to fit the observational data. In this paper, we present a new halo model for the anisotropies of the CIB using only four parameters. Our model connects the mass accretion on the dark matter halos to the star formation rate. Using this model, we find that the halo mass with the maximum efficiency for converting the accreted baryons into stars is $\log_{10}M_\mathrm{max} = {12.94}^{+0.02}_{-0.02} \: M_\odot$, consistent with other studies. Accounting for the mass loss through stellar evolution, we find, for an intermediate age galaxy, that the star formation efficiency defined as $M_\star(z)/M_b(z)$ is equal to 0.19 and 0.21 at redshift 0.1 and 2 respectively, in good agreement with the values obtained by previous studies. It is the first time that a CIB model is used to fit Planck and Herschel CIB power spectra simultaneously. However we find that the large angular scale Planck and Herschel data are not fully compatible with the small scale Herschel data (for $\ell>3000$). The CIB is expected to be correlated with the tSZ signal of the galaxy clusters. Using this halo model for the CIB, and a halo model for the tSZ with a single parameter, we also provide a consistent framework to calculate this cross-correlation, which requires no extra parameter. The CIB$\times$tSZ correlation is found to be higher when inferred with a combination of two widely spaced frequency channels (e.g. 143x857 GHz). The CIB, tSZ, and CIB$\times$tSZ act as foregrounds while measuring the kSZ power spectrum from the CMB power spectrum and need to be removed. Due to its simplistic nature and low number of parameters, the halo model formalism presented here is quite useful for such an analysis to measure the kSZ power spectrum accurately.

We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA Band-7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 16. per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (7 show a counterpart only in the mid-IR and no HST or near-IR identification, while 2 are detected as [CII] emitters at z=5). The 7 HST and near-IR dark galaxies with mid-IR counterpart contribute for about 15 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3.

We present the first [CII] 158 $\mu$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we identified 8 out of 14 detections as [CII] lines at $z\sim 5$, and two as CO transitions at lower redshifts. The remaining 4 lines have an elusive identification in the available catalogs and we considered them as [CII] candidates. We used the 8 confirmed [CII] and the 4 [CII] candidates to build one of the first [CII] LFs at $z\sim 5$. We found that 11 out of these 12 sources have a redshift very similar to that of the ALPINE target in the same pointing, suggesting the presence of overdensities around the targets. Therefore, we split the sample in two (a "clustered" and "field" sub-sample) according to their redshift separation and built two separate LFs. Our estimates suggest that there could be an evolution of the [CII] LF between $z \sim 5$ and $z \sim 0$. By converting the [CII] luminosity to star formation rate we evaluated the cosmic star formation rate density (SFRD) at $z\sim 5$. The clustered sample results in a SFRD $\sim 10$ times higher than previous measurements from UV-selected galaxies. On the other hand, from the field sample (likely representing the average galaxy population) we derived a SFRD $\sim 1.6$ higher compared to current estimates from UV surveys but compatible within the errors. Because of the large uncertainties, observations of larger samples are necessary to better constrain the SFRD at $z\sim 5$. This study represents one of the first efforts aimed at characterizing the demography of [CII] emitters at $z\sim 5$ using a mm-selection of galaxies.

We present the [CII]158$\mu$m line luminosity functions (LFs) at $z\sim 4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A} < -20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [C II] detections and 43 are upper limits. With this UV-selected sample, we set stringent lower bounds to the volume density of [C II] emitters at $z\sim4-6$. Our derived LFs are consistent with the $z\sim0$ [C II] LF at $\leq 10^9L_\odot$, but show a hint of an excess of [C II] emitters at the luminous end ($>10^9L_\odot$) compared to $z\sim0$. This excess is consistent with a potential population of [C II]-bright UV-faint sources that are largely excluded from the ALPINE target selection. To further constrain such sources, in a companion analysis, Loiacono et al. (2020) estimate the [C II] LF based on serendipitous [C II] sources in the ALPINE maps. Their resulting single data point is consistent with the excess we observe at high-L. Combined, these studies suggest possible redshift evolution in the [C II] LFs at high luminosity between $z\sim 5$ and $0$. This conclusion is further supported by the inferred [C II] LFs from the far-IR and CO LFs at $z\sim4-6$. Integrating our estimate on the best-fit [C II] LF, we set a constraint on molecular gas mass density at $z\sim4-6$, with $\rho_{mol} \sim (2-7)\times10^7M_\odot$ Mpc$^{-3}$. This is broadly consistent with previous studies, although the uncertainties span more than an order of magnitude. Finally, we find that available model predictions tend to underestimate the number densities of [C II] emitters at $z\sim5$.

The ALPINE-ALMA large program targets the [CII] 158 $\mu$m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z=4.4 and z=5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details of the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates, we measured the conversion factor from the ALMA 158 $\mu$m rest-frame dust continuum luminosity to the total infrared luminosity (L$_{\rm IR}$) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L$_{\rm IR}$ of 4.4$\times 10^{11}$ L$_\odot$. We also detected 57 additional continuum sources in our ALMA pointings. They are at lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5$\pm$0.2. We measured the 850 $\mu$m number counts between 0.35 and 3.5 mJy, improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40 % of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8$\times$10$^8$ L$_\odot$ and their median full width at half maximum is 252 km/s. After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR-L$_{\rm [CII]}$ relations of De Looze et al. (2014) and Lagache et al. (2018).

The [CII] 158 micron line is one of the strongest IR emission lines, which has been shown to trace the star-formation rate (SFR) of galaxies in the nearby Universe and up to $z \sim 2$. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at $4.4 < z< 5.9$, provides a new opportunity to examine this question with the first statistical dataset. Using the ALPINE data and earlier measurements from the literature we examine the relation between the [CII] luminosity and the SFR over the entire redshift range from $z \sim 4-8$. ALPINE galaxies, which are both detected in [CII] and dust continuum, show a good agreement with the local L([CII])-SFR relation. Galaxies undetected in the continuum with ALMA are found to be over-luminous in [CII], when the UV SFR is used. After accounting for dust-obscured star formation, by an amount SFR(IR)$\approx$SFR(UV) on average, which results from two different stacking methods and SED fitting, the ALPINE galaxies show an L([CII])-SFR relation comparable to the local one. When [CII] non-detections are taken into account, the slope may be marginally steeper at high-z, although this is still somewhat uncertain. When compared in a homogeneous manner, the $z>6 $ [CII] measurements (detections and upper limits) do not behave very differently from the $z \sim 4-6$ data. We find a weak dependence of L([CII])/SFR on the Lyman-alpha equivalent width. Finally, we find that the ratio L([CII])/LIR $\sim (1-3) \times 10^{-3}$ for the ALPINE sources, comparable to that of "normal" galaxies at lower redshift. Our analysis, which includes the largest sample ($\sim 150$ galaxies) of [CII] measurements at $z>4$ available so far, suggests no or little evolution of the L([CII])-SFR relation over the last 13 Gyr of cosmic time.

We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2.

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues.

We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at redshift ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame $158\,\rm{\mu m}$ obtained with the ALMA large program ALPINE. We study the connection between the UV spectral slope ($\beta$), stellar mass ($M_{\star}$), and infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$). Twenty-three galaxies are individually detected in the continuum at $>3.5\,\sigma$ significance. We perform a stacking analysis using both detections and non-detections to study the average dust attenuation properties at z~4.4-5.8. The individual detections and stacks show that the IRX-$\beta$ relation at z~5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts (z<4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud (SMC). This systematic change of the IRX-$\beta$ relation as a function of redshift suggests an evolution of dust attenuation properties at $z>4$. Similarly, we find that our galaxies have lower IRX values up to 1 dex on average at fixed mass compared to previously studied IRX-$M_{\star}$ relations at $z\lesssim4$, albeit with significant scatter. This implies a lower obscured fraction of star-formation than at lower redshifts. Our results suggest that dust properties of UV-selected star forming galaxies at $z\gtrsim4$ are characterised by (i) a steeper attenuation curve than at $z\lesssim4$, and (ii) a rapidly decreasing dust obscured fraction of star formation as a function of redshift. Nevertheless, even among this UV-selected sample, massive galaxies ($\log M_{\star}/M_\odot > 10$) at z~5-6 already exhibit an obscured fraction of star formation of $\sim45\%$, indicating a rapid build-up of dust during the epoch of reionization.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

The main limiting factor of cosmological analyses based on thermal Sunyaev-Zel'dovich (SZ) cluster statistics comes from the bias and systematic uncertainties that affect the estimates of the mass of galaxy clusters. High-angular resolution SZ observations at high redshift are needed to study a potential redshift or morphology dependence of both the mean pressure profile and of the mass-observable scaling relation used in SZ cosmological analyses. The NIKA2 camera is a new generation continuum instrument installed at the IRAM 30-m telescope. With a large field of view, a high angular resolution and a high-sensitivity, the NIKA2 camera has unique SZ mapping capabilities. In this paper, we present the NIKA2 SZ large program, aiming at observing a large sample of clusters at redshifts between 0.5 and 0.9, and the characterization of the first cluster oberved with NIKA2.

Clusters of galaxies are used to map the large-scale structures in the universe and as probe of universe evolution. They can be observed through the Sunyaev-Zel'dovich (SZ) effect. At this respect the spectro-imaging at low resolution frequency is an important tool, today, for the study of cluster of galaxies. We have developed KISS (KIDs-Interferometer-Spectrum-Survey), a spectrometric imager dedicated to the secondary anisotropies of the Cosmic Microwave Background (CMB). The multi-frequency approach permits to improve the component separation with respect to predecessor experiments. In this paper, firstly, we provide a description of the scientific context and the state of the art of SZ observations. Secondly, we describe the KISS instrument. Finally, we show preliminary results of the ongoing commissioning campaign.

NIKA2 is a dual-band millimetric continuum camera of 2900 Kinetic Inductance Detectors (KID), operating at $150$ and $260\,\rm{GHz}$, installed at the IRAM 30-meter telescope. We present the performance assessment of NIKA2 after one year of observation using a dedicated point-source calibration method, referred to as the \emph{baseline} method. Using a large data set acquired between January 2017 and February 2018 that span the whole range of observing elevations and atmospheric conditions encountered at the IRAM 30-m telescope, we test the stability of the performance parameters. We report an instantaneous field of view (FOV) of 6.5' in diameter, filled with an average fraction of $84\%$ and $90\%$ of valid detectors at $150$ and $260\,\rm{GHz}$, respectively. The beam pattern is characterized by a FWHM of $17.6'' \pm 0.1''$ and $11.1''\pm 0.2''$, and a beam efficiency of $77\% \pm 2\%$ and $55\% \pm 3\%$ at $150$ and $260\,\rm{GHz}$, respectively. The rms calibration uncertainties are about $3\%$ at $150\,\rm{GHz}$ and $6\%$ at $260\,\rm{GHz}$. The absolute calibration uncertainties are of $5\%$ and the systematic calibration uncertainties evaluated at the IRAM 30-m reference Winter observing conditions are below $1\%$ in both channels. The noise equivalent flux density (NEFD) at $150$ and $260\,\rm{GHz}$ are of $9 \pm 1\, \rm{mJy}\cdot s^{1/2}$ and $30 \pm 3\, \rm{mJy}\cdot s^{1/2}$. This state-of-the-art performance confers NIKA2 with mapping speeds of $1388 \pm 174$ and $111 \pm 11 \,\rm{arcmin}^2\cdot \rm{mJy}^{-2}\cdot \rm{h}^{-1}$ at $150$ and $260\,\rm{GHz}$. With these unique capabilities of fast dual-band mapping at high (better that 18'') angular resolution, NIKA2 is providing an unprecedented view of the millimetre Universe.

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

The NIKA2 polarization channel at 260 GHz (1.15 mm) has been proposed primarily to observe galactic star-forming regions and probe the critical scales between 0.01-0.05 pc at which magnetic field lines may channel the matter of interstellar filaments into growing dense cores. The NIKA2 polarimeter consists of a room temperature continuously rotating multi-mesh HWP and a cold polarizer that separates the two orthogonal polarizations onto two 260 GHz KIDs arrays. We describe in this paper the preliminary results obtained during the most recent commissioning campaign performed in December 2018. We concentrate here on the analysis of the extended sources, while the observation of compact sources is presented in a companion paper [12]. We present preliminary NIKA2 polarization maps of the Crab nebula. We find that the integrated polarization intensity flux measured by NIKA2 is consistent with expectations.In terms of polarization angle, we are still limited by systematic uncertainties that will be further investigated in the forthcoming commissioning campaigns.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

Mapping millimetre continuum emission has become a key issue in modern multi-wavelength astrophysics. In particular, spectrum-imaging at low frequency resolution is an asset for characterizing the clusters of galaxies via the Sunyaev Zeldovich (SZ) effect. In this context, we have built a ground-based spectrum-imager named KIDs Interferometer Spectrum Survey (KISS). This instrument is based on two 316-pixel arrays of Kinetic Inductance Detectors (KID) cooled to 150 mK by a custom dilution refrigerator-based cryostat. By using Ti-Al and Al absorbers, we can cover a wide frequency range between 80 and 300 GHz. In order to preserve a large instantaneous Field of View (FoV) 1 degree the spectrometer is based on a Fourier Transform interferometer. This represents a technological challenge due to the fast scanning speed that is needed to overcome the effects of background atmospheric fluctuations. KISS is installed at the QUIJOTE 2.25 m telescope in Tenerife since February 2019 and is currently in its commissioning phase. In this proceeding we present an overview of the instrument and the latest results.

One of the main goals of Cosmology is to search for the imprint of primordial gravitational waves in the CMB polarisation field, to probe inflationary theories. One of the obstacles toward the detection of the primordial signal is to extract the B-mode polarisation from astrophysical contaminations. We present a complete analysis of extragalactic foreground contamination due to polarised emission of radio and dusty star-forming galaxies. We update or use up-to-date models that are validated using the most recent measurements. We predict the flux limit (confusion noise) for the future CMB space or balloon experiments (IDS, PIPER, SPIDER, LiteBIRD, PICO), as well as ground-based experiments (C-BASS, NEXT-BASS, QUIJOTE, AdvACTPOL, BICEP3+Keck, BICEPArray, CLASS, SO, SPT3G, S4). Telescope aperture size (and frequency) is the main characteristic impacting the level of confusion noise. Using the flux limits and assuming constant polarisation fractions for radio and dusty galaxies, we compute the B-mode power spectra of the three extragalactic foregrounds (radio source shot noise, dusty galaxy shot noise and clustering), discuss their relative levels and compare their amplitudes to that of the primordial tensor modes parametrized by the tensor-to-scalar ratio r. At the reionization bump (l=5), contamination by extragalactic foregrounds is negligible. At the recombination peak (l=80), while the contamination is much lower than the targeted sensitivity on r for large-aperture telescopes, it is at comparable level for some of the medium- and small-aperture telescope experiments. For example, the contamination is at the level of the 68 per cent confidence level uncertainty on the primordial r for the LiteBIRD and PICO space experiments. Finally we also provide some useful unit conversion factors and give some predictions for the SPICA B-BOP experiment. Abridged

The new NIKA2 camera at the IRAM 30m radiotelescope was used to observe three known debris disks in order to constrain the SED of their dust emission in the millimeter wavelength domain. We have found that the spectral index between the two NIKA2 bands (1mm and 2mm) is consistent with the Rayleigh-Jeans regime (lambda^{-2}), unlike the steeper spectra (lambda^{-3}) measured in the submillimeter-wavelength domain for two of the three disks $-$ around the stars Vega and HD107146. We provide a succesful proof of concept to model this spectral inversion in using two populations of dust grains, those smaller and those larger than a grain radius a0 of 0.5mm. This is obtained in breaking the slope of the size distribution and the functional form of the absorption coefficient of the standard model at a0. The third disk - around the star HR8799 - does not exhibit this spectral inversion but is also the youngest.

We present here the polarization system of the NIKA camera and give a summary of the main results obtained and performed studies on Orion and the Crab nebula. The polarization system was equipped with a room temperature continuously rotating multi-mesh half wave plate and a grid polarizer facing the NIKA cryostat window. NIKA even though less sensitive than NIKA2 had polarization capability in both 1 and 2 millimiter bands. NIKA polarization observations demonstrated the ability of such a technology in detecting the polarization of different targets, compact and extended sources like the Crab nebula and Orion Molecular Cloud region OMC-1. These measurements together with the developed techniques to deal with systematics, opened the way to the current observations of NIKA2 in polarization that will provide important advances in the studies of galactic and extra-galactic emission and magnetic fields.

The steps taken to tailor to NIKA2 observations the Scanamorphos algorithm (initially developed to subtract low-frequency noise from Herschel on-the-fly observations) are described, focussing on the consequences of the different instrument architecture and observation strategy. The method, making the most extensive use of the redundancy built in the multi-scan coverage with large arrays of a given region of the sky, is applicable to extended sources, while the pipeline is so far optimized for compact sources. An example of application is given. A related tool to build consistent broadband SEDs from 60 microns to 2 mm, combining Herschel and NIKA2 data, has also been developed. Its main task is to process the data least affected by low-frequency noise and coverage limitations (i.e. the Herschel data) through the same transfer function as the NIKA2 data, simulating the same scan geometry and applying the same noise and atmospheric signal as extracted from the 1 mm and 2 mm data.

Clusters of galaxies, the largest bound objects in the Universe, constitute a cosmological probe of choice, which is sensitive to both dark matter and dark energy. Within this framework, the Sunyaev-Zel'dovich (SZ) effect has opened a new window for the detection of clusters of galaxies and for the characterization of their physical properties such as mass, pressure and temperature. NIKA, a KID-based dual band camera installed at the IRAM 30-m telescope, was particularly well adapted in terms of frequency, angular resolution, field-of-view and sensitivity, for the mapping of the thermal and kinetic SZ effect in high-redshift clusters. In this paper, we present the NIKA cluster sample and a review of the main results obtained via the measurement of the SZ effect on those clusters: reconstruction of the cluster radial pressure profile, mass, temperature and velocity.

Understanding where and when the mass of stars is determined is one of the fundamental, mostly unsolved, questions in astronomy. Here, we present the first results of GASTON, the Galactic Star Formation with NIKA2 large programme on the IRAM 30m telescope, that aims to identify new populations of low-brightness sources to tackle the question of stellar mass determination across all masses. In this paper, we focus on the high-mass star formation part of the project, for which we map a $\sim2$ deg$^2$ region of the Galactic plane around $l=24^\circ$ in both 1.2 mm and 2.0 mm continuum. Half-way through the project, we reach a sensitivity of 3.7 mJy/beam at 1.2mm. Even though larger than our target sensitivity of 2 mJy, the current sensitivity already allows the identification of a new population of cold, compact sources that remained undetected in any (sub-)mm Galactic plane survey so far. In fact, about 25% of the $\sim 1600$ compact sources identified in the 1.2 mm GASTON image are new detections. We present a preliminary analysis of the physical properties of the GASTON sources as a function of their evolutionary stage, arguing for a potential evolution of the mass distribution of these sources with time.

We present the results of the analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ based on high angular resolution NIKA2 Sunyaev-Zel'dovich (SZ) and $Chandra$ X-ray data. This multi-wavelength analysis enables us to estimate the shape of the temperature profile with unprecedented precision at this redshift and to obtain a map of the gas entropy distribution averaged along the line of sight. The comparison between the cluster morphological properties observed in the NIKA2 and $Chandra$ maps together with the analysis of the entropy map allows us to conclude that MOO J1142$+$1527 is an on-going merger hosting a cool-core at the position of the X-ray peak. This work demonstrates how the addition of spatially-resolved SZ observations to low signal-to-noise X-ray data can bring valuable insights on the intracluster medium thermodynamic properties at $z>1$.

High-resolution mapping of the hot gas in galaxy clusters is a key tool for cluster-based cosmological analyses. Taking advantage of the NIKA2 millimeter camera operated at the IRAM 30-m telescope, the NIKA2 SZ Large Program seeks to get a high-resolution follow-up of 45 galaxy clusters covering a wide mass range at high redshift in order to re-calibrate some of the tools needed for the cosmological exploitation of SZ surveys. We present the second cluster analysis of this program, targeting one of the faintest sources of the sample in order to tackle the difficulties in data reduction for such faint, low-SNR clusters. In this study, the main challenge is the precise estimation of the contamination by sub-millimetric point sources, which greatly affects the tSZ map of the cluster. We account for this contamination by performing a joint fit of the SZ signal and of the flux density of the compact sources. A prior knowledge of these fluxes is given by the adjustment of the SED of each source using data from both NIKA2 and the \textit{Herschel} satellite. The first results are very promising and demonstrate the possibility to estimate thermodynamic properties with NIKA2, even in a compact cluster heavily contaminated by point sources.

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.

The characterization of the Intra-Cluster Medium (ICM) properties of high-redshift galaxy clusters is fundamental to our understanding of large-scale structure formation processes. We present the results of a multi-wavelength analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ discovered as part of the Massive and Distant Clusters of WISE Survey (MaDCoWS). This analysis is based on high angular resolution $Chandra$ X-ray and NIKA2 Sunyaev-Zel'dovich (SZ) data. Although the X-ray data have only about 1700 counts, we are able to determine the ICM thermodynamic radial profiles, namely temperature, entropy, and hydrostatic mass. These have been obtained with unprecedented precision at this redshift and up to $0.7R_{500}$, thanks to the combination of high-resolution X-ray and SZ data. The comparison between the galaxy distribution mapped in infrared by $Spitzer$ and the morphological properties of the ICM derived from the combined analysis of the $Chandra$ and NIKA2 data leads us to the conclusion that the cluster is an on-going merger. We measure the hydrostatic mass profile of the cluster in four angular sectors centered on the large-scale X-ray centroid. This allows us to estimate a systematic uncertainty on the cluster total mass that characterizes both the impact of the observed deviations from spherical symmetry and of the core dynamics on the mass profile. We further combine the X-ray and SZ data at the pixel level to obtain maps of the temperature and entropy distributions averaged along the line of sight. We find a relatively low entropy core at the position of the X-ray peak and high temperature regions located on its south and west sides. The increase in ICM temperature at the location of the SZ peak is expected given the merger dynamics. (abridged)

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

The cosmic infrared background (CIB) is a powerful probe of large-scale structure across a very large redshift range, and consists of unresolved redshifted infrared emission from dusty galaxies. It can be used to study the astrophysics of galaxies, the star formation history of the universe, and the connection between dark and luminous matter. It can furthermore be used as a tracer of the large-scale structure and thus assist in de-lensing of the cosmic microwave background. The major difficulty in its use lies in obtaining accurate and unbiased large-scale CIB images that are cleaned of the contamination by Galactic dust. We used data on neutral atomic hydrogen from the recently-released HI4PI Survey to create template maps of Galactic dust, allowing us to remove this component from the Planck intensity maps from 353 to 857 GHz for approximately $25\%$ of the sky. This allows us to constrain the CIB power spectrum down to $\ell\gtrsim 70$. We present these CIB maps and the various processing and validation steps that we have performed to ensure their quality, as well as a comparison with previous studies. All our data products are made publicly available at https://doi.org/10.7910/DVN/8A1SR3, thereby enabling the community to investigate a wide range of questions related to the universe's large-scale structure.

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.

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $\tau$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $\tau$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $\Lambda$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$\sigma$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged)

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,\sigma$ to $9\,\sigma$. Combined with temperature, lensing is detected at $40\,\sigma$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{\Lambda CDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $\sigma_8 \Omega_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,\sigma$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $\sigma_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $\Omega_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $\sigma_8$ to better than $1\,\%$ precision, finding $\sigma_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

Analysis of the Planck 2018 data set indicates that the statistical properties of the cosmic microwave background (CMB) temperature anisotropies are in excellent agreement with previous studies using the 2013 and 2015 data releases. In particular, they are consistent with the Gaussian predictions of the $\Lambda$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on large angular scales. The novelty of the current study, however, lies in being a first attempt at a comprehensive analysis of the statistics of the polarization signal over all angular scales, using either maps of the Stokes parameters, $Q$ and $U$, or the $E$-mode signal derived from these using a new methodology (which we describe in an appendix). Although remarkable progress has been made in reducing the systematic effects that contaminated the 2015 polarization maps on large angular scales, it is still the case that residual systematics (and our ability to simulate them) can limit some tests of non-Gaussianity and isotropy. However, a detailed set of null tests applied to the maps indicates that these issues do not dominate the analysis on intermediate and large angular scales (i.e., $\ell \lesssim 400$). In this regime, no unambiguous detections of cosmological non-Gaussianity, or of anomalies corresponding to those seen in temperature, are claimed. Notably, the stacking of CMB polarization signals centred on the positions of temperature hot and cold spots exhibits excellent agreement with the $\Lambda$CDM cosmological model, and also gives a clear indication of how Planck provides state-of-the-art measurements of CMB temperature and polarization on degree scales.

The scatter in the relationship between the strength of [CII] 158$\mu$m emission and the star formation rate at high-redshift has been the source of much recent interest. Although the relationship is well-established locally, several intensely star-forming galaxies have been found whose [CII] 158$\mu$m emission is either weak, absent or spatially offset from the young stars. Here we present new ALMA data for the two most distant, gravitationally-lensed and spectroscopically-confirmed galaxies, A2744\_YD4 at $z=$8.38 and MACS1149\_JD1 at $z=$9.11, both of which reveal intense [OIII] 88$\mu$m emission. In both cases we provide stringent upper limits on the presence of [CII] 158$\mu$m with respect to [OIII] 88$\mu$m. We review possible explanations for this apparent redshift-dependent [CII] deficit in the context of our recent hydrodynamical simulations. Our results highlight the importance of using several emission line diagnostics with ALMA to investigate the nature of the interstellar medium in early galaxies.

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.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Galaxy clusters at high redshift are key targets for understanding matter assembly in the early Universe, yet they are challenging to locate. A sample of >2000 high-z candidate structures has been found using Planck's all-sky submm maps, and a sub-set of 234 have been followed up with Herschel-SPIRE, which showed that the emission can be attributed to large overdensities of dusty star-forming galaxies. In order to resolve and characterise the individual galaxies we targeted the eight brightest SPIRE sources in the centre of the Planck peak PLCK G073.4-57.5 using ALMA at 1.3 mm, and complemented these observations with data from IRAC, WIRCam J,K, and SCUBA-2. We detected a total of 18 millimetre galaxies brighter than 0.3 mJy in 2.4 arcmin^2. The ALMA source density is 8-30 times higher than average background estimates and larger than seen in typical 'proto-cluster' fields. We were able to match all but one of the ALMA sources to a NIR counterpart. The most significant (four) SCUBA-2 sources are not included in the ALMA pointings, but we find an 8sigma stacking detection of the ALMA sources in the SCUBA-2 map at 850 um. We derive photo-z, L_IR, SFR, stellar mass, T_dust, M_dust for all of the ALMA galaxies; the photo-zs identify two groups each of five sources, at z~1.5 and 2.4. The two groups show two 'red sequences' (i.e. similar NIR [3.6]-[4.5] colours and different J-K colours). The majority of the ALMA-detected galaxies are on the SFR versus stellar mass main sequence, and half of the sample is more massive than the characteristic stellar mass at the corresponding redshift. Serendipitous CO line detections in two of the galaxies appear to match their photometric redshifts at z~1.54. We performed an analysis of star-formation efficiencies and CO- and mm-continuum-derived gas fractions of our ALMA sources, combined with a sample of 1<z<3 cluster and proto-cluster members.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $\psi$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $\psi$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas.

Observations of the high-redshift Universe using the 21 cm line of neutral hydrogen and complimentary emission lines from the first galaxies promise to open a new door for our understanding of the epoch of reionization. We present predictions for the [C II] 158-micron line and H I 21 cm emission from redshifts z=6--9 using high-dynamic-range cosmological simulations combined with semi-analytical models. We find that the CONCERTO experiment should be able to detect the large scale power spectrum of [C II] emission to redshifts of up to z=8 (signal-to-noise ratio ~ 1 at k = 0.1 h/cMpc with 1500 hr of integration). A Stage II experiment similar to CCAT-p should be able to detect [C II] from even higher redshifts to high significance for similar integration times (signal-to-noise ratio of ~50 at k = 0.2 h/cMpc at z=6--9). We study the possibility of combining such future [C II] measurements with 21 cm measurements using LOFAR and SKA to measure the [C II]-21cm cross power spectra, and find that a Stage II experiment should be able to measure the cross-power spectrum for k < 1 h/cMpc to signal-to-noise ratio of better than 10. We discuss the capability of such measurements to constrain astrophysical parameters relevant to reionization and show that a measurement of the [C II]-21cm cross power spectrum helps break the degeneracy between the mass and brightness of ionizing sources.

Dust is a crucial component of the interstellar medium of galaxies. The presence of dust strongly affects the light produced by stars within a galaxy. As these photons are our main information vector to explore the stellar mass assembly and therefore understand a galaxy's evolution, modeling the luminous properties of galaxies and taking into account the impact of the dust is a fundamental challenge for semi-analytical models.We present the complete prescription of dust attenuation implemented in the new semi-analytical model: G.A.S. This model is based on a two-phase medium originating from a physically motivated turbulent model of gas structuring (G.A.S. I paper). Dust impact is treated by taking into account three dust components: Polycyclic Aromatic Hydrocarbons, Very Small Grains, and Big Grains. All three components evolve in both a diffuse and a fragmented/dense gas phase. Each phase has its own stars, dust content and geometry. Dust content evolves according to the metallicity of it associated phase.The G.A.S. model is used to predict both the UV and the IR luminosity functions from $z=9.0$ to $z=0.1$. Our two-phase ISM prescription catches very well the evolution of UV and IR luminosity functions. We note a small overproduction of the IR luminosity at low redshift ($z<0.5$). We also focus on the Infrared-Excess (IRX) and explore its dependency with the stellar mass, UV slope, stellar age, metallicity and slope of the attenuation curves. Our model predicts large scatters for relations based on IRX, especially for the IRX-$\beta$ relation. Our analysis reveals that the slope of the attenuation curve is more driven by absolute attenuation in the FUV band than by disk inclination. We confirm that the age of the stellar population and the slope of the attenuation curve can both shift galaxies below the fiducial star-birth relation in the IRX-$\beta$ diagram.

We present a linear clustering model of cosmic infrared background (CIB) anisotropies at large scales that is used to measure the cosmic star formation rate density up to redshift 6, the effective bias of the CIB and the mass of dark-matter halos hosting dusty star-forming galaxies. This is achieved using the Planck CIB auto- and cross-power spectra (between different frequencies) and CIBxCMB lensing cross-spectra measurements, as well as external constraints (e.g. on the CIB mean brightness). We recovered an obscured star formation history which agrees well with the values derived from infrared deep surveys and we confirm that the obscured star formation dominates the unobscured one up to at least z=4. The obscured and unobscured star formation rate densities are compatible at $1\sigma$ at z=5. We also determined the evolution of the effective bias of the galaxies emitting the CIB and found a rapid increase from $\sim$0.8 at z$=$0 to $\sim$8 at z$=$4. At 2$<$z$<$4, this effective bias is similar to that of galaxies at the knee of the mass functions and submillimeter galaxies. This effective bias is the weighted average of the true bias with the corresponding emissivity of the galaxies. The halo mass corresponding to this bias is thus not exactly the mass contributing the most to the star formation density. Correcting for this, we obtained a value of log(M$_h$/M$_{\odot}$)=12.77$_{-0.125}^{+0.128}$ for the mass of the typical dark matter halo contributing to the CIB at z=2. Finally, we also computed using a Fisher matrix analysis how the uncertainties on the cosmological parameters affect the recovered CIB model parameters and find that the effect is negligible.

There is a lack of large samples of spectroscopically confirmed clusters and protoclusters at high redshifts, $z>$1.5. Discovering and characterizing distant (proto-)clusters is important for yielding insights into the formation of large-scale structure and on the physical processes responsible for regulating star-formation in galaxies in dense environments. The Spitzer Planck Herschel Infrared Cluster (SPHerIC) survey was initiated to identify these characteristically faint and dust-reddened sources during the epoch of their early assembly. We present Spitzer IRAC observations of 82 galaxy (proto-)cluster candidates at 1.3<$z_p$<3.0 that were vetted in a two step process: (1) using Planck to select by color those sources with the highest star-formation rates, and (2) using Herschel at higher resolution to separate out the individual red sources. The addition of the Spitzer data enables efficient detection of the central and massive brightest red cluster galaxies (BRCGs). We find that BRCGs are associated with highly significant, extended and crowded regions of IRAC sources which are more overdense than the field. This result corroborates our hypothesis that BRCGs within the Planck - Herschel sources trace some of the densest and actively star-forming proto-clusters in the early Universe. On the basis of a richness-mass proxy relation, we obtain an estimate of their mean masses which suggests our sample consists of some of the most massive clusters at z$\approx$2 and are the likely progenitors of the most massive clusters observed today.

Galaxies are often used as tracers of the large scale structure (LSS) to measure the Integrated Sachs-Wolfe effect (ISW) by cross-correlating the galaxy survey maps with the Cosmic Microwave Background (CMB) map. We use the Cosmic Infrared Background (CIB) as a tracer of the LSS to perform a theoretical CIB-CMB cross-correlation to measure the ISW for different Planck HFI frequencies. We discuss the detectability of this ISW signal using a Signal-to-noise ratio analysis and find that the ISW detected this way can provide us with the highest SNR for a single tracer ranging from 5 to 6.7 (maximum being for 857 GHz) with the CIB and CMB maps extracted over the whole sky. A Fisher matrix analysis showed that this measurement of the ISW can improve the constraints on the cosmological parameters; especially the equation of state of the dark energy $w$ by $\sim 47\%$. Performing a more realistic analysis including the galactic dust residuals in the CIB maps over realistic sky fractions shows that the dust power spectra dominate over the CIB power spectra at $\ell < 100$ and ISW can't be detected with high SNR. We perform the cross-correlation on the existing CIB-CMB maps over $\sim 11\%$ of the sky in the southern hemisphere and find that the ISW is not detected with the existing CIB maps over such small sky fractions.

We explore the use of the cosmic infrared background as a tracer of the LSS for cross-correlating with the CMB and exploit the ISW. We use the improved linear CIB model of Maniyar et al (2018) and derive the theoretical CIBxISW cross-correlation for Planck HFI frequencies 217, 353, 545 and 857 GHz and IRAS 3000 GHz. We predict a positive cross-correlation between the CIB and CMB whose amplitude decreases rapidly at small scales. We perform a signal-to-noise ratio (SNR) analysis on this cross-correlation. In the ideal case of the cross-correlation obtained over 70% (40%) of the sky with no residual contaminants (e.g. galactic dust) in maps, the SNR ranges from 4.2-5.6 (3.2-4.3) with the highest for 857 GHz. A Fisher matrix analysis shows that an ISW signal detected with such high SNR on the 40% sky can improve the constraints on the cosmological parameters considerably; constraints on the equation of state of the dark energy are improved by 80%. We then perform a more realistic analysis with the effect of residual galactic dust contamination in CIB maps. We calculate the dust power spectra for several frequencies and sky fractions which dominate over CIB at lower multipoles we are interested in. Considering conservative 10% residual level of galactic dust in the CIB power spectra, we find that the SNR drops drastically making ISW detection difficult. To check the capability of current maps to detect ISW via this method, we measure the cross-correlation of the CIB and the CMB Planck maps on so called GASS field covering an area of 11% in the southern hemisphere. We find that with such a small sky fraction and dust residuals present in CIB maps, we do not detect ISW signal and the measured signal is consistent with zero. In order not to degrade the SNR for the ISW measurement by more than 10% on the 40% sky, we find that the dust needs to be cleaned up to 0.01% level on the power spectrum.

NIKA2 is a dual-band millimetric camera of thousands of Kinetic Inductance Detectors (KID) installed at the IRAM 30-meter telescope in the Spanish Sierra Nevada. The instrument commissioning was completed in September 2017, and NIKA2 is now open to the scientific community and will operate for the next decade. NIKA2 has well-adapted instrumental design and performance to produce high-resolution maps of the thermal Sunyaev-Zel'dovich (SZ) effect toward intermediate and high redshift galaxy clusters. Moreover, it benefits from a guaranteed time large program dedicated to mapping a representative sample of galaxy clusters via SZ and that includes X-ray follow-ups. The main expected outputs of the SZ large program are the constraints on the redshift evolution of the pressure profile and the mass-observable relation. The first SZ mapping of a galaxy cluster with NIKA2 was produced, as part of the SZ large program. We found a sizable impact of the intracluster medium dynamics on the integrated SZ observables. This shows NIKA2 capabilities for the precise characterisation of the mass-observable relation that is required for accurate cosmology with galaxy clusters.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes

We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow closely those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3 degree regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of beta_d = 1.55 +/- 0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of beta_s = -3.1 +/- 0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects.

This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous 2015 release. They enabled the first significant measurement of the reionization optical depth parameter using HFI data. This paper presents an extensive analysis of systematic effects, including the use of simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved. Calibration, based on the CMB dipole, is now extremely accurate and in the frequency range 100 to 353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than $0.35\mu$K, an accuracy of order $10^{-4}$. This is a major legacy from the HFI for future CMB experiments. The removal of bandpass leakage has been improved by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. Simulations reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect. Using these simulations, we measure and correct the small frequency calibration bias induced by this systematic effect at the $10^{-4}$ level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the $10^{-3}$ level.

The Crab nebula is a supernova remnant exhibiting a highly polarized synchrotron radiation at radio and millimeter wavelengths. It is the brightest source in the microwave sky with an extension of 7 by 5 arcminutes and commonly used as a standard candle for any experiment which aims at measuring the polarization of the sky. Though its spectral energy distribution has been well characterized in total intensity, polarization data are still lacking at millimetre wavelengths. We report in this paper high resolution (18 arcsec FWHM) observations of the Crab nebula in total intensity and linear polarization at 150 GHz with the NIKA camera. NIKA, operated at the IRAM 30 m telescope from 2012 to 2015, is a camera made of Lumped Element Kinetic Inductance Detectors (LEKIDs) observing the sky at 150 and 260 GHz. From these observations we are able to reconstruct the spatial distribution of the polarization degree and angle of the Crab nebula, which is found to be compatible with previous observations at lower and higher frequencies. Averaging across the source and using other existing data sets we find that the Crab nebula polarization angle is consistent with being constant over a wide range of frequencies with a value of -87.7$^\circ$ +- 0.3 in Galactic coordinates. We also present the first estimation of the Crab nebula spectral energy distribution polarized flux in a wide frequency range: 30-353 GHz. Assuming a single power law emission model we find that the polarization spectral index $\beta_{pol}$ = - 0.347 +- 0.026 is compatible with the intensity spectral index $\beta$ = - 0.323 +- 0.001.

The complete characterization of the pressure profile of high-redshift galaxy clusters, from their core to their outskirts, is a major issue for the study of the formation of large-scale structures. It is essential to constrain a potential redshift evolution of both the slope and scatter of the mass-observable scaling relations used in cosmology studies based on cluster statistics. In this paper, we present the first thermal Sunyaev-Zel'dovich (tSZ) mapping of a cluster from the sample of the NIKA2 SZ large program that aims at constraining the redshift evolution of cluster pressure profiles and the tSZ-mass scaling relation. We have observed the galaxy cluster PSZ2 G144.83+25.11 at redshift $z=0.58$ with the NIKA2 camera, a dual-band (150 and 260 GHz) instrument operated at the IRAM 30-meter telescope. We identify a thermal pressure excess in the south-west region of PSZ2 G144.83+25.11 and a high redshift sub-millimeter point source that affect the intracluster medium (ICM) morphology of the cluster. The NIKA2 data are used jointly with tSZ data acquired by the MUSTANG, Bolocam and $Planck$ experiments in order to non-parametrically set the best constraints on the electronic pressure distribution from the cluster core ($\rm{R} \sim 0.02 \rm{R_{500}}$) to its outskirts ($\rm{R} \sim 3 \rm{R_{500}} $). We investigate the impact of the over-pressure region on the shape of the pressure profile and on the constraints on the integrated Compton parameter $\rm{Y_{500}}$. A hydrostatic mass analysis is also performed by combining the tSZ-constrained pressure profile with the deprojected electronic density profile from XMM-$Newton$. This allows us to conclude that the estimates of $\rm{Y_{500}}$ and $\rm{M_{500}}$ obtained from the analysis with and without masking the disturbed ICM region differ by 65 and 79% respectively. (abridged)

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise ratio of 6, and is more than a factor of 2 larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, $(1-b)$. In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as an independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those from the Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of $(1-b)$ the tension is mild, only a little over one standard deviation, while it remains substantial ($3.7\,\sigma$) for the largest estimated value. We also examine constraints on extensions to the base flat $\Lambda$CDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base $\Lambda$CDM model.

The thermal emission of dust grains is a powerful tool for probing cold, dense regions of molecular gas in the ISM, and so constraining dust properties is key to obtaining accurate measurements of dust mass and temperature. By placing constraints on the dust emissivity spectral index, beta, towards two star-forming infrared dark clouds, SDC18.888 and SDC24.489, we evaluate the role of mass concentration in the associated star-formation activity. We exploit the simultaneous 1.2mm and 2.0mm imaging capability of NIKA on the IRAM 30m telescope to construct maps of beta for both clouds, and by incorporating Herschel observations, we create H2 column density maps with 13" resolution. While we find no significant systematic radial variations around the most massive clumps in either cloud on >0.1 pc scales, their mean beta values are significantly different, with beta = 2.07 +/- 0.09 (rand) +/- 0.25 (syst) for SDC18.888 and beta = 1.71 +/- 0.09 (rand) +/- 0.25 (syst) for SDC24.489. These differences could be a consequence of the very different environments in which both clouds lie, and we suggest that the proximity of SDC18.888 to the W39 HII region may raise beta on scales of 1 pc. We also find that the mass in SDC24.489 is more centrally concentrated and circularly symmetric than in SDC18.888, and is consistent with a scenario in which spherical globally-collapsing clouds concentrate a higher fraction of their mass into a single core than elongated clouds that will more easily fragment, distributing their mass into many cores. We demonstrate that beta variations towards interstellar clouds can be robustly constrained with high-SNR NIKA observations, providing more accurate estimates of their masses. The methods presented here will be applied to the Galactic Star Formation with NIKA2 (GASTON) large programme, extending our analysis to a statistically significant sample of star-forming clouds.

The New IRAM KID Arrays 2 (NIKA2) consortium has just finished installing and commissioning a millimetre camera on the IRAM 30 m telescope. It is a dual-band camera operating with three frequency multiplexed kilo-pixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK, designed to observe the intensity and polarisation of the sky at 260 and 150 GHz (1.15 and 2 mm). NIKA2 is today an IRAM resident instrument for millimetre astronomy, such as Intra Cluster Medium from intermediate to distant clusters and so for the follow-up of Planck satellite detected clusters, high redshift sources and quasars, early stages of star formation and nearby galaxies emission. We present an overview of the instrument performance as it has been evaluated at the end of the commissioning phase.

In the last decade, it has become clear that the dust-enshrouded star formation contributes significantly to early galaxy evolution. Detection of dust is therefore essential in determining the properties of galaxies in the high-redshift universe. This requires observations at the (sub-)millimeter wavelengths. Unfortunately, sensitivity and background confusion of single dish observations on the one hand, and mapping efficiency of interferometers on the other hand, pose unique challenges to observers. One promising route to overcome these difficulties is intensity mapping of fluctuations which exploits the confusion-limited regime and measures the collective light emission from all sources, including unresolved faint galaxies. We discuss in this contribution how 2D and 3D intensity mapping can measure the dusty star formation at high redshift, through the Cosmic Infrared Background (2D) and [CII] fine structure transition (3D) anisotropies.

Sub-structures in the hot gas of galaxy clusters are related to their formation history and to the astrophysical processes at play in the intracluster medium (ICM). The thermal Sunyaev-Zel'dovich (tSZ) effect is directly sensitive to the line-of-sight integrated ICM pressure, and is thus particularly adapted to study ICM sub-structures. We apply structure-enhancement filtering algorithms to high resolution tSZ observations of distant clusters, in order to search for pressure discontinuities, compressions, as well as secondary peaks in the ICM. The same filters are applied to synthetic tSZ images extracted from RHAPSODY-G hydrodynamic simulations, in order to better interpret the extracted features. We also study the noise propagation trough the filters and quantify the impact of systematic effects, point source residuals being identified as the dominant potential contaminant. In 3 of our 6 NIKA clusters we identify features at high S/N that show clear evidence for merger events. In MACSJ0717 (z=0.55), three strong pressure gradients are observed on the E, SE and W sectors, and two main peaks in the pressure distribution are identified. We observe a lack of tSZ compact structure in the cool-core cluster PSZ1G045.85 (z=0.61), and a tSZ gradient ridge dominates in the SE. In the highest z cluster, CLJ1227 (z=0.89), we detect a ~45" (360 kpc) long ridge pressure gradient associated with a secondary pressure peak in the W region. Our results show that current tSZ facilities have now reached the angular resolution and sensitivity to allow an exploration of the details of pressure sub-structures in clusters, even at high z. This opens the possibility to quantify the impact of the dynamical state on the relation between the tSZ signal and the mass of clusters, which is important when using tSZ clusters to test cosmological models. This work also marks the first NIKA cluster sample data release.

Millimeter-wave continuum astronomy is today an indispensable tool for both general Astrophysics studies and Cosmology. General purpose, large field-of-view instruments are needed to map the sky at intermediate angular scales not accessible by the high-resolution interferometers and by the coarse angular resolution space-borne or ground-based surveys. These instruments have to be installed at the focal plane of the largest single-dish telescopes. In this context, we have constructed and deployed a multi-thousands pixels dual-band (150 and 260 GHz, respectively 2mm and 1.15mm wavelengths) camera to image an instantaneous field-of-view of 6.5arc-min and configurable to map the linear polarization at 260GHz. We are providing a detailed description of this instrument, named NIKA2 (New IRAM KID Arrays 2), in particular focusing on the cryogenics, the optics, the focal plane arrays based on Kinetic Inductance Detectors (KID) and the readout electronics. We are presenting the performance measured on the sky during the commissioning runs that took place between October 2015 and April 2017 at the 30-meter IRAM (Institut of Millimetric Radio Astronomy) telescope at Pico Veleta. NIKA2 has been successfully deployed and commissioned, performing in-line with the ambitious expectations. In particular, NIKA2 exhibits FWHM angular resolutions of around 11 and 17.5 arc-seconds at respectively 260 and 150GHz. The NEFD (Noise Equivalent Flux Densities) demonstrated on the maps are, at these two respective frequencies, 33 and 8 mJy*sqrt(s). A first successful science verification run has been achieved in April 2017. The instrument is currently offered to the astronomical community during the coming winter and will remain available for at least the next ten years.

Galaxy clusters constitute a major cosmological probe. However, Planck 2015 results have shown a weak tension between CMB-derived and cluster-derived cosmological parameters. This tension might be due to poor knowledge of the cluster mass and observable relationship. As for now, arcmin resolution Sunyaev-Zeldovich (SZ) observations ({\it e.g.} SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters ($z<0.2$). For high redshift clusters ($z>0.5$) high resolution and high sensitivity SZ observations are needed. With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain) is particularly well adapted for these observations. The NIKA2 SZ observation program will map a large sample of clusters (50) at redshifts between 0.5 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.}

The [CII] fine structure transition at 158 microns is the dominant cooling line of cool interstellar gas, and is the brightest of emission lines from star forming galaxies from FIR through meter wavelengths. With the advent of ALMA and NOEMA, capable of detecting [CII]-line emission in high-redshift galaxies, there has been a growing interest in using the [CII] line as a probe of the physical conditions of the gas in galaxies, and as a SFR indicator at z>4. In this paper, we use a semi-analytical model of galaxy evolution (G.A.S.) combined with the code CLOUDY to predict the [CII] luminosity of a large number of galaxies at 4< z<8. At such high redshift, the CMB represents a strong background and we discuss its effects on the luminosity of the [CII] line. We study the LCII-SFR and LCII-Zg relations and show that they do not strongly evolve with redshift from z=4 and to z=8. Galaxies with higher [CII] luminosities tend to have higher metallicities and higher star formation rates but the correlations are very broad, with a scatter of about 0.5 dex for LCII-SFR. Our model reproduces the LCII-SFR relations observed in high-redshift star-forming galaxies, with [CII] luminosities lower than expected from local LCII-SFR relations. Accordingly, the local observed LCII-SFR relation does not apply at high-z. Our model naturally produces the [CII] deficit, which appears to be strongly correlated with the intensity of the radiation field in our simulated galaxies. We then predict the [CII] luminosity function, and show that it has a power law form in the range of LCII probed by the model with a slope alpha=1. The slope is not evolving from z=4 to z=8 but the number density of [CII]-emitters decreases by a factor of 20x. We discuss our predictions in the context of current observational estimates on both the differential and cumulative luminosity functions.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be $n_\mathrm{s} = 0.968 \pm 0.006$ and tightly constrain its scale dependence to $d n_s/d \ln k =-0.003 \pm 0.007$ when combined with the Planck lensing likelihood. When the high-$\ell$ polarization data is included, the results are consistent and uncertainties are reduced. The upper bound on the tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode polarization constraint $r< 0.12$ (95% CL) obtained from a joint BICEP2/Keck Array and Planck analysis. These results imply that $V(\phi) \propto \phi^2$ and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as $R^2$ inflation. Three independent methods reconstructing the primordial power spectrum are investigated. The Planck data are consistent with adiabatic primordial perturbations. We investigate inflationary models producing an anisotropic modulation of the primordial curvature power spectrum as well as generalized models of inflation not governed by a scalar field with a canonical kinetic term. The 2015 results are consistent with the 2013 analysis based on the nominal mission data.

The development of precision cosmology with clusters of galaxies requires high-angular resolution Sunyaev-Zel'dovich (SZ) observations. As for now, arcmin resolution SZ observations (e.g. SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters (z<0.2). With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain), will bring valuable information in the field of SZ imaging of clusters of galaxies. The NIKA2 SZ observation program will allow us to observe a large sample of clusters (50) at redshifts between 0.4 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.

Follow-up observations at high-angular resolution of submillimeter galaxies showed that the single-dish sources are comprised of a blend of several galaxies. Consequently, number counts derived from low and high angular resolution observations are in disagreement. This demonstrates the importance of resolution effects and the need to have realistic simulations to explore them. We built a new 2deg^2 simulation of the extragalactic sky from the far-infrared to the submillimeter. It is based on an updated version of the two star-formation mode galaxy evolution model. Using global galaxy properties, we use the abundance matching technique to populate a dark-matter lightcone and thus simulate the clustering. We produce maps and extract the sources, and show that the limited angular resolution of single-dish instruments have a strong impact on (sub)millimeter continuum observations. Taking into account these resolution effects, we are reproducing a large set of observables, including number counts, redshift distributions, and cosmic infrared background power spectra. Our simulation describes consistently the number counts from single-dish telescopes and interferometers. In particular, at 350 and 500 um, we find that number counts measured by Herschel between 5 and 50 mJy are biased towards high values by a factor 2, and that redshift distributions are biased towards low z. We also show that the clustering has an important impact on the Herschel pixel histogram used to derive number counts from P(D) analysis. Finally, we demonstrate that the large number density of red Herschel sources found in observations but not in models could be an observational artifact caused by the combination of noise, resolution effects, and the steepness of color and flux density distributions. Our simulation (SIDES) is available at http://cesam.lam.fr/sides

We propose a method to map the temperature distribution of the hot gas in galaxy clusters that uses resolved images of the thermal Sunyaev-Zel'dovich (tSZ) effect in combination with X-ray data. Application to images from the New IRAM KIDs Array (NIKA) and XMM-Newton allows us to measure and determine the spatial distribution of the gas temperature in the merging cluster MACS J0717.5+3745, at $z=0.55$. Despite the complexity of the target object, we find a good morphological agreement between the temperature maps derived from X-ray spectroscopy only -- using XMM-Newton ($T_{\rm XMM}$) and Chandra ($T_{\rm CXO}$) -- and the new gas-mass-weighted tSZ+X-ray imaging method ($T_{\rm SZ}$). We correlate the temperatures from tSZ+X-ray imaging and those from X-ray spectroscopy alone and find that $T_{\rm SZ}$ is higher than $T_{\rm XMM}$ and lower than $T_{\rm CXO}$ by $\sim 10\%$ in both cases. Our results are limited by uncertainties in the geometry of the cluster gas, contamination from kinetic SZ ($\sim 10\%$), and the absolute calibration of the tSZ map ($7\%$). Investigation using a larger sample of clusters would help minimise these effects.

Context: In the past decade, sensitive, resolved Sunyaev-Zel'dovich (SZ) studies of galaxy clusters have become common. Whereas many previous SZ studies have parameterized the pressure profiles of galaxy clusters, non-parametric reconstructions will provide insights into the thermodynamic state of the intracluster medium (ICM). Aims: We seek to recover the non-parametric pressure profiles of the high redshift ($z=0.89$) galaxy cluster CLJ 1226.9+3332 as inferred from SZ data from the MUSTANG, NIKA, Bolocam, and Planck instruments, which all probe different angular scales. Methods: Our non-parametric algorithm makes use of logarithmic interpolation, which under the assumption of ellipsoidal symmetry is analytically integrable. For MUSTANG, NIKA, and Bolocam we derive a non-parametric pressure profile independently and find good agreement among the instruments. In particular, we find that the non-parametric profiles are consistent with a fitted gNFW profile. Given the ability of Planck to constrain the total signal, we include a prior on the integrated Compton Y parameter as determined by Planck. Results: For a given instrument, constraints on the pressure profile diminish rapidly beyond the field of view. The overlap in spatial scales probed by these four datasets is therefore critical in checking for consistency between instruments. By using multiple instruments, our analysis of CLJ 1226.9+3332 covers a large radial range, from the central regions to the cluster outskirts: $0.05 R_{500} < r < 1.1 R_{500}$. This is a wider range of spatial scales than is typical recovered by SZ instruments. Similar analyses will be possible with the new generation of SZ instruments such as NIKA2 and MUSTANG2.

The characterization of dust properties in the interstellar medium (ISM) is key for star formation. Mass estimates are crucial to determine gravitational collapse conditions for the birth of new stellar objects in molecular clouds. However, most of these estimates rely on dust models that need further observational constraints from clouds to prestellar and protostellar cores. We present results of a study of dust emissivity changes based on mm-continuum data obtained with the NIKA camera at the IRAM-30m telescope. Observing dust emission at 1.15 mm and 2 mm allows us to constrain the dust emissivity index ($\beta$) in the Rayleigh-Jeans tail of the dust spectral energy distribution (SED) far from its peak emission, where the contribution of other parameters (i.e. dust temperature) is important. Focusing on the Taurus molecular cloud, a low-mass star-forming regions in the Gould Belt, we analyze the emission properties of several distinct objects in the B213 filament: three prestellar cores, two Class-0/I protostellar cores and one Class-II object. By means of the ratio of the two NIKA channel-maps, we show that in the Rayleigh-Jeans approximation the dust emissivity index varies among the objects. For one prestellar and two protostellar cores, we produce a robust study using Herschel data to constrain the dust temperature of the sources. By using the Abel transform inversion technique we get accurate radial $\beta$ profiles. We find systematic spatial variations of $\beta$ in the protostellar cores that is not observed in the prestellar core. While in the former case $\beta$ decreases toward the center, in the latter it remains constant. Moreover, $\beta$ appears anticorrelated with the dust temperature. We discuss the implication of these results in terms of dust grain evolution between pre- and protostellar cores.

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.

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 six parameters of the standard $\Lambda$CDM model have best-fit values derived from the Planck temperature power spectrum that are shifted somewhat from the best-fit values derived from WMAP data. These shifts are driven by features in the Planck temperature power spectrum at angular scales that had never before been measured to cosmic-variance level precision. We investigate these shifts to determine whether they are within the range of expectation and to understand their origin in the data. Taking our parameter set to be the optical depth of the reionized intergalactic medium $\tau$, the baryon density $\omega_{\rm b}$, the matter density $\omega_{\rm m}$, the angular size of the sound horizon $\theta_*$, the spectral index of the primordial power spectrum, $n_{\rm s}$, and $A_{\rm s}e^{-2\tau}$ (where $A_{\rm s}$ is the amplitude of the primordial power spectrum), we examine the change in best-fit values between a WMAP-like large angular-scale data set (with multipole moment $\ell<800$ in the Planck temperature power spectrum) and an all angular-scale data set ($\ell<2500$ Planck temperature power spectrum), each with a prior on $\tau$ of $0.07\pm0.02$. We find that the shifts, in units of the 1$\sigma$ expected dispersion for each parameter, are $\{\Delta \tau, \Delta A_{\rm s} e^{-2\tau}, \Delta n_{\rm s}, \Delta \omega_{\rm m}, \Delta \omega_{\rm b}, \Delta \theta_*\} = \{-1.7, -2.2, 1.2, -2.0, 1.1, 0.9\}$, with a $\chi^2$ value of 8.0. We find that this $\chi^2$ value is exceeded in 15% of our simulated data sets, and that a parameter deviates by more than 2.2$\sigma$ in 9% of simulated data sets, meaning that the shifts are not unusually large. Comparing $\ell<800$ instead to $\ell>800$, or splitting at a different multipole, yields similar results. We examine the $\ell<800$ model residuals in the $\ell>800$ power spectrum data and find that the features there... [abridged]

We explore the detection, with upcoming spectroscopic surveys, of three-dimensional power spectra of emission line fluctuations produced in different phases of the Interstellar Medium (ISM) by ionized carbon, ionized nitrogen and neutral oxygen at redshift z>4. The emission line [CII] from ionized carbon at 157.7 micron, and multiple emission lines from carbon monoxide, are the main targets of planned ground-based surveys, and an important foreground for future space-based surveys like the Primordial Inflation Explorer (PIXIE). However, the oxygen [OI] (145.5 micron) line, and the nitrogen [NII] (121.9 micron and 205.2 micron) lines, might be detected in correlation with [CII] with reasonable signal-to-noise ratio (SNR). These lines are important coolants of both the neutral and the ionized medium, and probe multiple phases of the ISM. We compute predictions of the three-dimensional power spectra for two surveys designed to target the [CII] line, showing that they have the required sensitivity to detect cross-power spectra with the [OI] line, and the [NII] lines with sufficient SNR. The importance of cross-correlating multiple lines is twofold. On the one hand, we will have multiple probes of the different phases of the ISM, which is key to understand the interplay between energetic sources, the gas and dust at high redshift. This kind of studies will be useful for a next-generation space observatory such as the NASA Far-IR Surveyor. On the other end, emission lines from external galaxies are an important foreground when measuring spectral distortions of the Cosmic Microwave Background spectrum with future space-based experiments like PIXIE; measuring fluctuations in the intensity mapping regime will help constraining the mean amplitude of these lines, and will allow us to better handle this important foreground.

The Draco nebula is a high Galactic latitude interstellar cloud likely to have been formed by the collision of a Galactic halo cloud entering the disk of the Milky Way. Such conditions are ideal to study the formation of cold and dense gas in colliding flows of warm gas. We present Herschel-SPIRE observations that reveal the fragmented structure of the interface between the infalling cloud and the Galactic layer. This front is characterized by a Rayleigh-Taylor instability structure. From the determination of the typical length of the periodic structure (2.2 pc) we estimated the gas kinematic viscosity and the turbulence dissipation scale (0.1 pc) that is compatible with that expected if ambipolar diffusion is the main mechanism of energy dissipation in the WNM. The small-scale structures of the nebula are typical of that seen in some molecular clouds. The gas density has a log-normal distribution with an average value of $10^3$ cm$^{-3}$. The size of the structures is 0.1-0.2 pc but this estimate is limited by the resolution of the observations. The mass ranges from 0.2 to 20 M$_{\odot}$ and the distribution of the more massive clumps follows a power law $dN/d\log(M) \sim M^{-1.4}$. We identify a mass-size relation with the same exponent as that found in GMCs ($M\sim L^{2.3}$) but only 15% of the mass of the cloud is in gravitationally bound structures. We conclude that the increase of pressure in the collision is strong enough to trigger the WNM-CNM transition caused by the interplay between turbulence and thermal instability as self-gravity is not dominating the dynamics.

We present SCUBA-2 follow-up of 61 candidate high-redshift Planck sources. Of these, 10 are confirmed strong gravitational lenses and comprise some of the brightest such submm sources on the observed sky, while 51 are candidate proto-cluster fields undergoing massive starburst events. With the accompanying Herschel-SPIRE observations and assuming an empirical dust temperature prior of $34^{+13}_{-9}$ K, we provide photometric redshift and far-IR luminosity estimates for 172 SCUBA-2-selected sources within these Planck overdensity fields. The redshift distribution of the sources peak between a redshift of 2 and 4, with one third of the sources having $S_{500}$/$S_{350} > 1$. For the majority of the sources, we find far-IR luminosities of approximately $10^{13}\,\mathrm{L}_\odot$, corresponding to star-formation rates of around $1000$ M$_\odot \mathrm{yr}^{-1}$. For $S_{850}>8$ mJy sources, we show that there is up to an order of magnitude increase in star-formation rate density and an increase in uncorrected number counts of $6$ for $S_{850}>8$ mJy when compared to typical cosmological survey fields. The sources detected with SCUBA-2 account for only approximately $5$ per cent of the Planck flux at 353 GHz, and thus many more fainter sources are expected in these fields.

Measurements of flux density are described for five planets, Mars, Jupiter, Saturn, Uranus, and Neptune, across the six Planck High Frequency Instrument frequency bands (100-857 GHz) and these are then compared with models and existing data. In our analysis, we have also included estimates of the brightness of Jupiter and Saturn at the three frequencies of the Planck Low Frequency Instrument (30, 44, and 70 GHz). The results provide constraints on the intrinsic brightness and the brightness time-variability of these planets. The majority of the planet flux density estimates are limited by systematic errors, but still yield better than 1% measurements in many cases. Applying data from Planck HFI, the Wilkinson Microwave Anisotropy Probe (WMAP), and the Atacama Cosmology Telescope (ACT) to a model that incorporates contributions from Saturn's rings to the planet's total flux density suggests a best fit value for the spectral index of Saturn's ring system of $\beta _\mathrm{ring} = 2.30\pm0.03$ over the 30-1000 GHz frequency range. The average ratio between the Planck-HFI measurements and the adopted model predictions for all five planets (excluding Jupiter observations for 353 GHz) is 0.997, 0.997, 1.018, and 1.032 for 100, 143, 217, and 353 GHz, respectively. Model predictions for planet thermodynamic temperatures are therefore consistent with the absolute calibration of Planck-HFI detectors at about the three-percent-level. We compare our measurements with published results from recent cosmic microwave background experiments. In particular, we observe that the flux densities measured by Planck HFI and WMAP agree to within 2%. These results allow experiments operating in the mm-wavelength range to cross-calibrate against Planck and improve models of radiative transport used in planetary science.

Measurement of the gas velocity distribution in galaxy clusters provides insight into the physics of mergers, through which large scale structures form in the Universe. Velocity estimates within the intracluster medium (ICM) can be obtained via the Sunyaev-Zel'dovich (SZ) effect, but its observation is challenging both in term of sensitivity requirement and control of systematic effects, including the removal of contaminants. In this paper we report resolved observations, at 150 and 260 GHz, of the SZ effect toward the triple merger MACS J0717.5+3745 (z=0.55), using data obtained with the NIKA camera at the IRAM 30m telescope. Assuming that the SZ signal is the sum of a thermal (tSZ) and a kinetic (kSZ) component and by combining the two NIKA bands, we extract for the first time a resolved map of the kSZ signal in a cluster. The kSZ signal is dominated by a dipolar structure that peaks at -5.1 and +3.4 sigma, corresponding to two subclusters moving respectively away and toward us and coincident with the cold dense X-ray core and a hot region undergoing a major merging event. We model the gas electron density and line-of-sight velocity of MACS J0717.5+3745 as four subclusters. Combining NIKA data with X-ray observations from XMM-Newton and Chandra, we fit this model to constrain the gas line-of-sight velocity of each component, and we also derive, for the first time, a velocity map from kSZ data (i.e. that is model-dependent). Our results are consistent with previous constraints on the merger velocities, and thanks to the high angular resolution of our data, we are able to resolve the structure of the gas velocity. Finally, we investigate possible contamination and systematic effects with a special care given to radio and submillimeter galaxies. Among the sources that we detect with NIKA, we find one which is likely to be a high redshift lensed submillimeter galaxy.

This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 lensing map. This cross-correlation yields a detection at $4\,\sigma$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,\sigma$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. The ISW signal allow us to detect $\Omega_\Lambda$ at more than $3\,\sigma$. This cross-correlation analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, polarization data is used to study the anomalously large ISW signal previously reported through the stacking of CMB features at the locations of known superstructures. We find that the current Planck polarization data do not exclude that this signal could be caused by the ISW effect. In addition, the stacking of the Planck lensing map on the locations of superstructures exhibits a positive cross-correlation with these large-scale structures. Finally, we have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the LSS traced by the 2MPZ survey by directly inverting the density field into the gravitational potential field.

We present spatially resolved ALMA [CII] observations of the bright (flux density S=400 mJy at 350 microns), gravitationally lensed, starburst galaxy PLCK G045.1+61.1 at z=3.427, the "Garnet". This source is part of our set of "Planck's Dusty GEMS", discovered with the Planck's all-sky survey. Two emission-line clouds with a relative velocity offset of ~600 km/s extend towards north-east and south-west, respectively, of a small, intensely star-forming clump with a star-formation intensity of 220 Msun/yr/kpc^2, akin to maximal starbursts. [CII] is also seen in absorption, with a redshift of +350 km/s relative to the brightest CO component. [CII] absorption has previously only been found in the Milky Way along sightlines toward bright high-mass star-forming regions, and this is the first detection in another galaxy. Similar to Galactic environments, the [CII] absorption feature is associated with [CI] emission, implying that this is diffuse gas shielded from the UV radiation of the clump, and likely at large distances from the clump. Since absorption can only be seen in front of a continuum source, the gas in this structure can definitely be attributed to gas flowing towards the clump. The absorber could be part of a cosmic filament or merger debris being accreted onto the galaxy. We discuss our results also in light of the on-going debate of the origin of the [CII] deficit in dusty star-forming galaxies.

We present new measurements of the power spectra of the cosmic infrared background (CIB) anisotropies using the Planck 2015 full-mission HFI data at 353, 545, and 857 GHz over 20000 square degrees. We use techniques similar to those applied for the cosmological analysis of Planck, subtracting dust emission at the power spectrum level. Our analysis gives stable solutions for the CIB power spectra with increasing sky coverage up to about 50% of the sky. These spectra agree well with Hi cleaned spectra from Planck measured on much smaller areas of sky with low Galactic dust emission. At 545 and 857 GHz our CIB spectra agree well with those measured from Herschel data. We find that the CIB spectra at l > 500 are well fitted by a power-law model for the clustered CIB, with a shallow index {\gamma}^cib = 0.53\pm0.02. This is consistent with the CIB results at 217 GHz from the cosmological parameter analysis of Planck. We show that a linear combination of the 545 and 857 GHz Planck maps is dominated by CIB fluctuations at multipoles l > 300.

The secondary cosmic microwave background (CMB) $B$-modes stem from the post-decoupling distortion of the polarization $E$-modes due to the gravitational lensing effect of large-scale structures. These lensing-induced $B$-modes constitute both a valuable probe of the dark matter distribution and an important contaminant for the extraction of the primary CMB $B$-modes from inflation. Planck provides accurate nearly all-sky measurements of both the polarization $E$-modes and the integrated mass distribution via the reconstruction of the CMB lensing potential. By combining these two data products, we have produced an all-sky template map of the lensing-induced $B$-modes using a real-space algorithm that minimizes the impact of sky masks. The cross-correlation of this template with an observed (primordial and secondary) $B$-mode map can be used to measure the lensing $B$-mode power spectrum at multipoles up to $2000$. In particular, when cross-correlating with the $B$-mode contribution directly derived from the Planck polarization maps, we obtain lensing-induced $B$-mode power spectrum measurement at a significance level of $12\,\sigma$, which agrees with the theoretical expectation derived from the Planck best-fit $\Lambda$CDM model. This unique nearly all-sky secondary $B$-mode template, which includes the lensing-induced information from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is delivered as part of the Planck 2015 public data release. It will be particularly useful for experiments searching for primordial $B$-modes, such as BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of the lensing-induced contribution to the measured total CMB $B$-modes.

We present the most significant measurement of the cosmic microwave background (CMB) lensing potential to date (at a level of 40 sigma), using temperature and polarization data from the Planck 2015 full-mission release. Using a polarization-only estimator we detect lensing at a significance of 5 sigma. We cross-check the accuracy of our measurement using the wide frequency coverage and complementarity of the temperature and polarization measurements. Public products based on this measurement include an estimate of the lensing potential over approximately 70% of the sky, an estimate of the lensing potential power spectrum in bandpowers for the multipole range 40<L<400 and an associated likelihood for cosmological parameter constraints. We find good agreement between our measurement of the lensing potential power spectrum and that found in the best-fitting LCDM model based on the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level measurement of the parameter combination $\sigma_8 \Omega_m^{0.25} = 0.591\pm 0.021$. We combine our determination of the lensing potential with the E-mode polarization also measured by Planck to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated with the B-modes observed directly by Planck at the expected level and with a statistical significance of 10 sigma, confirming Planck's sensitivity to this known sky signal. We also correlate our lensing potential estimate with the large-scale temperature anisotropies, detecting a cross-correlation at the 3 sigma level, as expected due to dark energy in the concordance LCDM model.

We investigate constraints on cosmic reionization extracted from the Planck cosmic microwave background (CMB) data. We combine the Planck CMB anisotropy data in temperature with the low-multipole polarization data to fit LCDM models with various parameterizations of the reionization history. We obtain a Thomson optical depth tau=0.058 +/- 0.012 for the commonly adopted instantaneous reionization model. This confirms, with only data from CMB anisotropies, the low value suggested by combining Planck 2015 results with other data sets and also reduces the uncertainties. We reconstruct the history of the ionization fraction using either a symmetric or an asymmetric model for the transition between the neutral and ionized phases. To determine better constraints on the duration of the reionization process, we also make use of measurements of the amplitude of the kinetic Sunyaev-Zeldovich (kSZ) effect using additional information from the high resolution Atacama Cosmology Telescope and South Pole Telescope experiments. The average redshift at which reionization occurs is found to lie between z=7.8 and 8.8, depending on the model of reionization adopted. Using kSZ constraints and a redshift-symmetric reionization model, we find an upper limit to the width of the reionization period of Dz < 2.8. In all cases, we find that the Universe is ionized at less than the 10% level at redshifts above z~10. This suggests that an early onset of reionization is strongly disfavoured by the Planck data. We show that this result also reduces the tension between CMB-based analyses and constraints from other astrophysical sources.

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $\Lambda$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5{\sigma} from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $\tau=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($\mu$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $\Lambda$CDM parameters consistent with those from temperature alone.

The characterization of the Galactic foregrounds has been shown to be the main obstacle in the challenging quest to detect primordial B-modes in the polarized microwave sky. We make use of the Planck-HFI 2015 data release at high frequencies to place new constraints on the properties of the polarized thermal dust emission at high Galactic latitudes. Here, we specifically study the spatial variability of the dust polarized spectral energy distribution, and its potential impact on the determination of the tensor-to-scalar ratio. We use the correlation ratio of the $C_\ell^{BB}$ angular power spectra between the 217- and 353-GHz channels as a tracer of these potential variations, computed on different high Galactic latitude regions, ranging from 80% to 20% of the sky. The new insight from Planck data is a departure of the correlation ratio from unity that cannot be attributed to a spurious decorrelation due to the cosmic microwave background, instrumental noise, or instrumental systematics. The effect is marginally detected on each region, but the statistical combination of all the regions gives more than 99% confidence for this variation in polarized dust properties. In addition, we show that the decorrelation increases when there is a decrease in the mean column density of the region of the sky being considered, and we propose a simple power-law empirical model for this dependence, which matches what is seen in the Planck data. We explore the effect that this measured decorrelation has on simulations of the BICEP2-Keck Array/Planck analysis and show that the 2015 constraints from those data still allow a decorrelation between the dust at 150 and 353GHz of the order of the one we measure. Finally we show that either spatial variation of the dust SED or of the dust polarization angle could produce decorrelations between 217- and 353-GHz data similar to those we observe in the data.

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The radio spectra peak at high frequencies and exhibit a variety of shapes. For a small set of low-z sources, we find a spectral upturn at high frequencies, indicating the presence of intrinsic cold dust. Variability can generally be approximated by achromatic variations, while sources with clear signatures of evolving shocks appear to be limited to the strongest outbursts.