Search

Quiescent galaxies over $34$, although our measurements are consistent within errors. This technique may be applied to other JWST surveys., Comment: Submitted to The Astrophysical Journal. 19 pages, 7 figures. Code to reproduce Figs. 2 - 4 can be found at https://github.com/alex-delavega/umap_jades_2025. Comments are welcome!

The role of internal and environmental factors in the star formation activity of galaxies is still a matter of debate, particularly at higher redshifts. Leveraging the most recent release of the COSMOS catalog, COSMOS2020, and density measurements from our previous study we disentangle the impact of environment and stellar mass on the star formation rate (SFR), and specific SFR (sSFR) of a sample of $\sim 210,000$ galaxies within redshift range $0.4< z < 4$ and present our findings in three cosmic epochs: 1) out to $z\sim 1$, the average SFR and sSFR decline at extremely dense environments and high mass end of the distribution which is mostly due to the presence of the massive quiescent population; 2) at $1

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides deep HST F275W and F435W imaging over four CANDELS fields (GOODS-N, GOODS-S, COSMOS, and EGS). We combine this newly acquired UV imaging with existing HST imaging from CANDELS as well as existing ancillary data to obtain robust photometric redshifts and reliable estimates for galaxy physical properties for over 150,000 galaxies in the $\sim$430 arcmin$^2$ UVCANDELS area. Here, we leverage the power of the new UV photometry to not only improve the photometric redshift measurements in these fields, but also constrain the full redshift probability distribution combining multiple redshift fitting tools. Furthermore, using the full UV-to-IR photometric dataset, we measure the galaxy physical properties by fitting templates from population synthesis models with two different parameterizations (flexible and fixed-form) of the star-formation histories (SFHs). Compared to the flexible SFH parametrization, we find that the fixed-form SFHs systematically underestimate the galaxy stellar masses, both at the low- ($\lesssim10^9 M_\odot$) and high- ($\gtrsim10^{10} M_\odot$) mass end, by as much as $\sim0.5$ dex. This underestimation is primarily due the limited ability of fixed-form SFH parameterization to simultaneously capture the chaotic nature of star-formation in these galaxies., Comment: 22 pages, 6 figures; accepted to ApJS; catalogs available via MAST

We spectroscopically confirm a new protocluster in the COSMOS field at $z$=2.24430 with Keck/MOSFIRE, dubbed CC2.2B, which is in the immediate vicinity of CC2.2A protocluster, originally presented in \cite{Darvish20}. CC2.2B and CC2.2A centroids are separated by $\sim$5.5 Mpc(angular) and $\sim$16 comoving Mpc(radial). CC2.2B and CC2.2A have similar properties, with CC2.2B having a line-of-sight velocity dispersion and estimated total mass of $\sigma_{los}$=693$\pm$65 km s$^{-1}$ and $M_{total}$=($\sim$2-3)$\times$10$^{14}$ $M_{\odot}$, respectively. These two similar overdensities are likely still in the merging process and will likely collapse into a more massive structure at lower redshifts. We combine CC2.2A and CC2.2B data to investigate the role of high-$z$ protocluster environments on the dynamics of star-forming (SF) galaxies compared to a similarly selected field sample. We find that on average, protocluster SF galaxies have $\sim$0.1 dex (at $\sim$1.8$\sigma$ significance) lower gas velocity dispersions, $\sim$0.2 dex (at $\sim$2.2$\sigma$ significance) lower dynamical masses, and $\sim$0.2 dex lower dynamical-to-stellar mass ratio than the field SF galaxies. We argue that galaxy harassment and galaxy-galaxy interactions can potentially explain these differences. We also find a factor of $\sim$2-3 lower scatter around the mean $\sigma$-$M_{*}$, $M_{dyn}$-$M_{*}$, and $M_{dyn}$/$M_{*}$ vs. $M_{*}$ relations for protocluster SF galaxies than the field. This could be due to a more uniform formation for protocluster galaxies than their field counterparts. Our results have potential implications for the physics of preprocessing in early environments.

The growing volume of data produced by large astronomical surveys necessitates the development of efficient analysis techniques capable of effectively managing high-dimensional datasets. This study addresses this need by demonstrating some applications of manifold learning and dimensionality reduction techniques, specifically the Self-Organizing Map (SOM), on the optical+NIR SED space of galaxies, with a focus on sample comparison, selection biases, and predictive power using a small subset. To this end, we utilize a large photometric sample from the five CANDELS fields and a subset with spectroscopic measurements from the KECK MOSDEF survey in two redshift bins at $z\sim1.5$ and $z\sim2.2$. We trained SOM with the photometric data and mapped the spectroscopic data onto it as our study case. We found that MOSDEF targets do not cover all SED shapes existing in the SOM. Our findings reveal that Active Galactic Nuclei (AGN) within the MOSDEF sample are mapped onto the more massive regions of the SOM, confirming previous studies and known selection biases towards higher-mass, less dusty galaxies. Furthermore, SOM were utilized to map measured spectroscopic features, examining the relationship between metallicity variations and galaxy mass. Our analysis confirmed that more massive galaxies exhibit lower [OIII]/H$\beta$ and [OIII]/[OII] ratios and higher H$\alpha$/H$\beta$ ratios, consistent with the known mass-metallicity relation. These findings highlight the effectiveness of SOM in analyzing and visualizing complex, multi-dimensional datasets, emphasizing their potential in data-driven astronomical studies.

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

We report the identification of two galaxy overdensities at $z\sim5.7$ in the sightline of the galaxy cluster Abell 2744. These overdensities consist of 25 and 17 member galaxies, spectroscopically confirmed with JWST NIRSpec/MSA and NIRCam/WFSS. Each overdensity has a total stellar mass of $\sim2\times10^{10} M_\odot$ and a star formation rate of $\sim200 M_\odot$/yr within a central region of radius $R=2$ Mpc (physical). The sensitive PRISM spectra allow us to identify six galaxies that show weak Ha+[NII] emissions within the overdensities ($27\pm6\%$), whereas the fraction of such galaxies is found significantly lower ($6\pm2\%$) in field samples of the equivalent redshift range. These weak emission line galaxies, dubbed as wELGs, exhibit a strong continuum break at $4000$AA rest-frame, a characteristic feature of evolved stellar populations. The high observed fraction of wELGs in the two overdensities is consistent with the idea that high-density environments are an ideal site where galaxies can accelerate their evolutionary pace compared to field analogs. Our study pinpoints an early onset of environmental effects, already important within one billion years after the Big Bang, and provides a complementary perspective on the emergence of quenched, massive galaxies at lower redshifts. Potential contributions from black hole accretion feedback to the reduction of star formation activity are discussed, but the connection to the local environments remains unclear., Comment: Submitted to ApJ

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 use the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields (UVCANDELS) to measure half-light radii in the rest-frame far-UV for $\sim$16,000 disk-like galaxies over $0.5\leq z \leq 3$. We compare these results to rest-frame optical sizes that we measure in a self-consistent way and find that the stellar mass-size relation of disk galaxies is steeper in the rest-frame UV than in the optical across our entire redshift range. We show that this is mainly driven by massive galaxies ($\gtrsim10^{10}$M$_\odot$), which we find to also be among the most dusty. Our results are consistent with the literature and have commonly been interpreted as evidence of inside-out growth wherein galaxies form their central structures first. However, they could also suggest that the centers of massive galaxies are more heavily attenuated than their outskirts. We distinguish between these scenarios by modeling and selecting galaxies at $z=2$ from the VELA simulation suite in a way that is consistent with UVCANDELS. We show that the effects of dust alone can account for the size differences we measure at $z=2$. This indicates that, at different wavelengths, size differences and the different slopes of the stellar mass-size relation do not constitute evidence for inside-out growth., Comment: Accepted for publication in ApJ. 22 pages, 12 figures, and 4 tables

To study the role of environment in galaxy evolution, we reconstruct the underlying density field of galaxies based on COSMOS2020 (The Farmer catalog) and provide the density catalog for a magnitude limited ($K_{s}<24.5$) sample of $\sim 210 \, k$ galaxies at $0.42$., Comment: 22 pages, 11 figures, 2 tables, Submitted to ApJ

UVCANDELS is a HST Cycle-26 Treasury Program awarded 164 orbits of primary ultraviolet (UV) F275W imaging and coordinated parallel optical F435W imaging in four CANDELS fields: GOODS-N, GOODS-S, EGS, and COSMOS, covering a total area of $\sim426$ arcmin$^2$. This is $\sim2.7$ times larger than the area covered by previous deep-field space UV data combined, reaching a depth of about 27 and 28 ABmag ($5\sigma$ in $0.2"$ apertures) for F275W and F435W, respectively. Along with the new photometric catalogs, we present an analysis of the rest-frame UV luminosity function (LF), relying on our UV-optimized aperture photometry method yielding a factor of $1.5\times$ increase than the H-isophot aperture photometry in the signal-to-noise ratios of galaxies in our F275W imaging. Using well tested photometric redshift measurements we identify 5810 galaxies at redshifts $0.6

We present the results of our Keck/DEIMOS spectroscopic follow-up of candidate galaxies of i-band-dropout protocluster candidate galaxies at $z\sim6$ in the COSMOS field. We securely detect Lyman-$\alpha$ emission lines in 14 of the 30 objects targeted, 10 of them being at $z=6$ with a signal-to-noise ratio of $5-20$, the remaining galaxies are either non-detections or interlopers with redshift too different from $z=6$ to be part of the protocluster. The 10 galaxies at $z\approx6$ make the protocluster one of the riches at $z>5$. The emission lines exhibit asymmetric profiles with high skewness values ranging from 2.87 to 31.75, with a median of 7.37. This asymmetry is consistent with them being Ly$\alpha$, resulting in a redshift range of $z=5.85-6.08$. Using the spectroscopic redshifts, we re-calculate the overdensity map for the COSMOS field and find the galaxies to be in a significant overdensity at the $4\sigma$ level, with a peak overdensity of $\delta=11.8$ (compared to the previous value of $\delta=9.2$). The protocluster galaxies have stellar masses derived from Bagpipes SED fits of $10^{8.29}-10^{10.28} \rm \,M_{\rm \odot}$ and star formation rates of $2-39\,\rm M_{\rm \odot}\rm\,yr^{-1}$, placing them on the main sequence at this epoch. Using a stellar-to-halo-mass relationship, we estimate the dark matter halo mass of the most massive halo in the protocluster to be $\sim 10^{12}\rm M_{\rm \odot}$. By comparison with halo mass evolution tracks from simulations, the protocluster is expected to evolve into a Virgo- or Coma-like cluster in the present day., Comment: 26 pages, 14 figues, 5 tables, main text is 16 pages, appendix is 10 pages, published in MNRAS

The connection between galaxies and dark matter halos is often quantified using the stellar mass-halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at $z \gtrsim 2$ can also constrain the SMHM relation via the galaxy autocorrelation function and through the cross-correlation between galaxies and Ly$\alpha$ absorption measured in transverse sightlines; however, such studies are very few and have produced some unexpected or inconclusive results. We use $\sim$3000 spectra of $z\sim2.5$ galaxies from the Lyman-alpha Tomography IMACS Survey (LATIS) to measure the galaxy-galaxy and galaxy-Ly$\alpha$ correlation functions in four bins of stellar mass spanning $10^{9.2} \lesssim M_* / M_{\odot} \lesssim 10^{10.5}$. Parallel analyses of the MultiDark N-body and ASTRID hydrodynamic cosmological simulations allow us to model the correlation functions, estimate covariance matrices, and infer halo masses. We find that results of the two methods are mutually consistent and are broadly in accord with standard SMHM relations. This consistency demonstrates that we are able to accurately measure and model Ly$\alpha$ transmission fluctuations $\delta_F$ in LATIS. We also show that the galaxy-Ly$\alpha$ cross-correlation, a free byproduct of optical spectroscopic galaxy surveys at these redshifts, can constrain halo masses with similar precision to galaxy-galaxy clustering., Comment: Accepted to ApJ

We present the stellar mass - stellar metallicity relation for 3491 star-forming galaxies at $2 \lesssim z \lesssim 3$ using rest-frame far-ultraviolet (FUV) spectra from the Ly$\alpha$ Tomography IMACS Survey (LATIS). We fit stellar population synthesis models from the Binary Population And Spectral Synthesis code (BPASS v$2.2.1$) to medium resolution (R $\sim 1000$) and high signal-to-noise ($>30$ per 100 km/s over a wavelength range of 1221 - 1800 \r{A}) composite spectra of galaxies in bins of stellar mass to determine their stellar metallicity, primarily tracing $\rm Fe/H$. We find a strong correlation between stellar mass and stellar metallicity, with stellar metallicity monotonically increasing with stellar mass at low masses and flattening at high masses ($M_* \gtrsim 10^{10.3} M_\odot$). Additionally, we compare our stellar metallicity measurements with the gas-phase oxygen abundance of galaxies at similar redshift and estimate the average $\rm [\alpha/Fe] \sim 0.6$. Such high $\alpha$-enhancement indicates that high-redshift galaxies have not yet undergone significant iron enrichment through Type Ia supernovae. Moreover, we utilize an analytic chemical evolution model to constrain the mass loading parameter of galactic winds as a function of stellar mass. We find that as the stellar mass increases, the mass loading parameter decreases. The parameter then flattens or reaches a turning point at around $M_* \sim 10^{10.5} M_\odot$. Our findings may signal the onset of black hole-driven outflows at $z \sim 2.5$ for galaxies with $M_* \gtrsim 10^{10.5} M_\odot$., Comment: 21 pages, 17 figures, 3 tables, accepted for publication in ApJ

The UltraViolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) survey is a Hubble Space Telescope (HST) Cycle-26 Treasury Program, allocated in total 164 orbits of primary Wide-Field Camera 3 Ultraviolet and Visible light F275W imaging with coordinated parallel Advanced Camera for Surveys F435W imaging, on four of the five premier extragalactic survey fields: GOODS-N, GOODS-S, EGS, and COSMOS. We introduce this survey by presenting a comprehensive analysis of the absolute escape fraction ($f_{\rm esc}^{\rm abs}$) of Lyman continuum (LyC) radiation through stacking the UV images of a population of star-forming galaxies with secure redshifts at $2.4\leq z\leq3.0$. Our stacking benefits from the catalogs of high-quality spectroscopic redshifts compiled from archival ground-based data and HST slitless spectroscopy, carefully vetted by dedicated visual inspection efforts. We develop a robust stacking method to apply to 10 samples of in total 56 galaxies, and perform detailed Monte Carlo (MC) simulations of the intergalactic medium (IGM) attenuation, to take into account the sample variance of the mean IGM transmission when measuring $f_{\rm esc}^{\rm abs}$. The full stack at $z\approx2.44$ from 28 galaxies places a stringent 1-$\sigma$ upper limit of $f_{\rm esc}^{\rm abs}\lesssim5\%$, whereas the full stack at $z\approx2.72$ of equal number of galaxies gives an upper limit of $f_{\rm esc}^{\rm abs}\lesssim26\%$ at 1-$\sigma$ confidence level. These new F275W and F435W imaging mosaics from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes (MAST)., Comment: 28 pages, 7 figures, and 4 tables. Accepted to ApJ

Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, basic demographics of clumps and their host galaxies remain uncertain. Using the HST/WFC3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (UVCANDELS), we detect and analyze giant star-forming clumps in galaxies at 0.5 $\leq$ z $\leq$ 1, connecting two epochs when clumps are common (at cosmic high-noon, z $\sim$ 2) and rare (in the local universe). We construct a clump sample whose rest-frame 1600 {\AA} luminosity is 3 times higher than the most luminous local HII regions (M$_{UV} \leq -$16 AB). In our sample, 35 $\pm$ 3$\%$ of low-mass galaxies (log[M$_{*}$/M$_{\odot}$] $<$ 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22 $\pm$ 3$\%$ and 22 $\pm$ 4$\%$ for intermediate (10 $\leq$ log[M$_{*}$/M$_{\odot}$] $\leq$ 10.5) and high-mass (log[M$_{*}$/M$_{\odot}$] $>$ 10.5) galaxies in agreement with previous studies. When compared to similar-mass non-clumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher SFRs and bluer rest-frame U-V colors, while high-mass clumpy SFGs tend to be larger than non-clumpy SFGs. However, clumpy and non-clumpy SFGs have similar S\'ersic index, indicating a similar underlying density profile. Furthermore, we investigate how UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher-specific SFR galaxies and are found closer to their hosts' galactic center., Comment: 21 pages, 13 figures, accepted for publication in ApJ

High-resolution imaging of galaxies in rest-frame UV has revealed the existence of giant star-forming clumps prevalent in high redshift galaxies. Studying these sub-structures provides important information about their formation and evolution and informs theoretical galaxy evolution models. We present a new method to identify clumps in galaxies' high-resolution rest-frame UV images. Using imaging data from CANDELS and UVCANDELS, we identify star-forming clumps in an HST/F160W$\leq 25$ AB mag sample of 6767 galaxies at $0.5\leq z\leq 3$ in four fields, GOODS-N, GOODS-S, EGS, and COSMOS. We use a low-pass band filter in Fourier space to reconstruct the background image of a galaxy and detect small-scale features (clumps) on the background-subtracted image. Clumpy galaxies are defined as those having at least one off-center clump that contributes a minimum of 10$\%$ of the galaxy's total rest-frame UV flux. We measure the fraction of clumpy galaxies ($\rm f_{clumpy}$) as a function of stellar mass, redshift, and galaxy environment. Our results indicate that $\rm f_{clumpy}$ increases with redshift, reaching $\sim 65\%$ at $z\sim 1.5$. We also find that $\rm f_{clumpy}$ in low-mass galaxies ($\rm 9.5\leq log(M_*/M_\odot)\leq 10$) is 10$\%$ higher compared to that of their high-mass counterparts ($\rm log(M_*/M_\odot)>10.5$). Moreover, we find no evidence of significant environmental dependence of $\rm f_{clumpy}$ for galaxies at the redshift range of this study. Our results suggest that the fragmentation of gas clouds under violent disk instability remains the primary driving mechanism for clump formation, and incidents common in dense environments, such as mergers, are not the dominant processes., Comment: 16 pages, 11 figures, 2 tables, accepted for publication in ApJ

We present the Texas Euclid Survey for Lyman-Alpha (TESLA), a spectroscopic survey in the 10 square degree of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Lyman-alpha emitters (LAEs) correlate with Lyman-alpha emission to understand the escape of Lyman alpha from galaxies at redshifts 2 -- 3.5. We present an analysis of 43 LAEs performed in the NEP field using early data from the TESLA survey. We use Subaru Hyper Suprime-Cam imaging in the grizy-bands, Spitzer/IRAC channels 1 and 2 from the Hawaii 20 square degree (H20) survey and spectra acquired by the Visible Integral-Field Replicable Unit Spectrograph (VIRUS) on the Hobby-Eberly Telescope. We perform spectral energy distribution (SED) fitting to compute the galaxy properties of 43 LAEs, and study correlations between stellar mass, star formation rate (SFR), and dust, to the Lyman-alpha rest-frame equivalent widths (EW). We uncover marginal (1 sigma significance) correlations between stellar mass and Lyman-alpha EW, and star formation rate (SFR) and Lyman-alpha EW, with a Spearman correlation coefficient of -0.$34_{-.14}^{+.17}$ and -0.$37_{-.14}^{+.16}$ respectively. We show that the Lyman-alpha distribution of the 43 LAEs is consistent with being drawn from an exponential distribution with an e-folding scale of 150 Angstrom. Once complete the TESLA survey will enable the study of ~ thousands of LAEs to explore correlations between galaxy properties and Lyman-alpha EW. The large sample size will allow the construction of a predictive model for the Lyman-alpha EW as a function of SED-derived galaxy properties, which could be used to improve Lyman-alpha based constraints on reionization.

The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply fourteen of the most commonly used SED fitting codes on samples from the CANDELS photometric catalogs at z~1 and z~3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust, and AGN models) on the derived stellar masses, SFRs, and A_V values. We then assess the difference among the codes on the SFR-stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (~0.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting., Comment: 25 pages, 11 figures. Accepted for publication in ApJ

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 present the first comprehensive release of photometric redshifts (photo-z's) from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) team. We use statistics based upon the Quantile-Quantile (Q--Q) plot to identify biases and signatures of underestimated or overestimated errors in photo-z probability density functions (PDFs) produced by six groups in the collaboration; correcting for these effects makes the resulting PDFs better match the statistical definition of a PDF. After correcting each group's PDF, we explore three methods of combining the different groups' PDFs for a given object into a consensus curve. Two of these methods are based on identifying the minimum f-divergence curve, i.e., the PDF that is closest in aggregate to the other PDFs in a set (analogous to the median of an array of numbers). We demonstrate that these techniques yield improved results using sets of spectroscopic redshifts independent of those used to optimize PDF modifications. The best photo-z PDFs and point estimates are achieved with the minimum f-divergence using the best 4 PDFs for each object (mFDa4) and the Hierarchical Bayesian (HB4) methods, respectively. The HB4 photo-z point estimates produced $\sigma_{\rm NMAD} = 0.0227/0.0189$ and $|\Delta z/(1+z)| > 0.15$ outlier fraction = 0.067/0.019 for spectroscopic and 3D-HST redshifts, respectively. Finally, we describe the structure and provide guidance for the use of the CANDELS photo-z catalogs, which are available at https://archive.stsci.edu/hlsp/candels., Comment: 35 pages, 19 figures, published in ApJ, data available at https://archive.stsci.edu/hlsp/candels

Observations and simulations of interacting galaxies and mergers in the local universe have shown that interactions can significantly enhance the star formation rates (SFR) and fueling of Active Galactic Nuclei (AGN). However, at higher redshift, some simulations suggest that the level of star formation enhancement induced by interactions is lower due to the higher gas fractions and already increased SFRs in these galaxies. To test this, we measure the SFR enhancement in a total of 2351 (1327) massive ($M_*>10^{10}M_\odot$) major ($1

We present a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with a desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wavebands for measuring the physical properties of galaxies in a Hawaii Two-0 (H20)- and UVISTA-like survey for a sample of $i<25$ AB mag galaxies. We find that with available $i$-band fluxes, $r$, $u$, IRAC/$ch2$ and $z$ bands provide most of the information regarding the redshift with importance decreasing from $r$-band to $z$-band. We also find that for the same sample, IRAC/$ch2$, $Y$, $r$ and $u$ bands are the most relevant bands in stellar mass measurements with decreasing order of importance. Investigating the inter-correlation between the bands, we train a model to predict UVISTA observations in near-IR from H20-like observations. We find that magnitudes in $YJH$ bands can be simulated/predicted with an accuracy of $1\sigma$ mag scatter $\lesssim 0.2$ for galaxies brighter than 24 AB mag in near-IR bands. One should note that these conclusions depend on the selection criteria of the sample. For any new sample of galaxies with a different selection, these results should be remeasured. Our results suggest that in the presence of a limited number of bands, a machine learning model trained over the population of observed galaxies with extensive spectral coverage outperforms template-fitting. Such a machine learning model maximally comprises the information acquired over available extensive surveys and breaks degeneracies in the parameter space of template-fitting inevitable in the presence of a few bands., Comment: 15 pages, 14 figures, accepted for publication in ApJ

The $\rm{H}\alpha$-to-UV luminosity ratio ($L(\rm H\alpha)/L(\rm UV)$) is often used to probe SFHs of star-forming galaxies and it is important to validate it against other proxies for burstiness. To address this issue, we present a statistical analysis of the resolved distribution of $\Sigma_{\rm{SFR}}$ as well as stellar age and their correlations with the globally measured $L(\rm H\alpha)/L(\rm UV)$ for a sample of 310 star-forming galaxies in two redshift bins of $1.37 < z < 1.70$ and $ 2.09 < z < 2.61$ observed by the MOSDEF survey. We use the multi-waveband CANDELS/3D-HST imaging of MOSDEF galaxies to construct $\Sigma_{\rm{SFR}}$ and stellar age maps. We analyze the composite rest-frame far-UV spectra of a subsample of MOSDEF targets obtained by the Keck/LRIS, which includes 124 star-forming galaxies (MOSDEF-LRIS) at redshifts $1.4 < z < 2.6$, to examine the average stellar population properties, and the strength of age-sensitive FUV spectral features in bins of $L(\rm H\alpha)/L(\rm UV)$. Our results show no significant evidence that individual galaxies with higher $L(\rm H\alpha)/L(\rm UV)$ are undergoing a burst of star formation based on the resolved distribution of $\Sigma_{\rm{SFR}}$ of individual star-forming galaxies. We segregate the sample into subsets with low and high $L(\rm H\alpha)/L(\rm UV)$. The high-$L(\rm H\alpha)/L(\rm UV)$ subset exhibits, on average, an age of $\log[\rm{Age/yr}]$ = 8.0, compared to $\log[\rm{Age/yr}]$ = 8.4 for the low-$L(\rm H\alpha)/L(\rm UV)$ galaxies, though the difference in age is significant at only the $2\sigma$ level. Furthermore, we find no variation in the strengths of Siiv$\lambda\lambda1393, 1402$ and Civ$\lambda\lambda1548, 1550$ P-Cygni features from massive stars between the two subsamples., Comment: 16 pages, 10 figures, published by the Monthly Notices of the Royal Astronomical Society

We explore how the fraction of quenched galaxies changes in groups of galaxies with respect to the distance to the center of the group, redshift, and stellar mass to determine the dominant process of environmental quenching in $0.2 < z < 0.8$ groups. We use new UV data from the UVCANDELS project in addition to existing multiband photometry to derive new galaxy physical properties of the group galaxies from the zCOSMOS 20k Group Catalog. Limiting our analysis to a complete sample of log$(M_*/M_{\odot})>10.56$ group galaxies we find that the probability of being quenched increases slowly with decreasing redshift, diverging from the stagnant field galaxy population. A corresponding analysis on how the probability of being quenched increases with time within groups suggests that the dominant environmental quenching process is characterized by slow ($\sim$Gyr) timescales. We find a quenching time of approximately $4.91^{+0.91}_{-1.47} $Gyrs, consistent with the slow processes of strangulation (Larson et al. 1980) and delayed-then-rapid quenching (Wetzel et al. 2013 arXiv:1206.3571v2 [astro-ph.CO]), although more data are needed to confirm this result.

Despite advances in observational data, theoretical models, and computational techniques to simulate key physical processes in the formation and evolution of galaxies, the stellar mass assembly of galaxies still remains an unsolved problem today. Optical spectroscopic measurements appear to show that the gas-phase metallicities of local ultra-luminous infrared galaxies (ULIRGs) are significantly lower than those of normal star-forming galaxies. This difference has resulted in the claim that ULIRGs are fueled by metal-poor gas accretion from the outskirts\cite{Mannucci10}. Here we report on a new set of gas-phase metallicity measurements making use of the far-infrared spectral lines of [O{\sc iii}]52 $\mu$m, [O{\sc iii}]88 $\mu$m, and [N{\sc iii}]57 $\mu$m instead of the usual optical lines. Photoionization models have resulted in a metallicity diagnostic based on these three lines that break the electron density degeneracy and reduce the scatter of the correlation significantly. Using new data from SOFIA and archival data from Herschel Space Observatory, we find that local ULIRGs lie on the mass-metallicity relation of star-forming galaxies and have metallicities comparable to other galaxies with similar stellar masses and star formation rates. The lack of a departure suggests that ULIRGs follow the same mass assembly mechanism as luminous star-forming galaxies and $\sim 0.3$ dex under-abundance in metallicities derived from optical lines is a result of heavily obscured metal-rich gas which has a negligible effect when using the FIR line diagnostics., Comment: Preprint submitted for publication 10/29/21

We study the mass-metallicity relation for 19 members of a spectroscopically-confirmed protocluster in the COSMOS field at $z=2.2$ (CC2.2), and compare it with that of 24 similarly selected field galaxies at the same redshift. Both samples are $\rm H\alpha$ emitting sources, chosen from the HiZELS narrow-band survey, with metallicities derived from $\rm N2\ (\frac{\rm [NII] \lambda 6584}{\rm H \alpha})$ line ratio. For the mass-matched samples of protocluster and field galaxies, we find that protocluster galaxies with $10^{9.9} \rm M_\odot \leq M_* \leq 10^{10.9} \rm M_\odot$ are metal deficient by $0.10 \pm 0.04$ dex ($2.5\sigma$ significance) compared to their coeval field galaxies. This metal deficiency is absent for low mass galaxies, $\rm M_* < 10^{9.9} \rm M_\odot$. Moreover, relying on both SED-derived and $\rm {H\alpha}$ (corrected for dust extinction based on $\rm {M_*}$) SFRs, we find no strong environmental dependence of SFR-$\rm {M_*}$ relation, however, we are not able to rule out the existence of small dependence due to inherent uncertainties in both SFR estimators. The existence of $2.5\sigma$ significant metal deficiency for massive protocluster galaxies favors a model in which funneling of the primordial cold gas through filaments dilutes the metal content of protoclusters at high redshifts ($z \gtrsim 2$). At these redshifts, gas reservoirs in filaments are dense enough to cool down rapidly and fall into the potential well of the protocluster to lower the gas-phase metallicity of galaxies. Moreover, part of this metal deficiency could be originated from galaxy interactions which are more prevalent in dense environments., Comment: 10 pages, 5 figures, 2 tables, accepted for publication in ApJ

Using the near-IR spectroscopy of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the role of local environment in the gas-phase metallicity of galaxies. The local environment measurements are derived from accurate and uniformly calculated photometric redshifts with well-calibrated probability distributions. Based on rest-frame optical emission lines, [NII]$\lambda6584$ and H$\alpha$, we measure gas-phase oxygen abundance of 167 galaxies at $1.37\leq z\leq1.7$ and 303 galaxies at $2.09\leq z\leq2.61$, located in diverse environments. We find that at $z\sim1.5$, the average metallicity of galaxies in overdensities with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.8}M_\odot$ is higher relative to their field counterparts by $0.094\pm0.051$, $0.068\pm0.028$ and $0.052\pm0.043$ dex, respectively. However, this metallicity enhancement does not exist at higher redshift, $z\sim2.3$, where, compared to the field galaxies, we find $0.056\pm0.043$, $0.056\pm0.028$ and $0.096\pm 0.034$ dex lower metallicity for galaxies in overdense environments with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.7}M_\odot$, respectively. Our results suggest that, at $1.37\leq z\leq2.61$, the variation of mass-metallicity relation with local environment is small ($<0.1$dex), and reverses at $z\sim2$. Our results support the hypothesis that, at the early stages of cluster formation, owing to efficient gas cooling, galaxies residing in overdensities host a higher fraction of pristine gas with prominent primordial gas accretion, which lowers their gas-phase metallicity compared to their coeval field galaxies. However, as the Universe evolves to lower redshifts ($z\lesssim2$), shock-heated gas in overdensities cannot cool down efficiently, and galaxies become metal-rich rapidly due to the suppression of pristine gas inflow and re-accretion of metal-enriched outflows in overdensities., Comment: 18 pages, 8 figures, 3 tables, accepted for publication in ApJ

We use a sample of star-forming field and protocluster galaxies at z=2.0-2.5 with Keck/MOSFIRE K-band spectra, a wealth of rest-frame UV photometry, and Spitzer/MIPS and Herschel/PACS observations, to dissect the relation between the ratio of IR to UV luminosity (IRX) versus UV slope ($\beta$) as a function of gas-phase metallicity (12+log(O/H)~8.2-8.7). We find no significant dependence of the IRX-$\beta$ trend on environment. However, we find that at a given $\beta$, IRX is highly correlated with metallicity, and less correlated with mass, age, and sSFR. We conclude that, of the physical properties tested here, metallicity is the primary physical cause of the IRX-$\beta$ scatter, and the IRX correlation with mass is presumably due to the mass dependence on metallicity. Our results indicate that the UV attenuation curve steepens with decreasing metallicity, and spans the full range of slope possibilities from a shallow Calzetti-type curve for galaxies with the highest metallicity in our sample (12+log(O/H)~8.6) to a steep SMC-like curve for those with 12+log(O/H)~8.3. Using a Calzetti (SMC) curve for the low (high) metallicity galaxies can lead to up to a factor of 3 overestimation (underestimation) of the UV attenuation and obscured SFR. We speculate that this change is due to different properties of dust grains present in the ISM of low- and high-metallicity galaxies., Comment: Accepted for publication in ApJL

Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the largest known sample of major spectroscopic galaxy pairs (2381 with $\Delta V <5000$ km s$^{-1}$) at $0.5

Using the CANDELS photometric catalogs for the HST/ACS and WFC3, we identified massive evolved galaxies at $3 < z < 4.5$, employing three different selection methods. We find the comoving number density of these objects to be $\sim 2 \times 10^{-5}$ and $8 \times 10^{-6}Mpc^{-3}$ after correction for completeness for two redshift bins centered at $z=3.4, 4.7$. We quantify a measure of how much confidence we should have for each candidate galaxy from different selections and what are the conservative error estimates propagated into our selection. Then we compare the evolution of the corresponding number densities and their stellar mass density with numerical simulations, semi-analytical models, and previous observational estimates, which shows slight tension at higher redshifts as the models tend to underestimate the number and mass densities. By estimating the average halo masses of the candidates ($M_h \approx 4.2, 1.9, 1.3 \times 10^{12} M_\odot$ for redshift bins centered at $z=3.4, 4.1, 4.7$), we find them to be consistent with halos that were efficient in turning baryons to stars and were relatively immune to the feedback effects and on the verge of transition into hot-mode accretion. This can suggest the relative cosmological starvation of the cold gas followed by an overconsumption phase in which the galaxy consumes the available cold gas rapidly as one of the possible drivers for the quenching of the massive evolved population at high redshift., Comment: 20 pages, 11 figures, Accepted for publication in ApJ

We report the spectroscopic confirmation of a new protocluster in the COSMOS field at $z$ $\sim$ 2.2, COSMOS Cluster 2.2 (CC2.2), originally identified as an overdensity of narrowband selected H$\alpha$ emitting candidates. With only two masks of Keck/MOSFIRE near-IR spectroscopy in both $H$ ($\sim$ 1.47-1.81 $\mu$m) and $K$ ($\sim$ 1.92-2.40 $\mu$m) bands ($\sim$ 1.5 hour each), we confirm 35 unique protocluster members with at least two emission lines detected with S/N $>$ 3. Combined with 12 extra members from the zCOSMOS-deep spectroscopic survey (47 in total), we estimate a mean redshift and a line-of-sight velocity dispersion of $z_{mean}$=2.23224 $\pm$ 0.00101 and $\sigma_{los}$=645 $\pm$ 69 km s$^{-1}$ for this protocluster, respectively. Assuming virialization and spherical symmetry for the system, we estimate a total mass of $M_{vir}$ $\sim$ $(1-2) \times$10$^{14}$ $M_{\odot}$ for the structure. We evaluate a number density enhancement of $\delta_{g}$ $\sim$ 7 for this system and we argue that the structure is likely not fully virialized at $z$ $\sim$ 2.2. However, in a spherical collapse model, $\delta_{g}$ is expected to grow to a linear matter enhancement of $\sim$ 1.9 by $z$=0, exceeding the collapse threshold of 1.69, and leading to a fully collapsed and virialized Coma-type structure with a total mass of $M_{dyn}$($z$=0) $\sim$ 9.2$\times$10$^{14}$ $M_{\odot}$ by now. This observationally efficient confirmation suggests that large narrowband emission-line galaxy surveys, when combined with ancillary photometric data, can be used to effectively trace the large-scale structure and protoclusters at a time when they are mostly dominated by star-forming galaxies.

The position of galaxies on the stellar mass, star formation rate plane with respect to the star-forming main sequence at each redshift is a convenient way to infer where the galaxy is in its evolution compared to the rest of the population. We use Hubble Space Telescope high resolution images in the GOODS-S field from the CANDELS survey and fit multi wavelength lights in resolution elements of galaxies with stellar population synthesis models. We then construct resolved kpc-scale stellar mass, star formation rate surface density curves for galaxies at z~1. Fitting these resolved main sequence curves with Schechter functions, we parameterize and explain the multi-wavelength structure of galaxies with three variables: phi*, alpha, and M*. For quenched galaxies below the main sequence, we find an average high mass slope (alpha) of the resolved main sequence curves to be ~ -0.4. The scatter of this slope is higher among the lower mass star forming galaxies and those above the main sequence compared to quenched galaxies, due to lack of an evolved bulge. Our findings agree well with an inside-out quenching of star-formation. We find that the knee of the Schechter fits (M*) for galaxies below the main sequence occurs at lower stellar mass surface densities compared to star forming galaxies, which hints at how far quenching has proceeded outwards., Comment: 5 pages, 3 figures, Accepted for publication by ApJL

We present a robust method, weighted von Mises kernel density estimation, along with boundary correction to reconstruct the underlying number density field of galaxies. We apply this method to galaxies brighter than $\rm HST/F160w\le 26$ AB mag at the redshift range of $0.4\leq z \leq 5$ in the five CANDELS fields (GOODS-N, GOODS-S, EGS, UDS, and COSMOS). We then use these measurements to explore the environmental dependence of the star formation activity of galaxies. We find strong evidence of environmental quenching for massive galaxies ($\rm M \gtrsim 10^{11} \rm {M}_\odot$) out to $z\sim 3.5$ such that an over-dense environment hosts $\gtrsim 20\%$ more massive quiescent galaxies compared to an under-dense region. We also find that environmental quenching efficiency grows with stellar mass and reaches $\sim 60\%$ for massive galaxies at $z\sim 0.5$. The environmental quenching is also more efficient in comparison to the stellar mass quenching for low mass galaxies ($\rm M \lesssim 10^{10} \rm {M}_\odot$) at low and intermediate redshifts ($z\lesssim 1.2$). Our findings concur thoroughly with the "over-consumption" quenching model where the termination of cool gas accretion (cosmological starvation) happens in an over-dense environment and the galaxy starts to consume its remaining gas reservoir in depletion time. The depletion time depends on the stellar mass and could explain the evolution of environmental quenching efficiency with the stellar mass., Comment: 17 pages, 10 figures, 3 tables, accepted for publication in the ApJ. Catalogs of the environment measurements and the density maps will be available with the final published paper in ApJ

Studying the internal processes of individual galaxies at kilo-parsec scales is crucial in enhancing our understanding of galaxy formation and evolution processes. In this work, we investigate the distribution of star formation rate (SFR), specific SFR (sSFR), and dust attenuation across individual galaxies for a sample of 32 galaxies selected from the MUSE-Wide Survey at 0.1 $< \textit{z} <$ 0.42 with a dynamic range in stellar masses between $10^{7.7}$ and $10^{10.3}$ M$_{\odot}$. We take advantage of the high spatial resolution of the MUSE integral-field spectrograph and measure reliable spatially resolved H$\alpha$ and H$\beta$ emission line maps for individual galaxies. We also derive resolved stellar mass, SFR and dust maps using pixel-by-pixel SED fitting on high resolution multi-band HST/ACS and HST/WFC3 data from the CANDELS survey. By combining these, we analyze the radial profile of various physical parameters across these galaxies. We observe a radial dependence in both stellar and nebular color excess profiles peaking at the inner regions of galaxies. We also find the color excess profiles to most strongly correlate with the integrated sSFRs of galaxies. The median sSFR$_{\mathrm{H}\alpha}$ radial profiles of galaxies in our sample show a 0.8 dex increase from the central regions outward. This increase compared to the almost flat median radial profile of sSFR$_{\mathrm{SED}}$, which traces longer timescales of star formation, is in favor of the inside-out quenching of star formation. We bring further evidence for this quenching scenario from the locus of different subregions of galaxies on the SFR-M$_{*}$ and sSFR-M$_{*}$ relations., Comment: 9 pages, 9 figures, accepted for publication in The Astrophysical Journal

With current and upcoming experiments such as WFIRST, Euclid and LSST, we can observe up to billions of galaxies. While such surveys cannot obtain spectra for all observed galaxies, they produce galaxy magnitudes in color filters. This data set behaves like a high-dimensional nonlinear surface, an excellent target for machine learning. In this work, we use a lightcone of semianalytic galaxies tuned to match CANDELS observations from Lu et al. (2014) to train a set of neural networks on a set of galaxy physical properties. We add realistic photometric noise and use trained neural networks to predict stellar masses and average star formation rates on real CANDELS galaxies, comparing our predictions to SED fitting results. On semianalytic galaxies, we are nearly competitive with template-fitting methods, with biases of $0.01$ dex for stellar mass, $0.09$ dex for star formation rate, and $0.04$ dex for metallicity. For the observed CANDELS data, our results are consistent with template fits on the same data at $0.15$ dex bias in $M_{\rm star}$ and $0.61$ dex bias in star formation rate. Some of the bias is driven by SED-fitting limitations, rather than limitations on the training set, and some is intrinsic to the neural network method. Further errors are likely caused by differences in noise properties between the semianalytic catalogs and data. Our results show that galaxy physical properties can in principle be measured with neural networks at a competitive degree of accuracy and precision to template-fitting methods., Comment: 19 pages, 10 figures, 6 tables. Accepted for publication in ApJ