Search

The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous QSOs whose SMBH must have had the fastest mass growth during the Universe first Gyr. In this paper we present the HYPERION sample and report on the first of the 3 years planned observations of the 2.4 Ms XMM-Newton Multi-Year Heritage program on which HYPERION is based. The goal of this program is to accurately characterize the X-ray nuclear properties of QSOs at the EoR. Through a joint X-ray spectral analysis of 10 sources, in the rest-frame $\sim2-50$ keV range, we report a steep average photon index ($\Gamma\sim2.4\pm0.1$). Absorption is not required. The average $\Gamma$ is inconsistent at $\geq4\sigma$ level with the canonical 1.8-2 value measured in QSO at z, Comment: 21 pages (including appendix), 12 figures, 4 tables. Accepted for pubblication in A&A

The outstanding mass growth of supermassive black holes (SMBHs) at the Reionisation Epoch and how it is related to the concurrent growth of their host galaxies, poses challenges to theoretical models aimed at explaining how these systems formed in short timescales (\dot M_{\rm BH}/M_{\rm BH}$, suggesting that AGN feedback might be efficiently acting to slow down the SMBH accretion, while the stellar mass assembly is still vigorously taking place in the host galaxy. In addition, we were also able to detect high-velocity emission on the red and blue sides of the [CII] emission line, that traces a gaseous outflow, and for the first time, we mapped a spatially-resolved water vapour disk through the H$_2$O v=0 $3_{(2,2)}-3_{(1,3)}$ emission line detected at $ν_{\rm obs} = 274.074$ GHz, whose kinematic properties and size are broadly consistent with those of the [CII] disk., 15 pages, 10 figures, 4 tables. Accepted in A&A

The population dynamics that followed the out of Africa expansion (OoA) and ultimately led to the formation of Oceanian, West and East Eurasian macro populations have long been debated 1–4. Furthermore, with the OoA being dated between 70 kya 1,4,5 and 65 kya 6 and the earliest splits between West and East Eurasian populations being inferred not earlier than 43 kya from modern DNA data 1,4,7, an additional question concerns the whereabouts of the early migrants out of Africa before those differentiations. Shedding light on these population dynamics may, in turn, provide clues to better understand cultural evolution in Eurasia between 50 kya and 35 kya, where the development of new technologies may be correlated to parallel independent evolution paths, to the arrival of new populations, or to long-term processes of cultural and biological exchanges. Here we jointly re-analyze Eurasian Paleolithic DNA available to date in light of material culture, and provide a comprehensive population model with minimal admixture events. Our integrated approach i) maintains Zlatý Kůň genetically as the most basal out of Africa human lineage sequenced to date, also in comparison to Oceanians and putatively links it with non-Mousterian material cultures documented in Europe 48-43 kya; ii) infers the presence of an OoA population Hub from which a major wave broadly associated with Initial Upper Paleolithic lithic industries emanated to populate West and East Eurasia before or around 45 kya, and of which Ust’Ishim, Bacho Kiro and Tianyuan were unadmixed descendants; iii) proposes a parsimonious placement of Oase1 as an individual related to Bacho Kiro who experienced additional Neanderthal introgression; and iv) explains the East/West Eurasian population split as a longer permanence of the latter in the OoA Hub, followed by a second population expansion (before 37 kya), broadly associated with Upper Paleolithic industries, that largely replaced pre-existing humans in Europe, and admixed with the previous wave to form Yana and Mal’ta in Siberia and, to a greater extent, GoyetQ116-1 in Belgium.

We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C+ at Early Times (ALPINE) survey - the first large multi-wavelength survey which aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 < z < 5.9 via the measurement of [CII]-emission at 158 micro-meter (64% at >3.5$\sigma$) and the surrounding far-infrared (FIR) continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes ($\beta$), and absorption line redshifts, as well as H-alpha emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased towards bluer colors ($\Delta\beta$ ~ 0.2). Using [CII] as tracer of the systemic redshift (confirmed for one galaxy at z=4.5 out of 118 for which we obtained optical [OII]$\lambda$3727A emission), we confirm red shifted Ly-alpha emission and blue shifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [CII] rest-frame we find that the absorption lines in galaxies with high specific SFR are more blue shifted, which could be indicative of stronger winds and outflows., Comment: 48 pages, 33 figures, 7 tables (including Appendix), accepted for publication in ApJS

The ALMA-ALPINE [CII] survey is aimed at characterizing the properties of a sample of normal star-forming galaxies (SFGs). The ALMA Large Program to INvestigate (ALPINE) features 118 galaxies observed in the [CII]-158 μm line and far infrared (FIR) continuum emission during the period of rapid mass assembly, right after the end of the HI reionization, at redshifts of 4 < z < 6. We present the survey science goals, the observational strategy, and the sample selection of the 118 galaxies observed with ALMA, with an average beam minor axis of about 0.85″, or ∼5 kpc at the median redshift of the survey. The properties of the sample are described, including spectroscopic redshifts derived from the UV-rest frame, stellar masses, and star-formation rates obtained from a spectral energy distribution (SED) fitting. The observed properties derived from the ALMA data are presented and discussed in terms of the overall detection rate in [CII] and FIR continuum, with the observed signal-to-noise distribution. The sample is representative of the SFG population in the main sequence at these redshifts. The overall detection rate in [CII] is 64% for a signal-to-noise ratio (S/N) threshold larger than 3.5 corresponding to a 95% purity (40% detection rate for S/N > 5). Based on a visual inspection of the [CII] data cubes together with the large wealth of ancillary data, we find a surprisingly wide range of galaxy types, including 40% that are mergers, 20% extended and dispersion-dominated, 13% compact, and 11% rotating discs, with the remaining 16% too faint to be classified. This diversity indicates that a wide array of physical processes must be at work at this epoch, first and foremost, those of galaxy mergers. This paper sets a reference sample for the gas distribution in normal SFGs at 4 < z < 6, a key epoch in galaxy assembly, which is ideally suited for studies with future facilities, such as the James Webb Space Telescope (JWST) and the Extremely Large Telescopes (ELTs).

At redshift $z>5$ the far-infrared (FIR) continuum spectra of main-sequence galaxies are sparsely sampled, often with a single data point. The dust temperature $T_{\rm d, SED}$ thus has to be assumed in the FIR continuum fitting. This introduces large uncertainties regarding the derived dust mass ($M_{\rm d}$), FIR luminosity, and obscured fraction of the star formation rate. These are crucial quantities to quantify the effect of dust obscuration in high-$z$ galaxies. To overcome observations limitations, we introduce a new method that combines dust continuum information with the overlying $\hbox{[C $\scriptstyle\rm II $]} 158\mu$m line emission. By breaking the $M_{\rm d} - T_{\rm d, SED}$ degeneracy, with our method, we can reliably constrain the dust temperature with a single observation at $158\mu$m. This method can be applied to all ALMA and NOEMA $\hbox{[C $\scriptstyle\rm II $]}$ observations and exploited in ALMA Large Programs such as ALPINE and REBELS targeting $\hbox{[C $\scriptstyle\rm II $]}$ emitters at high-$z$. We also provide a physical interpretation of the empirical relation recently found between $molecular$ gas mass and $\hbox{[C $\scriptstyle\rm II $]}$ luminosity. We derive an analogous relation linking the $total$ gas surface density and $\hbox{[C $\scriptstyle\rm II $]}$ surface brightness. By combining the two, we predict the cosmic evolution of the surface density ratio $\Sigma_{\rm H_2} / \Sigma_{\rm gas}$. We find that $\Sigma_{\rm H_2} / \Sigma_{\rm gas}$ slowly increases with redshift, which is compatible with current observations at $0 < z < 4$., Comment: 14 pages, 5 figures, submitted to MNRAS

We present the first attempt to detect outflows from galaxies approaching the Epoch of Reionization (EoR) using a sample of 9 star-forming ($\rm SFR=31\pm 20~M_{\odot}~yr^{-1}$) $z\sim 5.5$ galaxies for which the [CII]158$\mu$m line has been previously obtained with ALMA. We first fit each line with a Gaussian function and compute the residuals by subtracting the best fitting model from the data. We combine the residuals of all sample galaxies and find that the total signal is characterised by a flux excess of $\sim 0.5$ mJy extended over $\sim 1000$ km~s$^{-1}$. Although we cannot exclude that part of this signal is due to emission from faint satellite galaxies, we show that the most probable explanation for the detected flux excess is the presence of broad wings in the [CII] lines, signatures of starburst-driven outflows. We infer an average outflow rate of $\rm \dot{M}=54\pm23~ M_{\odot}~yr^{-1}$, providing a loading factor $\eta=\rm \dot{M}/SFR=1.7\pm1.3$ in agreement with observed local starbursts. Our interpretation is consistent with outcomes from zoomed hydro-simulations of {\it Dahlia}, a $z\sim 6$ galaxy ($\rm SFR\sim 100~\rm M_{\odot}~yr^{-1}$) whose feedback-regulated star formation results into an outflow rate $\rm \dot{M}\sim 30~ M_{\odot}~yr^{-1}$. The quality of the ALMA data is not sufficient for a detailed analysis of the [CII] line profile in individual galaxies. Nevertheless, our results suggest that starburst-driven outflows are in place in the EoR and provide useful indications for future ALMA campaigns. Deeper observations of the [CII] line in this sample are required to better characterise feedback at high-$z$ and to understand the role of outflows in shaping early galaxy formation.

IR spectroscopy in the range 12–230 μm with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes governing the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope and the upcoming Extremely Large Telescopes at shorter wavelengths and the Atacama Large Millimeter Array at longer wavelengths. The SPICA, with its 2.5-m telescope actively cooled to below 8 K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields, and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei, sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. Active galactic nuclei and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. The SPICA’s large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z ~ 6.

We present zoom-in, AMR, high-resolution ($\simeq 30$ pc) simulations of high-redshift ($z \simeq 6$) galaxies with the aim of characterizing their internal properties and interstellar medium. Among other features, we adopt a star formation model based on a physically-sound molecular hydrogen prescription, and introduce a novel scheme for supernova feedback, stellar winds and dust-mediated radiation pressure. In the zoom-in simulation the target halo hosts "Dahlia", a galaxy with a stellar mass $M_*=1.6\times 10^{10}$M$_\odot$, representative of a typical $z\sim 6$ Lyman Break Galaxy. Dahlia has a total H2 mass of $10^{8.5}$M$_\odot$, that is mainly concentrated in a disk-like structure of effective radius $\simeq 0.6$ kpc and scale height $\simeq 200$ pc. Frequent mergers drive fresh gas towards the center of the disk, sustaining a star formation rate per unit area of $\simeq 15 $M$_\odot$ yr$^{-1}$ kpc$^{-2}$. The disk is composed by dense ($n \gtrsim 25$ cm$^{-3}$), metal-rich ($Z \simeq 0.5 $ Z$_\odot$) gas, that is pressure-supported by radiation. We compute the $158\mu$m [CII] emission arising from {Dahlia}, and find that $\simeq 95\%$ of the total [CII] luminosity ($L_{[CII]}\simeq10^{7.5}$ L$_\odot$) arises from the H2 disk. Although $30\%$ of the CII mass is transported out of the disk by outflows, such gas negligibly contributes to [CII] emission, due to its low density ($n \lesssim 10$ cm$^{-3}$) and metallicity ($Z\lesssim 10^{-1}$Z$_\odot$). Dahlia is under-luminous with respect to the local [CII]-SFR relation; however, its luminosity is consistent with upper limits derived for most $z\sim6$ galaxies., Comment: 21 pages, 11 figures, 3 tables

We report the detection of [CII]158um emission from a system of three closely-separated sources in the COSMOS field at z~4.56, as part of the ALMA Large Program to INvestigate CII at Early times (ALPINE). The two dominant sources are closely associated, both spatially (1.6"~11kpc) and in velocity (~100km/s), while the third source is slightly more distant (2.8"~18kpc, ~300km/s). The second strongest source features a slight velocity gradient, while no significant velocity gradient is seen in the other two sources. Using the observed [CII] luminosities, we derive a total log(SFR_[CII]/[Msol/year])=2.8+/-0.2, which may be split into contributions of 59%, 31%, and 10% from the central, east, and west sources, respectively. Comparison of these [CII] detections to recent zoom-in cosmological simulations suggests an ongoing major merger. We are thus witnessing a system in a major phase of mass build-up by merging, including an on-going major merger and an upcoming minor merger, which is expected to end up in a single massive galaxy by z~2.5., ERC Advanced Grant 695671 “QUENCH"

We present new ALMA observations of the [O iii]88 μm line and high angular resolution observations of the [C ii]158 μm line in a normal star forming galaxy at z = 7.1. Previous [C ii] observations of this galaxy had detected [C ii] emission consistent with the Lyα redshift but spatially slightly offset relative to the optical (UV-rest frame) emission. The new [C ii] observations reveal that the [C ii] emission is partly clumpy and partly diffuse on scales larger than about 1 kpc. [O iii] emission is also detected at high significance, offset relative to the optical counterpart in the same direction as the [C ii] clumps, but mostly not overlapping with the bulk of the [C ii] emission. The offset between different emission components (optical/UV and different far-IR tracers) is similar to that which is observed in much more powerful starbursts at high redshift. We show that the [O iii] emitting clump cannot be explained in terms of diffuse gas excited by the UV radiation emitted by the optical galaxy, but it requires excitation by in-situ (slightly dust obscured) star formation, at a rate of about 7 M⊙ yr-1. Within 20 kpc from the optical galaxy the ALMA data reveal two additional [O iii] emitting systems, which must be star forming companions. We discuss that the complex properties revealed by ALMA in the z ~ 7.1 galaxy are consistent with expectations by recent models and cosmological simulations, in which differential dust extinction, differential excitation and different metal enrichment levels, associated with different subsystems assembling a galaxy, are responsible for the various appearance of the system when observed with distinct tracers.

We present Atacama Large Millimiter/submillimiter Array (ALMA) observations of eight highly excited CO (J$_{\rm up}>8$) lines and continuum emission in two $z\sim6$ quasars: SDSS J231038.88+185519.7 (hereafter J2310), for which CO(8-7), CO(9-8), and CO(17-16) lines have been observed, and ULAS J131911.29+095951.4 (J1319), observed in the CO(14-13), CO(17-16) and CO(19-18) lines. The continuum emission of both quasars arises from a compact region ($< 0.9$ kpc). By assuming a modified black-body law, we estimate dust masses of Log$(M_{\rm dust}/M_{\odot})=8.75\pm0.07$ and Log$(M_{\rm dust}/M_{\odot})=8.8\pm0.2$ and dust temperatures of $T_{\rm dust}=76\pm3~{\rm K}$ and $T_{\rm dust}=66^{+15}_{-10}~{\rm K}$, respectively for J2310 and J1319. Only CO(8-7) and CO(9-8) in J2310 are detected, while $3\sigma$ upper limits on luminosities are reported for the other lines of both quasars. The CO line luminosities and upper limits measured in J2310 and J1319 are consistent with those observed in local AGN and starburst galaxies, and other $z\sim 6$ quasars, except for SDSS J1148+5251 (J1148), the only quasar at $z=6.4$ with a previous CO(17-16) line detection. By computing the CO SLEDs normalised to the CO(6-5) line and FIR luminosities for J2310, J1319, and J1149, we conclude that different gas heating mechanisms (X-ray radiation and/or shocks) may explain the different CO luminosities observed in these $z\sim6$ quasar. Future J$_{\rm up}>8$ CO observations will be crucial to understand the processes responsible for molecular gas excitation in luminous high-$z$ quasars.

We report the serendipitous discovery of a bright galaxy (Gal-A) observed as part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). While this galaxy is detected both in line and continuum emission in ALMA Band 7, it is completely dark in UV/optical filters and only presents a marginal detection in the UltraVISTA Ks band. We discuss the nature of the observed ALMA line, i.e. whether the emission comes from [CII] at z~4.6, or from high-J CO transitions at z~2.2. In the first case we find a [CII]-to-FIR luminosity ratio of log(L_[CII]/L_FIR)=-2.5, consistent with the average value for local star-forming galaxies (SFGs); in the second case, instead, the source would lie outside of the empirical relations between L_CO and L_FIR found in the literature. At both redshifts, we derive the star-formation rate (SFR) from the ALMA continuum, and the stellar mass (M*) by using stellar population synthesis models as input for LePHARE spectral energy distribution (SED) fitting. Exploiting our results, we believe that Gal-A is a Main-Sequence (MS), dusty SFG at z=4.6 (i.e. [CII] emitter) with log(SFR/[M/yr])~1.4 and log(M*/M)~9.7. This work underlines the crucial role of the ALPINE survey in making a census of this class of objects, in order to unveil their contribution to the global star-formation rate density (SFRD) of the Universe at the end of the Reionisation epoch., Comment: 13 pages, 6 figures. Submitted to MNRAS

At redshift z > 5, the far-infrared (FIR) continuum spectra of main-sequence galaxies are sparsely sampled, often with a single data point. The dust temperature Td,SED, thus has to be assumed in the FIR continuum fitting. This introduces large uncertainties regarding the derived dust mass (Md), FIR luminosity, and obscured fraction of the star formation rate. These are crucial quantities to quantify the effect of dust obscuration in high-z galaxies. To overcome observation limitations, we introduce a new method that combines dust continuum information with the overlying [C $\scriptstyle \rm II$] 158 µm line emission. By breaking the Md–Td,SED degeneracy, with our method, we can reliably constrain the dust temperature with a single observation at 158 µm. This method can be applied to all Atacama Large Millimeter Array (ALMA) and NOEMA [C $\scriptstyle \rm II$] observations, and exploited in ALMA Large Programs such as ALPINE and REBELS targeting [C $\scriptstyle \rm II$] emitters at high-z. We also provide a physical interpretation of the empirical relation recently found between molecular gas mass and [C $\scriptstyle \rm II$] luminosity. We derive an analogous relation linking the total gas surface density and [C $\scriptstyle \rm II$] surface brightness. By combining the two, we predict the cosmic evolution of the surface density ratio $\Sigma _{\rm H_2} / \Sigma _{\rm gas}$. We find that $\Sigma _{\rm H_2} / \Sigma _{\rm gas}$ slowly increases with redshift, which is compatible with current observations at 0 < z < 4.

The intensity mapping of the [CII] 158um line redshifted to the sub-mm window is a promising probe of the z>4 star formation and its spatial distribution into the large-scale structure. To prepare the first-generation experiments (e.g., CONCERTO), we need realistic simulations of the sub-mm extragalactic sky in spectroscopy. We present a new version of the SIDES simulation including the main sub-mm lines around 1mm (CO, [CII], [CI]). This approach successfully reproduces the observed line luminosity functions. We then use our simulation to generate CONCERTO-like cubes (125-305GHz) and forecast the power spectra of the fluctuations caused by the various astrophysical components at those frequencies. Depending on our assumptions on the relation between star formation rate and [CII] luminosity, and the star formation history, our predictions of the z~6 [CII] power spectrum vary by two orders of magnitude. This highlights how uncertain the predictions are and how important future measurements will be to improve our understanding of this early epoch. SIDES can reproduce the CO shot noise recently measured at ~100 GHz by the mmIME experiment. Finally, we compare the contribution of the different astrophysical components at various redshift to the power spectra. The continuum is by far the brightest, by a factor of 3 to 100 depending on the frequency. At 300GHz, the CO foreground power spectrum is higher than the [CII] one for our base scenario. At lower frequency, the contrast between [CII] and extragalactic foregrounds is even worse. Masking the known galaxies from deep surveys should allow to reduce the foregrounds to 20% of the [CII] power spectrum up to z~6.5. However, this masking method will not be sufficient at higher redshifts. The code and the products of our simulation are released publicly and can be used for both intensity mapping experiments and sub-mm continuum and line surveys., 28 pages, 19 figures, 2 tables, accepted by A&A, the code and the products are available at https://cesamsi.lam.fr/instance/sides/home

We study the photoevaporation of Jeans-unstable molecular clumps by isotropic FUV (6 eV $< {\rm h}��$ < 13.6 eV) radiation, through 3D radiative transfer hydrodynamical simulations implementing a non-equilibrium chemical network that includes the formation and dissociation of H$_2$. We run a set of simulations considering different clump masses ($M=10-200$ M$_\odot$) and impinging fluxes ($G_0=2\times 10^3-8\times 10^4$ in Habing units). In the initial phase, the radiation sweeps the clump as an R-type dissociation front, reducing the H$_2$ mass by a factor 40-90%. Then, a weak ($\mathcal{M}\simeq 2$) shock develops and travels towards the centre of the clump, which collapses while loosing mass from its surface. All considered clumps remain gravitationally unstable even if radiation rips off most of the clump mass, showing that external FUV radiation is not able to stop clump collapse. However, the FUV intensity regulates the final H$_2$ mass available for star formation: for example, for $G_0 < 10^4$ more than 10% of the initial clump mass survives. Finally, for massive clumps ($\sim 100$ M$_\odot$) the H$_2$ mass increases by 25-50% during the collapse, mostly because of the rapid density growth that implies a more efficient H$_2$ self-shielding., 17 pages, 16 figures; added references, corrected typos

We present new estimates of AGN accretion and star-formation luminosity in galaxies obtained for the local 12-$\mu$m sample of Seyfert galaxies (12MGS), by performing a detailed broad-band spectral energy distribution (SED) decomposition including the emission of stars, dust heated by star formation and a possible AGN dusty torus. Thanks to the availability of data from the X-rays to the sub-millimetre, we constrain and test the contribution of the stellar, AGN and star-formation components to the SEDs. The availability of Spitzer-IRS low resolution mid-infrared (mid-IR) spectra is crucial to constrain the dusty torus component at its peak wavelengths. The results of SED-fitting are also tested against the available information in other bands: the reconstructed AGN bolometric luminosity is compared to those derived from X-rays and from the high excitation IR lines tracing AGN activity like [Ne V] and [O IV]. The IR luminosity due to star-formation (SF) and the intrinsic AGN bolometric luminosity are shown to be strongly related to the IR line luminosity. Variations of these relations with different AGN fractions are investigated, showing that the relation dispersions are mainly due to different AGN relative contribution within the galaxy. Extrapolating these local relations between line and SF or AGN luminosities to higher redshifts, by means of recent Herschel galaxy evolution results, we then obtain mid- and far-IR line luminosity functions useful to estimate how many star-forming galaxies and AGN we expect to detect in the different lines at different redshifts and luminosities with future IR facilities (e.g., JWST, SPICA)., Comment: 27 pages, 13 figures, 7 tables. MNRAS in press

We present the first attempt to detect outflows from galaxies approaching the Epoch of Reionization (EoR) using a sample of 9 star-forming ($\rm SFR=31\pm 20~M_{\odot}~yr^{-1}$) $z\sim 5.5$ galaxies for which the [CII]158$\mu$m line has been previously obtained with ALMA. We first fit each line with a Gaussian function and compute the residuals by subtracting the best fitting model from the data. We combine the residuals of all sample galaxies and find that the total signal is characterised by a flux excess of $\sim 0.5$ mJy extended over $\sim 1000$ km~s$^{-1}$. Although we cannot exclude that part of this signal is due to emission from faint satellite galaxies, we show that the most probable explanation for the detected flux excess is the presence of broad wings in the [CII] lines, signatures of starburst-driven outflows. We infer an average outflow rate of $\rm \dot{M}=54\pm23~ M_{\odot}~yr^{-1}$, providing a loading factor $\eta=\rm \dot{M}/SFR=1.7\pm1.3$ in agreement with observed local starbursts. Our interpretation is consistent with outcomes from zoomed hydro-simulations of {\it Dahlia}, a $z\sim 6$ galaxy ($\rm SFR\sim 100~\rm M_{\odot}~yr^{-1}$) whose feedback-regulated star formation results into an outflow rate $\rm \dot{M}\sim 30~ M_{\odot}~yr^{-1}$. The quality of the ALMA data is not sufficient for a detailed analysis of the [CII] line profile in individual galaxies. Nevertheless, our results suggest that starburst-driven outflows are in place in the EoR and provide useful indications for future ALMA campaigns. Deeper observations of the [CII] line in this sample are required to better characterise feedback at high-$z$ and to understand the role of outflows in shaping early galaxy formation., Comment: 10 pages, 7 figures, revised version to match the one accepted for publication in MNRAS

We present spectroscopic follow-up observations of CR7 with ALMA, targeted at constraining the infrared (IR) continuum and [CII]$_{158 \mu \rm m}$ line-emission at high spatial resolution matched to the HST/WFC3 imaging. CR7 is a luminous Ly$\alpha$ emitting galaxy at $z=6.6$ that consists of three separated UV-continuum components. Our observations reveal several well-separated components of [CII] emission. The two most luminous components in [CII] coincide with the brightest UV components (A and B), blue-shifted by $\approx 150$ km s$^{-1}$ with respect to the peak of Ly$\alpha$ emission. Other [CII] components are observed close to UV clumps B and C and are blue-shifted by $\approx300$ and $\approx80$ km s$^{-1}$ with respect to the systemic redshift. We do not detect FIR continuum emission due to dust with a 3$\sigma$ limiting luminosity L$_{\rm IR} (T_d = 35 \rm \, K) < 3.1\times10^{10}$ L$_{\odot}$. This allows us to mitigate uncertainties in the dust-corrected SFR and derive SFRs for the three UV clumps A, B and C of 28, 5 and 7 M$_{\odot}$ yr$^{-1}$. All clumps have [CII] luminosities consistent within the scatter observed in the local relation between SFR and L$_{\rm [CII]}$, implying that strong Ly$\alpha$ emission does not necessarily anti-correlate with [CII] luminosity. Combining our measurements with the literature, we show that galaxies with blue UV slopes have weaker [CII] emission at fixed SFR, potentially due to their lower metallicities and/or higher photoionisation. Comparison with hydrodynamical simulations suggests that CR7's clumps have metallicities of $0.1, Comment: 17 pages, accepted for publication in ApJ

We derive the CO luminosity function (LF) for different rotational transitions (i.e. (1-0), (3-2), (5-4)) starting from the Herschel LF by Gruppioni et al. and using appropriate $L'_{\rm CO} - L_{\rm IR}$ conversions for different galaxy classes. Our predicted LFs fit the data so far available at $z\approx0$ and $2$. We compare our results with those obtained by semi-analytical models (SAMs): while we find a good agreement over the whole range of luminosities at $z\approx0$, at $z\approx1$ and $z\approx2$ the tension between our LFs and SAMs in the faint and bright ends increases. We finally discuss the contribution of luminous AGN ($L_{X}>10^{44}\,\rm{erg\,s^{-1}}$) to the bright end of the CO LF concluding that they are too rare to reproduce the actual CO luminosity function at $z\approx2$., 5 pages, 3 figures, accepted for publication in MNRAS Letters

We present new ALMA observations of the [OIII]88$\mu$m line and high angular resolution observations of the [CII]158$\mu$m line in a normal star forming galaxy at z$=$7.1. Previous [CII] observations of this galaxy had detected [CII] emission consistent with the Ly$\alpha$ redshift but spatially slightly offset relative to the optical (UV-rest frame) emission. The new [CII] observations reveal that the [CII] emission is partly clumpy and partly diffuse on scales larger than about 1kpc. [OIII] emission is also detected at high significance, offset relative to the optical counterpart in the same direction as the [CII] clumps, but mostly not overlapping with the bulk of the [CII] emission. The offset between different emission components (optical/UV and different far-IR tracers) is similar to what observed in much more powerful starbursts at high redshift. We show that the [OIII] emitting clump cannot be explained in terms of diffuse gas excited by the UV radiation emitted by the optical galaxy, but it requires excitation by in-situ (slightly dust obscured) star formation, at a rate of about 7 M$_{\odot}$/yr. Within 20 kpc from the optical galaxy the ALMA data reveal two additional [OIII] emitting systems, which must be star forming companions. We discuss that the complex properties revealed by ALMA in the z$\sim$7.1 galaxy are consistent with expectations by recent models and cosmological simulations, in which differential dust extinction, differential excitation and different metal enrichment levels, associated with different subsystems assembling a galaxy, are responsible for the different appearance of the system when observed with different tracers., Comment: 12 pages, 8 figure, 4 tables

By combining high resolution, radiative transfer cosmological simulations of z~6 galaxies with a sub-grid multi-phase model of their interstellar medium we derive the expected intensity of several far infrared (FIR) emission lines ([C II] 158 micron, [O I] 63 micron, and [N II] 122 micron) for different values of the gas metallicity, Z. For Z = Z_sun the [C II] spectrum is very complex due to the presence of several emitting clumps of individual size < 3 kpc; the peak is displaced from the galaxy center by ~100 km/s. While the [O I] spectrum is also similarly displaced, the [N II] line comes predominantly from the central ionized regions of the galaxy. When integrated over ~500 km/s, the [C II] line flux is 185 mJy km/s; 95% of such flux originates from the cold (T~250 K) H I phase, and only 5% from the warm (T~5000 K) neutral medium. The [O I] and [N II] fluxes are ~6 and ~90 times lower than the [C II] one, respectively. By comparing our results with observations of Himiko, the most extended and luminous Lyman Alpha Emitter (LAE) at z = 6.6, we find that the gas metallicity in this source must be sub-solar. We conclude that the [C II] line from z~6 galaxies is detectable by the ALMA full array in 1.9 < t_ON < 7.7 hr observing time, depending on Z.

We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization., Comment: Accepted for publication in The Astrophysical Journal Letters. 6 pages, 5 figures

We present Atacama Large Millimiter/submillimiter Array (ALMA) observations of eight highly excited CO (J$_{\rm up}>8$) lines and continuum emission in two $z\sim6$ quasars: SDSS J231038.88+185519.7 (hereafter J2310), for which CO(8-7), CO(9-8), and CO(17-16) lines have been observed, and ULAS J131911.29+095951.4 (J1319), observed in the CO(14-13), CO(17-16) and CO(19-18) lines. The continuum emission of both quasars arises from a compact region ($< 0.9$ kpc). By assuming a modified black-body law, we estimate dust masses of Log$(M_{\rm dust}/M_{\odot})=8.75\pm0.07$ and Log$(M_{\rm dust}/M_{\odot})=8.8\pm0.2$ and dust temperatures of $T_{\rm dust}=76\pm3~{\rm K}$ and $T_{\rm dust}=66^{+15}_{-10}~{\rm K}$, respectively for J2310 and J1319. Only CO(8-7) and CO(9-8) in J2310 are detected, while $3��$ upper limits on luminosities are reported for the other lines of both quasars. The CO line luminosities and upper limits measured in J2310 and J1319 are consistent with those observed in local AGN and starburst galaxies, and other $z\sim 6$ quasars, except for SDSS J1148+5251 (J1148), the only quasar at $z=6.4$ with a previous CO(17-16) line detection. By computing the CO SLEDs normalised to the CO(6-5) line and FIR luminosities for J2310, J1319, and J1149, we conclude that different gas heating mechanisms (X-ray radiation and/or shocks) may explain the different CO luminosities observed in these $z\sim6$ quasar. Future J$_{\rm up}>8$ CO observations will be crucial to understand the processes responsible for molecular gas excitation in luminous high-$z$ quasars., 13 pages, 9 figures, 3 tables; MNRAS journal