Luminosity Function Research Articles

A Steep Decline in the Galaxy Space Density beyond Redshift 9 in the CANUCS UV Luminosity Function

We present a new sample of 158 galaxies at redshift z > 7.5 selected from deep James Webb Space Telescope (JWST) NIRCam imaging of five widely separated sight lines in the CANUCS survey. Two-thirds of the pointings and 80% of the galaxies are covered by 12–14 NIRCam filters, including seven to nine medium bands, providing accurate photometric redshifts and robustness against low-redshift interlopers. A sample of 28 galaxies at z > 7.5 with spectroscopic redshifts shows a low systematic offset and scatter in the difference between photometric and spectroscopic redshifts. We derive the galaxy UV luminosity function at redshifts 8–12, finding a slightly higher normalization than previously seen with the Hubble Space Telescope at redshifts 8–10. We observe a steeper decline in the galaxy space density from z = 8 to 12 than found by most JWST Cycle 1 studies. In particular, we find only eight galaxies at z > 10 and none at z > 12.5, with no z > 10 galaxies brighter than F277W AB = 28 or M UV = −20 in our unmasked, delensed survey area of 53.4 arcmin2. We attribute the lack of bright z > 10 galaxies in CANUCS compared to GLASS and CEERS to intrinsic variance in the galaxy density along different sight lines. The evolution in the CANUCS luminosity function between z = 8 and 12 is comparable to that predicted by simulations that assume a standard star formation efficiency without invoking any special adjustments.

SPICE: the connection between cosmic reionization and stellar feedback in the first galaxies

ABSTRACT We present SPICE, a new suite of radiation-hydrodynamic, cosmological simulations targeting the epoch of reionization. The goal of these simulations is to systematically probe a variety of stellar feedback models, including ‘bursty’ and ‘smooth’ forms of supernova energy injection, as well as poorly explored physical scenarios such as hypernova explosions and radiation pressure on dust. We show that even subtle differences in the behaviour of supernova feedback drive profound differences in reionization histories, with burstier forms of feedback causing earlier reionization. However, we also find that some global galaxy properties, such as the dust-attenuated luminosity functions and star formation main sequence, remain degenerate between models. In particular, we show that stellar feedback and its strength determine the morphological mix of galaxies emerging by $z \, = \, 5$ and that the reionization history is inextricably connected to intrinsic properties such as galaxy kinematics and morphology. While star-forming, massive discs are prevalent if supernova feedback is ‘smooth’, ‘bursty’ feedback preferentially generates dispersion-dominated systems. Different modes of feedback produce different strengths of outflows, altering the interstellar/circumgalactic medium in different ways, and in turn strongly affecting the escape of Lyman continuum (LyC) photons. We establish a correlation between galaxy morphology and LyC escape fraction, revealing that dispersion-dominated systems have escape fractions 10–50 times higher than their rotation-dominated counterparts at all redshifts. At the same intrinsic luminosity, dispersion-dominated systems should thus preferentially generate large H ii regions as compared to their rotation-dominated counterparts. Since dispersion-dominated systems are more prevalent if stellar feedback is more explosive, reionization occurs earlier in our simulation with burstier feedback. We argue that statistical samples of post-reionization galaxy morphologies (using both stellar and gaseous components) probed with telescopes such as JWST, ALMA, and MUSE can constrain stellar feedback at z > 5 and models of cosmic reionization.

Reconnaissance with JWST of the J-region Asymptotic Giant Branch in Distance Ladder Galaxies: From Irregular Luminosity Functions to Approximation of the Hubble Constant

We study stars in the J-regions of the asymptotic giant branch (JAGB) of near-infrared color–magnitude diagrams in the maser host NGC 4258 and four hosts of six Type Ia supernovae (SNe Ia): NGC 1448, NGC 1559, NGC 5584, and NGC 5643. These clumps of stars are readily apparent near 1.0 < F150W − F277W < 1.5 and m F150W = 22–25 mag with James Webb Space Telescope NIRCam photometry. Various methods have been proposed to assign an apparent reference magnitude to this recently proposed standard candle, including the mode, median, sigma-clipped mean, or a modeled luminosity function parameter. We test the consistency of these by measuring intrahost variations, finding differences of up to ∼0.2 mag that significantly exceed statistical uncertainties. Brightness differences appear intrinsic, and are further amplified by the nonuniform shape of the JAGB luminosity function, also apparent in the LMC and SMC. We follow a “many methods” approach to measure consistently JAGB magnitudes and distance moduli to the SN Ia host sample calibrated by NGC 4258. We find broad agreement with distance moduli measured from Cepheids, tip of the red giant branch, and Miras. However, the SN host mean distance modulus estimated via the JAGB method necessary to estimate H 0 differs by ∼0.19 mag among the above definitions, the result of different levels of luminosity function asymmetry. The methods yield a full range of 71−78 km s−1 Mpc−1, i.e., a fiducial result of H 0 = 74.7 ± 2.1(stat) ± 2.3(sys, ±3.1 if combined in quadrature) km s−1 Mpc−1, with systematic errors limited by the differences in methods. Future work may seek to standardize and refine this promising tool further, making it more competitive with established distance indicators.

Compatibility of JWST results with exotic halos

The James Webb Space Telescope (JWST) is unveiling astounding results about the first few hundred million years of life of the Universe, delivering images of galaxies at very high redshifts. Here, we develop a UV luminosity function model for high-redshift galaxies, considering parameters such as the stellar formation rate, dust extinction, and halo mass function. Calibration of this luminosity function model using UV luminosity data at redshifts z=4-7 yields optimal parameter values. Testing the model against data at higher redshifts reveals successful accommodation of the data at z=8-9, but challenges emerge at z≃13. Our findings suggest a negligible role of dust extinction at the highest redshifts, prompting a modification of the stellar formation rate to incorporate a larger fraction of luminous objects per massive halo, consistently with similar recent studies. This effect could be attributed to mundane explanations such as unknown evolution of standard astrophysics at high redshift or to the existence of exotic objects at high redshift. We comment on this latter possibility.

Euclid view of the dusty star-forming galaxies at ≳ detected in wide area submillimetre surveys

ABSTRACT We investigate the constraints provided by the Euclid space observatory on the physical properties of dusty star-forming galaxies (DSFGs) at $z\gtrsim 1.5$ detected in wide area submillimetre surveys with Herschel. We adopt a physical model for the high-z progenitors of spheroidal galaxies, which form the bulk of DSFGs at $z\gtrsim 1.5$. We improve the model by combining the output of the equations of the model with a formalism for the spectral energy distribution (SED). After optimizing the SED parameters to reproduce the measured infrared luminosity function and number counts of DFSGs, we simulated a sample of DSFGs over 100 $\hbox{deg}^2$ and then applied a $5\, \sigma$ detection limit of $37\,$ mJy at 250 $\mu$m. We estimated the redshifts from the Euclid data and then fitted the Euclid+Herschel photometry with the code CIGALE to extract the physical parameters. We found that 100 per cent of the Herschel galaxies are detected in all 4 Euclid bands above $3\, \sigma$. For 87 per cent of these sources the accuracy on $1 + z$ is better than 15 per cent. The sample comprises mostly massive, i.e. $\log (M_{\star }/{\rm M}_{\odot })\sim 10.5{\!-\!}12.9$, highly star forming, i.e. $\log (\hbox{SFR}/{\rm M}_{\odot }\hbox{yr}^{-1})\sim 1.5{\!-\!}4$, dusty, i.e. $\log (M_{\rm dust}/{\rm M}_{\odot })\sim 7.5{\!-\!}9.9$, galaxies. The measured stellar masses have a dispersion of 0.19 dex around the true value, thus showing that Euclid will provide reliable stellar mass estimates for the majority of the bright DSFGs at $z\gtrsim 1.5$ detected by Herschel. We also explored the effect of complementing the Euclid photometry with that from the Vera C. Rubin Observatory/LSST.

NSCs from groups to clusters: A catalogue of dwarf galaxies in the Shapley Supercluster and the role of environment in galaxy nucleation

Abstract Nuclear star clusters (NSCs) are dense star clusters located at the centre of galaxies spanning a wide range of masses and morphologies. Analysing NSC occupation statistics in different environments provides an invaluable window into investigating early conditions of high-density star formation and mass assembly in clusters and group galaxies. We use HST/ACS deep imaging to obtain a catalogue of dwarf galaxies in two galaxy clusters in the Shapley Supercluster: the central cluster Abell 3558 and the northern Abell 1736a. The Shapley region is an ideal laboratory to study nucleation as it stands as the highest mass concentration in the nearby Universe. We investigate the NSC occurrence in quiescent dwarf galaxies as faint as MI = −10 mag and compare it with all other environments where nucleation data is available. We use galaxy cluster/group halo mass as a proxy for the environment and employ a Bayesian logistic regression framework to model the nucleation fraction (fn) as a function of galaxy luminosity and environment. We find a notably high fn in Abell 3558: at MI ≈ −13.1 mag, half the galaxies in the cluster host NSCs. This is higher than in the Virgo and Fornax clusters but comparable to the Coma Cluster. On the other hand, the fn in Abell 1736a is relatively lower, comparable to groups in the Local Volume. We find that the probability of nucleation varies with galaxy luminosity remarkably similarly in galaxy clusters. These results reinforce previous findings of the important role of the environment in NSC formation/growth.

The High-mass X-Ray Binary Luminosity Functions of Dwarf Galaxies

Drawing from the Chandra archive and using a carefully selected set of nearby dwarf galaxies, we present a calibrated high-mass X-ray binary (HMXB) luminosity function in the low-mass galaxy regime and search for an already hinted at dependence on metallicity. Our study introduces a new sample of local dwarf galaxies (D < 12.5 Mpc and M * < 5 × 109 M ⊙), expanding the specific star formation rates (sSFR) and gas-phase metallicities probed in previous investigations. Our analysis of the observed X-ray luminosity function indicates a shallower power-law slope for the dwarf galaxy HMXB population. In our study, we focus on dwarf galaxies that are more representative in terms of sSFR compared to prior work. In this regime, the HMXB luminosity function exhibits significant stochastic sampling at high luminosities. This likely accounts for the pronounced scatter observed in the galaxy-integrated HMXB population’s L X/SFR versus metallicity for our galaxy sample. Our calibration is necessary to understand the active galactic nuclei content of low-mass galaxies identified in current and future X-ray survey fields and has implications for binary population synthesis models, as well as X-ray-driven cosmic heating in the early Universe.

ALMA survey of a massive node of the Cosmic Web at z ∼ 3

Submillimeter surveys toward overdense regions in the early Universe are essential for uncovering the obscured star formation and the cold gas content of assembling galaxies within massive dark matter halos. In this work, we present deep ALMA mosaic observations covering an area of ∼2′×2′ around MUSE Quasar Nebula 01 (MQN01), one of the largest and brightest Ly-α emitting nebulae discovered thus far; it surrounds a radio-quiet quasar at z ≃ 3.25. Our observations target the 1.2 and the 3 mm dust continuum as well as the carbon monoxide CO(4–3) transition in galaxies in the vicinity of the quasar. We identify a robust sample of 11 CO-line-emitting galaxies (including a closely separated quasar companion) that lie within ±4000 km s−1 of the quasar systemic redshift. A fraction of these objects were missed in previous deep rest-frame optical/UV surveys, which highlights the critical role of (sub)millimeter imaging. We also detect a total of 11 sources revealed in the dust continuum at 1.2 mm; six of them have either high-fidelity spectroscopic redshift information from rest-frame UV metal absorptions or the CO(4–3) line that places them in the same narrow redshift range. A comparison of the CO luminosity function and 1.2 mm number count density with those of the general fields points to a galaxy overdensity of δ &gt; 10. We find evidence of a systematic flattening at the bright end of the CO luminosity function with respect to the trend measured in blank fields. Our findings reveal that galaxies in dense regions at z ∼ 3 are more massive and significantly richer in molecular gas than galaxies in fields, which enables a faster and accelerated assembly. This is the first in a series of studies aimed at characterizing one of the densest regions of the Universe found so far at z &gt; 3.

A Candidate Supermassive Black Hole in a Gravitationally Lensed Galaxy at Z ≈ 10

While supermassive black holes (SMBHs) are widely observed in the nearby and distant Universe, their origin remains debated with two viable formation scenarios with light and heavy seeds. In the light seeding model, the seed of the first SMBHs form from the collapse of massive stars with masses of 10–100 M ⊙, while the heavy seeding model posits the formation of 104–5 M ⊙ seeds from direct collapse. The detection of SMBHs at redshifts z ≳ 10, edging closer to their formation epoch, provides critical observational discrimination between these scenarios. Here, we focus on the JWST-detected galaxy, GHZ 9, at z ≈ 10 that is lensed by the foreground cluster, A2744. Based on 2.1 Ms deep Chandra observations, we detect a candidate X-ray active galactic nucleus (AGN), which is spatially coincident with the high-redshift galaxy, GHZ 9. The SMBH candidate is inferred to have a bolometric luminosity of (1.0−0.4+0.5)×1046ergs−1 , which corresponds to a black hole (BH) mass of (8.0−3.2+3.7)×107M⊙ assuming Eddington-limited accretion. This extreme mass at such an early cosmic epoch suggests the heavy seed origin for this BH candidate. Based on the Chandra and JWST discoveries of extremely high-redshift quasars, we have constructed the first simple AGN luminosity function extending to z ≈ 10. Comparison of this luminosity function with theoretical models indicates an overabundant z ≈ 10 SMBH population, consistent with a higher-than-expected seed formation efficiency.

The universal power spectrum of quasars in optical wavelengths

Context. The optical variability of quasars is one of the few windows through which we can explore the behaviour of accretion discs around supermassive black holes. Aims. We aim to establish the dependence of variability properties, such as characteristic timescales and the variability amplitude, on basic quasar parameters such as black hole mass and the accretion rate, controlling for the rest-frame wavelength of emission. Methods. Using large catalogues of quasars, we selected the g-band light curves for 4770 objects from the Zwicky Transient Facility archive. All the selected objects fall into a narrow redshift bin, 0.6 &lt; z &lt; 0.7, but cover a wide range of accretion rates in Eddington units (REdd) and black hole masses (M). We grouped these objects into 26 independent bins according to these parameters, calculated low-resolution g-band variability power spectra for each of these bins, and approximated the power spectra with a simple analytic model that features a break at a timescale, tb. Results. We find a clear dependence of the break timescale, tb, on REdd, on top of the known dependence of tb on the black hole mass, M. In our fits, tb ∝ M0.65 − 0.55REdd0.35−0.3, where the ranges in the exponents correspond to the best-fitting parameters of different power spectrum models. This mass dependence is slightly steeper than that found in other studies. Scaling tb to the orbital timescale of the innermost stable circular orbit (ISCO), tISCO, results approximately in tb/tISCO ∝ (REdd/M)0.35. In the standard thin disc model, (REdd/M) ∝ Tmax4, where Tmax is the maximum disc temperature, so that tb/tISCO appears to scale approximately with the maximum temperature of the disc to a small power. The observed values of tb are ∼10 longer than the orbital timescale at the light-weighted average radius of the disc region emitting in the (observer frame) g-band. The different scaling of the break frequency with M and REdd shows that the shape of the variability power spectrum cannot be solely a function of the quasar luminosity, even for a single rest-frame wavelength. Finally, the best-fitting models have slopes above the break in the range between −2.5 and −3. A slope of −2, as in the damped random walk models, fits the data significantly worse.

An X-Ray Census of Active Galactic Nuclei in the Virgo and Fornax Clusters of Galaxies with SRG/eROSITA

We present a uniform and sensitive X-ray census of active galactic nuclei (AGNs) in the two nearest galaxy clusters, Virgo and Fornax, utilizing the newly released X-ray source catalogs from the first all-sky scan of Spectrum-Roentgen-Gamma/eROSITA. A total of 50 and 10 X-ray sources are found positionally coincident with the nuclei of member galaxies in Virgo and Fornax, respectively, down to a 0.2–2.3 keV luminosity of ∼1039 erg s−1 and reaching out to a projected distance well beyond the virial radius of both clusters. The majority of the nuclear X-ray sources are newly identified. There is weak evidence that the nuclear X-ray sources are preferentially found in late-type hosts. Several hosts are dwarf galaxies with a stellar mass below ∼109 M ⊙. We find that contamination by nonnuclear X-ray emission can be neglected in most cases, indicating the dominance of a genuine AGN. In the meantime, no nuclear X-ray source exhibits a luminosity higher than a few times 1041 erg s−1, which might be owing to a steep intrinsic luminosity function. The X-ray AGN occupation rate is only ∼3% in both clusters, apparently much lower than that in field galaxies inferred from previous X-ray studies. Both aspects suggest that the cluster environment effectively suppresses AGN activity. The findings of this census have important implications for the interplay between galaxies and their central massive black holes in cluster environments.

A too-many-dwarf-galaxy-satellites problem in the M 83 group

Dwarf galaxies in groups of galaxies provide excellent test cases for models of structure formation. This led to a so-called small-scale crisis, including the famous missing-satellites and too-big-to-fail problems. It was suggested that these two problems can be resolved by introducing baryonic physics to cosmological simulations. We tested the nearby grand spiral M 83 – a Milky Way sibling – to determine whether its number of dwarf galaxy companions is compatible with today’s Λ cold dark matter model using two methods: with cosmological simulations that include baryons and with theoretical predictions from the subhalo mass function. By employing distance measurements, we recovered a list of confirmed dwarf galaxies within 330 kpc of M 83 down to a magnitude of MV = −10. We find that both the state-of-the-art hydrodynamical cosmological simulation Illustris-TNG50 and theoretical predictions agree with the number of confirmed satellites around M 83 at the bright end of the luminosity function (&gt; 108 solar masses) but underestimate it at the faint end (down to 106 solar masses) at more than 3σ and 5σ levels, respectively. This indicates a too-many-satellites problem for M 83 in the Λ cold dark matter model. The actual degree of tension with cosmological models is underestimated because the number of observed satellites is incomplete due to the high contamination of spurious stars and Galactic cirrus.

The neutrino background from non-jetted active galactic nuclei

Aims. We calculate the contribution to the neutrino background from the non-jetted active galactic nuclei (AGN) population following the recent IceCube association of TeV neutrinos with NGC 1068. Methods. We exploited our robust knowledge of the AGN X-ray luminosity function and evolution and converted it to the neutrino band by using NGC 1068 as a benchmark, together with a theoretically motivated neutrino spectrum. Results. The resulting neutrino background up to redshift 5 does not violate either the IceCube diffuse flux or the upper bounds for non-jetted AGN, although barely so. This is consistent with a scenario in which the latter class makes a substantial contribution mostly below 1 PeV, while jetted AGN, that is, blazars, dominate above this energy, in intriguing agreement with the dip in the neutrino data at ∼300 TeV. More and better IceCube data on Seyfert galaxies will allow us to constrain the fraction of neutrino emitters among non-jetted AGN.

Fast forward modelling of galaxy spatial and statistical distributions

A forward modelling approach provides simple, fast and realistic simulations of galaxy surveys, without a complex underlying model. For this purpose, galaxy clustering needs to be simulated accurately, both for the usage of clustering as its own probe and to control systematics. We present a forward model to simulate galaxy surveys, where we extend the Ultra-Fast Image Generator to include galaxy clustering. We use the distribution functions of the galaxy properties, derived from a forward model adjusted to observations. This population model jointly describes the luminosity functions, sizes, ellipticities, SEDs and apparent magnitudes. To simulate the positions of galaxies, we then use a two-parameter relation between galaxies and halos with Subhalo Abundance Matching (SHAM). We simulate the halos and subhalos using the fast PINOCCHIO code, and a method to extract the surviving subhalos from the merger history. Our simulations contain a red and a blue galaxy population, for which we build a SHAM model based on star formation quenching. For central galaxies, mass quenching is controlled with the parameter Mlimit, with blue galaxies residing in smaller halos. For satellite galaxies, environmental quenching is implemented with the parameter tquench, where blue galaxies occupy only recently merged subhalos. We build and test our model by comparing to imaging data from the Dark Energy Survey Year 1. To ensure completeness in our simulations, we consider the brightest galaxies with i < 20. We find statistical agreement between our simulations and the data for two-point correlation functions on medium to large scales. Our model provides constraints on the two SHAM parameters Mlimit and tquench and offers great prospects for the quick generation of galaxy mock catalogues, optimized to agree with observations.

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy ( flim,line,5σ≈1.2× 10−18 erg s−1 cm−2) to measure metallicities and star formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.5 < z < 4). In parallel, NGDEEP targets the HUDF-Par2 parallel field with NIRCam ( mlim,5σ=30.6−30.9 ) to discover galaxies to z > 12, constraining the slope of the faint end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST Advanced Camera for Surveys optical imaging in the sky, F435W in the HUDF ( mlim,F435W=29.6 ) and F814W in HUDF-Par2 ( mlim,F814W=30 ), making this a premier HST+JWST deep field. As a treasury survey, NGDEEP data are public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep-field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products.

A JWST Project on 47 Tucanae. Overview, Photometry, and Early Spectroscopic Results of M Dwarfs and Observations of Brown Dwarfs* *The main data set presented in this article was obtained from the Mikulski Archive for Space Telescopes (MAST) at the Space Telescope Science Institute. The specific observations analyzed can be accessed via doi:10.17909/6t82-4360.

James Webb Space Telescope (JWST) observations have been demonstrated to be efficient in detecting multiple stellar populations in globular clusters (GCs) in the low-mass regime of M dwarfs. We present an overview, and first results, of different projects that can be explored by using the JWST observations gathered under program GO2560 for 47 Tucanae, the first program entirely devoted to the investigation of multiple populations in very-low-mass stars, which includes spectroscopic data for the faintest GC stars for which spectra are available. Our color–magnitude diagram (CMD) shows some substructures for ultracool stars, including gaps and breaks in slope. In particular, we observe both a gap and a minimum in the F322W2 luminosity function less than 1 mag apart, and discuss which it could be associated with the H-burning limit. We detect stars fainter than this minimum, very likely brown dwarfs. We corroborate the ubiquity of the multiple populations across different masses, from ∼0.1 M ⊙ up to red giants (∼0.8 M ⊙). The oxygen range inferred for the M dwarfs, both from the CMD and from the spectra of two M dwarfs associated with different populations, is similar to that observed for giants. We have not detected any difference between the fractions of stars in distinct populations across stellar masses ≳ 0.1 M ⊙. This work demonstrates the JWST's capability in uncovering multiple populations within M dwarfs and illustrates the possibility to analyze very-low-mass stars in GCs approaching the H-burning limit and the brown-dwarf sequence.

EPOCHS. II. The Ultraviolet Luminosity Function from 7.5 < z < 13.5 Using 180 arcmin2 of Deep, Blank Fields from the PEARLS Survey and Public JWST Data

We present an analysis of the ultraviolet luminosity function (UV LF) and star formation rate density of distant galaxies (7.5 < z < 13.5) in the “blank” fields of the Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) survey combined with Early Release Science data from the CEERS, GLASS, and NGDEEP surveys/fields and the first data release of JADES. We use strict quality cuts on EAZY photometric redshifts to obtain a reliable selection and characterization of high-redshift (z > 6.5) galaxies from a consistently processed set of deep, near-infrared imaging. Within an area of 180 arcmin2, we identify 1046 candidate galaxies at redshifts z > 6.5 and we use this sample to study the UV LF in four redshift bins between 7.5 < z < 13.5. The measured number density of galaxies at z = 8 and z = 9 matches those of past observations undertaken by the Hubble Space Telescope (HST). Our z = 10.5 measurements lie between early James Webb Space Telescope (JWST) results and past HST results, indicating cosmic variance may be the cause of previous high density measurements. However, the number densities of UV-luminous galaxies at z = 12.5 are high compared to predictions from simulations. When examining the star formation rate density of galaxies at this period, our observations are still largely consistent with a constant star formation efficiency, are slightly lower than previous early estimations using JWST, and support galaxy driven reionization at z ≤ 8.

COSMOS-Web: Intrinsically Luminous z ≳ 10 Galaxy Candidates Test Early Stellar Mass Assembly

We report the discovery of 15 exceptionally luminous 10 ≲ z ≲ 14 candidate galaxies discovered in the first 0.28 deg2 of JWST/NIRCam imaging from the COSMOS-Web survey. These sources span rest-frame UV magnitudes of −20.5 > M UV > −22, and thus constitute the most intrinsically luminous z ≳ 10 candidates identified by JWST to date. Selected via NIRCam imaging, 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 ≳ 10 sources and 3/15 (20%) likely low-redshift contaminants. Three of our z ∼ 12 candidates push the limits of early stellar mass assembly: they have estimated stellar masses ∼ 5 × 109 M ⊙, implying an effective stellar baryon fraction of ϵ ⋆ ∼ 0.2−0.5, where ϵ ⋆ ≡ M ⋆/(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 ⋆ ∼ 1010 M ⊙ galaxies relative to M ⋆ ∼ 109 M ⊙—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 UV luminosity function from a double power law to a Schechter function at z ≈ 8. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understand how, and if, such early massive galaxies push the limits of galaxy formation in the Lambda cold dark matter paradigm.

Mapping the Growth of Supermassive Black Holes as a Function of Galaxy Stellar Mass and Redshift

The growth of supermassive black holes is strongly linked to their galaxies. It has been shown that the population mean black hole accretion rate ( BHAR¯ ) primarily correlates with the galaxy stellar mass (M ⋆) and redshift for the general galaxy population. This work aims to provide the best measurements of BHAR¯ as a function of M ⋆ and redshift over ranges of 109.5 < M ⋆ < 1012 M ⊙ and z < 4. We compile an unprecedentedly large sample with 8000 active galactic nuclei (AGNs) and 1.3 million normal galaxies from nine high-quality survey fields following a wedding cake design. We further develop a semiparametric Bayesian method that can reasonably estimate BHAR¯ and the corresponding uncertainties, even for sparsely populated regions in the parameter space. BHAR¯ is constrained by X-ray surveys sampling the AGN accretion power and UV-to-infrared multiwavelength surveys sampling the galaxy population. Our results can independently predict the X-ray luminosity function (XLF) from the galaxy stellar mass function (SMF), and the prediction is consistent with the observed XLF. We also try adding external constraints from the observed SMF and XLF. We further measure BHAR¯ for star-forming and quiescent galaxies and show that star-forming BHAR¯ is generally larger than or at least comparable to the quiescent BHAR¯ .

JWST NIRCam Photometry: A Study of Globular Clusters Surrounding Bright Elliptical Galaxy VV 191a at z = 0.0513

James Webb Space Telescope NIRCam images have revealed 154 reliable globular cluster (GC) candidates around the z = 0.0513 elliptical galaxy VV 191a after subtracting 34 likely interlopers from background galaxies inside our search area. NIRCam broadband observations are made at 0.9–4.5 μm using the F090W, F150W, F356W, and F444W filters. Using point-spread-function-matched photometry, the data are analyzed to present color–magnitude diagrams and color distributions that suggest a relatively uniform population of GCs, except for small fractions of reddest (5%–8%) and bluest (2%–4%) outliers. GC models in the F090W versus (F090W–F150W) diagram fit the NIRCam data well and show that the majority of GCs detected have a mass of ∼106.5 M ⊙, with metallicities [Fe/H] spanning the typical range expected for GCs (−2.5 ≲ [Fe/H]≲ 0.5). However, the models predict ∼0.3–0.4 mag bluer (F356W–F444W) colors than the NIRCam data for a reasonable range of GC ages, metallicities, and reddening. Although our data do not quite reach the luminosity function turnover, the measured luminosity function is consistent with previous measurements, suggesting an estimated peak at m AB ∼ −9.4 ± 0.2 mag in the F090W filter.