Search

We explore the evolving relationship between galaxies and their dark matter halos from $z \sim 0.1$ to $z \sim 12$ using mass-limited angular clustering measurements in the 0.54 deg$^2$ of the COSMOS-Web survey. This study provides the first measurements of the mass-limited two-point correlation function at $z \ge 10$ and a consistent analysis spanning 13.4 Gyr of cosmic history, setting new benchmarks for future simulations and models. Using a halo occupation distribution (HOD) framework, we derive characteristic halo masses and the stellar-to-halo mass relationship (SHMR) across redshifts and stellar mass bins. Our results first indicate that HOD models fit data at $z \ge 2.5$ best when incorporating a non-linear scale-dependent halo bias, boosting clustering at non-linear scales (r = 10-100 kpc). We find that galaxies at z > 10.5 with $\log(M_\star / M_\odot) \ge 8.85$ are hosted by halos with $M_{\rm h} \sim 10^{10.5}\,M_\odot$, achieving a star formation efficiency (SFE) $M_\star / (f_b M_{\rm h}) $ up to 1 dex higher than at $z \le 1$. The high galaxy bias at $z \ge 8$ suggests that these galaxies reside in massive halos with intrinsic high SFE. Our SHMR evolves significantly with redshift, starting high at $z \ge 10.5$, decreasing until $z \sim 2 - 3$, then increasing again until the present. Current simulations fail to reproduce both massive high-$z$ galaxies and this evolution, while semi-empirical models linking SFE to halo mass, accretion rates, and redshift align with our findings. We propose that $z > 8$ galaxies experience bursty star formation without significant feedback altering their growth, driving the rapid growth of massive galaxies observed by JWST. Over time, increasing feedback efficiency and exponential halo growth suppress star formation. At $z \sim 2 - 3$ and after, halo growth slows down while star formation continues, supported by gas reservoirs in halos., Comment: Submitted to Astronomy & Astrophysics

Galaxy groups with $M_{tot} \lesssim 10^{14}$ $M_\odot$ and up to a few tens of members are the most common galaxy environment, marking the transition between field and massive clusters. Identifying groups plays a crucial role in understanding structure formation and galaxy evolution. Modern deep surveys allow us to build well-characterized samples of groups up to the regime where structures were taking shape. We aimed to build the largest deep catalog of galaxy groups to date over the COSMOS-Web field effective area of 0.45 deg$^2$, leveraging the deep high quality data of the new COSMOS-Web photometric catalog resulted from the James Webb Space Telescope observations of the COSMOS-Web field. We performed the group search with the AMICO algorithm, a linear matched filter based on an analytical model for the group signal. AMICO has already been tested in wide and deep field surveys, including COSMOS data up to $z=2$. In this work, we tested the algorithm performances at even higher redshift and searched for protocluster cores at $z>2$. We compiled a list of known protoclusters in COSMOS at $2 \leq z \leq 3.7$, matched them with our detections and studied the clustering of the detected cores. We estimated purity and completeness of our sample by creating data-driven mocks with the SinFoniA code and linked signal-to-noise to purity. We detected 1678 groups in the COSMOS-Web field up to $z=3.7$, including lists of members extending nearly two magnitudes deeper than the previous AMICO-COSMOS catalog. 756 groups were detected with purity of 80\%. More than 500 groups have their redshift confirmed by assigning spectroscopic counterparts. This group catalog offers a unique opportunity to explore galaxy evolution in different environments spanning $\sim$12 Gyr and to study groups, from the least rich population to the formation of the most massive clusters., Comment: 17 pages, 14 figures, 2 tables (including appendices), submitted to Astronomy & Astrophysics (A&A)

Star-forming clumps have been found to significantly influence the star formation of gas-rich $z>1$ galaxies. Using public data from JWST/NIRCam (COSMOS-Web) and ALMA (FMOS-COSMOS), we study a sample of 32 massive ($>10^{10.5}\,\rm M_{\odot}$) main-sequence galaxies at $z_{\rm spec}\sim1.5$ with $\sim0.3\,\rm kpc$ resolution. We create composite morphological models consisting of bulge, disk, and clumps to fully 'deconstruct' the galaxy images. With the resulting measurements of the flux and size of these components, we find the following: (I)The combined contribution of clumps is $1-30\%$ towards the net star formation rate (SFR) of the host while contributing $1-20\%$ to its stellar mass. The clumps show a correlation between their stellar mass and SFR, but have an increased specific-SFR (sSFR) relative to the star-forming main sequence, with offsets ranging from $0\lesssim\Delta\log\rm sSFR\lesssim 0.4$. They feature star formation surface densities of $10^{-2}-10^{2}\,\rm M_{\odot}/yr/kpc^{2}$, consistent with values observed in local star-forming and starburst galaxies. (II)The clumps span a large range of characteristic sizes ($r_{e}\sim0.1-1\,\rm kpc$) and stellar masses ($\sim 10^{8.0-9.5}\,\rm M_{\odot}$). We estimate a mass-size relation ($r_{e}\propto\rm M_{\star}^{\,0.52\pm0.07}$) along with a stellar mass function (slope, $\alpha=-1.85\pm 0.19$), both suggesting a hierarchical nature similar to that expected in star-forming regions in local galaxies. (III)Our measurements agree with the properties of stellar clumps in $z\gtrsim1$ lensed systems, bridging the gap between lensed and unlensed studies by detecting structures at sub-kpc scales.(IV)Clumps are found to be preferentially located along spiral features visible primarily in the residual rest-frame near-IR images. In conclusion, we present an observation-based, coherent picture of star-forming clumps at $z>1$., Comment: Published in MNRAS; 22 pages, 24 figures

Recent JWST observations have revealed the prevalence of spiral structures at $z > 1$. Unlike in the local Universe, the origin and the consequence of spirals at this epoch remain unexplored. We use public JWST/NIRCam data from the COSMOS-Web survey to map spiral structures in eight massive ($> 10^{10.5}\,\rm M_{\odot}$) star-forming galaxies at $z_{\rm spec} \sim 1.5$. We present a method for systematically quantifying spiral arms at $z>1$, enabling direct measurements of flux distributions. Using rest-frame near-IR images, we construct morphological models accurately tracing spiral arms. We detect offsets ($\sim 0.2 - 0.8\,\rm kpc$) between the rest-frame optical and near-IR flux distributions across most arms. Drawing parallels to the local Universe, we conclude that these offsets reflect the presence of density waves. For nine out of eighteen arms, the offsets indicate spiral shocks triggered by density waves. Five arms have offsets in the opposite direction and are likely associated with tidal interactions. For the remaining cases with no detected offsets, we suggest that stochastic 'clumpy' star formation is the primary driver of their formation. In conclusion, we find a multi-faceted nature of spiral arms at $z > 1$, similar to that in the local Universe., Comment: Accepted for publication in ApJ Letters; 12 pages, 6 figures

``Little Red Dots'' (LRDs) are an abundant high-redshift population newly discovered by the James Webb Space Telescope (JWST). They are characterized by a red color in the rest-frame optical band, compact morphology, and broad Balmer emission lines (${\rm FWHM} \gtrsim 1000~{\rm km\,s^{-1}}$) that suggest an AGN nature. Using a method of pixel-by-pixel color selection and relaxing the compactness criteria, we identify three of the first dual LRD candidates in the COSMOS-Web survey with projected separations of $0.\!\!^{\prime\prime}2-0.\!\!^{\prime\prime}4$ (1-2 pkpc at their photometric redshifts). A comparison between existing LRD samples and mock data reveals that the projected separations of these dual LRD candidates are unlikely to result from chance projections of objects at different redshifts. In one case (CW-B5-15958), the dual LRD includes two bright sources ($m_{\rm F444W}=24.3$ and $24.8$) with characteristic V-shape spectral energy distribution (SEDs) and photometric redshifts consistent with each other. We find that CW-B5-15958 has a faint off-centered component and a companion galaxy. In the other two dual systems, the brighter LRD exhibits a V-shape SED, while the fainter LRD ($m_{\rm F444W}\gtrsim26$) is undetected in both F115W and F150W. These discoveries suggest that the angular auto-correlation function (ACF) of LRDs exhibits a significant excess ($\sim3\times10^2$ times) on sub-arcsec (kilo-parsec) separations compared to the extrapolation of a power-law ACF of JWST-found AGNs measured over $10^{\prime\prime}-100^{\prime\prime}$. Follow-up spectroscopic confirmation of their redshifts and the construction of a larger sample are essential to advance our understanding of the evolution of supermassive black holes and the importance of mergers in the early universe., Comment: 12 pages, 5 figures, and 1 table. Comments are welcome

We report the discovery of an ultra-massive grand-design red spiral galaxy, named Zh\'ul\'ong (Torch Dragon), at $z_{\rm phot} = 5.2^{+0.3}_{-0.2}$ in the JWST PANORAMIC survey, identified as the most distant bulge+disk galaxy with spiral arms known to date. Zh\'ul\'ong displays an extraordinary combination of properties: 1) a classical bulge centered in a large, face-on exponential stellar disk (half-light radius of $R_{\rm e} = 3.7 \pm 0.1 \, \mathrm{kpc}$), with spiral arms extending across 19 kpc; 2) a clear transition from the red, quiescent core ($F150W-F444W=3.1$ mag) with high stellar mass surface density ($\log(\Sigma M_{\star}/M_{\odot} \, \mathrm{kpc}^{-2}) = 9.91_{-0.09}^{+0.11}$) to the star-forming outer regions, as revealed by spatially resolved SED analysis, which indicates significant inside-out galaxy growth; 3) an extremely high stellar mass at its redshift, with $\log (M_{\star}/M_{\odot})=11.03_{-0.08}^{+0.10}$ comparable to the Milky Way, and an implied baryon-to-star conversion efficiency ($\epsilon \sim 0.3$) that is 1.5 times higher than even the most efficient galaxies at later epochs; 4) despite an active disk, a relatively modest overall star formation rate ($\mathrm{SFR} =66_{-46}^{+89} ~M_{\odot} \, \mathrm{yr}^{-1}$), which is $>$0.5 dex below the star formation main sequence at $z \sim 5.2$ and $>$10 times lower than ultra-massive dusty galaxies at $z=5-6$. Altogether, Zh\'ul\'ong shows that mature galaxies emerged much earlier than expected in the first billion years after the Big Bang through rapid galaxy formation and morphological evolution. Our finding offers key constraints for models of massive galaxy formation and the origin of spiral structures in the early universe., Comment: 11 pages, 5 main figures, 2 extended figures, 1 table; submitted to A&A

Optically dark dusty star-forming galaxies (DSFGs) play an essential role in massive galaxy formation at early cosmic time, however their nature remains elusive. Here we present a detailed case study of all the baryonic components of a $z=4.821$ DSFG, XS55. Selected from the ultra-deep COSMOS-XS 3GHz map with a red SCUBA-2 450$\mu$m/850$\mu$m colour, XS55 was followed up with ALMA 3mm line scans and spectroscopically confirmed to be at $z=4.821$ via detections of the CO(5-4) and [CI](1-0) lines. JWST/NIRCam imaging reveals that XS55 is a F150W-dropout with red F277W/F444W colour, and a complex morphology: a compact central component embedded in an extended structure with a likely companion. XS55 is tentatively detected in X-rays with both Chandra and XMM-Newton, suggesting an active galactic nucleus (AGN) nature. By fitting a panchromatic SED spanning NIR to radio wavelengths, we revealed that XS55 is a massive main-sequence galaxy with a stellar mass of $M_\ast=(5\pm1)\times10^{10}\,{\rm M_\odot}$ and a star formation rate of ${\rm SFR}=540\pm177~{\rm M_\odot\,yr^{-1}}$. The dust of XS55 is optically thick in the far infrared (FIR) with a surprisingly cold dust temperature of $T_{\rm dust}=33\pm2\,{\rm K}$, making XS55 one of the coldest DSFGs at $z>4$ known to date. This work unveils the nature of a radio-selected F150W-dropout, suggesting the existence of a population of DSFGs hosting active black holes embedded in optically thick dust., Comment: 11 pages, 7 figures, accepted in Astronomy & Astrophysics

We present a detailed analysis of EELG1002: a $z = 0.8275$ EELG identified within archival Gemini/GMOS spectroscopy as part of the COSMOS Spectroscopic Archive. Combining GMOS spectra and available multi-wavelength photometry, we find EELG1002 is a low-mass ($10^{7 - 8}$ M$_\odot$), compact ($\sim 530$ pc), and bursty star-forming galaxy with mass doubling timescales of $\sim 5 - 15$ Myr. EELG1002 has record-breaking rest-frame [OIII]+H$\beta$ EW of $\sim 2800 - 3700$\AA~which is $\sim 16 - 35 \times$ higher than typical $z \sim 0.8$ [OIII] emitters with similar stellar mass and even higher than typical $z > 5$ galaxies. We find no clear evidence of an AGN suggesting the emission lines are star formation driven. EELG1002 is chemically unevolved (direct $T_e$; $12+\log_{10} (\textrm{O/H}) \sim 7.5$ consistent with $z > 5$ galaxies at fixed stellar mass) and may be undergoing a first intense, bursty star formation phase analogous to conditions expected of galaxies in the early Universe. We find evidence for a highly energetic ISM ([OIII]/[OII] $\sim 11$) and hard ionizing radiation field (elevated [NeIII]/[OII] at fixed [OIII]/[OII]). Coupled with its compact, metal-poor, and actively star-forming nature, EELG1002 is found to efficiently produce ionizing photons with $\xi_{ion} \sim 10^{25.70 - 25.75}$ erg$^{-1}$ Hz and may have $\sim 10 - 20\%$ LyC escape fraction suggesting such sources may be important reionization-era analogs. We find dynamical mass of $\sim 10^9$ M$_\odot$ suggesting copious amounts of gas to support intense star-formation activity as also suggested by analogs identified in Illustris-TNG. EELG1002 may be an ideal low-$z$ laboratory of galaxies in the early Universe and demonstrates how archival datasets can support high-$z$ science and next-generation surveys planned with \textit{Euclid} and \textit{Roman}., Comment: 27 pages, 9 Figures, 8 Tables. Submitted to ApJ. Abstract is abridged. Comments are welcomed. All code can be accessed via our github repository (https://github.com/akhostov/EELG1002)

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^{12} 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 17$ 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 $\lesssim 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: Accepted to MNRAS

Photometric redshifts are widely used in studies of dusty star-forming galaxies (DSFGs), but catastrophic photo-$z$ failure can undermine all redshift-dependent results. Here we report the spectroscopic redshift confirmation of COSBO-7, a strongly lensed DSFG in the COSMOS-PRIMER field. Recently, a photometric redshift solution of $z\gtrsim7.0$ was reported for COSBO-7 based on ten bands of {\it James Webb} Space Telescope (JWST) NIRCam and MIRI imaging data. This $z$ value was favored by four independent spectral energy distribution (SED) fitting codes, and the result provided an appealing candidate for the most distant massive DSFG known to date. This photo-$z$ solution was also supported by a single line detection in Atacama Large Millimeter Array (ALMA) Band 3 consistent with CO(7-6) at $z=7.46$. However, our new ALMA observations robustly detect two lines in Band 6 identified as CO(7-6) and [CI](2-1) at $z_{\rm spec}=2.625$, and thus the Band 3 line as CO(3-2). These three robust line detections decidedly place COSBO-7 at $z=2.625$, refuting the photo-$z$ solution. We derive physical parameters by fitting near-infrared(NIR)-to-millimeter(mm) photometry and lens modeling, revealing that COSBO-7 is a main sequence galaxy. We examine possible reasons for this photo-$z$ failure and attribute it to (1) the likely underestimation of photometric uncertainties at 0.9\,$\mu$m and 1.15 \,$\mu$m; and (2) the lack of photometry at wavelengths beyond 20\,$\mu$m. Notably, we recover a bona fide $z_{\rm phot}\sim 2.3$ by including the existing MIPS $24\,\mu$m photometry, demonstrating the critical importance of mid-infrared (MIR) data in bolstering photo-$z$ measurements. This work highlights a common challenge in modeling the SEDs of DSFGs, and provides a cautionary tale regarding the reliability of photometric redshifts as well as pseudo-spectroscopic redshifts based on single line detection., Comment: Published in A&A

JWST has revealed a population of compact and extremely red galaxies at $z>4$, which likely host active galactic nuclei (AGN). We present a sample of 434 ``little red dots'' (LRDs), selected from the 0.54 deg$^2$ COSMOS-Web survey. We fit galaxy and AGN SED models to derive redshifts and physical properties; the sample spans $z\sim5$-$9$ after removing brown dwarf contaminants. We consider two extreme physical scenarios: either LRDs are all AGN, and their continuum emission is dominated by the accretion disk, or they are all compact star-forming galaxies, and their continuum is dominated by stars. If LRDs are AGN-dominated, our sample exhibits bolometric luminosities $\sim10^{45-47}$ erg\,s$^{-1}$, spanning the gap between JWST AGN in the literature and bright, rare quasars. We derive a bolometric luminosity function (LF) $\sim100$ times the (UV-selected) quasar LF, implying a non-evolving black hole accretion density of $\sim10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ from $z\sim2$-$9$. By contrast, if LRDs are dominated by star formation, we derive stellar masses $\sim10^{8.5-10}\,M_\odot$. MIRI/F770W is key to deriving accurate stellar masses; without it, we derive a mass function inconsistent with $\Lambda$CDM. The median stellar mass profile is broadly consistent with the maximal stellar mass surface densities seen in the nearby universe, though the most massive $\sim50$\% of objects exceed this limit, requiring substantial AGN contribution to the continuum. Nevertheless, stacking all available X-ray, mid-IR, far-IR/sub-mm, and radio data yields non-detections. Whether dominated by dusty AGN, compact star-formation, or both, the high masses/luminosities and remarkable abundance of LRDs implies a dominant mode of early galaxy/SMBH growth., Comment: 27 pages, 13 figures. Submitted to ApJ, comments welcome! Data access at https://github.com/hollisakins/akins24_cw

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

SPHEREx, the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and ices Explorer, is a NASA MIDEX mission planned for launch in 2024. SPHEREx will carry out the first all-sky spectral survey at wavelengths between 0.75 micron and 5 micron with spectral resolving power ~40 between 0.75 and 3.8 micron and ~120 between 3.8 and 5 micron At the end of its two-year mission, SPHEREx will provide 0.75-to-5 micron spectra of each 6.2"x6.2" pixel on the sky - 14 billion spectra in all. This paper updates an earlier description of SPHEREx presenting changes made during the mission's Preliminary Design Phase, including a discussion of instrument integration and test and a summary of the data processing, analysis, and distribution plans.

We explore the morphological features and star formation activities of [OII] emitters in the COSMOS UltraDeep field at $z \sim 1.5$ using JWST NIRCam data from the COSMOS-Web survey and Subaru Hyper Suprime-Cam. We also report the discovery of large filamentary structures traced by [OII] emitters, surrounding an extremely overdense core with a galaxy number density $\sim11\times$ higher than the field average. These structures span over 50 cMpc, underscoring their large scale in the cosmic web at this epoch. After matching the stellar mass distributions, the core galaxies show a higher frequency of disturbances (50$\%$ $ \pm$ 9$\%$) than those in outskirts (41$\%$ $ \pm$ 9$\%$) and the field (21$\%$ $ \pm$ 5$\%$), indicative of more frequent mergers and interactions in the innermost $\lesssim1.5 $ arcmin region. Additionally, we observe that specific star formation rates are elevated in denser environments. A Kolmogorov-Smirnov (KS) test comparing the distribution of specific star formation rates of core and field galaxies yields a $\textit{p}$-value of 0.02, suggesting an enhancement of star-formation activity driven by the dense environment. Our findings underscore the environmental impact on galaxy evolution during a pivotal cosmic epoch and set the stage for further investigation with the increasing larger data from upcoming surveys., Comment: 13 pages, 4 figures; accepted for publication in The Astrophysical Journal Letters

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

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

We report the identification of 64 new candidates of compact galaxies, potentially hosting faint quasars with bolometric luminosities of $L_\mathrm{bol} = 10^{43}$--10$^{46}$ erg s$^{-1}$, residing in the reionization epoch within the redshift range of $6 \lesssim z \lesssim 8$. These candidates were selected by harnessing the rich multiband datasets provided by the emerging JWST-driven extragalactic surveys, focusing on COSMOS-Web, as well as JADES, UNCOVER, CEERS, and PRIMER. Our search strategy includes two stages: applying stringent photometric cuts to catalog-level data and detailed spectral energy distribution fitting. These techniques effectively isolate the quasar candidates while mitigating contamination from low-redshift interlopers, such as brown dwarfs and nearby galaxies. The selected candidates indicate physical traits compatible with low-luminosity active galactic nuclei, likely hosting $\approx10^5$--$10^7~M_\odot$ supermassive black holes (SMBHs) living in galaxies with stellar masses of $\approx10^8$--$10^{10}~M_\odot$. The SMBHs selected in this study, on average, exhibit an elevated mass compared to their hosts, with the mass ratio distribution slightly higher than those of galaxies in the local Universe. As with other high-$z$ studies, this is at least in part due to the selection method for these quasars. An extensive Monte Carlo analysis provides compelling evidence that heavy black hole seeds from the direct collapse scenario appear to be the preferred pathway to mature this specific subset of SMBHs by $z\approx7$. This work underscores the significance of further spectroscopic observations, as the quasar candidates presented here offer exceptional opportunities to delve into the nature of the earliest galaxies and SMBHs that formed during cosmic infancy., Comment: Accepted for publication in the Astronomy & Astrophysics journal. 19 pages, 10 figures, and 4 tables. We welcome comments from the reader

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

We present observations and analysis of the starburst, PACS-819, at z=1.45 ($M_*=10^{10.7}$ M$_{ \odot}$), using high-resolution ($0^{\prime \prime}.1$; 0.8 kpc) ALMA and multi-wavelength JWST images from the COSMOS-Web program. Dissimilar to HST/ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such an intense starburst. Through dynamical modeling of the CO J=5--4 emission with ALMA, PACS-819 is rotation-dominated thus has a disk-like nature. However, kinematic anomalies in CO and asymmetric features in the bluer JWST bands (e.g., F150W) support a more disturbed nature likely due to interactions. The JWST imaging further enables us to map the distribution of stellar mass and dust attenuation, thus clarifying the relationships between different structural components, not discernable in the previous HST images. The CO J = 5 -- 4 and FIR dust continuum emission are co-spatial with a heavily-obscured starbursting core (<1 kpc) which is partially surrounded by much less obscured star-forming structures including a prominent arc, possibly a tidally-distorted dwarf galaxy, and a clump, either a sign of an ongoing violent disk instability or a recently accreted low-mass satellite. With spatially-resolved maps, we find a high molecular gas fraction in the central area reaching $\sim3$ ($M_{\text{gas}}$/$M_*$) and short depletion times ($M_{\text{gas}}/SFR\sim$ 120 Myrs) across the entire system. These observations provide insights into the complex nature of starbursts in the distant universe and underscore the wealth of complementary information from high-resolution observations with both ALMA and JWST., Comment: 21 pages(including appendix), 13 figures, accepted by 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

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

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

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

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

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

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