This paper presents cosmological results based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. Our results are in very good agreement with the 2013 analysis of the Planck nominal-mission temperature data, but with increased precision. The temperature and polarization power spectra are consistent with the standard spatially-flat 6-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations (denoted "base ΛCDM" in this paper). From the Planck temperature data combined with Planck lensing, for this cosmology we find a Hubble constant, H 0 = (67.8 ± 0.9) km s −1 Mpc −1 , a matter density parameter Ω m = 0.308 ± 0.012, and a tilted scalar spectral index with n s = 0.968 ± 0.006, consistent with the 2013 analysis. Note that in this abstract we quote 68% confidence limits on measured parameters and 95% upper limits on other parameters. We present the first results of polarization measurements with the Low Frequency Instrument at large angular scales. Combined with the Planck temperature and lensing data, these measurements give a reionization optical depth of τ = 0.066 ± 0.016, corresponding to a reionization redshift of z re = 8.8+1.7 −1.4 . These results are consistent with those from WMAP polarization measurements cleaned for dust emission using 353-GHz polarization maps from the High Frequency Instrument. We find no evidence for any departure from base ΛCDM in the neutrino sector of the theory; for example, combining Planck observations with other astrophysical data we find N eff = 3.15 ± 0.23 for the effective number of relativistic degrees of freedom, consistent with the value N eff = 3.046 of the Standard Model of particle physics. The sum of neutrino masses is constrained to m ν < 0.23 eV. The spatial curvature of our Universe is found to be very close to zero, with |Ω K | < 0.005. Adding a tensor component as a single-parameter extension to base ΛCDM we find an upper limit on the tensor-to-scalar ratio of r 0.002 < 0.11, consistent with the Planck 2013 results and consistent with the B-mode polarization constraints from a joint analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP B-mode data to our analysis leads to a tighter constraint of r 0.002 < 0.09 and disfavours inflationary models with a V(φ) ∝ φ 2 potential. The addition of Planck polarization data leads to strong constraints on deviations from a purely adiabatic spectrum of fluctuations. We find no evidence for any contribution from isocurvature perturbations or from cosmic defects. Combining Planck data with other astrophysical data, including Type Ia supernovae, the equation of state of dark energy is constrained to w = −1.006 ± 0.045, consistent with the expected Corresponding author: G. Efstathiou, e-mail: gpe@ast.cam.ac.ukArticle published by EDP Sciences A13, page 1 of 63 A&A 594, A13 (2016) value for a cosmological constant. The standard big bang nucleosynthesis predictions for the helium and deuterium abundanc...

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ∼ 1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40 − 8 + 8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 M ⊙ . An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ∼ 40 Mpc ) less than 11 hours after the merger by the One-Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ∼10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ∼ 9 and ∼ 16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC 4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta.

New far-infrared and submillimeter data are used to solidify and to extend to long wavelengths the empirical calibration of the infrared spectral energy distribution (SED) of normal star-forming galaxies. As was found by Dale et al. (2001), a single parameter family, characterized by fν (60µm) fν (100µm) , is adequate to describe the range of normal galaxy spectral energy distributions observed by IRAS and ISO from 3 to 100 µm. However, predictions based on the first generation models at longer wavelengths (122 to 850 µm) are increasingly overluminous compared to the data for smaller fν (60µm) fν (100µm) , or alternatively, for weaker global interstellar radiation fields. After slightly modifying the far-infrared/submillimeter dust emissivity in those models as a function of the radiation field intensity to better match the long wavelength data, a suite of SEDs from 3 µm to 20 cm in wavelength is presented. Results from relevant applications are also discussed, including submillimeter-based photometric redshift indicators, the infrared energy budget and simple formulae for recovering the bolometric infrared luminosity, and dust mass estimates in galaxies. Regarding the latter, since galaxy infrared SEDs are not well-described by single blackbody curves, the usual methods of estimating dust masses can be grossly inadequate. The improved model presented herein is used to provide a more accurate relation between infrared luminosity and dust mass.

With the goal of investigating the degree to which the mid-infrared emission traces the star formation rate (SFR), we analyze Spitzer 8 µm and 24 µm