The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works.

Characterizing the variability of the extragalactic sources used for calibration in the Atacama Large Millimeter/Sub-millimeter Array (ALMA) is key to assess the flux scale uncertainty of science observations. To this end, we model the variability of 39 quasars which have been used by ALMA as secondary flux calibrators using continuous time stochastic processes. This formalism is specially adapted to the multi-frequency, quasi-periodic sampling which characterizes the calibration monitoring of ALMA. We find that simple mixtures of Ornstein-Uhlenbeck processes can describe well the flux and spectral index variability of these sources for Bands 3 and 7 (91.5 and 103.5, and 343.5 GHz, respectively). The spectral shape of the calibrators are characterized by negative spectral indices, mostly between $-0.35$ and $-0.80$, and with additional concavity. The model provides forecasts, interpolations, and uncertainty estimations for the observed fluxes that depend on the intrinsic variability of the source. These can be of practical use for the ALMA data calibrator survey and data quality assurance.

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 Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

The bright [CI] 1-0 and [CI] 2-1 lines of atomic carbon are becoming more and more widely employed tracers of the cold neutral gas in high-redshift galaxies. Here we present observations of these lines in the 11 galaxies of the set of Planck's Dusty GEMS, the brightest gravitationally lensed galaxies on the extragalactic submillimeter sky probed by the Planck satellite. We have [CI] 1-0 measurements for seven, and [CI] 2-1 measurements for eight galaxies, including four galaxies where both lines are measured. We use our observations to constrain the gas excitation mechanism, excitation temperatures, optical depths, atomic carbon and molecular gas masses, and carbon abundances. Ratios of L_CI/L_ FIR are similar to those found in the local Universe, and suggest that the total cooling budget through atomic carbon has not strongly changed in the last 12 Gyr. Both lines are optically thin and trace 1 - 6 x 10^7 M_sun of atomic carbon. Carbon abundance ratios with H_2, X_CI, are between 2.5 and 4 x 10^-5, for a "ULIRG" CO-to-H_2 conversion factor of alpha_CO=0.8 M_sun/ [K km s^-1 pc^2]. Ratios of molecular gas masses derived from [CI] 1-0 and CO agree within the measurement uncertainties for five galaxies, and to better than a factor of 2 for another two with [CI] 1-0 measurements, after taking CO excitation carefully into account. This does not support the idea that intense, high-redshift starburst galaxies host large quantities of "CO-dark" gas. These results also support the common assumptions underlying most molecular gas mass estimates made for massive, dusty, high-redshift starburst galaxies, although the good agreement between the masses obtained with both tracers cannot be taken as an independent confirmation of either alpha_CO or X_CI.

We present an extensive CO emission-line survey of the Planck's dusty GEMS, a small set of 11 strongly lensed dusty star-forming galaxies at z = 2-4 discovered with Planck and Herschel satellites, using EMIR on the IRAM 30-m telescope. We detected a total of 45 CO rotational lines from Jup=3 to Jup=11, and up to eight transitions per source, allowing a detailed analysis of the gas excitation and interstellar medium conditions within these extremely bright, vigorous starbursts. We applied radiative transfer models using the large velocity gradient approach to infer the spatially-averaged molecular gas densities, $n_{H_2}$~10$^{2.6-4.1}$ cm$^{-3}$, and kinetic temperatures, $T_k$~30-1000 K. In five sources, we find evidence of two distinct gas phases with different properties and model their CO ladder with two excitation components. The warm (70-320 K) and dense gas reservoirs in these galaxies are highly excited, while the cooler (15-60 K) and more extended low-excitation components cover a range of gas densities. In two sources, the latter is associated with diffuse Milky Way-like gas phases, which provides evidence that a significant fraction of the total gas masses of dusty starburst galaxies can be embedded in cool, low-density reservoirs. Finally, we show that the CO line luminosity ratios are consistent with those predicted by models of photon-dominated regions and disfavor scenarios of gas clouds irradiated by intense X-ray fields from active galactic nuclei. By combining CO, [CI] and [CII] line diagnostics, we obtain average PDR gas densities significantly higher than in normal star-forming galaxies at low-redshift, and far-ultraviolet radiation fields at least 100 times more intense than in the Milky Way. These spatially-averaged conditions are consistent with those in high-redshift SMGs and in a range of low-redshift environments.

We report the discovery of a molecular wind signature from a massive intensely star-forming clump of a few $10^9$ Msun, in the strongly gravitationally lensed submillimeter galaxy "the Emerald" (PLCK_G165.7+49.0) at z=2.236. The Emerald is amongst the brightest high-redshift galaxies on the submillimeter sky, and was initially discovered with the Planck satellite. The system contains two magnificient structures with projected lengths of 28.5" and 21" formed by multiple, near-infrared arcs, falling behind a massive galaxy cluster at z=0.35, as well as an adjacent filament that has so far escaped discovery in other wavebands. We used HST/WFC3 and CFHT optical and near-infrared imaging together with IRAM and SMA interferometry of the CO(4-3) line and 850 $\mu$m dust emission to characterize the foreground lensing mass distribution, construct a lens model with Lenstool, and calculate gravitational magnification factors between 20 and 50 in most of the source. The majority of the star formation takes place within two massive star-forming clumps which are marginally gravitationally bound and embedded in a $9 \times 10^{10}$ Msun, fragmented disk with 20% gas fraction. One of the clumps shows a pronounced blue wing in the CO(4-3) line profile, which we interpret as a wind signature. The mass outflow rates are high enough for us to suspect that the clump might become unbound within a few tens of Myr, unless the outflowing gas can be replenished by gas accretion from the surrounding disk. The velocity offset of -200 km s$^{-1}$ is above the escape velocity of the clump, but not that of the galaxy overall, suggesting that much of this material might ultimately rain back onto the galaxy and contribute to fueling subsequent star formation.

[abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, enhancing their measured emission while increasing spatial resolution. Derived number counts depend, however, on mass reconstruction models that properly constrain these clusters. We estimate the 1.1 mm number counts along the line of sight of three galaxy clusters, i.e. Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223, which are part of the ALMA Frontier Fields Survey. We perform detailed simulations to correct these counts for lensing effects. We use several publicly available lensing models for the galaxy clusters to derive the intrinsic flux densities of our sources. We perform Monte Carlo simulations of the number counts for a detailed treatment of the uncertainties in the magnifications and adopted source redshifts. We find an overall agreement among the number counts derived for the different lens models, despite their systematic variations regarding source magnifications and effective areas. Our number counts span ~2.5 dex in demagnified flux density, from several mJy down to tens of uJy. Our number counts are consistent with recent estimates from deep ALMA observations at a 3$\sigma$ level. Below $\approx$ 0.1 mJy, however, our cumulative counts are lower by $\approx$ 1 dex, suggesting a flattening in the number counts. In our deepest ALMA mosaic, we estimate number counts for intrinsic flux densities $\approx$ 4 times fainter than the rms level. This highlights the potential of probing the sub-10 uJy population in larger samples of galaxy cluster fields with deeper ALMA observations.

The use of galaxy clusters as precision cosmological probes relies on an accurate determination of their masses. However, inferring the relationship between cluster mass and observables from direct observations is difficult and prone to sample selection biases. In this work, we use weak lensing as the best possible proxy for cluster mass to calibrate the Sunyaev-Zel'dovich (SZ) effect measurements from the APEX-SZ experiment. For a well-defined (ROSAT) X-ray complete cluster sample, we calibrate the integrated Comptonization parameter, $Y_{\rm SZ}$, to the weak-lensing derived total cluster mass, $M_{500}$. We employ a novel Bayesian approach to account for the selection effects by jointly fitting both the SZ Comptonization, $Y_{\rm SZ}\text{--}M_{500}$, and the X-ray luminosity, $L_{\rm x}\text{--}M_{500}$, scaling relations. We also account for a possible correlation between the intrinsic (log-normal) scatter of $L_{\rm x}$ and $Y_{\rm SZ}$ at fixed mass. We find the corresponding correlation coefficient to be $r= 0.47_{-0.35}^{+0.24}$, and at the current precision level our constraints on the scaling relations are consistent with previous works. For our APEX-SZ sample, we find that ignoring the covariance between the SZ and X-ray observables biases the normalization of the $Y_{\rm SZ}\text{--}M_{500}$ scaling high by $1\text{--}2\sigma$ and the slope low by $\sim 1\sigma$, even when the SZ effect plays no role in the sample selection. We conclude that for higher-precision data and larger cluster samples, as anticipated from on-going and near-future cluster cosmology experiments, similar biases (due to intrinsic covariances of cluster observables) in the scaling relations will dominate the cosmological error budget if not accounted for correctly.

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.

Most sub-mm emission line studies of galaxies to date have targeted sources with known redshifts where the frequencies of the lines are well constrained. Recent blind line scans circumvent the spectroscopic redshift requirement, which could represent a selection bias. Our aim is to detect emission lines present in continuum oriented observations. The detection of such lines provides spectroscopic redshift and yields properties of the galaxies. We perform a search for emission lines in the ALMA observations of five Frontier Fields clusters and assess the reliability of our detection by associating line candidates with detected galaxies in deep near-infrared imaging. We find 26 significant emission lines candidates, with observed line fluxes between 0.2-4.6 Jy km s$^{-1}$ and velocity dispersions (FWHM) of 25-600 km s$^{-1}$. Nine of these candidates lie nearby to near-infrared sources, boosting their reliability; in six cases the observed line frequency and strength are consistent with expectations given the photometric redshift and properties of the galaxy counterparts. We present redshift identifications, magnifications and molecular gas estimates for the galaxies with identified lines. We show that two of these candidates likely originate from starburst galaxies, one of which is a jellyfish galaxy, while another two are consistent with being main sequence galaxies based in their depletion times. This work highlights the degree to which serendipitous emission lines can be discovered in large mosaic continuum observations when deep ancillary data are available. The low number of high-significance line detections, however, confirms that such surveys are not as optimal as blind line scans. We stress that Monte Carlo simulations should be used to assess the line detections significances, since using the negative noise suffers from stochasticity and incurs larger uncertainties.

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.

We present an analysis of high-resolution ALMA interferometry of CO(4-3) line emission and dust continuum in the "Ruby" (PLCK_G244.8+54.9), a bright, gravitationally lensed galaxy at z = 3.0 discovered with the Planck all-sky survey. The Ruby is the brightest of Planck's Dusty GEMS, a sample of 11 of the brightest gravitationally lensed high-redshift galaxies on the extragalactic sub-mm sky. We resolve the high-surface-brightness continuum and CO line emission of the Ruby in several extended clumps along a partial, nearly circular Einstein ring with 1.4" diameter around a massive galaxy at z = 1.5. Local star-formation intensities are up to 2000 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, amongst the highest observed at high redshift, and clearly in the range of maximal starbursts. Gas-mass surface densities are a few $\times$ 10$^4$ M$_{\odot}$ pc$^{-2}$. The Ruby lies at, and in part even above, the starburst sequence in the Schmidt-Kennicutt diagram, and at the limit expected for star formation that is self-regulated through the kinetic energy injection from radiation pressure, stellar winds, and supernovae. We show that these processes can also inject sufficient kinetic energy and momentum into the gas to explain the turbulent line widths, which are consistent with marginally gravitationally bound molecular clouds embedded in a critically Toomre-stable disk. The star-formation efficiency is in the range 1-10% per free-fall time, consistent with the notion that the pressure balance that sets the local star-formation law in the Milky Way may well be universal out to the highest star-formation intensities. AGN feedback is not necessary to regulate the star formation in the Ruby, in agreement with the absence of a bright AGN component in the infrared and radio regimes.

[abridged] The Hubble and Spitzer Space Telescope surveys of the Frontier Fields (FF) provide extremely deep images around six massive, strong-lensing clusters of galaxies. The ALMA FF survey aims to cover the same fields at 1.1mm, with maps reaching (unlensed) sensitivities of $<$70$\mu$Jy, in order to explore the properties of background dusty star-forming galaxies. We report on the multi-wavelength photometric analysis of all 12 significantly detected ($>$5$\sigma$) sources in the first three FF clusters observed by ALMA, based on data from Hubble and Spitzer, the VLT and Herschel. We measure the total photometry in all available bands and determine the photometric redshifts and the physical properties of the counterparts via SED-fitting. In particular, we carefully estimate the FIR photometry using 1.1mm priors to limit the misidentification of blended FIR counterparts, which strongly affect some flux estimates in previous FIR catalogs. We identify robust near-infrared (NIR) counterparts for all 11 sources with K$_s$ detection, the majority of which are quite red, with eight having $F814W-K_s\gtrsim 4$ and five having $F160W-[4.5]\gtrsim3$. From the FIR point of view, all our objects have $z_{phot}$$\sim$1--3, whereas based on the optical SED one object prefers a high-$z$ solution ($z\geq\ $7). Five objects among our sample have spectroscopic redshifts from the GLASS survey for which we can reproduce their SEDs with existing templates. This verification confirms the validity of our photometric redshift methodology. The mean redshift of our sample is $z_{phot}$=1.99$\pm$0.27. All 1.1mm selected objects are massive (10.0$<\log[M_{\star}(M_{\odot})]<$ 11.5), with high star formation rates ($<\log[SFR(M_{\odot}/yr)]> \approx$1.6) and high dust contents (8.1 $<\log[M_{dust} (M_{\odot})]<$8.8), consistent with previous ALMA surveys.

We study the properties of the foreground galaxy of the Ruby, the brightest gravitationally lensed high-redshift galaxy on the sub-millimeter sky as probed by the Planck satellite, and part of our sample of Planck's Dusty GEMS. The Ruby consists of an Einstein ring of 1.4" diameter at z = 3.005 observed with ALMA at 0.1" resolution, centered on a faint, red, massive lensing galaxy seen with HST/WFC3, which itself has an exceptionally high redshift, z = 1.525 $\pm$ 0.001, as confirmed with VLT/X-Shooter spectroscopy. Here we focus on the properties of the lens and the lensing model obtained with LENSTOOL. The rest-frame optical morphology of this system is strongly dominated by the lens, while the Ruby itself is highly obscured, and contributes less than 10% to the photometry out to the K band. The foreground galaxy has a lensing mass of (3.70 $\pm$ 0.35) $\times$ 10$^{11}$ M$_{\odot}$. Magnification factors are between 7 and 38 for individual clumps forming two image families along the Einstein ring. We present a decomposition of the foreground and background sources in the WFC3 images, and stellar population synthesis modeling with a range of star-formation histories for Chabrier and Salpeter initial mass functions (IMFs). Only the stellar mass range obtained with the latter agrees well with the lensing mass. This is consistent with the bottom-heavy IMFs of massive high-redshift galaxies expected from detailed studies of the stellar masses and mass profiles of their low-redshift descendants, and from models of turbulent gas fragmentation. This may be the first direct constraint on the IMF in a lens at z = 1.5, which is not a cluster central galaxy.

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.

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.

Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion.

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

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

Although infrared (IR) overall dust emission from clusters of galaxies has been statistically detected using data from the Infrared Astronomical Satellite (IRAS), it has not been possible to sample the spectral energy distribution (SED) of this emission over its peak, and thus to break the degeneracy between dust temperature and mass. By complementing the IRAS spectral coverage with Planck satellite data from 100 to 857 GHz, we provide new constraints on the IR spectrum of thermal dust emission in clusters of galaxies. We achieve this by using a stacking approach for a sample of several hundred objects from the Planck cluster sample; this procedure averages out fluctuations from the IR sky, allowing us to reach a significant detection of the faint cluster contribution. We also use the large frequency range probed by Planck, together with component-separation techniques, to remove the contamination from both cosmic microwave background anisotropies and the thermal Sunyaev-Zeldovich effect (tSZ) signal, which dominate below 353 GHz. By excluding dominant spurious signals or systematic effects, averaged detections are reported at frequencies between 353 and 5000 GHz. We confirm the presence of dust in clusters of galaxies at low and intermediate redshifts, yielding an SED with a shape similar to that of the Milky Way. Planck's beam does not allow us to investigate the detailed spatial distribution of this emission (e.g., whether it comes from intergalactic dust or simply the dust content of the cluster galaxies), but the radial distribution of the emission appears to follow that of the stacked SZ signal, and thus the extent of the clusters. The recovered SED allows us to constrain the dust mass responsible for the signal, as well as its temperature. We additionally explore the evolution of the IR emission as a function of cluster mass and redshift.

The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two sub-catalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these catalogues covers most of the sky and allows the user to produce subsamples at higher reliabilities than the target 80% integral reliability of the catalogue. The PCCS2E contains sources detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and ERCSC catalogues.

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

The all-sky Planck survey in 9 frequency bands was used to search for emission from all 274 known Galactic supernova remnants. Of these, 16 were detected in at least two Planck frequencies. The radio-through-microwave spectral energy distributions were compiled to determine the emission mechanism for microwave emission. In only one case, IC 443, is there high-frequency emission clearly from dust associated with the supernova remnant.In all cases, the low-frequency emission is from synchrotron radiation. A single power law, as predicted for a population of relativistic particles with energy distribution that extends continuously to high energies, is evident for many sources, including the Crab and PKS 1209-51/52. A decrease in flux density relative to the extrapolation of radio emission is evident in several sources. Their spectral energy distributions can be approximated as broken power laws, $S_\nu\propto\nu^{-\alpha}$, with the spectral index, alpha, increasing by 0.5-1 above a break frequency in the range 10-60 GHz. The break could be due to synchrotron losses.

The Virgo cluster is the largest Sunyaev-Zeldovich (SZ) source in the sky, both in terms of angular size and total integrated flux. Planck's wide angular scale and frequency coverage, together with its high sensitivity, allow a detailed study of this large object through the SZ effect. Virgo is well resolved by Planck, showing an elongated structure, which correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find a good agreement between the SZ signal (or Compton paranmeter, y_c) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Due to its proximity to us, the gas beyond the virial radius can be studied with unprecedented sensitivity by integrating the SZ signal over tens of square degrees. We study the signal in the outskirts of Virgo and compare it with analytical models and a constrained simulation of the environment of Virgo. Planck data suggest that significant amounts of low-density plasma surround Virgo out to twice the virial radius. We find the SZ signal in the outskirts of Virgo to be consistent with a simple model that extrapolates the inferred pressure at lower radii while assuming that the temperature stays in the keV range beyond the virial radius. The observed signal is also consistent with simulations and points to a shallow pressure profile in the outskirts of the cluster. This reservoir of gas at large radii can be linked with the hottest phase of the elusive warm/hot intergalactic medium. Taking the lack of symmetry of Virgo into account, we find that a prolate model is favoured by the combination of SZ and X-ray data, in agreement with predictions.

By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present a significant detection of baryons participating in large-scale bulk flows around central galaxies (CGs) at redshift $z\approx 0.1$. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the CGC (Central Galaxy Catalogue) samples extracted from Sloan Digital Sky Survey (DR7) data. For the foreground-cleaned maps, we find $1.8$-$2.5\sigma$ detections of the kSZ signal, which are consistent with the kSZ evidence found in individual Planck raw frequency maps, although lower than found in the WMAP-9yr W band ($3.3\sigma$). We further reconstruct the peculiar velocity field from the CG density field, and compute for the first time the cross-correlation function between kSZ temperature fluctuations and estimates of CG radial peculiar velocities. This correlation function yields a $3.0$-$3.7$$\sigma$ detection of the peculiar motion of extended gas on Mpc scales, in flows correlated up to distances of 80-100 $h^{-1}$ Mpc. Both the pairwise momentum estimates and kSZ temperature-velocity field correlation find evidence for kSZ signatures out to apertures of 8 arcmin and beyond, corresponding to a physical radius of $> 1$ Mpc, more than twice the mean virial radius of halos. This is consistent with the predictions from hydro simulations that most of the baryons are outside the virialized halos. We fit a simple model, in which the temperature-velocity cross-correlation is proportional to the signal seen in a semi-analytic model built upon N-body simulations, and interpret the proportionality constant as an "effective" optical depth to Thomson scattering. We find $\tau_T=(1.4\pm0.5)\times 10^{-4}$; the simplest interpretation of this measurement is that much of the gas is in a diffuse phase, which contributes little signal to X-ray or thermal SZ observations.

The quest for a B-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between E-modes and B-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder, we identify 259 filaments at high Galactic latitude, with lengths larger or equal to 2\deg\ (corresponding to 3.5 pc in length for a typical distance of 100 pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes I, Q, U, E, and B, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be 11 %. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the E and B asymmetry and the $C_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck 353 GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter.

We use Planck data to detect the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) effect and the infrared emission from the galaxies that make up the the cosmic infrared background (CIB). We first perform a stacking analysis towards Planck-confirmed galaxy clusters. We detect infrared emission produced by dusty galaxies inside these clusters and demonstrate that the infrared emission is about 50% more extended than the tSZ effect. Modelling the emission with a Navarro--Frenk--White profile, we find that the radial profile concentration parameter is $c_{500} = 1.00^{+0.18}_{-0.15}$. This indicates that infrared galaxies in the outskirts of clusters have higher infrared flux than cluster-core galaxies. We also study the cross-correlation between tSZ and CIB anisotropies, following three alternative approaches based on power spectrum analyses: (i) using a catalogue of confirmed clusters detected in Planck data; (ii) using an all-sky tSZ map built from Planck frequency maps; and (iii) using cross-spectra between Planck frequency maps. With the three different methods, we detect the tSZ-CIB cross-power spectrum at significance levels of (i) 6 $\sigma$, (ii) 3 $\sigma$, and (iii) 4 $\sigma$. We model the tSZ-CIB cross-correlation signature and compare predictions with the measurements. The amplitude of the cross-correlation relative to the fiducial model is $A_{\rm tSZ-CIB}= 1.2\pm0.3$. This result is consistent with predictions for the tSZ-CIB cross-correlation assuming the best-fit cosmological model from Planck 2015 results along with the tSZ and CIB scaling relations.

We present the 8th Full Focal Plane simulation set (FFP8), deployed in support of the Planck 2015 results. FFP8 consists of 10 fiducial mission realizations reduced to 18144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising $10^4$ mission realizations reduced to about $10^6$ maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects; remaining subdominant effects will be included in future updates. Generated at a cost of some 25 million CPU-hours spread across multiple high-performance-computing (HPC) platforms, FFP8 is used for the validation and verification of analysis algorithms, as well as their implementations, and for removing biases from and quantifying uncertainties in the results of analyses of the real data.

The Andromeda Galaxy (M31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiral arms and sub-features. We examine the morphology of this long-wavelength dust emission as seen by Planck, including a study of its outermost spiral arms, and investigate the dust heating mechanism across M31. We find that dust dominating the longer wavelength emission ($\gtrsim 0.3\,$mm) is heated by the diffuse stellar population (as traced by 3.6$\,\mu$m emission), with the dust dominating the shorter wavelength emission heated by a mix of the old stellar population and star-forming regions (as traced by 24$\,\mu$m emission). We also fit spectral energy distributions (SEDs) for individual 5' pixels and quantify the dust properties across the galaxy, taking into account these different heating mechanisms, finding that there is a linear decrease in temperature with galactocentric distance for dust heated by the old stellar population, as would be expected, with temperatures ranging from around 22$\,$K in the nucleus to 14$\,$K outside of the 10$\,$kpc ring. Finally, we measure the integrated spectrum of the whole galaxy, which we find to be well-fitted with a global dust temperature of ($18.2\pm1.0$)$\,$K with a spectral index of $1.62\pm0.11$ (assuming a single modified blackbody), and a significant amount of free-free emission at intermediate frequencies of 20-60$\,$GHz, which corresponds to a star formation rate of around $0.12$M$_\odot\,$yr$^{-1}$. We find a $2.3\,\sigma$ detection of the presence of spinning dust emission, with a 30$\,$GHz amplitude of $0.7\pm0.3\,$Jy, which is in line with expectations from our Galaxy.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

The Planck High Frequency Instrument (HFI) has observed the full sky at six frequencies (100, 143, 217, 353, 545, and 857 GHz) in intensity and at four frequencies in linear polarization (100, 143, 217, and 353 GHz). In order to obtain sky maps, the time-ordered information (TOI) containing the detector and pointing samples must be processed and the angular response must be assessed. The full mission TOI is included in the Planck 2015 release. This paper describes the HFI TOI and beam processing for the 2015 release. HFI calibration and map-making are described in a companion paper. The main pipeline has been modified since the last release (2013 nominal mission in intensity only), by including a correction for the non-linearity of the warm readout and by improving the model of the bolometer time response. The beam processing is an essential tool that derives the angular response used in all the Planck science papers and we report an improvement in the effective beam window function uncertainty of more than a factor 10 relative to the 2013 release. Noise correlations introduced by pipeline filtering function are assessed using dedicated simulations. Angular cross-power spectra using datasets that are decorrelated in time are immune to the main systematic effects.

This paper describes the processing applied to the Planck High Frequency Instrument (HFI) cleaned, time-ordered information to produce photometrically calibrated maps in temperature and (for the first time) in polarization. The data from the entire 2.5 year HFI mission include almost five independent full-sky surveys. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To obtain the best accuracy on the calibration over such a large range, two different photometric calibration schemes have been used. The 545 and 857 GHz data are calibrated using models of planetary atmospheric emission. The lower frequencies (from 100 to 353 GHz) are calibrated using the time-variable cosmological microwave background dipole, which we call the "orbital dipole". This source of calibration only depends on the satellite velocity with respect to the solar system. Using a CMB temperature of 2.7255 +/- 0.0006 K, it permits an independent measurement of the amplitude of the CMB solar dipole (3364.3 +/- 1.5 \mu K) which is approximatively 1\sigma\ higher than the WMAP measurement with a direction that is consistent between both experiments. We describe the pipeline used to produce the maps of intensity and linear polarization from the HFI timelines, and the scheme used to set the zero level of the maps a posteriori. We also summarize the noise characteristics of the HFI maps in the 2015 Planck data release and present some null tests to assess their quality. Finally, we discuss the major systematic effects and in particular the leakage induced by flux mismatch between the detectors that leads to spurious polarization signal.

We present observations and analysis of a sample of 123 galaxy clusters from the 2013 Planck catalogue of Sunyaev-Zel'dovich sources with the Arcminute Microkelvin Imager (AMI), a ground-based radio interferometer. AMI provides an independent measurement with higher angular resolution, 3 arcmin compared to the Planck beams of 5-10 arcmin. The AMI observations thus provide validation of the cluster detections, improved positional estimates, and a consistency check on the fitted 'size' ($\theta_{s}$) and 'flux' ($Y_{\rm tot}$) parameters in the Generalised Navarro, Frenk and White (GNFW) model. We detect 99 of the clusters. We use the AMI positional estimates to check the positional estimates and error-bars produced by the Planck algorithms PowellSnakes and MMF3. We find that $Y_{\rm tot}$ values as measured by AMI are biased downwards with respect to the Planck constraints, especially for high Planck-SNR clusters. We perform simulations to show that this can be explained by deviation from the 'universal' pressure profile shape used to model the clusters. We show that AMI data can constrain the $\alpha$ and $\beta$ parameters describing the shape of the profile in the GNFW model for individual clusters provided careful attention is paid to the degeneracies between parameters, but one requires information on a wider range of angular scales than are present in AMI data alone to correctly constrain all parameters simultaneously.

We present an analysis of 11 bright far-IR/submm sources discovered through a combination of the Planck survey and follow-up Herschel-SPIRE imaging. Each source has a redshift z=2.2-3.6 obtained through a blind redshift search with EMIR at the IRAM 30-m telescope. Interferometry obtained at IRAM and the SMA, and optical/near-infrared imaging obtained at the CFHT and the VLT reveal morphologies consistent with strongly gravitationally lensed sources. Additional photometry was obtained with JCMT/SCUBA-2 and IRAM/GISMO at 850 um and 2 mm, respectively. All objects are bright, isolated point sources in the 18 arcsec beam of SPIRE at 250 um, with spectral energy distributions peaking either near the 350 um or the 500 um bands of SPIRE, and with apparent far-infrared luminosities of up to 3x10^14 L_sun. Their morphologies and sizes, CO line widths and luminosities, dust temperatures, and far-infrared luminosities provide additional empirical evidence that these are strongly gravitationally lensed high-redshift galaxies. We discuss their dust masses and temperatures, and use additional WISE 22-um photometry and template fitting to rule out a significant contribution of AGN heating to the total infrared luminosity. Six sources are detected in FIRST at 1.4 GHz. Four have flux densities brighter than expected from the local far-infrared-radio correlation, but in the range previously found for high-z submm galaxies, one has a deficit of FIR emission, and 6 are consistent with the local correlation. The global dust-to-gas ratios and star-formation efficiencies of our sources are predominantly in the range expected from massive, metal-rich, intense, high-redshift starbursts. An extensive multi-wavelength follow-up programme is being carried out to further characterize these sources and the intense star-formation within them.

Planck data when combined with ancillary data provide a unique opportunity to separate the diffuse emission components of the inner Galaxy. The purpose of the paper is to elucidate the morphology of the various emission components in the strong star-formation region lying inside the solar radius and to clarify the relationship between the various components. The region of the Galactic plane covered is l=300-0-60deg where star-formation is highest and the emission is strong enough to make meaningful component separation. The latitude widths in this longitude range lie between 1deg and 2deg, which correspond to FWHM z-widths of 100-200pc at a typical distance of 6kpc. The four emission components studied here are synchrotron, free-free, anomalous microwave emission (AME), and thermal (vibrational) dust emission. These components are identified by constructing spectral energy distributions (SEDs) at positions along the Galactic plane using the wide frequency coverage of Planck (28.4-857GHz) in combination with low-frequency radio data at 0.408-2.3GHz plus WMAP data at 23-94GHz, along with far-infrared (FIR) data from DIRBE and IRAS. The free-free component is determined from radio recombination line (RRL) data. AME is found to be comparable in brightness to the free-free emission on the Galactic plane in the frequency range 20-40GHz with a width in latitude similar to that of the thermal dust; it comprises 45+/-1% of the total 28.4GHz emission in the longitude range l=300-0-60deg. The free-free component is the narrowest, reflecting the fact that it is produced by current star-formation as traced by the narrow distribution of OB stars. It is the dominant emission on the plane between 60 and 100GHz. RRLs from this ionized gas are used to assess its distance, leading to a free-free z-width of FWHM ~100pc...(abridged)

Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We make use of the Planck 353 GHz I, Q, and U Stokes maps as dust templates, and cross-correlate them with the Planck and WMAP data at 12 frequencies from 23 to 353 GHz, over circular patches with 10 degree radius. The cross-correlation analysis is performed for both intensity and polarization data in a consistent manner. We use a mask that focuses our analysis on the diffuse interstellar medium at intermediate Galactic latitudes. We determine the spectral indices of dust emission in intensity and polarization between 100 and 353 GHz, for each sky-patch. The mean values, $1.59\pm0.02$ for polarization and $1.51\pm0.01$ for intensity, for a mean dust temperature of 19.6 K, are close, but significantly different ($3.6\,\sigma$). We determine the mean spectral energy distribution (SED) of the microwave emission, correlated with the 353 GHz dust templates, by averaging the results of the correlation over all sky-patches. We find that the mean SED increases for decreasing frequencies at $\nu < 60$ GHz, for both intensity and polarization. The rise of the polarization SED towards low frequencies may be accounted for by a synchrotron component correlated with dust, with no need for any polarization of the anomalous microwave emission. We use a spectral model to separate the synchrotron and dust polarization and to characterize the spectral dependence of the dust polarization fraction. The polarization fraction ($p$) of the dust emission decreases by $(21\pm6)$ % from 353 to 70 GHz. The decrease of $p$ could indicate differences in polarization efficiency among components of interstellar dust (e.g., carbon versus silicate grains). Our observational results provide inputs to quantify and optimize the separation between Galactic and cosmological polarization.

We present thermal Sunyaev-Zel'dovich effect (SZE) measurements for 42 galaxy clusters observed at 150 GHz with the APEX-SZ experiment. For each cluster, we model the pressure profile and calculate the integrated Comptonization $Y$ to estimate the total thermal energy of the intracluster medium (ICM). We compare the measured $Y$ values to X-ray observables of the ICM from the literature (cluster gas mass $M_{\rm{gas}}$, temperature $T_X$, and $Y_X =M_{\rm{gas}}T_X$) that relate to total cluster mass. We measure power law scaling relations, including an intrinsic scatter, between the SZE and X-ray observables for three subsamples within the set of 42 clusters that have uniform X-ray analysis in the literature. We observe that differences between these X-ray analyses introduce significant variability into the measured scaling relations, particularly affecting the normalization. For all three subsamples, we find results consistent with a self-similar model of cluster evolution dominated by gravitational effects. Comparing to predictions from numerical simulations, these scaling relations prefer models that include cooling and feedback in the ICM. Lastly, we measure an intrinsic scatter of $\sim28$ per cent in the $Y-Y_X\,$ scaling relation for all three subsamples.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5\deg$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $\mu$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$\Lambda$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,\sigma$ significance.

[Abridged] We use the Planck all-sky submm and mm maps to search for rare sources distinguished by extreme brightness, a few hundreds of mJy, and their potential for being situated at high redshift. These "cold" Planck sources, selected using the High Frequency Instrument (HFI) directly from the maps and from the Planck Catalogue of Compact Sources (PCCS), all satisfy the criterion of having their rest-frame far-infrared peak redshifted to the frequency range 353 and 857 GHz. This colour-selection favours galaxies in the redshift range z=2-4, which we consider as cold peaks in the cosmic infrared background (CIB). We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500um sources, in comparison to reference SPIRE fields. About 94% of the SPIRE sources in the Planck fields are consistent with being overdensities of galaxies peaking at 350um. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z>2.2. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10sigma. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z~2. Under the Td=35K assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4x10^12 Lsun, yielding star formation rates of typically 700 Msun.yr^-1. If the observed overdensities are actual gravitationally-bound structures, the total total star formation rates reaches 7x10^3 Msun.yr^-1. Taken together, these sources show the signatures of high-z (z>$) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

The role of the magnetic field in the formation of the filamentary structures observed in the interstellar medium (ISM) is a debated topic. The Planck all-sky maps of linearly polarized emission from dust at 353GHz provide the required combination of imaging and statistics to study the correlation between the structures of the Galactic magnetic field and of interstellar matter, both in the diffuse ISM and in molecular clouds. The data reveal structures, or ridges, in the intensity map with counterparts in the Stokes Q and/or U maps. We focus on structures at intermediate and high Galactic latitudes with column density from $10^{20}$ to $10^{22}$ cm$^{-2}$. We measure the magnetic field orientation on the plane of the sky from the polarization data, and present an algorithm to estimate the orientation of the ridges from the dust intensity map. We use analytical models to account for projection effects. Comparing polarization angles on and off the structures, we estimate the mean ratio between the strengths of the turbulent and mean components of the magnetic field to be between 0.6 and 1.0, with a preferred value of 0.8. We find that the ridges are preferentially aligned with the magnetic field measured on the structures. This trend becomes more striking for increasing polarization fraction and decreasing column density. We interpret the increase of alignment with polarization fraction as a consequence of projections effects. The decrease of alignment for high column density is not due to a loss of correlation between the structures and the geometry of the magnetic field. In molecular complexes, we observe structures perpendicular to the magnetic field, which cannot be accounted for by projection effects. We discuss our results in the context of models and MHD simulations, which describe the formation of structures in the magnetized ISM.

Shortened abstract: Observations of the nearby Chamaeleon clouds in gamma rays with the Fermi Large Area Telescope and in thermal dust emission with Planck and IRAS have been used with the HI and CO radio data to (i) map the gas column densities in the different phases and at the dark neutral medium (DNM) transition between the HI-bright and CO-bright media; (ii) constrain the CO-to-$H_2$ conversion factor, $X_{CO}$; (iii) probe the dust properties per gas nucleon in each gas phase and spatially across the clouds. We have separated clouds in velocity in HI and CO emission and modelled the 0.4-100 GeV intensity, the dust optical depth at 353 GHz, the thermal radiance of the large grains, and an estimate of the dust extinction empirically corrected for the starlight intensity, $A_{VQ}$. The gamma-ray emissivity spectra confirm that the GeV-TeV cosmic rays uniformly permeate all gas phases up to the CO cores. The dust and cosmic rays reveal large amounts of DNM gas, with comparable spatial distributions and twice as much mass as in the CO-bright clouds. We give constraints on the HI-DNM-CO transitions and CO-dark $H_2$ fractions for separate clouds. The corrected extinction provides the best fit to the total gas traced by the gamma rays, but we find evidence for a rise in $A_{VQ}/N_H$ and a steep rise in opacity, with increasing $N_H$ and $H_2$ fraction, and with decreasing dust temperature. We observe less variations for the specific power of the grains, except for a decline by half in the CO cores. This combined information suggests grain evolution. The gamma rays and dust radiance yield consistent $X_{CO}$ estimates near $0.7\times10^{20}$ cm$^{-2}$ (K km/s)$^{-1}$. The other dust tracers yield biased values because of the grain opacity rise in the CO clouds. These results also confirm a factor of 2 difference between $X_{CO}$ estimates at pc and kpc scales.

The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust, bringing new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, $p_V$ , and the optical depth in the V band to the star, $\tau_V$. Toward these stars we measure the submillimetre polarized intensity, $P_S$, and total intensity, $I_S$, in the Planck 353 GHz channel. For those lines of sight through the diffuse interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, $R_{S/V} = (P_S/I_S)/(p_V/\tau_V)$ and $R_{P/p} = P_S/p_V$ , the latter focusing directly on the polarization properties of the aligned grain population alone. We find $R_{S/V}$ = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and $R_{P/p}$ = 5.4 MJy sr$^{-1}$, with uncertainties 0.2 and 0.3 MJy sr$^{-1}$, respectively. Our estimate of $R_{S/V}$ is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse interstellar medium. However, our estimate of $R_{P/p}$ is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, $R_{P/p}$, indicates that changes to the optical properties in the models of the aligned grain population are required.

Planck has mapped the microwave sky in nine frequency bands between 30 and 857 GHz in temperature and seven bands between 30 and 353 GHz in polarization. In this paper we consider the problem of diffuse astrophysical component separation, and process these maps within a Bayesian framework to derive a consistent set of full-sky astrophysical component maps. For the temperature analysis, we combine the Planck observations with the 9-year WMAP sky maps and the Haslam et al. 408 MHz map to derive a joint model of CMB, synchrotron, free-free, spinning dust, CO, line emission in the 94 and 100 GHz channels, and thermal dust emission. Full-sky maps are provided with angular resolutions varying between 7.5 arcmin and 1 deg. Global parameters (monopoles, dipoles, relative calibration, and bandpass errors) are fitted jointly with the sky model, and best-fit values are tabulated. For polarization, the model includes CMB, synchrotron, and thermal dust emission. These models provide excellent fits to the observed data, with rms temperature residuals smaller than 4 uK over 93% of the sky for all Planck frequencies up to 353 GHz, and fractional errors smaller than 1% in the remaining 7% of the sky. The main limitations of the temperature model at the lower frequencies are degeneracies among the spinning dust, free-free, and synchrotron components; additional observations from external low-frequency experiments will be essential to break these. The main limitations of the temperature model at the higher frequencies are uncertainties in the 545 and 857 GHz calibration and zero-points. For polarization, the main outstanding issues are instrumental systematics in the 100-353 GHz bands on large angular scales in the form of temperature-to-polarization leakage, uncertainties in the analog-to-digital conversion, and very long time constant corrections, all of which are expected to improve in the near future.

We present foreground-reduced CMB maps derived from the full Planck data set in both temperature and polarization. Compared to the corresponding Planck 2013 temperature sky maps, the total data volume is larger by a factor of 3.2 for frequencies between 30 and 70 GHz, and by 1.9 for frequencies between 100 and 857 GHz. In addition, systematic errors in the forms of temperature-to-polarization leakage, analogue-to-digital conversion uncertainties, and very long time constant errors have been dramatically reduced, to the extent that the cosmological polarization signal may now be robustly recovered on angular scales $\ell\gtrsim40$. On the very largest scales, instrumental systematic residuals are still non-negligible compared to the expected cosmological signal, and modes with $\ell < 20$ are accordingly suppressed in the current polarization maps by high-pass filtering. As in 2013, four different CMB component separation algorithms are applied to these observations, providing a measure of stability with respect to algorithmic and modelling choices. The resulting polarization maps have rms instrumental noise ranging between 0.21 and 0.27$\,\mu\textrm{K}$ averaged over 55 arcmin pixels, and between 4.5 and 6.1$\,\mu\textrm{K}$ averaged over 3.4 arcmin pixels. The cosmological parameters derived from the analysis of temperature power spectra are in agreement at the $1\sigma$ level with the Planck 2015 likelihood. Unresolved mismatches between the noise properties of the data and simulations prevent a satisfactory description of the higher-order statistical properties of the polarization maps. Thus, the primary applications of these polarization maps are those that do not require massive simulations for accurate estimation of uncertainties, for instance estimation of cross-spectra and cross-correlations, or stacking analyses.