Showing 1–50 of 120 results for author: Ciesla, L

NOEMA formIng Cluster survEy (NICE): Characterizing eight massive galaxy groups at $1.5 < z < 4$ in the COSMOS field

Authors: Nikolaj B. Sillassen, Shuowen Jin, Georgios E. Magdis, Emanuele Daddi, Tao Wang, Shiying Lu, Hanwen Sun, Vinod Arumugam, Daizhong Liu, Malte Brinch, Chiara D'Eugenio, Raphael Gobat, Carlos Gómez-Guijarro, Michael Rich, Eva Schinnerer, Veronica Strazzullo, Qinghua Tan, Francesco Valentino, Yijun Wang, Mengyuan Xiao, Luwenjia Zhou, David Blánquez-Sesé, Zheng Cai, Yanmei Chen, Laure Ciesla , et al. (19 additional authors not shown)

Abstract: The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are c… ▽ More The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are confirmed with ALMA, and one is confirmed by H$α$ from Subaru/FMOS. We constructed the integrated FIR SEDs for the eight groups, obtaining total IR SFR $=260-1300~{\rm M_\odot}$~yr$^{-1}$. We adopted six methods to estimate the dark matter masses, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass density. We found the radial stellar mass density are consistent with a NFW profile, supporting that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are $\log(M_{\rm h}/{\rm M_\odot})=12.8-13.7$ with uncertainty of 0.3 dex. From halo mass estimates, we derive baryonic accretion rate ${\rm BAR}=(1-8)\times10^{3}\,{\rm M_{\odot}/yr}$ for this sample. We find a quasi-linear correlation between the integrated SFR/BAR and the theoretical halo mass limit for cold streams, $M_{\rm stream}/M_{\rm h}$, with ${\rm SFR/BAR}=10^{-0.46\pm0.22}\left({M_{\rm stream}/M_{\rm h}}\right)^{0.71\pm0.16}$ with a scatter of $0.40\,{\rm dex}$. Further, we compare halo masses and stellar masses with simulations, and find all structures are consistent with being progenitors of $M_{\rm h}(z=0)>10^{14}\,{\rm M_{\odot}}$ galaxy clusters, and the most massive central galaxies have stellar masses consistent with brightest cluster galaxies (BCGs) progenitors in the TNG300 simulation. The results strongly suggest these structures are forming massive galaxy clusters via baryonic and dark matter accretion. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 44 pages (27pp appendix), 32 figures, 18 tables, accepted for publication in A&A

A$^3$COSMOS: Measuring the cosmic dust-attenuated star formation rate density at $4 < z < 5$

Authors: Benjamin Magnelli, Sylvia Adscheid, Tsan-Ming Wang, Laure Ciesla, Emanuele Daddi, Ivan Delvecchio, David Elbaz, Yoshinobu Fudamoto, Shuma Fukushima, Maximilien Franco, Carlos Gómez-Guijarro, Carlotta Gruppioni, Eric F. Jiménez-Andrade, Daizhong Liu, Pascal Oesch, Eva Schinnerer, Alberto Traina

Abstract: [Abridged] In recent years, conflicting results have provided an uncertain view of the dust-attenuated properties of $z>4$ star-forming galaxies (SFGs). To solve this, we used the deepest data publicly available in COSMOS to build a mass-complete ($>10^{9.5}\,M_{\odot}$) sample of SFGs at $4<z<5$ and measured their dust-attenuated properties by stacking all archival ALMA band 6 and 7 observations… ▽ More [Abridged] In recent years, conflicting results have provided an uncertain view of the dust-attenuated properties of $z>4$ star-forming galaxies (SFGs). To solve this, we used the deepest data publicly available in COSMOS to build a mass-complete ($>10^{9.5}\,M_{\odot}$) sample of SFGs at $4<z<5$ and measured their dust-attenuated properties by stacking all archival ALMA band 6 and 7 observations available. Combining this information with their rest-frame ultraviolet emission from the COSMOS2020 catalog, we constrained the IRX ($\equiv L_{\rm IR}/L_{\rm UV}$)--$β_{\rm UV}$, IRX--$M_\ast$, and SFR--$M_\ast$ relations at $z\sim4.5$. Finally, using these relations and the stellar mass function of SFGs at $z\sim4.5$, we inferred the unattenuated and dust-attenuated SFRD at this epoch. SFGs at $z\sim4.5$ follow an IRX--$β_{\rm UV}$ relation that is consistent with that of local starbursts, while they follow a steeper IRX--$M_\ast$ relation than observed locally. The grain properties of dust in these SFGs seems thus similar to those in local starbursts but its mass and geometry result in lower attenuation in low-mass SFGs. SFGs at $z\sim4.5$ lie on a linear SFR--$M_\ast$ relation, whose normalization varies by 0.3 dex, when we exclude or include from our stacks the ALMA primary targets. The cosmic SFRD$(>M_\ast)$ converges at $M_\ast<10^{9}\,M_\odot$ and is dominated by SFGs with $M_\ast\sim10^{9.5-10.5}\,M_\odot$. The fraction of the cosmic SFRD that is attenuated by dust, ${\rm SFRD}_{\rm IR}(>M_\ast)/ {\rm SFRD}(>M_\ast)$, is $90\pm4\%$ for $M_\ast\,=\,10^{10}\,M_\odot$, $68\pm10\%$ for $M_\ast=10^{8.9}\,M_\odot$ (i.e., $0.03\times M^\star$; $M^\star$ being the characteristic stellar mass of SFGs) and this value converges to $60\pm10\%$ for $M_\ast=10^{8}\,M_\odot$. Even at this early epoch, the fraction of the cosmic SFRD that is attenuated by dust remains thus significant. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: Accepted for publication in A&A; 14 pages and 8 figures

Disentangling the co-evolution of galaxies and supermassive black holes with PRIMA

Authors: L. Bisigello, C. Gruppioni, A. Bolatto, L. Ciesla, A. Pope, L. Armus, L., J. D. Smith, R. Somerville, L. Y. A. Yung, R. J. Wright, C. M. Bradford, J. Glenn, A. Feltre

Abstract: The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-I… ▽ More The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-IR provide unmatched diagnostic power to determine the properties of gas and dust, measure gas-phase metallicities and map cold galactic outflows in even the most obscured galaxies. By combining multi-band photometric surveys with low and high-resolution far-IR spectroscopy, the PRobe far-Infrared Mission for Astrophysics (PRIMA), a concept for a far-IR, 1.8m-diameter, cryogenically cooled observatory, will revolutionise the field of galaxy evolution by taking advantage of this IR toolkit to find and study dusty galaxies across galactic time. In this work, we make use of the phenomenological simulation SPRITZ and the Santa Cruz semi-analytical model to describe how a moderately deep multi-band PRIMA photometric survey can easily reach beyond previous IR missions to detect and study galaxies down to $10^{11}\,L_{\odot}$ beyond cosmic noon and at least up to z=4, even in the absence of gravitational lensing. By decomposing the spectral energy distribution (SED) of these photometrically selected galaxies, we show that PRIMA can be used to accurately measure the relative AGN power, the mass fraction contributed by polycyclic aromatic hydrocarbon (PAH) and the total IR luminosity. At the same time, spectroscopic follow up with PRIMA will allow to trace both the star formation and black hole accretion rates (SFR, BHAR), the gas phase metallicities and the mass outflow rates of cold gas in hundreds to thousands of individual galaxies to z=2. △ Less

Submitted 3 June, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: 16 pages, 17 figures, accepted for publication in A&A

Overcoming Confusion Noise with Hyperspectral Imaging from PRIMAger

Authors: James M. S. Donnellan, Seb J. Oliver, Matthieu Bethermin, Longji Bing, Alberto Bolatto, Charles M. Bradford, Denis Burgarella, Laure Ciesla, Jason Glenn, Alexandra Pope, Stephen Serjeant, Raphael Shirley, JD T. Smith, Chris Sorrell

Abstract: The PRobe far-Infrared Mission for Astrophysics (PRIMA) concept aims to perform mapping with spectral coverage and sensitivities inaccessible to previous FIR space telescopes. PRIMA's imaging instrument, PRIMAger, provides unique hyperspectral imaging simultaneously covering 25-235 $μ$m. We synthesise images representing a deep, 1500 hr deg$^{-2}$ PRIMAger survey, with realistic instrumental and c… ▽ More The PRobe far-Infrared Mission for Astrophysics (PRIMA) concept aims to perform mapping with spectral coverage and sensitivities inaccessible to previous FIR space telescopes. PRIMA's imaging instrument, PRIMAger, provides unique hyperspectral imaging simultaneously covering 25-235 $μ$m. We synthesise images representing a deep, 1500 hr deg$^{-2}$ PRIMAger survey, with realistic instrumental and confusion noise. We demonstrate that we can construct catalogues of galaxies with a high purity ($>95$ per cent) at a source density of 42k deg$^{-2}$ using PRIMAger data alone. Using the XID+ deblending tool we show that we measure fluxes with an accuracy better than 20 per cent to flux levels of 0.16, 0.80, 9.7 and 15 mJy at 47.4, 79.7, 172, 235 $μ$m respectively. These are a factor of $\sim$2 and $\sim$3 fainter than the classical confusion limits for 72-96 $μ$m and 126-235 $μ$m, respectively. At $1.5 \leq z \leq 2$, we detect and accurately measure fluxes in 8-10 of the 10 channels covering 47-235 $μ$m for sources with $2 \leq$ log(SFR) $\leq 2.5$, a 0.5 dex improvement on what might be expected from the classical confusion limit. Recognising that PRIMager will operate in a context where high quality data will be available at other wavelengths, we investigate the benefits of introducing additional prior information. We show that by introducing even weak prior flux information when employing a higher source density catalogue (more than one source per beam) we can obtain accurate fluxes an order of magnitude below the classical confusion limit for 96-235 $μ$m. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 14 pages, 11 figures

Confusion of extragalactic sources in the far infrared: a baseline assessment of the performance of PRIMAger in intensity and polarization

Authors: Matthieu Béthermin, Alberto D. Bolatto, François Boulanger, Charles M. Bradford, Denis Burgarella, Laure Ciesla, James Donnellan, Brandon S. Hensley, Jason Glenn, Guilaine Lagache, Enrique Lopez-Rodriguez, Seb Oliver, Alexandra Pope, Marc Sauvage

Abstract: Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in in… ▽ More Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in intensity and polarization. We aim to provide predictions of the confusion level and its consequences for future surveys. We produced simulated PRIMAger maps affected only by the confusion noise using the simulated infrared extragalactic sky (SIDES) semi-empirical simulation. We then estimated the confusion limit in these maps and extracted the sources using a basic blind extractor. By comparing the input galaxy catalog and the extracted source catalog, we derived various performance metrics as completeness, purity, and the accuracy of various measurements. In intensity, we predict that the confusion limit increases rapidly with increasing wavelength. The confusion limit in polarization is more than 100x lower. The measured flux density is dominated by the brightest galaxy in the beam, but other objects also contribute at longer wavelength (~30% at 235 $μ$m). We also show that galaxy clustering has a mild impact on confusion in intensity (up to 25%), while it is negligible in polarization. In intensity, a basic blind extraction will be sufficient to detect galaxies at the knee of the luminosity function up to z~3 and 10$^{11}$ M$_\odot$ main-sequence galaxies up to z~5. In polarization for a conservative sensitivity, we expect ~8 000 detections up to z=2.5 opening a totally new window on the high-z dust polarization. Finally, we show that intensity surveys at short wavelength and polarization surveys at long wavelength tend to reach confusion at similar depth. There is thus a strong synergy. △ Less

Submitted 28 April, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

Comments: 21 pages, 15 figures, 4 tables, submitted to A&A, minor update on the names of the sensitivity scenarii (for consistency with other papers)

CEERS: Increasing Scatter along the Star-Forming Main Sequence Indicates Early Galaxies Form in Bursts

Authors: Justin W. Cole, Casey Papovich, Steven L. Finkelstein, Micaela B. Bagley, Mark Dickinson, Kartheik G. Iyer, L. Y. Aaron Yung, Laure Ciesla, Ricardo O. Amorin, Pablo Arrabal Haro, Rachana Bhatawdekar, Antonello Calabro, Nikko J. Cleri, Alexander de la Vega, Avishai Dekel, Ryan Endsley, Eric Gawiser, Mauro Giavalisco, Nimish P. Hathi, Michaela Hirschmann, Benne W. Holwerda, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Ray A. Lucas, Sara Mascia , et al. (7 additional authors not shown)

Abstract: We present the star-formation-rate -- stellar-mass (SFR-M$_\ast$) relation for galaxies in the CEERS survey at $4.5\leq z\leq 12$. We model the \jwst\ and \hst\ rest-UV and rest-optical photometry of galaxies with flexible star-formation histories (SFHs) using \bagpipes. We consider SFRs averaged from the SFHs over 10~Myr (\sfrten) and 100~Myr (\sfrcen), where the photometry probes SFRs on these t… ▽ More We present the star-formation-rate -- stellar-mass (SFR-M$_\ast$) relation for galaxies in the CEERS survey at $4.5\leq z\leq 12$. We model the \jwst\ and \hst\ rest-UV and rest-optical photometry of galaxies with flexible star-formation histories (SFHs) using \bagpipes. We consider SFRs averaged from the SFHs over 10~Myr (\sfrten) and 100~Myr (\sfrcen), where the photometry probes SFRs on these timescales, effectively tracing nebular emission lines in the rest-optical (on $\sim10$~Myr timescales) and the UV/optical continuum (on $\sim100$ Myr timescales). We measure the slope, normalization and intrinsic scatter of the SFR-M$_\ast$ relation, taking into account the uncertainty and the covariance of galaxy SFRs and $M_\ast$. From $z\sim 5-9$ there is larger scatter in the $\sfrten-M_\ast$ relation, with $σ(\log \sfrcen)=0.4$~dex, compared to the $\sfrcen-M_\ast$ relation, with $σ(\log \sfrten)=0.1$~dex. This scatter increases with redshift and increasing stellar mass, at least out to $z\sim 7$. These results can be explained if galaxies at higher redshift experience an increase in star-formation variability and form primarily in short, active periods, followed by a lull in star formation (i.e. ``napping'' phases). We see a significant trend in the ratio $R_\mathrm{SFR}=\log(\sfrten/\sfrcen)$ in which, on average, $R_\mathrm{SFR}$ decreases with increasing stellar mass and increasing redshift. This yields a star-formation ``duty cycle'' of $\sim40\%$ for galaxies with $\log M_\ast/M_\odot\geq 9.3$, at $z\sim5$, declining to $\sim20\%$ at $z\sim9$. Galaxies also experience longer lulls in star formation at higher redshift and at higher stellar mass, such that galaxies transition from periods of higher SFR variability at $z\gtrsim~6$ to smoother SFR evolution at $z\lesssim~4.5$. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 28 pages, 12 figures, 2 Appendix figures

MAGIC: Muse gAlaxy Groups In Cosmos -- A survey to probe the impact of environment on galaxy evolution over the last 8 Gyr

Authors: B. Epinat, T. Contini, W. Mercier, L. Ciesla, B. C. Lemaux, S. D. Johnson, J. Richard, J. Brinchmann, L. A. Boogaard, D. Carton, L. Michel-Dansac, R. Bacon, D. Krajnovic, H. Finley, I. Schroetter, E. Ventou, V. Abril-Melgarejo, A. Boselli, N. F. Bouché, W. Kollatschny, K. Kovac, M. Paalvast, G. Soucail, T. Urrutia, P. M. Weilbacher

Abstract: We introduce the MUSE gAlaxy Groups in COSMOS (MAGIC) survey, which was built to study the impact of environment on galaxy evolution over the last 8 Gyr. It consists of 17 MUSE fields targeting 14 massive structures at intermediate redshift ($0.3<z<0.8$) in the COSMOS area. We securely measured the redshifts for 1419 sources and identified 76 galaxy pairs and 67 groups of at least 3 members using… ▽ More We introduce the MUSE gAlaxy Groups in COSMOS (MAGIC) survey, which was built to study the impact of environment on galaxy evolution over the last 8 Gyr. It consists of 17 MUSE fields targeting 14 massive structures at intermediate redshift ($0.3<z<0.8$) in the COSMOS area. We securely measured the redshifts for 1419 sources and identified 76 galaxy pairs and 67 groups of at least 3 members using a friends-of-friends algorithm. The environment of galaxies is quantified from group properties, as well as from global and local density estimators. The MAGIC survey has increased the number of objects with a secure spectroscopic redshift over its footprint by a factor of about 5. Most of the new redshifts have apparent magnitudes in the $z^{++}$ band $z_{app}^{++}>21.5$. The spectroscopic redshift completeness is high: in the redshift range of [OII] emitters ($0.25 \le z < 1.5$), where most of the groups are found, it globally reaches a maximum of 80% down to $z_{app}^{++}=25.9$, and locally decreases from $\sim 100$% to $\sim50$% in magnitude bins from $z_{app}^{++}=23-24$ to $z_{app}^{++}=25.5$. We find that the fraction of quiescent galaxies increases with local density and with the time spent in groups. A morphological dichotomy is also found between bulge-dominated quiescent and disk-dominated star-forming galaxies. As environment gets denser, the peak of the stellar mass distribution shifts towards $M_*>10^{10}~M_\odot$, and the fraction of galaxies with $M_*<10^9~M_\odot$ decreases significantly, even for star-forming galaxies. We also highlight peculiar features such as close groups, extended nebulae, and a gravitational arc. Our results suggest that galaxies are preprocessed in groups of increasing mass before entering rich groups and clusters. We publicly release two catalogs containing the properties of galaxies and groups, respectively. △ Less

Submitted 28 March, 2024; v1 submitted 1 December, 2023; originally announced December 2023.

Comments: 27 pages, 22 figures, accepted for publication in A&A

Journal ref: A&A, 683, A205 (2024)

The COSMOS-Web ring: in-depth characterization of an Einstein ring lensing system at z~2

Authors: W. Mercier, M. Shuntov, R. Gavazzi, J. W. Nightingale, R. Arango, O. Ilbert, A. Amvrosiadis, L. Ciesla, C. Casey, S. Jin, A. L. Faisst, I. T. Andika, N. E. Drakos, A. Enia, M. Franco, S. Gillman, G. Gozaliasl, C. C. Hayward, M. Huertas-Company, J. S. Kartaltepe, A. M. Koekemoer, C. Laigle, D. Le Borgne, G. Magdis, G. Mahler , et al. (12 additional authors not shown)

Abstract: Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED f… ▽ More Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED fitting codes. Using JWST/NIRCam images, we also produce two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. The lens is a very massive and quiescent (sSFR < 10^(-13) yr-1) elliptical galaxy at z = 2.02 \pm 0.02 with a total mass Mtot(<thetaE) = (3.66 \pm 0.36) x 10^11 Msun and a stellar mass M* = (1.37 \pm 0.14) x 10^11 Msun. Compared to SHMRs from the literature, we find that the total mass is consistent with the presence of a DM halo of mass Mh = 1.09^(+1.46)_(-0.57) x 10^13 Msun. In addition, the background source is a M* = (1.26 \pm 0.17) x 10^10 Msun star-forming galaxy (SFR=(78 \pm 15) Msun/yr) at z = 5.48 \pm 0.06. Its reconstructed morphology shows two components with different colors. Dust attenuation values from SED fitting and nearby detections in the FIR also suggest it could be partially dust-obscured. Conclusions. We find the lens at z=2. Its total, stellar, and DM halo masses are consistent within the Einstein ring, so we do not need any unexpected changes in our description of the lens (e.g. change its IMF or include a non-negligible gas contribution). The most likely solution for the lensed source is at z = 5.5. Its reconstructed morphology is complex and highly wavelength dependent, possibly because it is a merger or a main sequence galaxy with a heterogeneous dust distribution. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 16 pages, submitted to A&A

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

Identification of a transition from stochastic to secular star formation around $z=9$ with JWST

Authors: L. Ciesla, D. Elbaz, O. Ilbert, V. Buat, B. Magnelli, D. Narayanan, E. Daddi, C. Gómez-Guijarro, R. Arango-Toro

Abstract: Star formation histories (SFH) of early (6$<z<$12) galaxies have been found to be highly stochastic in both simulations and observations, while at $z\lesssim$6 the presence of a main sequence (MS) of star-forming galaxies imply secular processes at play. In this work, we aim at characterising the SFH variability of early galaxies as a function of their stellar mass and redshift. We use the JADES p… ▽ More Star formation histories (SFH) of early (6$<z<$12) galaxies have been found to be highly stochastic in both simulations and observations, while at $z\lesssim$6 the presence of a main sequence (MS) of star-forming galaxies imply secular processes at play. In this work, we aim at characterising the SFH variability of early galaxies as a function of their stellar mass and redshift. We use the JADES public catalogue and derive the physical properties of the galaxies as well as their SFH using the spectral energy distribution modelling code CIGALE. To this aim, we implement a non-parametric SFH with a flat prior allowing for as much stochasticity as possible. We use the SFR gradient, an indicator of the movement of galaxies on the SFR-$M_\ast$ plane, linked to the recent SFH of galaxies. This dynamical approach of the relation between the SFR and stellar mass allows us to show that, at $z>9$, 87% of massive galaxies, ($\log(M_\ast/M_\odot)\gtrsim$9), have SFR gradients consistent with a stochastic star-formation activity during the last 100 Myr, while this fraction drops to 15% at $z<7$. On the other hand, we see an increasing fraction of galaxies with a star-formation activity following a common stream on the SFR-$M_\ast$ plane with cosmic time, indicating that a secular mode of star-formation is emerging. We place our results in the context of the observed excess of UV emission as probed by the UV luminosity function at $z\gtrsim10$, by estimating $σ_{UV}$, the dispersion of the UV absolute magnitude distribution, to be of the order of 1.2mag and compare it with predictions from the literature. In conclusion, we find a transition of star-formation mode happening around $z\sim9$: Galaxies with stochastic SFHs dominates at $z\gtrsim9$, although this level of stochasticity is too low to reach those invoked by recent models to reproduce the observed UV luminosity function. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: Submitted to A&A

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

Stellar angular momentum of disk galaxies at z = 0.7 in the MAGIC survey I. Impact of the environment

Authors: W. Mercier, B. Epinat, T. Contini, D. Krajnović, L. Ciesla, B. C. Lemaux, V. Abril-Melgarejo, L. Boogaard, D. Pelliccia

Abstract: Aims: At intermediate redshift, galaxy groups/clusters are thought to impact galaxies (e.g. their angular momentum). We investigate whether the environment has an impact on the galaxies' angular momentum and identify underlying driving physical mechanisms. Methods: We derive robust estimates of the stellar angular momentum using Hubble Space Telescope (HST) images combined with spatially resolve… ▽ More Aims: At intermediate redshift, galaxy groups/clusters are thought to impact galaxies (e.g. their angular momentum). We investigate whether the environment has an impact on the galaxies' angular momentum and identify underlying driving physical mechanisms. Methods: We derive robust estimates of the stellar angular momentum using Hubble Space Telescope (HST) images combined with spatially resolved ionised gas kinematics from the Multi-Unit Spectroscopic Explorer (MUSE) for a sample of ~200 galaxies in groups and in the field at z~0.7 drawn from the MAGIC survey. Using various environmental tracers, we study the position of the galaxies in the the angular momentum-stellar mass (Fall) relation as a function of environment. Results: We measure a 0.12 dex (2sigma significant) depletion of angular momentum for low-mass galaxies (M* < 10^10 Msun) in groups with respect to the field. Massive galaxies located in dense environments have less angular momentum than expected from the low-mass Fall relation but, without a comparable field sample, we cannot infer whether this effect is mass- or environmentally-driven. Furthermore, massive galaxies are found in the centre of the structures and have low systemic velocities. The observed depletion of angular momentum at low mass does not appear linked with the strength of the over-density around the galaxies but it is strongly correlated with the galaxies' systemic velocity normalised by the dispersion of their host group and with their ionised gas velocity dispersion. Conclusions: Group galaxies seem depleted in angular momentum, especially at low mass. Our results suggest that this depletion might be induced by physical mechanisms that scale with the systemic velocity of the galaxies (e.g. stripping or merging) and that such mechanism might be responsible for enhancing the velocity dispersion of the gas as galaxies lose angular momentum. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 23 pages, including 5 appendices, 14 figures, accepted version before language correction

Journal ref: A&A 677, A143 (2023)

The relation of cosmic environment and morphology with the star formation and stellar populations of AGN and non-AGN galaxies

Authors: G. Mountrichas, G. Yang, V. Buat, B. Darvish, M. Boquien, Q. Ni, D. Burgarella, L. Ciesla

Abstract: In this work, we study the relation of cosmic environment and morphology with the star-formation (SF) and the stellar population of galaxies. Most importantly, we examine if this relation differs for systems with active and non-active supermassive black holes. For that purpose, we use 551 X-ray detected active galactic nuclei (AGN) and 16,917 non-AGN galaxies in the COSMOS-Legacy survey, for which… ▽ More In this work, we study the relation of cosmic environment and morphology with the star-formation (SF) and the stellar population of galaxies. Most importantly, we examine if this relation differs for systems with active and non-active supermassive black holes. For that purpose, we use 551 X-ray detected active galactic nuclei (AGN) and 16,917 non-AGN galaxies in the COSMOS-Legacy survey, for which the surface-density field measurements are available. The sources lie at redshift of $\rm 0.3<z<1.2$, probe X-ray luminosities of $\rm 42<log\,[L_{X,2-10keV}(erg\,s^{-1})]<44$ and have stellar masses, $\rm 10.5<log\,[M_*(M_\odot)]<11.5$. Our results show that isolated AGN (field) have lower SFR compared to non AGN, at all L$_X$ spanned by our sample. However, in denser environments (filaments, clusters), moderate L$_X$ AGN ($\rm log\,[L_{X,2-10keV}(erg\,s^{-1})]>43$) and non-AGN galaxies have similar SFR. We, also, examine the stellar populations and the morphology of the sources in different cosmic fields. For the same morphological type, non-AGN galaxies tend to have older stellar populations and are less likely to have undergone a recent burst in denser environments compared to their field counterparts. The differences in the stellar populations with the density field are, mainly, driven by quiescent systems. Moreover, low L$_X$ AGN present negligible variations of their stellar populations, in all cosmic environments, whereas moderate L$_X$ AGN have, on average, younger stellar populations and are more likely to have undergone a recent burst, in high density fields. Finally, in the case of non-AGN galaxies, the fraction of bulge-dominated (BD) systems increases with the density field, while BD AGN are scarce in denser environments. Our results are consistent with a scenario in which a common mechanism, such as mergers, triggers both the SF and the AGN activity. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A. 10 pages, 5 figures, 3 tables

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

Probing the Timescale of the 1.4 GHz Radio emissions as a Star formation tracer

Authors: R. C. Arango-Toro, L. Ciesla, O. Ilbert, B. Magnelli, E. F. Jiménez-Andrade, V. Buat

Abstract: Radio used as a star formation rate (SFR) tracer presents enormous advantages by being unaffected by dust and radio sources being pinpointed at the sub-arc-second level. The interpretation of the low frequency 1.4 GHz luminosity is hampered by the difficulty in modeling the cosmic ray paths in the interstellar medium, and their interactions with the magnetic field. In this work, we compare the SFR… ▽ More Radio used as a star formation rate (SFR) tracer presents enormous advantages by being unaffected by dust and radio sources being pinpointed at the sub-arc-second level. The interpretation of the low frequency 1.4 GHz luminosity is hampered by the difficulty in modeling the cosmic ray paths in the interstellar medium, and their interactions with the magnetic field. In this work, we compare the SFR derived from radio observations, and the ones derived from spectral energy distribution (SED) modeling. We aim at better understand the behavior of the SFR radio tracer, with a specific emphasis on the link with star-formation histories. We used the SED modeling code Code Investigating GALaxy Emission, CIGALE, with a non-parametric star formation history model (SFH) and fit the data over the wavelength range from the ultraviolet (UV) up to the mid-infrared (mid-IR). We interpret the difference between radio and SED-based SFR tracers in the light of recent gradients in the derived SFH. To validate the robustness of the results, we checked for any remaining active galaxy nuclei (AGN) contribution and tested the impact of our SFH modeling approach. Approximately 27% our galaxies present a radio SFR (SFR$_{\rm radio}$) at least ten times larger than the instantaneous SFR from SED-fitting (SFR$_{\rm SED}$). This trend affects primarily the galaxies that show a declining SFH activity over the last 300 Myr. Both SFR indicators converge toward a consistent value, when the SFHs are averaged over a period larger than 150 Myr to derive SFR$_{\rm SED}$. Although the radio at low frequency 1.4 GHz is a good tracer of the star formation activity of galaxies with constant or increasing SFH, our results indicate that this is not the case for galaxies that are quenching. Our analysis suggests that the star formation time sensitivity of the radio low frequency could be longer than 150 Myr. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 10 pages, 10 figures

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

The slippery slope of dust attenuation curves: Correlation of dust attenuation laws with star-to-dust compactness up to z = 4

Authors: M. Hamed, K. Małek, V. Buat, Junais, L. Ciesla, D. Donevski, G. Riccio, M. Figueira

Abstract: Aims. We investigate dust attenuation of 122 heavily dust-obscured galaxies detected with the Atacama Large Millimeter Array (ALMA) and Herschel in the COSMOS field. We search for correlations between dust attenuation recipes and the variation of physical parameters, mainly the effective radii of galaxies, their star formation rates (SFR), and stellar masses, and aim to understand which of the com… ▽ More Aims. We investigate dust attenuation of 122 heavily dust-obscured galaxies detected with the Atacama Large Millimeter Array (ALMA) and Herschel in the COSMOS field. We search for correlations between dust attenuation recipes and the variation of physical parameters, mainly the effective radii of galaxies, their star formation rates (SFR), and stellar masses, and aim to understand which of the commonly used laws best describes dust attenuation in dusty star-forming galaxies at high redshift. Methods. We make use of the extensive photometric coverage of the COSMOS data combined with highly-resolved dust continuum maps from ALMA. We use CIGALE to estimate various physical properties of these dusty objects, mainly their SFR, their stellar masses and their attenuation. We infer galaxy effective radii (Re) using GALFIT in the Y band of HSC and ALMA continuum maps. We use these radii to investigate the relative compactness of the dust continuum and the extension of the rest-frame UV/optical Re(y)/Re(ALMA). Results. We find that the physical parameters calculated from our models strongly depend on the assumption of dust attenuation curve. As expected, the most impacted parameter is the stellar mass, which leads to a change in the "starburstiness" of the objects. We find that taking into account the relative compactness of star-to-dust emission prior to SED fitting is crucial, especially when studying dust attenuation of dusty star-forming galaxies. Shallower attenuation curves did not show a clear preference of compactness with attenuation, while the Calzetti attenuation curve preferred comparable spatial extent of unattenuated stellar light and dust emission. The evolution of the Re(UV)/Re(ALMA) ratio with redshift peeks around the cosmic noon in our sample of DSFGs, showing that this compactness is correlated with the cosmic SFR density of these dusty sources. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 18 pages, 13 figures, accepted for publication in A&A. Abstract abridged for arXiv submission

Journal ref: A&A 674, A99 (2023)

JWST CEERS probes the role of stellar mass and morphology in obscuring galaxies

Authors: Carlos Gómez-Guijarro, Benjamin Magnelli, David Elbaz, Stijn Wuyts, Emanuele Daddi, Aurélien Le Bail, Mauro Giavalisco, Mark Dickinson, Pablo G. Pérez-González, Pablo Arrabal Haro, Micaela B. Bagley, Laura Bisigello, Véronique Buat, Denis Burgarella, Antonello Calabrò, Caitlin M. Casey, Yingjie Cheng, Laure Ciesla, Avishai Dekel, Henry C. Ferguson, Steven L. Finkelstein, Maximilien Franco, Norman A. Grogin, Benne W. Holwerda, Shuowen Jin , et al. (16 additional authors not shown)

Abstract: In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JW… ▽ More In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs. △ Less

Submitted 4 September, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: Accepted by A&A. 13 pages, 9 figures

Journal ref: A&A 677, A34 (2023)

A z=1.85 galaxy group in CEERS: evolved, dustless, massive intra-halo light and a brightest group galaxy in the making

Authors: Rosemary T. Coogan, Emanuele Daddi, Aurélien Le Bail, David Elbaz, Mark Dickinson, Mauro Giavalisco, Carlos Gómez-Guijarro, Alexander de la Vega, Micaela Bagley, Steven L. Finkelstein, Maximilien Franco, Asantha R. Cooray, Peter Behroozi, Laura Bisigello, Caitlin M. Casey, Laure Ciesla, Paola Dimauro, Alexis Finoguenov, Anton M. Koekemoer, Ray A. Lucas, Pablo G. Pérez-González, L. Y. Aaron Yung, Pablo Arrabal Haro, Jeyhan S. Kartaltepe, Shardha Jogee , et al. (3 additional authors not shown)

Abstract: We present CEERS JWST/NIRCam imaging of a massive galaxy group at z=1.85, to explore the early JWST view on massive group formation in the distant Universe. The group contains >16 members (including 6 spectros. confirmations) down to log10(Mstar/Msun)=8.5, including the brightest group galaxy (BGG) in the process of actively assembling at this redshift. The BGG is comprised of multiple merging com… ▽ More We present CEERS JWST/NIRCam imaging of a massive galaxy group at z=1.85, to explore the early JWST view on massive group formation in the distant Universe. The group contains >16 members (including 6 spectros. confirmations) down to log10(Mstar/Msun)=8.5, including the brightest group galaxy (BGG) in the process of actively assembling at this redshift. The BGG is comprised of multiple merging components extending ~3.6" (30kpc) across the sky. The BGG contributes 69% of the group's total galactic stellar mass, with one of the merging components containing 76% of the total mass of the BGG and a SFR>1810Msun/yr. Most importantly, we detect intra-halo light (IHL) in several HST and JWST/NIRCam bands, allowing us to construct a state-of-the-art rest-frame UV-NIR Spectral Energy Distribution of the IHL for the first time at this high redshift. This allows stellar population characterisation of both the IHL and member galaxies, as well as the morphology distribution of group galaxies vs. their star-formation activity when coupled with Herschel data. We create a stacked image of the IHL, giving us a sensitivity to extended emission of 28.5 mag/arcsec2 at rest-frame 1um. We find that the IHL is extremely dust poor (Av~0), containing an evolved stellar population of log10(t50/yr)=8.8, corresponding to a formation epoch for 50% of the stellar material 0.63Gyr before z=1.85. There is no evidence of ongoing star-formation in the IHL. The IHL in this group at z=1.85 contributes ~10% of the total stellar mass, comparable with what is observed in local clusters. This suggests that the evolution of the IHL fraction is more self-similar with redshift than predicted by some models, challenging our understanding of IHL formation during the assembly of high-redshift clusters. JWST is unveiling a new side of group formation at this redshift, which will evolve into Virgo-like structures in the local Universe. △ Less

Submitted 7 June, 2023; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 14 pages + appendix, 11 figures, 4 tables. Accepted to A&A on 15th May 2023

Journal ref: A&A 677, A3 (2023)

ALMA FIR View of Ultra High-redshift Galaxy Candidates at $z\sim$ 11-17: Blue Monsters or Low-$z$ Red Interlopers?

Authors: Seiji Fujimoto, Steven L. Finkelstein, Denis Burgarella, Chris L. Carilli, Véronique Buat, Caitlin M. Casey, Laure Ciesla, Sandro Tacchella, Jorge A. Zavala, Gabriel Brammer, Yoshinobu Fudamoto, Masami Ouchi, Francesco Valentino, M. C. Cooper, Mark Dickinson, Maximilien Franco, Mauro Giavalisco, Taylor A. Hutchison, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Takashi Kojima, Rebecca L. Larson, Eric J. Murphy, Casey Papovich, Pablo G. Pérez-González , et al. (28 additional authors not shown)

Abstract: We present ALMA Band~7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866~$μ$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr… ▽ More We present ALMA Band~7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866~$μ$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr$^{-1}$. We detect a 5.1$σ$ line feature at $338.726\pm0.007$~GHz exactly coinciding with the JWST source position, with a 2\% likelihood of the signal being spurious. The most likely line identification would be [OIII]52$μ$m at $z=16.01$ or [CII]158$μ$m at $z=4.61$, whose line luminosities do not violate the non-detection of the dust continuum in both cases. Together with three other $z\gtrsim$ 11--13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high-$z$ solution at $z\sim$11--17 is favored in every candidate as a very blue (UV continuum slope of $\simeq-2.3$) and luminous ($M_{\rm UV}\simeq[-$24:$-21]$) system. Still, we find in several candidates that reasonable SED fits ($Δ$ $χ^{2}\lesssim4$) are reproduced by type-II quasar and/or quiescent galaxy templates with strong emission lines at $z\sim3$--5, where such populations predicted from their luminosity functions and EW([OIII]+H$β$) distributions are abundant in survey volumes used for the identification of the $z\sim$11--17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high-$z$ solutions, lower-$z$ possibilities are not completely ruled out in several of the $z\sim$11--17 candidates, indicating the need to consider the relative surface densities of the lower-$z$ contaminants in the ultra high-$z$ galaxy search. △ Less

Submitted 26 July, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 25 pages, 8 figures, 7 tables. ApJ in press

GOODS-ALMA 2.0: Last gigayear star formation histories of the so-called starbursts within the main sequence

Authors: L. Ciesla, C. Gómez-Guijarro, V. Buat, D. Elbaz, S. Jin, M. Béthermin, E. Daddi, M. Franco, H. Inami, G. Magdis, B. Magnelli

Abstract: Recently, a population of compact main sequence (MS) galaxies exhibiting starburst-like properties have been identified in the GOODS-ALMA blind survey at 1.1mm. Several evolution scenarios were proposed to explain their particular physical properties (e.g., compact size, low gas content, short depletion time). In this work, we aim at studying the star formation history (SFH) of the GOODS-ALMA gala… ▽ More Recently, a population of compact main sequence (MS) galaxies exhibiting starburst-like properties have been identified in the GOODS-ALMA blind survey at 1.1mm. Several evolution scenarios were proposed to explain their particular physical properties (e.g., compact size, low gas content, short depletion time). In this work, we aim at studying the star formation history (SFH) of the GOODS-ALMA galaxies to understand if the so-called ``starburst (SB) in the MS'' galaxies exhibit a different star formation activity over the last Gyr compared to MS galaxies that could explain their specificity. We use the CIGALE SED modelling code to which we add non-parametric SFHs. To compare quantitatively the recent SFH of the galaxies, we define a parameter, the star formation rate (SFR) gradient that provides the angle showing the direction that a galaxy has followed in the SFR vs stellar mass plane over a given period. We show that ``SB in the MS'' have positive or weak negative gradients over the last 100, 300, and 1000 Myr, at odds with a scenario where these galaxies would be transitioning from the SB region at the end of a strong starburst phase. Normal GOODS-ALMA galaxies and ``SB in the MS'' have the same SFR gradients distributions meaning that they have similar recent SFH, despite their different properties (compactness, low depletion time). The ``SBs in the MS'' manage to maintain a star-formation activity allowing them to stay within the MS. This points toward a diversity of galaxies within a complex MS. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: Submitted to A&A

Journal ref: A&A 672, A191 (2023)

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).XIV. The main sequence relation in a rich environment down to M_star ~ 10^6 Mo

Authors: A. Boselli, M. Fossati, J. Roediger, M. Boquien, M. Fumagalli, M. Balogh, S. Boissier, J. Braine, L. Ciesla, P. Côté, J. C. Cuillandre, L. Ferrarese, G. Gavazzi, S. Gwyn, Junais, G. Hensler, A. Longobardi, M. Sun

Abstract: Using a compilation of Halpha fluxes for 384 star forming galaxies detected during the VESTIGE survey, we study several important scaling relations for a complete sample of galaxies in a rich environment. The extraordinary sensitivity of the data allows us to sample the whole dynamic range of the Halpha luminosity function, from massive (M*~10^11 Mo) to dwarf systems (M*~10^6 Mo). This extends pre… ▽ More Using a compilation of Halpha fluxes for 384 star forming galaxies detected during the VESTIGE survey, we study several important scaling relations for a complete sample of galaxies in a rich environment. The extraordinary sensitivity of the data allows us to sample the whole dynamic range of the Halpha luminosity function, from massive (M*~10^11 Mo) to dwarf systems (M*~10^6 Mo). This extends previous works to a dynamic range in stellar mass and star formation rate (10^-4<SFR<10 Mo yr^-1) never explored so far. The main sequence (MS) relation derived for all star forming galaxies within one virial radius of the Virgo cluster has a slope comparable to that observed in other nearby samples of isolated objects, but has a dispersion ~3 times larger. The dispersion is tightly connected to the available amount of HI gas, with gas-poor systems located far below objects of similar stellar mass but with a normal HI content. When measured on unperturbed galaxies with a normal HI gas content, the relation has a slope a=0.92, an intercept b=-1.57, and a scatter ~0.40. We compare these observational results to the prediction of models. The observed scatter in the MS relation can be reproduced only after a violent and active stripping process such as ram-pressure that removes gas from the disc and quenches star formation on short (<1 Gyr) timescales. This rules out milder processes such as starvation. This interpretation is also consistent with the position of galaxies of different star formation activity and gas content within the phase-space diagram. We also show that the star forming regions formed in the stripped material outside perturbed galaxies are located well above the MS relation drawn by unperturbed systems. These HII regions, which might be at the origin of compact sources typical in rich environments, are living a starburst phase lasting only <50 Myr, later becoming quiescent systems. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: Accepted for publication on A&A

Journal ref: A&A 669, A73 (2023)

The hidden side of cosmic star formation at z > 3: Bridging optically-dark and Lyman break galaxies with GOODS-ALMA

Authors: Mengyuan Xiao, David Elbaz, Carlos Gómez-Guijarro, Lucas Leroy, Longji Bing, Emanuele Daddi, Benjamin Magnelli, Maximilien Franco, Luwenjia Zhou, Mark Dickinson, Tao Wang, Wiphu Rujopakarn, Georgios E. Magdis, Ezequiel Treister, Hanae Inami, Ricardo Demarco, Mark T. Sargent, Xinwen Shu, Jeyhan S. Kartaltepe, David M. Alexander, Matthieu Béthermin, Frederic Bournaud, Laure Ciesla, Henry C. Ferguson, Steven L. Finkelstein , et al. (15 additional authors not shown)

Abstract: Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order… ▽ More Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift. △ Less

Submitted 10 February, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

Journal ref: A&A 672, A18 (2023)

Dusty Starbursts Masquerading as Ultra-high Redshift Galaxies in JWST CEERS Observations

Authors: Jorge A. Zavala, Veronique Buat, Caitlin M. Casey, Denis Burgarella, Steven L. Finkelstein, Micaela B. Bagley, Laure Ciesla, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Maximilien Franco, E. F. Jim'enez-Andrade, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Aurélien Le Bail, E. J. Murphy, Casey Papovich, Sandro Tacchella, Stephen M. Wilkins, Itziar Aretxaga, Peter Behroozi, Jaclyn B. Champagne, Adriano Fontana, Mauro Giavalisco, Andrea Grazian , et al. (99 additional authors not shown)

Abstract: Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may als… ▽ More Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations. △ Less

Submitted 30 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: Published in The Astrophysical Journal Letters (updated to match the published version)

Journal ref: ApJL 943 L9 (2023)

A Long Time Ago in a Galaxy Far, Far Away: A Candidate z ~ 12 Galaxy in Early JWST CEERS Imaging

Authors: Steven L. Finkelstein, Micaela B. Bagley, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Jeyhan S. Kartaltepe, Casey Papovich, Denis Burgarella, Dale D. Kocevski, Marc Huertas-Company, Kartheik G. Iyer, Rebecca L. Larson, Pablo G. Pérez-González, Caitlin Rose, Sandro Tacchella, Stephen M. Wilkins, Katherine Chworowsky, Aubrey Medrano, Alexa M. Morales, Rachel S. Somerville, L. Y. Aaron Yung, Adriano Fontana, Mauro Giavalisco, Andrea Grazian, Norman A. Grogin , et al. (95 additional authors not shown)

Abstract: We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging f… ▽ More We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang. △ Less

Submitted 7 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

The Physical Conditions of Emission-Line Galaxies at Cosmic Dawn from JWST/NIRSpec Spectroscopy in the SMACS 0723 Early Release Observations

Authors: Jonathan R. Trump, Pablo Arrabal Haro, Raymond C. Simons, Bren E. Backhaus, Ricardo O. Amorín, Mark Dickinson, Vital Fernández, Casey Papovich, David C. Nicholls, Lisa J. Kewley, Samantha W. Brunker, John J. Salzer, Stephen M. Wilkins, Omar Almaini, Micaela B. Bagley, Danielle A. Berg, Rachana Bhatawdekar, Laura Bisigello, Véronique Buat, Denis Burgarella, Antonello Calabrò, Caitlin M. Casey, Laure Ciesla, Nikko J. Cleri, Justin W. Cole , et al. (39 additional authors not shown)

Abstract: We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wav… ▽ More We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain \textit{absolute} spectrophotometry of the current version of the reductions. Compared to $z\sim3$ galaxies in the literature, the $z>5$ galaxies have similar [OIII]$λ$5008/H$β$ ratios, similar [OIII]$λ$4364/H$γ$ ratios, and higher ($\sim$0.5 dex) [NeIII]$λ$3870/[OII]$λ$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$λ$3870/[OII]$λ$3728, [OIII]$λ$4364/H$γ$, and [OIII]$λ$5008/H$β$ emission-line ratios are consistent with an interstellar medium that has very high ionization ($\log(Q) \simeq 8-9$, units of cm~s$^{-1}$), low metallicity ($Z/Z_\odot \lesssim 0.2$), and very high pressure ($\log(P/k) \simeq 8-9$, units of cm$^{-3}$). The combination of [OIII]$λ$4364/H$γ$ and [OIII]$λ$(4960+5008)/H$β$ line ratios indicate very high electron temperatures of $4.1<\log(T_e/{\rm K})<4.4$, further implying metallicities of $Z/Z_\odot \lesssim 0.2$ with the application of low-redshift calibrations for ``$T_e$-based'' metallicities. These observations represent a tantalizing new view of the physical conditions of the interstellar medium in galaxies at cosmic dawn. △ Less

Submitted 19 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

Spatial disconnection between stellar and dust emissions: the test of the Antennae Galaxies (Arp 244)

Authors: L. -M. Seillé, V. Buat, W. Haddad, A. Boselli, M. Boquien, L. Ciesla, Y. Roehlly, D. Burgarella

Abstract: The detection with of the Atacama Large Millimeter Array (ALMA) of dust-rich high redshift galaxies whose cold dust emission is spatially disconnected from the ultraviolet emission bears a challenge for modelling their spectral energy distributions (SED) with codes based on an energy budget between the stellar and dust components. We test the validity of energy balance modelling on a nearby resolv… ▽ More The detection with of the Atacama Large Millimeter Array (ALMA) of dust-rich high redshift galaxies whose cold dust emission is spatially disconnected from the ultraviolet emission bears a challenge for modelling their spectral energy distributions (SED) with codes based on an energy budget between the stellar and dust components. We test the validity of energy balance modelling on a nearby resolved galaxy with vastly different ultraviolet and infrared spatial distributions and infer what information can be reliably retrieved from the analysis of the full spectral energy distribution. We use 15 broadband images of the Antennae Galaxies ranging from far-ultraviolet to far-infrared and divide Arp 244 into 58 square ~1 kpc$^2$ regions. We fit the data with CIGALE to determine the star formation rate, stellar mass and dust attenuation of each region. We compare these quantities for the addition of the 58 regions to the ones obtained for Arp 244 as a whole and find that both estimates are consistent within one sigma. We present the spatial distribution of these physical parameters as well as the shape of the attenuation curve across the Antennae Galaxies . We also observe a flattening of the attenuation curves with increasing attenuation and dust surface density in agreement with the predictions of hydrodynamical simulations coupled with radiative transfer modelling. △ Less

Submitted 16 July, 2022; originally announced July 2022.

Comments: 13 pages, 9 figures, Accepted for publication in Astronomy & Astrophysics

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

Starbursts with suppressed velocity dispersion revealed in a forming cluster at z=2.51

Authors: Mengyuan Xiao, Tao Wang, David Elbaz, Daisuke Iono, Xing Lu, Longji Bing, Emanuele Daddi, Benjamin Magnelli, Carlos Gómez-Guijarro, Frederic Bournaud, Qiusheng Gu, Shuowen Jin, Francesco Valentino, Anita Zanella, Raphael Gobat, Sergio Martin, Gabriel Brammer, Kotaro Kohno, Corentin Schreiber, Laure Ciesla, Xiaoling Yu, Koryo Okumura

Abstract: One of the most prominent features of galaxy clusters is the presence of a dominant population of massive ellipticals in their cores. Stellar archaeology suggests that these gigantic beasts assembled most of their stars in the early Universe via starbursts. However, the role of dense environments and their detailed physical mechanisms in triggering starburst activities remain unknown. Here we repo… ▽ More One of the most prominent features of galaxy clusters is the presence of a dominant population of massive ellipticals in their cores. Stellar archaeology suggests that these gigantic beasts assembled most of their stars in the early Universe via starbursts. However, the role of dense environments and their detailed physical mechanisms in triggering starburst activities remain unknown. Here we report spatially resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO $J= 3-2$ emission line, with a resolution of about 2.5 kiloparsecs, toward a forming galaxy cluster core with starburst galaxies at $z=2.51$. In contrast to starburst galaxies in the field often associated with galaxy mergers or highly turbulent gaseous disks, our observations show that the two starbursts in the cluster exhibit dynamically cold (rotation-dominated) gas-rich disks. Their gas disks have extremely low velocity dispersion ($σ_{\mathrm{0}} \sim 20-30$ km s$^{-1}$), which is three times lower than their field counterparts at similar redshifts. The high gas fraction and suppressed velocity dispersion yield gravitationally unstable gas disks, which enables highly efficient star formation. The suppressed velocity dispersion, likely induced by the accretion of corotating and coplanar cold gas, might serve as an essential avenue to trigger starbursts in massive halos at high redshifts. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 19 pages, 14 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 664, A63 (2022)

Comparison of the star formation of X-ray selected active galactic nuclei in eFEDS with star-forming galaxies

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla, K. Malek, R. Shirley

Abstract: We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with… ▽ More We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with $\sim 17000$ sources. Both samples have the same photometric coverage, from optical to far-infrared. We construct the spectral energy distributions (SEDs) of all sources and fit them using the CIGALE code, utilizing the same templates and parametric grid for both samples. We account for the mass incompleteness and exclude quiescent systems from both datasets. These allow us to compare the SFR of the two populations in a uniform manner, minimising systematic effects. Based on our analysis, AGN at low and moderate L$_X$ ($\rm L_{X,2-10keV} < 10^{44}\,ergs^{-1}$), have SFR that is lower, or at most, equal to that of star-forming galaxies, in agreement with previous studies. The large number of luminous X-ray AGN available in the eFEDS dataset, enable us to expand our investigations at higher L$_X$ to test previous, tentative results. At $\rm L_{X,2-10keV} > 10^{44.2}\,ergs^{-1}$, the SFR of AGN appears enhanced, by $\sim 30\%$, compared to that of star-forming sources, for systems with stellar mass, $\rm 10.5 < log\,[M_*(M_\odot)] < 11.5$, confirming indications found in previous studies. The most massive sources $\rm log\,[M_*(M_\odot)] > 11.5$, present a flat SFR$_{norm}$-L$_X$ relation up to $\rm L_{X,2-10keV} \sim 10^{44.5}\,ergs^{-1}$, with SFR similar to that of star-forming galaxies. However, at higher L$_X$ ($\rm L_{X,2-10keV} \sim 10^{45}\,ergs^{-1}$), we find indications that the SFR of these massive AGN hosts may be enhanced compared to that of non-AGN systems. △ Less

Submitted 5 May, 2022; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 14 pages, 7 figures, 2 tables

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

Fitting AGN/galaxy X-ray-to-radio SEDs with CIGALE and improvement of the code

Authors: Guang Yang, Médéric Boquien, William N. Brandt, Véronique Buat, Denis Burgarella, Laure Ciesla, Bret D. Lehmer, Katarzyna E. Małek, George Mountrichas, Casey Papovich, Estelle Pons, Marko Stalevski, Patrice Theulé, Shifu Zhu

Abstract: Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Re… ▽ More Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of CIGALE (named X-CIGALE) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply X-CIGALE to different samples, including COSMOS spectroscopic type 2 AGNs, CDF-S X-ray detected normal galaxies, SDSS quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of X-CIGALE and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as $L_X(θ) \propto 1+\cos θ$, where $θ$ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of CIGALE, and publicly release this new version as CIGALE v2022.0 on https://cigale.lam.fr △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: ApJ accepted. The code is available at https://cigale.lam.fr/news/

GOODS-ALMA 2.0: Starbursts in the main sequence reveal compact star formation regulating galaxy evolution prequenching

Authors: C. Gómez-Guijarro, D. Elbaz, M. Xiao, V. I. Kokorev, G. E. Magdis, B. Magnelli, E. Daddi, F. Valentino, M. T. Sargent, M. Dickinson, M. Béthermin, M. Franco, A. Pope, B. S. Kalita, L. Ciesla, R. Demarco, H. Inami, W. Rujopakarn, X. Shu, T. Wang, L. Zhou, D. M. Alexander, F. Bournaud, R. Chary, H. C. Ferguson , et al. (16 additional authors not shown)

Abstract: Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived… ▽ More Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Comments: Accepted by A&A. 26 pages, 13 figures

Journal ref: A&A 659, A196 (2022)

The Evolving Interstellar Medium of Star-Forming Galaxies, as traced by $\texttt{Stardust}$

Authors: Vasily Kokorev, Georgios Magdis, Iary Davidzon, Gabriel Brammer, Francesco Valentino, Emanuele Daddi, Laure Ciesla, Daizhong Liu, Shuowen Jin, Isabella Cortzen, Ivan Delvecchio, Clara Giménez-Arteaga, Carlos Gómez-Guijarro, Mark Sargent, Sune Toft, John R. Weaver

Abstract: We analyse the far-infrared properties of $\sim$ 5,000 star-forming galaxies at $z<4.5$, drawn from the deepest, super-deblended catalogues in the GOODS-N and COSMOS fields. We develop a novel panchromatic SED fitting algorithm, $\texttt{Stardust}$, that models the emission from stars, AGN, and infrared dust emission, without relying on energy balance assumptions. Our code provides robust estimate… ▽ More We analyse the far-infrared properties of $\sim$ 5,000 star-forming galaxies at $z<4.5$, drawn from the deepest, super-deblended catalogues in the GOODS-N and COSMOS fields. We develop a novel panchromatic SED fitting algorithm, $\texttt{Stardust}$, that models the emission from stars, AGN, and infrared dust emission, without relying on energy balance assumptions. Our code provides robust estimates of the UV-optical and FIR physical parameters, such as the stellar mass ($M_*$), dust mass ($M_{\rm dust}$), infrared luminosities ($L_{\rm IR}$) arising from AGN and star formation activity, and the average intensity of the interstellar radiation field ($\langle U \rangle$). Through a set of simulations we quantify the completeness of our data in terms of $M_{\rm dust}$, $L_{\rm IR}$ and $\langle U \rangle$, and subsequently characterise the distribution and evolution of these parameters with redshift. We focus on the dust-to-stellar mass ratio ($f_{\rm dust}$), which we parametrise as a function of cosmic age, stellar mass, and specific star formation rate. The $f_{\rm dust}$ is found to increase by a factor of 10 from $z=0$ to $z=2$ and appears to remain flat at higher$-z$, mirroring the evolution of the gas fraction. We also find a growing fraction of warm to cold dust with increasing distance from the main sequence, indicative of more intense interstellar radiation fields, higher star formation efficiencies and more compact star forming regions for starburst galaxies. Finally, we construct the dust mass functions (DMF) of star-forming galaxies up to $z=1$ by transforming the stellar mass function to DMF through the scaling relations derived here. The evolution of $f_{\rm dust}$ and the recovered DMFs are in good agreement with the theoretical predictions of the Horizon-AGN and IllustrisTNG simulations. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 39 pages, 21 figures, 3 tables. Accepted to ApJ July 10, 2021

Polar dust obscuration in broad-line active galaxies from the XMM-XXL field

Authors: V. Buat, G. Mountrichas, G. Yang, M. Boquien, Y. Roehlly, D. Burgarella, M. Stalevski, L. Ciesla, P. Theulé

Abstract: Dust is observed in the polar regions of nearby AGN and it is known to contribute substantially to their mid-IR emission and to the obscuration of their UV to optical emission. We selected a sample of 1275 BLAGN in the XMM-XXL field, with optical to infrared photometric data. These AGN are seen along their polar direction and we expect a maximal impact of dust located around the poles when it is p… ▽ More Dust is observed in the polar regions of nearby AGN and it is known to contribute substantially to their mid-IR emission and to the obscuration of their UV to optical emission. We selected a sample of 1275 BLAGN in the XMM-XXL field, with optical to infrared photometric data. These AGN are seen along their polar direction and we expect a maximal impact of dust located around the poles when it is present. We used X-CIGALE, which introduces a dust component to account for obscuration along the polar directions, modeled as a foreground screen, and an extinction curve that is chosen as it steepens significantly at short wavelengths or is much grayer. By comparing the results of different fits, we are able to define subsamples of sources with positive statistical evidence in favor of or against polar obscuration and described using the gray or steep extinction curve. We find a similar fraction of sources with positive evidence for and against polar dust. Applying statistical corrections, we estimate that half of our sample could contain polar dust and among them, 60% exhibit a steep extinction curve and 40% a flat extinction curve; although these latter percentages are found to depend on the adopted extinction curves. The obscuration in the V-band is not found to correlate with the X-ray column density, while A_V/N_H ratios span a large range of values and higher dust temperatures are found with the flat, rather than with the steep extinction curve. Ignoring this polar dust component in the fit of the spectral energy distribution of these composite systems leads to an overestimation of the stellar contribution. A single fit with a polar dust component described with an SMC extinction curve efficiently overcomes this issue but it fails at identifying all the AGN with polar dust obscuration. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: 12 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

GOODS-ALMA 2.0: Source catalog, number counts, and prevailing compact sizes in 1.1 mm galaxies

Authors: C. Gómez-Guijarro, D. Elbaz, M. Xiao, M. Béthermin, M. Franco, B. Magnelli, E. Daddi, M. Dickinson, R. Demarco, H. Inami, W. Rujopakarn, G. E. Magdis, X. Shu, R. Chary, L. Zhou, D. M. Alexander, F. Bournaud, L. Ciesla, H. C. Ferguson, S. L. Finkelstein, M. Giavalisco, D. Iono, S. Juneau, J. S. Kartaltepe, G. Lagache , et al. (14 additional authors not shown)

Abstract: Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas… ▽ More Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $σ= 68.4μ$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks. △ Less

Submitted 21 December, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: Accepted by A&A. 30 pages, 18 figures

Journal ref: A&A 658, A43 (2022)

The role of AGN and its obscuration on the position of the host galaxy relative to Main Sequence

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla, K. Malek, R. Shirley

Abstract: We use X-ray Active Galactic Nuclei (AGN) observed by the Chandra X-ray Observatory within the 9.3 deg$^2$ Bo$\rm \ddot{o}$tes field of the NDWFS to study whether there is a correlation between X-ray luminosity (L$_X$) and star formation rate (SFR) of the host galaxy, at $\rm 0.5<z<2.0$, with respect to the position of the galaxy to the main sequence (SFR$_{norm}$). About half of the sources in th… ▽ More We use X-ray Active Galactic Nuclei (AGN) observed by the Chandra X-ray Observatory within the 9.3 deg$^2$ Bo$\rm \ddot{o}$tes field of the NDWFS to study whether there is a correlation between X-ray luminosity (L$_X$) and star formation rate (SFR) of the host galaxy, at $\rm 0.5<z<2.0$, with respect to the position of the galaxy to the main sequence (SFR$_{norm}$). About half of the sources in the X-ray sample have spectroscopic redshifts. We also construct a reference galaxy catalogue. For both datasets, we use photometric data from optical to the far infrared, compiled by the HELP project and apply spectral energy distribution (SED) fitting, using the X-CIGALE code. We exclude quiescent sources from both the X-ray and the reference samples. We also account for the mass completeness of our dataset, in different redshifts bins. Our analysis highlights the importance of studying the SFR-L$_X$ relation, in a uniform manner, taking into account the systematics and selection effects. Our results suggest that, in less massive galaxies ($\rm log\,[M_*(M_\odot)] \sim 11$), AGN enhances the SFR of the host galaxy by $\sim 50\%$ compared to non AGN systems. A flat relation is observed for the most massive galaxies. SFR$_{norm}$ does not evolve with redshift. The results, although tentative, are consistent with a scenario in which, in less massive systems, both AGN and star formation (SF) are fed by cold gas, supplied by a merger event. In more massive galaxies, the flat relation could be explained by a different SMBH fuelling mechanism that is decoupled from the star formation of the host galaxy (e.g. hot diffuse gas). Finally, we compare the host galaxy properties of X-ray absorbed and unabsorbed sources. Our results show no difference which suggests that X-ray absorption is not linked with the properties of the galaxy. △ Less

Submitted 20 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A

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

Investigating the delay between dust radiation and star-formation in local and distant quenching galaxies

Authors: Laure Ciesla, Veronique Buat, Mederic Boquien, Alessandro Boselli, David Elbaz, Gregoire Aufort

Abstract: We investigate the timescale with which the IR luminosity decreases after a complete and rapid quenching of star formation using observations of local and high-redshift galaxies. From SED modelling, we derive the time since quenching of a subsample of 14 galaxies from the Herschel Reference Survey suffering from ram-pressure stripping due to the environment of the Virgo cluster and of a subsample… ▽ More We investigate the timescale with which the IR luminosity decreases after a complete and rapid quenching of star formation using observations of local and high-redshift galaxies. From SED modelling, we derive the time since quenching of a subsample of 14 galaxies from the Herschel Reference Survey suffering from ram-pressure stripping due to the environment of the Virgo cluster and of a subsample of 7 rapidly quenched COSMOS galaxies selected through a state-of-the-art statistical method already tested on the determination of galaxies' star formation history. Three out of the 7 COSMOS galaxies have an optical spectra with no emission line, confirming their quenched nature. Present physical properties of the two samples are obtained as well as the past L$_{IR}$ of these galaxies, just before their quenching, from the long-term SFH properties. This past L$_{IR}$ is shown to be consistent with the L$_{IR}$ of reference samples of normally star-forming galaxies with same $M_*$ and $z$ than each of our quenched galaxies. We put constraints on the present to past L$_{IR}$ ratio as a function of quenching time. The two samples probe different dynamical ranges in terms of quenching age with the HRS galaxies exhibiting longer timescales (0.2-3\,Gyr) compared to the COSMOS one ($<100$\,Myr). Assuming an exponential decrease of the L$_{IR}$ after quenching, the COSMOS quenched galaxies are consistent with short e-folding times less than a couple of hundreds of Myr while the properties of the HRS quenched galaxies are compatible with timescales of several hundreds of Myr. For the HRS sample, this result is consistent with ram pressure stripping due to the environment. For the COSMOS sample, different quenching processes are acting on short to intermediate timescales. Processes such as galaxy mergers, disk instabilities or environmental effects can produce such strong star formation variability. △ Less

Submitted 31 May, 2021; originally announced June 2021.

Comments: 14 pages, 14 figures, accepted for publication in A&A

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

Multiwavelength dissection of a massive heavily dust-obscured galaxy and its blue companion at z $\sim $2

Authors: M. Hamed, L. Ciesla, M. Béthermin, K. Małek, E. Daddi, M. T. Sargent, R. Gobat

Abstract: In this work we study a system of two galaxies, Astarte and Adonis, at z $\sim $2 when the Universe was undergoing its peak of star formation activity. Astarte is a dusty star-forming galaxy at the massive-end of the main sequence (MS) and Adonis is a less-massive, bright in ultraviolet (UV), companion galaxy with an optical spectroscopic redshift. We analyse the physical properties of this system… ▽ More In this work we study a system of two galaxies, Astarte and Adonis, at z $\sim $2 when the Universe was undergoing its peak of star formation activity. Astarte is a dusty star-forming galaxy at the massive-end of the main sequence (MS) and Adonis is a less-massive, bright in ultraviolet (UV), companion galaxy with an optical spectroscopic redshift. We analyse the physical properties of this system, and probe the gas mass of Astarte with its ALMA CO emission, to investigate whether this ultra-massive galaxy is quenching or not. We use CIGALE - a spectral energy distribution modeling code - to derive the key physical properties of Astarte and Adonis, mainly their star formation rates (SFRs), stellar masses, and dust luminosities. We inspect the variation of the physical parameters depending on the assumed dust attenuation law. We also estimate the molecular gas mass of Astarte from its CO emission, using different $α_{CO}$ and transition ratios ($r_{31}$) and discuss the implication of the various assumptions on the gas mass derivation. We find that Astarte exhibits a MS-like star formation activity, while Adonis is undergoing a strong starburst (SB) phase. The molecular gas mass of Astarte is far below the gas fraction of typical star-forming galaxies at z=2. This low gas content and high SFR, result in a depletion time of $0.22\pm0.07$ Gyrs, slightly shorter than what is expected for a MS galaxy at this redshift. The CO luminosity versus the total IR luminosity suggests a MS-like activity assuming a galactic conversion factor and a low transition ratio. The SFR of Astarte is of the same order using different attenuation laws, unlike its stellar mass that increases using shallow attenuation laws. We discuss these properties and suggest that Astarte might be experiencing a recent decrease of star formation activity and is quenching through the MS following a SB epoch. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 10 pages, 7 figures, accepted for publication in A&A. Abstract abridged for arXiv submission

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

X-ray flux in the SED modelling: An application of X-CIGALE in the XMM-XXL field

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla

Abstract: X-CIGALE, built upon the spectral energy distribution (SED) code of CIGALE, implements important new features: The code accounts for obscuring material in the polars of the AGN and has the ability to fit X-ray fluxes. In this work, we use ~2500 spectroscopic, X-ray AGN from the XMM-XXL-North field and examine the improvements the new features bring in the SED modelling analysis. Based on our resul… ▽ More X-CIGALE, built upon the spectral energy distribution (SED) code of CIGALE, implements important new features: The code accounts for obscuring material in the polars of the AGN and has the ability to fit X-ray fluxes. In this work, we use ~2500 spectroscopic, X-ray AGN from the XMM-XXL-North field and examine the improvements the new features bring in the SED modelling analysis. Based on our results, X-CIGALE successfully connects the X-ray with the UV luminosity in the whole range spanned by our sample (log Lx(2-10 keV) = (42 - 46) erg/s). The addition of the new features globally improves the efficiency of X-CIGALE in the estimation and characterization of the AGN component. The classification into type 1 and type 2 based on their inclination angle is improved, especially at redshifts lower than 1. The statistical significance of the AGN fraction, fracAGN, measurements is increased, in particular for luminous X-ray sources (log LX > 45 erg/s). These conclusions hold under the condition that (mid-) IR photometry is available in the SED fitting process. The addition of polar dust increases the AGN fraction and the efficiency of the SED decomposition to detect AGN among X-ray selected sources. X-CIGALE estimates a strong AGN (fracAGN > 0.3) in more than 90% of the infrared selected AGN and 75% of X-ray detected AGN not selected by IR colour criteria. The latter drops to ~50% when polar dust is not included. The ability of X-CIGALE to include X-ray information in the SED fitting process can be instrumental in the optimal exploitation of the wealth of data that current (eROSITA) and future (ATHENA) missions will provide us. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: A&A accepted, 19 pages, 20 figures, 5 tables

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

GOODS-ALMA: Optically dark ALMA galaxies shed light on a cluster in formation at z = 3.5

Authors: L. Zhou, D. Elbaz, M. Franco, B. Magnelli, C. Schreiber, T. Wang, L. Ciesla, E. Daddi, M. Dickinson, N. Nagar, G. Magdis, D. M. Alexander, M. Béthermin, R. Demarco, J. Mullaney, F. Bournaud, H. Ferguson, S. L. Finkelstein, M. Giavalisco, H. Inami, D. Iono, S. Juneau, G. Lagache, H. Messias, K. Motohara , et al. (8 additional authors not shown)

Abstract: In this paper we study the properties of the six optically dark galaxies detected in the 69 arcmin^2 GOODS-ALMA 1.1mm continuum survey. While none of them are listed in the deepest H-band based CANDELS catalog in the GOODS-South field down to H=28.16AB, we were able to de-blend two of them from their bright neighbor and measure an $H$-band flux for them. We note that AGS4 and AGS15 have H=25.23, 2… ▽ More In this paper we study the properties of the six optically dark galaxies detected in the 69 arcmin^2 GOODS-ALMA 1.1mm continuum survey. While none of them are listed in the deepest H-band based CANDELS catalog in the GOODS-South field down to H=28.16AB, we were able to de-blend two of them from their bright neighbor and measure an $H$-band flux for them. We note that AGS4 and AGS15 have H=25.23, 27.11AB respectively. Their extreme proximity (0.50", 0.27") to a bright optical source and their extreme faintness prevented them from being included in the H-band catalog. We present the spectroscopic scan follow-up of five of the six sources with ALMA band 4. All are detected in the 2mm continuum with signal-to-noise ratios higher than eight. One emission line is detected in AGS4 (ν_{obs} =151.44GHz with a S/N=8.58) and AGS17 (ν_{obs} =154.78GHz with a S/N=10.23), which we interpret in both cases as being due to the CO(6-5) line at z^{AGS4}_{spec}=3.556 and z^{AGS4}_{spec}=3.467, respectively. These redshifts match both the probability distribution of the photometric redshifts derived from the UV to near-infrared spectral energy distributions (SEDs) and the far-infrared SEDs for typical dust temperatures of galaxies at these redshifts. We present evidence that nearly 70% (4/6 of galaxies) of the optically dark galaxies belong to the same overdensity of galaxies at z~3.5. overdensity The most massive one, AGS24 (M_{\star} = 10^{11.32^{+0.02}_{-0.19}} M_{\odot}), is the most massive galaxy without an active galactic nucleus (AGN) at $z$\,>\,3 in the GOODS-ALMA field. It falls in the very center of the peak of the galaxy surface density, which suggests that the surrounding overdensity is a proto-cluster in the process of virialization and that AGS24 is the candidate progenitor of the future brightest cluster galaxy (BCG). △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: 17 pages, 13 figures, 3 tables, accepted by A&A

Journal ref: A&A 642, A155 (2020)

Space Project for Astrophysical and Cosmological Exploration (SPACE), an ESA stand-alone mission and a possible contribution to the Origins Space Telescope

Authors: Denis Burgarella, Andrew Bunker, Rychard Bouwens, Laurent Pagani, Jose Afonso, Hakim Atek, Marc Audard, Sylvie Cabrit, Karina Caputi, Laure Ciesla, Christopher Conselice, Asantha Cooray, Giovanni Cresci, Mirko Curti, Jose Miguel Rodriguez Espinosa, Marc Ferrari, Chiaki Kobayashi, Nadege Lagarde, Jesus Gallego Maestro, Roberto Maiolino, Katarzyna Malek, Filippo Mannucci, Julien Montillaud, Pascal Oesch, Chris Pearson , et al. (8 additional authors not shown)

Abstract: We propose a new mission called Space Project for Astrophysical and Cosmological Exploration (SPACE) as part on the ESA long term planning Voyage 2050 programme. SPACE will study galaxy evolution at the earliest times, with the key goals of charting the formation of the heavy elements, measuring the evolution of the galaxy luminosity function, tracing the build-up of stellar mass in galaxies over… ▽ More We propose a new mission called Space Project for Astrophysical and Cosmological Exploration (SPACE) as part on the ESA long term planning Voyage 2050 programme. SPACE will study galaxy evolution at the earliest times, with the key goals of charting the formation of the heavy elements, measuring the evolution of the galaxy luminosity function, tracing the build-up of stellar mass in galaxies over cosmic time, and finding the first super-massive black holes (SMBHs) to form. The mission will exploit a unique region of the parameter space, between the narrow ultra-deep surveys with HST and JWST, and shallow wide-field surveys such as Roman Space Telescope and EUCLID, and should yield by far the largest sample of any current or planned mission of very high redshift galaxies at z > 10 which are sufficiently bright for detailed follow-up spectroscopy. Crucially, we propose a wide-field spectroscopic near-IR + mid-IR capability which will greatly enhance our understanding of the first galaxies by detecting and identifying a statistical sample of the first galaxies and the first SMBH, and to chart the metal enrichment history of galaxies in the early Universe - potentially finding signatures of the very first stars to form from metal-free primordial gas. The wide-field and wavelength range of SPACE will also provide us a unique opportunity to study star formation by performing a wide survey of the Milky Way in the near-IR + mid-IR. This science project can be enabled either by a stand-alone ESA-led M mission or by an instrument for an L mission (with ESA and/or NASA, JAXA and other international space agencies) with a wide-field (sub-)millimetre capability at wavelength > 500 microns. △ Less

Submitted 5 August, 2020; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: Submitted to Experimental Astronomy as one of the proposals for the European Space Agency Voyage 2050 Programme

GOODS-ALMA: The slow downfall of star-formation in $z$ = 2-3 massive galaxies

Authors: M. Franco, D. Elbaz, L. Zhou, B. Magnelli, C. Schreiber, L. Ciesla, M. Dickinson, N. Nagar, G. Magdis, D. M. Alexander, M. Béthermin, R. Demarco, E. Daddi, T. Wang, J. Mullaney, M. Sargent, H. Inami, X. Shu, F. Bournaud, R. Chary, R. T. Coogan, H. Ferguson, S. L. Finkelstein, M. Giavalisco, C. Gómez-Guijarro , et al. (12 additional authors not shown)

Abstract: We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA a… ▽ More We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2. △ Less

Submitted 7 October, 2020; v1 submitted 6 May, 2020; originally announced May 2020.

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

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

GOODS-ALMA: Using IRAC and VLA to probe fainter millimeter galaxies

Authors: M. Franco, D. Elbaz, L. Zhou, B. Magnelli, C. Schreiber, L. Ciesla, M. Dickinson, N. Nagar, G. Magdis, D. M. Alexander, M. Béthermin, R. Demarco, E. Daddi, T. Wang, J. Mullaney, H. Inami, X. Shu, F. Bournaud, R. Chary, R. T. Coogan, H. Ferguson, S. L. Finkelstein, M. Giavalisco, C. Gómez-Guijarro, D. Iono , et al. (11 additional authors not shown)

Abstract: In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $σ$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $μ$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$σ$ at which the numbe… ▽ More In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $σ$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $μ$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$σ$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$σ$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues. △ Less

Submitted 7 October, 2020; v1 submitted 6 May, 2020; originally announced May 2020.

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

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

Constraining the recent star formation history of galaxies : an Approximate Bayesian Computation approach

Authors: G. Aufort, L. Ciesla, P. Pudlo, V. Buat

Abstract: [Abridged] Although galaxies are found to follow a tight relation between their star formation rate and stellar mass, they are expected to exhibit complex star formation histories (SFH), with short-term fluctuations. The goal of this pilot study is to present a method that will identify galaxies that are undergoing a strong variation of star formation activity in the last tens to hundreds Myr. In… ▽ More [Abridged] Although galaxies are found to follow a tight relation between their star formation rate and stellar mass, they are expected to exhibit complex star formation histories (SFH), with short-term fluctuations. The goal of this pilot study is to present a method that will identify galaxies that are undergoing a strong variation of star formation activity in the last tens to hundreds Myr. In other words, the proposed method will determine whether a variation in the last few hundreds of Myr of the SFH is needed to properly model the SED rather than a smooth normal SFH. To do so, we analyze a sample of COSMOS galaxies using high signal-to-noise ratio broad band photometry. We apply Approximate Bayesian Computation, a state-of-the-art statistical method to perform model choice, associated to machine learning algorithms to provide the probability that a flexible SFH is preferred based on the observed flux density ratios of galaxies. We present the method and test it on a sample of simulated SEDs. The input information fed to the algorithm is a set of broadband UV to NIR (rest-frame) flux ratios for each galaxy. The method has an error rate of 21% in recovering the right SFH and is sensitive to SFR variations larger than 1 dex. A more traditional SED fitting method using CIGALE is tested to achieve the same goal, based on fits comparisons through Bayesian Information Criterion but the best error rate obtained is higher, 28%. We apply our new method to the COSMOS galaxies sample. The stellar mass distribution of galaxies with a strong to decisive evidence against the smooth delayed-$τ$ SFH peaks at lower M* compared to galaxies where the smooth delayed-$τ$ SFH is preferred. We discuss the fact that this result does not come from any bias due to our training. Finally, we argue that flexible SFHs are needed to be able to cover that largest SFR-M* parameter space possible. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Journal ref: A&A 635, A136 (2020)

X-CIGALE: fitting AGN/galaxy SEDs from X-ray to infrared

Authors: Guang Yang, M. Boquien, V. Buat, D. Burgarella, L. Ciesla, F. Duras, M. Stalevski, W. N. Brandt, C. Papovich

Abstract: CIGALE is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of CIGALE is not able to fit X-ray data, which often provide unique insights into AGN intrinsic power. We develop a new X-ray module for CIGALE, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of CIGALE from UV-to-IR w… ▽ More CIGALE is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of CIGALE is not able to fit X-ray data, which often provide unique insights into AGN intrinsic power. We develop a new X-ray module for CIGALE, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of CIGALE from UV-to-IR wavelengths. We implement a modern clumpy two-phase torus model, SKIRTOR. To account for moderately extincted type 1 AGNs, we implement polar-dust extinction. We publicly release the source code (named X-CIGALE). We test X-CIGALE with X-ray detected AGNs in SDSS, COSMOS, and AKARI-NEP. The fitting quality (as indicated by reduced $χ^2$) is good in general, indicating that X-CIGALE is capable of modelling the observed SED from X-ray to IR. We discuss constrainability and degeneracy of model parameters in the fitting of AKARI-NEP, for which excellent mid-IR photometric coverage is available. We also test fitting a sample of AKARI-NEP galaxies for which only X-ray upper limits are available from Chandra observations, and find that the upper limit can effectively constrain the AGN SED contribution for some systems. Finally, using X-CIGALE, we assess the ability of Athena to constrain the AGN activity in future extragalactic studies. △ Less

Submitted 22 January, 2020; originally announced January 2020.

Comments: Published in MNRAS. The source code is temporarily available at https://gitlab.lam.fr/gyang/cigale/tree/xray, but will be eventually integrated into CIGALE (https://cigale.lam.fr/) in its next release

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 491, Issue 1, January 2020, Pages 740-757

A hyper luminous starburst at z=4.72 magnified by a lensing galaxy pair at z=1.48

Authors: L. Ciesla, M. Béthermin, E. Daddi, J. Richard, T. Diaz-Santos, M. T. Sargent, D. Elbaz, M. Boquien, T. Wang, C. Schreiber, C. Yang, J. Zabl, M. Fraser, M. Aravena, R. J. Assef, A. J. Baker, A. Beelen, A. Boselli, F. Bournaud, D. Burgarella, V. Charmandaris, P. Côté, B. Epinat, L. Ferrarese, R. Gobat , et al. (1 additional authors not shown)

Abstract: [Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm… ▽ More [Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: 14 pages, 15 figures, accepted for publication in A&A

Journal ref: A&A 635, A27 (2020)

BlueMUSE: Project Overview and Science Cases

Authors: Johan Richard, Roland Bacon, Jérémy Blaizot, Samuel Boissier, Alessandro Boselli, NicolasBouché, Jarle Brinchmann, Norberto Castro, Laure Ciesla, Paul Crowther, Emanuele Daddi, Stefan Dreizler, Pierre-Alain Duc, David Elbaz, Benoit Epinat, Chris Evans, Matteo Fossati, Michele Fumagalli, Miriam Garcia, Thibault Garel, Matthew Hayes, Angela Adamo, Artemio Herrero, Emmanuel Hugot, Andrew Humphrey , et al. (37 additional authors not shown)

Abstract: We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases… ▽ More We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena. △ Less

Submitted 9 November, 2021; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: 60 pages, 22 figures, minor updates

Revealing the dust attenuation properties on resolved scales in NGC628 with SWIFT UVOT data

Authors: Marjorie Decleir, Ilse De Looze, Médéric Boquien, Maarten Baes, Sam Verstocken, Daniela Calzetti, Laure Ciesla, Jacopo Fritz, Rob Kennicutt, Angelos Nersesian, Mathew Page

Abstract: Understanding how dust attenuation laws vary between and within galaxies is a key question if we want to reliably measure the physical properties of galaxies at both global and local scales. To shed new light on this question, we present a detailed study of the slope and bump strength of the attenuation law in the nearby spiral galaxy NGC628 at the resolved spatial scale of 325 pc. To do so, we ha… ▽ More Understanding how dust attenuation laws vary between and within galaxies is a key question if we want to reliably measure the physical properties of galaxies at both global and local scales. To shed new light on this question, we present a detailed study of the slope and bump strength of the attenuation law in the nearby spiral galaxy NGC628 at the resolved spatial scale of 325 pc. To do so, we have modelled a broad multi-wavelength dataset from the ultraviolet (UV) to the infrared (IR) with the state-of-the-art SED fitting code CIGALE, including SWIFT UVOT data for which we have developed a new optimized reduction pipeline. We find that the median dust attenuation curve of NGC628 is fairly steep, but not as steep as the SMC curve, and has a sub-MW-type UV bump. We observe intriguing variations within the galaxy, with regions of high $A_V$ exhibiting a shallower attenuation curve. We argue that the flattening of the curve is due to a dominance of absorption-over-scattering events at higher $A_V$. No trend between the bump strength and the IRAC 8.0 micron emission was found. However, this does not necessarily rule out PAHs as the main contributors to the UV bump. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 26 pages, 16 figures, 4 tables, accepted for publication in MNRAS

Exploring the halo occupation of AGN using dark-matter cosmological simulations

Authors: A. Georgakakis, J. Comparat, A. Merloni, L. Ciesla, J. Aird, A. Finoguenov

Abstract: A semi-empirical model is presented that describes the distribution of Active Galactic Nuclei (AGN) on the cosmic web. It populates dark-matter halos in N-body simulations (MultiDark) with galaxy stellar masses using empirical relations based on abundance matching techniques, and then paints accretion events on these galaxies using state-of-the-art measurements of the AGN occupation of galaxies. T… ▽ More A semi-empirical model is presented that describes the distribution of Active Galactic Nuclei (AGN) on the cosmic web. It populates dark-matter halos in N-body simulations (MultiDark) with galaxy stellar masses using empirical relations based on abundance matching techniques, and then paints accretion events on these galaxies using state-of-the-art measurements of the AGN occupation of galaxies. The explicit assumption is that the large-scale distribution of AGN is independent of the physics of black-hole fueling. The model is shown to be consistent with current measurements of the two-point correlation function of AGN samples. It is then used to make inferences on the halo occupation of the AGN population. Mock AGN are found in halos with a broad distribution of masses with a mode of $\approx 10^{12}\,h^{-1} \, M_{\odot}$ and a tail extending to cluster-size halos. The clustering properties of the model AGN depend only weakly on accretion luminosity and redshift. The fraction of satellite AGN in the model increases steeply toward more massive halos, in contrast with some recent observational results. This discrepancy, if confirmed, could point to a dependence of the halo occupation of AGN on the physics of black-hole fueling. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: 22 pages, submitted to MNRAS, light-cones and correlation functions available at https://doi.org/10.5281/zenodo.2114283

CIGALE: a python Code Investigating GALaxy Emission

Authors: M. Boquien, D. Burgarella, Y. Roehlly, V. Buat, L. Ciesla, D. Corre, A. K. Inoue, H. Salas

Abstract: Context. Measuring how the physical properties of galaxies change across cosmic times is essential to understand galaxy formation and evolution. With the advent of numerous ground-based and space-borne instruments launched over the past few decades we now have exquisite multi-wavelength observations of galaxies from the FUV to the radio domain. To tap into this mine of data and obtain new insight… ▽ More Context. Measuring how the physical properties of galaxies change across cosmic times is essential to understand galaxy formation and evolution. With the advent of numerous ground-based and space-borne instruments launched over the past few decades we now have exquisite multi-wavelength observations of galaxies from the FUV to the radio domain. To tap into this mine of data and obtain new insight into the formation and evolution of galaxies, it is essential that we are able to extract information from their SED. Aims. We present a completely new implementation of CIGALE. Written in python, its main aims are to easily and efficiently model the FUV to radio spectrum of galaxies and estimate their physical properties such as star formation rate, attenuation, dust luminosity, stellar mass, and many other physical quantities. Methods. To compute the spectral models, CIGALE builds composite stellar populations from simple stellar populations combined with highly flexible star formation histories, calculates the emission from gas ionised by massive stars, and attenuates both the stars and the ionised gas with a highly flexible attenuation curve. Based on an energy balance principle, the absorbed energy is then re-emitted by the dust in the mid- and far-infrared domains while thermal and non-thermal components are also included, extending the spectrum far into the radio range. A large grid of models is then fitted to the data and the physical properties are estimated through the analysis of the likelihood distribution. Results. CIGALE is a versatile and easy-to-use tool that makes full use of the architecture of multi-core computers, building grids of millions of models and analysing samples of thousands of galaxies, both at high speed. Beyond fitting the SEDs of galaxies and parameter estimations, it can also be used as a model-generation tool or serve as a library to build new applications. △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: 34 pages, 14 figures. Accepted for publication in A&A. The code is freely available at https://cigale.lam.fr . Comments welcome

Journal ref: A&A 622, A103 (2019)

Revealing environmental dependence of molecular gas content in a distant X-ray cluster at z=2.51

Authors: Tao Wang, David Elbaz, Emanuele Daddi, Daizhong Liu, Tadayuki Kodama, Ichi Tanaka, Corentin Schreiber, Anita Zanella, Francesco Valentino, Mark Sargent, Kotaro Kohno, Mengyuan Xiao, Maurilio Pannella, Laure Ciesla, Raphael Gobat, Yusei Koyama

Abstract: We present a census of the molecular gas properties of galaxies in the most distant known X-ray cluster, CLJ1001, at z=2.51, using deep observations of CO(1-0) with JVLA. In total 14 cluster members with $M_{*} > 10^{10.5} M_{\odot}$ are detected, including all the massive star-forming members within the virial radius, providing the largest galaxy sample in a single cluster at $z > 2$ with CO(1-0)… ▽ More We present a census of the molecular gas properties of galaxies in the most distant known X-ray cluster, CLJ1001, at z=2.51, using deep observations of CO(1-0) with JVLA. In total 14 cluster members with $M_{*} > 10^{10.5} M_{\odot}$ are detected, including all the massive star-forming members within the virial radius, providing the largest galaxy sample in a single cluster at $z > 2$ with CO(1-0) measurements. We find a large variety in the gas content of these cluster galaxies, which is correlated with their relative positions (or accretion states), with those closer to the cluster core being increasingly gas-poor. Moreover, despite their low gas content, the galaxies in the cluster center exhibit an elevated star formation efficiency (SFE=SFR/$M_{\rm gas}$) compared to field galaxies, suggesting that the suppression on the SFR is likely delayed compared to that on the gas content. Their gas depletion time is around $t_{\rm dep} \sim 400$ Myrs, comparable to the cluster dynamical time. This implies that they will likely consume all their gas within a single orbit around the cluster center, and form a passive cluster core by $z\sim2$. This result is one of the first direct pieces of evidence for the influence of environment on the gas reservoirs and SFE of $z > 2$ cluster galaxies, thereby providing new insights into the rapid formation and quenching of the most massive galaxies in the early universe. △ Less

Submitted 24 October, 2018; originally announced October 2018.

Comments: accepted for publication in ApJ Letters

The integrated properties of the CALIFA galaxies: Model-derived galaxy parameters and quenching of star formation

Authors: T. Bitsakis, S. F. Sanchez, L. Ciesla, P. Bonfini, V. Charmandaris, B. Cervantes Sodi, A. Maragkoudakis, T. Diaz-Santos, A. Zezas

Abstract: We present a study of the integrated properties of the 835 galaxies in the CALIFA survey. To derive the main physical parameters of the galaxies we have fitted their UV-to-IR spectral energy distributions (SED) with sets of theoretical models using CIGALE. We perform a comparison of the integrated galaxy parameters derived from multi-band SED fitting with those obtained from modelling the Integral… ▽ More We present a study of the integrated properties of the 835 galaxies in the CALIFA survey. To derive the main physical parameters of the galaxies we have fitted their UV-to-IR spectral energy distributions (SED) with sets of theoretical models using CIGALE. We perform a comparison of the integrated galaxy parameters derived from multi-band SED fitting with those obtained from modelling the Integral Field Unit (IFU) spectra and show the clear advantage of using the SED-derived star formation rates (SFR). A detailed analysis of galaxies in the SFR/Mstar plane as a function of their properties reveals that quenching of star formation is caused by a combination of gas deficiency and the inefficiency of the existing gas to form new stars. Exploring the plausible mechanisms that could produce this effect, we find a strong correlation with galaxy morphology and the build-up of central bulge. On the other hand, the presence of AGN and/or a stellar bar, as well as the local environment have only temporal effects on the current star formation, a result also consistent with their model-derived star formation histories. △ Less

Submitted 12 October, 2018; originally announced October 2018.

Comments: Accepted for publication in MNRAS

Comprehensive comparison of models for spectral energy distributions from 0.1 micron to 1 mm of nearby star-forming galaxies

Authors: L. K. Hunt, I. De Looze, M. Boquien, R. Nikutta, A. Rossi, S. Bianchi, D. A. Dale, G. L. Granato, R. C. Kennicutt, L. Silva, L. Ciesla, M. Relano, S. Viaene, B. Brandl, D. Calzetti, K. V. Croxall, B. T. Draine, M. Galametz, K. D. Gordon, B. A. Groves, G. Helou, R. Herrera-Camus, J. L. Hinz, J. Koda, S. Salim , et al. (4 additional authors not shown)

Abstract: We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of… ▽ More We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR. △ Less

Submitted 13 November, 2018; v1 submitted 11 September, 2018; originally announced September 2018.

Comments: 42 pages, 23 figures, 8 tables, 2 appendices: updated to v2, A&A, accepted

Journal ref: A&A 621, A51 (2019)

The bivariate luminosity and mass functions of the local HRS galaxy sample. The stellar, dust, gas mass functions

Authors: P. Andreani, A. Boselli, L. Ciesla, R. Vio, L. Cortese, V. Buat, Y. Miyamoto

Abstract: We discuss the results of the relationships between the K-band and stellar mass, far-infrared luminosities, star formation rate, dust and gas masses of nearby galaxies computing the bivariate K-band Luminosity Function (BLF) and bivariate K-band Mass Function (BMF) of the Herschel Reference Survey (HRS), a volume-limited sample with full wavelength coverage. We derive the BLFs and BMFs from the K-… ▽ More We discuss the results of the relationships between the K-band and stellar mass, far-infrared luminosities, star formation rate, dust and gas masses of nearby galaxies computing the bivariate K-band Luminosity Function (BLF) and bivariate K-band Mass Function (BMF) of the Herschel Reference Survey (HRS), a volume-limited sample with full wavelength coverage. We derive the BLFs and BMFs from the K-band and stellar mass, far-infrared luminosities, star formation rate, dust and gas masses cumulative distributions using a copula method which is outlined in detail. The use of the bivariate computed taking into account the upper limits allows us to derive on a more solid statistical ground the relationship between the observed physical quantities. The analysis shows that the behaviour of the morphological (optically selected) subsamples is quite different. A statistically meaningful result can be obtained over the whole HRS sample only from the relationship between the K-band and the stellar mass, while for the remaining physical quantities (dust and gas masses, far-IR luminosity and star formation rate), the analysis is distinct for late-type (LT) and early-type galaxies (ETG). However, the number of ETGs is small to perform a robust statistical analysis, and in most of the case results are discussed only for the LTG subsample. The Luminosity and Mass Functions (LFs, MFs) of LTGs are generally dependent on the K-band and the various dependencies are discussed in detail. We are able to derive the corresponding LFs and MFs and compare them with those computed with other samples. Our statistical analysis allows us to characterise the HRS, that, although non homogeneously selected and partially biased towards low IR luminosities, may be considered as representative of the local LT galaxy population. △ Less

Submitted 7 June, 2018; originally announced June 2018.

Comments: A&A, accepted

Journal ref: A&A 617, A33 (2018)

Identification of galaxies that experienced a recent major drop of star formation

Authors: L. Ciesla, D. Elbaz, C. Schreiber, E. Daddi, T. Wang

Abstract: [abridged] Building upon a previous study, we define a method to blindly identify galaxies that underwent, and may still be undergoing, a fast downfall of their star-formation activity, that is, a more than 80% drop in star-formation rate (SFR) occurring in less than 500 Myr. Modeling galaxies' SED with a delayed-$τ$ star formation history (SFH), with and without allowing an instantaneous SFR drop… ▽ More [abridged] Building upon a previous study, we define a method to blindly identify galaxies that underwent, and may still be undergoing, a fast downfall of their star-formation activity, that is, a more than 80% drop in star-formation rate (SFR) occurring in less than 500 Myr. Modeling galaxies' SED with a delayed-$τ$ star formation history (SFH), with and without allowing an instantaneous SFR drop within the last hundreds Myr, we isolate 102 candidates out of a subsample of 6,680 galaxies classified as star-forming from the UVJ criterion in the ZFOURGE catalogues. These galaxies are mostly located in the lower part of the SFR-M$_*$ main sequence (MS) and extend up to a factor 100 below it. They also lie close to the limit between the passive and active regions on the UVJ diagram, indicating that they are in a transition phase. We show that the selected candidates have different physical properties compared to galaxies with similar UVJ colors, namely, lower star-formation rates and different stellar masses. The morphology of the candidates show no preference for a particular type. Among the 102 candidates, only 4 show signs of an AGN activity (from X-ray luminosity or UV-IR SED fitting decomposition). This low fraction of AGNs among the candidates implies that AGN activity may not be the main driver of the recent downfall, although timescale differences and duty cycle must be taken into account. We finally attempt to recover the past position of these galaxies on the SFR-M$_*$ plane, before the downfall of their star-formation and show that some of them were in the starburst region before and are back on the MS. These candidates constitute a promising sample that need more investigation in order to understand the different mechanisms at the origin of the star formation decrease of the Universe since $z$$\sim$2. △ Less

Submitted 27 March, 2018; originally announced March 2018.

Comments: Accepted for publication in A&A