Showing 1–50 of 86 results for author: Vallini, L

HYPERION. Shedding light on the first luminous quasars: A correlation between UV disc winds and X-ray continuum

Authors: A. Tortosa, L. Zappacosta, E. Piconcelli, M. Bischetti, C. Done, G. Miniutti, I. Saccheo, G. Vietri, A. Bongiorno, M. Brusa, S. Carniani, I. V. Chilingarian, F. Civano, S. Cristiani, V. D'Odorico, M. Elvis, X. Fan, C. Feruglio, F. Fiore, S. Gallerani, E. Giallongo, R. Gilli, A. Grazian, M. Guainazzi, F. Haardt , et al. (19 additional authors not shown)

Abstract: One of the main open questions in the field of luminous ($L_{\rm bol}>10^{47}\,\rm erg\,s^{-1}$) quasars (QSOs) at $z \gtrsim 6$ is the rapid formation ($< 1\,$Gyr) of their supermassive black holes (SMBHs). For this work we analysed the relation between the X-ray properties and other properties describing the physics and growth of both the accretion disc and the SMBH in QSOs at the Epoch of Reion… ▽ More One of the main open questions in the field of luminous ($L_{\rm bol}>10^{47}\,\rm erg\,s^{-1}$) quasars (QSOs) at $z \gtrsim 6$ is the rapid formation ($< 1\,$Gyr) of their supermassive black holes (SMBHs). For this work we analysed the relation between the X-ray properties and other properties describing the physics and growth of both the accretion disc and the SMBH in QSOs at the Epoch of Reionization (EoR). The sample consists of 21 $z>6$ QSOs, which includes 16 sources from the rapidly grown QSOs from the HYPERION sample and five other luminous QSOs with available high-quality archival X-ray data. We discovered a strong and statistically significant ($>3σ$) relation between the X-ray continuum photon index ($Γ$) and the $\rm C\,IV$ disc wind velocity ($v_{\rm C\,IV}$) in $z>6$ luminous QSOs, whereby the higher the $v_{\rm C\,IV}$, the steeper the $Γ$. This relation suggests a link between the disc-corona configuration and the kinematics of disc winds. Furthermore, we find evidence at $>2-3σ$ level that $Γ$ and $v_{\rm C\,IV}$ are correlated to the growth rate history of the SMBH. Although additional data are needed to confirm it, this result may suggest that, in luminous $z>6$ QSOs, the SMBH predominantly grows via fast accretion rather than via initial high seed BH mass. △ Less

Submitted 16 October, 2024; originally announced October 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

The mass-metallicity relation as a ruler for galaxy evolution: insights from the James Webb Space Telescope

Authors: A. Pallottini, A. Ferrara, S. Gallerani, L. Sommovigo, S. Carniani, L. Vallini, M. Kohandel, G. Venturi

Abstract: Galaxy evolution emerges from the balance between cosmic gas accretion, fueling star formation, and supernova (SN) feedback, regulating the metal enrichment. Hence, the stellar mass ($M_*$) - gas metallicity relation (MZR) is key to understand the physics of galaxies. High-quality JWST data enable accurate measurements of the MZR up to redshift z=10. Our aims are to understand the observed MZR, it… ▽ More Galaxy evolution emerges from the balance between cosmic gas accretion, fueling star formation, and supernova (SN) feedback, regulating the metal enrichment. Hence, the stellar mass ($M_*$) - gas metallicity relation (MZR) is key to understand the physics of galaxies. High-quality JWST data enable accurate measurements of the MZR up to redshift z=10. Our aims are to understand the observed MZR, its connection with the star formation rate (SFR), the role played by SFR stochasticity, and how it is regulated by SN feedback. We compare the MZR from the JADES, CEERS, and UNCOVER surveys, which comprise about 180 galaxies at $z=3-10$ with $10^6<M_*/M_\odot<10^{10}$, with 200 galaxies from the SERRA cosmological simulations. To interpret the MZR, we develop a minimal model for galaxy evolution that includes: cosmic accretion modulated with an amplitude $A_{100}$ on 100 Myr; a time delay $t_d$ between SFR and SN; SN-driven outflows with a varying mass loading factor $ε_{SN}$. Using our minimal model, we find the observed mean MZR is reproduced by weak outflows ($ε_{SN}=1/4$), in line with findings from JADES. Matching the observed MZR dispersion requires $t_d=20$ Myr and a $A_{100}=1/3$ modulation of the accretion rate. Successful models have low stochasticity ($σ_{SFR}=0.2$), yielding a MZR dispersion of $σ_{Z}=0.2$. Such values are close but lower than SERRA predictions ($σ_{SFR}=0.24$, $σ_{Z}=0.3$), clarifying why SERRA show no clear MZR trend and some tension with the observations. As the MZR is very sensitive to SFR stochasticity, models predicting high r.m.s. values ($σ_{SFR}=0.5$) result in a ``chemical chaos'' (i.e. $σ_{Z}=1.4$), virtually destroying the MZR. As a consequence, invoking a highly stochastic SFR ($σ_{SFR}=0.8$) to explain the overabundance of bright, super-early galaxies leads to inconsistencies with the observed MZR. △ Less

Submitted 31 July, 2024; originally announced August 2024.

Comments: 11 pages and 6 figures in the main text, submitted to A&A

A hidden AGN powering bright [O III] nebulae in a protocluster core at $z=4.5$ revealed by JWST

Authors: M. Solimano, J. González-López, M. Aravena, B. Alcalde Pampliega, R. J. Assef, M. Béthermin, M. Boquien, S. Bovino, C. M. Casey, P. Cassata, E. da Cunha, R. L. Davies, I. De Looze, X. Ding, T. Díaz-Santos, A. L. Faisst, A. Ferrara, D. B. Fisher, N. M. Förster-Schreiber, S. Fujimoto, M. Ginolfi, C. Gruppioni, L. Guaita, N. Hathi, R. Herrera-Camus , et al. (26 additional authors not shown)

Abstract: We present new JWST/NIRSpec IFU observations of the J1000+0234 system at $z=4.54$, the dense core of a galaxy protocluster hosting a massive, dusty star forming galaxy (DSFG) with a low luminosity radio counterpart. The new data reveals two extended, high equivalent width (EW$_0 > 1000$ Å) nebulae at each side of the DSFG disk along its minor axis (namely O3-N and O3-S). On one hand, O3-N's spectr… ▽ More We present new JWST/NIRSpec IFU observations of the J1000+0234 system at $z=4.54$, the dense core of a galaxy protocluster hosting a massive, dusty star forming galaxy (DSFG) with a low luminosity radio counterpart. The new data reveals two extended, high equivalent width (EW$_0 > 1000$ Å) nebulae at each side of the DSFG disk along its minor axis (namely O3-N and O3-S). On one hand, O3-N's spectrum shows a prominent FWHM $\sim1300$ km s$^{-1}$ broad and blueshifted component, suggesting an outflow origin. On the other hand, O3-S stretches over parsec and has a velocity gradient that spans $800$ km s$^{-1}$ but no evidence of a broad component. Both sources, however, seem to be powered at least partially by an active galactic nucleus (AGN), so we classify them as extended emission-line regions (EELRs). The strongest evidence comes from the detection of the high-ionization [Ne V] $\lambda3427$ line toward O3-N, which paired with the non-detection of hard X-rays implies an obscuring column density above the Compton-thick regime. In O3-S, the [Ne V] line is not detected, but we measure a He II well above the expectation for star formation. We interpret this as O3-S being externally irradiated by the AGN, akin to the famous Hanny's Voorwerp object in the local Universe. In addition, more classical line ratio diagnostics (e.g. [O III]/H$β$ vs [N II]/H$α$) put the DSFG itself in the AGN region of the diagrams, and hence the most probable host of the AGN. These results showcase the ability of JWST of unveiling highly obscured AGN at high redshifts. △ Less

Submitted 17 July, 2024; originally announced July 2024.

Comments: 5 pages, 4 figures plus 5 appendices (incl. 3 extra figures and one table). Submitted to A&A on July 17th 2024

A$^3$COSMOS: the dust mass function and dust mass density at $0.5<z<6$

Authors: A. Traina, B. Magnelli, C. Gruppioni, I. Delvecchio, M. Parente, F. Calura, L. Bisigello, A. Feltre, F. Pozzi, L. Vallini

Abstract: Context. Although dust in galaxies represents only a few percent of the total baryonic mass, it plays a crucial role in the physical processes occurring in galaxies. Studying the dust content of galaxies, particularly at high$-z$, is therefore crucial to understand the link between dust production, obscured star formation and the build-up of galaxy stellar mass. Aims. To study the dust propertie… ▽ More Context. Although dust in galaxies represents only a few percent of the total baryonic mass, it plays a crucial role in the physical processes occurring in galaxies. Studying the dust content of galaxies, particularly at high$-z$, is therefore crucial to understand the link between dust production, obscured star formation and the build-up of galaxy stellar mass. Aims. To study the dust properties (mass and temperature) of the largest Atacama Large Millimeter/submillimeter Array (ALMA)-selected sample of star-forming galaxies available from the archive (A$^3$COSMOS) and derive the dust mass function and dust mass density of galaxies from $z=0.5\,-\,6$. Methods. We performed spectral energy distribution (SED) fitting with the CIGALE code to constrain the dust mass and temperature of the A$^3$COSMOS galaxy sample, thanks to the UV-to-near-infrared photometric coverage of each galaxies combined with the ALMA (and Herschel when available) coverage of the Rayleigh-Jeans tail of their dust-continuum emission. We then computed and fitted the dust mass function by combining the A$^3$COSMOS and state-of-the-art {\it Herschel} samples, in order to obtain the best estimate of the integrated dust mass density up to $z \sim 6$. Results. Galaxies in \a3 have dust masses between $\sim 10^8$ and $\sim 10^{9.5}$ M$_{\odot}$. From the SED fitting, we were also able to derive a dust temperature, finding that the distribution of the dust temperature peaks at $\sim 30-35$K. The dust mass function at $z=0.5\,-\,6$ evolves with an increase of $M^*$ and decrease of the number density ($Φ^*$) and is in good agreement with literature estimates. The dust mass density shows a smooth decrease in its evolution from $z \sim 0.5$ to $z \sim 6$, which is steeper than what is found by models at $z \gtrsim 2$. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 11 pages, 9 figures, Accepted for publication on A&A

Journal ref: A&A 690, A84 (2024)

Gas conditions of a star-formation selected sample in the first billion years

Authors: Tom J. L. C. Bakx, Hiddo S. B. Algera, Bram Venemans, Laura Sommovigo, Seiji Fujimoto, Stefano Carniani, Masato Hagimoto, Takuya Hashimoto, Akio K. Inoue, Dragan Salak, Stephen Serjeant, Livia Vallini, Stephen Eales, Andrea Ferrara, Yoshinobu Fudamoto, Chihiro Imamura, Shigeki Inoue, Kirsten K. Knudsen, Hiroshi Matsuo, Yuma Sugahara, Yoichi Tamura, Akio Taniguchi, Satoshi Yamanaka

Abstract: We present Atacama Large Millimetre/submillimetre Array (ALMA) observations of the [O$_{\rm III}$] 88 $μ$m emission of a sample of thirteen galaxies at $z$ = 6 to 7.6 selected as [C$_{\rm II}$]-emitting companion sources of quasars. To disentangle the origins of the luminous Oxygen line in the $z$ > 6 Universe, we looked at emission-line galaxies that are selected through an excellent star-formati… ▽ More We present Atacama Large Millimetre/submillimetre Array (ALMA) observations of the [O$_{\rm III}$] 88 $μ$m emission of a sample of thirteen galaxies at $z$ = 6 to 7.6 selected as [C$_{\rm II}$]-emitting companion sources of quasars. To disentangle the origins of the luminous Oxygen line in the $z$ > 6 Universe, we looked at emission-line galaxies that are selected through an excellent star-formation tracer [C$_{\rm II}$] with star-formation rates between 9 and 162 M$_{\odot}$/yr. Direct observations reveal [O$_{\rm III}$] emission in just a single galaxy (L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ = 2.3), and a stacked image shows no [O$_{\rm III}$] detection, providing deep upper limits on the L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ ratios in the $z > 6$ Universe (L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ < 1.2 at 3$σ$). While the fidelity of this sample is high, no obvious optical/near-infrared counterpart is seen in the JWST imaging available for four galaxies. Additionally accounting for low-redshift CO emitters, line stacking shows that our sample-wide result remains robust: The enhanced L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ reported in the first billion years of the Universe is likely due to the selection towards bright, blue Lyman-break galaxies with high surface star-formation rates or young stellar populations. The deep upper limit on the rest-frame 90 $μ$m continuum emission (< 141 $μ$Jy at 3$σ$), implies a low average dust temperature (T$_{\rm dust}$ < 30K) and high dust mass (M$_{\rm dust}$ ~ 10$^8$ M$_{\odot}$). As more normal galaxies are explored in the early Universe, synergy between JWST and ALMA is fundamental to further investigate the ISM properties of the a broad range of samples of high-$z$ galaxies. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 20 pages; 13 figures; accepted for publication in MNRAS

The ALPINE-ALMA [CII] Survey: Dust emission effective radius up to 3 kpc in the Early Universe

Authors: F. Pozzi, F. Calura, Q. D'Amato, M. Gavarente, M. Bethermin, M. Boquien, V. Casasola, A. Cimatti, R. Cochrane, M. Dessauges-Zavadsky, A. Enia, F. Esposito, A. L. Faisst, R. Gilli, M. Ginolfi, R. Gobat, C. Gruppioni, C. C. Hayward, E. Ibar, A. M. Koekemoer, B. C. Lemaux, G. E. Magdis, J. Molina, M. Talia, L. Vallini , et al. (2 additional authors not shown)

Abstract: Measurements of the size of dust continuum emission are an important tool for constraining the spatial extent of star formation and hence the build-up of stellar mass. Compact dust emission has generally been observed at Cosmic Noon (z~2-3). However, at earlier epochs, toward the end of the Reionization (z~4-6), only the sizes of a handful of IR-bright galaxies have been measured. In this work, we… ▽ More Measurements of the size of dust continuum emission are an important tool for constraining the spatial extent of star formation and hence the build-up of stellar mass. Compact dust emission has generally been observed at Cosmic Noon (z~2-3). However, at earlier epochs, toward the end of the Reionization (z~4-6), only the sizes of a handful of IR-bright galaxies have been measured. In this work, we derive the dust emission sizes of main-sequence galaxies at z~5 from the ALPINE survey. We measure the dust effective radius r_e,FIR in the uv-plane in Band 7 of ALMA for seven ALPINE galaxies with resolved emission and we compare it with rest-frame UV and [CII]158$μ$m measurements. We study the r_e,FIR-L_IR scaling relation by considering our dust size measurements and all the data in literature at z~4-6. Finally, we compare our size measurements with predictions from simulations. The dust emission in the selected ALPINE galaxies is rather extended (r_e,FIR~1.5-3 kpc), similar to [CII]158 um but a factor of ~2 larger than the rest-frame UV emission. Putting together all the measurements at z~5, spanning 2 decades in luminosity from L_IR ~ 10^11 L_sun to L_IR ~ 10^13 L_sun, the data highlight a steeply increasing trend of the r_e,FIR-L_IR relation at L_IR< 10^12 L_sun, followed by a downturn and a decreasing trend at brighter luminosities. Finally, simulations that extend up to the stellar masses of the ALPINE galaxies considered in the present work predict a sub-set of galaxies (~25% at 10^10 M_sun < M_star < 10^11 M_sun) with sizes as large as those measured. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 8 pages, Accepted by A&A

Journal ref: A&A 686, A187 (2024)

AGN feedback in the Local Universe: multiphase outflow of the Seyfert galaxy NGC 5506

Authors: Federico Esposito, Almudena Alonso-Herrero, Santiago García-Burillo, Viviana Casasola, Françoise Combes, Daniele Dallacasa, Richard Davies, Ismael García-Bernete, Begoña García-Lorenzo, Laura Hermosa Muñoz, Luis Peralta de Arriba, Miguel Pereira-Santaella, Francesca Pozzi, Cristina Ramos Almeida, Thomas Taro Shimizu, Livia Vallini, Enrica Bellocchi, Omaira González-Martín, Erin K. S. Hicks, Sebastian Hönig, Alvaro Labiano, Nancy A. Levenson, Claudio Ricci, David J. Rosario

Abstract: We present new optical GTC/MEGARA seeing-limited (0.9") integral-field observations of NGC 5506, together with ALMA observations of the CO(3-2) transition at a 0.2" (25 pc) resolution. NGC 5506 is a luminous (bolometric luminosity of $\sim 10^{44}$ erg/s) nearby (26 Mpc) Seyfert galaxy, part of the Galaxy Activity, Torus, and Outflow Survey (GATOS). We modelled the CO(3-2) kinematics with 3D-Barol… ▽ More We present new optical GTC/MEGARA seeing-limited (0.9") integral-field observations of NGC 5506, together with ALMA observations of the CO(3-2) transition at a 0.2" (25 pc) resolution. NGC 5506 is a luminous (bolometric luminosity of $\sim 10^{44}$ erg/s) nearby (26 Mpc) Seyfert galaxy, part of the Galaxy Activity, Torus, and Outflow Survey (GATOS). We modelled the CO(3-2) kinematics with 3D-Barolo, revealing a rotating and outflowing cold gas ring within the central 1.2 kpc. We derived an integrated cold molecular gas mass outflow rate for the ring of 8 M$_{\odot}$/yr. We fitted the optical emission lines with a maximum of two Gaussian components to separate rotation from non-circular motions. We detected high [OIII]$λ$5007 projected velocities (up to 1000 km/s) at the active galactic nucleus (AGN) position, decreasing with radius to an average 330 km/s around 350 pc. We also modelled the [OIII] gas kinematics with a non-parametric method, estimating the ionisation parameter and electron density in every spaxel, from which we derived an ionised mass outflow rate of 0.076 M$_{\odot}$/yr within the central 1.2 kpc. Regions of high CO(3-2) velocity dispersion, extending to projected distances of 350 pc from the AGN, appear to be the result from the interaction of the AGN wind with molecular gas in the galaxy's disc. Additionally, we find the ionised outflow to spatially correlate with radio and soft X-ray emission in the central kiloparsec. We conclude that the effects of AGN feedback in NGC 5506 manifest as a large-scale ionised wind interacting with the molecular disc, resulting in outflows extending to radial distances of 610 pc △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 24 pages, 30 figures, accepted for publication by Astronomy & Astrophysics

Gas-phase metallicity gradients in galaxies at $z \sim 6-8$

Authors: G. Venturi, S. Carniani, E. Parlanti, M. Kohandel, M. Curti, A. Pallottini, L. Vallini, S. Arribas, A. J. Bunker, A. J. Cameron, M. Castellano, A. Ferrara, A. Fontana, S. Gallerani, V. Gelli, R. Maiolino, E. Ntormousi, C. Pacifici, L. Pentericci, S. Salvadori, E. Vanzella

Abstract: The study of gas-phase metallicity and its spatial distribution at high redshift is crucial to understand the processes that shaped the growth and evolution of galaxies in the early Universe. Here we study the spatially resolved metallicity in three systems at $z\sim6-8$, namely A2744-YD4, BDF-3299, and COSMOS24108, with JWST NIRSpec IFU low-resolution ($R\sim100$) spectroscopic observations. Thes… ▽ More The study of gas-phase metallicity and its spatial distribution at high redshift is crucial to understand the processes that shaped the growth and evolution of galaxies in the early Universe. Here we study the spatially resolved metallicity in three systems at $z\sim6-8$, namely A2744-YD4, BDF-3299, and COSMOS24108, with JWST NIRSpec IFU low-resolution ($R\sim100$) spectroscopic observations. These are among the highest-$z$ sources in which metallicity gradients have been probed so far. Each of these systems hosts several spatial components in the process of merging within a few kpc, identified from the rest-frame UV and optical stellar continuum and ionised gas emission line maps. The sources have heterogeneous properties, with stellar masses log($M_*/M_\odot) \sim 7.6-9.3$, star formation rates (SFRs) $\sim1-15$ $M_\odot$ yr$^{-1}$, and gas-phase metallicities 12+log(O/H) $\sim 7.7-8.3$, which exhibit a large scatter within each system. Their properties are generally consistent with those of the highest-$z$ samples to date ($z\sim3-10$), though the sources in A2744-YD4 and COSMOS24108 are at the high end of the mass-metallicity relation (MZR) defined by the $z\sim3-10$ sources. Moreover, the targets in this work follow the predicted slope of the MZR at $z\sim 6-8$ from most cosmological simulations. The gas-phase metallicity gradients are consistent with being flat in the main sources of each system. Flat metallicity gradients are thought to arise from gas mixing processes on galaxy scales, such as mergers or galactic outflows and SN winds driven by intense stellar feedback, which wash out any gradient formed in the galaxy. The existence of flat gradients at $z\sim6-8$ sets also important constraints on cosmological simulations and chemical evolution models, whose predictions on the cosmic evolution of metallicity gradients differ significantly, but are mostly limited to $z<3$ so far. △ Less

Submitted 10 September, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: 15 pages (+5 appendix); 7 figures (+8 in appendix); in production on A&A

Unveiling [CII] clumps in a lensed star-forming galaxy at z ~ 3.4

Authors: A. Zanella, E. Iani, M. Dessauges-Zavadsky, J. Richard, C. De Breuck, J. Vernet, M. Kohandel, F. Arrigoni Battaia, A. Bolamperti, F. Calura, C. -C. Chen, T. Devereaux, A. Ferrara, V. Mainieri, A. Pallottini, G. Rodighiero, L. Vallini, E. Vanzella

Abstract: Observations at UV and optical wavelengths have revealed that galaxies at z~1-4 host star-forming regions, dubbed "clumps", which are believed to form due to the fragmentation of gravitationally unstable, gas-rich disks. However, the detection of the parent molecular clouds that give birth to such clumps is still possible only in a minority of galaxies, mostly at z~1. We investigated the [CII] and… ▽ More Observations at UV and optical wavelengths have revealed that galaxies at z~1-4 host star-forming regions, dubbed "clumps", which are believed to form due to the fragmentation of gravitationally unstable, gas-rich disks. However, the detection of the parent molecular clouds that give birth to such clumps is still possible only in a minority of galaxies, mostly at z~1. We investigated the [CII] and dust morphology of a z~3.4 lensed galaxy hosting four clumps detected in the UV continuum. We aimed to observe the [CII] emission of individual clumps that, unlike the UV, is not affected by dust extinction, to probe their nature and cold gas content. We conducted ALMA observations probing scales down to ~300 pc and detected three [CII] clumps. One (dubbed "NE") coincides with the brightest UV clump, while the other two ("SW" and "C") are not detected in the UV continuum. We do not detect the dust continuum. We converted the [CII] luminosity of individual clumps into molecular gas mass and found Mmol~10^8 Msun. By complementing it with the star formation rate (SFR) estimate from the UV continuum, we estimated the gas depletion time (tdep) of clumps and investigated their location in the Schmidt-Kennicutt plane. While the NE clump has a short tdep=0.16 Gyr, comparable with high-redshift starbursts, the SW and C clumps instead have longer tdep>0.65 Gyr and are likely probing the initial phases of star formation. The lack of dust continuum detection is consistent with the blue UV continuum slope estimated for this galaxy (beta~-2.5) and it indicates that dust inhomogeneities do not significantly affect the detection of UV clumps in this target. We pushed the observation of the cold gas content of individual clumps up to z~3.4 and showed that the [C II] line emission is a promising tracer of molecular clouds at high redshift, allowing the detection of clumps with a large range of depletion times. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A

The ALPINE-ALMA [C II] survey: Characterisation of Spatial Offsets in Main-Sequence Galaxies at $z \sim$ 4-6

Authors: Meghana Killi, Michele Ginolfi, Gergö Popping, Darach Watson, Giovanni Zamorani, Brian C. Lemaux, Seiji Fujimoto, Andreas Faisst, Matthieu Bethermin, Michael Romano, Yoshinobu Fudamoto, Sandro Bardelli, Médéric Boquien, Stefano Carniani, Miroslava Dessauges-Zavadsky, Carlotta Gruppioni, Nimish Hathi, Eduardo Ibar, Gareth C. Jones, Anton M. Koekemoer, Ivanna Langan, Hugo Méndez-Hernández, Yuma Sugahara, Livia Vallini, Daniela Vergani

Abstract: Galaxy morphology is shaped by stellar activity, feedback, gas and dust properties, and interactions with surroundings, and can therefore provide insight into these processes. In this paper, we study the spatial offsets between stellar and interstellar medium emission in a sample of 54 main-sequence star-forming galaxies at $z\sim4-6$ observed with the Atacama Large Millimeter/submillimeter Array… ▽ More Galaxy morphology is shaped by stellar activity, feedback, gas and dust properties, and interactions with surroundings, and can therefore provide insight into these processes. In this paper, we study the spatial offsets between stellar and interstellar medium emission in a sample of 54 main-sequence star-forming galaxies at $z\sim4-6$ observed with the Atacama Large Millimeter/submillimeter Array (ALMA) and drawn from the ALMA Large Program to INvestigate C$^+$ at Early times (ALPINE). We find no significant spatial offset for the majority ($\sim$ 70 percent) of galaxies in the sample among any combination of [C II], far-infrared continuum, optical, and ultraviolet emission. However, a fraction of the sample ($\sim$ 30 percent) shows offsets larger than the median by more than 3$σ$ significance (compared to the uncertainty on the offsets), especially between [C II] and ultraviolet emission. We find that these significant offsets are of the order of $\sim$0.5-0.7 arcsec, corresponding to $\sim$3.5-4.5 kiloparsecs. The offsets could be caused by a complex dust geometry, strong feedback from stars and active galactic nuclei, large-scale gas inflow and outflow, or a combination of these phenomena. However, our current analysis does not definitively constrain the origin. Future, higher resolution ALMA and JWST observations may help resolve the ambiguity. Regardless, since there exist at least some galaxies that display such large offsets, galaxy models and spectral energy distribution fitting codes cannot assume co-spatial emission in all main-sequence galaxies, and must take into account that the observed emission across wavelengths may be spatially segregated. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 20 pages, 8 figures. Submitted to MNRAS; second revision

Dark progenitors and massive descendants: A first ALMA perspective on Radio-Selected NIRdark galaxies in the COSMOS field

Authors: Fabrizio Gentile, Margherita Talia, Emanuele Daddi, Marika Giulietti, Andrea Lapi, Marcella Massardi, Francesca Pozzi, Giovanni Zamorani, Meriem Behiri, Andrea Enia, Matthieu Bethermin, Daniele Dallacasa, Ivan Delvecchio, Andreas L. Faisst, Carlotta Gruppioni, Federica Loiacono, Alberto Traina, Mattia Vaccari, Livia Vallini, Cristian Vignali, Vernesa Smolcic, Andrea Cimatti

Abstract: We present the first spectroscopic ALMA follow-up for a pilot sample of nine Radio-Selected NIRdark galaxies in the COSMOS field. These sources were initially selected as radio-detected sources (S(3GHz)>12.65 uJy), lacking an optical/NIR counterpart in the COSMOS2015 catalog (Ks>24.7 mag), with just three of them subsequently detected in the deeper COSMOS2020. Several studies highlighted how this… ▽ More We present the first spectroscopic ALMA follow-up for a pilot sample of nine Radio-Selected NIRdark galaxies in the COSMOS field. These sources were initially selected as radio-detected sources (S(3GHz)>12.65 uJy), lacking an optical/NIR counterpart in the COSMOS2015 catalog (Ks>24.7 mag), with just three of them subsequently detected in the deeper COSMOS2020. Several studies highlighted how this selection could provide a population of highly dust-obscured, massive, and star-bursting galaxies. With these new ALMA observations, we assess the spectroscopic redshifts of this pilot sample of sources and improve the quality of the physical properties estimated through SED-fitting. Moreover, we measure the quantity of molecular gas present inside these galaxies and forecast their potential evolutionary path, finding that the RS-NIRdark galaxies could represent a likely population of high-z progenitors of the massive and passive galaxies discovered at z~3. Finally, we present some initial constraints on the kinematics of the ISM within the analyzed galaxies, reporting a high fraction (~55%) of double-peaked lines that can be interpreted as the signature of a rotating structure in our targets or with the presence of major mergers in our sample. Our results presented in this paper showcase the scientific potential of (sub)mm observations for this elusive population of galaxies and highlight the potential contribution of these sources in the evolution of the massive and passive galaxies at high-z. △ Less

Submitted 13 May, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 18 pages, 9 figures (+2 in the appendix), 6 tables. Accepted by A&A

Journal ref: A&A 687, A288 (2024)

Modelling molecular clouds and CO excitation in AGN-host galaxies

Authors: Federico Esposito, Livia Vallini, Francesca Pozzi, Viviana Casasola, Almudena Alonso-Herrero, Santiago García-Burillo, Roberto Decarli, Francesco Calura, Cristian Vignali, Matilde Mingozzi, Carlotta Gruppioni, Dhrubojyoti Sengupta

Abstract: We present a new physically-motivated model for estimating the molecular line emission in active galaxies. The model takes into account (i) the internal density structure of giant molecular clouds (GMCs), (ii) the heating associated both to stars and to the active galactic nuclei (AGN), respectively producing photodissociation regions (PDRs) and X-ray dominated regions (XDRs) within the GMCs, and… ▽ More We present a new physically-motivated model for estimating the molecular line emission in active galaxies. The model takes into account (i) the internal density structure of giant molecular clouds (GMCs), (ii) the heating associated both to stars and to the active galactic nuclei (AGN), respectively producing photodissociation regions (PDRs) and X-ray dominated regions (XDRs) within the GMCs, and (iii) the mass distribution of GMCs within the galaxy volume. The model needs, as input parameters, the radial profiles of molecular mass, far-UV flux and X-ray flux for a given galaxy, and it has two free parameters: the CO-to-H2 conversion factor $α_{CO}$, and the X-ray attenuation column density $N_H$. We test this model on a sample of 24 local ($z \leq 0.06$) AGN-host galaxies, simulating their carbon monoxide spectral line energy distribution (CO SLED). We compare the results with the available observations and calculate, for each galaxy, the best ($α_{CO}$, $N_H$) with a Markov chain Monte Carlo algorithm, finding values consistent with those present in the literature. We find a median $α_{CO} = 4.8$ M$_{\odot}$ (K km s$^{-1}$ pc$^{2}$)$^{-1}$ for our sample. In all the modelled galaxies, we find the XDR component of the CO SLED to dominate the CO luminosity from $J_{\text{upp}} \geq 4$. We conclude that, once a detailed distribution of molecular gas density is taken into account, PDR emission at mid-/high-$J$ becomes negligible with respect to XDR. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 19 pages, 36 figures, resubmitted to MNRAS after addressing referee's comments

SERENADE II: An ALMA Multi-Band Dust-Continuum Analysis of 28 Galaxies at $5<z<8$ and the Physical Origin of the Dust Temperature Evolution

Authors: Ikki Mitsuhashi, Yuichi Harikane, Franz E. Bauer, Tom Bakx, Andrea Ferrara, Seiji Fujimoto, Takuya Hashimoto, Akio K. Inoue, Kazushi Iwasawa, Yuri Nishimura, Masatoshi Imanishi, Yoshiaki Ono, Toshiki Saito, Yuma Sugahara, Hideki Umehata, Livia Vallini, Tao Wang

Abstract: We present an analysis of ALMA multi-band dust-continuum observations for 28 spectroscopically-confirmed bright Lyman-break galaxies at $5<z<8$. Our sample consists of 11 galaxies at $z\sim6$ newly observed in our ALMA program, which substantially increases the number of $5<z<8$ galaxies with both rest-frame 88 and 158 $μ{\rm m}$ continuum observations, allowing us to simultaneously measure the IR… ▽ More We present an analysis of ALMA multi-band dust-continuum observations for 28 spectroscopically-confirmed bright Lyman-break galaxies at $5<z<8$. Our sample consists of 11 galaxies at $z\sim6$ newly observed in our ALMA program, which substantially increases the number of $5<z<8$ galaxies with both rest-frame 88 and 158 $μ{\rm m}$ continuum observations, allowing us to simultaneously measure the IR luminosity and dust temperature for a statistical sample of $z\gtrsim5$ galaxies for the first time. We derive the relationship between the UV slope ($β_{\rm UV}$) and infrared excess (IRX) for the $z\sim6$ galaxies, and find a shallower IRX-$β_{\rm UV}$ relation compared to the previous results at $z\sim2$--4. Based on the IRX-$β_{\rm UV}$ relation consistent with our results and the $β_{\rm UV}$-$M_{\rm UV}$ relation including fainter galaxies in the literature, we find a limited contribution of the dust-obscured star formation to the total SFR density, $\sim30\%$ at $z\sim6$. Our measurements of the dust temperature at $z\sim6-7$, $T_{\rm dust}=40.9_{-9.1}^{+10.0}\,{\rm K}$ on average, supports a gentle increase of $T_{\rm dust}$ from $z=0$ to $z\sim6$--7. Using an analytic model with parameters consistent with recent {\it{JWST}} results, we discuss that the observed redshift evolution of the dust temperature can be reproduced by an $\sim0.6\,{\rm dex}$ increase in the gas depletion timescale and $\sim0.4\,{\rm dex}$ decrease of the metallicity. The variety of $T_{\rm dust}$ observed at high redshifts can also be naturally explained by scatters around the star-formation main sequence and average mass-metallicity relation, including an extremely high dust temperature of $T_{\rm dust}>80\,{\rm K}$ observed in a galaxy at $z=8.3$. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: Submitted to ApJ

The ALMA-ALPINE [CII] survey: Kennicutt-Schmidt relation in four massive main-sequence galaxies at z~4.5

Authors: M. Béthermin, C. Accard, C. Guillaume, M. Dessauges-Zavadsky, E. Ibar, P. Cassata, T. Devereaux, A. Faisst, J. Freundlich, G. C. Jones, K. Kraljic, H. Algera, R. O. Amorin, S. Bardelli, M. Boquien, V. Buat, E. Donghia, Y. Dubois, A. Ferrara, Y. Fudamoto, M. Ginolfi, P. Guillard, M. Giavalisco, C. Gruppioni, G. Gururajan , et al. (18 additional authors not shown)

Abstract: The Kennicutt-Schmidt (KS) relation between the gas and the star formation rate (SFR) surface density ($Σ_{\rm gas}$-$Σ_{\rm SFR}$) is essential to understand star formation processes in galaxies. So far, it has been measured up to z~2.5 in main-sequence galaxies. In this letter, we aim to put constraints at z~4.5 using a sample of four massive main-sequence galaxies observed by ALMA at high resol… ▽ More The Kennicutt-Schmidt (KS) relation between the gas and the star formation rate (SFR) surface density ($Σ_{\rm gas}$-$Σ_{\rm SFR}$) is essential to understand star formation processes in galaxies. So far, it has been measured up to z~2.5 in main-sequence galaxies. In this letter, we aim to put constraints at z~4.5 using a sample of four massive main-sequence galaxies observed by ALMA at high resolution. We obtained ~0.3"-resolution [CII] and continuum maps of our objects, which we then converted into gas and obscured SFR surface density maps. In addition, we produced unobscured SFR surface density maps by convolving Hubble ancillary data in the rest-frame UV. We then derived the average $Σ_{\rm SFR}$ in various $Σ_{\rm gas}$ bins, and estimated the uncertainties using a Monte Carlo sampling. Our galaxy sample follows the KS relation measured in main-sequence galaxies at lower redshift and is slightly lower than predictions from simulations. Our data points probe the high end both in terms of $Σ_{\rm gas}$ and $Σ_{\rm gas}$, and gas depletion timescales (285-843 Myr) remain similar to z~2 objects. However, three of our objects are clearly morphologically disturbed, and we could have expected shorter gas depletion timescales (~100 Myr) similar to merger-driven starbursts at lower redshifts. This suggests that the mechanisms triggering starbursts at high redshift may be different than in the low- and intermediate-z Universe. △ Less

Submitted 17 November, 2023; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: 9 pages, 7 figures, 2 tables, accepted by A&A (letter)

Studying [CII] Emission in Low-mass Galaxies at z ~ 7

Authors: Kelsey Glazer, Marusa Bradac, Ryan L. Sanders, Seiji Fujimoto, Patricia Bolan, Andrea Ferrara, Victoria Strait, Tucker Jones, Brian C. Lemaux, Livia Vallini, Russell Ryan

Abstract: We report on a $\rm{[CII]}_{158μ\rm{m}}$ search using the Atacama Large Millimeter/submillimeter Array (ALMA) on three lensed, confirmed {\lya} emitting galaxies at $z \sim 7$. Our targets are ultra-violet (UV) faint systems with stellar masses on the order of $M_{*} \sim 10^{9} M_{\odot}$. We detect a single [CII] line emission ($4σ$) from the brightest ($L \sim 2.4 \times 10^{10}L_{\odot}$) gala… ▽ More We report on a $\rm{[CII]}_{158μ\rm{m}}$ search using the Atacama Large Millimeter/submillimeter Array (ALMA) on three lensed, confirmed {\lya} emitting galaxies at $z \sim 7$. Our targets are ultra-violet (UV) faint systems with stellar masses on the order of $M_{*} \sim 10^{9} M_{\odot}$. We detect a single [CII] line emission ($4σ$) from the brightest ($L \sim 2.4 \times 10^{10}L_{\odot}$) galaxy in our sample, MACS0454-1251. We determine a systemic redshift ($z_{\rm{[CII]}} = 6.3151 \pm 0.0005$) for MACS0454-1251 and measure a {\lya} velocity offset of $Δv \approx 300 \pm 70 \rm{km\,s}^{-1}$. The remaining two galaxies we detect no {\ct} but provide $3 σ$ upper limits on their {\ct} line luminosities which we use to investigate the $L_{\textrm{[CII]}} - \rm{SFR}$ relation. Overall our single {\ct} detection shows agreement with the relation for dwarf and local starburst galaxies. Our [CII] deficient galaxies could potentially be exhibiting low metallicities ($Z<Z_{\odot}$). Another possible explanation for weaker [CII] emission could be strong feedback from star formation disrupting molecular clouds. We do not detect continuum emission in any of the sources, placing upper limits on their dust masses. Assuming a single dust temperature of $T_{d}=35 \rm{K}$ dust masses ($M_{\rm{dust}}$) range from $< 4.8 \times 10^{7} M_{\odot} $ to $2.3 \times 10^{8} M_{\odot}$. Collectively, our results suggest faint reionization era sources could be metal poor and/or could have strong feedback suppressing [CII] emission. △ Less

Submitted 30 April, 2024; v1 submitted 20 September, 2023; originally announced September 2023.

Comments: 7 pages, 4 figures accepted for publication to MNRAS

Spatially resolved Kennicutt-Schmidt relation at z~7 and its connection with the interstellar medium properties

Authors: Livia Vallini, Joris Witstok, Laura Sommovigo, Andrea Pallottini, Andrea Ferrara, Stefano Carniani, Mahsa Kohandel, Renske Smit, Simona Gallerani, Carlotta Gruppioni

Abstract: We exploit moderately resolved [OIII], [CII] and dust continuum ALMA observations to derive the gas density ($n$), the gas-phase metallicity ($Z$) and the deviation from the Kennicutt-Schmidt (KS) relation ($κ_s$) on ~sub-kpc scales in the interstellar medium (ISM) of five bright Lyman Break Galaxies at the Epoch of Reionization ($z\approx 7$). To do so, we use GLAM, a state-of-art, physically mot… ▽ More We exploit moderately resolved [OIII], [CII] and dust continuum ALMA observations to derive the gas density ($n$), the gas-phase metallicity ($Z$) and the deviation from the Kennicutt-Schmidt (KS) relation ($κ_s$) on ~sub-kpc scales in the interstellar medium (ISM) of five bright Lyman Break Galaxies at the Epoch of Reionization ($z\approx 7$). To do so, we use GLAM, a state-of-art, physically motivated Bayesian model that links the [CII] and [OIII] surface brightness ($Σ_{\rm [CII]}$, $Σ_{\rm [OIII]}$) and the SFR surface density ($Σ_{\rm SFR}$) to $n$, $κ_s$, and $Z$. All five sources are characterized by a central starbursting region, where the $Σ_{\rm gas}$ vs $Σ_{\rm SFR}$ align ~10x above the KS relation ($κ_s\approx10$). This translates into gas depletion times in the range $t_{\rm dep}\approx 80-250$ Myr. The inner starbursting centers are characterized by higher gas density ($\log (n/{\rm cm^{-3}}) \approx 2.5-3.0$) and higher metallicity ($\log (Z/Z_{\odot}) \approx -0.5$) than the galaxy outskirts. We derive marginally negative radial metallicity gradients ($\nabla \log Z \approx -0.03 \pm 0.07$dex/kpc), and a dust temperature ($T_d\approx$32-38 K) that anticorrelates with the gas depletion time. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 13 pages, 7 figures. Resubmitted to MNRAS after addressing the referee comments

HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). A new regime for the X-ray nuclear properties of the first quasars

Authors: L. Zappacosta, E. Piconcelli, F. Fiore, I. Saccheo, R. Valiante, C. Vignali, F. Vito, M. Volonteri, M. Bischetti, A. Comastri, C. Done, M. Elvis, E. Giallongo, F. La Franca, G. Lanzuisi, M. Laurenti, G. Miniutti, A. Bongiorno, M. Brusa, F. Civano, S. Carniani, V. D'Odorico, C. Feruglio, S. Gallerani, R. Gilli , et al. (18 additional authors not shown)

Abstract: 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 Q… ▽ More 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 ($Γ\sim2.4\pm0.1$). Absorption is not required. The average $Γ$ is inconsistent at $\geq4σ$ level with the canonical 1.8-2 value measured in QSO at z<6. This spectral slope is also much steeper than that reported in lower-z QSOs with similar luminosity or accretion rate, thus suggesting a genuine redshift evolution. Alternatively, we can interpret this result as the presence of an unusually low-energy cutoff $E_{cut}\sim20$ keV on a standard $Γ=1.9$ power-law. We also report on mild indications that HYPERION QSOs show higher soft X-ray emission at 2 keV compared to the UV one at 2500A than expected by lower-z luminous AGN. We speculate that a redshift-dependent coupling between the corona and accretion disc or intrinsically different coronal properties may account for the steep spectral slopes, especially in the presence of powerful winds. The reported slopes, if confirmed at lower luminosities, may have an important impact on future X-ray AGN studies in the early Universe. △ Less

Submitted 19 July, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 21 pages (including appendix), 12 figures, 4 tables. Accepted for pubblication in A&A

Journal ref: A&A 678, A201 (2023)

Optical and mid-infrared line emission in nearby Seyfert galaxies

Authors: A. Feltre, C. Gruppioni, L. Marchetti, A. Mahoro, F. Salvestrini, M. Mignoli, L. Bisigello, F. Calura, S. Charlot, J. Chevallard, E. Romero-Colmenero, E. Curtis-Lake, I. Delvecchio, O. L. Dors, M. Hirschmann, T. Jarrett, S. Marchesi, M. E. Moloko, A. Plat, F. Pozzi, R. Sefako, A. Traina, M. Vaccari, P. Väisänen, L. Vallini , et al. (2 additional authors not shown)

Abstract: Line ratio diagnostics provide valuable clues on the source of ionizing radiation in galaxies with intense black hole accretion and starbursting events, such as local Seyfert or galaxies at the peak of the star formation history. We aim to provide a reference joint optical and mid-IR analysis for studying AGN identification via line ratios and testing predictions from photoionization models. We ob… ▽ More Line ratio diagnostics provide valuable clues on the source of ionizing radiation in galaxies with intense black hole accretion and starbursting events, such as local Seyfert or galaxies at the peak of the star formation history. We aim to provide a reference joint optical and mid-IR analysis for studying AGN identification via line ratios and testing predictions from photoionization models. We obtained homogenous optical spectra with the Southern Africa Large Telescope for 42 Seyfert galaxies with Spitzer/IRS spectroscopy and X-ray to mid-IR multiband data available. After confirming the power of the main optical ([OIII]) and mid-IR ([NeV], [OIV], [NeIII]) emission lines in tracing AGN activity, we explore diagrams based on ratios of optical and mid-IR lines by exploiting photoionization models of different ionizing sources (AGN, star formation and shocks). We find that pure AGN photoionization models are good at reproducing observations of Seyfert galaxies with an AGN fractional contribution to the mid-IR (5-40 micron) emission larger than 50 per cent. For targets with a lower AGN contribution these same models do not fully reproduce the observed mid-IR line ratios. Mid-IR ratios like [NeV]/[NeII], [OIV]/[NeII] and [NeIII]/[NeII] show a dependence on the AGN fractional contribution to the mid-IR unlike optical line ratios. An additional source of ionization, either from star formation or radiative shocks, can help explain the observations in the mid-IR. Among combinations of optical and mid-IR diagnostics in line ratio diagrams, only those involving the [OI]/Halpha ratio are promising diagnostics for simultaneously unraveling the relative role of AGN, star formation and, shocks. A proper identification of the dominant ionizing source would require the exploitation of analysis tools based on advanced statistical techniques as well as spatially resolved data. △ Less

Submitted 19 May, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 31 pages, 15 figures, 2 tables. Accepted for publication in A&A

Interpreting ALMA non-detections of JWST super-early galaxies

Authors: M. Kohandel, A. Ferrara, A. Pallottini, L. Vallini, L. Sommovigo, F. Ziparo

Abstract: Recent attempts to detect [OIII] 88$μ$m emission from super-early ($z>10$) galaxy candidates observed by JWST have been unsuccessful. By using zoom-in simulations, we show that these galaxies are faint, and mostly fall below the local metal-poor $\rm [OIII]-SFR$ relation as a result of their low ionization parameter, $U_{\rm ion}\lesssim 10^{-3}$. Such low $U_{\rm ion}$ values are found in galaxie… ▽ More Recent attempts to detect [OIII] 88$μ$m emission from super-early ($z>10$) galaxy candidates observed by JWST have been unsuccessful. By using zoom-in simulations, we show that these galaxies are faint, and mostly fall below the local metal-poor $\rm [OIII]-SFR$ relation as a result of their low ionization parameter, $U_{\rm ion}\lesssim 10^{-3}$. Such low $U_{\rm ion}$ values are found in galaxies that are in an early assembly stage, and whose stars are still embedded in high-density natal clouds. However, the most luminous galaxy in our sample ($\rm{log}[L_{\rm{[OIII]}}/L_\odot] = 8.4$, $U_{\rm ion} \approx 0.1$) could be detected by ALMA in only $2.8$ hrs. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 6 pages, 3 figures, Submitted to MNRAS (comments are welcome)

CONCERTO: Simulating the CO, [CII], and [CI] line emission of galaxies in a 117 $\rm deg^2$ field and the impact of field-to-field variance

Authors: A. Gkogkou, M. Béthermin, G. Lagache, M. Van Cuyck, E. Jullo, M. Aravena, A. Beelen, A. Benoit, J. Bounmy, M. Calvo, A. Catalano, S. Cora, D. Croton, S. de la Torre, A. Fasano, A. Ferrara, J. Goupy, C. Hoarau, W. Hu, T. Ishiyama, K. K. Knudsen, J. -C. Lambert, J. F. Macías-Pérez, J. Marpaud, G. Mellema , et al. (7 additional authors not shown)

Abstract: In the submm regime, spectral line scans and line intensity mapping (LIM) are new promising probes for the cold gas content and star formation rate of galaxies across cosmic time. However, both of these two measurements suffer from field-to-field variance. We study the effect of field-to-field variance on the predicted CO and [CII] power spectra from future LIM experiments such as CONCERTO, as wel… ▽ More In the submm regime, spectral line scans and line intensity mapping (LIM) are new promising probes for the cold gas content and star formation rate of galaxies across cosmic time. However, both of these two measurements suffer from field-to-field variance. We study the effect of field-to-field variance on the predicted CO and [CII] power spectra from future LIM experiments such as CONCERTO, as well as on the line luminosity functions (LFs) and the cosmic molecular gas mass density that are currently derived from spectral line scans. We combined a 117 $\rm deg^2$ dark matter lightcone from the Uchuu cosmological simulation with the simulated infrared dusty extragalactic sky (SIDES) approach. We find that in order to constrain the CO LF with an uncertainty below 20%, we need survey sizes of at least 0.1 $\rm deg^2$. Furthermore, accounting for the field-to-field variance using only the Poisson variance can underestimate the total variance by up to 80%. The lower the luminosity is and the larger the survey size is, the higher the level of underestimate. At $z$<3, the impact of field-to-field variance on the cosmic molecular gas density can be as high as 40% for the 4.6 arcmin$^2$ field, but drops below 10% for areas larger than 0.2 deg$^2$. However, at $z>3$ the variance decreases more slowly with survey size and for example drops below 10% for 1 deg$^2$ fields. Finally, we find that the CO and [CII] LIM power spectra can vary by up to 50% in $\rm 1 deg^2$ fields. This limits the accuracy of the constraints provided by the first 1 deg$^2$ surveys. The level of the shot noise power is always dominated by the sources that are just below the detection thresholds. We provide an analytical formula to estimate the field-to-field variance of current or future LIM experiments. The code and the full SIDES-Uchuu products (catalogs, cubes, and maps) are publicly available. △ Less

Submitted 5 December, 2022; originally announced December 2022.

The ALPINE-ALMA [CII] survey: double stellar population and AGN activity in a galaxy at $z\sim5.5$

Authors: L. Barchiesi, M. Dessauges-Zavadsky, C. Vignali, F. Pozzi, R. Marques-Chaves, A. Feltre, A. Faisst, M. Béthermin, P. Cassata, S. Charlot, Y. Fudamoto, M. Ginolfi, E. Ibar, G. C. Jones, M. Romano, D. Schaerer, L. Vallini, E. Vanzella, L. Yan

Abstract: GDS J033218.92-275302.7 (here GS-14) is a $z\sim5.5$ galaxy detected in [CII] as part of the ALPINE survey with unusual UV spectral features that have been interpreted as signatures of either a double stellar population or of an active galactic nucleus (AGN). We exploited the multi-wavelength coverage of GS-14 to investigate the properties and the origin of its emission. We performed UV-to-NIR SED… ▽ More GDS J033218.92-275302.7 (here GS-14) is a $z\sim5.5$ galaxy detected in [CII] as part of the ALPINE survey with unusual UV spectral features that have been interpreted as signatures of either a double stellar population or of an active galactic nucleus (AGN). We exploited the multi-wavelength coverage of GS-14 to investigate the properties and the origin of its emission. We performed UV-to-NIR SED-fitting, with single/double stellar population and/or AGN component. We analyzed the VIMOS spectrum, which shows highly-ionized emission lines (Ovi, Nv, and Niv). The line properties have been compared with those observed in galaxies and AGN, and with the predictions from radiation transfer models for star-forming galaxies, AGN, and shocks. The SED-fitting provides a total stellar mass of $M_*=(4 \pm 1) \times 10^{10} M_\odot$, an age of the main stellar population of $\sim670 Myr$ and a recent short (8 Myr) burst of star formation (SF) of $\sim 90 M_\odot yr^{-1}$. The Nv line has a characteristic P-Cygni profile, which suggests a $\sim 3 Myr$ old population of stars with a mass of $\sim 5 \times10^{7} M_\odot$. The Nv profile also shows evidence for an additional component of nebular emission. The comparison of the line ratios with theoretical models allows us to associate the emission with SF or AGN, but the strong radiation field required to ionize the Ovi is more commonly related to AGN activity. We found evidence for an old and already evolved stellar population at $z\sim 5.5$ and show that the galaxy is experiencing a second short burst of SF. In addition, GS-14 carries signatures of obscured AGN activity. The AGN could be responsible for the short depletion time of this galaxy, thus making GS-14 one of the two ALPINE sources with hints of an active nucleus and an interesting target for future follow-ups. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: Accepted in A&A, 10 pages, 4 figures

Journal ref: A&A 675, A30 (2023)

ALMA Observations for CO Emission from Luminous Lyman-break Galaxies at $z=6.0293$-$6.2037$

Authors: Yoshiaki Ono, Seiji Fujimoto, Yuichi Harikane, Masami Ouchi, Livia Vallini, Andrea Ferrara, Takatoshi Shibuya, Andrea Pallottini, Akio K. Inoue, Masatoshi Imanishi, Kazuhiro Shimasaku, Takuya Hashimoto, Chien-Hsiu Lee, Yuma Sugahara, Yoichi Tamura, Kotaro Kohno, Malte Schramm

Abstract: We present our new Atacama Large Millimeter/submillimeter Array (ALMA) observations targeting CO(6-5) emission from three luminous Lyman break galaxies (LBGs) at $z_{\rm spec} = 6.0293$-$6.2037$ found in the Subaru/Hyper Suprime-Cam survey, whose [OIII]$88μ$m and [CII]$158μ$m emission have been detected with ALMA. We find a marginal detection of the CO(6-5) line from one of our LBGs, J0235-0532, a… ▽ More We present our new Atacama Large Millimeter/submillimeter Array (ALMA) observations targeting CO(6-5) emission from three luminous Lyman break galaxies (LBGs) at $z_{\rm spec} = 6.0293$-$6.2037$ found in the Subaru/Hyper Suprime-Cam survey, whose [OIII]$88μ$m and [CII]$158μ$m emission have been detected with ALMA. We find a marginal detection of the CO(6-5) line from one of our LBGs, J0235-0532, at the $\simeq 4 σ$ significance level and obtain upper limits for the other two LBGs, J1211-0118 and J0217-0208. Our $z=6$ luminous LBGs are consistent with the previously found correlation between the CO luminosity and the infrared luminosity. The unique ensemble of the multiple far-infrared emission lines and underlying continuum fed to a photodissociation region model reveal that J0235-0532 has a relatively high hydrogen nucleus density that is comparable to those of low-$z$ (U)LIRGs, quasars, and Galactic star-forming regions with high $n_{\rm H}$ values, while the other two LBGs have lower $n_{\rm H}$ consistent with local star-forming galaxies. By carefully taking account of various uncertainties, we obtain total gas mass and gas surface density constraints from their CO luminosity measurements. We find that J0235-0532 locates below the Kennicutt-Schmidt (KS) relation, comparable to the previously CO(2-1) detected $z=5.7$ LBG, HZ10. Combined with previous results for dusty starbursts at similar redshifts, the KS relation at $z=5$-$6$ is on average consistent with the local one. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: 33 pages, 12 figures, accepted for publication in ApJ

Spritz is sparkling: simulated CO and [CII] luminosities

Authors: L. Bisigello, L. Vallini, C. Gruppioni, F. Esposito, F. Calura, I. Delvecchio, A. Feltre, F. Pozzi, G. Rodighiero

Abstract: We present a new prediction of the luminosity functions of the [CII] line at 158 $μ$m, of the CO lines from J=0 to J=24, and of the molecular gas mass density up to z=10, using the Spectro-Photometric Realisations of Infrared-selected Targets at all-z (SPRITZ) simulation (Bisigello et al. 2021). We update the state-of-the-art phenomenological simulation SPRITZ to include both the CO ($J\leq24$) an… ▽ More We present a new prediction of the luminosity functions of the [CII] line at 158 $μ$m, of the CO lines from J=0 to J=24, and of the molecular gas mass density up to z=10, using the Spectro-Photometric Realisations of Infrared-selected Targets at all-z (SPRITZ) simulation (Bisigello et al. 2021). We update the state-of-the-art phenomenological simulation SPRITZ to include both the CO ($J\leq24$) and [CII] line luminosities. This has been performed using different empirical and theoretical relations to convert the total infrared luminosity (or star formation rate) to [CII] or CO luminosity. The resulting line luminosity functions have been compared for validation with a large set of observations available in the literature. We then used the derived CO and [CII] line luminosities to estimate the molecular gas mass density and compare it with available observations. The CO and [CII] luminosity functions presented here are well in agreement with all the available observations. In particular, the best results for [CII] are obtained deriving the [CII] luminosity directly from the star formation rate, but considering a dependence of this relation on the gas metallicity. For all the CO luminosity functions, the estimates favoured by the data are derived considering different relations, depending on the ionisation mechanism dominating each galaxy, i.e. star formation or active galactic nuclei, and, moreover, deriving the $J\geq4$ CO lines directly from the [CII] luminosity. However, further data are necessary to fully discriminate between models. Finally, the best agreement with observations of the molecular gas mass density are derived by converting the [CII] luminosity to H2 mass, using a [CII]-to-H2 conversion ~130 $\rm M_{\odot}/{\rm L}_{\odot}$. All the line luminosity functions, useful for planning and interpreting future observations, are made publicly available. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 18 pages, 16 figures, 8 tables, accepted for publication in A&A

Journal ref: A&A 666, A193 (2022)

The black hole and host galaxy growth in an isolated $z\sim 6$ QSO observed with ALMA

Authors: R. Tripodi, C. Feruglio, F. Fiore, M. Bischetti, V. D'Odorico, S. Carniani, S. Cristiani, S. Gallerani, R. Maiolino, A. Marconi, A. Pallottini, E. Piconcelli, L. Vallini, T. Zana

Abstract: 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 (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass (… ▽ More 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 (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass ($M_{\rm dyn}$) plane, we compare the star formation rate (SFR), derived from the accurate estimate of the dust temperature and the dust mass ($T_{\rm dust}, M_{\rm dust}$), with the BH accretion rate. To this aim, we analysed a deep, $900$ pc resolution ALMA observation of the sub-mm continuum, [CII] and H$_2$O of the $z\sim 6$ QSO J2310+1855, enabling a detailed study of dust properties and cold gas kinematics. We performed an accurate SED analysis obtaining a dust temperature of $T_{\rm dust} = 71$ K and a dust mass of $M_{\rm dust}= 4.4 \times 10^8\ \rm M_{\odot}$. The implied AGN-corrected SFR is $1240 \ \rm M_{\odot}yr^{-1}$, a factor of 2 smaller than previously reported for this QSO. We derived the best estimate of the dynamical mass $M_{\rm dyn} = 5.2\times 10^{10}\ \rm M_{\odot}$ within $r = 1.7$ kpc, based on a dynamical model of the system. We found that ${\rm SFR}/M_{\rm dyn}>\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. △ Less

Submitted 30 July, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 15 pages, 10 figures, 4 tables. Accepted in A&A

Journal ref: A&A 665, A107 (2022)

The interstellar medium of high-redshift galaxies: Gathering clues from C III] and [C II] lines

Authors: Vladan Markov, Stefano Carniani, Livia Vallini, Andrea Ferrara, Andrea Pallottini, Roberto Maiolino, Simona Gallerani, Laura Pentericci

Abstract: A tight relation between [C II] line luminosity and the star formation rate (SFR) has been observed for local galaxies. At high redshift (z > 5), galaxies instead deviate downwards from the local $Σ$_[C II] - $Σ$_SFR relation. This deviation might be caused by different interstellar medium (ISM) properties in galaxies at early epochs. To test this hypothesis, we combined the [C II] and SFR data wi… ▽ More A tight relation between [C II] line luminosity and the star formation rate (SFR) has been observed for local galaxies. At high redshift (z > 5), galaxies instead deviate downwards from the local $Σ$_[C II] - $Σ$_SFR relation. This deviation might be caused by different interstellar medium (ISM) properties in galaxies at early epochs. To test this hypothesis, we combined the [C II] and SFR data with C III] line observations and our physical models. We additionally investigated how ISM properties, such as burstiness, $κ_s$, total gas density, $n$, and metallicity, $Z$, affect the deviation from the $Σ$_[C II] - $Σ$_SFR relation in these sources. We present the VLT/X-SHOOTER observations targeting the C III] $λ1909$ line emission in three galaxies at 5.5 < z < 7. We include X-SHOOTER and VLT/MUSE archival data of eight galaxies at 2 < z < 7, and eleven star-forming systems at 6 < z < 7.5, with either C III] or [C II] detection reported in the literature. We detected C III] $λλ1907, 1909$ line emission in HZ10 and we derived the intrinsic, integrated flux of the C III] $λ1909$ line. We constrained the ISM properties for our sample of galaxies, $κ_s$ , $n$, and $Z$, by applying our physically motivated model based on the MCMC algorithm. For the most part, high-z star-forming galaxies show subsolar metallicities. The majority of the sources have $log(κ_s) > 1$, that is, they overshoot the Kennicutt-Schmidt (KS) relation by about one order of magnitude, implying that the whole KS relation might be shifted upwards at early times. Furthermore, all the high-z galaxies of our sample lie below the $Σ$_[C II] - $Σ$_SFR local relation. The total gas density, $n$, shows the strongest correlation with the deviation from the local $Σ$_[C II] - $Σ$_SFR relation. △ Less

Submitted 17 June, 2022; v1 submitted 7 June, 2022; originally announced June 2022.

Comments: 16 pages, 13 figures, to be published in A&A

Journal ref: A&A 663, A172 (2022)

CONCERTO: High-fidelity simulation of millimeter line emissions of galaxies and [CII] intensity mapping

Authors: M. Bethermin, A. Gkogkou, M. Van Cuyck, G. Lagache, A. Beelen, M. Aravena, A. Benoit, J. Bounmy, M. Calvo, A. Catalano, B. de Batz de Trenquelleon, C. De Breuck, A. Fasano, A. Ferrara, J. Goupy, C. Hoarau, C. Horellou, W. Hu, A. Julia, K. Knudsen, J. -C. Lambert, J. Macias-Perez, J. Marpaud, A. Monfardini, A. Pallottini , et al. (5 additional authors not shown)

Abstract: 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-m… ▽ More 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. △ Less

Submitted 24 August, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: 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

Journal ref: A&A 667, A156 (2022)

Dynamical characterization of galaxies up to $z \sim 7$

Authors: F. Rizzo, M. Kohandel, A. Pallottini, A. Zanella, A. Ferrara, L. Vallini, S. Toft

Abstract: The characterization of the dynamical state of galaxies up to z~7 is crucial for constraining the mechanisms driving the mass assembly in the early Universe. However, it is unclear whether the data quality of current and future observations is sufficient to perform a solid dynamical analysis. This paper defines the angular resolution and S/N required for a robust characterization of the dynamical… ▽ More The characterization of the dynamical state of galaxies up to z~7 is crucial for constraining the mechanisms driving the mass assembly in the early Universe. However, it is unclear whether the data quality of current and future observations is sufficient to perform a solid dynamical analysis. This paper defines the angular resolution and S/N required for a robust characterization of the dynamical state of galaxies up to the EoR. The final aim is to help design spatially-resolved surveys targeting emission lines of primeval galaxies. We investigate the [CII]-158um emission from z~6-7 LBGs from the SERRA cosmological simulation, covering a range of dynamical states: from disks to major mergers. We create ALMA mock observations with various data quality and apply the kinematic classification methods used in the literature. These tests allow us to quantify the performances of such methods as a function of angular resolution and S/N. We find that barely-resolved observations do not allow the correct dynamical characterization of a galaxy, resulting in the misclassification of all disks in our sample. However, even when using spatially-resolved observations with data quality typical of high-z galaxies, the standard kinematic classification methods, based on the analysis of the moment maps, fail to distinguish a merger from a disk. The high angular resolution and S/N needed to apply these standard methods successfully can be achieved with current data only for a handful of bright galaxies. We propose a new classification method, called PVsplit, that quantifies the asymmetries and morphological features in position-velocity diagrams using three empirical parameters. We test PVsplit on our mock data concluding that it can predict whether a galaxy is a disk or a merger provided that S/N $\gtrsim10$, and the major axis is covered by $\gtrsim3$ independent resolution elements. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: Submitted to Astronomy and Astrophysics (A&A) Journal. Comments are welcome

Journal ref: A&A 667, A5 (2022)

Photodissociation and X-Ray Dominated Regions

Authors: Mark G. Wolfire, Livia Vallini, Mélanie Chevance

Abstract: The radiation from stars and active galactic nuclei (AGN) creates photodissociation regions (PDRs) and X-ray dominated regions (XDRs), where the chemistry or heating are dominated by far-ultraviolet (FUV) radiation or X-ray radiation, respectively. PDRs include a wide range of environments from the diffuse interstellar medium to dense star-forming regions. XDRs are found in the center of galaxies… ▽ More The radiation from stars and active galactic nuclei (AGN) creates photodissociation regions (PDRs) and X-ray dominated regions (XDRs), where the chemistry or heating are dominated by far-ultraviolet (FUV) radiation or X-ray radiation, respectively. PDRs include a wide range of environments from the diffuse interstellar medium to dense star-forming regions. XDRs are found in the center of galaxies hosting AGN, in protostellar disks, and in the vicinity of X-ray binaries. In this review, we describe the dominant thermal, chemical, and radiation transfer processes in PDRs and XDRs, as well as a brief description of models and their use to analyze observations. We then present recent results from Milky Way, nearby extragalactic, and high-redshift observations. Several important results are: $\bullet$ Velocity resolved PDR lines reveal the kinematics of the neutral atomic gas and provide constraints on the stellar feedback process. Their interpretation is, however, in dispute as observations suggest a prominent role for stellar winds while they are much less important in theoretical models. $\bullet$ A significant fraction of molecular mass resides in CO-dark gas especially in low-metallicity/highly irradiated environments. $\bullet$ The CO ladder and CI/CII ratios can determine if FUV or X-rays dominate the ISM heating of extragalactic sources. $\bullet$ With ALMA, PDR and XDR tracers are now routinely detected on galactic scales over cosmic time. This makes it possible to link the star formation history of the Universe to the evolution of the physical and chemical properties of the gas. △ Less

Submitted 15 February, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: Submitted to Annual Reviews of Astronomy and Astrophysics, Volume 60, 2022

AGN impact on the molecular gas in galactic centers as probed by CO lines

Authors: Federico Esposito, Livia Vallini, Francesca Pozzi, Viviana Casasola, Matilde Mingozzi, Cristian Vignali, Carlotta Gruppioni, Francesco Salvestrini

Abstract: We present a detailed analysis of the X-ray, infrared, and carbon monoxide (CO) emission for a sample of 35 local ($z \leq 0.15$), active ($L_X \geq 10^{42}$ erg s$^{-1}$) galaxies. Our goal is to infer the contribution of far-ultraviolet (FUV) radiation from star formation (SF), and X-ray radiation from the active galactic nuclei (AGN), respectively producing photodissociation regions (PDRs) and… ▽ More We present a detailed analysis of the X-ray, infrared, and carbon monoxide (CO) emission for a sample of 35 local ($z \leq 0.15$), active ($L_X \geq 10^{42}$ erg s$^{-1}$) galaxies. Our goal is to infer the contribution of far-ultraviolet (FUV) radiation from star formation (SF), and X-ray radiation from the active galactic nuclei (AGN), respectively producing photodissociation regions (PDRs) and X-ray dominated regions (XDRs), to the molecular gas heating. To this aim, we exploit the CO spectral line energy distribution (CO SLED) as traced by Herschel, complemented with data from single-dish telescopes for the low-J lines, and high-resolution ALMA images of the mid-J CO emitting region. By comparing our results to the Schmidt-Kennicutt relation, we find no evidence for AGN influence on the cold and low-density gas on kpc-scales. On nuclear (r = 250 pc) scales, we find weak correlations between the CO line ratios and either the FUV or X-ray fluxes: this may indicate that neither SF nor AGN radiation dominates the gas excitation, at least at r = 250 pc. From a comparison of the CO line ratios with PDR and XDR models, we find that PDRs can reproduce observations only in presence of extremely high gas densities ($n > 10^5$ cm$^{-3}$). In the XDR case, instead, the models suggest moderate densities ($n \approx 10^{2-4}$ cm$^{-3}$). We conclude that a mix of the two mechanisms (PDR for the mid-J, XDR or possibly shocks for the high-J) is necessary to explain the observed CO excitation in active galaxies. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: 25 pages, 19 figures, accepted for publication in MNRAS

A survey of high-$z$ galaxies: SERRA simulations

Authors: A. Pallottini, A. Ferrara, S. Gallerani, C. Behrens, M. Kohandel, S. Carniani, L. Vallini, S. Salvadori, V. Gelli, L. Sommovigo, V. D'Odorico, F. Di Mascia, E. Pizzati

Abstract: We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/ga… ▽ More We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/gas/dust mass, metallicity) of high-redshift $6 \lesssim z \lesssim 15$ galaxies. This flagship paper focuses on the $z=7.7$ sub-sample, including 202 galaxies with stellar mass $10^7 M_\odot \lesssim M_\star \lesssim 5\times 10^{10}M_\odot$, and specific star formation ranging from ${\rm sSFR} \sim 100\,{\rm Gyr}^{-1}$ in young, low-mass galaxies to $\sim 10\,{\rm Gyr}^{-1}$ for older, massive ones. At this redshift, SERRA galaxies are typically bursty, i.e. they are located above the Schmidt-Kennicutt relation by a factor $κ_s = 3.03^{+4.9}_{-1.8}$, consistent with recent findings for [OIII] and [CII] emitters at high-$z$. They also show relatively large ${\rm IRX} = L_{\rm FIR}/L_{\rm UV}$ values as a result of their compact/clumpy morphology effectively blocking the stellar UV luminosity. Note that this conclusion might be affected by insufficient spatial resolution at the molecular cloud level. We confirm that early galaxies lie on the standard $\rm [CII]-SFR$ relation; their observed $L_{\rm [OIII]}/L_{\rm [CII]} \simeq 1-10$ ratios can be reproduced by a part of the SERRA galaxies without the need of a top-heavy IMF and/or anomalous C/O abundances. [OI] line intensities are similar to local ones, making ALMA high-$z$ detections challenging but feasible ($\sim 6\,\rm hr$ for a SFR of $50\,M_\odot\,{\rm yr}^{-1}$). △ Less

Submitted 4 May, 2022; v1 submitted 7 January, 2022; originally announced January 2022.

Comments: 23 pages, 14 figures, accepted by MNRAS; part of the data used for this study is available at the website http://cosmology.sns.it/data_access.html

High [OIII]/[CII] surface brightness ratios trace early starburst galaxies

Authors: L. Vallini, A. Ferrara, A. Pallottini, S. Carniani, S. Gallerani

Abstract: We study the impact of deviations from the Kennicutt-Schmidt relation (quantified by the `burstiness' parameter $κ_s$), gas metallicity ($Z$), and density ($n$) on the observed [OIII]88$μ$m/[CII]158$μ$m surface brightness ratios ($Σ_{[OIII]}/Σ_{[CII]}$) in nine galaxies at $z\approx6-9$. We first discuss possible biases in the measured $Σ_{[OIII]}/Σ_{[CII]}$ ratios by comparing the data with zoom-… ▽ More We study the impact of deviations from the Kennicutt-Schmidt relation (quantified by the `burstiness' parameter $κ_s$), gas metallicity ($Z$), and density ($n$) on the observed [OIII]88$μ$m/[CII]158$μ$m surface brightness ratios ($Σ_{[OIII]}/Σ_{[CII]}$) in nine galaxies at $z\approx6-9$. We first discuss possible biases in the measured $Σ_{[OIII]}/Σ_{[CII]}$ ratios by comparing the data with zoom-in cosmological simulations, and then use a Markov Chain Monte Carlo algorithm to derive the best fit values of ($κ_s, Z, n$). We find that (i) the strongest dependence of $Σ_{[OIII]}/Σ_{[CII]}$ is on $κ_s$; (ii) high ratios identify starburst galaxies with short gas depletion times ($t_{dep}=6-49\,\rm Myr$); (iii) a secondary dependence on density is found, with $Σ_{[OIII]}/Σ_{[CII]}$ anticorrelating with $n$ as a result of the lower [OIII] critical density, (iv) the ratio only weakly depends on $Z$. The nine galaxies are significantly enriched (Z=0.2-0.5 $Z_\odot$), and dense ($n=10^{1-3} {\rm cm}^{-3}$). This lends further support to the starburst scenario in which a rapid enrichment of the interstellar medium is expected. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: accepted for publication in MNRAS, 12 pages, 8 figures

The ALPINE-ALMA [CII] survey: Dust mass budget in the early Universe

Authors: F. Pozzi, F. Calura, Y. Fudamoto, M. Dessauges-Zavadsky, C. Gruppioni, M. Talia, G. Zamorani, M. Bethermin, A. Cimatti, A. Enia, Y. Khusanova, R. Decarli, O. Le Fevre, P. Capak, P. Cassata, A. L. Faisst, L. Yan, D. Schaerer, J. Silverman, S. Bardelli, M. Boquien, A. Enia, D. Narayanan, M. Ginolfi, N. P. Hathi , et al. (11 additional authors not shown)

Abstract: The dust content of normal galaxies and the dust mass density (DMD) at high-z (z>4) are unconstrained given the source confusion and the sensitivity limitations of previous observations. The ALMA Large Program to INvestigate [CII] at Early Times (ALPINE), which targeted 118 UV-selected star-forming galaxies at 4.4<z<5.9, provides a new opportunity to tackle this issue for the first time with a sta… ▽ More The dust content of normal galaxies and the dust mass density (DMD) at high-z (z>4) are unconstrained given the source confusion and the sensitivity limitations of previous observations. The ALMA Large Program to INvestigate [CII] at Early Times (ALPINE), which targeted 118 UV-selected star-forming galaxies at 4.4<z<5.9, provides a new opportunity to tackle this issue for the first time with a statistically robust dataset. We have exploited the rest-frame far-infrared (FIR) fluxes of the 23 continuum individually detected galaxies and stacks of continuum images to measure the dust content of the 118 UV-selected ALPINE galaxies. We have focused on the dust scaling relations and, by comparing them with predictions from chemical evolution models, we have probed the evolutionary stage of UV-selected galaxies at high-z. By using the observed correlation between the UV-luminosity and the dust mass, we have estimated the DMD of UV-selected galaxies at z~5, weighting the galaxies by means of the UV-luminosity function (UVLF). The derived DMD has been compared with the value we have estimated from the 10 ALPINE far-IR continuum blindly detected galaxies at the redshift of the ALPINE targets. The comparison of the observed dust scaling relations with chemical evolution models suggests that ALPINE galaxies are not likely progenitors of disc galaxies, but of intermediate and low mass proto-spheroids, resulting in present-day bulges of spiral or elliptical galaxies. Interestingly, this conclusion is in line with the independent morphological analysis, that shows that the majority (~70\%) of the dust-continuum detected galaxies have a disturbed morphology. The DMD obtained at z~5 from UV-selected sources is ~30% of the value obtained from blind far-IR selected sources, showing that the UV-selection misses the most dust-rich, UV-obscured galaxies. △ Less

Submitted 31 May, 2021; originally announced May 2021.

Comments: 16 pages, 6 figures. Accepted for publication in A&A

Journal ref: A&A 653, A84 (2021)

Dust temperature in ALMA $\hbox{[C $\scriptstyle\rm II $]}$-detected high-$z$ galaxies

Authors: L. Sommovigo, A. Ferrara, S. Carniani, A. Zanella, A. Pallottini, S. Gallerani, L. Vallini

Abstract: 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 qua… ▽ More 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μ$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μ$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 $Σ_{\rm H_2} / Σ_{\rm gas}$. We find that $Σ_{\rm H_2} / Σ_{\rm gas}$ slowly increases with redshift, which is compatible with current observations at $0 < z < 4$. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 14 pages, 5 figures, submitted to MNRAS

Velocity dispersion in the interstellar medium of early galaxies

Authors: M. Kohandel, A. Pallottini, A. Ferrara, S. Carniani, S. Gallerani, L. Vallini, A. Zanella, C. Behrens

Abstract: We study the structure of spatially resolved, line-of-sight velocity dispersion for galaxies in the Epoch of Reionization (EoR) traced by [CII] $158μ\rm{m}$ line emission. Our laboratory is a simulated prototypical Lyman-break galaxy, "Freesia", part of the SERRA suite. The analysis encompasses the redshift range 6 < z < 8, when Freesia is in a very active assembling phase. We build velocity dispe… ▽ More We study the structure of spatially resolved, line-of-sight velocity dispersion for galaxies in the Epoch of Reionization (EoR) traced by [CII] $158μ\rm{m}$ line emission. Our laboratory is a simulated prototypical Lyman-break galaxy, "Freesia", part of the SERRA suite. The analysis encompasses the redshift range 6 < z < 8, when Freesia is in a very active assembling phase. We build velocity dispersion maps for three dynamically distinct evolutionary stages (Spiral Disk at z=7.4, Merger at z=8.0, and Disturbed Disk at z=6.5) using [CII] hyperspectral data cubes. We find that, at a high spatial resolution of 0.005" ($\simeq 30 pc$), the luminosity-weighted average velocity dispersion is $σ_{\rm{CII}}$~23-38 km/s with the highest value belonging to the highly-structured Disturbed Disk stage. Low resolution observations tend to overestimate $σ_{\rm CII}$ values due to beam smearing effects that depend on the specific galaxy structure. For an angular resolution of 0.02" (0.1"), the average velocity dispersion is 16-34% (52-115%) larger than the actual one. The [CII] emitting gas in Freesia has a Toomre parameter $\mathcal{Q}$~0.2 and a rotational-to-dispersion ratio of $v_{\rm c}/σ$~ 7 similar to that observed in z=2-3 galaxies. The primary energy source for the velocity dispersion is due to gravitational processes, such as merging/accretion events; energy input from stellar feedback is generally subdominant (< 10%). Finally, we find that the resolved $σ_{\rm{CII}} - Σ_{\rm SFR}$ relation is relatively flat for $0.02<Σ_{\rm SFR}/{{\rm M}_{\odot}} \mathrm{yr}^{-1} {\mathrm kpc}^{-2} < 30$, with the majority of data lying on the derived analytical relation $σ\propto Σ_{\rm SFR}^{5/7}$. At high SFR, the increased contribution from stellar feedback steepens the relation, and $σ_{\rm{CII}}$ rises slightly. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS; 16 pages (plus appendix), 8 figures, 3 tables

Missing [CII] emission from early galaxies

Authors: S. Carniani, A. Ferrara, R. Maiolino, M. Castellano, S. Gallerani, A. Fontana, M. Kohandel, A. Lupi, A. Pallottini, L. Pentericci, L. Vallini, E. Vanzella

Abstract: ALMA observations have revealed that [CII] 158$μ$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of n… ▽ More ALMA observations have revealed that [CII] 158$μ$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of nine z>6 galaxies observed in both [CII] and [OIII]. After performing several uv-tapering experiments to optimize the identification of extended line emission, we detect [CII] emission in the whole sample, with an extent systematically larger than the [CII] emission. Next, we use interferometric simulations to study the effect of SBD on the line luminosity estimate. About 40% of the extended [CII] component might be missed at an angular resolution of 0.8$^{\prime\prime}$, implying that $L_{\rm [CII]}$ is underestimated by a factor $\approx2$ in data at low (<7) signal-to-noise ratio . By combining these results, we conclude that $L_{\rm [CII]}$ of z>6 galaxies lies, on average, slightly below the local $L_{\rm [CII]}-SFR$ relation ($Δ^{z=6-9}=-0.07\pm0.3$), but within the intrinsic dispersion of the relation. SBD correction also yields $L_{\rm [OIII]}/L_{\rm [CII]}<10$, i.e. more in line with current hydrodynamical simulations. △ Less

Submitted 7 December, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

The ALPINE-ALMA [CII] Survey: The nature, luminosity function and star formation history of dusty galaxies up to z~6

Authors: C. Gruppioni, M. Bethermin, F. Loiacono, O. Le Fevre, P. Capak, P. Cassata, A. L. Faisst, D. Schaerer, J. Silverman, L. Yan, S. Bardelli, M. Boquien, R. Carraro, A. Cimatti, M. Dessauges-Zavadsky, M. Ginolfi, S. Fujimoto, N. P. Hathi, G. C. Jones, Y. Khusanova, A. M. Koekemoer, G. Lagache, B. C. Lemaux, P. Oesch, F. Pozzi , et al. (8 additional authors not shown)

Abstract: We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all… ▽ More We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 14 per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (six show a counterpart only in the mid-IR and no HST or near-IR identification, while two are detected as [CII] emitters at z=5). The six HST and near-IR dark galaxies with mid-IR counterpart contribute for about 17 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3. △ Less

Submitted 30 September, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 26 pages, 18 figures. Accepted for publication in A&A

Journal ref: A&A 643, A8 (2020)

The ALPINE-ALMA [C II] survey: the luminosity function of serendipitous [C II] line emitters at $z\sim 5$

Authors: Federica Loiacono, Roberto Decarli, Carlotta Gruppioni, Margherita Talia, Andrea Cimatti, Gianni Zamorani, Francesca Pozzi, Lin Yan, Brian C. Lemaux, Dominik A. Riechers, Olivier Le Fèvre, Matthieu Béthermin, Peter Capak, Paolo Cassata, Andreas Faisst, Daniel Schaerer, John D. Silverman, Sandro Bardelli, Médéric Boquien, Sandra Burkutean, Miroslava Dessauges-Zavadsky, Yoshinobu Fudamoto, Seiji Fujimoto, Michele Ginolfi, Nimish P. Hathi , et al. (10 additional authors not shown)

Abstract: We present the first [CII] 158 $μ$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we… ▽ More We present the first [CII] 158 $μ$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we identified 8 out of 14 detections as [CII] lines at $z\sim 5$, and two as CO transitions at lower redshifts. The remaining 4 lines have an elusive identification in the available catalogs and we considered them as [CII] candidates. We used the 8 confirmed [CII] and the 4 [CII] candidates to build one of the first [CII] LFs at $z\sim 5$. We found that 11 out of these 12 sources have a redshift very similar to that of the ALPINE target in the same pointing, suggesting the presence of overdensities around the targets. Therefore, we split the sample in two (a "clustered" and "field" sub-sample) according to their redshift separation and built two separate LFs. Our estimates suggest that there could be an evolution of the [CII] LF between $z \sim 5$ and $z \sim 0$. By converting the [CII] luminosity to star formation rate we evaluated the cosmic star formation rate density (SFRD) at $z\sim 5$. The clustered sample results in a SFRD $\sim 10$ times higher than previous measurements from UV-selected galaxies. On the other hand, from the field sample (likely representing the average galaxy population) we derived a SFRD $\sim 1.6$ higher compared to current estimates from UV surveys but compatible within the errors. Because of the large uncertainties, observations of larger samples are necessary to better constrain the SFRD at $z\sim 5$. This study represents one of the first efforts aimed at characterizing the demography of [CII] emitters at $z\sim 5$ using a mm-selection of galaxies. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: 19 pages, 12 figures; submitted to Astronomy & Astrophysics

Journal ref: A&A 646, A76 (2021)

The ALPINE-ALMA [C II] Survey: [C II]158micron Emission Line Luminosity Functions at $z \sim 4-6$

Authors: Lin Yan, A. Sajina, F. Loiacono, G. Lagache, M. Bèthermin, A. Faisst, M. Ginolfi, O. Le Fèvre, C. Gruppioni, P. L. Capak, P. Cassata, D. Schaerer, J. D. Silverman, S. Bardelli, M. Dessauges-Zavadsky, A. Cimatti, N. P. Hathi, B. C. Lemaux, E. Ibar, G. C. Jones, A. M. Koekemoer, P. A. Oesch, M. Talia, F. Pozzi, D. A. Riechers , et al. (6 additional authors not shown)

Abstract: We present the [CII]158$μ$m line luminosity functions (LFs) at $z\sim4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A}<-20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [CII] detections and 43 are upper limits. This is by far the largest sample of [CII] detections which allows us to set constr… ▽ More We present the [CII]158$μ$m line luminosity functions (LFs) at $z\sim4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A}<-20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [CII] detections and 43 are upper limits. This is by far the largest sample of [CII] detections which allows us to set constraints to the volume density of [CII] emitters at $z\sim4-6$. But because this is a UV-selected sample, we are missing [CII]-bright but UV-faint sources making our constraints strict lower limits. Our derived LFs are statistically consistent with the $z\sim0$ [CII] LF at $10^{8.25} - 10^{9.75}L_\odot$. We compare our results with the upper limits of the [CII] LF derived from serendipitous sources in the ALPINE maps (Loiacono et al. 2020). We also infer the [CII] LFs based on published far-IR and CO LFs at $z\sim4-6$. Combining our robust lower limits with these additional estimates, we set further constraints to the true number density of [CII] emitters at $z\sim 4 - 6$. These additional LF estimates are largely above our LF at $L_{[CII]}>10^9L_{\odot}$, suggesting that UV-faint but [CII]-bright sources likely make a significant contributions to the [CII] emitter volume density. When we include all the LF estimates, we find that available model predictions underestimate the number densities of [CII] emitters at $z\sim4-6$. Finally, we set a constraint on the molecular gas mass density at $z\sim4-6$, with $ρ_{mol} \sim (2-7)\times10^7M_\odot$\,Mpc$^{-3}$. This is broadly consistent with previous studies. △ Less

Submitted 29 October, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 17 pages, 9 Figures, Update to match with the published version. Accepted for the publication in ApJ

The ALPINE-ALMA [CII] survey: Molecular gas budget in the Early Universe as traced by [CII]

Authors: M. Dessauges-Zavadsky, M. Ginolfi, F. Pozzi, M. Béthermin, O. Le Fèvre, S. Fujimoto, J. D. Silverman, G. C. Jones, D. Schaerer, A. L. Faisst, Y. Khusanova, Y. Fudamoto, P. Cassata, F. Loiacono, P. L. Capak, L. Yan, R. Amorin, S. Bardelli, M. Boquien, A. Cimatti, C. Gruppioni, N. P. Hathi, E. Ibar, A. M. Koekemoer, B. C. Lemaux , et al. (10 additional authors not shown)

Abstract: The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Lar… ▽ More The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. The agreement between molecular gas masses derived from [CII] luminosity, dynamical mass, and rest-frame 850um luminosity, extrapolated from the rest-frame 158um continuum, supports [CII] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z=0 to z=5.9, which reaches a mean value of (4.6+/-0.8)x10^8 yr at z~5.5, only a factor of 2-3 shorter than in present-day galaxies. This suggests a mild enhancement of star formation efficiency toward high redshifts, unless the molecular gas fraction significantly increases. Our estimates show that the rise in molecular gas fraction as reported previously, flattens off above z~3.7 to achieve a mean value of 63%+/-3 over z=4.4-5.9. This redshift evolution of the gas fraction is in line with the one of the specific star formation rate. We use multi-epoch abundance matching to follow the gas fraction evolution over cosmic time of progenitors of z=0 Milky Way-like galaxies in 10^13 Msun halos and of more massive z=0 galaxies in 10^14 Msun halos. Interestingly, the former progenitors show a monotonic decrease of the gas fraction with cosmic time, while the latter show a constant gas fraction from z=5.9 to z~2 and a decrease at z<2. We discuss three possible effects, namely outflows, halt of gas supplying, and over-efficient star formation, which may jointly contribute to the gas fraction plateau of the latter massive galaxies. △ Less

Submitted 22 April, 2020; originally announced April 2020.

Comments: 16 pages, 9 figures, 2 tables. Submitted to Astronomy & Astrophysics

Journal ref: A&A 643, A5 (2020)

The ALPINE-ALMA [CII] survey: Dust attenuation properties and obscured star formation at z~4.4-5.8

Authors: Yoshinobu Fudamoto, P. A. Oesch, A. Faisst, M. Bethermin, M. Ginolfi, Y. Khusanova, F. Loiacono, O. Le Fevre, P. Capak, D. Schaerer, J. D. Silverman, P. Cassata, L. Yan, R. Amorin, S. Bardelli, M. Boquien, A. Cimatti, M. Dessauges-Zavadsky, S. Fujimoto, C. Gruppioni, N. P. Hathi, E. Ibar, G. C. Jones, A. M. Koekemoer, G. Lagache , et al. (12 additional authors not shown)

Abstract: We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at a redshift of ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame 158μm obtained with the Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE). We study the connection between the ultrav… ▽ More We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at a redshift of ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame 158μm obtained with the Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE). We study the connection between the ultraviolet (UV) spectral slope ($β$), stellar mass (M_*), and infrared excess (IRX=L_IR/L_UV). Twenty-three galaxies are individually detected in the continuum at >3.5 sigma significance. We perform a stacking analysis using both detections and nondetections to study the average dust attenuation properties at z~4.4-5.8. The individual detections and stacks show that the IRX-$β$ relation at z~5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts (z<4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud (SMC). This systematic change of the IRX-$β$ relation as a function of redshift suggests an evolution of dust attenuation properties at z>4. Similarly, we find that our galaxies have lower IRX values, up to 1dex on average, at a fixed mass compared to previously studied IRX-M_* relations at z<4, albeit with significant scatter. This implies a lower obscured fraction of star formation than at lower redshifts. Our results suggest that dust properties of UV-selected star forming galaxies at z>4 are characterised by (i) a steeper attenuation curve than at z<4, and (ii) a rapidly decreasing dust obscured fraction of star formation as a function of redshift. Nevertheless, even among this UV-selected sample, massive galaxies (log M_*/$M_\odot$>10) at z~5-6 already exhibit an obscured fraction of star formation of ~45%, indicating a rapid build-up of dust during the epoch of reionization. △ Less

Submitted 6 October, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

Comments: 12 pages, 6 figures. accepted for publication in Astronomy & Astrophysics. Minor revisions

Journal ref: A&A 643, A4 (2020)

Predicting FIR lines from simulated galaxies

Authors: Alessandro Lupi, Andrea Pallottini, Andrea Ferrara, Stefano Bovino, Stefano Carniani, Livia Vallini

Abstract: Far-infrared (FIR) emission lines are a powerful tool to investigate the properties of the interstellar medium, especially in high-redshift galaxies, where ALMA observations have provided unprecedented information. Interpreting such data with state-of-the-art cosmological simulations post-processed with CLOUDY, has provided insights on the internal structure and gas dynamics of these systems. Howe… ▽ More Far-infrared (FIR) emission lines are a powerful tool to investigate the properties of the interstellar medium, especially in high-redshift galaxies, where ALMA observations have provided unprecedented information. Interpreting such data with state-of-the-art cosmological simulations post-processed with CLOUDY, has provided insights on the internal structure and gas dynamics of these systems. However, no detailed investigation of the consistency and uncertainties of this kind of analysis has been performed to date. Here, we compare different approaches to estimate FIR line emission from state-of-the-art cosmological simulations, either with CLOUDY or with on-the-fly non-equilibrium chemistry. We find that [CII]$_{158μ}$ predictions are robust to the model variations we explored. [OI] emission lines, that typically trace colder and denser gas relative to [CII]$_{158μ}$, are instead model-dependent, as these lines are strongly affected by the thermodynamic state of the gas and non-equilibrium photoionisation effects. For the same reasons, [OI] lines represent an excellent tool to constrain emission models, hence future observations targeting these lines will be crucial. △ Less

Submitted 25 September, 2020; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: 18 pages, 10 figures, 2 tables, matches the published version (MNRAS, 496, 5160-5175)

Star formation law in the EoR from [CII] and CIII] lines

Authors: L. Vallini, A. Ferrara, A. Pallottini, S. Carniani, S. Gallerani

Abstract: We present a novel method to simultaneously characterise the star formation law and the interstellar medium properties of galaxies in the Epoch of Reionization (EoR) through the combination of [CII] 158$μ$m (and its known relation with star formation rate) and CIII]$λ$1909Å emission line data. The method, based on a Markov Chain Monte Carlo algorithm, allows to determine the target galaxy average… ▽ More We present a novel method to simultaneously characterise the star formation law and the interstellar medium properties of galaxies in the Epoch of Reionization (EoR) through the combination of [CII] 158$μ$m (and its known relation with star formation rate) and CIII]$λ$1909Å emission line data. The method, based on a Markov Chain Monte Carlo algorithm, allows to determine the target galaxy average density, $n$, gas metallicity, $Z$, and burstiness parameter, $κ_s$, quantifying deviations from the Kennicutt-Schmidt relation. As an application, we consider COS-3018 (z=6.854), the only EoR Lyman Break Galaxy so far detected in both [CII] and CIII]. We show that COS-3018 is a moderate starburst ($κ_s \approx 3$), with $Z\approx 0.4\, Z_{\odot}$, and $n \approx 500\, {\rm cm^{-3}}$. Our method will be optimally applied to joint ALMA and JWST targets. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: 6 pages, 1 figure, Accepted for publication in MNRAS Letters

The ALPINE-ALMA [CII] Survey: Size of Individual Star-Forming Galaxies at z=4-6 and their Extended Halo Structure

Authors: S. Fujimoto, J. D. Silverman, Matthieu Bethermin, M. Ginolfi, G. C. Jones, O. Le Fèvre, M. Dessauges-Zavadsky, W. Rujopakarn, A. L. Faisst, Y. Fudamoto, P. Cassata, L. Morselli, R. Maiolino, D. Schaerer, P. Capak, L. Yan, L. Vallini, S. Toft, F. Loiacono, G. Zamorani, M. Talia, D. Narayanan, N. P. Hathi, B. C. Lemaux, Médéric Boquien , et al. (9 additional authors not shown)

Abstract: We present the physical extent of [CII] 158um line-emitting gas from 46 star-forming galaxies at z=4-6 from the ALMA Large Program to INvestigate CII at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [CII] line (r_e,[CII]) for individual galaxies and compare them with the rest-frame ultra-violet (UV) continuum (r_e,UV) from Hubble Space Telescope image… ▽ More We present the physical extent of [CII] 158um line-emitting gas from 46 star-forming galaxies at z=4-6 from the ALMA Large Program to INvestigate CII at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [CII] line (r_e,[CII]) for individual galaxies and compare them with the rest-frame ultra-violet (UV) continuum (r_e,UV) from Hubble Space Telescope images. The effective radius r_e,[CII] exceeds r_e,UV by factors of ~2-3 and the ratio of r_e,[CII]/r_e,UV increases as a function of M_star. We do not find strong evidence that [CII] line, the rest-frame UV, and FIR continuum are always displaced over ~ 1-kpc scale from each other. We identify 30% of isolated ALPINE sources as having an extended [CII] component over 10-kpc scales detected at 4.1$σ$-10.9$σ$ beyond the size of rest-frame UV and far-infrared (FIR) continuum. One object has tentative rotating features up to ~10-kpc, where the 3D model fit shows the rotating [CII]-gas disk spread over 4 times larger than the rest-frame UV-emitting region. Galaxies with the extended [CII] line structure have high star-formation rate (SFR), stellar mass (M_star), low Lya equivalent-width, and more blue-shifted (red-shifted) rest-frame UV metal absorption (Lya line), as compared to galaxies without such extended [CII] structures. Although we cannot rule out the possibility that a selection bias towards luminous objects may be responsible for such trends, the star-formation driven outflow also explains all these trends. Deeper observations are essential to test whether the extended [CII] line structures are ubiquitous to high-z star-forming galaxies. △ Less

Submitted 18 May, 2020; v1 submitted 28 February, 2020; originally announced March 2020.

Comments: 21 pages, 11 figures, ApJ in press

The ALPINE-ALMA [CII] survey. Little to no evolution in the [CII]-SFR relation over the last 13 Gyr

Authors: D. Schaerer, M. Ginolfi, M. Bethermin, Y. Fudamoto, P. A. Oesch, O. Le Fevre, A. Faisst, P. Capak, P. Cassata, J. D. Silverman, Lin Yan, G. C. Jones, R. Amorin, S. Bardelli, M. Boquien, A. Cimatti, M. Dessauges-Zavadsky, M. Giavalisco, N. P. Hathi, S. Fujimoto, E. Ibar, A. Koekemoer, G. Lagache, B. C. Lemaux, F. Loiacono , et al. (12 additional authors not shown)

Abstract: The [CII] 158 micron line is one of the strongest IR emission lines, which has been shown to trace the star-formation rate (SFR) of galaxies in the nearby Universe and up to $z \sim 2$. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at $4.4 < z< 5.9$, provides a new opportunity to examine this… ▽ More The [CII] 158 micron line is one of the strongest IR emission lines, which has been shown to trace the star-formation rate (SFR) of galaxies in the nearby Universe and up to $z \sim 2$. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at $4.4 < z< 5.9$, provides a new opportunity to examine this question with the first statistical dataset. Using the ALPINE data and earlier measurements from the literature we examine the relation between the [CII] luminosity and the SFR over the entire redshift range from $z \sim 4-8$. ALPINE galaxies, which are both detected in [CII] and dust continuum, show a good agreement with the local L([CII])-SFR relation. Galaxies undetected in the continuum with ALMA are found to be over-luminous in [CII], when the UV SFR is used. After accounting for dust-obscured star formation, by an amount SFR(IR)$\approx$SFR(UV) on average, which results from two different stacking methods and SED fitting, the ALPINE galaxies show an L([CII])-SFR relation comparable to the local one. When [CII] non-detections are taken into account, the slope may be marginally steeper at high-z, although this is still somewhat uncertain. When compared in a homogeneous manner, the $z>6 $ [CII] measurements (detections and upper limits) do not behave very differently from the $z \sim 4-6$ data. We find a weak dependence of L([CII])/SFR on the Lyman-alpha equivalent width. Finally, we find that the ratio L([CII])/LIR $\sim (1-3) \times 10^{-3}$ for the ALPINE sources, comparable to that of "normal" galaxies at lower redshift. Our analysis, which includes the largest sample ($\sim 150$ galaxies) of [CII] measurements at $z>4$ available so far, suggests no or little evolution of the L([CII])-SFR relation over the last 13 Gyr of cosmic time. △ Less

Submitted 25 May, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: 11 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics. Minor revisions, see e.g. Table A.1

Journal ref: A&A 643, A3 (2020)

The ALPINE-ALMA [CII] Survey: data processing, catalogs, and statistical source properties

Authors: M. Bethermin, Y. Fudamoto, M. Ginolfi, F. Loiacono, Y. Khusanova, P. L. Capak, P. Cassata, A. Faisst, O. Le Fevre, D. Schaerer, J. D. Silverman, L. Yan, R. Amorin, S. Bardelli, M. Boquien, A. Cimatti, I. Davidzon, M. Dessauges-Zavadsky, S. Fujimoto, C. Gruppioni, N. P. Hathi, E. Ibar, G. C. Jones, A. M. Koekemoer, G. Lagache , et al. (11 additional authors not shown)

Abstract: The ALPINE-ALMA large program targets the [CII] 158 $μ$m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z=4.4 and z=5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details of the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive ac… ▽ More The ALPINE-ALMA large program targets the [CII] 158 $μ$m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z=4.4 and z=5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details of the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates, we measured the conversion factor from the ALMA 158 $μ$m rest-frame dust continuum luminosity to the total infrared luminosity (L$_{\rm IR}$) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L$_{\rm IR}$ of 4.4$\times 10^{11}$ L$_\odot$. We also detected 57 additional continuum sources in our ALMA pointings. They are at lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5$\pm$0.2. We measured the 850 $μ$m number counts between 0.35 and 3.5 mJy, improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40 % of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8$\times$10$^8$ L$_\odot$ and their median full width at half maximum is 252 km/s. After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR-L$_{\rm [CII]}$ relations of De Looze et al. (2014) and Lagache et al. (2018). △ Less

Submitted 16 July, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: 44 pages, 29 figures, 14 tables, accepted by A&A, ALPINE products released at https://cesam.lam.fr/a2c2s/, luminosities in table B2 corrected

Journal ref: A&A 643, A2 (2020)

The ALPINE-ALMA [CII] Survey: On the nature of an extremely obscured serendipitous galaxy

Authors: Michael Romano, P. Cassata, L. Morselli, B. C. Lemaux, M. Bethermin, P. Capak, A. Faisst, O. Le Fevre, D. Schaerer, J. Silverman, L. Yan, S. Bardelli, M. Boquien, A. Cimatti, M. Dessauges-Zavadsky, A. Enia, Y. Fudamoto, S. Fujimoto, M. Ginolfi, C. Gruppioni, N. P. Hathi, E. Ibar, G. C. Jones, A. M. Koekemoer, F. Loiacono , et al. (9 additional authors not shown)

Abstract: 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.… ▽ More 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. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: 13 pages, 6 figures. Submitted to MNRAS

Outflows and extended [CII] halos in high redshift galaxies

Authors: Elia Pizzati, Andrea Ferrara, Andrea Pallottini, Simona Gallerani, Livia Vallini, Davide Decataldo, Seiji Fujimoto

Abstract: Recent stacked ALMA observations have revealed that normal, star-forming galaxies at $z\approx 6$ are surrounded by extended ($\approx 10\,\mathrm{kpc}$) [CII] emitting halos which are not predicted by the most advanced, zoom-in simulations. We present a model in which these halos are the result of supernova-driven cooling outflows. Our model contains two free parameters, the outflow mass loading… ▽ More Recent stacked ALMA observations have revealed that normal, star-forming galaxies at $z\approx 6$ are surrounded by extended ($\approx 10\,\mathrm{kpc}$) [CII] emitting halos which are not predicted by the most advanced, zoom-in simulations. We present a model in which these halos are the result of supernova-driven cooling outflows. Our model contains two free parameters, the outflow mass loading factor, $η$, and the parent galaxy dark matter halo circular velocity, $v_c$. The outflow model successfully matches the observed [CII] surface brightness profile if $η= 3.20 \pm 0.10$ and $v_c = 170 \pm 10{\,\rm km\,s^{-1}}$, corresponding to a dynamical mass of $\approx 10^{11}\, \mathrm{M}_\odot$. The predicted outflow rate and velocity range are $128 \pm 5 \,\mathrm{M}_\odot {\rm yr}^{-1}$ and $300-500 {\,\rm km\,s^{-1}}$, respectively. We conclude that: (a) extended halos can be produced by cooling outflows; (b) the large $η$ value is marginally consistent with starburst-driven outflows, but it might indicate additional energy input from AGN; (c) the presence of [CII] halos requires an ionizing photon escape fraction from galaxies $f_{\rm esc} \ll 1$. The model can be readily applied also to individual high-$z$ galaxies, as those observed, e.g., by the ALMA ALPINE survey now becoming available. △ Less

Submitted 24 April, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

Comments: 13 pages, 7 figures; accepted for publication in MNRAS

The ALPINE-ALMA [CII] Survey: Multi-Wavelength Ancillary Data and Basic Physical Measurements

Authors: A. L. Faisst, D. Schaerer, B. C. Lemaux, P. A. Oesch, Y. Fudamoto, P. Cassata, M. Bethermin, P. L. Capak, O. Le Fevre, J. D. Silverman, L. Yan, M. Ginolfi, A. M. Koekemoer, L. Morselli, R. Amorin, S. Bardelli, M. Boquien, G. Brammer, A. Cimatti, M. Dessauges-Zavadsky, S. Fujimoto, C. Gruppioni, N. P. Hathi, S. Hemmati, E. Ibar , et al. (15 additional authors not shown)

Abstract: 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$σ$) and the surround… ▽ More 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$σ$) 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 ($β$), 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 ($Δβ$ ~ 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]$λ$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. △ Less

Submitted 12 March, 2020; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: 48 pages, 33 figures, 7 tables (including Appendix), accepted for publication in ApJS

Large Population of ALMA Galaxies at z>6 with Very High [OIII]88um to [CII]158um Flux Ratios: Evidence of Extremely High Ionization Parameter or PDR Deficit?

Authors: Yuichi Harikane, Masami Ouchi, Akio K. Inoue, Yoshiki Matsuoka, Yoichi Tamura, Tom Bakx, Seiji Fujimoto, Kana Moriwaki, Yoshiaki Ono, Tohru Nagao, Ken-ichi Tadaki, Takashi Kojima, Takatoshi Shibuya, Eiichi Egami, Andrea Ferrara, Simona Gallerani, Takuya Hashimoto, Kotaro Kohno, Yuichi Matsuda, Hiroshi Matsuo, Andrea Pallottini, Yuma Sugahara, Livia Vallini

Abstract: We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to esti… ▽ More We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find. △ Less

Submitted 24 May, 2020; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: 20 pages, 18 figures, Accepted for publication in ApJ

Impact of X-rays on CO emission from high-z galaxies

Authors: L. Vallini, A. G. G. M. Tielens, A. Pallottini, S. Gallerani, C. Gruppioni, S. Carniani, F. Pozzi, M. Talia

Abstract: We study the impact of active galactic nuclei (AGN) on the CO Spectral Line Energy Distribution (SLED) of high-$z$ galaxies. In particular, we want to assess if the CO SLED can be used as a probe of AGN activity. To this purpose, we develop a semi-analytical model that takes into account the radiative transfer and the clumpy structure of giant molecular clouds where the CO lines are excited, their… ▽ More We study the impact of active galactic nuclei (AGN) on the CO Spectral Line Energy Distribution (SLED) of high-$z$ galaxies. In particular, we want to assess if the CO SLED can be used as a probe of AGN activity. To this purpose, we develop a semi-analytical model that takes into account the radiative transfer and the clumpy structure of giant molecular clouds where the CO lines are excited, their distribution in the galaxy disk, and the torus obscuration of the AGN radiation. We study the joint effect on the CO SLED excitation of (i) the X-ray luminosity from the AGN ($L_{X}$), (ii) the size of the molecular disk, (iii) the inclination angle ($Ω$) of the torus with respect to the molecular disk, and (iv) the GMC mean density. We also discuss the possibility of an enhanced Cosmic Ray Ionization Rate (CRIR). We find that the X-ray Dominated Region (XDR) generated by the AGN in every case enhances the CO SLED for $J>5$, with increasing excitation of high-$J$ CO lines for increasing X-ray luminosity. Because high-$z$ galaxies are compact, the XDR region typically encloses the whole disk, thus its effect can be more important with respect to lower redshift objects. The impact of the XDR can be disentangled from an enhanced CRIR either if $L_X>10^{44} \rm \,erg\, s^{-1}$, or if $Ω\geq 60^{\circ}$. We finally provide predictions on the CO(7-6)/[CII] and CO(17-16)/[CII] ratios as a function of $L_X$, which can be relevant for ALMA follow up of galaxies and quasars previously detected in [CII]. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: accepted for publication in MNRAS