(abridged) We present ALMA observations of hydrogen fluoride, HF J=1-0, H20 (220-211), and the 1.2 THz rest-frame continuum emission from the z=4.7 system BR1202-0725. BR1202-0725 is a galaxy group consisting of a QSO, a sub-millimeter galaxy (SMG), and two Ly-alpha emitters. We detected HF in emission in the QSO and possibly in absorption in the SMG, while water is detected in emission in both the QSO and SMG. The QSO is the most luminous HF emitter yet found and has the same ratio of HF emission line to infrared luminosity as a sample of local AGN and the Orion Bar. This consistency covers about 10 orders-of-magnitude in infrared luminosity, L_IR. Based on the conclusions of a study of HF emission in the Orion Bar and modeling, the HF emission in the QSO is either excited by collisions with electrons and H2 in molecular plasmas irradiated by the AGN and intense star formation or predominately by collisions with H2, with a modest contribution from electrons, in a relatively high temperature (~120 K), dense (~10^5 cm^-3) medium. Although HF should be an excellent tracer of molecular outflows, we do not find strong evidence for outflows in HF in either the QSO or the SMG. From a putative absorption feature in HF in the SMG, we estimate an upper limit on the outflow rate, dM/dt_outflow <~45 M_sun/yr. The ratio of the outflow rate to the star formation rate is <5% for the SMG. The broadness of the H2O line in the SMG, FWHM~1020 km/s, may suggest that either the gas on large scales (>4 kpc) is significantly more disturbed and turbulent due either to interactions and mass exchange with the other members of the group or to the dissipation of the energy of the intense star formation or both. The lack of significant molecular outflows in either source may imply that much of the energy from the intense star formation and AGN activity in this pair is being dissipated in their ISM.

With the spatial resolution of the Atacama Large Millimetre Array (ALMA), dusty galaxies in the distant Universe typically appear as single, compact blobs of dust emission, with a median half-light radius, $\approx$ 1 kpc. Occasionally, strong gravitational lensing by foreground galaxies or galaxy clusters has probed spatial scales 1-2 orders of magnitude smaller, often revealing late-stage mergers, sometimes with tantalising hints of sub-structure. One lensed galaxy in particular, the Cosmic Eyelash at $z=$ 2.3, has been cited extensively as an example of where the interstellar medium exhibits obvious, pronounced clumps, on a spatial scale of $\approx$ 100 pc. Seven orders of magnitude more luminous than giant molecular clouds in the local Universe, these features are presented as circumstantial evidence that the blue clumps observed in many $z\sim$ 2-3 galaxies are important sites of ongoing star formation, with significant masses of gas and stars. Here, we present data from ALMA which reveal that the dust continuum of the Cosmic Eyelash is in fact smooth and can be reproduced using two S\'ersic profiles with effective radii, 1.2 and 4.4 kpc, with no evidence of significant star-forming clumps down to a spatial scale of $\approx$ 80 pc and a star-formation rate of $<$ 3 M$_\odot$ yr$^{-1}$.

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.

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.

Understanding the physics of how stars form is a highly-prioritized goal of modern Astrophysics, in part because star formation is linked to both galactic dynamics on large scales and to the formation of planets on small scales. It is well-known that stars form from the gravitational collapse of molecular clouds, which are in turn formed out of the turbulent interstellar medium. Star formation is highly inefficient, with one of the likely culprits being the regulation against gravitational collapse provided by magnetic fields. Measurement of the polarized emission from interstellar dust grains, which are partially aligned with the magnetic field, provides a key tool for understanding the role these fields play in the star formation process. Over the past decade, much progress has been made by the most recent generation of polarimeters operating over a range of wavelengths (from the far-infrared through the millimeter part of the spectrum) and over a range of angular resolutions (from less than an arcsecond through fractions of a degree). Future developments in instrument sensitivity for ground-based, airborne, and space-borne polarimeters operating over range of spatial scales are critical for enabling revolutionary steps forward in our understanding of the magnetized turbulence from which stars are formed.

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.

We present Keck/KCWI integral field spectrograph observations of the complex system surrounding SMM J02399$-$0136 (a lensed $z=2.8$ sub-mm galaxy), including an associated Ly$\rm \alpha$ nebula, a dust-obscured, broad-absorption-line quasar, and neighboring galaxies. At a 3$\sigma$ surface brightness contour of 1.6$\times$ 10$^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, the Ly$\rm \alpha$ nebula extends over 17 arcsec ($\gtrsim 140$ physical kpc) and has a total Ly$\rm \alpha$ luminosity of $2.5 \times 10^{44} \, \rm erg \, s^{-1}$ (uncorrected for lensing). The nebula exhibits a kinematic shear of $\sim$ 1000 km s$^{-1}$ over 100 pkpc with lowest velocities east of SMM J02399$-$0136 and increasing to the southwest. We also discover a bright, Ly$\rm \alpha$ emitter, separated spatially and kinematically from the nebula, at a projected separation of $\approx$60 kpc from the quasar. This source has no clear central counterpart in deep $Hubble \, Space \, Telescope$ imaging, giving an intrinsic Ly$\rm \alpha$ rest-frame equivalent width greater than 312 $ \rm \mathring{A}$ (5-$\sigma$). We argue that this 'dark cloud' is illuminated by the quasar with a UV flux that is orders of magnitude brighter than the emission along our sightline. This result confirms statistical inferences that luminous quasars at $z>2$ emit UV radiation anisotropically. Future KCWI observations of other lines, e.g. Ly$\rm \beta$, He$\rm _{II}$, C$\rm _{IV}$, etc, and with polarimetry will further reveal the origin of the Ly$\rm \alpha$ nebula and nature of the dark cloud.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

Aims. The goal of the paper is to present a detailed study of the propagation of low velocity (5 to 25 km s-1) stationary molecular shocks in environments illuminated by an external ultraviolet (UV) radiation field. In particular, we intend to show how the structure, dynamics, energetics, and chemical properties of shocks are modified by UV photons and to estimate how efficiently shocks can produce line emission. Methods. We implemented several key physico-chemical processes in the Paris-Durham shock code to improve the treatment of the radiative transfer and its impact on dust and gas particles. We propose a new integration algorithm to find the steady-state solutions of magnetohydrodynamics equations in a range of parameters in which the fluid evolves from a supersonic to a subsonic regime. We explored the resulting code over a wide range of physical conditions, which encompass diffuse interstellar clouds and hot and dense photon-dominated regions (PDR). Results. We find that C-type shock conditions cease to exist as soon as G0 > 0.2 (nH/cm-3)^1/2. Such conditions trigger the emergence of another category of stationary solutions, called C*-type and CJ-type shocks, in which the shocked gas is momentarily subsonic along its trajectory. These solutions are shown to be unique for a given set of physical conditions and correspond to dissipative structures in which the gas is heated up to temperatures comprised between those found in C-type and adiabatic J-type shocks. High temperatures combined with the ambient UV field favour the production or excitation of a few molecular species to the detriment of others, hence leading to specific spectroscopic tracers such as rovibrational lines of H2 and rotational lines of CH+. Unexpectedly, the rotational lines of CH+ may carry as much as several percent of the shock kinetic energy.

We utilize the multi-wavelength data of M33 to study the origin of turbulence in its interstellar medium. We find that the HI turbulent energy surface density inside 8 kpc is $\sim1-3~\times~10^{46}$ erg pc$^{-2}$, and has no strong dependence on galactocentric radius because of the lack of variation in HI surface density and HI velocity dispersion. Then, we consider the energies injected by supernovae (SNe), the magneto-rotational instability (MRI), and the gravity-driven turbulence from accreted materials as the sources of turbulent energy. For a constant dissipation time of turbulence, the SNe energy can maintain turbulence inside $\sim 4$ kpc radius (equivalent to $\sim 0.5~R_{25}$), while the MRI energy is always smaller than the turbulent energy within 8 kpc radius. However, when we let the dissipation time to be equal to the crossing time of turbulence across the HI scale-height, the SNe energy is enough to maintain turbulence out to 7 kpc radius, and the sum of SNe and MRI energies is able to maintain turbulence out to 8 kpc radius. Due to lack of constraint in the mass accretion rate through the disk of M33, we can not rule out the accretion driven turbulence as a possible source of energy. Furthermore, by resolving individual Giant Molecular Clouds in M33, we also show that the SNe energy can maintain turbulence within individual molecular clouds with $\sim 1\%$ of coupling efficiency. This result strengthens the proposition that stellar feedback is an important source of energy to maintain turbulence in nearby galaxies.

The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization. We use new Planck maps to characterize Galactic dust emission as a foreground to the CMB polarization. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for six nested sky regions covering from 24 to 71 % of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra. The TE correlation and E/B power asymmetry extend to low multipoles that were not included in earlier Planck polarization papers. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we determine the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component separation procedure required for detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified blackbody emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean spectral index for dust polarization is $\beta_{\rm d}^{P} = 1.53\pm0.02 $. By fitting multi-frequency cross-spectra, we examine the correlation of the dust polarization maps across frequency. We find no evidence for decorrelation. If the Planck limit for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then decorrelation might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio $r\simeq0.01$ at the recombination peak.

Using the ${\it Planck}$ full-mission data, we present a detection of the temperature (and therefore velocity) dispersion due to the kinetic Sunyaev-Zeldovich (kSZ) effect from clusters of galaxies. To suppress the primary CMB and instrumental noise we derive a matched filter and then convolve it with the ${\it Planck}$ foreground-cleaned `${\tt 2D-ILC\,}$' maps. By using the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we determine the normalized ${\it rms}$ dispersion of the temperature fluctuations at the positions of clusters, finding that this shows excess variance compared with the noise expectation. We then build an unbiased statistical estimator of the signal, determining that the normalized mean temperature dispersion of $1526$ clusters is $\langle \left(\Delta T/T \right)^{2} \rangle = (1.64 \pm 0.48) \times 10^{-11}$. However, comparison with analytic calculations and simulations suggest that around $0.7\,\sigma$ of this result is due to cluster lensing rather than the kSZ effect. By correcting this, the temperature dispersion is measured to be $\langle \left(\Delta T/T \right)^{2} \rangle = (1.35 \pm 0.48) \times 10^{-11}$, which gives a detection at the $2.8\,\sigma$ level. We further convert uniform-weight temperature dispersion into a measurement of the line-of-sight velocity dispersion, by using estimates of the optical depth of each cluster (which introduces additional uncertainty into the estimate). We find that the velocity dispersion is $\langle v^{2} \rangle =(123\,000 \pm 71\,000)\,({\rm km}\,{\rm s}^{-1})^{2}$, which is consistent with findings from other large-scale structure studies, and provides direct evidence of statistical homogeneity on scales of $600\,h^{-1}{\rm Mpc}$. Our study shows the promise of using cross-correlations of the kSZ effect with large-scale structure in order to constrain the growth of structure.

Cosmic rays are predominantly accelerated in shocks associated with star formation such as supernova remnants and stellar wind bubbles, so the cosmic-ray flux and thus cosmic-ray ionization rate, $\zeta_{\rm H}$, should correlate with the star-formation rate in a galaxy. Submillimeter bright galaxies (SMGs) are some of the most prolific star forming galaxies in the Universe, and gravitationally lensed SMGs provide bright continuum sources suitable for absorption line studies. Abundances of OH$^+$ and H$_2$O$^+$ are useful for inferring $\zeta_{\rm H}$ when combined with chemical models, and have been used for this purpose within the Milky Way. At redshifts $z\gtrsim2$ transitions out of the ground rotational states of OH$^+$ and H$_2$O$^+$ are observable with ALMA, and we present observations of both molecules in absorption toward the lensed SMGs SMM J2135$-$0102 and SDP 17b. These detections enable an exploration of $\zeta_{\rm H}$ in galaxies with extreme star formation and high supernova rates, both of which should significantly enhance cosmic-ray production. The observed OH$^+$ and H$_2$O$^+$ absorption is thought to arise in massive, extended halos of cool, diffuse gas that surround these galaxies. Using a chemical model designed to focus on the reaction network important to both species, we infer cosmic-ray ionization rates of $\zeta_{\rm H}\sim10^{-16}$-$10^{-14}$ s$^{-1}$ in these extended gaseous halos. As our estimates come from gas that is far away from the sites of cosmic-ray acceleration, they imply that cosmic-ray ionization rates in the compact regions where star formation occurs in these galaxies are orders of magnitude higher.

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.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

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.

[abridged] The interstellar medium is now widely recognized to display features ascribable to magnetized turbulence. With the public release of Planck data and the current balloon-borne and ground-based experiments, the growing amount of data tracing the polarized thermal emission from Galactic dust in the submillimetre provides choice diagnostics to constrain the properties of this magnetized turbulence. We aim to constrain these properties in a statistical way, focusing in particular on the power spectral index of the turbulent component of the interstellar magnetic field in a diffuse molecular cloud, the Polaris Flare. We present an analysis framework which is based on simulating polarized thermal dust emission maps using model dust density (proportional to gas density) and magnetic field cubes, integrated along the line of sight, and comparing these statistically to actual data. The model fields are derived from fBm processes, which allow a precise control of their one- and two-point statistics. We explore the nine-dimensional parameter space of these models through a MCMC analysis, which yields best-fitting parameters and associated uncertainties. We find that the power spectrum of the turbulent component of the magnetic field in the Polaris Flare molecular cloud scales with wavenumber as a power law with a spectral index $2.8\pm 0.2$. It complements a uniform field whose norm in the POS is approximately twice the norm of the fluctuations of the turbulent component. The density field is well represented by a log-normally distributed field with a mean gas density $40\,\mathrm{cm}^{-3}$ and a power spectrum with as spectral index $1.7^{+0.4}_{-0.3}$. The agreement between the Planck data and the simulated maps for these best-fitting parameters is quantified by a $\chi^2$ value that is only slightly larger than unity.

As the next step toward an improved large scale Galactic magnetic field model, we present a simple comparison of polarised synchrotron and thermal dust emission on the Galactic plane. We find that the field configuration in our previous model that reproduces the polarised synchrotron is not compatible with the WMAP 94 GHz polarised emission data. In particular, the high degree of dust polarisation in the outer Galaxy (90deg < l < 270deg) implies that the fields in the dust-emitting regions are more ordered than the average of synchrotron-emitting regions. This new dust information allows us to constrain the spatial mixing of the coherent and random magnetic field components in the outer Galaxy. The inner Galaxy differs in polarisation degree and apparently requires a more complicated scenario than our current model. In the scenario that each interstellar component (including fields and now dust) follows a spiral arm modulation, as observed in external galaxies, the changing degree of ordering of the fields in dust-emitting regions may imply that the dust arms and the field component arms are shifted as a varying function of Galacto-centric radius. We discuss the implications for how the spiral arm compression affects the various components of the magnetised interstellar medium but conclude that improved data such as that expected from the Planck satellite will be required for a thorough analysis.

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.

Using the PACS and SPIRE spectrometers on-board Herschel, we obtained observations of the Taffy galaxies (UGC 12914/12915) and bridge. The Taffy system is believed to be the result of a face-on collision between two gas-rich galaxies, in which the stellar disks passed through each other, but the gas was dispersed into a massive H I and molecular bridge between them. Emission is detected and mapped in both galaxies and the bridge in the [C II]157.7 $\mu$m and [O I]63.2 $\mu$m fine-structure lines. Additionally, SPIRE FTS spectroscopy detects the [C I] $^3$P$_2$$\rightarrow$$^3$P$_1$(809.3 GHz) and [C I] $^3$P$_1$$\rightarrow$$3$P$_0$(492.2 GHz) neutral carbon lines, and weakly detects high-J CO transitions in the bridge. These results indicate that the bridge is composed of a warm multi-phase medium consistent with shock and turbulent heating. Despite low star formation rates in the bridge, the [C II] emission appears to be enhanced, reaching [C II]/FIR ratios of 3.3% in parts of the bridge. Both the [C II] and [O I] lines show broad intrinsic multi-component profiles, similar to those seen in previous CO 1-0 and H I observations. The [C II] emission shares similar line profiles with both the double-peaked H I profiles and shares a high-velocity component with single-peaked CO profiles in the bridge, suggesting that the [C II] emission originates in both the neutral and molecular phases. We show that it is feasible that a combination of turbulently heated H$_2$ and high column-density H I, resulting from the galaxy collision, is responsible for the enhanced [C II] emission.

Starburst galaxies at the peak of cosmic star formation are among the most extreme starforming engines in the universe, producing stars over ~100 Myr. The star formation rates of these galaxies, which exceed 100 $M_\odot$ per year, require large reservoirs of cold molecular gas to be delivered to their cores, despite strong feedback from stars or active galactic nuclei. Starburst galaxies are therefore ideal targets to unravel the critical interplay between this feedback and the growth of a galaxy. The methylidyne cation, CH$^+$, is a most useful molecule for such studies because it cannot form in cold gas without supra-thermal energy input, so its presence highlights dissipation of mechanical energy or strong UV irradiation. Here, we report the detection of CH$^+$(J=1-0) emission and absorption lines in the spectra of six lensed starburst galaxies at redshifts z~2.5. This line has such a high critical density for excitation that it is emitted only in very dense ($>10^5$ cm$^{-3}$) gas, and is absorbed in low-density gas. We find that the CH$^+$ emission lines, which are broader than 1000 km s$^{-1}$, originate in dense shock waves powered by hot galactic winds. The CH$^+$ absorption lines reveal highly turbulent reservoirs of cool ($T\sim 100$K), low-density gas, extending far outside (>10 kpc) the starburst cores (radii <1 kpc). We show that the galactic winds sustain turbulence in the 10 kpc-scale environments of the starburst cores, processing these environments into multi-phase, gravitationally bound reservoirs. However, the mass outflow rates are found to be insufficient to balance the star formation rates. Another mass input is therefore required for these reservoirs, which could be provided by on-going mergers or cold stream accretion. Our results suggest that galactic feedback, coupled jointly to turbulence and gravity, extends the starburst phase instead of quenching it.

The z=0.89 molecular absorber toward PKS1830-211 provides us with the opportunity to probe the chemical and physical properties of the interstellar medium in the disk of a galaxy at a look-back time of half the present age of the Universe. Recent ALMA observations of hydrides have unveiled the multi-phase composition of this source's interstellar medium along two absorbing sightlines. Here, we report ALMA observations of CH+ and SH+, and of their 13C- and 34S- isotopologues, as potential tracers of energetic processes in the interstellar medium. CH+ and 13CH+ are detected toward both images of PKS1830-211, CH+ showing the deepest and broadest absorption among all species observed so far. The [CH+]/[13CH+] abundance ratio is ~100 in the south-west line of sight. [...] Toward the north-east image, we find an even larger value of [CH+]/[13CH+], 146 +/- 43, although with a large uncertainty. This sightline intercepts the absorber at a larger galactocentric radius than the southwestern one, where material might be less processed in stellar nucleosynthesis. In contrast to CH+ and its 13C isotopologue, SH+ and 34SH+ are only detected on the south-west sightline. These are the first detections of extragalactic SH+ and interstellar 34SH+. The spectroscopic parameters of SH+ are reevaluated and improved rest frequencies of 34SH+ are obtained. The [CH+]/[SH+] column density ratios show a large difference between the two lines of sight: ~25 and >600 toward the SW and NE image, respectively. We are not able to shed light on the formation process of CH+ and SH+ with these data, but the differences in the two sightlines toward PKS1830-211 suggest that their absorptions arise from gas with molecular fraction gtrsim 10%, with SH+ tracing significantly higher molecular fractions than CH+.

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.

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 report the detection of OH+ and H2O+ in the z=0.89 absorber toward the lensed quasar PKS1830-211. The abundance ratio of OH+ and H2O+ is used to quantify the molecular hydrogen fraction (fH2) and the cosmic-ray ionization rate of atomic hydrogen (zH) along two lines of sight, located at ~2 kpc and ~4 kpc to either side of the absorber's center. The molecular fraction decreases outwards, from ~0.04 to ~0.02, comparable to values measured in the Milky Way at similar galactocentric radii. For zH, we find values of ~2x10^-14 s^-1 and ~3x10^-15 s^-1, respectively, which are slightly higher than in the Milky Way at comparable galactocentric radii, possibly due to a higher average star formation activity in the z=0.89 absorber. The ALMA observations of OH+, H2O+, and other hydrides toward PKS1830-211 reveal the multi-phase composition of the absorbing gas. Taking the column density ratios along the southwest and northeast lines of sight as a proxy of molecular fraction, we classify the species ArH+, OH+, H2Cl+, H2O+, CH, and HF as tracing gases increasingly more molecular. Incidentally, our data allow us to improve the accuracy of H2O+ rest frequencies and thus refine the spectroscopic parameters.

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.

Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular power spectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above $b = \pm 20{\deg}$. We find that the dust temperature is $T = (19.4 \pm 1.3)$ K and the dust spectral index is $\beta = 1.6 \pm 0.1$ averaged over the whole sky, while $T = (19.4 \pm 1.5)$ K and $\beta = 1.6 \pm 0.2$ on 21 % of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60 % of the sky at Galactic latitudes $|b| > 20{\deg}$. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.

The origin of striations aligned along the local magnetic field direction in the translucent envelope of the Taurus molecular cloud is examined with new observations of 12CO and 13CO J=2-1 emission obtained with the 10~m submillimeter telescope of the Arizona Radio Observatory. These data identify a periodic pattern of excess blue and redshifted emission that is responsible for the striations. For both 12CO and 13CO, spatial variations of the J=2-1 to J=1-0 line ratio are small and are not spatially correlated with the striation locations. A medium comprised of unresolved CO emitting substructures (cells) with a beam area filling factor less than unity at any velocity is required to explain the average line ratios and brightness temperatures. We propose that the striations result from the modulation of velocities and the beam filling factor of the cells as a result of either the Kelvin-Helmholtz instability or magnetosonic waves propagating through the envelope of the Taurus molecular cloud. Both processes are likely common features in molecular clouds that are sub-Alfvenic and may explain low column density, cirrus-like features similarly aligned with the magnetic field observed throughout the interstellar medium in far-infrared surveys of dust emission.

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.

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

Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field orientation directly in vicinities with an area equivalent to that subtended by the Planck effective beam at 353 GHz (10') and using the second-order structure functions of the field orientation angles. We find that the average dispersion of the starlight-inferred field orientations within 10'-diameter vicinities is less than 20 deg, and that at these scales the mean field orientation is on average within 5 deg of that inferred from the submillimetre polarization observations in the considered regions. We also find that the dispersion of starlight polarization orientations and the polarization fractions within these vicinities are well reproduced by a Gaussian model of the turbulent structure of the magnetic field, in agreement with the findings reported by the Planck collaboration at scales greater than 10' and for comparable column densities. At scales greater than 10', we find differences of up to 14.7 deg between the second-order structure functions obtained from starlight and submillimetre polarization observations in the same positions in the plane of the sky, but comparison with a Gaussian model of the turbulent structure of the magnetic field indicates that these differences are small and are consistent with the difference in angular resolution between both techniques.

This paper describes the identification, modelling, and removal of previously unexplained systematic effects in the polarization data of the Planck High Frequency Instrument (HFI) on large angular scales, including new mapmaking and calibration procedures, new and more complete end-to-end simulations, and a set of robust internal consistency checks on the resulting maps. These maps, at 100, 143, 217, and 353 GHz, are early versions of those that will be released in final form later in 2016. The improvements allow us to determine the cosmic reionization optical depth $\tau$ using, for the first time, the low-multipole $EE$ data from HFI, reducing significantly the central value and uncertainty, and hence the upper limit. Two different likelihood procedures are used to constrain $\tau$ from two estimators of the CMB $E$- and $B$-mode angular power spectra at 100 and 143 GHz, after debiasing the spectra from a small remaining systematic contamination. These all give fully consistent results. A further consistency test is performed using cross-correlations derived from the Low Frequency Instrument maps of the Planck 2015 data release and the new HFI data. For this purpose, end-to-end analyses of systematic effects from the two instruments are used to demonstrate the near independence of their dominant systematic error residuals. The tightest result comes from the HFI-based $\tau$ posterior distribution using the maximum likelihood power spectrum estimator from $EE$ data only, giving a value $0.055\pm 0.009$. In a companion paper these results are discussed in the context of the best-fit Planck $\Lambda$CDM cosmological model and recent models of reionization.

The CH$^+$ ion is a key species in the initial steps of interstellar carbon chemistry. Its formation in diverse environments where it is observed is not well understood, however, because the main production pathway is so endothermic (4280 K) that it is unlikely to proceed at the typical temperatures of molecular clouds. We investigation CH$^+$ formation with the first velocity-resolved spectral mapping of the CH$^+$ $J=1-0, 2-1$ rotational transitions, three sets of CH $\Lambda$-doubled triplet lines, $^{12}$C$^+$ and $^{13}$C$^+$, and CH$_3$OH 835~GHz E-symmetry Q branch transitions, obtained with Herschel/HIFI over $\approx$12 arcmin$^2$ centered on the Orion BN/KL source. We present the spatial morphologies and kinematics, cloud boundary conditions, excitation temperatures, column densities, and $^{12}$C$^+$ optical depths. Emission from C$^+$, CH$^+$, and CH is indicated to arise in the diluted gas, outside of the explosive, dense BN/KL outflow. Our models show that UV-irradiation provides favorable conditions for steady-state production of CH$^+$ in this environment. Surprisingly, no spatial or kinematic correspondences of these species are found with H$_2$ S(1) emission tracing shocked gas in the outflow. We propose that C$^+$ is being consumed by rapid production of CO to explain the lack of C$^+$ and CH$^+$ in the outflow, and that fluorescence provides the reservoir of H$_2$ excited to higher ro-vibrational and rotational levels. Hence, in star-forming environments containing sources of shocks and strong UV radiation, a description of CH$^+$ formation and excitation conditions is incomplete without including the important --- possibly dominant --- role of UV irradiation.

We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($\psi$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $\psi$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $\psi$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS.

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 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 present the first 7.5'x11.5' velocity-resolved map of the [CII]158um line toward the Orion molecular cloud-1 (OMC-1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41alpha hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [CII] luminosity (~85%) is from the extended, FUV-illuminated face of the cloud G_0>500, n_H>5x10^3 cm^-3) and from dense PDRs (G_0~10^4, n_H~10^5 cm^-3) at the interface between OMC-1 and the HII region surrounding the Trapezium cluster. Around 15% of the [CII] emission arises from a different gas component without CO counterpart. The [CII] excitation, PDR gas turbulence, line opacity (from [13CII]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the [CII]/FIR and FIR/M_Gas ratios and show that [CII]/FIR decreases from the extended cloud component (10^-2-10^-3) to the more opaque star-forming cores (10^-3-10^-4). The lowest values are reminiscent of the "[CII] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing [CII]/FIR ratio correlates better with the column density of dust through the molecular cloud than with FIR/M_Gas. We conclude that the [CII] emitting column relative to the total dust column along each line of sight is responsible for the observed [CII]/FIR variations through the cloud.

Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impact of an expanding HII region on the morphology of the field. We derive an analytical solution for the magnetic field, assumed to evolve from an initially uniform configuration following the expansion of ionized gas and the formation of a shell of swept-up ISM. From the RM data we estimate a mean value of the line-of-sight component of the magnetic field of about 3microG (towards the observer) in the Rosette Nebula, for a uniform electron density of about 12cm-3. The dust shell that surrounds the Rosette HII region is clearly observed in the Planck intensity map at 353 GHz, with a polarization signal significantly different from that of the local background when considered as a whole. The Planck observations constrain the plane-of-the-sky orientation of the magnetic field in the Rosette's parent molecular cloud to be mostly aligned with the large-scale field along the Galactic plane. The Planck data are compared with the analytical model, which predicts the mean polarization properties of a spherical and uniform dust shell for a given orientation of the field. This comparison leads to an upper limit of about 45degr on the angle between the line of sight and the magnetic field in the Rosette complex, for an assumed intrinsic dust polarization fraction of 4%. This field direction can reproduce the RM values detected in the ionized region if the magnetic field strength in the Monoceros molecular cloud is in the range 6.5--9microG.

Within ten nearby (d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, $N_H$. The selected regions, covering several degrees in size, are analyzed at an effective angular resolution of 10' FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed in bins of column density using the novel statistical tool called "Histogram of Relative Orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing $N_H$, from preferentially parallel or having no preferred orientation to preferentially perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfv\'enic or sub-Alfv\'enic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

The Planck mission, thanks to its large frequency range and all-sky coverage, has a unique potential for systematically detecting the brightest, and rarest, submillimetre sources on the sky, including distant objects in the high-redshift Universe traced by their dust emission. A novel method, based on a component-separation procedure using a combination of Planck and IRAS data, has been applied to select the most luminous cold submm sources with spectral energy distributions peaking between 353 and 857GHz at 5' resolution. A total of 2151 Planck high-z source candidates (the PHZ) have been detected in the cleanest 26% of the sky, with flux density at 545GHz above 500mJy. Embedded in the cosmic infrared background close to the confusion limit, these high-z candidates exhibit colder colours than their surroundings, consistent with redshifts z>2, assuming a dust temperature of 35K and a spectral index of 1.5. First follow-up observations obtained from optical to submm have confirmed that this list consists of two distinct populations. A small fraction (around 3%) of the sources have been identified as strongly gravitationally lensed star-forming galaxies, which are amongst the brightest submm lensed objects (with flux density at 545GHz ranging from 350mJy up to 1Jy) at redshift 2 to 4. However, the vast majority of the PHZ sources appear as overdensities of dusty star-forming galaxies, having colours consistent with z>2, and may be considered as proto-cluster candidates. The PHZ provides an original sample, complementary to the Planck Sunyaev-Zeldovich Catalogue; by extending the population of the virialized massive galaxy clusters to a population of sources at z>1.5, the PHZ may contain the progenitors of today's clusters. Hence the PHZ opens a new window on the study of the early ages of structure formation, and the understanding of the intensively star-forming phase at high-z.

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.

Planck observations at 353GHz provide the first fully-sampled maps of the polarized dust emission towards interstellar filaments and their backgrounds. The polarization data provide insight on the structure of their magnetic field (B). We present the polarization maps of three nearby star forming filament of moderate column density (NH~10^22cm^-2): Musca, B211, and L1506. We use the spatial information to separate Stokes I, Q, and U of the filaments from those of their backgrounds, an essential step to measure the intrinsic polarization fraction (p) and angle (psi) of each emission component. We find that the polarization angles in the three filaments (psi_fil) are coherent along their lengths and not the same as in their backgrounds (psi_bg). The differences between psi_fil and psi_bg are 12deg, 6deg, and 54deg for Musca, B211, and L1506, respectively. These differences for Musca and L1506 are larger than the dispersions of psi, both along the filaments and in their backgrounds. The observed changes of psi are direct evidence for variations of the orientation of the plane of the sky (POS) projection of the B-field. As in previous studies, we find a decrease of several percent of p with NH. We show that the drop in p cannot be explained by random fluctuations of the orientation of B within the filaments because they are too small (sigma_psi<10deg). We recognize the degeneracy between dust alignment efficiency and the structure of B in causing variations in p, but we argue that the decrease of p from the backgrounds to the filaments results in part from depolarization associated with the 3D structure of B: both its orientation in the POS and with respect to the POS. We do not resolve the inner structure of the filaments, but at the smallest scales accessible with Planck (~0.2pc), the observed changes of psi and p hold information on the B-field structure within filaments.

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 Atacama Large Millimeter Array (ALMA) observations of two high-redshift systems (SMMJ02399-0136 and the Cloverleaf QSO) in their rest-frame 122 micron continuum (~650 GHz or ~450 micron on-sky) and [NII] 122 micron line emission. The continuum observations with a synthesized beam of ~0.25" resolve both sources and recover the expected flux. The Cloverleaf is resolved into a partial Einstein ring, while the SMMJ02399-0136 is unambiguously separated into two components; an AGN associated point source and an extend region at the location of a previously identified dusty starburst. We detect the [NII] line in both systems, though significantly weaker than our previous detections made with the 1st generation z(Redshift) and Early Universe Spectrometer. We show that this discrepancy is mostly explained if the line flux is resolved out due to significantly more extended emission and longer ALMA baselines than expected. Based on the ALMA observations we determine that greater than 75% of the total [NII] line flux in each source is produced via star formation. We use the [NII] line flux that is recovered by ALMA to constrain the N/H abundance, ionized gas mass, hydrogen ionizing photon rate, and star formation rate. In SMMJ02399-0136 we discover it contains a significant amount (~1000 solar masses per year) of unobscured star formation in addition to its dusty starburst and argue that SMMJ02399-0136 may be similar to the Antennae Galaxies (Arp 244) locally. In total these observations provide a new look at two well-studied systems while demonstrating the power and challenges of Band-9 ALMA observations of high-z systems.

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.

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.

We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2\Pi_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H$_2$S 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N(SH)/N(H$_2$S) in the range 1.1 - 3.0, N(CS)/N(H$_2$S) in the range 0.32 - 0.61, and N(SO)/N(H$_2$S) in the range 0.08 - 0.30. The observed SH/H$_2$ ratios - in the range (0.5-2.6) $\times 10^{-8}$ - indicate that SH (and other sulphur-bearing molecules) account for << 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H$_2$S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity drift, relative to H$_2$, of SH molecules produced by the dissociative recombination of H$_3$S$^+$ - reduce this discrepancy to a factor ~ 3.

The dissipation of kinetic and magnetic energy in the interstellar medium (ISM) can proceed through viscous, Ohmic or ambipolar diffusion (AD). It occurs at very small scales compared to the scales at which energy is presumed to be injected. This localized heating may impact the ISM evolution but also its chemistry, thus providing observable features. Here, we perform 3D spectral simulations of decaying magnetohydrodynamic turbulence including the effects of AD. We find that the AD heating power spectrum peaks at scales in the inertial range, due to a strong alignment of the magnetic and current vectors in the dissipative range. AD affects much greater scales than the AD scale predicted by dimensional analysis. We find that energy dissipation is highly concentrated on thin sheets. Its probability density function follows a lognormal law with a power-law tail which hints at intermittency, a property which we quantify by use of structure function exponents. Finally, we extract structures of high dissipation, defined as connected sets of points where the total dissipation is most intense and we measure the scaling exponents of their geometric and dynamical characteristics: the inclusion of AD favours small sizes in the dissipative range.

We present the Planck Catalogue of Galactic Cold Clumps (PGCC), an all-sky catalogue of Galactic cold clump candidates detected by Planck. This catalogue is the full version of the Early Cold Core (ECC) catalogue, which was made available in 2011 with the Early Release Compact Source Catalogue (ERCSC) and contained 915 high S/N sources. It is based on the Planck 48 months mission data that are currently being released to the astronomical community. The PGCC catalogue is an observational catalogue consisting exclusively of Galactic cold sources. The three highest Planck bands (857, 545, 353 GHz) have been combined with IRAS data at 3 THz to perform a multi-frequency detection of sources colder than their local environment. After rejection of possible extragalactic contaminants, the PGCC catalogue contains 13188 Galactic sources spread across the whole sky, i.e., from the Galactic plane to high latitudes, following the spatial distribution of the main molecular cloud complexes. The median temperature of PGCC sources lies between 13 and 14.5 K, depending on the quality of the flux density measurements, with a temperature ranging from 5.8 to 20 K after removing sources with the 1% largest temperature estimates. Using seven independent methods, reliable distance estimates have been obtained for 5574 sources, which allows us to derive their physical properties such as their mass, physical size, mean density and luminosity. The PGCC sources are located mainly in the solar neighbourhood, up to a distance of 10.5 kpc towards the Galactic centre, and range from low-mass cores to large molecular clouds. Because of this diversity and because the PGCC catalogue contains sources in very different environments, the catalogue is useful to investigate the evolution from molecular clouds to cores. Finally, the catalogue also includes 54 additional sources located in the SMC and LMC.