We present a systematic search for wide-separation (Einstein radius >1.5"), galaxy-scale strong lenses in the 30 000 sq.deg of the Pan-STARRS 3pi survey on the Northern sky. With long time delays of a few days to weeks, such systems are particularly well suited for catching strongly lensed supernovae with spatially-resolved multiple images and open new perspectives on early-phase supernova spectroscopy and cosmography. We produce a set of realistic simulations by painting lensed COSMOS sources on Pan-STARRS image cutouts of lens luminous red galaxies with known redshift and velocity dispersion from SDSS. First of all, we compute the photometry of mock lenses in gri bands and apply a simple catalog-level neural network to identify a sample of 1050207 galaxies with similar colors and magnitudes as the mocks. Secondly, we train a convolutional neural network (CNN) on Pan-STARRS gri image cutouts to classify this sample and obtain sets of 105760 and 12382 lens candidates with scores pCNN>0.5 and >0.9, respectively. Extensive tests show that CNN performances rely heavily on the design of lens simulations and choice of negative examples for training, but little on the network architecture. Finally, we visually inspect all galaxies with pCNN>0.9 to assemble a final set of 330 high-quality newly-discovered lens candidates while recovering 23 published systems. For a subset, SDSS spectroscopy on the lens central regions proves our method correctly identifies lens LRGs at z~0.1-0.7. Five spectra also show robust signatures of high-redshift background sources and Pan-STARRS imaging confirms one of them as a quadruply-imaged red source at z_s = 1.185 strongly lensed by a foreground LRG at z_d = 0.3155. In the future, we expect that the efficient and automated two-step classification method presented in this paper will be applicable to the deeper gri stacks from the LSST with minor adjustments.

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

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

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

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

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

We present spatially resolved ALMA [CII] observations of the bright (flux density S=400 mJy at 350 microns), gravitationally lensed, starburst galaxy PLCK G045.1+61.1 at z=3.427, the "Garnet". This source is part of our set of "Planck's Dusty GEMS", discovered with the Planck's all-sky survey. Two emission-line clouds with a relative velocity offset of ~600 km/s extend towards north-east and south-west, respectively, of a small, intensely star-forming clump with a star-formation intensity of 220 Msun/yr/kpc^2, akin to maximal starbursts. [CII] is also seen in absorption, with a redshift of +350 km/s relative to the brightest CO component. [CII] absorption has previously only been found in the Milky Way along sightlines toward bright high-mass star-forming regions, and this is the first detection in another galaxy. Similar to Galactic environments, the [CII] absorption feature is associated with [CI] emission, implying that this is diffuse gas shielded from the UV radiation of the clump, and likely at large distances from the clump. Since absorption can only be seen in front of a continuum source, the gas in this structure can definitely be attributed to gas flowing towards the clump. The absorber could be part of a cosmic filament or merger debris being accreted onto the galaxy. We discuss our results also in light of the on-going debate of the origin of the [CII] deficit in dusty star-forming galaxies.

(abridged) We report the discovery of PHz G95.5-61.6, a complex structure detected in emission in the Planck all-sky survey that corresponds to two over-densities of high-redshift galaxies. This is the first source from the Planck catalogue of high-z candidates that has been completely characterised with follow-up observations from the optical to the sub-millimetre domain. Herschel/SPIRE observations at 250, 350 and 500 microns reveal the existence of five sources producing a 500 microns emission excess that spatially corresponds to the candidate proto-clusters discovered by Planck. Further observations at CFHT in the optical bands (g and i) and in the near infrared (J, H and K_s), plus mid infrared observations with IRAC/Spitzer (at 3.6 and 4.5 microns) confirm that the sub-mm red excess is associated with an over-density of colour-selected galaxies. Follow-up spectroscopy of 13 galaxies with VLT/X-Shooter establishes the existence of two high-z structures: one at z~1.7 (three confirmed member galaxies), the other at z~2.0 (six confirmed members). This double structure is also seen in the photometric redshift analysis of a sample of 127 galaxies located inside a circular region of 1'-radius containing the five Herschel/SPIRE sources, where we found a double-peaked excess of galaxies at z~1.7 and z~2.0 with respect to the surrounding region. These results suggest that PHz G95.5-61.6 corresponds to two accreting nodes, not physically linked to one another, embedded in the large scale structure of the Universe at z~2 and along the same line-of-sight. In conclusion, the data, methods and results illustrated in this pilot project confirm that Planck data can be used to detect the emission from clustered, dusty star forming galaxies at high-z, and, thus, to pierce through the early growth of cluster-scale structures.

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.

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