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

Recent observations reveal a rapid dust build-up in high-redshift galaxies (z > 4), challenging current models of galaxy formation. While our understanding of dust production and destruction in the interstellar medium (ISM) is advancing, probing baryonic processes in the early Universe remains a complex task. We characterize the evolution of 98 z~5 star-forming galaxies observed as part of the ALPINE survey by constraining the physical processes underpinning the gas and dust production, consumption, and destruction in their ISM. We make use of chemical evolution models to simultaneously reproduce the observed dust and gas content. For each galaxy, we estimate initial gas mass, inflows and outflows, and efficiencies of dust growth and destruction. We test the models with the canonical Chabrier and top-heavy initial mass functions (IMFs), with the latter enabling rapid dust production on shorter timescales. Our models successfully reproduce gas and dust content in older galaxies (> 600 Myr) regardless of the IMF, with Type II SNe as the primary dust source and no dust growth in ISM with moderate inflow of primordial gas. In case of intermediate-age galaxies (300 - 600 Myr), we reproduce the gas and dust content through Type II SNe and dust growth in ISM, though we observe an over-prediction of dust mass in older galaxies, potentially indicating an unaccounted dust destruction mechanism and/or an overestimation of the observed dust masses. The number of young galaxies (< 300 Myr) reproduced, increases for models assuming top-heavy IMF but with maximal prescriptions of dust production. Galactic outflows are necessary to reproduce observed gas and dust masses. The Chabrier IMF models reproduce 65% of galaxies, while top-heavy IMF models improve this to 93%, easing tensions with observations. Upcoming JWST data will refine these models by resolving degeneracies in intrinsic galaxy properties., Comment: Accepted for publication in A&A

The James Webb Space Telescope (JWST) is opening new frontiers of transient discovery and follow-up at high-redshift. Here we present the discovery of a spectroscopically confirmed Type Ia supernova (SN Ia; SN 2023aeax) at z=2.15 with JWST, with cadenced NIRCam observations that enable multi-band light curve fitting. SN 2023aeax lands at the edge of traditional low-z cosmology color cuts because of its blue color (peak rest-frame B-V~-0.3), but we still apply a fiducial standardization approach with the BayeSN model and find that the SN 2023aeax luminosity distance measurement is in agreement (~0.1sigma) with LambdaCDM. SN 2023aeax is only the second spectroscopically confirmed SN Ia in the dark matter-dominated Universe at z>2 (the other is SN 2023adsy), giving it rare leverage to constrain any potential evolution in SN Ia standardized luminosities. Similar to SN 2023adsy (B-V~0.8), SN 2023aeax has a fairly extreme (but opposite) color, which may be due to the small sample size or a secondary factor, such as host galaxy properties. Nevertheless, the SN 2023aeax spectrum is well-represented by normal low-z SN Ia spectra and we find no definitive evolution in SN Ia standardization with redshift. Still, the first two spectroscopically confirmed z>2 SNe Ia have peculiar colors and combine for a ~1sigma distance slope relative to LambdaCDM, the implications of which require a larger sample and dedicated host galaxy observations to investigate., Comment: Submitted to ApJL. arXiv admin note: text overlap with arXiv:2406.05089

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)

In this work, we devise a model for estimating UV and optical line emission (i.e., CIII] $1909$A, H$\beta$, [OIII] $5007$A, H$\alpha$, [NII] $6583$A) tracing HII regions in the interstellar medium (ISM) of galaxies at $z\sim4-6$ from the ALMA Large Programme ALPINE. The aim is to investigate the impact of binary stars in the stellar population along with an abrupt quenching in the Star Formation History (SFH) on line emission. This is crucial for understanding the ISM's properties in early galaxies and identifying new star formation tracers in high-$z$ galaxies. The model simulates HII+Photodissociation Region (PDR) complexes through radiative transfer in 1D slabs, characterized by gas density ($n$), ionisation parameter ($U$), and metallicity ($Z$). It considers: (a) heating from star formation (SF), simulated with Starburst99 and BPASS to quantify binary stars impact; (b) constant, exponentially declining, and quenched SFH scenarios. For each galaxy, we select theoretical ratios from CLOUDY models between [CII] line emission, tracing PDRs, and nebular lines from HII regions, using these to derive expected optical/UV lines from observed [CII]. We find binary stars strongly impact line emission post-quenching, keeping UV photon flux higher for longer, maintaining free electron temperature and ionised column density in HII regions up to 5 Myr after quenching. We constrain ISM properties of our subsample, finding a low ionisation parameter $\log U{\approx}-3.8\pm 0.2$ and moderate/high densities $\log(n/\rm cm^{-3}){\approx}2.9\pm 0.6$. Finally, we derive UV/optical line luminosities-SFR relations for different burstiness parameters ($k_s$). In the fiducial BPASS model, relations show negligible SFH dependence but depend on $k_s$, while in the SB99 case, dependence is on SFH. We propose their use for characterising the burstiness of high-$z$ galaxies., Comment: 14 pages, 8 figures. Accepted for publication in A&A

We study the stellar mass function (SMF) and the co-evolution with dark matter halos via abundance matching in the largest redshift range to date $0.25$, we find increased abundances of massive (log$\, M_{\star}/M_{\odot}>10.5$) implying integrated star formation efficiencies (SFE) $\epsilon_{\star}\equiv M_{\star}\, f_{\rm b}^{-1} M_{\rm halo}^{-1} \gtrsim 0.5$. We find a flattening of the SMF at the high-mass end that is better described by a double power law at $z>5.5$. At $z \lesssim 5.5$ it transitions to a Schechter law which coincides with the emergence of the first massive quiescent galaxies in the Universe. We trace the cosmic stellar mass density (SMD) and infer the star formation rate density (SFRD), which at $z>7.5$ agrees remarkably with recent \JWST{} UV luminosity function-derived estimates. However, at $z \lesssim 3.5$, we find significant tension ($\sim 0.3$ dex) with the cosmic star formation (SF) history from instantaneous SF measures, the causes of which remain poorly understood. We infer the stellar-to-halo mass relation (SHMR) and the SFE from abundance matching out to $z=12$, finding a non-monotonic evolution. The SFE has the characteristic strong dependence with mass in the range of $0.02 - 0.2$, and mildly decreases at the low mass end out to $z\sim3.5$. At $z\sim3.5$ the SFE increases sharply from $\sim 0.1$ to approach high SFE of $0.8-1$ by $z\sim 10$ for log$(M_{\rm h}/M_{\odot})\approx11.5$, albeit with large uncertainties. Finally, we use the SHMR to track the SFE and stellar mass growth throughout the halo history and find that they do not grow at the same rate -- from the earliest times up until $z\sim3.5$ the halo growth rate outpaces galaxy assembly, but at $z>3.5$ halo growth stagnates and accumulated gas reservoirs keep the SF going and galaxies outpace halos.

The stellar mass-star formation rate (M$_\star$ - SFR) plane is essential for distinguishing galaxy populations, but how galaxies move within this plane over cosmic time remains unclear. This study aims to describe galaxy migrations in the M$_\star$ - SFR plane by reconstructing star formation histories (SFHs) for a sample of galaxies out to redshift $z < 4$. This provides insights into the physical processes driving star formation. We use data from the COSMOS field, selecting 299131 galaxies at $z < 4$ with COSMOS-Web NIRCam data (m$_\mathrm{F444W} < 27$) over 0.54 deg$^2$. Using the SED modeling code CIGALE with non-parametric SFHs, we derive physical properties and migration vectors for these galaxies. These vectors describe the direction and velocity of evolutionary paths across the M$_\star$ - SFR plane. To assess the accuracy of these vectors, we compare them to results from the Horizon-AGN simulation. Galaxies within the main sequence show low migration amplitudes and dispersed movement directions, indicating oscillation within the main sequence. Most progenitors were already on the main sequence a billion years earlier. Starburst galaxies assembled half their mass in the last 350 Myr and originated from the main sequence. Passive galaxies show uniformly declining SFHs and include massive galaxies already in the passive region at $z = 3.5-4$, with increasing density over time. Using reconstructed SFHs up to $z < 4$, we propose a coherent picture of how galaxies migrate over cosmic time in the M$_\star$ - SFR plane, highlighting the connection between major phases in the star-formation history.

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

Parametric strong lensing studies of galaxy clusters often display "misleading features". This is the case in the galaxy cluster Abell 370. Using strong lensing techniques, it has been described parametrically by a four dark matter clumps model and galaxy scale perturbers, as well as a significant external shear component, which physical origin remains a challenge. The dark matter distribution features a mass clump with no stellar counterpart and a significant offset between one of the dark matter clumps and its associated stellar counterpart. In this paper, based on BUFFALO data, we begin by revisiting this mass model. We find a four dark matter clumps solution which does not require any external shear and provides a slightly better RMS compared to previous models. Investigating further this new solution, we present a class of models which can accurately reproduce the strong lensing data, but whose parameters for the dark matter component are poorly constrained. We then develop a model where each large scale dark matter component must be associated with a stellar counterpart. This model is unable to reproduce the observational constraints with an RMS smaller than 2.3", and the parameters describing this dark matter component remain poorly constrained. Examining the total projected mass maps, we find a good agreement between the total mass and the stellar distribution, both being bimodal. We interpret the "misleading features" of the four dark matter clumps mass model and the failure of the three dark matter clumps mass model as being symptomatic of the lack of realism of a parametric description of the dark matter distribution, and encourage caution and criticism on the outputs of parametric strong lensing modelling. We briefly discuss the implications of our results for using Abell 370 as a gravitational telescope., Comment: Abstract largely abridged - A&A accepted on 25/10/2024

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

The unique characteristics of the brightest group galaxies (BGGs) link the evolutionary continuum between galaxies like the Milky Way and more massive BCGs in dense clusters. This study investigates the stellar properties of BGGs over cosmic time (z = 0.08-1.30), extending our previous work (Gozaliasl et al. 2016, 2018; Paper I and Paper II). We analyze data of 246 BGGs from our X-ray galaxy group catalog in the COSMOS field, examining stellar age, mass, star formation rate (SFR), specific SFR (sSFR), and halo mass. Comparisons are made with Millennium and Magneticum simulations. We explore the variation of stellar properties with the projected offset from the X-ray peak or host halo center. Using a mock galaxy catalog, we evaluated the accuracy of SED-derived stellar ages, finding a mean absolute error of about one Gyr. Observed BGG age distributions show a bias towards younger ages compared to semi-analytical models and the Magneticum simulation. Our analysis of stellar age versus mass reveals trends with a positive slope, suggesting complex evolutionary pathways across redshifts. We observe a negative correlation between stellar age and SFR across all redshift ranges. Using a cosmic-time-dependent main sequence framework, we identify star-forming BGGs, finding that about 20% of BGGs in the local universe exhibit star-forming characteristics, increasing to 50% at $z=1.0$. Our findings support an inside-out formation scenario for BGGs, where older stellar populations are near the X-ray peak and younger populations at larger offsets indicate ongoing star formation. The distribution of stellar ages for lower-mass BGGs ($10^{10-11} M_\odot$) deviates from constant ages predicted by models, highlighting current models' limitations in capturing galaxies' complex star formation histories., Comment: 24 Pages, 45 figures

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

Galaxy protoclusters are sites of rapid growth, with a high density of massive galaxies driving elevated rates of star formation and accretion onto supermassive black holes. Here, we present new JWST/NIRSpec IFU observations of the J1000+0234 group at $z=4.54$, a dense region of a protocluster hosting a massive, dusty star forming galaxy (DSFG). The new data reveal two extended, high-equivalent-width (EW$_0>1000\r{A}$) [O III] nebulae that appear at both sides of the DSFG along its minor axis (namely O3-N and O3-S). On one hand, the spectrum of O3-N shows a broad and blueshifted component with a full width at half maximum (FWHM) of 1300 km/s, suggesting an outflow origin. On the other hand, O3-S stretches over 8.6 kpc, and has a velocity gradient that spans 800 km/s, but shows no evidence of a broad component. However, both sources seem to be powered by an active galactic nucleus (AGN), so we classified them as extended emission-line regions (EELRs). The strongest evidence comes from the detection of the high-ionization [Ne V] $\lambda 3427$ line toward O3-N, which paired with the lack of hard X-rays implies an obscuring column density above the Compton-thick regime. The [Ne V] line is not detected in O3-S, but we measure a He II $\lambda 4687$/H$\beta$=0.25, which is well above the expectation for star formation. Despite the remarkable alignment of O3-N and O3-S with two radio sources, we do not find evidence of shocks from a radio jet that could be powering the EELRs. We interpret this as O3-S being externally irradiated by the AGN, akin to the famous Hanny's Voorwerp object in the local Universe. In addition, classical line ratio diagnostics (e.g., [O III]/H$\beta$ vs [N II]/H$\alpha$) put the DSFG itself in the AGN region of the diagrams, and therefore suggest it to be the most probable AGN host. These results showcase the ability of JWST to unveil obscured AGN at high redshifts., Comment: 13 pages, 10 figures, and one table. Accepted for publication in A&A on November 12th 2024

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

We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z = 4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find an increased fraction of interacting galaxies in the starburst regime at all redshifts out to z = 4. This increase correlates with star formation efficiency (SFE), but not with gas fraction. However, the correlation is weak (and only significant out to z = 2), which could be explained by the short duration of SFE increase during interaction. In addition, we find that isolated disk galaxies make up a significant fraction of the starburst population. The fraction of such galaxies with star-forming clumps ("clumpy disks") is significantly increased compared to the main-sequence disk population. Furthermore, this fraction directly correlates with SFE. This is direct observational evidence for a long-term increase of SFE maintained due to disk instabilities, contributing to the majority of starburst galaxies in our sample and hence to substantial mass growth in these systems. This result could also be of importance for explaining the growth of the most massive galaxies at z > 6., Comment: 16 pages, 6 figures. Accepted by ApJ

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.

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

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

The recent discovery of a galaxy at z=7.3 with undetected optical emission lines and a blue UV to optical continuum ratio in JWST spectroscopy is surprising and needs to be explained physically. Here, we explore two possibilities that could cause such a seemingly quiescent ~5e8 $\rm M_\odot$ galaxy in the early Universe: (i) stochastic variations in the star formation history (SFH) and (ii) the effect of spatially varying dust attenuation on the measured line and continuum emission properties. Both scenarios can play at the same time to amplify the effect. A stochastic star formation model (similar to realistic SFHs from hydrodynamical simulations of similar-mass galaxies) can create such observed properties if star formation is fast-varying with a correlation time of <150 Myrs given a reasonable burst amplitude of ~0.6 dex. The total time spent in this state is less than 20 Myrs and the likelihood of such a state to occur over 500 Myrs at z=7 is ~50% (consistent with current observations). On the other hand, we show that a spectrum with blue UV continuum and lack of emission lines can be reproduced by a blue+red composite spectrum. The UV continuum is emitted from dust-free density bounded HII regions (blue component), while the red component is a dust-obscured starburst with weakened emission lines due to strong differential dust attenuation between stellar and nebular emission. Future resolving far-infrared observations with ALMA will shed light on the latter scenario., Comment: 8 pages, 5 figures. Accepted for publication in ApJ

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

Context: From redshift 6 to redshift $\approx$ 4 galaxies grow rapidly from low mass galaxies towards the more mature massive galaxies we see at the cosmic noon. Growth via gas accretion and mergers undoubtedly shape this evolution - however, there currently exists much uncertainty over the contribution of each of these processes to the overall evolution of galaxies. Furthermore, previous characterisations of the morphology of galaxies in the molecular gas phase has been limited by the coarse resolution of previous observations. Aims: The goal of this paper is to derive the morpho-kinematic properties of 3 main-sequence systems at $z\sim4.5$, drawn from the ALPINE survey, using brand new high-resolution ALMA data in band 7. The objects were previously characterised as one merger with three components, and and two dispersion-dominated galaxies. Methods: We use intensity and velocity maps, position-velocity diagrams and radial profiles of [CII], in combination with dust continuum maps, to analyse the morphology and kinematics of the 3 systems.} Results: In general, we find that the high-resolution ALMA data reveal more complex morpho-kinematic properties. We identify in one galaxy interaction-induced clumps, showing the profound effect that mergers have on the molecular gas in galaxies, consistent with what is suggested in recent simulations. A galaxy that was previously classified as dispersion dominated turned out to show two bright [CII] emission regions, that could either be merging galaxies or massive star-forming regions within the galaxy itself. The high resolution data for the other dispersion dominated object also revealed clumps of [CII] that were not previously identified. Within the sample, we might also detect star-formation powered outflows (or outflows from Active Galactic Nuclei) which appear to be fuelling diffuse gas regions and enriching the circumgalactic medium., Comment: Submitted to A&A, 13 pages

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

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

While space-borne optical and near-infrared facilities have succeeded in delivering a precise and spatially resolved picture of our Universe, their small survey area is known to under-represent the true diversity of galaxy populations. Ground-based surveys have reached comparable depths but at lower spatial resolution, resulting in source confusion that hampers accurate photometry extractions. What once was limited to the infrared regime has now begun to challenge ground-based ultra-deep surveys, affecting detection and photometry alike. Failing to address these challenges will mean forfeiting a representative view into the distant Universe. We introduce The Farmer: an automated, reproducible profile-fitting photometry package that pairs a library of smooth parametric models from The Tractor (Lang et al. 2016) with a decision tree that determines the best-fit model in concert with neighboring sources. Photometry is measured by fitting the models on other bands leaving brightness free to vary. The resulting photometric measurements are naturally total, and no aperture corrections are required. Supporting diagnostics (e.g. $\chi^2$) enable measurement validation. As fitting models is relatively time intensive, The Farmer is built with high-performance computing routines. We benchmark The Farmer on a set of realistic COSMOS-like images and find accurate photometry, number counts, and galaxy shapes. The Farmer is already being utilized to produce catalogs for several large-area deep extragalactic surveys where it has been shown to tackle some of the most challenging optical and near-infrared data available, with the promise of extending to other ultra-deep surveys expected in the near future. The Farmer is available to download from GitHub and Zenodo., Comment: 30 pages, 17 figures, accepted for publication in ApJS. The Farmer software is publicly accessible on Github at https://github.com/astroweaver/the_farmer

Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED fitting codes. Using JWST/NIRCam images, we also produce two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. The lens is a very massive and quiescent (sSFR < 10^(-13) yr-1) elliptical galaxy at z = 2.02 \pm 0.02 with a total mass Mtot(

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

In order to study transient phenomena in the Universe, existing and forthcoming imaging surveys are covering wide areas of sky repeatedly over time, with a range of cadences, point spread functions, and depths. We describe here a framework that allows an efficient search for different types of time-varying astrophysical phenomena in current and future, large data repositories. We first present a methodology to generate and store key survey parameters that enable researchers to determine if a survey, or a combination of surveys, allows specific time-variable astrophysical phenomena to be discovered. To facilitate further exploration of sources in regions of interest, we then generate a few sample metrics that capture the essential brightness characteristics of a sky pixel at a specific wavelength. Together, we refer to these as "annotated coadds". The techniques presented here for WISE/NEOWISE-R data are sensitive to 10 percent brightness variations at around 12th Vega magnitude at 4.5 microns wavelength. Application of the technique to ZTF data also enabled the detection of 0.5 mag variability at 20 AB mag in the r-band. We demonstrate the capabilities of these metrics for different classes of sources: high proper-motion stars, periodic variable stars, and supernovae, and find that each metric has its advantages depending on the nature of variability. We also present a data structure which will ease the search for temporally varying phenomena in future surveys., Comment: 26 pages, 11 figures, PASP in press

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