Tracing the reionization epoch with ALMA: [CII] emission in z~7 galaxies
L. Pentericci, S.Carniani, M. Castellano, A. Fontana, R. Maiolino, L. Guaita, E. Vanzella, A. Grazian, P. Santini, H. Yan, S. Cristiani, C. Conselice, M. Giavalisco, N. Hathi, A. Koekemoer
DD RAFT VERSION S EPTEMBER
1, 2016
Preprint typeset using L A TEX style emulateapj v. 5/2/11
TRACING THE REIONIZATION EPOCH WITH ALMA: [CII] EMISSION IN Z ∼ L. P
ENTERICCI , S.C ARNIANI , M. C
ASTELLANO , A. F ONTANA , R. M AIOLINO , L. G
UAITA , E. V ANZELLA , A. G RAZIAN , P.S ANTINI , H. Y AN , S. C RISTIANI , C. C ONSELICE , M. G IAVALISCO , N. H ATHI AND
A. K
OEKEMOER Draft version September 1, 2016
ABSTRACTWe present new results on [CII]158 µ m emission from four galaxies in the reionization epoch. These galaxieswere previously confirmed to be at redshifts between 6.6 and 7.15 from the presence of the Ly α emission linein their spectra. The Ly α emission line is redshifted by 100-200 km s − compared to the systemic redshiftgiven by the [CII] line. These velocity offsets are smaller than what is observed in z ∼ α and hence somewhat alleviate the need for a very neutral IGM at z ∼ α emitters observed at this epoch. The galaxies show [CII] emission with L[CII]=0 . − . × L (cid:12) : theseluminosities place them consistently below the SFR-L[CII] relation observed for low redshift star forming andmetal poor galaxies and also below z=5.5 Lyman break galaxies with similar star formation rates. We argue thatprevious undetections of [CII] in z ∼ α emitters and therefore probably metal poor systems, while ourgalaxies are more representative of the general high redshift star forming population. Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: formation INTRODUCTION
We have just entered an exciting era when cosmic mi-crowave background observations can be directly comparedto observations of the first galaxies. Reionization is thoughtto begin at z ∼ −
15 and be completed by z ∼ z > z ∼
10 primar-ily using the Lyman Break technique but their spectroscopicconfirmation has been difficult. At present very few galaxiesat z > INAF, Osservatorio Astronomico di Roma, via Frascati 33, 00040Monteporzio, Italy Kavli Institute for Cosmology, University of Cambridge, MadingleyRoad, Cambridge CB3 0HA, UK Cavendish Laboratory, University of Cambridge, 19 J. J. ThomsonAve., Cambridge CB3 0HE, UK INAF, Osservatorio Astronomico di Bologna, via Ranzani 1, I-40127Bologna, Italy Department of Physics and Astronomy, University of Missouri,Columbia, MO, USA INAF, Osservatorio Astronomico di Trieste, Via G. B. Tiepolo, 11,I-34143 Trieste, Italy University of Nottingham, School of Physics and Astronomy, Not-tingham NG7 2RD, UK Astronomy Department, University of Massachusetts, Amherst, MA01003, USA 0000-0002-7831-8751 Aix Marseille UniversitÃl’, CNRS, LAM (Laboratoired’Astrophysique de Marseille Space Telescope Science Institute, 3700 San Martin Drive, Baltimore,MD 21208, USA verse was not yet complete at z (cid:39)
7. This marked decrease ofLy α emission at z ∼ z ∼ ∆ χ ∼ . ∼ − M (cid:12) yr − . However dust is surely present even if insmall amounts, at such early epochs (Watson et al. 2015) andcan strongly suppress the UV-continuum. Accounting for thismay easily raise the inferred SFR to several tens M (cid:12) yr − .ALMA can play a key role in settling the above is-sues and assess the nature of high redshift galaxies. The[CII]157.74 µ m line, the strongest FIR emission lines in starforming galaxies (accounting for 0.1–1% of their bolometricluminosity), is accessible by ALMA even at the highest red-shifts probed to date. [CII] can provide SFR estimates that arenot biased by dust extinction, although they might depend onthe metallicity. In addition this line allows us to accuratelymeasure the systemic redshift of the galaxies. The recent[CII] line detections of Lyman Break galaxies (LBGs) andLy α emitters (LAEs) at z ∼ z = 6 galaxy A383-5.1.At variance with the above picture, the few observations at-tempted on z ∼ z ∼ z ∼ a r X i v : . [ a s t r o - ph . GA ] A ug constraints on the non-detection of the [CII] line, whereas thethird one showed a detection from a region that is not cen-tered on the galaxy. None of the z ∼ z ∼ . z ∼ z ∼ z ∼
7. Throughout the paper, we adopt a cosmologywith Ω Λ = 0 . Ω M = 0 . H =70 km/s/Mpc. Magnitudesare in the AB system. SAMPLE SELECTION AND OBSERVATIONS
We have recently completed CANDELSz7, an ESO spec-troscopic large program using FORS2 at the VLT. Our goalis to systematically study the Ly α emission in galaxies fromz=5.5 to z=7.2 selected in the CANDELS fields (Grogin et al.2011, Koekemoer et al. 2011). With the addition of previ-ous data we have assembled a sample of >
120 LBGs at z ∼ α line, with EW in some cases as low as 5 Å.From this sample we have selected the seven brightest galax-ies with: (1) precise redshifts between 6.6 and 7.2 from Ly α emission; (2) SFR > − M (cid:12) yr − based on UV emission,assuming the Kennicutt (1998) calibration with no dust cor-rection. While these LBGs all have Ly α emission, the major-ity of them would not be selected by the usual LAE criteriabecause the lines only have modest EW.During Cycle-3 we obtained observations for four of theseven galaxies approved in program 2015.1.01105.S. In Ta-ble 1 we show their optical properties. NTTDF6345 was partof one of our early studies (Pentericci et al. 2011), whilethe other galaxies were confirmed by CANDELSz7. COS-MOS13679 has been independently confirmed by Stark et al.(2016). In Table 1 the coordinates reported are the HST H-band centroids in all cases except for NTTDF6345, which wasdetected with HAWK-I in Y-band (Castellano et al. 2010).For NTTDF6345 we had previously obtained lower S/N datain Cycle-2 (prog. 2013.1.01031.S ).Objects were observed with ALMA in band 6 which hasbaselines between 14.7 and 376.9 meters, and provides a min-imal resolution of 0.9”. The center of the band was set at theredshift corresponding to the peak of the Ly α emission. Sen-sitivity was set to reach ∼ log(L[CII])=7.5 (with L[CII] in L (cid:12) ).Observations were performed in Frequency Division Mode.Out of the four spectral windows, SPW1 was centred on theexpected frequency of the [CII] line in the Upper Side Band.This spectral band was set to a spectral resolution of 10 kms − . SPW0 was located on the continuum next to SPW1 (onthe higher frequency side), while SPW2 and SPW3 were lo-cated in the Lower Side Band to sample the continuum.The ALMA observations were carried out from April 2015 toMarch 2016. The number of antennas ranged from 36 to 46.The precipitable water vapour during the observations rangedbetween 0.82 and 3.0 mm. The phases were centred at thepositions reported in Table 1. The data were reduced withthe Common Astronomy Software Application (CASA) andthe final images were produced using the CLEAN task. Thecontinuum image of each target was extracted using all theline-free channels of the four spectral windows, while spec-tral cubes were generated from the SPW1 dataset. F l u x [ m J y ] COSMOS13679 F l u x [ m J y ] NTTDF6345 F l u x [ m J y ] UDS16291
500 0 500 1000Velocity [km/s]0.40.20.00.20.40.60.81.0 F l u x [ m J y ] COSMOS24108
Figure 1.
S. Spectra of the four galaxies in the region of the [CII] emissionline. The velocities are shown with respect to that inferred from the Ly α which is marked with the dashed vertical lines. Table 1
Galaxies optical and spectroscopic properties
ID RA Dec redshift SFR Ly α EW M UV M SED ref M (cid:12) yr − Å 10 M (cid:12) COSMOS13679 150.0990366 2.3436268 7.1453 23.9 15 -21.46 3.0 1,2NTTDF6345 181.4039006 -7.7561900 6.701 25.0 15 -21.57 – 3UDS16291 34.3561430 -5.1856263 6.6381 15.8 6 -20.97 0.6 1COSMOS24108 150.1972224 2.4786508 6.6294 29 27 -21.67 3.9 1
Note . — Redshift reference: 1 Pentericci et al. in preparation; 2 Stark et al. (2016); 3 Pentericci et al. (2011) RESULTS [CII] line detections
In Table 2 we present our results: a [CII] emission line isdetected for the three sources observed only in Cycle-3, witha S/N=4.5. For NTTDF6345 an emission line is observed sep-arately in the Cycle-2 and Cycle-3 data-sets with a S/N of 4.1and 5.6 respectively, and a S/N>6 is obtained in the combi-nation. Although the S/N is modest, the spatial coincidence(or very close spatial association) between the [CII] and thenear-IR counterparts, and the consistent small shift with re-spect to the redshift determined from the Ly α emission allargue for the reality of the detections. In addition in three outof four cases the detections are spatially resolved. In Figure1 for each source we present the spectrum of the [CII] regionwith a rebinning of 40 km s − for NTTDF6345 and 20 km s − for the other sources: the vertical dotted line at v=0 km s − indicates the redshift determined from the Ly α emission. Themaps of the line emission are shown in Figure 2, extractedwith a spectral width of 440 km s − for NTTDF6345, 280 kms − for COSMOS24108, 100 km s − for UDS16291 and 120km s − for COSMOS13679. The black crosses indicates thecentroid of the Y-band image for NTTDF6345 and of the HSTH-band images for the other galaxies.Thermal far-infrared continuum is not detected in any ofthe galaxies: the limits on the total IR-luminosity convertinto limits on the dust obscured SFR dust that are reported inTable 2 (assuming a Kennicutt (1998) relation with a Salpeterinitial mass function).
Offsets between [CII] and rest-frame UV position
In Table 2 we report the offset between the near-IR coor-dinates and the ALMA detections. In two cases the shiftsare consistent with the ALMA astrometric uncertainty (0.1-0.15 (cid:48)(cid:48) ), while for COSMOS24108 and COSMOS13679 theyare larger. Dunlop et al. (2016) recently noted that the HSTand ALMA astrometry of the HUDF field presented both asystematic shift of 0.25” to the south and a random shift ofup to 0.5 (cid:48)(cid:48) . Spatial offsets of up to 0.5 (cid:48)(cid:48) are clearly evidentin most of the LBGs observed by Capak et al. (2015) andfor the z ∼ ∼ . (cid:48)(cid:48) ), while in the other galaxy CLM1, the [CII] emissionis co-spatial with the UV continuum.To further investigate this issue we looked at brightserendipitous sources detected in the continuum band in ourfields. We find at least one source per field and we measurea shift between 0.1 (cid:48)(cid:48) and 0.6 (cid:48)(cid:48) in random directions betweenthe HST (or HAWK-I) counterpart and the ALMA detection.Given the depth of the CANDELS images ( H ∼
27 at 5 σ ) it isunlikely that the ALMA detections are associated to other ob-jects undetected in the H-band. Note that for NTTDF6345 where we only have HAWK-I images, the shift is negligi-ble both for the serendipitous source and the LBG. We con-clude that there are still substantial uncertainties in the rel-ative ALMA-HST astrometry. In summary, for UDS16291and NTTDF6345 the [CII] detections are centered at the posi-tion of the near-IR sources, and for COSMOS13679 the offset(0.4”) is within the range of those reported in the literature andmeasured for the serendipitous sources: for these three galax-ies we conclude that the [CII] emission comes from the sameregion as the bulk of the SF observed in the near-IR images.For COSMOS24108 the offset is slightly larger than the rangereported in the literature, and while the [CII] emission is al-most certainly associated to the source, it could actually comefrom an external region of the galaxy, not coincident with thebulk of the star formation. A similar case was already ob-served in BDF3299 at z=7.109 (Maiolino et al. 2015): forthis galaxy we concluded that the [CII] emission arises froman external accreting/satellite clump of neutral gas, in agree-ment with recent models of galaxy formation (Vallini et al.2015). While the identification with another transition froma foreground galaxy, which by chance happens to be at 0.8 (cid:48)(cid:48) from COSMOS24108, is unlikely, we cannot completely dis-card the possibility of a spurious detection given the low S/N. DISCUSSION [CII] - Ly α velocity shifts The [CII] line traces the systemic redshift of the sources,unlike the Ly α line which is typically red-shifted by up toseveral hundreds km s − (Erb et al. 2014, Trainor et al. 2016)consistent with the presence of outflowing gas, although thefinal observed Ly α profiles depend on many factors such asgeometry, gas covering fraction, dust and IGM ionizationstate. For our galaxies the velocity shifts are not very large,of the order of 100-200 km s − , smaller than those reported at z ∼ − and av-erage values of 400 km s − for LBGs with M UV < −
21. Themean shift is also lower than those reported by Willott et al.(2015), 430 and 275 km s − respectively for their two z ∼ . The small shifts in our galaxies are particularlysignificant given that in general objects with low Ly α emis-sion have larger velocity offsets. A similar tentative evidencefor smaller Ly α velocities at z ∼ z ∼ λ α compared to the systemic redshift isvery relevant when interpreting the line visibility during thereionization epoch, in the presence of a partly neutral IGM For the Capak et al. (2015) sources this comparison is not possible sincetheir UV redshifts are not from Ly α emission Table 2
Galaxies ALMA properties ID λ [CII]) [CII] flux FWHM [CII] S/N rms cont t int ∆ pos ∆ v M dyn ∗ sin ( i ) SFR dust mm Jy km/s km/s µ Jybeam − s arcsec km s − M (cid:12) M (cid:12) yr − COSMOS13679 1.28426 ± × − ±
35 4.5 14 2782 0.4 135 <2.5 <6.2NTTDF6345 1.2143 ± × − ±
70 6.1 16 2087 0.0 110 <18 <5.7UDS16291 1.20438 ± × − ±
15 4.5 20 2117 0.1 110 <0.7 <6.6COSMOS24108 1.20249 ± × − ±
40 4.5 18 2177 0.8 240 <6.5 <6.2
Figure 2.
Maps of the [CII] emission for COSMOS13679 (upper-left), NTTDF6345 (upper-right), UDS16291(lower-left) and COSMOS24108 (lower-right).Contours are -2,2,3,4 σ level. The crosses indicate the position of the near-IR detections. In the insets we show the near-IR images with ALMA contoursoverlayed. (Dijkstra et al. 2011). Smaller velocity offsets imply that Ly α is closer to resonance and more easily quenched by a neutralIGM. If at redshift 6 and 7 the offsets are as large as thosefound for lower redshift LBGs, the IGM must be very neutralto produce the drop in Ly α fraction that is observed betweenthese two epochs (Pentericci et al. 2014). On the other hand,if the velocity offsets were much smaller, as our observationsindicate, the drop in the Ly α visibility could be produced byan IGM that is still substantially more ionized (Mesinger etal. 2015). SFR-L[CII] relation
In Figure 3 we show the SFR-L[CII] relation for our galax-ies and previously observed sources. We plot COSMOS24108with a different symbol because it is not certain whether its[CII] emission is from the main galaxy or just from a clumpin its outskirt, in which case the point would shift ∼ dex tothe left. We remark that the SFR for our sources as well asprevious z ∼ SFR dust are very low, at least for our galaxies. The SFR for the Capaket al. (2015) sample include both the UV and dust obscuredcontribution. Finally for the Willot et al. sample we plot SED-derived SFR. Our galaxies are a factor of 2-3 less luminous inL[CII] than z ∼ z ∼ α emitters or were LBGs but their spectrashowed a Ly α emission with very high EW (typically >
40 Å).The only exception z8-GND-5296 at z=7.5 which has modestLy α emission, but in this case the L[CII] limit is very shallow(Schaerer et al. 2015). Our four new sources have Ly α emis-sion with low EW (Table 1). The Ly α emission strength isknown to depend on the presence on dust and possibly metal-licity (Raiter et al. 2010): although the derivation of metal-licity is not easy, several studies indicate that LAEs are moremetal poor galaxies compared to the rest of the LBG popula-tion (e.g. Song et al. 2014). Metallicity plays an importantrole in shaping the SFR-L[CII] relation: in Figure 3 we showdifferent metallicity dependent relations produced by a recentstudy of Vallini et al. (2015) based on high-resolution, ra-diative transfer cosmological simulations. The different lineshave Z=0.05,0.1 and 0.2 Z (cid:12) respectively, and for metallic-ity ≤ . z ∼ Dynamical masses
Assuming the sources have ordered motions, we can esti-mate the dynamical masses based on the [CII] velocity dis-persion. We follow the method described in Wang et al.(2013) who approximate M dyn = 1 . × V cir D , where V cir is the circular velocity in km s − , D is the size in kpc, V cir =1 . σ [ CII ] / sin ( i ), and i , the disk inclination angle. Oursources are only marginally resolved so we assume an upperlimit on D from the resolution of 0.9 (cid:48)(cid:48) . We obtain massesbetween 0 . − × M (cid:12) . We can compare these valuesto the total stellar masses determined by a classical SED fit-ting of the multi-wavelength photometry for the 3 CANDELSsources. We use the photometry from the CANDELS cat-alogs, including deep IRAC data essential for sampling therest-frame optical emission at these redshifts. We include thecontribution from nebular emission lines, which can stronglycontaminate the IRAC bands. The stellar masses are reportedin Table 1 and are very similar to the dynamical masses, espe- Figure 3.
S. The blue circles represent the galaxies of the present study withthe empty circle indicating COSMOS24108 for which the association be-tween the [CII] emission and the optical galaxy is uncertain. Red circles areprevious z ∼ ∼ ∼ Z = 0 . Z (cid:12) , 0 . Z (cid:12) and 0 . Z (cid:12) respectively. cially considering the uncertainty due to the unknown incli-nation angle i . CONCLUSIONS
We have presented ALMA [CII] emission line detections infour z ∼ ∼ α emission compared to z ∼ α emission at z ∼
7. These results can help us tune futureALMA observations of high-z LBGs where no Ly α emissionis visible but where precise photometric redshifts exist, as inthe CANDELS fields. Given that in the reionization epochmost galaxies do not show Ly α emission, this means that [CII]can be an efficient alternative to derive the redshift for the ma-jority of the galaxy population. While larger samples of tar-gets with different luminosities and Ly α properties are neededto put our conclusions on firmer grounds, additional diagnos-tic can be used to interpret the physical conditions in theseobjects. For example in case of photo-ionization feedback[NII] at 205 µ m can be much stronger than [CII] (Pavesi et al.2016). Alternatively the [OII] line at 88m µ could be brighterthan [CII] in chemically unevolved systems (Cormier et al.2015) as recently observed in a dust-poor galaxy at z=7.2 (In-oue et al. 2016). REFERENCESCapak, P. L., Carilli, C., Jones, G., et al. 2015, Nature, 522, 455Castellano, M., Fontana, A., Paris, D., et al. 2010, A&A, 524, A28