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

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

Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ~ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critical photoelectric heating mechanism in interstellar gas, and [C II] 158{\mu}m fine-structure emission, a principal coolant. We present ALMA Band 9 observations targeting [C II] in six z ~ 2 star-forming galaxies with strong Spitzer IRS detections of PAH emission. All six galaxies are detected in dust continuum and marginally resolved. We compare the properties of PAH and [C II] emission, and constrain their relationship as a function of total infrared luminosity (LIR) and IR surface density. [C II] emission is detected in one galaxy at high signal-to-noise (34{\sigma}), and we place a secure upper limit on a second source. The rest of our sample are not detected in [C II] likely due to redshift uncertainties and narrow ALMA bandpass windows. Our results are consistent with the deficit in [C II]/LIR and PAH/LIR observed in the literature. However, the ratio of [C II] to PAH emission at z ~ 2 is possibly much lower than what is observed in nearby dusty star-forming galaxies. This could be the result of enhanced cooling via [O I] at high-z, hotter gas and dust temperatures, and/or a reduction in the photoelectric efficiency, in which the coupling between interstellar radiation and gas heating is diminished.

Cosmology and galaxy evolution studies with LSST, Euclid, and WFIRST, will require accurate redshifts for the detected galaxies. In this study, we present improved photometric redshift estimates for galaxies using a template library that populates three-color space and is constrained by HST/CANDELS photometry. For the training sample, we use a sample of galaxies which have median photometric redshifts derived by combining the results of different redshift-fitting codes. It allows us to train on a large, unbiased galaxy sample having deep, unconfused photometry at optical-to-mid infrared wavelengths. Galaxies in the training sample are assigned to cubes in three-dimensional color space, V-H, I-J, and z-H. We then derive the best-fit spectral energy distributions (SEDs) of the training sample at the fixed CANDELS median photometric redshifts to construct the new template library for each individual color cube (i.e. color-cube-based template library). We derive photometric redshifts (photo-z) of our target galaxies in the GOODS-S and -N fields using our new color-cube-based template library. We find that our method yields $\sigma_{NMAD}$ of 0.026 and an outlier fraction of 6% using only photometry in the LSST and Euclid/WFIRST bands. This is an improvement of 10% on $\sigma_{NMAD}$ and a reduction in outlier fraction of 13% compared to other techniques. In particular, we improve the photo-z precision by about 30% at 2 < z < 3. We also assess photo-z improvements by including K or mid-infrared bands to the ugrizY JH photometry. Our color-cube-based template library is a powerful tool to constrain photometric redshifts for future large surveys.

Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

We report the discovery of a mid-infrared variable AGN which is hosted by an ultraluminous infrared galaxy (ULIRG) candidate in the Sloan Stripe 82 field. \textit{WISE} \textit{J030654.88+010833.6} is a red, extended galaxy, which we estimate to be at a photometric redshift of 0.28 $\leq$ z $\leq$ 0.31, based on its optical and near-infrared spectral energy distribution (SED). The factor of two variability over 8 years seen in the \textit{WISE} 3.4 and 4.6 $\mu$m wavelength channels is not clearly correlated with optical variability in archival data. Based on our estimation of the physical parameters of the host galaxy, \textit{J030654.88+010833.6} is possibly a composite AGN/starburst ULIRG in a phase where high star formation $\sim$ 70 M$_{\odot}$ year$^{-1}$ is occurring. Our estimate of the black hole mass to stellar mass ratio also appears to be consistent with that of broad-line AGN in the local universe. The long-term variability of \textit{J030654.88+010833.6} as seen in the \textit{WISE} \textit{W1} and \textit{W2} light curves is likely due to variations in the accretion rate, with the energy being reprocessed by dust in the vicinity of the AGN.

We take advantage of new 100-hour Spitzer/IRAC observations available for MACS0647-JD, a strongly lensed $z\approx11$ galaxy candidate, to provide improved constraints on its physical properties. Probing the physical properties of galaxies at $z>8$ is challenging due to the inherent faintness of such sources and the limited wavelength coverage available. Thanks to the high $\approx$2-6$\times$ lensing magnification of the multiple images of MACS0647-JD, we can use the sensitive Spitzer/IRAC data to probe the rest-frame optical fluxes of MACS0647-JD and investigate its physical properties including the age and the stellar mass. In deriving Spitzer/IRAC fluxes for MACS0647-JD, great care is taken in coping with the impact of three bright ($\approx$8-16 mag) stars in our field to ensure robust results. Assuming a constant star formation rate, the age, stellar mass, and rest-frame UV slope we estimate for MACS0647-JD based on a stack of the photometry are log$_{10}$(age/yr) = 8.6$^{+0.1}_{-2.1}$, log$_{10}$(M$_{*}$/M$_{\odot}$) = 9.1$^{+0.2}_{-1.4}$, and $\beta = -$1.3$\pm$0.6, respectively. We compare our results with expectations from the EAGLE simulation and find that MACS0647-JD has properties consistent with corresponding to the most massive and rapidly star-forming galaxies in the simulation. We also find that its radius, 105$\pm$28 pc, is a factor of $\approx$2 smaller than the mean size in a separate simulation project DRAGONS. Interestingly enough, the observed size is similar to the small sizes seen in very low-luminosity $z\approx$6-10 galaxies behind lensing clusters.

The Cosmic Dawn Intensity Mapper (CDIM) will transform our understanding of the era of reionization when the Universe formed the first stars and galaxies, and UV photons ionized the neutral medium. CDIM goes beyond the capabilities of upcoming facilities by carrying out wide area spectro-imaging surveys, providing redshifts of galaxies and quasars during reionization as well as spectral lines that carry crucial information on their physical properties. CDIM will make use of unprecedented sensitivity to surface brightness to measure the intensity fluctuations of reionization on large-scales to provide a valuable and complementary dataset to 21-cm experiments. The baseline mission concept is an 83-cm infrared telescope equipped with a focal plane of 24 \times 20482 detectors capable of R = 300 spectro-imaging observations over the wavelength range of 0.75 to 7.5 {\mu}m using Linear Variable Filters (LVFs). CDIM provides a large field of view of 7.8 deg2 allowing efficient wide area surveys, and instead of moving instrumental components, spectroscopic mapping is obtained through a shift-and-stare strategy through spacecraft operations. CDIM design and capabilities focus on the needs of detecting faint galaxies and quasars during reionization and intensity fluctuation measurements of key spectral lines, including Lyman-{\alpha} and H{\alpha} radiation from the first stars and galaxies. The design is low risk, carries significant science and engineering margins, and makes use of technologies with high technical readiness level for space observations.

Galaxies form and evolve in the context of their local and large-scale environments. Their baryonic content that we observe with imaging and spectroscopy is intimately connected to the properties of their dark matter halos, and to their location in the "cosmic web" of large-scale structure. Very large spectroscopic surveys of the local universe (e.g., SDSS and GAMA) measure galaxy positions (location within large-scale structure), statistical clustering (a direct constraint on dark matter halo masses), and spectral features (measuring physical conditions of the gas and stars within the galaxies, as well as internal velocities). Deep surveys with the James Webb Space Telescope (JWST) will revolutionize spectroscopic measurements of redshifts and spectral properties for galaxies out to the epoch of reionization, but with numerical statistics and over cosmic volumes that are too small to map large-scale structure and to constrain halo properties via clustering. Here, we consider advances in understanding galaxy evolution that would be enabled by very large spectroscopic surveys at high redshifts: very large numbers of galaxies (outstanding statistics) over large co-moving volumes (large-scale structure on all scales) over broad redshift ranges (evolution over most of cosmic history). The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

Astrophysical measurements away from the 1 AU orbit of Earth can enable several astrophysical science cases that are challenging or impossible to perform from Earthbound platforms, including: building a detailed understanding of the extragalactic background light throughout the electromagnetic spectrum; measurements of the properties of dust and ice in the inner and outer solar system; determinations of the mass of planets and stellar remnants far from luminous stars using gravitational microlensing; and stable time-domain astronomy. Though potentially transformative for astrophysics, opportunities to fly instrumentation capable of these measurements are rare, and a mission to the distant solar system that includes instrumentation expressly designed to perform astrophysical science, or even one primarily for a different purpose but capable of precise astronomical investigation, has not yet been flown. In this White Paper, we describe the science motivations for this kind of measurement, and advocate for future flight opportunities that permit intersectional collaboration and cooperation to make these science investigations a reality.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

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 an exploration of the expected detection of the earliest Active Galactic Nuclei (AGN) in the Universe from state-of-art galaxy formation and evolution semi-analytic models and hydro-dynamical simulations. We estimate the number and radiative characteristics of Super Massive Black Holes (SMBHs) at $z\geq 6$, a redshift range that will be intensively explored by the next generation of telescopes, in particular in the radio through the Square Kilometre Array (SKA) and at high energies with ESA's Athena X-ray Observatory. We find that Athena will be able to observe over 5000 AGN/$\rm deg^2$ at the Epoch of Re-ionization (EoR), $6\leq z \leq 10$. Similarly, for the same redshift range the models/simulations suggest that SKA will detect at least 400 AGN/$\rm deg^2$. Additionally, we stress the importance of the volume of the simulation box as well as the initial physical conditions of the models/simulations on their effect on the luminosity functions (LFs) and the creation of the most massive SMBHs that we currently observe at the EoR. Furthermore, following the evolution of the accretion mode of the SMBHs in each model/simulation, we show that, while the quasar dominates over the radio mode at the EoR, detection at radio wavelengths still reaches significant numbers even at the highest redshifts. Finally, we present the effect that the radiative efficiency has on the LFs by comparing results produced with a constant value for the radiative efficiency and more complex calculations based on the spin of each SMBH.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way.

We report our Spitzer Space Telescope observations and detections of the binary neutron star merger GW170817. At 4.5um, GW170817 is detected at 21.9 mag AB at +43 days and 23.9 mag AB at +74 days after merger. At 3.6um, GW170817 is not detected to a limit of 23.2 mag AB at +43 days and 23.1 mag AB at +74 days. Our detections constitute the latest and reddest constraints on the kilonova/macronova emission and composition of heavy elements. The 4.5um luminosity at this late phase cannot be explained by elements exclusively from the first abundance peak of the r-process. Moreover, the steep decline in the Spitzer band, with a power-law index of 3.4 +/- 0.2, can be explained by a few of the heaviest isotopes in the third abundance peak with half-life around 14 days dominating the luminosity (e.g. 140Ba, 143Pr, 147Nd, 156Eu, 191Os, 223Ra, 225Ra, 233Pa, 234Th) or a model with lower deposition efficiency. This data offers evidence that the heaviest elements in the second and third r-process abundance peak were indeed synthesized. Our conclusion is verified by both analytics and network simulations and robust despite intricacies and uncertainties in the nuclear physics. Future observations with Spitzer and James Webb Space Telescope will further illuminate the relative abundance of the synthesized heavy elements.

The outer solar system provides a unique, quiet vantage point from which to observe the universe around us, where measurements could enable several niche astrophysical science cases that are too difficult to perform near Earth. NASA's New Horizons mission comprises an instrument package that provides imaging capability from ultraviolet (UV) to near-infrared (near-IR) wavelengths with moderate spectral resolution located beyond the orbit of Pluto. A carefully designed survey with New Horizons can optimize the use of expendable propellant and the limited data telemetry bandwidth to allow several measurements, including a detailed understanding of the cosmic extragalactic background light; studies of the local and extragalactic UV background; measurements of the properties of dust and ice in the outer solar system; confirmation and characterization of transiting exoplanets; determinations of the mass of dark objects using gravitational microlensing; and rapid follow-up of transient events. New Horizons is currently in an extended mission designed to focus on Kuiper Belt science that will conclude in 2021. The astrophysics community has a unique, generational opportunity to use this mission for astronomical observation at heliocentric distances beyond 50 au in the next decade. In this paper, we discuss the potential science cases for such an extended mission, and provide an initial assessment of the most important operational requirements and observation strategies it would require. We conclude that New Horizons is capable of transformative science, and that it would make a valuable and unique asset for astrophysical science that is unlikely to be replicated in the near future.

We investigate the value of the near-infrared imaging from upcoming surveys for constraining the ellipticities of galaxies. We select galaxies between 0.5 < z < 3 that are brighter than expected Euclid sensitivity limits from the GOODS-S and N fields in CANDELS. The co-added CANDELS/HST V+I and J+H images are degraded in resolution and sensitivity to simulate Euclid-quality optical and near-infrared (NIR) images. We then run GALFIT on these simulated images and find that optical and NIR provide similar performance in measuring galaxy ellipticities at redshifts 0.5 < z < 3. At z > 1.0, the NIR-selected source density is higher by a factor of 1.4 and therefore the standard error in NIR-derived ellipticities is about 30% smaller, implying a more precise ellipticity measurement. The good performance of the NIR is mainly because galaxies have an intrinsically smoother light distribution in the NIR bands than in the optical, the latter tracing the clumpy star-forming regions. In addition, the NIR bands have a higher surface brightness per pixel than the optical images, while being less affected by dust attenuation. Despite the worse spatial sampling and resolution of Euclid NIR compared to optical, the NIR approach yields equivalent or more precise galaxy ellipticity measurements. If systematics that affect shape such as dithering strategy and point spread function undersampling can be mitigated, inclusion of the NIR can improve galaxy ellipticity measurements over all redshifts. This is particularly important for upcoming weak lensing surveys, such as with Euclid and WFIRST.

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

We report the detection of linearly polarized emission at 53 and 89 $\mu$m, from the radio-loud active galactic nucleus (AGN) Cygnus A using HAWC+ onboard SOFIA. We measure a highly polarized core of $11\pm3$% and $9\pm2$% with a position angle (P.A.) of polarization of $43\pm8^{\circ}$ and $39\pm7^{\circ}$ at 53 and 89 $\mu$m, respectively. We find (1) a synchrotron dominated core with a flat spectrum ($+0.21\pm0.05$) and a turn-over at $543\pm120$ $\mu$m, which implies synchrotron emission is insignificant in the infrared (IR), and (2) a $2-500$ $\mu$m bump peaking at $\sim40$ $\mu$m described by a blackbody component with color temperature of $107\pm9$ K. The polarized SED has the same shape as the IR bump of the total flux SED. We observe a change in the P.A. of polarization of $\sim20^{\circ}$ from 2 to 89 $\mu$m, which suggests a change of polarization mechanisms. The ultraviolet, optical and near-IR polarization has been convincingly attributed to scattering by polar dust, consistent with the usual torus scenario, though this scattered component can only be directly observed from the core in the near-IR. By contrast, the gradual rotation by $\sim20^{\circ}$ towards the far-IR, and the near-perfect match between the total and polarized IR bumps, indicate that dust emission from aligned dust grains becomes dominant at $10-100$ $\mu$m, with a large polarization of 10% at a nearly constant P.A. This result suggests that a coherent dusty and magnetic field structure dominates the $10-100$ $\mu$m emission around the AGN.

We conduct a stacking analysis using 1.4 GHz NRAO VLA Sky Survey (NVSS) detections and Planck all-sky maps to estimate the differential source counts down to the few 100 $\mu$Jy level at 30, 44, 70 and 100 GHz. Consequently, we are able to measure the integrated extragalactic background light from discrete sources at these frequencies. By integrating down to a 1.4 GHz flux density of $\approx$2$ \mu$Jy, we measure integrated, extragalactic brightness temperatures from discrete sources of $105.63\pm10.56 $mK, $21.76\pm3.09 \mu$K, $8.80\pm0.95 \mu$K, $2.59\pm0.27 \mu$K, and $1.15\pm0.10 \mu$k at 1.4, 30, 44, 70, and 100 GHz, respectively. Our measurement at 1.4 GHz is slightly larger than previous measurements, most likely due to using NVSS data compared to older interferometric data in the literature, but still remains a factor of $\approx$4.5 below that required to account for the excess extragalactic sky brightness measured at 1.4 GHz by ARCADE 2. The fit to ARCADE 2 total extragalactic sky brightness measurements is also a factor of $\approx$8.6, 6.6, 6.2, and 4.9 times brighter than what we estimate from discrete sources at 30, 44, 70 and 100 GHz, respectively. The extragalactic sky spectrum (i.e., $T_{\rm b} \propto \nu^{\beta}$) from discrete sources appears to flatten with increasing frequency, having a spectral index of $\beta=-2.82\pm0.06$ between 1.4 and 30 GHz and $\beta=-2.39\pm0.12$ between 30 and 100 GHz. We believe that the spectral flattening most likely arises from a combination of Gigahertz-peaked sources and the spectral hardening of radio-detected sources at higher frequencies, particularly at faint flux densities. However, the precise origin of a hard component of energetic electrons responsible for the emission remains unclear.

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.

We characterize the influence that inclination has on the shape and normalization in average dust attenuation curves derived from a sample of ~10,000 local star-forming galaxies. To do this we utilize aperture-matched multi-wavelength data from the Galaxy Evolution Explorer, the Sloan Digital Sky Survey, the United Kingdom Infrared Telescope, and the Two Micron All-Sky Survey. We separate our sample into groups according to axial ratio (b/a) and obtain an estimate of their average total-to-selective attenuation $k(\lambda)$. The attenuation curves are found to be shallower at UV wavelengths with increasing inclination, whereas the shape at longer wavelengths remains unchanged. The highest inclination subpopulation (b/a<0.42) exhibits a NUV excess in its average selective attenuation, which, if interpreted as a 2175\AA\ feature, is best fit with a bump strength of 17-26% of the MW value. No excess is apparent in the average attenuation curve of lower inclination galaxies. The differential reddening between the stellar continuum and ionized gas is found to decrease with increasing inclination. We find that higher inclination galaxies have slightly higher values of $R_V$, although this is poorly constrained given the uncertainties. We outline possible explanations for these trends within a two component dust model (dense cloud+diffuse dust) and find that they can be naturally explained if carriers of the 2175\AA\ feature are preferentially destroyed in star-forming regions (UV-bright regions).

We present Gemini and Keck spectroscopic redshifts and velocity dispersions for twenty clusters detected via the Sunyaev-Zel'dovich (SZ) effect by the Planck space mission, with estimated masses in the range $2.3 \times 10^{14} M_{\odot} < M < 9.4 \times 10^{14} M_{\odot}$. Cluster members were selected for spectroscopic follow-up with Palomar, Gemini and Keck optical and (in some cases) infrared imaging. Seven cluster redshifts were measured for the first time with this observing campaign, including one of the most distant Planck clusters confirmed to date, at $z=0.782\pm0.010$, PSZ2 G085.95+25.23. The spectroscopic redshift catalogs of members of each confirmed cluster are included as on-line tables. We show the galaxy redshift distributions and measure the cluster velocity dispersions. The cluster velocity dispersions obtained in this paper were used in a companion paper to measure the Planck mass bias and to constrain the cluster velocity bias.

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

We present 1-second cadence observations of M32 (NGC221) with the CHIMERA instrument at the Hale 200-inch telescope of the Palomar Observatory. Using field stars as a baseline for relative photometry, we are able to construct a light curve of the nucleus in the g-prime and r-prime band with 1sigma=36 milli-mag photometric stability. We derive a temporal power spectrum for the nucleus and find no evidence for a time-variable signal above the noise as would be expected if the nuclear black hole were accreting gas. Thus, we are unable to constrain the spin of the black hole although future work will use this powerful instrument to target more actively accreting black holes. Given the black hole mass of (2.5+/-0.5)*10^6 Msun inferred from stellar kinematics, the absence of a contribution from a nuclear time-variable signal places an upper limit on the accretion rate which is 4.6*10^{-8} of the Eddington rate, a factor of two more stringent than past upper limits from HST. The low mass of the black hole despite the high stellar density suggests that the gas liberated by stellar interactions was primarily at early cosmic times when the low-mass black hole had a small Eddington luminosity. This is at least partly driven by a top-heavy stellar initial mass function at early cosmic times which is an efficient producer of stellar mass black holes. The implication is that supermassive black holes likely arise from seeds formed through the coalescence of 3-100 Msun mass black holes that then accrete gas produced through stellar interaction processes.

We study the mass-richness relation of 116 spectroscopically-confirmed massive clusters at $0.4 < z < 2$ by mining the $Spitzer$ archive. We homogeneously measure the richness at 4.5$\mu$m for our cluster sample within a fixed aperture of $2^{\prime}$ radius and above a fixed brightness threshold, making appropriate corrections for both background galaxies and foreground stars. We have two subsamples, those which have a) literature X-ray luminosities and b) literature Sunyaev-Zeldovich effect masses. For the X-ray subsample we re-derive masses adopting the most recent calibrations. We then calibrate an empirical mass-richness relation for the combined sample spanning more than one decade in cluster mass and find the associated uncertainties in mass at fixed richness to be $\pm 0.25$ dex. We study the dependance of the scatter of this relation with galaxy concentration, defined as the ratio between richness measured within an aperture radius of 1 and 2 arcminutes. We find that at fixed aperture radius the scatter increases for clusters with higher concentrations. We study the dependance of our richness estimates with depth of the [4.5]$\mu$m imaging data and find that reaching a depth of at least [4.5]= 21 AB mag is sufficient to derive reasonable mass estimates. We discuss the possible extension of our method to the mid-infrared $WISE$ all-sky survey data, and the application of our results to the $Euclid$ mission. This technique makes richness-based cluster mass estimates available for large samples of clusters at very low observational cost.

We measure the Planck cluster mass bias using dynamical mass measurements based on velocity dispersions of a subsample of 17 Planck-detected clusters. The velocity dispersions were calculated using redshifts determined from spectra obtained at Gemini observatory with the GMOS multi-object spectrograph. We correct our estimates for effects due to finite aperture, Eddington bias and correlated scatter between velocity dispersion and the Planck mass proxy. The result for the mass bias parameter, $(1-b)$, depends on the value of the galaxy velocity bias $b_v$ adopted from simulations: $(1-b)=(0.51\pm0.09) b_v^3$. Using a velocity bias of $b_v=1.08$ from Munari et al., we obtain $(1-b)=0.64\pm 0.11$, i.e, an error of 17% on the mass bias measurement with 17 clusters. This mass bias value is consistent with most previous weak lensing determinations. It lies within $1\sigma$ of the value needed to reconcile the Planck cluster counts with the Planck primary CMB constraints. We emphasize that uncertainty in the velocity bias severely hampers precision measurements of the mass bias using velocity dispersions. On the other hand, when we fix the Planck mass bias using the constraints from Penna-Lima et al., based on weak lensing measurements, we obtain a positive velocity bias $b_v \gtrsim 0.9$ at $3\sigma$.

We characterize the near-infrared (NIR) dust attenuation for a sample of ~5500 local (z<0.1) star-forming galaxies and obtain an estimate of their average total-to-selective attenuation $k(\lambda)$. We utilize data from the United Kingdom Infrared Telescope (UKIRT) and the Two Micron All-Sky Survey (2MASS), which is combined with previously measured UV-optical data for these galaxies. The average attenuation curve is slightly lower in the far-UV than local starburst galaxies, by roughly 15%, but appears similar at longer wavelengths with a total-to-selective normalization at V-band of $R_V=3.67\substack{+0.44 \\ -0.35}$. Under the assumption of energy balance, the total attenuated energy inferred from this curve is found to be broadly consistent with the observed infrared dust emission ($L_{\rm{TIR}}$) in a small sample of local galaxies for which far-IR measurements are available. However, the significant scatter in this quantity among the sample may reflect large variations in the attenuation properties of individual galaxies. We also derive the attenuation curve for sub-populations of the main sample, separated according to mean stellar population age (via $D_n4000$), specific star formation rate, stellar mass, and metallicity, and find that they show only tentative trends with low significance, at least over the range which is probed by our sample. These results indicate that a single curve is reasonable for applications seeking to broadly characterize large samples of galaxies in the local Universe, while applications to individual galaxies would yield large uncertainties and is not recommended.

We directly detect dust emission in an optically-detected, multiply-imaged galaxy lensed by the Frontier Fields cluster MACSJ0717.5+3745. We detect two images of the same galaxy at 1.1mm with the AzTEC camera on the Large Millimeter Telescope leaving no ambiguity in the counterpart identification. This galaxy, MACS071_Az9, is at z>4 and the strong lensing model (mu=7.5) allows us to calculate an intrinsic IR luminosity of 9.7e10 Lsun and an obscured star formation rate of 14.6 +/- 4.5 Msun/yr. The unobscured star formation rate from the UV is only 4.1 +/- 0.3 Msun/yr which means the total star formation rate (18.7 +/- 4.5 Msun/yr) is dominated (75-80%) by the obscured component. With an intrinsic stellar mass of only 6.9e9Msun, MACS0717_Az9 is one of only a handful of z>4 galaxies at these lower masses that is detected in dust emission. This galaxy lies close to the estimated star formation sequence at this epoch. However, it does not lie on the dust obscuration relation (IRX-beta) for local starburst galaxies and is instead consistent with the Small Magellanic Cloud (SMC) attenuation law. This remarkable lower mass galaxy showing signs of both low metallicity and high dust content may challenge our picture of dust production in the early Universe.

(Abridged) Our sensitive ($\sigma_{\rm n} = 572\,{\rm nJy\,beam}^{-1}$), high-resolution (FWHM $\theta_{1/2} = 220\,{\rm mas} \approx2\mathrm{\,kpc~at~}z\gtrsim1$) 10$\,$GHz image covering a single Karl G.~Jansky Very Large Array (VLA) primary beam (FWHM $\Theta_{1/2} = 4.25'$) in the GOODS-N field contains 32 sources with $S_{\rm p}\gtrsim2\,\mu{\rm Jy~beam}^{-1}$ and optical and/or near-infrared (OIR) counterparts. Most are about as large as the star-forming regions that power them. Their median FWHM major axis is $\langle\theta_{\rm M} \rangle=167\pm32\,{\rm mas} \approx 1.2\pm0.28\,{\rm kpc}$ with rms scatter 91 was $\approx$ 0.79 kpc. In units of the effective radius $r_{\rm e}$ that encloses half their flux, these radio sizes are $\langle r_{\rm e}\rangle = 69\pm13{\rm mas} \approx \pm114\mathrm{\,pc}$ and have rms scatter $38\mathrm{\,mas}\approx324\mathrm{\,pc}$. These sizes are smaller than those measured at lower radio frequencies, but agree with dust emission sizes measured at mm/sub-mm wavelengths and extinction-corrected H$\alpha$ sizes. We made a low-resolution ($\theta_{1/2}=1.0"$) image with $\approx10\times$ better brightness sensitivity to detect extended sources and measure matched-resolution spectral indices $\alpha_{1.4}^{10}$. It contains 6 new sources with $S_{\rm p}\gtrsim3.9\,\mu{\rm Jy~beam}^{-1}$ and OIR counterparts. The median redshift of all 38 sources is $1.24\pm0.15$. The 19 sources with 1.4$\,$GHz counterparts have median spectral index $-0.74\pm0.10$ with rms scatter $0.35$. Including upper limits on $\alpha$ for sources not detected at 1.4$\,$GHz flattens the median to $\gtrsim-0.61$, suggesting that the $\mu$Jy radio sources at higher redshifts, and hence selected at higher rest-frame frequencies, may have flatter spectra. If the non-thermal spectral index is -0.85, the median thermal fraction at rest-frame frequency 20$\,$GHz is $\gtrsim$48%.

We review the uncertainties in high-z star-formation rate (SFR) measures and the constraints that one obtains from high-z gamma-ray burst (GRB) rates on them. We show that at the present time, the GRB rates per unit star-formation at z>3 are higher than at lower redshift. We also compare metallicity predictions made using a hierarchical model of cosmic chemical evolution based on two recently proposed SFRs, one based on the observed galaxy luminosity function at high redshift and one based on the GRB rate and find that within the considerable scatter in metal abundance measures, they both are consistent with the data. Analyzing the ensemble of different measurements together, we conclude that despite metallicity biases, GRBs may be a less biased probe of star-formation at z>3 than at z<2. There is likely to be a common origin to the high GRB rate per unit star-formation and the high observed Lyman-continuum production rate in high redshift galaxies and that this may be due to a relative overabundance of stars with mass $>$25Msun which are likely GRB progenitors. We also find that to reconcile these measurements with the Thomson scattering cross section of cosmic microwave background (CMB) photons measured by Planck, the escape fraction of Lyman-continuum photons from galaxies must be low, about 15 percent or less and that the clumping factor of the IGM is likely to be small, about 3. Finally, we demonstrate that GRBs are unique probes of metallicity evolution in low-mass galaxy samples and that GRB hosts likely lost a significant fraction of metals to the intergalactic medium (IGM) due to feedback processes such as stellar winds and supernovae.

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

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

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

The Planck Early Release Compact Source Catalogue (ERCSC) includes nine lists of highly reliable sources, individually extracted at each of the nine Planck frequency channels. To facilitate the study of the Planck sources, especially their spectral behaviour across the radio/infrared frequencies, we provide a "bandmerged" catalogue of the ERCSC sources. This catalogue consists of 15191 entries, with 79 sources detected in all nine frequency channels of Planck and 6818 sources detected in only one channel. We describe the bandmerging algorithm, including the various steps used to disentangle sources in confused regions. The multi-frequency matching allows us to develop spectral energy distributions of sources between 30 and 857 GHz, in particular across the 100 GHz band, where the energetically important CO J=1->0 line enters the Planck bandpass. We find ~3-5sigma evidence for contribution to the 100 GHz intensity from foreground CO along the line of sight to 147 sources with |b|>30 deg. The median excess contribution is 4.5+/-0.9 percent of their measured 100 GHz flux density which cannot be explained by calibration or beam uncertainties. This translates to 0.5+/-0.1 K km s^{-1} of CO which must be clumped on the scale of the Planck 100 GHz beam, i.e., ~10 arcmin. If this is due to a population of low mass (~15 Msun) molecular gas clumps, the total mass in these clumps may be more than 2000 Msun. Further, high-spatial-resolution, ground-based observations of the high-latitude sky will help shed light on the origin of this diffuse, clumpy CO emission.

Cosmic Dawn Intensity Mapper is a "Probe Class" mission concept for reionization studies of the universe. It will be capable of spectroscopic imaging observations between 0.7 to 6-7 microns in the near-Infrared. The primary observational objective is pioneering observations of spectral emission lines of interest throughout the cosmic history, but especially from the first generation of distant, faint galaxies when the universe was less than 800 million years old. With spectro-imaging capabilities, using a set of linear variable filters (LVFs), CDIM will produce a three-dimensional tomographic view of the epoch of reionization (EoR). CDIM will also study galaxy formation over more than 90% of the cosmic history and will move the astronomical community from broad-band astronomical imaging to low-resolution (R=200-300) spectro-imaging of the universe.

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 the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with telescopes at the Canary Islands observatories as part of the general optical follow-up programme undertaken by the Planck collaboration. In total, 78 SZ sources are discussed. Deep-imaging observations were obtained for most of these sources; spectroscopic observations in either in long-slit or multi-object modes were obtained for many. We effectively used 37.5 clear nights. We found optical counterparts for 73 of the 78 candidates. This sample includes 53 spectroscopic redshift determinations, 20 of them obtained with a multi-object spectroscopic mode. The sample contains new redshifts for 27 Planck clusters that were not included in the first Planck SZ source catalogue (PSZ1).

We introduce the Swift Gamma-Ray Burst Host Galaxy Legacy Survey ("SHOALS"), a multi-observatory high-redshift galaxy survey targeting the largest unbiased sample of long-duration gamma-ray burst hosts yet assembled (119 in total). We describe the motivations of the survey and the development of our selection criteria, including an assessment of the impact of various observability metrics on the success rate of afterglow-based redshift measurement. We briefly outline our host-galaxy observational program, consisting of deep Spitzer/IRAC imaging of every field supplemented by similarly-deep, multi-color optical/NIR photometry, plus spectroscopy of events without pre-existing redshifts. Our optimized selection cuts combined with host-galaxy follow-up have so far enabled redshift measurements for 110 targets (92%) and placed upper limits on all but one of the remainder. About 20% of GRBs in the sample are heavily dust-obscured, and at most 2% originate from z>5.5. Using this sample we estimate the redshift-dependent GRB rate density, showing it to peak at z~2.5 and fall by about an order of magnitude towards low (z=0) redshift, while declining more gradually towards high (z~7) redshift. This behavior is consistent with a progenitor whose formation efficiency varies modestly over cosmic history. Our survey will permit the most detailed examination to date of the connection between the GRB host population and general star-forming galaxies, directly measure evolution in the host population over cosmic time and discern its causes, and provide new constraints on the fraction of cosmic star-formation occurring in undetectable galaxies at all redshifts.

We present rest-frame NIR luminosities and stellar masses for a large and uniformly-selected population of GRB host galaxies using deep Spitzer Space Telescope imaging of 119 targets from the Swift GRB Host Galaxy Legacy Survey spanning 0.03 < z < 6.3, and determine the effects of galaxy evolution and chemical enrichment on the mass distribution of the GRB host population across cosmic history. We find strong evolution in the host luminosity distribution between z~0.5 (median absolute NIR AB magnitude ~ -18.5, corresponding to M* ~ 3x10^8 M_sun and z~1.5), but negligible variation between z~1.5 and z~5 (median magnitude ~ -21.2, corresponding to M* ~ 5x10^9 M_sun). Dust-obscured GRBs dominate the massive host population but are only rarely seen associated with low-mass hosts, indicating that massive star-forming galaxies are universally and (to some extent) homogeneously dusty at high-redshift while low-mass star-forming galaxies retain little dust in their ISM. Comparing our luminosity distributions to field surveys and measurements of the high-z mass-metallicity relation, our results have good consistency with a model in which the GRB rate per unit star-formation is constant in galaxies with gas-phase metallicity below approximately the Solar value but heavily suppressed in more metal-rich environments. This model also naturally explains the previously-reported "excess" in the GRB rate beyond z>2; metals stifle GRB production in most galaxies at z<1.5 but have only minor impact at higher redshifts. The metallicity threshold we infer is much higher than predicted by single-star models and favors a binary progenitor. Our observations also constrain the fraction of cosmic star-formation in low-mass galaxies undetectable to Spitzer to be a small minority at most redshifts (~10% at z~2, ~25% at z~3, and ~50% at z=3.5-6.0).

The dust attenuation for a sample of $\sim$10000 local ($z\lesssim0.1$) star forming galaxies is constrained as a function of their physical properties. We utilize aperture-matched multi-wavelength data available from the Galaxy Evolution Explorer (GALEX) and the Sloan Digital Sky Survey (SDSS) to ensure that regions of comparable size in each galaxy are being analyzed. We follow the method of Calzetti et al. (1994) and characterize the dust attenuation through the UV power-law index, $\beta$, and the dust optical depth, which is quantified using the difference in Balmer emission line optical depth, $\tau_B^l=\tau_{\mathrm{H}\beta} - \tau_{\mathrm{H}\alpha}$. The observed linear relationship between $\beta$ and $\tau_B^l$ is similar to the local starburst relation, but the large scatter ($\sigma_{\mathrm{int}}=0.44$) suggests there is significant variation in the local Universe. We derive a selective attenuation curve over the range $1250\mathrm{\AA}<\lambda<8320\mathrm{\AA}$ and find that a single attenuation curve is effective for characterizing the majority of galaxies in our sample. This curve has a slightly lower selective attenuation in the UV compared to previously determined curves. We do not see evidence to suggest that a 2175\AA\ feature is significant in the average attenuation curve. Significant positive correlations are seen between the amount of UV and optical reddening and galaxy metallicity, mass, star formation rate (SFR), and SFR surface density. This provides a potential tool for gauging attenuation where the stellar population is unresolved, such as at high-$z$.

The fine tuning of parameters required to reproduce our present day Universe suggests that our Universe may simply be a region within an eternally inflating super-region. Many other regions beyond our observable Universe would exist with each such region governed by a different set of physical parameters. Collision between these regions, if they occur, should leave signatures of anisotropy in the cosmic microwave background (CMB) but have not been seen. We analyze the spectral properties of masked, foreground-cleaned maps between 100 and 545 GHz constructed from the Planck dataset. Four distinct $\sim2-4\arcdeg$ regions associated with CMB cold spots show anomalously strong 143 GHz emission but no correspondingly strong emission at either 100 or 217 GHz. The signal to noise of this 143 GHz residual emission is at the $\gtrsim$6$\sigma$ level which reduces to $3.2-5.4\sigma$ after subtraction of remaining synchrotron/free-free foregrounds. We assess different mechanisms for this residual emission and conclude that although there is a 30\% probability that noise fluctuations may cause foregrounds to fall within 3$\sigma$ of the excess, there is less than a 0.5\% probability that foregrounds can explain all the excess. A plausible explanation is that the collision of our Universe with an alternate Universe whose baryon to photon ratio is a factor of $\sim$4500 larger than ours, could produce enhanced Hydrogen Paschen-series emission at the epoch of recombination. Future spectral mapping and deeper observations at 100 and 217 GHz are needed to mitigate systematics arising from unknown Galactic foregrounds and to confirm this unusual hypothesis.

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.

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

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

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

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

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We use a Bayesian software package to analyze CARMA-8 data towards 19 unconfirmed Planck SZ-cluster candidates from Rodriguez-Gonzalvez et al. (2015), that are associated with significant overdensities in WISE. We used two cluster parameterizations, one based on a (fixed shape) generalized-NFW pressure profile and another based on a beta-gas-density profile (with varying shape parameters) to obtain parameter estimates for the nine CARMA-8 SZ-detected clusters. We find our sample is comprised of massive, Y_{500}=0.0010 \pm 0.0015 arcmin^2, relatively compact, theta_{500}= 3.9 \pm 2.0 arcmin systems. Results from the beta model show that our cluster candidates exhibit a heterogeneous set of brightness-temperature profiles. Comparison of Planck and CARMA-8 measurements showed good agreement in Y_{500} and an absence of obvious biases. We estimated the total cluster mass M_{500} as a function of z for one of the systems; at the preferred photometric redshift of 0.5, the derived mass, M_{500} \approx 0.8 \pm 0.2 \times 10^{15} Msun. Spectroscopic Keck/MOSFIRE data confirmed a galaxy member of one of our cluster candidates to be at z=0.565. Applying a Planck prior in Y_{500} to the CARMA-8 results reduces uncertainties for both parameters by a factor >4, relative to the independent Planck or CARMA-8 measurements. We here demonstrate a powerful technique to find massive clusters at intermediate z \gtrsim 0.5 redshifts using a cross-correlation between Planck and WISE data, with high-resolution follow-up with CARMA-8. We also use the combined capabilities of Planck and CARMA-8 to obtain a dramatic reduction by a factor of several, in parameter uncertainties.

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

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

We present the results of deep \chandra\ imaging of the central region of the Extended Groth Strip, the AEGIS-X Deep (AEGIS-XD) survey. When combined with previous \chandra\ observations of a wider area of the strip, AEGIS-X Wide (AEGIS-XW; Laird et~al. 2009), these provide data to a nominal exposure depth of 800ks in the three central ACIS-I fields, a region of approximately $0.29$~deg$^{2}$. This is currently the third deepest X-ray survey in existence, a factor $\sim 2-3$ shallower than the Chandra Deep Fields (CDFs) but over an area $\sim 3$ times greater than each CDF. We present a catalogue of 937 point sources detected in the deep \chandra\ observations. We present identifications of our X-ray sources from deep ground-based, Spitzer, GALEX and HST imaging. Using a likelihood ratio analysis, we associate multi band counterparts for 929/937 of our X-ray sources, with an estimated 95~\% reliability, making the identification completeness approximately 94~\% in a statistical sense. Reliable spectroscopic redshifts for 353 of our X-ray sources are provided predominantly from Keck (DEEP2/3) and MMT Hectospec, so the current spectroscopic completeness is $\sim 38$~per cent. For the remainder of the X-ray sources, we compute photometric redshifts based on multi-band photometry in up to 35 bands from the UV to mid-IR. Particular attention is given to the fact that the vast majority the X-ray sources are AGN and require hybrid templates. Our photometric redshifts have mean accuracy of $\sigma=0.04$ and an outlier fraction of approximately 5\%, reaching $\sigma=0.03$ with less than 4\% outliers in the area covered by CANDELS . The X-ray, multi-wavelength photometry and redshift catalogues are made publicly available.

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 all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

This paper describes the mapmaking procedure applied to Planck LFI (Low Frequency Instrument) data. The mapmaking step takes as input the calibrated timelines and pointing information. The main products are sky maps of $I,Q$, and $U$ Stokes components. For the first time, we present polarization maps at LFI frequencies. The mapmaking algorithm is based on a destriping technique, enhanced with a noise prior. The Galactic region is masked to reduce errors arising from bandpass mismatch and high signal gradients. We apply horn-uniform radiometer weights to reduce effects of beam shape mismatch. The algorithm is the same as used for the 2013 release, apart from small changes in parameter settings. We validate the procedure through simulations. Special emphasis is put on the control of systematics, which is particularly important for accurate polarization analysis. We also produce low-resolution versions of the maps, and corresponding noise covariance matrices. These serve as input in later analysis steps and parameter estimation. The noise covariance matrices are validated through noise Monte Carlo simulations. The residual noise in the map products is characterized through analysis of half-ring maps, noise covariance matrices, and simulations.

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.

We measure the star formation rates (SFRs) of massive ($M_{\star}>10^{10.1}M_{\odot}$) early-type galaxies (ETGs) in a sample of 11 high-redshift ($1.0 < z < 1.5$) galaxy clusters drawn from the IRAC Shallow Cluster Survey (ISCS). We identify ETGs visually from Hubble Space Telescope imaging and select likely cluster members as having either an appropriate spectroscopic redshift or red sequence color. Mid-infrared SFRs are measured using Spitzer 24 $\mu$m data for isolated cluster galaxies for which contamination by neighbors, and active galactic nuclei, can be ruled out. Cluster ETGs show enhanced specific star formation rates (sSFRs) compared to cluster galaxies in the local Universe, but have sSFRs more than four times lower than that of field ETGs at $1 < z < 1.5$. Relative to the late-type cluster population, isolated ETGs show substantially quenched mean SFRs, yet still contribute 12% of the overall star formation activity measured in $1 < z < 1.5$ clusters. We find that new ETGs are likely being formed in ISCS clusters; the fraction of cluster galaxies identified as ETGs increases from 34% to 56% from $z \sim 1.5 \rightarrow 1.25$. While the fraction of cluster ETGs that are highly star-forming ($\textrm{SFR}\geq26\ M_{\odot}$ yr$^{-1}$) drops from 27% to 10% over the same period, their sSFRs are roughly constant. All these factors taken together suggest that, particularly at $z\gtrsim1.25$, the events that created these distant cluster ETGs$-$likely mergers, at least among the most massive$-$were both recent and gas-rich.

The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_\ell$. For the dust, they are well described by power laws in $\ell$ with exponents $\alpha^{EE,BB}=-2.42\pm0.02$. The amplitudes of the polarization $C_\ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $\beta_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no "clean" windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $\ell(\ell+1)C_\ell^{BB}/(2\pi)$ of $1.32\times10^{-2}\mu$K$_{CMB}^2$ over the $40<\ell<120$ range; the statistical uncertainty is $\pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $\ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky.

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (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 present 1-2 arcmin spatial resolution CARMA-8 31-GHz observations towards 19 unconfirmed Planck cluster candidates, selected to have significant galaxy overdensities from the WISE early data release and thought to be at z>1 from the WISE colors of the putative brightest cluster galaxy (BCG). We find a Sunyaev-Zeldovich (SZ) detection in the CARMA-8 data towards 9 candidate clusters, where one detection is considered tentative. For each cluster candidate we present CARMA-8 maps, a study of their radio-source environment and we assess the reliability of the SZ detection. The CARMA SZ detections appear to be SZ-bright, with the mean, primary-beam-corrected peak flux density of the decrement being -2.9mJy/beam with a standard deviation of 0.8, and are typically offset from the Planck position by approximately 80 arcsec. Using archival imaging data in the vicinity of the CARMA SZ centroids, we present evidence that one cluster matches Abell 586-a known z~0.2 cluster; four candidate clusters are likely to have 0.3<z<0.7; and, for the remaining 4, the redshift information is inconclusive. We also argue that the sensitivity limits resulting from the cross-correlation between Planck and WISE makes it challenging to use our selection criterion to identify clusters at z > 1.

We present rest-frame far-infrared (FIR) and optical observations of the host galaxy of GRB090423 at z=8.23 from the Atacama Large Millimeter Array (ALMA) and the Spitzer Space Telescope, respectively. The host remains undetected to 3-sigma limits of Fnu(222 GHz)<33 microJy and Fnu(3.6 micron)<81 nJy. The FIR limit is about 20 times fainter than the luminosity of the local ULIRG Arp220, and comparable to the local starburst M82. Comparing to model spectral energy distributions we place a limit on the IR luminosity of L_IR(8-1000 micron)<3e10 Lsun, corresponding to a limit on the obscured star formation rate of SFR_IR<5 Msun/yr; for comparison, the limit on the unobscured star formation rate from Hubble Space Telescope rest-frame UV observations is SFR_UV<1 Msun/yr. We also place a limit on the host galaxy stellar mass of <5e7 Msun (for a stellar population age of 100 Myr and constant star formation rate). Finally, we compare our millimeter observations to those of field galaxies at z>4 (Lyman break galaxies, Ly-alpha emitters, and submillimeter galaxies), and find that our limit on the FIR luminosity is the most constraining to date, although the field galaxies have much larger rest-frame UV/optical luminosities than the host of GRB090423 by virtue of their selection techniques. We conclude that GRB host galaxies at z>4, especially those with measured interstellar medium metallicities from afterglow spectroscopy, are an attractive sample for future ALMA studies of high redshift obscured star formation.

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

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.

This paper presents the large-scale polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse large-scale maps of dust polarization fraction and polarization direction, while taking account of noise bias and possible systematic effects. We find that the maximum observed dust polarization fraction is high (pmax > 18%), in particular in some of the intermediate dust column density (AV < 1mag) regions. There is a systematic decrease in the dust polarization fraction with increasing dust column density, and we interpret the features of this correlation in light of both radiative grain alignment predictions and fluctuations in the magnetic field orientation. We also characterize the spatial structure of the polarization angle using the angle dispersion function and find that, in nearby fields at intermediate latitudes, the polarization angle is ordered over extended areas that are separated by filamentary structures, which appear as interfaces where the magnetic field sky projection rotates abruptly without apparent variations in the dust column density. The polarization fraction is found to be anti-correlated with the dispersion of the polarization angle, implying that the variations are likely due to fluctuations in the 3D magnetic field orientation along the line of sight sampling the diffuse interstellar medium.We also compare the dust emission with the polarized synchrotron emission measured with the Planck LFI, with low-frequency radio data, and with Faraday rotation measurements of extragalactic sources. The two polarized components are globally similar in structure along the plane and notably in the Fan and North Polar Spur regions. A detailed comparison of these three tracers shows, however, that dust and cosmic rays generally sample different parts of the line of sight and confirms that much of the variation observed in the Planck data is due to the 3D structure of the magnetic field.