Search

The current paradigm for the co-evolution of galaxies and their supermassive black holes postulates that dust-obscured active galactic nuclei (AGNs) represent a transitional phase towards a more luminous and unobscured state. However, our understanding of dusty AGNs and their host galaxies at early cosmic times is inadequate due to observational limitations. Here, we present JWST observations of CID-931, an X-ray-detected AGN at a spectroscopic redshift of $z_{\rm spec}=4.91$. Multiband NIRCam imaging from the COSMOS-Web program reveals an unresolved red core, similar to JWST-discovered dusty AGNs. Strikingly, the red core is surrounded by at least eight massive star-forming clumps spread over $1.\!\!^{\prime\prime}6 \approx 10~{\rm kpc}$, each of which has a stellar mass of $10^9-10^{10}M_\odot$ and $\sim0.1-1~{\rm kpc}$ in radius. The whole system amounts to $10^{11}M_\odot$ in stellar mass, higher than typical star-forming galaxies at the same epoch. In this system, gas inflows and/or complex merger events may trigger clump formation and AGN activity thus leading to the rapid formation of a massive galaxy hosting a supermassive black hole. Future follow-up observations will provide new insights into the evolution of the galaxy-black hole relationship during such transitional phases in the early universe., Comment: 13 pages, 10 figures, submitted to PASJ

A multitude of JWST studies reveal a surprising over-abundance of over-massive accreting super-massive blackholes (SMBHs) -- leading to a deepening tension between theory and observation in the first billion years of cosmic time. Across X-ray to infrared wavelengths, models built off of pre-JWST predictions fail to easily reproduce observed AGN signatures (or lack thereof), driving uncertainty around the true nature of these sources. Using a sample of JWST AGN identified via their broadened Halpha emission and covered by the deepest X-ray surveys, we find neither any measurable X-ray emission nor any detection of high-ionization emission lines frequently associated with accreting SMBHs. We propose that these sources are accreting at or beyond the Eddington limit, which reduces the need for efficient production of heavy SMBH seeds at cosmic dawn. Using a theoretical model of super-Eddington accretion, we can produce the observed relative dearth of both X-ray and ultraviolet emission, as well as the high Balmer decrements, without the need for significant dust attenuation. This work indicates that super-Eddington accretion is easily achieved through-out the early Universe, and further study is required to determine what environments are required to trigger this mode of black hole growth., Comment: Submitted

One of the greatest challenges in galaxy evolution over the last decade has been constraining the prevalence of heavily dust-obscured galaxies in the early Universe. At $z>3$, these galaxies are increasingly rare, and difficult to identify as they are interspersed among the more numerous dust-obscured galaxy population at $z=1-3$, making efforts to secure confident spectroscopic redshifts expensive, and sometimes unsuccessful. In this work, we present the Extended Mapping Obscuration to Reionization with ALMA (Ex-MORA) Survey -- a 2mm blank-field survey in the COSMOS-Web field, and the largest ever ALMA blank-field survey to-date covering 577 arcmin$^2$. Ex-MORA is an expansion of the MORA survey designed to identify primarily $z>3$ dusty, star-forming galaxies while simultaneously filtering out the more numerous $z<3$ population by leveraging the very negative $K$-correction at observed-frame 2mm. We identify 37 significant ($>$5$\sigma$) sources, 33 of which are robust thermal dust emitters. We measure a median redshift of $\langle z \rangle = 3.6^{+0.1}_{-0.2}$, with two-thirds of the sample at $z>3$, and just under half at $z>4$, demonstrating the overall success of the 2mm-selection technique. The integrated $z>3$ volume density of Ex-MORA sources is $\sim1-3\times10^{-5}$ Mpc$^{-3}$, consistent with other surveys of infrared luminous galaxies at similar epochs. We also find that techniques using rest-frame optical emission (or lack thereof) to identify $z>3$ heavily dust-obscured galaxies miss at least half of Ex-MORA galaxies. This supports the idea that the dusty galaxy population is heterogeneous, and that synergies across observatories spanning multiple energy regimes are critical to understanding their formation and evolution at $z>3$., Comment: 13 pages, 5 figures, submitted to ApJ; fully reduced mosaic will be shared upon publication

JWST has enabled detecting and spatially resolving the heavily dust-attenuated stellar populations of sub-millimeter galaxies, revealing detail that was previously inaccessible. In this work we construct a sample of 289 sub-millimeter galaxies with detailed joint ALMA and JWST constraints in the COSMOS field. Sources are originally selected using the SCUBA-2 instrument and have archival ALMA observations from various programs. Their JWST NIRCam imaging is from COSMOS-Web and PRIMER. We extract multi-wavelength photometry in a manner that leverages the unprecedented near-infrared spatial resolution of JWST, and fit the data with spectral energy distribution models to derive photometric redshifts, stellar masses, star-formation rates and optical attenuation. The sample has an average z=2.6, A_V=2.5, SFR=270 and log(M*)=11.1. There are 81 (30%) galaxies that have no previous optical/near-infrared detections, including 75% of the z>4 sub-sample (n=28). The faintest observed near-infrared sources have the highest redshifts and largest A_V=4. In a preliminary morphology analysis we find that ~10% of our sample exhibit spiral arms and 5% host stellar bars, with one candidate bar found at z>3. Finally, we find that the clustering of JWST galaxies within 10 arcseconds of a sub-mm galaxy is a factor of 2 greater than what is expected based on either random clustering or the distribution of sources around any red galaxy irrespective of a sub-mm detection., Comment: 37 pages (15 for RGBs + references), 14 figures, submitted to ApJ

We present new free-form and hybrid mass reconstructions of the galaxy cluster lens MACS J0416.1-2403 at $z=0.396$ using the lens inversion method GRALE. The reconstructions use 237 spectroscopically confirmed multiple images from Bergamini et. al. 2023 as the main input. Our primary model reconstructs images to a positional accuracy of 0.191", thus representing one of the most precise reconstructions of this lens to date. Our models find broad agreement with previous reconstructions, and identify two $\sim 10^{11} M_{\odot}$ light-unaffiliated substructures. We focus on two highly magnified arcs: Spock and Mothra. Our model features a unique critical curve structure around the Spock arc with 2 crossings. This structure enables sufficient magnification across this arc to potentially explain the large number of transients as microlensing events of supergiant stars. Additionally, we develop a model of the millilens substructure expected to be magnifying Mothra, which may be a binary pair of supergiants with $\mu \sim 6000$. This model accounts for flexibility in the millilens position while preserving the observed flux and minimizing image position displacements along the Mothra arc. We constrain the millilens mass and core radius to $\lesssim 10^6 M_{\odot}$ and $\lesssim 12$ pc, respectively, which would render it one of the smallest and most compact substructures constrained by lensing. If the millilens is dominated by wave dark matter, the axion mass is constrained to be $< 3.0 \times 10^{-21}$ eV. Further monitoring of this lens with JWST will uncover more transients, permitting tighter constraints on the structure surrounding these two arcs., Comment: Submitted. Comments Welcome

Recent stellar chemical abundance measurements of a handful of $z\sim2$ quiescent galaxies have suggested these galaxies exhibit a remarkably strong $\alpha$-enhancement compared to their local and intermediate redshift counterparts. This apparent chemical evolution following quenching suggests that even the innermost regions of massive early-type galaxies may have experienced substantial mixing of stars in mergers, challenging a purely inside-out growth model. However, larger samples are needed to determine whether a high $\alpha$-enhancement ([Mg/Fe] $\approx 0.5$) is common in $z \sim 2$ quiescent galaxies, and a comparative analysis is needed to determine whether it is consistently inferred using different stellar population synthesis models. We report age and stellar chemical abundance measurements for a sample of four gravitationally lensed quiescent galaxies at $z\sim2.1-2.65$ based on Magellan/FIRE spectroscopy. For three of these galaxies we constrain the $\alpha$-enhancement, and in two cases we measure high values comparable to earlier results when the spectra are analyzed consistently. We also find that the choice of modeling approach can exert a significant effect on the measured abundances. This model dependence can be partly, but not entirely, explained by the complex abundance patterns of $\alpha$ elements in galaxies, which has been observed at lower redshifts and in one $z \sim 2$ quiescent galaxy. Our investigation highlights the importance of independently varying abundance of $\alpha$ elements when fitting the spectra of such galaxies. Observations with JWST will soon deliver precise and spatially resolved abundances of these and other quiescent galaxies at cosmic noon, opening a new window into their evolution., Comment: 10 pages, submitted to The Astrophysical Journal

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

Our knowledge of relations between supermassive black holes and their host galaxies at $z\gtrsim1$ is still limited, even though being actively sought out to $z\sim6$. Here, we use the high resolution and sensitivity of JWST to measure the host galaxy properties for 107 X-ray-selected type-I AGNs at $0.68

With their high angular resolutions of 30--100 mas, large fields of view, and complex optical systems, imagers on next-generation optical/near-infrared space observatories, such as the Near-Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST), present both new opportunities for science and also new challenges for empirical point spread function (PSF) characterization. In this context, we introduce ShOpt, a new PSF fitting tool developed in Julia and designed to bridge the advanced features of PIFF (PSFs in the Full Field of View) with the computational efficiency of PSFEx (PSF Extractor). Along with ShOpt, we propose a suite of non-parametric statistics suitable for evaluating PSF fit quality in space-based imaging. Our study benchmarks ShOpt against the established PSF fitters PSFEx and PIFF using real and simulated COSMOS-Web Survey imaging. We assess their respective PSF model fidelity with our proposed diagnostic statistics and investigate their computational efficiencies, focusing on their processing speed relative to the complexity and size of the PSF models. We find that ShOpt can already achieve PSF model fidelity comparable to PSFEx and PIFF while maintaining competitive processing speeds, constructing PSF models for large NIRCam mosaics within minutes., Comment: 53 pages, 27 figures, submitted to Astronomical Journal. Same as v2 with a corrected typo, which contained revisions from v1

In the last decades, an increasing scientific interest has been growing in the elusive population of "dark" (i.e. lacking an optical/NIR counterpart) Dusty Star-Forming Galaxies (DSFGs). Although extremely promising for their likely contribution to the cosmic Star Formation Rate Density and for their possible role in the evolution of the first massive and passive galaxies around $z\sim3$, the difficulty in selecting statistically significant samples of dark DSFGs is limiting their scientific potentialities. This work presents the first panchromatic study of a sample of 263 Radio-Selected NIRdark galaxies discovered in the COSMOS field following the procedure by Talia+21. These sources are selected as radio-bright galaxies (S(3GHz)>12.65 uJy) with no counterpart in the NIR-selected COSMOS2020 catalog (Ks > 25.5 mag). For these sources, we build a new photometric catalog including accurate photometry from the optical to the radio obtained with a new deblending pipeline (PhoEBO: Photometry Extractor for Blended Objects). We employ this catalog to estimate the photo-zs and the physical properties of the galaxies through an SED-fitting procedure performed with two different codes (Magphys and Cigale). Finally, we estimate the AGN contamination in our sample by performing a series of complementary tests. The high values of the median extinction (Av ~ 4) and star formation rate (SFR ~ 500 Msun/yr) confirm the likely DSFG nature of the RS-NIRdark galaxies. The median photo-z (z~3) and the presence of a significant tail of high-z candidates (z>4.5) suggest that these sources are important contributors to the cosmic SFRD and the evolutionary path of galaxies at high redshifts., Comment: 22 pages, 9 figures, 3 tables. Accepted for publication in ApJ

Several large JWST blank field observing programs have not yet discovered the first galaxies expected to form at $15 \leq z \leq 20$. This has motivated the search for more effective survey strategies that will be able to effectively probe this redshift range. Here, we explore the use of gravitationally lensed cluster fields, that have historically been the most effective discovery tool with HST. In this paper, we analyze the effectiveness of the most massive galaxy clusters that provide the highest median magnification factor within a single JWST NIRCam module in uncovering this population. The results of exploiting these lensing clusters to break the $z > 15$ barrier are compared against the results from large area, blank field surveys such as JADES and CEERS in order to determine the most effective survey strategy for JWST. We report that the fields containing massive foreground galaxy clusters specifically chosen to occupy the largest fraction of a single NIRCam module with high magnification factors in the source plane, whilst containing all multiple images in the image plane within a single module provide the highest probability of both probing the $15 \leq z \leq 20$ regime, as well as discovering the highest redshift galaxy possible with JWST. We also find that using multiple massive clusters in exchange for shallower survey depths is a more time efficient method of probing the $z > 15$ regime., Comment: 12 pages, 5 figures

We present the first results from the Web Epoch of Reionization Lyman-$\alpha$ Survey (WERLS), a spectroscopic survey of Lyman-$\alpha$ emission using Keck I/MOSFIRE and LRIS. WERLS targets bright ($J<26$) galaxy candidates with photometric redshifts of $5.5\lesssim z \lesssim 8$ selected from pre-JWST imaging embedded in the Epoch of Reionization (EoR) within three JWST deep fields: CEERS, PRIMER, and COSMOS-Web. Here, we report 11 $z\sim7-8$ Lyman-$\alpha$ emitters (LAEs; 3 secure and 8 tentative candidates) detected in the first five nights of WERLS MOSFIRE data. We estimate our observed LAE yield is $\sim13$%, broadly consistent with expectations assuming some loss from redshift uncertainty, contamination from sky OH lines, and that the Universe is approximately half-ionized at this epoch, whereby observable Lyman-$\alpha$ emission is unlikely for galaxies embedded in a neutral intergalactic medium. Our targets are selected to be UV-bright, and span a range of absolute UV magnitudes with $-23.1 < M_{\text{UV}} < -19.8$. With two LAEs detected at $z=7.68$, we also consider the possibility of an ionized bubble at this redshift. Future synergistic Keck+JWST efforts will provide a powerful tool for pinpointing beacons of reionization and mapping the large scale distribution of mass relative to the ionization state of the Universe., Comment: 27 pages, 8 figures; ApJ submitted

About 12 billion years ago, the Universe was first experiencing light again after the dark ages, and galaxies filled the environment with stars, metals and dust. How efficient was this process? How fast did these primordial galaxies form stars and dust? We can answer these questions by tracing the Star Formation Rate Density (SFRD) back to its widely unknown high redshift tail, traditionally observed in the Near-InfraRed (NIR), Optical and UV bands. Thus, the objects with a high amount of dust were missing. We aim to fill this knowledge gap by studying Radio Selected NIR-dark (\textit{RS-NIRdark}) sources, i.e. sources not having a counterpart at UV-to-NIR wavelengths. We widen the sample by Talia et al. (2021) from 197 to 272 objects in the COSMic evolution Survey (COSMOS) field, including also photometrically contaminated sources, previously excluded. Another important step forward consists in the visual inspection of each source in the bands from u* to MIPS-24$\mu$m. According to their "environment" in the different bands, we are able to highlight different cases of study and calibrate an appropriate photometric procedure for the objects affected by confusion issues. We estimate that the contribution of RS-NIRdark to the Cosmic SFRD at 3$<$z$<$5 is $\sim$10--25$\%$ of that based on UV-selected galaxies.

In this letter, we report the discovery of the highest redshift, heavily obscured, radio-loud AGN candidate selected using JWST NIRCam/MIRI, mid-IR, sub-mm, and radio imaging in the COSMOS-Web field. Using multi-frequency radio observations and mid-IR photometry, we identify a powerful, radio-loud (RL), growing supermassive black hole (SMBH) with significant spectral steepening of the radio SED ($f_{1.28 \mathrm{GHz}} \sim 2$ mJy, $q_{24\mu m} = -1.1$, $\alpha_{1.28-3\mathrm{GHz}}=-1.2$, $\Delta \alpha = -0.4$). In conjunction with ALMA, deep ground-based observations, ancillary space-based data, and the unprecedented resolution and sensitivity of JWST, we find no evidence of AGN contribution to the UV/optical/NIR data and thus infer heavy amounts of obscuration (N$_{\mathrm{H}} > 10^{23}$ cm$^{-2}$). Using the wealth of deep UV to sub-mm photometric data, we report a singular solution photo-z of $z_\mathrm{phot}$ = 7.7$^{+0.4}_{-0.3}$ and estimate an extremely massive host-galaxy ($\log M_{\star} = 11.4 -12\,\mathrm{M}_{\odot}$) hosting a powerful, growing SMBH (L$_{\mathrm{Bol}} = 4-12 \times 10^{46}$ erg s$^{-1}$). This source represents the furthest known obscured RL AGN candidate, and its level of obscuration aligns with the most representative but observationally scarce population of AGN at these epochs., Comment: Accepted to ApJL

We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ candidate galaxies discovered in the first 0.28 deg$^2$ of JWST/NIRCam imaging from the COSMOS-Web Survey. These sources span rest-frame UV magnitudes of $-20.5>M_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble ACS/F814W, deep ground-based observations corroborate their detection and help significantly constrain their photometric redshifts. We analyze their spectral energy distributions using multiple open-source codes and evaluate the probability of low-redshift solutions; we conclude that 12/15 (80%) are likely genuine $z\gtrsim10$ sources and 3/15 (20%) likely low-redshift contaminants. Three of our $z\sim12$ candidates push the limits of early stellar mass assembly: they have estimated stellar masses $\sim5\times10^{9}\,M_\odot$, implying an effective stellar baryon fraction of $\epsilon_{\star}\sim0.2-0.5$, where $\epsilon_{\star}\equiv M_{\star}/(f_{b}M_{halo})$. The assembly of such stellar reservoirs is made possible due to rapid, burst-driven star formation on timescales $<$100\,Myr where the star-formation rate may far outpace the growth of the underlying dark matter halos. This is supported by the similar volume densities inferred for $M_\star\sim10^{10}\,M_\odot$ galaxies relative to $M_\star\sim10^{9}\,M_\odot$ -- both about $10^{-6}$ Mpc$^{-3}$ -- implying they live in halos of comparable mass. At such high redshifts, the duty cycle for starbursts would be of order unity, which could cause the observed change in the shape of the UVLF from a double powerlaw to Schechter at $z\approx8$. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understanding how, and if, such early massive galaxies push the limits of galaxy formation in $\Lambda$CDM., Comment: 30 pages, 9 figures; ApJ submitted

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

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

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

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

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

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

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

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

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

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

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

Large photometric surveys provide a rich source of observations of quiescent galaxies, including a surprisingly large population at z>1. However, identifying large, but clean, samples of quiescent galaxies has proven difficult because of their near-degeneracy with interlopers such as dusty, star-forming galaxies. We describe a new technique for selecting quiescent galaxies based upon t-distributed stochastic neighbor embedding (t-SNE), an unsupervised machine learning algorithm for dimensionality reduction. This t-SNE selection provides an improvement both over UVJ, removing interlopers which otherwise would pass color selection, and over photometric template fitting, more strongly towards high redshift. Due to the similarity between the colors of high- and low-redshift quiescent galaxies, under our assumptions t-SNE outperforms template fitting in 63% of trials at redshifts where a large training sample already exists. It also may be able to select quiescent galaxies more efficiently at higher redshifts than the training sample., Comment: ApJ, in press

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters., Comment: Accepted ApJS; MAST archive will be live concurrent with publication

We perform a comprehensive analysis of the redshift evolution of the rest-UV spectra of star-forming galaxies out to z~5. We combine new z~5 measurements of HI Ly$\alpha$ and low- and high-ionization interstellar metal absorption features with comparable measurements at z~2-4. We measure the equivalent widths of interstellar absorption features using stacked spectra in bins of Ly$\alpha$ equivalent width, performing corrections to Ly$\alpha$ strengths based on a model for the transmission of the intergalactic medium. We find a strong correlation between decreasing low-ionization absorption strength and increasing Ly$\alpha$ emission strength over the redshift range z~2-5, suggesting that both of these quantities are fundamentally linked to neutral gas covering fraction. At the highest Ly$\alpha$ equivalent widths, we observe evolution at $z\sim5$ towards greater Ly$\alpha$ emission strength at fixed low-ionization absorption strength. If we interpret the non-evolving relationship of Ly$\alpha$ emission strength and low-ionization line strength at z~2-4 as primarily reflecting the radiative transfer of Ly$\alpha$ photons, this evolution at z~5 suggests a higher intrinsic production rate of Ly$\alpha$ photons than at lower redshift. Our conclusion is supported by the joint evolution of the relationships among Ly$\alpha$ emission strength, interstellar absorption strength, and dust reddening. We perform additional analysis in bins of stellar mass, star-formation rate, UV luminosity, and age, examining how the relationships between galaxy properties and Ly$\alpha$ emission evolve towards higher redshift. We conclude that increasing intrinsic Ly$\alpha$ photon production and strong detection of nebular CIV emission (signaling lower metallicity) at z~5 indicate an elevated ionized photon production efficiency ($\xi_{\rm ion}$)., Comment: 20 pages, 17 figures, accepted for publication in MNRAS

We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers., Comment: Published in the Astrophysical Journal letters. Fixed an error in dynamical mass

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves., Comment: 20 pages (+appendix), 11 figures, accepted for publication in the Astrophysical Journal (ApJ). This work comes with a companion letter by Tanaka, Valentino, Toft et al. 2019, ApJL, 885, L34

The redshift range z=4-6 marks a transition phase between primordial and mature galaxy formation in which galaxies considerably increase their stellar mass, metallicity, and dust content. The study of galaxies in this redshift range is therefore important to understand early galaxy formation and the fate of galaxies at later times. Here, we investigate the burstiness of the recent star-formation history (SFH) of 221 $z\sim4.5$ main-sequence galaxies at log(M) > 9.7 by comparing their ultra-violet (UV) continuum, H$\alpha$ luminosity, and H$\alpha$ equivalent-width (EW). The H$\alpha$ properties are derived from the Spitzer [3.6$\mu$m]-[4.5$\mu$m] broad-band color, thereby properly taking into account model and photometric uncertainties. We find a significant scatter between H$\alpha$ and UV-derived luminosities and star-formation rates (SFRs). About half of the galaxies show a significant excess in H$\alpha$ compared to expectations from a constant smooth SFH. We also find a tentative anti-correlation between H$\alpha$ EW and stellar mass, ranging from 1000$\r{A}$ at log(M) < 10 to below 100$\r{A}$ at log(M) > 11. Consulting models suggests that most $z\sim4.5$ galaxies had a burst of star-formation within the last 50 Myrs, increasing their SFRs by a factor of > 5. The most massive galaxies on the other hand might decrease their SFRs, and may be transitioning to a quiescent stage by z=4. We identify differential dust attenuation (f) between stars and nebular regions as the main contributor to the uncertainty. With local galaxies selected by increasing H$\alpha$ EW (reaching values similar to high-z galaxies), we predict that f approaches unity at $z>4$ consistent with the extrapolation of measurements out to z=2., Comment: 31 pages (including appendix), 21 figures, accepted for publication in ApJ

Machine-learning (ML) algorithms will play a crucial role in studying the large datasets delivered by new facilities over the next decade and beyond. Here, we investigate the capabilities and limits of such methods in finding galaxies with brightness-variable active galactic nuclei (AGN). Specifically, we focus on an unsupervised method based on self-organizing maps (SOM) that we apply to a set of nonparametric variability estimators. This technique allows us to maintain domain knowledge and systematics control while using all the advantages of ML. Using simulated light curves that match the noise properties of observations, we verify the potential of this algorithm in identifying variable light curves. We then apply our method to a sample of ~8300 WISE color-selected AGN candidates in Stripe 82, in which we have identified variable light curves by visual inspection. We find that with ML we can identify these variable classified AGN with a purity of 86% and a completeness of 66%, a performance that is comparable to that of more commonly used supervised deep-learning neural networks. The advantage of the SOM framework is that it enables not only a robust identification of variable light curves in a given dataset, but it is also a tool to investigate correlations between physical parameters in multi-dimensional space - such as the link between AGN variability and the properties of their host galaxies. Finally, we note that our method can be applied to any time-sampled light curve (e.g., supernovae, exoplanets, pulsars, and other transient events)., Comment: 9 pages, 4 figures, accepted for publication in ApJL

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$., Comment: 17 pages, 11 figures, Accepted for Publication in ApJ

Using the VLA and ALMA, we have obtained CO(2-1), [C II], [N II] line emission and multiple dust continuum measurements in a sample of "normal" galaxies at $z=5-6$. We report the highest redshift detection of low-$J$ CO emission from a Lyman Break Galaxy, at $z\sim5.7$. The CO line luminosity implies a massive molecular gas reservoir of $(1.3\pm0.3)(\alpha_{\rm CO}/4.5\,M_\odot$ (K km s$^{-1}$ pc$^2)^{-1})\times10^{11}\,M_\odot$, suggesting low star formation efficiency, with a gas depletion timescale of order $\sim$1 Gyr. This efficiency is much lower than traditionally observed in $z\gtrsim5$ starbursts, indicating that star forming conditions in Main Sequence galaxies at $z\sim6$ may be comparable to those of normal galaxies probed up to $z\sim3$ to-date, but with rising gas fractions across the entire redshift range. We also obtain a deep CO upper limit for a Main Sequence galaxy at $z\sim5.3$ with $\sim3$ times lower SFR, perhaps implying a high $\alpha_{\rm CO}$ conversion factor, as typically found in low metallicity galaxies. For a sample including both CO targets, we also find faint [N II] 205$\,\mu$m emission relative to [C II] in all but the most IR-luminous "normal" galaxies at $z=5-6$, implying more intense or harder radiation fields in the ionized gas relative to lower redshift. These radiation properties suggest that low metallicity may be common in typical $\sim$10$^{10}\,M_\odot$ galaxies at $z=5-6$. While a fraction of Main Sequence star formation in the first billion years may take place in conditions not dissimilar to lower redshift, lower metallicity may affect the remainder of the population., Comment: 17 pages, 9 figures, 3 tables, accepted by the Astrophysical Journal

The advent of new sub-millimeter observational facilities has stimulated the desire to model the sub-mm line emission of galaxies within cosmological galaxy formation models. This is typically done in post-processing by applying sub-resolution recipes to describe the properties of the unresolved interstellar medium. At the same time, while there is freedom in how one implements these sorts of recipes, the impact of various choices has yet to be systematically explored in simulations. In this paper, we do just that. We combine a semi-analytic model of galaxy formation with chemical equilibrium networks and numerical radiative transfer models and explore how different choices for the sub-grid modeling affect the predicted CO, [CI], and [CII] emission of galaxies. We find that a key component for a successful model includes a molecular cloud mass-size relation and scaling for the ultraviolet and cosmic ray radiation field that varies based on local ISM properties. Our most successful model adopts a Plummer profile for the radial density distribution of gas within molecular clouds. On the other hand, different assumptions for the clumping of gas within molecular clouds and changes in the slope of the molecular cloud mass distribution function hardly affect the CO, [CI], and [CII] luminosities of galaxies. At fixed star-formation rate the [CII]-SFR ratio of galaxies scales inversely with the pressure acting on molecular clouds, increasing the molecular clouds density and hence decreasing the importance of [CII] line cooling. Overall we find that it is essential that a wide range of sub-mm emission lines arising in vastly different phases of the ISM are used as model constraints in order to limit the freedom in sub-grid choices and we present the details of the most successful model variant., Comment: accepted for publication in MNRAS

We report the serendipitous discovery of a dusty, starbursting galaxy at $z=5.667$ (hereafter called CRLE) in close physical association with the "normal" main-sequence galaxy HZ10 at $z=5.654$. CRLE was identified by detection of [CII], [NII] and CO(2-1) line emission, making it the highest redshift, most luminous starburst in the COSMOS field. This massive, dusty galaxy appears to be forming stars at a rate of at least 1500$\,M_\odot$ yr$^{-1}$ in a compact region only $\sim3$ kpc in diameter. The dynamical and dust emission properties of CRLE suggest an ongoing merger driving the starburst, in a potentially intermediate stage relative to other known dusty galaxies at the same epoch. The ratio of [CII] to [NII] may suggest that an important ($\sim15\%$) contribution to the [CII] emission comes from a diffuse ionized gas component, which could be more extended than the dense, starbursting gas. CRLE appears to be located in a significant galaxy overdensity at the same redshift, potentially associated with a large-scale cosmic structure recently identified in a Lyman Alpha Emitter survey. This overdensity suggests that CRLE and HZ10 reside in a protocluster environment, offering the tantalizing opportunity to study the effect of a massive starburst on protocluster star formation. Our findings support the interpretation that a significant fraction of the earliest galaxy formation may occur from the inside out, within the central regions of the most massive halos, while rapidly evolving into the massive galaxy clusters observed in the local Universe., Comment: 16 pages, 9 figures, 4 tables, final version to appear on ApJ (accepted May 19, 2018)

We trace the specific star formation rate (sSFR) of massive star-forming galaxies ($\gtrsim\!10^{10}\,\mathcal{M}_\odot$) from $z\sim2$ to 7. Our method is substantially different from previous analyses, as it does not rely on direct estimates of star formation rate, but on the differential evolution of the galaxy stellar mass function (SMF). We show the reliability of this approach by means of semi-analytical and hydrodynamical cosmological simulations. We then apply it to real data, using the SMFs derived in the COSMOS and CANDELS fields. We find that the sSFR is proportional to $(1+z)^{1.1\pm0.2}$ at $z>2$, in agreement with other observations but in tension with the steeper evolution predicted by simulations from $z\sim4$ to 2. We investigate the impact of several sources of observational bias, which however cannot account for this discrepancy. Although the SMF of high-redshift galaxies is still affected by significant errors, we show that future large-area surveys will substantially reduce them, making our method an effective tool to probe the massive end of the main sequence of star-forming galaxies., Comment: ApJ accepted

We investigate Lya, [OIII]5007, Ha, and [CII]158um emission from 1124 galaxies at z=4.9-7.0. Our sample is composed of 1092 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 34 galaxies at z=5.148-7.508 with deep ALMA [CII]158um data in the literature. Fluxes of strong rest-frame optical lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, we find that the rest-frame Ha equivalent width and the Lya escape fraction f_Lya positively correlate with the rest-frame Lya equivalent width EW^0_Lya. The f_Lya-EW^0_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5 significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW^0_Lya~0-30 A, and then decreases out to EW^0_Lya~130 A. We also identify an anti-correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW^0_Lya at the >99% confidence level. We carefully investigate physical origins of the correlations with stellar-synthesis and photoionization models, and find that a simple anti-correlation between EW_Lya^0 and metallicity explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study, indicating detections of metal-poor (~0.03 Zo) galaxies with EW^0_Lya~200 A., Comment: 24 pages, 21 figures, accepted in ApJ

We present an empirical parameterization of the [NII]/H$\alpha$ flux ratio as a function of stellar mass and redshift valid at 0 < z < 2.7 and 8.5 < log(M) < 11.0. This description can easily be applied to (i) simulations for modeling [NII]$\lambda6584$ line emission, (ii) deblend [NII] and H$\alpha$ in current low-resolution grism and narrow-band observations to derive intrinsic H$\alpha$ fluxes, and (iii) to reliably forecast the number counts of H$\alpha$ emission-line galaxies for future surveys, such as those planned for Euclid and the Wide Field Infrared Survey Telescope (WFIRST). Our model combines the evolution of the locus on the Baldwin, Phillips & Terlevich (BPT) diagram measured in spectroscopic data out to z = 2.5 with the strong dependence of [NII]/H$\alpha$ on stellar mass and [OIII]/H$\beta$ observed in local galaxy samples. We find large variations in the [NII]/H$\alpha$ flux ratio at a fixed redshift due to its dependency on stellar mass; hence, the assumption of a constant [NII] flux contamination fraction can lead to a significant under- or overestimate of H$\alpha$ luminosities. Specifically, measurements of the intrinsic H$\alpha$ luminosity function derived from current low-resolution grism spectroscopy assuming a constant 29% contamination of [NII] can be overestimated by factors of ~8 at log(L) > 43.0 for galaxies at redshifts z = 1.5. This has implications for the prediction of H$\alpha$ emitters for Euclid and WFIRST. We also study the impact of blended H$\alpha$ and [NII] on the accuracy of measured spectroscopic redshifts., Comment: 18 pages, 11 figures, 1 table, accepted for publication in ApJ

Recent studies have found a significant evolution and scatter in the IRX-$\beta$ relation at z > 4, suggesting different dust properties of these galaxies. The total far-infrared (FIR) luminosity is key for this analysis but poorly constrained in normal (main-sequence) star-forming z > 5 galaxies where often only one single FIR point is available. To better inform estimates of the FIR luminosity, we construct a sample of local galaxies and three low-redshift analogs of z > 5 systems. The trends in this sample suggest that normal high-redshift galaxies have a warmer infrared (IR) SED compared to average z < 4 galaxies that are used as prior in these studies. The blue-shifted peak and mid-IR excess emission could be explained by a combination of a larger fraction of the metal-poor inter-stellar medium (ISM) being optically thin to ultra-violet (UV) light and a stronger UV radiation field due to high star formation densities. Assuming a maximally warm IR SED suggests 0.6 dex increased total FIR luminosities, which removes some tension between dust attenuation models and observations of the IRX-$\beta$ relation at z > 5. Despite this, some galaxies still fall below the minimum IRX-$\beta$ relation derived with standard dust cloud models. We propose that radiation pressure in these highly star-forming galaxies causes a spatial offset between dust clouds and young star-forming regions within the lifetime of O/B stars. These offsets change the radiation balance and create viewing-angle effects that can change UV colors at fixed IRX. We provide a modified model that can explain the location of these galaxies on the IRX-$\beta$ diagram., Comment: 34 pages including Appendix, 18 figures (4 in Appendix), 5 tables (1 in Appendix). Accepted for publication in ApJ

We examine the rest-frame ultra-violet (UV) properties of 10 [CII]$\lambda158\,{\rm \mu m}$$-$detected galaxies at $z\sim5.5$ in COSMOS using new HST/WFC3 near-infrared imaging. Together with pre-existing $158\,{\rm \mu m}-$continuum and [CII] line measurements by ALMA, we study their dust attenuation properties on the IRX-$\beta$ diagram, which connects the total dust emission ($\propto$ IRX=log($L_{FIR}/L_{1600}$)) to the line-of-sight dust column ($\propto\beta$). We find systematically bluer UV continuum spectral slopes ($\beta$) compared to previous low-resolution ground-based measurements, which relieves some of the tension between models of dust attenuation and observations at high redshifts. While most of the galaxies are consistent with local starburst or Small Magellanic cloud like dust properties, we find galaxies with low IRX values and a large range in $\beta$ that cannot be explained by models of a uniform dust distribution well mixed with stars. A stacking analysis of Keck/DEIMOS optical spectra indicates that these galaxies are metal-poor with young stellar populations which could significantly alter their spatial dust distribution., Comment: 9 pages, 3 figures, 1 table. Accepted for publication in The Astrophysical Journal