Planck 's Dusty GEMS

N. Nesvadba, R. Kneissl, R. Cañameras, F. Boone, E. Falgarone, B. Frye, M.

Gerin, S. Koenig, Guilaine Lagache, E. Le Floc'h, et al.

► To cite this version:

N. Nesvadba, R. Kneissl, R. Cañameras, F. Boone, E. Falgarone, et al. Planck 's Dusty GEMS. Astronomy & Astrophysics - A&A, 2016, 593, pp.L2. 10.1051/0004-6361/201629037. hal-02431829

HAL Id: hal-02431829

https://hal.science/hal-02431829v1

Submitted on 13 Nov 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
LETTER TO THE EDITOR

**Planck's Dusty GEMS**

**II. Extended [CII] emission and absorption in the Garnet at z = 3.4**
**seen with ALMA***

N. Nesvadba¹,²,³,⁴, R. Kneissl⁵,⁶, R. Cañamasas¹,²,³,⁴, F. Boone⁷, E. Falgarone⁸, B. Frye⁹, M. Gerin⁸,
S. Koenig¹⁰, G. Lagache¹¹, E. Le Floc'h¹², S. Malhotra¹³, and D. Scott¹⁴

*(Affiliations can be found after the references)*

Received 1 June 2016 / Accepted 20 July 2016

**ABSTRACT**

We present spatially resolved ALMA [CII] observations of the bright (flux density *S*<sub>350</sub> = 400 mJy at 350 µm), 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<sup>-1</sup> extend towards north-east and south-west, respectively, of a small, intensely star-forming clump with a star-formation intensity of 220 M<sub>☉</sub> yr<sup>-1</sup> kpc<sup>-2</sup>, akin to maximal starbursts. [CII] is also seen in absorption, with a redshift of +350 km s<sup>-1</sup> 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.

**Key words.** galaxies: high-redshift – Galaxy: formation – galaxies: ISM – galaxies: starburst – galaxies: kinematics and dynamics – infrared: galaxies

**1. Introduction**

The bright 158 µm line of singly ionized carbon, [CII]158, is one of the most versatile tracers of the interstellar gas in star-forming galaxies. With a low ionization potential, C+ is a probe of the cold neutral gas in galaxies, and can be associated with intensely star-forming environments (e.g., Stacey et al. 2010; Rigopoulou et al. 2014) and diffuse gas (e.g. Langer et al. 2010; Gerin et al. 2015). It is the main coolant of the cold neutral medium (Bennett et al. 1994), and the most luminous line of gas heated by UV photons from star formation over wide ranges of density and UV intensity (e.g., Goldsmith et al. 2012; Le Petit et al. 2006; Kaufman et al. 1999), but can also be bright in shocked gas (Appleton et al. 2013). The diversity of environments and gas conditions probed by [CII] (Rawle et al. 2014; Boone et al. 2014; Schaerer et al. 2015; Knudsen et al. 2016) entail that many empirical properties of [CII] are not yet very well understood, especially at high redshift.

Here we present spatially resolved ALMA cycle 2 observations of [CII] in a strongly gravitationally lensed dusty starburst galaxy at *z* = 3.427, G045.1+61.1 ("the Garnet"), which was discovered using the *Planck* all-sky survey, and was subsequently confirmed with SPIRE imaging (*Planck Collaboration Int. XXVII 2015*), and through ground-based observations (*Cañamasas et al. 2015*, C15 hereafter). The Garnet consists of four counter-images seen behind a small group of galaxies at a spectroscopic redshift of *z* = 0.56 (C15), which together reach a peak flux density of 400 mJy at 350 µm

as seen with SPIRE. The spectral energy distribution (SED) of the dust is consistent with star formation, without obvious signs of AGN contamination (C15). Here we focus on the brightest of the four counter-images, which contributes 46% of the total flux at 850 µm seen with the SMA (C15). The SMA receivers over 90% of the total flux measured with SCUBA-2 at the same wavelength. Our ALMA data resolve the [CII] emission into two clouds at different velocities, around an intensely star-forming clump, against which [CII] is seen in absorption. To our knowledge, no such absorption line has previously been found outside the Milky Way, where it is seen against massive star-forming clouds and is a tracer of the diffuse interstellar medium (ISM, Gerin et al. 2015). Throughout the paper we adopt a flat $H_0 = 70$ km s$^{-1}$ Mpc$^{-1}$ concordance cosmology with $\Omega_M = 0.3$ and $\Omega_\Lambda = 0.7$.

**2. Observations and data reduction**

ALMA observed our Cycle 2 project with 37 antennae on 6 June, 2015 and obtained 27 min of band 8 data in the C34-3 configuration, tuned to the [CII] line (428.775 GHz sky frequency) in one spectral window. Three more windows were placed on the line-free continuum. The intrinsic channel width was 15.6 MHz in all spectral windows, yielding a total bandwidth of 1.856 GHz per window after flagging the edge channels. This corresponds to about 1300 km s⁻¹ at *z* = 3.4, sufficient to cover all [CII] emission, assuming a line profile similar to CO(4–3), which we used as reference. One spectral window had four times higher noise levels on average, which is due to the nearby telluric oxygen line at 424.8 GHz and was ultimately not used in the analysis.

<footnote>* Based on data obtained with ALMA in program 2013.1.01230.S, and with EMIR on the IRAM 30 m telescope in program 223-13.</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Left to right: [CII] channel maps centered on −500 km s⁻¹, systemic velocity, +100 km s⁻¹, and +350 km s⁻¹, respectively, all relative to z = 3.427 and with a channel width of 87 km s⁻¹. Contours show the continuum in the line-free spectral windows and are given for 3, 5, 10, 15, and 20σ. Negative contours are −3σ. The inset shows the SMA 850 µm image with all counter images; the gray circle is the primary beam of ALMA. The other panels show image A. The ALMA beam size in the lower left corner of each map corresponds to 150 pc at z = 3.4 (for µ = 20).</figcaption></figure>

<table><caption>**Table 1.** Emission-line properties for the redshifted and blueshifted EELRs and the two line components seen against the clump.</caption><thead><tr><th>Component</th><th>μ I<sub>CII</sub><br>[Jy km s<sup>-1</sup>]</th><th>Velocity<br>[km s<sup>-1</sup>]</th><th>FWHM<br>[km s<sup>-1</sup>]</th><th>μ L'<sub>CII</sub><br>[10<sup>10</sup> K km s<sup>-1</sup> pc<sup>2</sup>]</th><th>μL<sub>CII</sub><br>[10<sup>9</sup> L<sub>☉</sub>]</th><th>μ</th></tr></thead><tbody><tr><td>Blue EELR</td><td>52.4 ± 0.6</td><td>-467 ± 15</td><td>202 ± 20</td><td>9.5 ± 0.7</td><td>21 ± 0.2</td><td>21</td></tr><tr><td>Red EELR</td><td>8.3 ± 0.6</td><td>97 ± 36</td><td>213 ± 50</td><td>1.5 ± 0.2</td><td>3.3 ± 0.2</td><td>22</td></tr><tr><td rowspan="2">Clump</td><td>4.8 ± 0.6</td><td>-418 ± 47</td><td>190 ± 66</td><td>0.87 ± 0.1</td><td>1.9 ± 0.2</td><td>10</td></tr><tr><td>8.73 ± 0.8</td><td>0 ± 35</td><td>274 ± 49</td><td>1.5 ± 0.1</td><td>3.5 ± 0.3</td><td>10</td></tr><tr><td>Absorber</td><td>-0.60 ± 0.15</td><td>343 ± 34</td><td>124 ± 25</td><td>...</td><td>...</td><td>10</td></tr></tbody></table>

**Notes.** For the absorption line, we list the depth τ of the flux-normalized spectrum instead of the flux.

We used the standard manual scripts for ALMA reduction with the Common Astronomy Software Application CASA, and CLEAN to construct the synthesized beam-deconvolved images of the frequency data cubes. The spectral channels were rebinned to a width of 50 km s⁻¹. The major and minor axis size of the beam is 0.41" × 0.27" along PA = 13.5°. The largest angular scale is 3.9", and the primary beam size 12". The rms is 1.3 mJy in each 50 km s⁻¹ channel, and 0.46 µJy in the averaged continuum image extracted from spectral windows 1 and 2.

In addition, we also used the CO(4–3) line from C15 and the [CI] line corresponding to the ³P₁⁻³P₀ transition of atomic carbon at 492.16 GHz in the rest frame, which we observed during the same campaign with EMIR on the IRAM 30 m telescope with a 20" beam. Our [CI] detection also confirmed the redshift of the Garnet. At z = 3.4266 ± 0.0005, the line falls at (111.183 ± 0.015) GHz. Observations were carried out on 2 March, 2014, under good and stable conditions, for a total observing time of 81 min. Observations and the data reduction are described in C15. The integrated apparent line flux of [CI](1–0) is (8.6 ± 1.6) Jy km s⁻¹, with a peak brightness temperature of 3.3 mK, and FWHM = (420 ± 120) km s⁻¹. The integrated [CI] luminosity is μL' = (3.3 ± 0.7) × 10¹¹ K km s⁻¹ pc², and we estimate a gas mass of (17.4 ± 4) × 10¹¹ μ⁻¹ M☉, a factor of 2 higher than previously estimated in C15 from CO(4–3), and broadly confirming our previous choice of a low CO-to-H₂ conversion factor in C15. All values are uncorrected for the gravitational magnification factor μ. The ³P₂⁻³P₁ [CI](2–1) line at 809 GHz is inaccessible from the ground at z = 3.43. We made this estimate using the relation reported in Papadopoulos et al. (2004) and Wagg et al. (2006), M<sup>[CI]</sup>(H₂) = 1375.8 D²<sub>1+z</sub> X<sub>[CI, -5]</sub>⁻¹ A<sub>10,7</sub>⁻¹ Q<sub>10</sub>⁻¹ S<sub>[CI]</sub>[M☉], where D² is the luminosity distance at redshift z, X<sub>[CI, -5]</sub> = 5 the abundance of atomic carbon in the ISM in units of 10⁻⁵, and A₁₀ the Einstein coefficient in units of 10⁻⁷ s⁻¹. Q₁₀ = 0.7 captures the population of

the upper and lower level of the transition (Papadopoulos et al. 2004)

### 3. Observed morphologies and spectrum

We used the line-free continuum image from spectral windows 1 and 2 to estimate the continuum flux density in each spaxel of the [CII] cube and subtracted this value from every spectral bin. We also corrected for the (insignificant) brightening of the dust continuum by 2% between the line and continuum windows expected for an apparent dust temperature, T<sub>d</sub> = 36 K (C15).

Figure 1 shows the channel maps of the two extended emission line regions (EELRs) as seen in the image plane, which are offset by 564 km s⁻¹ from each other (redshifted and blueshifted channels, respectively), and centered on a small continuum clump with high surface brightness. The integrated spectral properties of each region are listed in Table 1. Each EELR extends about ~2" radially from the continuum source, but the southwestern region is much wider in tangential direction, 2.2" instead of 0.9" measured for the northeastern region (all sizes are 3σ isophotal sizes). This might either represent differences in the intrinsic gas properties or in the gravitational magnification. In either case, the strong sudden velocity offset, which is larger than the full width at half maximum (FWHM) of the individual clouds, shows that the [CII] lines do not probe large-scale rotational motion, but are two kinematically separated clouds.

We used Lenstool (Kneib et al. 1996; Jullo et al. 2007) to construct a gravitational lens model from the positions of the four images seen in Fig. 2 of C15, and calculated the magnification of each image, assuming that all four are multiple images of the same region. The average magnification at the position of the bright star-forming clump is μ = 10, while μ = 21 and μ = 22 for the blueshifted and redshifted gas, respectively. We
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Integrated spectrum in a 1.2" aperture centered on the bright continuum source. The spectrum is shown in black with 87 km s⁻¹ channels, blue lines are individual Gaussian fits, the red line plots the combined fit. The continuum (dashed black line) was measured independently from the same aperture in the line-free continuum spectral windows. The dashed green line shows the integrated [CI] 1–0 profile.</figcaption></figure>

will present a more detailed analysis of the lensing configuration when our scheduled *Hubble* Space Telescope imaging has been taken (GO-14223/PI Frye). The present analysis does not strongly depend on the details of the lens modeling, because the most important physical parameters (surface brightnesses, velocities) are not directly affected by the lensing.

The continuum clump is marginally resolved, with a FWHM size of $0.66 \pm 0.01'' \times 0.51 \pm 0.01''$ along the major and minor axis, respectively, compared to a beam size of $FWHM = 0.41 \times 0.27''$. The continuum flux density of the clump within $2\times$ the FWHM beam size is $42.9 \mu^{-1}$ mJy, corresponding to that expected for a modified blackbody with $\beta = 2.0$ and $T = 36$ K (C15). The clump alone contributes 28% of the total FIR luminosity of G045.1+61.1. Another 18% comes from faint, diffuse emission associated with the two emission-line regions (Fig. 1), the remaining flux comes from the other lensed images shown in the inset of Fig. 1. Adopting a conversion from $L_{FIR}$, integrated between 8 $\mu$m and 1000 $\mu$m, and star formation rate, SFR, of SFR [$M_{\odot} yr^{-1}$] = $4.5 \times 10^{-37}$ W (Kennicutt 1989), and the apparent size of the clump of $0.66'' \times 0.51''$, this implies an average star-formation intensity of $220 M_{\odot} yr^{-1} kpc^{-2}$, which is in the range of maximal starbursts (Elmegreen 1999).

In Fig. 1 we also show the spectrum extracted from an aperture with 1.2" diameter centered on the bright continuum source. As expected from the strong velocity jump in this region (also visible in the channel maps), we see both velocity components, one at +100 km s⁻¹ from the systemic redshift of $z = 3.427$ measured at the [CI] peak and one at -500 km s⁻¹. Line emission is faint near the bright star-forming clump; most flux comes from the extended blueshifted and redshifted EELRs (Fig. 1 and Table 1). In addition, we also see an absorption feature around velocities of +350 km s⁻¹. The depth of the absorption trough was measured from the normalized spectrum to be $-0.60 \pm 0.15$, and other fit parameters are listed in Table 1.

We examined whether this feature might be an observational artifact. The line is resolved into at least five spectral channels and does not depend on the specific choice of aperture size or CLEAN parameters used in the reduction. It is also apparent in

the dirty maps. We do not detect any strongly negative regions in the maps, and our field does not contain bright sources that could cause strong negative signals. The spectra in the continuum spectral windows are flat, and neither the check source nor the phase calibrator show a corresponding amplitude drop. We therefore conclude that there is no evidence that the absorption-line feature in G045.1+61.1 is spurious.

## 4. Astrophysical nature of [CII] in the Garnet

Our ALMA data show a rich environment with two kinematically offset emission-line clouds around a bright star-forming clump, and another component is seen in absorption. For a magnification factor of $\mu \sim 10$, the FIR luminosity of the clump implies $\sim 400 M_{\odot} yr^{-1}$ of star formation, an order of magnitude greater than in massive star clusters in the Milky Way and nearby galaxies ($0.1 - \text{few} \times 10 M_{\odot} yr^{-1}$, Larsen & Richtler 2000).

We used the PDR models of Kaufman et al. (1999) to derive the average gas properties from the [CII] luminosities obtained with ALMA, the CO(4–3) luminosity of C15, and the [CI] 1–0 luminosity (Sect. 2). From the ratios $L_{CII}/L_{CI} = 72.5 \pm 13$ and $L_{CII}/L_{CO} = 23 \pm 1$, we estimate that the gas in the Garnet is on average exposed to a radiation field 200 times greater than in the solar neighborhood, and has a density of about $10^4 \text{ cm}^{-3}$. To derive these line ratios, we corrected the [CI] 1–0 and CO(4–3) luminosities by a factor of 0.46 to take into account that the fainter images are also in the beam, and we assumed that the gas probed by the CO line emission is optically thick.

[CII] emission extends over at least 1200 km s⁻¹ in the Garnet, over a small area of $\le 1$ kpc in the source plane (assuming $\mu \ge 10$ as suggested by our preliminary lensing analysis, see Sect. 3). It is therefore interesting that the Garnet has only redshifted, and no blueshifted [CII] absorption. Unlike EELRs, where blue- and redshifts may either indicate outflows or inflows from the host galaxy (because the gas may be in front of the galaxy or behind it), redshifted absorption is always found in front of the emitter, and therefore a unique signature of an inflow, not an outflow. The [CI] emission line shows a distinct wing at the same velocity as the [CII] absorption. This is also characteristic of [CII] absorption in the Milky Way. Gerin et al. (2015) have shown that both lines can simultaneously arise from gas with similar conditions. By analogy with the Milky Way, we also assume that the line we see consists of multiple deep narrower absorption components that sample a velocity range of 120 km s⁻¹. This might be the intrinsic velocity range of multiple clouds in a filament or merger debris. Even in the Milky Way, where average velocity dispersions in the ISM are about 10 km s⁻¹, the total width over which absorption is found is known to be up to 70–80 km s⁻¹.

We followed Gerin et al. (2015) to estimate a column density of the absorbing gas. With their Eq. (1) originally taken from Goldsmith et al. (2012) to derive a $C^+$ column density, $N(C^+) = 1.4 \times 10^{17} \int \tau dv[\text{cm}^{-2}]$, we find $N(C^+) = 8.2 \times 10^{18} \text{ cm}^{-2}$, corresponding to $N(\text{HI}) = 5.9 \times 10^{22} \text{ cm}^{-2}$ when adopting the Galactic carbon abundance of $1.4 \times 10^{-4}$. This choice is reasonable for massive dusty starburst galaxies at high redshift, which have gas-phase metallicities comparable to those in the Milky Way (Takata et al. 2006; Nesvadba et al. 2007), but might be too high for infall of more pristine gas. For solar metallicity, this column density is about a factor of 2 higher than the most extreme values found in the Milky Way, and is accordingly more suitable for lower metallicity gas. The implied total column density of neutral gas is comparable to those estimated for $H_2$ from
the CO emission-line surface brightness (Swinbank et al. 2011; Cañamas et al., in prep.) and also plausible here, given our gas mass estimate from [CI].

It is impossible to estimate a total mass from an absorption line, therefore we constrained the total mass within this region from the flux in the red wing of [CI], that coincides in velocity with the [CII] absorption (Fig. 1). The wing contains 26% of the [CI] 1–0 flux, and $3.2 \times 10^{10} M_{\odot}$ in mass (for $\mu \sim 10$), comparable to the total molecular gas content of high-z star-bursts (e.g., Tacconi et al. 2008). Carbon is easily ionized (e.g., Goldsmith et al. 2012), and the absorber must therefore be effectively shielded from the intense UV radiation of the clump, and cannot be associated with intense star formation. The absorber is most likely located at a large distance from the clump, potentially several kpc. A velocity offset of 350 km s$^{-1}$ is also too large for the gas to be associated with the clump itself, and is comparable to velocity dispersions of massive low-redshift galaxies (Bernardi et al. 2005), suggesting the gas is bound to the host galaxy. The Garnet provides a rare opportunity to study gas in a high-redshift galaxy outside of intense star-forming regions, and possibly even infalling gas from a cosmic filament, merger debris or satellite accretion.

The blue- and redshifted EELRs extend from a massive star-forming clump, which could be probing a wind. We used the recent empirical analysis of Heckman et al. (2015) of momentum-driven winds in intensely star-forming low-redshift galaxies to show that this scenario is unlikely, and alternative explanations, such as overlapping gas clouds in a star-forming region or merger, are more realistic. Heckman et al. (2015) argued that winds will escape if the momentum input from star formation is $\ge 10\times$ a critical value that depends on the column density, $N_{\text{C}}$, and circular velocity, $v_{\text{c}}$, of the galaxy. We used the 400 km s$^{-1}$ velocity of the absorbing gas to approximate the circular velocity of the Garnet at large distances.

Using Eq. (3) of Heckman et al. (2015), we found $p_{\text{crit}} = 2 \times 10^{36}$ dyn, for a fiducial cloud distance of 100 pc from the star-burst. The critical momentum increases linearly with increasing distance. We also adopted $N_{\text{C}} = 1 \times 10^{23} \text{ cm}^{-2}$, $v_{\text{c}} = 400 \text{ km s}^{-1}$, and used the mass per H$_2$ molecule, that is $3.24 \times 10^{-24} \text{ g}$. $N_{\text{C}}$ is set by the gas mass surface density estimates of 0.4 and $3.9 \times 10^9 M_{\odot} \text{ kpc}^{-2}$ implied by the [CI] line profile, assuming that [CI] and [CII] have the same morphology. Assuming (again following Heckman et al. 2015) that the combined momentum flux from ram and radiation pressure in the star-forming clump is $4.8 \times 10^{33} \text{ SFR dyn}$, we find a total momentum input from star formation of $1.9 \times 10^{36}$ dyn for SFR = 400 $M_{\odot} \text{ yr}^{-1}$. Following Heckman et al. (2015), this may be enough to balance gravity, but not to form a wind that escapes, the main difference to low-redshift galaxies being the high column densities. Of course, this does not rule out the presence of fainter more energetic wind components, as reported for unlensed galaxies (Maiolino et al. 2012; Cicone et al. 2015).

## 5. Implications for high-z galaxies and conclusions

It is also interesting to investigate how our sources relate to dusty star-forming galaxies at high redshifts in the field, where we have more comprehensive knowledge of the global properties, but lack the detailed spatially resolved information that only strong lensing can provide. Infrared-selected, intensely star-forming galaxies at low redshift show a marked deficit in [CII] luminosity relative to the total FIR luminosity, compared to galaxies with lower star formation rates and less extinction; this is

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Ratio of [CII] to FIR luminosity as a function of FIR luminosity for the Garnet and for several samples of low and high-redshift galaxies in the literature. The FIR luminosities have been corrected for average magnification factors of 10, 21, and 21 at the position of the total source, and the blueshifted gas and the continuum source, respectively.</figcaption></figure>

a trend that is rarely found at high redshift (Rigopoulou et al. 2014; Malhotra et al. 2016). We show in Fig. 3 where the Garnet lies relative to other populations of low- and high-redshift galaxies. The absorption does not lower the continuum luminosity, which is derived from the integrated flux between 8 and $1000 \ \mu\text{m}$ in the rest-frame. The blue EELR falls into the region of the diagram, in which distant and nearby star-forming galaxies, and local ULIRGs overlap, whereas the aperture centered on the clump lies near the extreme end of the local ULIRG sample. This highlights the fact that unresolved observations of high-redshift galaxies sample wide ranges of intrinsic line ratios. The integrated spectrum of the Garnet lies very close to the blue EELR, suggesting that the global [CII]/FIR ratios are dominated by the global spatial distribution and sizes of gas clouds within star-forming regions, and do not reflect the properties of individual star clusters.

These ALMA data of G045.1+61.1 at $z = 3.43$, the first of Planck's Dusty GEMS, shows that [CII] at high redshift probes a wide range of environments, from the dense gas irradiated by UV photons from maximally star-forming clumps to the diffuse interstellar gas at large distances seen in absorption. Although we found signatures of [CII] spanning at least $1200 \text{ km s}^{-1}$, we did not detect strong evidence for a wind that would regulate star formation by removing the gas reservoirs, as is often suggested in starburst-driven wind scenarios. This might indicate that broad wings of [CII] and other lines in high-redshift galaxies probe a much richer phenomenology than previously thought. Deeper higher-resolution observations spanning a wider spectral range will open an interesting new window to study the cold neutral gas in intense high-redshift starburst galaxies, including the ambient gas outside of the intense star-forming regions.

**Acknowledgements.** We thank the ALMA staff for carrying out the observations, and the team at the ARC in Grenoble for their help with preparing the data taking. We also thank the referee for suggestions that helped improve the paper. NPHN wishes to acknowledge funding through the JAO within their visiting scientists program, and is grateful for their hospitality during her stay as a science visitor in Vitacura.

## References

* Appleton, P. N., Guillard, P., Boulanger, F., et al. 2013, ApJ, 777, 66

* Bennett, C. L., Fixsen, D. J., Hinshaw, G., et al. 1994, ApJ, 434, 587

* Bernardi, M., Sheth, R. K., Nichol, R. C., Schneider, D. P., & Brinkmann, J. 2005, AJ, 129, 61
Boone, F., Combes, F., Krips, M., et al. 2014, in SF2A-2014: Proc. of the Annual meeting of the French Society of Astronomy and Astrophysics, eds. J. Ballet, F. Martins, F. Bournaud, R. Monier, & C. Reylé, 381

Cañameras, R., Nesvadba, N. P. H., Guery, D., et al. 2015, A&A, 581, A105

Cicone, C., Maiolino, R., Gallerani, S., et al. 2015, A&A, 574, A14

Elmegreen, B. G. 1999, ApJ, 517, 103

Gerin, M., Ruaud, M., Goicoechea, J. R., et al. 2015, A&A, 573, A30

Goldsmith, P. F., Langer, W. D., Pineda, J. L., & Velusamy, T. 2012, ApJS, 203, 13

Heckman, T. M., Alexandroff, R. M., Borthakur, S., Overzier, R., & Leitherer, C. 2015, ApJ, 809, 147

Jullo, E., Kneib, J.-P., Limousin, M., et al. 2007, New J. Phys., 9, 447

Kaufman, M. J., Wolfire, M. G., Hollenbach, D. J., & Luhman, M. L. 1999, ApJ, 527, 795

Kennicutt, Jr., R. C. 1989, ApJ, 344, 685

Kneib, J.-P., Ellis, R. S., Smail, I., Couch, W. J., & Sharples, R. M. 1996, ApJ, 471, 643

Knudsen, K. K., Richard, J., Kneib, J.-P., et al. 2016, MNRAS, 462, L6

Langer, W. D., Velusamy, T., Pineda, J. L., et al. 2010, A&A, 521, L17

Larsen, S. S., & Richtler, T. 2000, A&A, 354, 836

Le Petit, F., Nehmé, C., Le Bourlot, J., & Roueff, E. 2006, ApJS, 164, 506

Maiolino, R., Gallerani, S., Neri, R., et al. 2012, MNRAS, 425, L66

Nesvadba, N. P. H., Lehnert, M. D., Genzel, R., et al. 2007, ApJ, 657, 725

Papadopoulos, P. P., Thi, W.-F., & Viti, S. 2004, MNRAS, 351, 147

Planck Collaboration Int. XXVII. 2015, A&A, 582, A30

Rawle, T. D., Egami, E., Bussmann, R. S., et al. 2014, ApJ, 783, 59

Rigopoulou, D., Hopwood, R., Magdis, G. E., et al. 2014, ApJ, 781, L15

Schäerer, D., Boone, F., Jones, T., et al. 2015, A&A, 576, L2

Stacey, G. J., Hailey-Dunsheath, S., Ferkinhoff, C., et al. 2010, ApJ, 724, 957

Swinbank, A. M., Papadopoulos, P. P., Cox, P., et al. 2011, ApJ, 742, 11

Tacconi, L. J., Genzel, R., Smail, I., et al. 2008, ApJ, 680, 246

Takata, T., Sekiguchi, K., Smail, I., et al. 2006, ApJ, 651, 713

Wagg, J., Wilner, D. J., Neri, R., Downes, D., & Wiklind, T. 2006, ApJ, 651, 46

¹ Institut d'Astrophysique Spatiale, Bât. 121, 91405 Orsay Cedex, France
e-mail: nicole.nesvadba@ias.u-psud.fr

² CNRS, France

³ Université Paris-Sud, France

⁴ Université Paris-Saclay, France

5 European Southern Observatory, ESO Vitacura, Alonso de Cordova 3107, Vitacura, Casilla 19001 Santiago, Chile

6 Atacama Large Millimeter/submillimeter Array, ALMA Santiago Central Offices, Alonso de Cordova 3107, Vitacura, Casilla 763-0355 Santiago, Chile

7 Université de Toulouse, UMS-OMP, IRAP, 31028 Toulouse Cedex 4, France

8 LERMA, UMR 8112 CNRS, École Normale Supérieure and Observatoire de Paris, 75000 Paris, France

9 Steward Observatory, University of Arizona, Tucson, AZ 85721, USA

10 Chalmers University of Technology, Onsala Space Observatory, 439 92 Onsala, Sweden

11 Aix Marseille Université, CNRS, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, 13388 Marseille, France

12 CEA-Saclay, 91191 Gif-sur-Yvette, France

13 School of Earth and Space Exploration, Arizona State University, Tempe, AZ 85287, USA

14 Department of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, 6658 British Columbia, Canada