We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution.

Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper.

The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) suite of instruments operated throughout the over two years of the Rosetta mission operations in the vicinity of comet 67P/Churyumov-Gerasimenko. It measured gas densities and composition throughout the comet's atmosphere, or coma. Here we present two-years' worth of measurements of the relative densities of the four major volatile species in the coma of the comet, H2O. CO2, CO and O2, by one of the ROSINA sub-systems called the Double Focusing Mass Spectrometer (DFMS). The absolute total gas densities were provided by the Comet Pressure Sensor (COPS), another ROSINA sub-system. DFMS is a very high mass resolution and high sensitivity mass spectrometer able to resolve at a tiny fraction of an atomic mass unit. We have analyzed the combined DFMS and COPS measurements using an inversion scheme based on spherical harmonics that solves for the distribution of potential surface activity of each species as the comet rotates, changing solar illumination, over short intervals and as the comet changes distance from the sun and orientation of its spin axis over long time intervals. We also use the surface boundary conditions derived from the inversion scheme to simulate the whole coma with our fully kinetic Direct Simulation Monte Carlo model and calculate the production rates of the four major species throughout the mission. We compare the derived production rates with revised remote sensing observations by the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) as well as with published observations from the Microwave Instrument for the Rosetta Orbiter (MIRO). Finally we use the variation of the surface production of the major species to calculate the total mass loss over the mission and, for different estimates of the dust/gas ratio, calculate the variation of surface loss over the nucleus.

The physical properties of cometary nuclei observed today relate to their complex history and help to constrain their formation and evolution. In this article, we review some of the main physical properties of cometary nuclei and focus in particular on the thermal, mechanical, structural and dielectric properties, emphasizing the progress made during the Rosetta mission. Comets have a low density of 480 $\pm$ 220 kg m-3 and a low permittivity of 1.9 - 2.0, consistent with a high porosity of 70 - 80 %, are weak with a very low global tensile strength $<$100 Pa, and have a low bulk thermal inertia of 0 - 60 J K-1 m-2 s-1/2 that allowed them to preserve highly volatiles species (e.g. CO, CO2, CH4, N2) into their interior since their formation. As revealed by 67P/Churyumov-Gerasimenko, the above physical properties vary across the nucleus, spatially at its surface but also with depth. The broad picture is that the bulk of the nucleus consists of a weakly bonded, rather homogeneous material that preserved primordial properties under a thin shell of processed material, and possibly covered by a granular material; this cover might in places reach a thickness of several meters. The properties of the top layer (the first meter) are not representative of that of the bulk nucleus. More globally, strong nucleus heterogeneities at a scale of a few meters are ruled out on 67P small lobe.

The MIRO instrument's remote sensing capability is integral in constraining water density, temperature and velocity fields in the coma of 67P/Churyumov-Gersimenko. Our aim is to quantify the contribution to the water density from all facets inside and outside the field-of-view (FOV) of MIRO, in both, nadir and limb geometries. This information is crucial for understanding the MIRO derived coma production rates and their relation to the nucleus characteristics, and inherent spatial resolution of the data. This study relies on a detailed 3D nucleus shape model, illumination conditions and the pointing information of the viewing geometry. With these parameters, we can evaluate the relative contribution of water density originating from facets directly inside the MIRO beam as well outside of the beam as a function of distance along the MIRO line-of-sight. We also calculate the ratio of in-beam versus out-of-beam number density. We demonstrate that despite the rather small MIRO field-of-view there is only a small fraction of molecules that originate from facets within the MIRO beam. This is true for nadir, but also translated into the limb observing geometry. The MIRO instrument cannot discriminate active from non-active regions directly from observations. This study also suggests that the beam averaged solar incidence angle, local time and mean normal vectors are not necessary related to molecules within the MIRO beam. These results also illustrate why the 1D spherical Haser model can be applied with relative success for analyzing the MIRO data (and generally any Rosetta measurements). The future possibilities of constraining gas activity distribution on the surface should use 3D codes extracting information from the MIRO spectral line shapes which contain additional information. The presented results are in fact applicable to all relevant instruments onboard Rosetta.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

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

We have explored the capabilities of dust extinction and $\gamma$ rays to probe the properties of the interstellar medium in the nearby anti-centre region. We have jointly modelled the $\gamma$-ray intensity and the stellar reddening, E(B-V) as a combination of H$_{\rm I}$-bright, CO-bright, and ionised gas components. The complementary information from dust reddening and $\gamma$ rays is used to reveal the dark gas not seen, or poorly traced, by H$_{\rm I}$, free-free, and $^{12}$CO emissions. We compare the total gas column densities, $N_{\rm{H}}$, derived from the $\gamma$ rays and stellar reddening with those inferred from a similar analysis (Remy et al. 2017) of $\gamma$ rays and of the optical depth of the thermal dust emission, $\tau_{353}$, at 353 GHz. We can therefore compare environmental variations in specific dust reddening, E(B-V)/$N_{\rm H}$, and in dust emission opacity (dust optical depth per gas nucleon), $\tau_{353}/N_{\rm{H}}$. Over the whole anti-centre region, we find an average E(B-V)/$N_{\rm H}$ ratio of $(2.02\pm0.48)\times$ $10^{-22}$~mag~cm$^2$, with maximum local variations of about $\pm30\%$ at variance with the two to six fold coincident increase seen in emission opacity as the gas column density increases. In the diffuse medium, the small variations in specific reddening, E(B-V)/$N_{\rm H}$ implies a rather uniform dust-to-gas mass ratio in the diffuse parts of the anti-centre clouds. The small amplitude of the E(B-V)/$N_{\rm H}$ variations with increasing $N_{\rm{H}}$ column density confirms that the large opacity $\tau_{353}/N_{\rm{H}}$ rise seen toward dense CO clouds is primarily due to changes in dust emissivity. The environmental changes are qualitatively compatible with model predictions based on mantle accretion on the grains and the formation of grain aggregates.

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths.

H$_{\rm I}$ 21-cm and $^{12}$CO 2.6-mm line emissions respectively trace the atomic and molecular gas phases, but they miss most of the opaque H$_{\rm I}$ and diffuse H$_2$ present in the Dark Neutral Medium (DNM) at the transition between the H$_{\rm I}$-bright and CO-bright regions. Jointly probing H$_{\rm I}$, CO, and DNM gas, we aim to constrain the threshold of the H$_{\rm I}$-H$_2$ transition. We also aim to measure gas mass fractions in the different phases and to test their relation to cloud properties. We have used dust optical depth measurements at 353 GHz, $\gamma$-ray maps at GeV energies, and H$_{\rm I}$ and CO line data to trace the gas column densities and map the DNM in nearby clouds toward the Galactic anticentre and Chamaeleon regions. The onset of detectable CO intensities varies by only a factor of 4 from cloud to cloud, between 0.6$ \times 10^{21}$ cm$^{-2}$ and 2.5$ \times 10^{21}$ cm$^{-2}$ in total gas column density. We observe larger H$_2$ column densities than linearly inferred from the CO intensities at $A_{\rm{V}}> 3$ magnitudes because of the large CO optical thickness: the additional H$_2$ mass in this regime represents on average 20$\%$ of the CO-inferred molecular mass. In the DNM envelopes, we find that the fraction of diffuse CO-dark H$_2$ in the molecular column densities decreases with increasing $A_{\rm{V}}$ in a cloud. For a half molecular DNM, the fraction decreases from more than 80$\%$ at 0.4 magnitude to less than 20$\%$ beyond 2 magnitudes. In mass, the DNM fraction varies with the cloud properties. The mass present in the DNM envelopes appears to scale with the molecular mass seen in CO as $M_{\rm H}^{\rm DNM} = (62 \pm 7) {M_{\rm H_2}^{\rm CO}}^{\,0.51 \pm 0.02}$ across two decades in mass. The phase transitions in these clouds show both common trends and environmental differences.

The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first identified from the Planck data as a region of particularly cold dust emission and was mapped at 100-500 $\mu$m as a part of Herschel key program Galactic Cold Cores. The observations reveal several velocity components in the field, with strongest line emission concentrated to velocity range $\sim [3,5]$ km $\rm s^{-1}$. The column density of molecular hydrogen along the filament varies from 1.0 to 2.3 $\times 10^{22}$ $\rm cm^{2}$. We have examined six cold clumps from the central part of the filament. The clumps have masses in the range $10 - 20$ $M_{\odot}$ ($\sim70M_{\odot}$ in total) and are close to or above the virial mass. Furthermore, the main filament is heavily fragmented and most of the the substructures have a mass lower than or close to the virial mass, suggesting that the filament is dispersing as a whole. Position-velocity maps of $\rm ^{12}CO$ and $\rm ^{13}CO$ lines indicate that at least one of the striations is kinematically connected to two of the clumps, potentially indicating mass accretion from the striation onto the main filament. We tentatively estimate the accretion rate to be $\dot{M}$ = $2.23 \times 10^{-6}$ $ M_{\odot} / \rm year$.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

We aim to explore the capabilities of dust emission and rays for probing the properties of the interstellar medium in the nearby anti-centre region, using gamma-ray observations with the Fermi Large Area Telescope (LAT), and the thermal dust optical depth inferred from Planck and IRAS observations. In particular, we aim at quantifying potential variations in cosmic-ray density and dust properties per gas nucleon across the different gas phases and different clouds, and at measuring the CO-to-H2 conversion factor, X$_{CO}$ , in different environments. We have separated six nearby anti-centre clouds that are coherent in velocities and distances, from the Galactic-disc background in HI 21-cm and $^{12}$CO 2.6-mm line emission. We have jointly modelled the gamma-ray intensity recorded between 0.4 and 100 GeV, and the dust optical depth at 353 GHz as a combination of HI-bright, CO-bright, and ionised gas components. The complementary information from dust emission and gamma rays was used to reveal the gas not seen, or poorly traced, by HI , free-free, and $^{12}$CO emissions, namely (i) the opaque HI and diffuse H$_2$ present in the Dark Neutral Medium at the atomic-molecular transition, and (ii) the dense H$_2$ to be added where $^{12}$CO lines saturate. The measured interstellar gamma-ray spectra support a uniform penetration of the cosmic rays with energies above a few GeV through the clouds. We find a gradual increase in grain opacity as the gas becomes more dense. The increase reaches a factor of four to six in the cold molecular regions that are well shielded from stellar radiation. Consequently, the X$_{CO}$ factor derived from dust is systematically larger by 30% to 130% than the gamma-ray estimate. We also evaluate the average gamma-ray X$_{CO}$ factorfor each cloud, and find that X$_{CO}$ tends to decrease from diffuse to more compact molecular clouds, as expected from theory.

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

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two main components: the central (Gaussian) region (core component), and the power-law like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and its environmental column density were compared with theoretical models describing the evolution of filaments under gravity-dominated conditions. The feasibility of a transition to supercritical state by accretion is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high Mline-core) filaments in dense environments (av\sim3mag) can become supercritical in timescales of t\sim1Myr by accreting mass at constant or decreasing width. The trend of increasing Mline-tot (Mline-core and Mline-wing), and ridge Av with background also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the association between dense environments and high Mline-core values, and argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations.

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

IRIS (InfraRed Imaging Spectrograph) is the diffraction-limited first light instrument for the Thirty Meter Telescope (TMT) that consists of a near-infrared (0.84 to 2.4 $\mu$m) imager and integral field spectrograph (IFS). The IFS makes use of a lenslet array and slicer for spatial sampling, which will be able to operate in 100's of different modes, including a combination of four plate scales from 4 milliarcseconds (mas) to 50 mas with a large range of filters and gratings. The imager will have a field of view of 34$\times$34 arcsec$^{2}$ with a plate scale of 4 mas with many selectable filters. We present the preliminary design of the data reduction system (DRS) for IRIS that need to address all of these observing modes. Reduction of IRIS data will have unique challenges since it will provide real-time reduction and analysis of the imaging and spectroscopic data during observational sequences, as well as advanced post-processing algorithms. The DRS will support three basic modes of operation of IRIS; reducing data from the imager, the lenslet IFS, and slicer IFS. The DRS will be written in Python, making use of open-source astronomical packages available. In addition to real-time data reduction, the DRS will utilize real-time visualization tools, providing astronomers with up-to-date evaluation of the target acquisition and data quality. The quicklook suite will include visualization tools for 1D, 2D, and 3D raw and reduced images. We discuss the overall requirements of the DRS and visualization tools, as well as necessary calibration data to achieve optimal data quality in order to exploit science cases across all cosmic distance scales.

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

The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role of environment in the formation of supercritical filaments. Filaments were extracted from fields at D<500pc with the getfilaments algorithm and characterised according to their column density profiles and intrinsic properties. Each profile was fitted with a beam-convolved Plummer-like function and quantified based on the relative contributions from the filament 'core', represented by a Gaussian, and 'wing' component, dominated by the power-law of the Plummer-like function. These parameters were examined for populations associated with different background levels. We find that filaments increase their core (Mcore) and wing (Mwing) contributions while increasing their total linear mass density (Mtot). Both components appear to be linked to the local environment, with filaments in higher backgrounds having systematically more massive Mcore and Mwing. This dependence on the environment supports an accretion-based model for filament evolution in the local neighbourhood (D<500pc). Structures located in the highest backgrounds develop the highest central Av, Mcore, and Mwing as Mtot increases with time, favoured by the local availability of material and the enhanced gravitational potential. Our results indicate that filaments acquiring a significantly massive central region with Mcore>Mcrit/2 may become supercritical and form stars. This translates into a need for filaments to become at least moderately self-gravitating in order to undergo localised star formation or become star-forming filaments.

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

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

There is accumulating evidence for the presence of complex molecules, including carbon-bearing and organic molecules, in the interstellar medium. Much of this evidence comes to us from studies of chemical composition, photo- and mass-spectroscopy in cometary, meteoritic and asteroid samples, indicating a need to better understand the surface chemistry of astrophysical objects. There is also considerable interest in the origins of life-forming and life-sustaining molecules on Earth. Here, we perform reactive molecular dynamics simulations to probe the formation of carbon-rich molecules and clusters on carbonaceous surfaces resembling dust grains and meteoroids. Our results show that large chains form on graphitic surfaces at low temperatures (100K - 500K) and smaller fullerene-like molecules form at higher temperatures (2000K - 3000K). The formation is faster on the surface than in the gas at low temperatures but slower at high temperatures as surface interactions prevent small clusters from coagulation. We find that for efficient formation of molecular complexity, mobility about the surface is important and helps to build larger carbon chains on the surface than in the gas phase at low temperatures. Finally, we show that the temperature of the surface strongly determines what kind of structures forms and that low turbulent environments are needed for efficient formation.

The Galactic Cold Cores project has made Herschel observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. The dust opacity spectral index beta and the dust colour temperature T are derived using Herschel and Planck data. The relation between beta and T is examined for the whole sample and inside individual fields. Based on IRAS and Planck data, the fields are characterised by a median colour temperature of 16.1 K and a median opacity spectral index of beta=1.84. We observe a clear T-beta anti-correlation. In Herschel observations, constrained at lower resolution by Planck data, the variations follow the column density structure and beta(FIR) can rise to ~2.2 in individual clumps. The Planck 217 GHz band shows a systematic excess that is consistent with a general flattening of the dust emission spectrum at millimetre wavelengths. When fitted separately below and above 700 um, the median spectral index values are beta(FIR) ~ 1.91 and beta(mm) ~ 1.66. The spectral index changes as a function of column density and wavelength. Beta variations are partly masked by temperature gradients and the changes in the intrinsic grain properties may be even greater.

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.

Our velocity relative to the rest frame of the cosmic microwave background (CMB) generates a dipole temperature anisotropy on the sky which has been well measured for more than 30 years, and has an accepted amplitude of v/c = 0.00123, or v = 369km/s. In addition to this signal generated by Doppler boosting of the CMB monopole, our motion also modulates and aberrates the CMB temperature fluctuations (as well as every other source of radiation at cosmological distances). This is an order 0.1% effect applied to fluctuations which are already one part in roughly one hundred thousand, so it is quite small. Nevertheless, it becomes detectable with the all-sky coverage, high angular resolution, and low noise levels of the Planck satellite. Here we report a first measurement of this velocity signature using the aberration and modulation effects on the CMB temperature anisotropies, finding a component in the known dipole direction, (l,b)=(264, 48) [deg], of 384km/s +- 78km/s (stat.) +- 115km/s (syst.). This is a significant confirmation of the expected velocity.

It is widely accepted that cosmic rays (CRs) up to at least PeV energies are Galactic in origin. Accelerated particles are injected into the interstellar medium where they propagate to the farthest reaches of the Milky Way, including a surrounding halo. The composition of CRs coming to the solar system can be measured directly and has been used to infer the details of CR propagation that are extrapolated to the whole Galaxy. In contrast, indirect methods, such as observations of gamma-ray emission from CR interactions with interstellar gas, have been employed to directly probe the CR densities in distant locations throughout the Galactic plane. In this article we use 73 months of data from the Fermi Large Area Telescope in the energy range between 300 MeV and 10 GeV to search for gamma-ray emission produced by CR interactions in several high- and intermediate-velocity clouds located at up to ~ 7 kpc above the Galactic plane. We achieve the first detection of intermediate-velocity clouds in gamma rays and set upper limits on the emission from the remaining targets, thereby tracing the distribution of CR nuclei in the halo for the first time. We find that the gamma-ray emissivity per H atom decreases with increasing distance from the plane at 97.5% confidence level. This corroborates the notion that CRs at the relevant energies originate in the Galactic disk. The emissivity of the upper intermediate-velocity Arch hints at a 50% decline of CR densities within 2 kpc from the plane. We compare our results to predictions of CR propagation models.

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.

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

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to $n_s = 0.9603 \pm 0.0073$, ruling out exact scale invariance at over 5 $\sigma$. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining $d n_s/d ln k = -0.0134 \pm 0.0090$. Several analyses dropping the slow-roll approximation are carried out, including detailed model comparison and inflationary potential reconstruction. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit $\chi^2$ by ~ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with bounds on $f_\mathrm{NL}$ measured by Planck. The fractional primordial contribution of CDM isocurvature modes in the curvaton and axion scenarios has upper bounds of 0.25% or 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelation can improve $\chi^2$ by approximatively 4 as a result of slightly lowering the theoretical prediction for the $\ell<40$ multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

The project Galactic Cold Cores has made Herschel observations of interstellar clouds where the Planck satellite survey has located cold and compact clumps. The sources range from starless clumps to protostellar cores. We examine 116 Herschel fields to estimate the submillimetre dust opacity and its variations. The submillimetre dust opacity was derived from Herschel data, and near-infrared observations of the reddening of background stars are converted into near-infrared optical depth. We studied the systematic errors affecting these parameters and used modelling to correct for the expected biases. The ratio of 250um and J band opacities is correlated with the cloud location and star formation activity. We find a median ratio of tau(250um)/tau(J)= (1.6+-0.2)*10^-3, which is more than three times the mean value in diffuse medium. Assuming a spectral index beta=1.8 instead of beta=2.0, the value would be lower by ~30%. No significant systematic variation is detected with Galactocentric distance or with Galactic height. The tau(250um)/tau(J) maps reveal six fields with clear increase of submillimetre opacity of up to tau(250um)/tau(J) ~ 4*10^-3. These are all nearby fields with spatially resolved clumps of high column density. We interpret the increase in the far-infrared opacity as a sign of grain growth in the densest and coldest regions of interstellar clouds.

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales ($\ell > 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim \theta \lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $\sigma_8(\Omega_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $\sigma_{8}$.

Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $\psi$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction $p_\mathrm{max}$ decreases with column density $N_\mathrm{H}$ in the more opaque fields with $N_\mathrm{H} > 10^{21}\,\mathrm{cm}^{-2}$; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical (MHD) turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction $p_\mathrm{max}$ in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function $\mathcal{S}$. [abridged]

This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

We discuss the methods employed to photometrically calibrate the data acquired by the Low Frequency Instrument on Planck. Our calibration is based on a combination of the Orbital Dipole plus the Solar Dipole, caused respectively by the motion of the Planck spacecraft with respect to the Sun and by motion of the Solar System with respect to the CMB rest frame. The latter provides a signal of a few mK with the same spectrum as the CMB anisotropies and is visible throughout the mission. In this data release we rely on the characterization of the Solar Dipole as measured by WMAP. We also present preliminary results (at 44GHz only) on the study of the Orbital Dipole, which agree with the WMAP value of the Solar System speed within our uncertainties. We compute the calibration constant for each radiometer roughly once per hour, in order to keep track of changes in the detectors' gain. Since non-idealities in the optical response of the beams proved to be important, we implemented a fast convolution algorithm which considers the full beam response in estimating the signal generated by the dipole. Moreover, in order to further reduce the impact of residual systematics due to sidelobes, we estimated time variations in the calibration constant of the 30GHz radiometers (the ones with the largest sidelobes) using the signal of a reference load. We have estimated the calibration accuracy in two ways: we have run a set of simulations to assess the impact of statistical errors and systematic effects in the instrument and in the calibration procedure, and we have performed a number of consistency checks on the data and on the brightness temperature of Jupiter. Calibration errors for this data release are expected to be about 0.6% at 44 and 70 GHz, and 0.8% at 30 GHz. (Abriged.)

On the arcminute angular scales probed by Planck, the CMB anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217GHz frequency bands with an overall significance of greater than 25sigma. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB last-scattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the best-fitting LCDM model for the Planck temperature power spectrum, showing that this measurement at z=1100 correctly predicts the properties of the lower-redshift, later-time structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter constraints; it improves CMB-alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z~2.

We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.7 to 4.6 arcmin. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12 and 39 microKelvin in HFI four lowest frequency channel (100--353 GHz) and 13 and 14 kJy/sr for the 545 and 857 GHz channels. Using the 143 GHz channel as a reference, these two high frequency channels are intercalibrated within 5% and the 353 GHz relative calibration is at the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50 < l <2500), are intercalibrated at better than 0.2 %.

Context. Using observations to deduce dust properties, grain size distribution, and physical conditions in molecular clouds is a highly degenerate problem. Aims. The coreshine phenomenon, a scattering process at 3.6 and 4.5 $\mu$m that dominates absorption, has revealed its ability to explore the densest parts of clouds. We want to use this effect to constrain the dust parameters. The goal is to investigate to what extent grain growth (at constant dust mass) inside molecular clouds is able to explain the coreshine observations. We aim to find dust models that can explain a sample of Spitzer coreshine data. We also look at the consistency with near-infrared data we obtained for a few clouds. Methods. We selected four regions with a very high occurrence of coreshine cases: Taurus-Perseus, Cepheus, Chameleon and L183/L134. We built a grid of dust models and investigated the key parameters to reproduce the general trend of surface bright- nesses and intensity ratios of both coreshine and near-infrared observations with the help of a 3D Monte-Carlo radiative transfer code. The grid parameters allow to investigate the effect of coagulation upon spherical grains up to 5 $\mu$m in size derived from the DustEm diffuse interstellar medium grains. Fluffiness (porosity or fractal degree), ices, and a handful of classical grain size distributions were also tested. We used the near- and mostly mid-infrared intensity ratios as strong discriminants between dust models. Results. The determination of the background field intensity at each wavelength is a key issue. In particular, an especially strong background field explains why we do not see coreshine in the Galactic plane at 3.6 and 4.5 $\mu$m. For starless cores, where detected, the observed 4.5 $\mu$m / 3.6 $\mu$m coreshine intensity ratio is always lower than $\sim$0.5 which is also what we find in the models for the Taurus-Perseus and L183 directions. Embedded sources can lead to higher fluxes (up to four times greater than the strongest starless core fluxes) and higher coreshine ratios (from 0.5 to 1.1 in our selected sample). Normal interstellar radiation field conditions are sufficient to find suitable grain models at all wavelengths for starless cores. The standard interstellar grains are not able to reproduce observations and, due to the multi-wavelength approach, only a few grain types meet the criteria set by the data. Porosity does not affect the flux ratios while the fractal dimension helps to explain coreshine ratios but does not seem able to reproduce near-infrared observations without a mix of other grain types. Conclusions. Combined near- and mid-infrared wavelengths confirm the potential to reveal the nature and size distribution of dust grains. Careful assessment of the environmental parameters (interstellar and background fields, embedded or nearby reddened sources) is required to validate this new diagnostic.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.