Search

We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ ($\langle z\rangle=10.0$), UV-magnitudes between M$_{\rm UV}$ = $-$21.2 and $-$19.5 (with $\langle $M$_{\rm UV}\rangle=-20.2$) and rest-frame UV slopes ($\langle \beta\rangle=-2.4$). These galaxies are, on average, more luminous than most $z\ge9$ candidates discovered by JWST so far in the literature, while exhibiting similar blue colors in their rest-frame UV. The rest-frame UV slopes derived from SED-fitting are blue ($\beta\sim$[$-$2.0, $-$2.7]) without reaching extremely blue values as reported in other recent studies at these redshifts. The blue color is consistent with models that suggest the underlying stellar population is not yet fully enriched in metals like similarly luminous galaxies in the lower redshift Universe. The derived stellar masses with $\langle \log_{\rm 10} ($M$_\star/$M$_\odot)\rangle\approx8-9$ are not in tension with the standard $\Lambda$CDM model and our measurement of the volume density of such UV luminous galaxies aligns well with previously measured values presented in the literature at $z\sim9-10$. Our sample of galaxies, although compact, are significantly resolved., Comment: Submitted to ApJ

We present the results of a systematic study of the rest-frame optical morphology of quiescent galaxies at $z \geq 3$ using the Near-Infrared Camera (NIRCam) onboard $JWST$. Based on a sample selected by $UVJ$ color or $NUVUVJ$ color, we focus on 26 quiescent galaxies with $9.8<\log{(M_\star/M_\odot)}<11.4$ at $2.810.3$ implies that our size - stellar mass relations are below those at lower redshifts, with the amplitude of $\sim0.6\, {\rm kpc}$ at $M_\star = 5\times 10^{10}\, M_\odot$. This value agrees with the extrapolation from the size evolution of quiescent galaxies at $z<3$ in the literature, implying that the size of quiescent galaxies increases monotonically from $z\sim3-5$. Our sample is mainly composed of galaxies with bulge-like structures according to their median S\'ersic index and axis ratio of $n\sim3-4$ and $q\sim0.6-0.8$, respectively. On the other hand, there is a trend of increasing fraction of galaxies with low S\'ersic index, suggesting $3

A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ($z\gtrsim1$) Universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution in X-ray selected AGN up to $z\sim2$. Here, we focus on a sample of three AGNs with their host galaxies exhibiting prominent dust lanes, potentially due to their edge-on alignment. These represent 27% (3 out of 11 with early NIRCam data) of the heavily obscured ($N_H>10^{23}$ cm$^{-2}$) AGN population. With limited signs of a central AGN in the optical and near-infrared, the NIRCam images are used to produce reddening maps $E(B-V)$ of the host galaxies. We compare the mean central value of $E(B-V)$ to the X-ray obscuring column density along the line-of-sight to the AGN ($N_H\sim10^{23-23.5}$ cm$^{-2}$). We find that the extinction due to the host galaxy is present ($0.6\lesssim E(B-V) \lesssim 0.9$; $1.9 \lesssim A_V \lesssim 2.8$) and significantly contributes to the X-ray obscuration at a level of $N_H\sim10^{22.5}$ cm$^{-2}$ assuming an SMC gas-to-dust ratio which amounts to $\lesssim$30% of the total obscuring column density. These early results, including three additional cases from CEERS, demonstrate the ability to resolve such dust structures with JWST and separate the different circumnuclear and galaxy-scale obscuring structures., Comment: 12 pages, 8 figures, Accepted for publication in ApJL

Understanding the star formation rate (SFR) variability and how it depends on physical properties of galaxies is important for developing and testing the theory of galaxy formation. We investigate how statistical measurements of the extragalactic background light (EBL) can shed light on this topic and complement traditional methods based on observations of individual galaxies. Using semi-empirical models of galaxy evolution and SFR indicators sensitive to different star formation timescales (e.g., H$\alpha$ and UV continuum luminosities), we show that the SFR variability, quantified by the joint probability distribution of the SFR indicators (i.e., the bivariate conditional luminosity function), can be characterized as a function of galaxy mass and redshift through the cross-correlation between deep, near-infrared maps of the EBL and galaxy distributions. As an example, we consider combining upcoming SPHEREx maps of the EBL with galaxy samples from Rubin/LSST. We demonstrate that their cross-correlation over a sky fraction of $f_\mathrm{sky}\sim0.5$ can constrain the joint SFR indicator distribution at high significance up to $z\sim2.5$ for mass-complete samples of galaxies down to $M_{*}\sim10^9\,M_{\odot}$. These constraints not only allow models of different SFR variability to be distinguished, but also provide unique opportunities to investigate physical mechanisms that require large number statistics such as environmental effects. The cross-correlations investigated illustrate the power of combining cosmological surveys to extract information inaccessible from each data set alone, while the large galaxy populations probed capture ensemble-averaged properties beyond the reach of targeted observations towards individual galaxies., Comment: 12 pages, 7 figures, MNRAS accepted

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

A growing number of far-infrared bright sources completely invisible in deep extragalactic optical surveys hint at an elusive population of z>4 dusty, star-forming galaxies. Cycle 1 JWST surveys are now detecting their rest-frame optical light, which provides key insight into their stellar properties and statistical constraints on the population as a whole. This work presents the JWST/NIRCam counterpart from the COSMOS-Web survey to a far-infrared SCUBA-2 and ALMA source, AzTECC71, which was previously undetected at wavelengths shorter than 850 microns. AzTECC71, amongst the reddest galaxies in COSMOS-Web with F277W - F444W~0.9, is undetected in NIRCam/F150W and F115W and fainter in F444W than other sub-millimeter galaxies identified in COSMOS-Web by 2-4 magnitudes. This is consistent with the system having both a lower stellar mass and higher redshift than the median dusty, star-forming galaxy. With deep ground- and space-based upper limits combined with detections in F277W, F444W and the far-IR including ALMA Band 6, we find a high probability (99%) that AzTECC71 is at z>4 with z_phot=5.7(+0.8,-0.7). This galaxy is massive (logM*/Msun~10.7) and IR-luminous (logLIR/Lsun~12.7), comparable to other optically-undetected but far-IR bright dusty, star-forming galaxies at z>4. This population of luminous, infrared galaxies at z>4 is largely unconstrained but comprises an important bridge between the most extreme dust-obscured galaxies and more typical high-redshift star-forming galaxies. If further far-IR-selected galaxies that drop out of the F150W filter in COSMOS-Web have redshifts z>4 like AzTECC71, then the volume density of such sources may be ~3-10x greater than previously estimated., Comment: 13 pages, 7 figures, submitted to ApJ

Utilizing the James Webb Space Telescope Early Release NIRCam Observations, we perform a weak-lensing analysis of the massive galaxy cluster SMACS J0723.3-7327 ($z=0.39$). We investigate the spatial variation of the PSF from the stars in the mosaic image. Our measurements show that the PSF for both modules has very small spatial and temporal variation with average complex ellipticity components of $e_1=0.007\pm0.001$ and $e_2=0.029\pm0.001$ in the observed north-up reference frame. We create PSF models through a principal component analysis of the stars and show that they properly account for the ellipticity of the PSF with residual shapes of $e_1=(0.3\pm3.5)\times10^{-4}$ and $e_2=(1.8\pm4.0)\times10^{-4}$. We select background galaxies by their photometric redshift and measure galaxy shapes by model fitting. Our weak-lensing source catalog achieves 215 galaxies arcmin$^{-2}$. We map the projected mass density of SMACSJ0723 and detect the cluster with a peak significance of $12.2\sigma$. The mass distribution is found to elongate in the east-west direction with an extension to the northeast edge of the field of view where a candidate substructure is found in the Chandra X-ray imaging. We fit the tangential shear with a Navarro-Frenk-White model and estimate the mass of the cluster to be $M_{500}=7.9\pm1.1\times10^{14}$ M$_{\odot}$ ($M_{200}=11.4\pm1.5\times10^{14}$ M$_\odot$ ), which agrees with existing mass estimates. Combining the multiwavelength evidence from literature with our weak-lensing analysis, we hypothesize that SMACSJ0723 is observed near first pericenter passage and we identify candidate radio relics., Comment: 14 pages, 11 figures, submitted to ApJ

We present the results of a systematic search for candidate quiescent galaxies in the distant Universe in eleven $JWST$ fields with publicly available observations collected during the first three months of operations and covering an effective sky area of $\sim145$ arcmin$^2$. We homogeneously reduce the new $JWST$ data and combine them with existing observations from the $Hubble\,Space\,Telescope$. We select a robust sample of $\sim80$ candidate quiescent and quenching galaxies at $3 < z < 5$ using two methods: (1) based on their rest-frame $UVJ$ colors, and (2) a novel quantitative approach based on Gaussian Mixture Modeling of the $NUV-U$, $U-V$, and $V-J$ rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ($M_\star\geq 10^{10.6} M_\odot$) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of $2-3$, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at $z>3$, as it could be expected for highly clustered massive systems. Importantly, $JWST$ enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available., Comment: 15 pages, 6 Figures + Appendix. Accepted for publication in ApJ on Feb, 9. Data release: - Reduced HST+JWST mosaics + photometric catalogs and Eazy-py modeling: https://erda.ku.dk/archives/7166d013c1ca1371aac3c57b9e73190d/published-archive.html - Supplementary material and tables: https://zenodo.org/record/7614908#.Y-4ZruzMLmE - MAST: https://doi.org/10.17909/g3nt-a370 - See also Gould et al. 2023

The quasar SDSS 0956+5128 exhibits three distinct velocity components with large offsets in emission: the systemic velocity of [O II], [O III] and [Ne III] narrow lines have redshift $z=0.7142$; broad Mg II line is shifted by $-1200$ km s$^{-1}$ with respect to the narrow lines; broad H$\alpha$, H$\beta$ lines are at $-4100$ km s$^{-1}$. We present new Hubble Space Telescope spectra of Ly$\alpha$ and C IV emission lines and high-resolution images of the quasar. The offsets of these lines are consistent with the velocity component of the Balmer emission, and the photometry in optical and near-infrared wavelengths does not show any signs of recent mergers in the host galaxy or irregularities in the location of the quasar. The data do not confirm predictions of the previous most-likely hypotheses involving a special orientation and morphology of the quasar disk, such as in the recoiling black hole scenario, neither it is consistent with accretion disk winds. Instead, based on the cumulative evidence, we propose a new scenario, in which the broad line region is in the state of outflow caused by a strong shock wave, with a supernova as a possible event for producing the shock ejecta., Comment: 15 pages, 4 figures, accepted for publication in ApJ

How galaxies form, assemble, and cease their star-formation is a central question within the modern landscape of galaxy evolution studies. These processes are indelibly imprinted on the galaxy stellar mass function (SMF). We present constraints on the shape and evolution of the SMF, the quiescent galaxy fraction, and the cosmic stellar mass density across 90% of the history of the Universe from $z=7.5\rightarrow0.2$ via the COSMOS survey. Now with deeper and more homogeneous near-infrared coverage exploited by the COSMOS2020 catalog, we leverage the large 1.27 deg$^{2}$ effective area to improve sample statistics and understand cosmic variance particularly for rare, massive galaxies and push to higher redshifts with greater confidence and mass completeness than previous studies. We divide the total stellar mass function into star-forming and quiescent sub-samples through $NUVrJ$ color-color selection. Measurements are then fitted with Schechter functions to infer the intrinsic SMF, the evolution of its key parameters, and the cosmic stellar mass density out to $z=7.5$. We find a smooth, monotonic evolution in the galaxy SMF since $z=7.5$, in agreement with previous studies. The number density of star-forming systems seems to have undergone remarkably consistent growth spanning four decades in stellar mass from $z=7.5\rightarrow2$ whereupon high-mass systems become predominantly quiescent (i.e. downsizing). An excess of massive systems at $z\sim2.5-5.5$ with strikingly red colors, some newly identified, increase the observed number densities to the point where the SMF cannot be reconciled with a Schechter function. Systematics including cosmic variance and/or AGN contamination are unlikely to fully explain this excess, and so we speculate that there may be contributions from dust-obscured objects similar to those found in FIR surveys. (abridged), Comment: 39 pages, 24 figures, accepted for publication in A&A. Data files containing key measurements are available for download: https://doi.org/10.5281/zenodo.7808832

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

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$\sigma$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$\sigma$ point source depths of $\sim$25.3-26.0 magnitudes. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($64$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery., Comment: 46 pages, 16 figures, ApJ accepted

ALMA observations of the long wavelength dust continuum are used to estimate the gas masses in a sample of 708 star-forming (SF) galaxies at z = 0.3 to 4.5. We determine the dependence of gas masses and star formation efficiencies (SFE=SFR per unit gass mass). We find that 70 percent of the increase in SFRs of the MS is due to the increased gas masses at earlier epochs while 30 percent is due to increased efficiency of SF. For galaxies above the MS this is reversed with 70 percent of the increased SFR relative to the MS being due to elevated SFEs. Thus, the major evolution of star formation activity at early epochs is driven by increased gas masses, while the starburst activity taking galaxies above the MS is due to enhanced triggering of star formation (likely due to galactic merging). The interstellar gas peaks at z = 2 and dominates the stellar mass down to z = 1.2. Accretion rates needed to maintain continuity of the MS evolution exceed 100 Msun per yr at z > 2. The galactic gas contents are likely the driving determinant for both the rise in SF and AGN activity from z = 5 to their peak at z = 2 and subsequent fall to lower z. We suggest that for self-gravitating clouds with supersonic turbulence, cloud collisions and the filamentary structure of the clouds regulate the star formation activity., Comment: AAS Russel Lecture June 2022 (N. Scoville), 18 pages. arXiv admin note: text overlap with arXiv:1702.04729

Recent ALMA large surveys unveiled the presence of significant dust continuum emission in star-forming galaxies at $z>4$. Unfortunately, such large programs -- i.e. ALPINE ($z\sim 5$) and REBELS ($z \sim 7$) -- only provide us with a single Far-Infrared (FIR) continuum data point for their individual targets. Therefore, high-$z$ galaxies FIR spectral energy densities (SEDs) remain mostly unconstrained, hinging on an assumption for their dust temperature ($T_{\rm d}$) in the SED fitting procedure. This introduces uncertainties in the inferred dust masses ($M_{\rm d }$), infrared luminosities ($L_{\rm IR}$), and obscured Star Formation Rate (SFR) fraction at $z > 4$. In this work we use a method that allows us to constrain $T_{\rm d}$ with a single band measurement by combining the $158\ \mathrm{\mu m}$ continuum information with the overlying [CII] emission line. We analyse the $21$ [CII] and FIR continuum detected $z\sim 5$ galaxies in ALPINE, finding a range of $T_{\rm d}=25-60\ \mathrm{K}$ and $M_{\rm d} = 0.6-25.1\ \times 10^{7}\ \mathrm{M_{\odot}}$. Given the measured stellar masses of ALPINE galaxies, the inferred dust yields are around $M_{\rm d}/M_{\star} = (0.2-8) \times 10^{-3}$, consistent with theoretical dust-production constraints. We find that $8$ out of $21$ ALPINE galaxies have $L_{\rm IR} \geq 10^{12}\ \mathrm{L_{\odot}}$, comparable to UltraLuminous IR Galaxies (ULIRGs). Relying on ultraviolet-to-optical SED fitting, the SFR was underestimated by up to $2$ orders of magnitude in $4$ of these $8$ ULIRGs-like galaxies. We conclude that these $4$ peculiar sources should be characterised by a two-phase interstellar medium structure with "spatially-segregated" FIR and ultraviolet emitting regions., Comment: Accepted for publication in MNRAS. 13 pages, 4 figures

We use the unprecedented resolution and depth of the JWST NIRCam Early Release Observations at 1-5 $\mu$m to study the stellar mass, age, and metallicity of compact star clusters in the neighborhood of the host galaxies in the SMACS J0723.3-7327 galaxy cluster field at z = 0.39. The measured colors of these star clusters show a similar distribution as quiescent galaxies at the same redshift, but are >3 magnitudes fainter than the current depths of wide-field galaxy survey. The star clusters are unresolved in the NIRCam/F150W data suggesting sizes smaller than 50pc. This is significantly smaller than star forming clumps or dwarf galaxies in local galaxies. From fitting their photometry with simple stellar population (SSP) models, we find stellar metallicities consistent with 0.2-0.3 $Z_{\odot}$ and ages of $1.5^{+0.5}_{-0.5}$ Gyrs. We rule out metallicities <0.2 $Z_{\odot}$ and solar/super-solar at 4$\sigma$ significance. Assuming mass-to-light ratios obtained from the best-fit SSPs, we estimate stellar masses of $2.4^{+3.0}_{-1.5}\times10^6$ M$_{\odot}$. These are between average masses of local globular clusters and dwarf galaxies. Our analysis suggests middle-aged globulars with relatively recent formation times at z=0.5-0.7, which could have been subsequently stripped away from their host galaxies due to interactions in the cluster environment, or formed in cold flows onto the cluster core. However, we cannot rule out these objects being compact cores of stripped dwarf galaxies., Comment: 20 pages, 8 figures, 2 tables. Accepted for publication in ApJ

We present the radio properties of 66 spectroscopically-confirmed normal star-forming galaxies (SFGs) at $4.4

A multiwavelength study of galaxies is important to understand their formation and evolution. Only in the recent past, thanks to the Atacama Large (Sub) Millimeter Array (ALMA), were we able to study the far-infrared (IR) properties of galaxies at high redshifts. In this article, we summarize recent research highlights and their significance to our understanding of early galaxy evolution from the ALPINE survey, a large program with ALMA to observe the dust continuum and 158um C+ emission of normal star-forming galaxies at z = 4-6. Combined with ancillary data at UV through near-IR wavelengths, ALPINE provides the currently largest multiwavelength sample of post-reionization galaxies and has advanced our understanding of (i) the demographics of C+ emission; (ii) the relation of star formation and C+ emission; (iii) the gas content; (iv) outflows and enrichment of the intergalactic medium; and (v) the kinematics, emergence of disks, and merger rates in galaxies at z > 4. ALPINE builds the basis for more detailed measurements with the next generation of telescopes, and places itself as an important post-reionization baseline sample to allow a continuous study of galaxies over 13 billion years of cosmic time., Comment: 29 pages, 6 figures; Invited article for the special issue "Recent Advances in Infrared Galaxies and AGN", edited by Anna Sajina and Asantha R. Cooray, in Universe

We present the average [CII] $158\,\rm{\mu m}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{\mu m}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4

We use the COSMOS2020 catalogue to measure the stellar-to-halo mass relation (SHMR) divided by central and satellite galaxies from $z=0.2$ to $z = 5.5$. Starting from accurate photometric redshifts we measure the near-infrared selected two-point angular correlation and stellar mass functions in ten redshift bins and fit them with an HOD-based model. At each redshift, we measure the ratio of stellar mass to halo mass, $M_*/M_h$, which shows the characteristic strong dependence of halo mass with a peak at $M_h^{\rm peak} \sim 10^{12}\, M_{\odot}$. Our results are in accordance with the scenario in which the peak of star-formation efficiency moves towards more massive halos at higher redshifts. We also measure the fraction of satellites as a function of stellar mass and redshift. For all stellar mass thresholds the satellite fraction decreases at higher redshifts. At a given redshift there is a higher fraction of low-mass satellites. The satellite contribution to the total stellar mass budget in halos becomes more important than centrals at halo masses of about $M_h > 10^{13} \, M_{\odot}$ and always stays below by peak, indicating that quenching mechanisms are present in massive halos that keep the star-formation efficiency low. Finally, we compare our results with three hydrodynamical simulations Horizon-AGN, Illustris-TNG-100 and EAGLE. We find that the most significant discrepancy is at the high mass end, where the simulations generally show that satellites have a higher contribution to the total stellar mass budget than the observations. This, together with the finding that the fraction of satellites is higher in the simulations, indicates that the feedback mechanisms acting in group-and cluster-scale halos appear to be less efficient in quenching the mass assembly of satellites, and/or that quenching occurs much later in the simulations.

We analyze stellar proper motions in the COSMOS field to assess the presence of bulk motions. At bright magnitudes (G-band 18.5--20.76 AB), we use the proper motions of 1,010 stars in the Gaia DR2 catalog. At the faint end, we computed proper motions of 11,519 point-like objects at i-band magnitudes 19--25 AB using Hubble ACS and Subaru HSC, which span two epochs about 11 years apart. In order to measure these proper motions with unprecedented accuracy at faint magnitudes, we developed a foundational set of astrometric tools which will be required for Joint Survey Processing (JSP) of data from the next generation of optical/infrared surveys. The astrometric grids of Hubble ACS and Subaru HSC mosaics were corrected at the catalog level, using proper motion-propagated and parallax-corrected Gaia DR2 sources. These astrometric corrections were verified using compact extragalactic sources. Upon comparison of our measured proper motions with Gaia DR2, we estimate the uncertainties in our measurements to be ~2--3 mas/yr per axis, down to 25.5 AB mag. We corrected proper motions for the mean motion of the Sun, and we find that late-type main-sequence stars predominantly in the thin disk in the COSMOS field have space velocities mainly towards the Galactic center. We detect candidate high-velocity (> 220 km/s) stars, 6 of them at ~0.4-6 kpc from the Gaia sample, and 5 of them at ~20 kpc from the faint star HSC and ACS sample. The sources from the faint star sample may be candidate halo members of the Sangarius stream., Comment: 33 pages; 22 figures; 6 tables; 3 appendices; 7 authors

The luminosity functions at z < 4 - 5 suggest that most galaxies have a relatively low stellar mass (logM_star = 10) and a low dust attenuation (A_FUV = 1.0). The physical properties of these objects are quite homogeneous. We used an approach where we combined their rest-frame far-infrared and submillimeter emissions and utilized the universe and the redshift as a spectrograph to increase the amount of information in a collective way. From a subsample of 27 ALMA-detected galaxies at z > 4.5, we built an infrared spectral energy distribution composite template. It was used to fit, with CIGALE, the 105 galaxies (detections and upper limits) in the sample from the FUV to the FIR. The derived physical parameters provide information to decipher the nature of the dust cycle and of the stellar populations in these galaxies. The derived IR composite template is consistent with the galaxies in the studied sample. A delayed star formation history with tau_main = 500 Myrs is slightly favored by the statistical analysis as compared to a delayed with a final burst or a continuous star formation history. The position of the sample in the star formation rate (SFR)- M_star diagram is consistent with previous papers. The redshift evolution of the log M_star versus A_FUV relation is in agreement with evolution in the redshift of this relation. This evolution is necessary to explain the cosmic evolution of the average dust attenuation of galaxies. Evolution is also observed in the L_dust/ L_FUV (IRX) versus UV slope beta_FUV diagram: younger galaxies have bluer beta_FUV. We modeled the shift of galaxies in the IRX versus the beta_FUV diagram with the mass-weighted age as a free parameter, and we provide an equation to make predictions., Comment: 40 pages, 13 figures

There is now ample evidence that dust is already present in abundance at high z. However, given the faintness of distant galaxies in the optical and the NIR, datasets are still limited and how the dust affects the emerging radiation of galaxies at very high redshift is not yet fully understood. Using the ALPINE survey, our objective is to quantify the dust attenuation properties in galaxies at z=4.4-5.5, and in particular the shape of their attenuation curve. Using the CIGALE code, we model the stellar populations and their interaction with the dust in order to measure some of the physical properties of a subsample of 23 main-sequence ALPINE galaxies. We find that the attenuation curves span a broad range of properties, from curves that are much steeper than the SMC extinction curve, to shallower than the starburst attenuation curve. The shape of the attenuation curves strongly depends on the V-band attenuation. Galaxies with the lowest attenuation also present the steepest curves. The steepness of such curves is probably the consequence of the combination of the intrinsic physical properties of the dust, the relative distribution of stars and dust in the interstellar medium, and the differential reddening. The broad range of attenuation curves found at z~5 shows that no single attenuation curve is appropriate for main sequence galaxies and that assuming a fixed curve can lead to large errors, for instance in the interpretation and use of the IRX-beta diagram, if SED modeling is not feasible. Great caution should be exercised when correcting high redshift galaxies for the presence of dust using the UV slope beta as it can affect the estimation of both SFR and stellar mass even at low V-band attenuation due to the steepness of the attenuation curve. However, when SED modeling can be used, the impact of the choice of the attenuation curve on the SFR and the stellar mass is limited., Comment: 20 pages, 14 figures. Accepted for publication in A&A

We investigate the strong-lensing cluster Abell 370 (A370) using a wide Integral Field Unit (IFU) spectroscopic mosaic from the Multi-Unit Spectroscopic Explorer (MUSE). IFU spectroscopy provides significant insight into the structure and mass content of galaxy clusters, yet IFU-based cluster studies focus almost exclusively on the central Einstein-radius region. Covering over 14 arcmin$^2$, the new MUSE mosaic extends significantly beyond the A370 Einstein radius, providing, for the first time, a detailed look at the cluster outskirts. Combining these data with wide-field, multi-band Hubble Space Telescope (HST) imaging from the BUFFALO project, we analyse the distribution of objects within the cluster and along the line of sight. Identifying 416 cluster galaxies, we use kinematics to trace the radial mass profile of the halo, providing a mass estimate independent from the lens model. We also measure radially-averaged properties of the cluster members, tracking their evolution as a function of infall. Thanks to the high spatial resolution of our data, we identify six cluster members acting as galaxy-galaxy lenses, which constrain localized mass distributions beyond the Einstein radius. Finally, taking advantage of MUSE's 3D capabilities, we detect and analyse multiple spatially extended overdensities outside of the cluster that influence lensing-derived halo mass estimates. We stress that much of this work is only possible thanks to the robust, extended IFU coverage, highlighting its importance even in less optically dense cluster regions. Overall, this work showcases the power of combining HST+MUSE, and serves as the initial step towards a larger and wider program targeting several clusters., Comment: 20 pages, 12 figures, 3 tables (including 1 in an appendix). Accepted in MNRAS. Data Release products available at https://astro.dur.ac.uk/~hbpn39/pilot-wings.html

The [CII] 158$~\mu$m emission line represents so far one of the most profitable tools for the investigation of the high-redshift galaxies in the early Universe. Being one of the brightest cooling lines in the rest-frame far-infrared regime of star-forming galaxies, it has been successfully exploited as a tracer of star-formation rate (SFR) in local sources. The picture is more complex at higher redshifts, where its usability in this context is still under investigation. Recent results from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey suggest that there is no (or weak) evolution of the L$\mathrm{_{[CII]}}$-SFR relation up to $z\sim6$ but their reliability is hampered by the presence of a large population of [CII] non-detected galaxies. In this work, we characterize the population of [CII] non-detections in ALPINE. By stacking their ALMA spectra, we obtain a signal detected at $\sim5.1\sigma$, resulting in a [CII] luminosity of $\mathrm{log(L_\mathrm{[CII]}}/\mathrm{L_{\odot}})$ $\sim7.8$. When combining this value with those from the [CII] detections, we find a $\mathrm{L_{[CII]}}$-SFR relation with a slope $b=1.14\pm0.11$, in agreement within the uncertainties both with the linear relation found in the local Universe, and with the previous findings from ALPINE at $z\sim5$. This suggests that the [CII] line can be considered a good tracer of star formation up to the distant Universe. Finally, we show that the galaxies of our sample that most deviate from the observed L$_\mathrm{[CII]}$-SFR relation could suffer from a less precise redshift estimation, perhaps artificially reducing their [CII] luminosity. In this respect, we claim that there is no evidence in favour of a deficit of [CII] content in high-z galaxies, in contrast with earlier studies., Comment: Accepted for publication in A&A. 10 pages, 5 figures

We present $10$ main-sequence ALPINE galaxies (log($M/M_{\odot}$) = 9.2-11.1 and ${\rm SFR}=23-190\,{\rm M_{\odot}\,yr^{-1}}$) at $z\sim4.5$ with optical [OII] measurements from Keck/MOSFIRE spectroscopy and Subaru/MOIRCS narrow-band imaging. This is the largest such multi-wavelength sample at these redshifts, combining various measurements in the ultra-violet, optical, and far-infrared including [CII]$_{158{\rm \mu m}}$ line emission and dust continuum from ALMA and H$\alpha$ emission from Spitzer photometry. For the first time, this unique sample allows us to analyze the relation between [OII] and total star-formation rate (SFR) and the interstellar medium (ISM) properties via [OII]/[CII] and [OII]/\halpha luminosity ratios at $z\sim4.5$. The [OII]$-$SFR relation at $z\sim4.5$ cannot be described using standard local descriptions, but is consistent with a metal-dependent relation assuming metallicities around $50\%$ solar. To explain the measured dust-corrected luminosity ratios of $L[OII]/L[CII] \sim 0.98^{+0.21}_{-0.22}$ and $L[OII]/LHa \sim -0.22^{+0.13}_{-0.15}$ for our sample, ionization parameters $\log(U)< -2$ and electron densities $\log(\rm n_e / {\rm [cm^{-3}]}) \sim 2.5-3$ are required. The former is consistent with galaxies at $z\sim2-3$, however lower than at $z>6$. The latter may be slightly higher than expected given the galaxies' specific SFR. The analysis of this pilot sample suggests that typical log($ M/M_{\odot})$ > 9 galaxies at $z\sim4.5$ to have broadly similar ISM properties as their descendants at $z\sim2$ and suggest a strong evolution of ISM properties since the Epoch of Reionization at $z>6$., Comment: 14 pages, 8 figures, Accepted for publication to MNRAS

We report CO(5$\rightarrow$4) and CO(6$\rightarrow$5) line observations in the dusty starbursting galaxy CRLE ($z = 5.667$) and the main-sequence (MS) galaxy HZ10 ($z = 5.654$) with the Northern Extended Millimeter Array (NOEMA). CRLE is the most luminous $z>5$ starburst in the COSMOS field and HZ10 is the most gas-rich "normal" galaxy currently known at $z>5$. We find line luminosities for CO(5$\rightarrow$4) and CO(6$\rightarrow$5) of (4.9 $\pm$ 0.5) and (3.8 $\pm$ 0.4) $\times$ 10$^{10}$ K km s$^{-1}$ pc$^{2}$ for CRLE and upper limits of $< 0.76$ and $< 0.60$ $\times$ 10$^{10}$ K km s$^{-1}$ pc$^{2}$ for HZ10, respectively. The CO excitation of CRLE appears comparable to other $z>5$ dusty star-forming galaxies (DSFGs). For HZ10, these line luminosity limits provide the first significant constraints of this kind for a MS galaxy at $z > 5$. We find the upper limit of $L'_{5\rightarrow4}/L'_{2\rightarrow1}$ in HZ10 could be similar to the average value for MS galaxies around $z\approx 1.5$, suggesting that MS galaxies with comparable gas excitation may already have existed one billion years after the Big Bang. For CRLE we determine the most likely values for the H$_2$ density, kinetic temperature and dust temperature based on excitation modeling of the CO line ladder. We also derive a total gas mass of $(7.1 \pm 1.3) \times 10^{10} M_\odot$. Our findings provide some of the currently most detailed constraints on the gas excitation that sets the conditions for star formation in a galaxy protocluster environment at $z > 5$., Comment: 14 pages, 6 figures, accepted to ApJ

ISCEA (Infrared Satellite for Cosmic Evolution Astrophysics) is a small astrophysics mission whose Science Goal is to discover how galaxies evolved in the cosmic web of dark matter at cosmic noon. Its Science Objective is to determine the history of star formation and its quenching in galaxies as a function of local density and stellar mass when the Universe was 3-5 Gyrs old (1.2 10 Mpc in each of 50 protocluster (cluster and cosmic web) fields at 1.2 < z < 2.1. ISCEA measures the star formation quenching factor in those fields, and galaxy kinematics with a precision < 50 km/s to deduce the 3D spatial distribution in each field. ISCEA will transform our understanding of galaxy evolution at cosmic noon. ISCEA is a small satellite observatory with a 30cm equivalent diameter aperture telescope with a FoV of 0.32 deg^2, and a multi-object spectrograph with a digital micro-mirror device (DMD) as its programmable slit mask. ISCEA will obtain spectra of 1000 galaxies simultaneously at an effective resolving power of R=1000, with 2.8"x2.8" slits, over the NIR wavelength range of 1.1 to 2.0 microns, a regime not accessible from the ground without large gaps in coverage. ISCEA will achieve a pointing accuracy of <= 2" FWHM over 200s. ISCEA will be launched into a Low Earth Orbit, with a prime mission of 2.5 years. ISCEA's space-qualification of DMDs opens a new window for spectroscopy from space, enabling revolutionary advances in astrophysics., Comment: 40 pages, 31 figures

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s., Comment: 15 pages with 11 figures, approved by Euclid Consortium Publication Board and submitted to Astronomy and Astrophysics. Data products will become available via the IRSA website once the paper is accepted. This paper is a companion to "COSMOS2020: A panchromatic view of the Universe to z~10 from two complementary catalogs" by John Weaver et al., which is being posted in parallel

The Cosmic Evolution Survey (COSMOS) has become a cornerstone of extragalactic astronomy. Since the last public catalog in 2015, a wealth of new imaging and spectroscopic data has been collected in the COSMOS field. This paper describes the collection, processing, and analysis of this new imaging data to produce a new reference photometric redshift catalog. Source detection and multi-wavelength photometry is performed for 1.7 million sources across the $2\,\mathrm{deg}^{2}$ of the COSMOS field, $\sim$966,000 of which are measured with all available broad-band data using both traditional aperture photometric methods and a new profile-fitting photometric extraction tool, The Farmer, which we have developed. A detailed comparison of the two resulting photometric catalogs is presented. Photometric redshifts are computed for all sources in each catalog utilizing two independent photometric redshift codes. Finally, a comparison is made between the performance of the photometric methodologies and of the redshift codes to demonstrate an exceptional degree of self-consistency in the resulting photometric redshifts. The $i<21$ sources have sub-percent photometric redshift accuracy and even the faintest sources at $25

Galaxy mergers are thought to be one of the main mechanisms of the mass assembly of galaxies. Recently, many works have suggested a possible increase in the fraction of major mergers in the early Universe, reviving the debate on which processes (e.g., cold accretion, star formation, mergers) most contribute to the mass build-up of galaxies through cosmic time. To estimate the importance of major mergers in this context, we make use of the new data collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE), which observed the [CII] 158 $\mu$m emission line from a sample of 75 main-sequence star-forming galaxies at 4.4 < z < 5.9. We used, for the first time, the morpho-kinematic information provided by the [CII] emission to obtain the fraction of major mergers ($f_{MM}$) at z~5. By adopting different prescriptions for the merger timescales ($T_{MM}$), we converted this fraction into the merger rate per galaxy ($R_{MM}$) and per volume ($\Gamma_{MM}$). We then combined our results with those at lower redshifts from the literature, computing the cosmic evolution of the merger fraction. This is described by a rapid increase from z~0 to higher redshifts, a peak at z~3, and a slow decrease towards earlier epochs. Depending on the timescale prescription used, this fraction translates into a merger rate ranging between ~0.1 and ~4.0 Gyr$^{-1}$ at z~5. Finally, we compare the specific star formation and star-formation rate density with the analogous quantities from major mergers. Our new ALPINE data reveal the presence of a significant merging activity in the early Universe. However, whether this population of mergers can provide a relevant contribution to the galaxy mass assembly at these redshifts and through the cosmic epochs is strongly dependent on the assumption of the merger timescale., Comment: Accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics

The dust content of normal galaxies and the dust mass density (DMD) at high-z (z>4) are unconstrained given the source confusion and the sensitivity limitations of previous observations. The ALMA Large Program to INvestigate [CII] at Early Times (ALPINE), which targeted 118 UV-selected star-forming galaxies at 4.4

The mean size ( effective radius $R_e$) of Massive Galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$) is observed to increase steadily with cosmic time. It is still unclear whether this trend originates from the size growth of individual galaxies (via, e.g., mergers and/or AGN feedback) or from the inclusion of larger galaxies entering the selection at later epochs (progenitor bias). We here build a data-driven, flexible theoretical framework to probe the structural evolution of MGs. We assign galaxies to dark matter haloes via stellar mass-halo mass (SMHM) relations with varying high-mass slopes and scatters $\sigma_{\rm SMHM}$ in stellar mass at fixed halo mass, and assign sizes to galaxies using an empirically-motivated, constant and linear relationship between $R_e$ and the host dark matter halo radius $R_h$. We find that: 1) the fast mean size growth of MGs is well reproduced independently of the shape of the input SMHM relation; 2) the numbers of compact MGs grow steadily until $z\gtrsim2$ and fall off at lower redshifts, suggesting a lesser role of progenitor bias at later epochs; 3) a time-independent scatter $\sigma_{\rm SMHM}$ is consistent with a scenario in which compact starforming MGs transition into quiescent MGs in a few $10^8$yr with a negligible structural evolution during the compact phase, while a scatter increasing at high redshift implies significant size growth during the starforming phase. A robust measurement of the size function of MGs at high redshift can set strong constraints on the scatter of the SMHM relation and, by extension, on models of galaxy evolution., Comment: Accepted for publication in MNRAS. Comments still welcome

We use the Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) strong lensing image catalog of the merging galaxy cluster Abell 370 to obtain a mass model using the free-form lens inversion algorithm GRALE. The improvement of the strong lensing data quality results in a lens plane rms of only 0.45 arcsec, about a factor of two lower than that of our existing HFF v4 reconstruction. We attribute the improvement to spectroscopic data and use of the full reprocessed HST mosaics. In our reconstructed mass model, we found indications of three distinct mass features in Abell 370: (i) a $\sim 35$ kpc offset between the northern BCG and the nearest mass peak, (ii) a $\sim 100$ kpc mass concentration of roughly critical density $\sim 250$ kpc east of the main cluster, and (iii) a probable filament-like structure passing N-S through the cluster. While (i) is present in some form in most publicly available reconstructions spanning the range of modeling techniques: parametric, hybrid, and free-form, (ii) and (iii) are recovered by only about half of the reconstructions. We tested our hypothesis on the presence of the filament-like structure by creating a synthetic cluster - Irtysh IIIc - mocking the situation of a cluster with external mass. We also computed the source plane magnification distributions. Using them we estimated the probabilities of magnifications in the source plane, and scrutinized their redshift dependence. Finally, we explored the lensing effects of Abell 370 on the luminosity functions of sources at $z_s=9.0$, finding it consistent with published results., Comment: 16 pages, 13 figures

While the kinematics of galaxies up to z~3 have been characterized in detail, only a handful of galaxies at high redshift (z>4) have been examined in such a way. The Atacama Large Millimeter/submillimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE) survey observed a statistically significant sample of 118 star-forming main sequence galaxies at z=4.4-5.9 in [CII]158um emission, increasing the number of such observations by nearly 10x. A preliminary qualitative classification of these sources revealed a diversity of kinematic types (i.e., rotators, mergers, and dispersion-dominated systems). In this work, we supplement the initial classification by applying quantitative analyses to the ALPINE data: a tilted ring model (TRM) fitting code (3DBarolo), a morphological classification (Gini-M20), and a set of disk identification criteria. Of the 75 [CII]-detected ALPINE galaxies, 29 are detected at sufficient significance and spatial resolution to allow for TRM fitting and the derivation of morphological and kinematic parameters. These 29 sources constitute a high-mass subset of the ALPINE sample (M_*>10^9.5Msol). We robustly classify 14 of these sources (six rotators, five mergers, and three dispersion-dominated systems); the remaining sources showing complex behaviour. By exploring the G-M20 of z>4 rest-frame FIR and [CII] data for the first time, we find that our 1"~6kpc resolution data alone are insufficient to separate galaxy types. We compare the rotation curves and dynamical mass profiles of the six ALPINE rotators to the two previously detected z~4-6 unlensed main sequence rotators, finding high rotational velocities (~50-250km/s) and a diversity of rotation curve shapes., Comment: 25 pages, 15 figures. Accepted for publication in MNRAS

The faint-end slope of the quasar luminosity function at z~6 and its implication on the role of quasars in reionizing the intergalactic medium at early times has been an outstanding problem for some time. The identification of faint high-redshift quasars with luminosities of <1e44.5 erg/s is challenging. They are rare (few per square-degree) and the separation of these unresolved quasars from late-type stars and compact star-forming galaxies is difficult from ground-based observations alone. In addition, source confusion becomes significant at >25mag, with 30% of sources having their flux contaminated by foreground objects when the seeing resolution is ~0.7". We mitigate these issues by performing a pixel-level joint processing of ground and space-based data from Subaru/HSC and HST/ACS. We create a deconfused catalog over the 1.64 square-degrees of the COSMOS field, after accounting for spatial varying PSFs and astrometric differences between the two datasets. We identify twelve low-luminosity (M_UV ~ -21 mag) z>6 quasar candidates through (i) their red color measured between ACS/F814W and HSC/i-band and (ii) their compactness in the space-based data. Non-detections of our candidates in Hubble DASH data argues against contamination from late-type stars. Our constraints on the faint end of the quasar luminosity function at z~6.4 suggests a negligibly small contribution to reionization compared to the star-forming galaxy population. The confirmation of our candidates and the evolution of number density with redshift could provide better insights into how supermassive galaxies grew in the first billion years of cosmic time., Comment: 24 pages (incl. Appendix), 21 figures, 1 table; accepted for publication in ApJ

We examine the Fundamental Plane (FP) and mass-to-light ratio ($M/L$) scaling relations using the largest sample of massive quiescent galaxies at $1.511.26$, subset of 8 quiescent galaxies at $z>2$, from Stockmann et al. (2020), we show that they cannot passively evolve to the local Coma cluster relation alone and must undergo significant structural evolution to mimic the sizes of local massive galaxies. The evolution of the FP and $M/L$ scaling relations, from $z=2$ to present-day, for this subset are consistent with passive aging of the stellar population and minor merger structural evolution into the most massive galaxies in the Coma cluster and other massive elliptical galaxies from the MASSIVE Survey. Modeling the luminosity evolution from minor merger added stellar populations favors a history of merging with "dry" quiescent galaxies., Comment: 17 pages, 6 figures

How and when did galaxies form and assemble their stars and stellar mass? The answer to these questions, so crucial to astrophysics and cosmology, requires the full reconstruction of the so called cosmic star formation rate density (SFRD), i.e. the evolution of the average star formation rate per unit volume of the universe. While the SFRD has been reliably traced back to 10-11 billion years ago, its evolution is still poorly constrained at earlier cosmic epochs, and its estimate is mainly based on galaxies luminous in the ultraviolet and with low obscuration by dust. This limited knowledge is largely due to the lack of an unbiased census of all types of star-forming galaxies in the early universe. We present a new approach to find dust-obscured star-forming galaxies based on their emission at radio wavelengths coupled with the lack of optical counterparts. Here, we present a sample of 197 galaxies selected with this method. These systems were missed by previous surveys at optical and near-infrared wavelengths, and 22 of them are at very high redshift (i.e. z > 4.5). The contribution of these elusive systems to the SFRD is substantial and can be as high as 40% of the previously known SFRD based on UV-luminous galaxies. The mere existence of such heavily obscured galaxies in the first two billion years after the Big Bang opens new avenues to investigate the early phases of galaxy formation and evolution, and to understand the links between these systems and the massive galaxies which ceased their star formation at later cosmic times., Comment: 15 pages, 6 figures. Accepted for publication in ApJ. Corrected a typo in Table 3. Minor changes to match the accepted version

The Euclid, Rubin/LSST and Roman (WFIRST) projects will undertake flagship optical/near-infrared surveys in the next decade. By mapping thousands of square degrees of sky and covering the electromagnetic spectrum between 0.3 and 2 microns with sub-arcsec resolution, these projects will detect several tens of billions of sources, enable a wide range of astrophysical investigations by the astronomical community and provide unprecedented constraints on the nature of dark energy and dark matter. The ultimate cosmological, astrophysical and time-domain science yield from these missions will require joint survey processing (JSP) functionality at the pixel level that is outside the scope of the individual survey projects. The JSP effort scoped here serves two high-level objectives: 1) provide precise concordance multi-wavelength images and catalogs over the entire sky area where these surveys overlap, which accounts for source confusion and mismatched isophotes, and 2) provide a science platform to analyze concordance images and catalogs to enable a wide range of astrophysical science goals to be formulated and addressed by the research community. For the cost of about 200WY, JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also allow for the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone., Comment: Final Report from March 2019 submitted to NASA, NSF, DOE, 41 pages. This is the extended version of the Astro2020 white papers

Star formation rate (SFR) measurements at z>4 have relied mostly on rest-frame far-ultraviolet (FUV) observations. The corrections for dust attenuation based on IRX-$\beta$ relation are highly uncertain and are still debated in the literature. Hence, rest-frame far-infrared (FIR) observations are necessary to constrain the dust-obscured component of the SFR. In this paper, we exploit the rest-frame FIR continuum observations collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) to directly constrain the obscured SFR in galaxies at 4.4

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

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

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

The infrared spectral energy distributions (SEDs) of main-sequence galaxies in the early universe (z > 4) is currently unconstrained as infrared continuum observations are time consuming and not feasible for large samples. We present Atacama Large Millimetre Array (ALMA) Band 8 observations of four main-sequence galaxies at z ~ 5.5 to study their infrared SED shape in detail. Our continuum data (rest-frame 110$\rm \mu m$, close to the peak of infrared emission) allows us to constrain luminosity weighted dust temperatures and total infrared luminosities. With data at longer wavelengths, we measure for the first time the emissivity index at these redshifts to provide more robust estimates of molecular gas masses based on dust continuum. The Band 8 observations of three out of four galaxies can only be reconciled with optically thin emission redward of rest-frame 100$\rm \mu m$. The derived dust peak temperatures at z ~ 5.5 (38$\pm$8K) are elevated compared to average local galaxies, however, 5-10K below what would be predicted from an extrapolation of the trend at $z<4$. This behaviour can be explained by decreasing dust abundance (or density) towards high redshifts, which would cause the infrared SED at the peak to be more optically thin, making hot dust more visible to the external observer. From the 850$\rm \mu m$ dust continuum, we derive molecular gas masses between $10^{10}$ and $10^{11}\,{\rm M_{\odot}}$ and gas fractions (gas over total mass) of 30-80% (gas depletion times of 100-220Myrs). All in all, our results provide a first measured benchmark SED to interpret future millimetre observations of normal, main-sequence galaxies in the early Universe., Comment: 12 pages, 8 Figures, 3 Table, Submitted to MNRAS

We present ALMA observations of a merging system at z ~ 4.57, observed as a part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. Combining ALMA [CII] 158 micron and far-infrared continuum data with multi-wavelength ancillary data we find that the system is composed of two massive (Mstar >~ 10^10 Msun) star-forming galaxies experiencing a major merger (stellar mass ratio r_mass ~ 0.9) at close spatial (~13 kpc; projected) and velocity (delta_v < 300 km/s) separations, and two additional faint narrow [CII]-emitting satellites. The overall system belongs to a larger-scale protocluster environment and is coincident to one of its overdensity peaks. ALMA reveals also the presence of [CII] emission arising from a circumgalactic gas structure, extending up to a diameter-scale of ~30 kpc. Our morpho-spectral decomposition analysis shows that about 50% of the total flux resides between the individual galaxy components, in a metal-enriched gaseous envelope characterized by a disturbed morphology and complex kinematics. Similarly to observations of shock-excited [CII] emitted from tidal tails in local groups, our results can be interpreted as a possible signature of interstellar gas stripped by strong gravitational interactions, with a possible contribution from material ejected by galactic outflows and emission triggered by star formation in small faint satellites. Our findings suggest that mergers could be an efficient mechanism of gas mixing in the circumgalactic medium around high-z galaxies, and thus play a key role in the galaxy baryon cycle at early epochs., Comment: 12 pages, 7 figures. Accepted for publication in A&A

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

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