Quiescent Galaxies Research Articles

The First Quenched Galaxies: When and How?

Many quiescent galaxies discovered in the early Universe by JWST raise fundamental questions on when and how these galaxies became and stayed quenched. Making use of the latest version of the semianalytic model GAEA that provides good agreement with the observed quenched fractions up to z ∼ 3, we make predictions for the expected fractions of quiescent galaxies up to z ∼ 7 and analyze the main quenching mechanism. We find that in a simulated box of 685 Mpc on a side, the first quenched massive (M ⋆ ∼ 1011 M ⊙), Milky Way–mass, and low-mass (M ⋆ ∼ 109.5 M ⊙) galaxies appear at z ∼ 4.5, z ∼ 6.2, and before z = 7, respectively. Most quenched galaxies identified at early redshifts remain quenched for more than 1 Gyr. Independently of galaxy stellar mass, the dominant quenching mechanism at high redshift is accretion disk feedback (quasar winds) from a central massive black hole, which is triggered by mergers in massive and Milky Way–mass galaxies and by disk instabilities in low-mass galaxies. Environmental stripping becomes increasingly more important at lower redshift.

Beyond the rings: Polar ring galaxy NGC 4262 and its globular cluster system

ABSTRACT In the context of the hierarchical model of galaxy evolution, polar ring galaxies (PRGs) are considered the intermediate phase between ongoing mergers and quiescent galaxies. This study explores the globular cluster system (GCS) and its properties in the nearest PRG, NGC4262, serving as a pilot investigation to study GCS in nearby PRGs. We utilize wide and deep-field observations of the CFHT as part of the NGVS to investigate the GCS of NGC4262. We presented the first optical image of NGC4262 with an optically faint ring component. The photometric analysis of the GCS displays a distinct colour bimodality. We estimate the total number of GCs for NGC4262 to be 266 $\pm$ 16 GCs with a specific frequency of 4.2 $\pm$ 0.8 and a specific mass of 0.23 $\pm$ 0.01, which is relatively high compared to other galaxies of similar mass and environmental conditions. The spatial and azimuthal distributions of subpopulations reveal strong evidence of previous interactions within the host galaxy. The colour distribution of GCS in NGC4262 shows a gradient of –0.05 $\pm$ 0.01 within 5.5$^\prime$, supporting the notion of past interactions and evolutionary transitions. PRG NGC4262 conforms to the overall trend of the GCS mass with respect to the halo mass. Furthermore, our investigation of the global scaling relations between GCS and host galaxy parameters provides further support for the hypothesis that PRGs are an intermediate phase connecting ongoing mergers and quiescent galaxies.

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.

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

Teaming up Radio and Sub-mm/FIR Observations to Probe Dusty Star-Forming Galaxies

In this paper, we investigate the benefits of teaming up data from the radio to the far-infrared (FIR) regime for the characterization of dusty star-forming galaxies (DSFGs). These galaxies are thought to be the star-forming progenitors of local massive quiescent galaxies and to play a pivotal role in the reconstruction of the cosmic star formation rate density up to high redshift. Due to their dust-enshrouded nature, DSFGs are often invisible in the near-infrared/optical/UV bands. Therefore, they necessitate observations at longer wavelengths, primarily the FIR band, where dust emission occurs, and the radio band, which is not affected by dust absorption. Combining data from these two spectral windows makes it possible to characterize even the dustiest objects, enabling the retrieval of information about their age, dust temperature, and star-formation status, and facilitates the differentiation between various galaxy populations that evolve throughout cosmic history. Despite the detection of faint radio sources being a challenging task, this study demonstrates that an effective strategy to build statistically relevant samples of DSFGs would be reaching deep sensitivities in the radio band, even restricted to smaller areas, and then combining these radio observations with FIR/submm data. Additionally, this paper quantifies the improvement in the spectral energy distribution (SED) reconstruction of DSFGs by incorporating ALMA band measurements, in particular, in its upgraded status thanks to the anticipated Wideband Sensitivity Upgrade.

The ALCHEMI Atlas: Principal Component Analysis Reveals Starburst Evolution in NGC 253

Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultrawide millimeter spectral scan survey from the Atacama Large Millimeter/submillimeter Array Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1.″6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH+ and the first interferometric images of C3H+, NO, and HCS+. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our data set is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as (i) young starburst tracers characterized by high-excitation transitions of HC3N and complex organic molecules versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing photodissociation regions, (ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH3OH, HNCO, HOCO+, and OCS) versus shocks from star formation-induced outflows (high-excitation transitions of SiO), as well as (iii) outflows showing emission from HOC+, CCH, H3O+, CO isotopologues, HCN, HCO+, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback.

X-ray stacking reveals average SMBH accretion properties of star-forming galaxies and their cosmic evolution over 4 ≲ z ≲ 7

ABSTRACT With an X-ray stacking analysis of $\simeq 12\, 000$ Lyman-break galaxies (LBGs) using the Chandra Legacy Survey image, we investigate average supermassive black hole (SMBH) accretion properties of star-forming galaxies (SFGs) at 4 ≲ z ≲ 7. Although no X-ray signal is detected in any stacked image, we obtain strong 3σ upper limits for the average black hole accretion rate (BHAR) as a function of star formation rate (SFR). At z ∼ 4 (5) where the stacked image is deeper, the 3σ BHAR upper limits per SFR are ∼1.5 (1.0) dex lower than the local black hole-to-stellar mass ratio, indicating that the SMBHs of SFGs in the inactive (BHAR $\lesssim 1 \, \mathrm{M}_\odot$ yr−1) phase are growing much more slowly than expected from simultaneous evolution. We obtain a similar result for BHAR per dark halo accretion rate. QSOs from the literature are found to have ∼1 dex higher SFRs and ≳ 2 dex higher BHARs than LBGs with the same dark halo mass. We also make a similar comparison for dusty starburst galaxies and quiescent galaxies from the literature. A duty-cycle corrected analysis shows that for a given dark halo, the SMBH mass increase in the QSO phase dominates over that in the much longer inactive phase. Finally, a comparison with the TNG300, TNG100, SIMBA100, and EAGLE100 simulations finds that they overshoot our BHAR upper limits by ≲ 1.5 dex, possibly implying that simulated SMBHs are too massive.

Remarkably Compact Quiescent Candidates at 3 < z < 5 in JWST-CEERS

In this Letter, we measure the rest-frame optical and near-infrared sizes of 10 quiescent candidates at 3 < z < 5, first reported by Carnall et al. We use James Webb Space Telescope Near-Infrared Camera F277W and F444W imaging obtained through the public CEERS Early Release Science program and imcascade, an astronomical fitting code that utilizes multi-Gaussian expansion, to carry out our size measurements. When compared to the extrapolation of rest-optical size–mass relations for quiescent galaxies at lower redshift, 8 out of 10 candidates in our sample (80%) are on average more compact by ∼40%. A total of 7 out of 10 candidates (70%) exhibit rest-frame infrared sizes ∼10% smaller than rest-frame optical sizes, indicative of negative color gradients. Two candidates (20%) have rest-frame infrared sizes ∼1.4× larger than rest-frame optical sizes; one of these candidates exhibits signs of ongoing or residual star formation, suggesting this galaxy may not be fully quenched. The remaining candidate is unresolved in both filters, which may indicate an active galactic nucleus. Strikingly, we observe three of the most massive galaxies in the sample (log(M ⋆/M ⊙) = 10.74–10.95) are extremely compact, with effective radii ∼0.7 kpc. Our findings provide no indication that the size evolution relation flattens out, and may indicate that the size evolution of quiescent galaxies is steeper than previously anticipated beyond z > 3.

The TNG50-SKIRT Atlas: Post-processing methodology and first data release

Galaxy morphology is a powerful diagnostic to assess the realism of cosmological hydrodynamical simulations. Determining the morphology of simulated galaxies requires the generation of synthetic images through 3D radiative transfer post-processing that properly accounts for different stellar populations and interstellar dust attenuation. We use the SKIRT code to generate the TNG50-SKIRT Atlas, a synthetic UV to near-infrared broadband image atlas for a complete stellar-mass selected sample of 1154 galaxies extracted from the TNG50 cosmological simulation at z = 0. The images have a high spatial resolution (100 pc) and a wide field of view (160 kpc). In addition to the dust-obscured images, we also release dust-free images and physical parameter property maps with matching characteristics. As a sanity check and preview application we discuss the UVJ diagram of the galaxy sample. We investigate the effect of dust attenuation on the UVJ diagram and find that it affects both the star-forming and the quiescent galaxy populations. The quiescent galaxy region is polluted by younger and star-forming highly inclined galaxies, while dust attenuation induces a separation in inclination of the star-forming galaxy population, with low-inclination galaxies remaining at the blue side of the diagram and high-inclination galaxies systematically moving towards the red side. This image atlas can be used for a variety of other applications, including galaxy morphology studies and the investigation of local scaling relations. We publicly release the images and parameter maps, and we invite the community to use them.

MAGIC: MUSE gAlaxy Groups In COSMOS – A survey to probe the impact of environment on galaxy evolution over the last 8 Gyr

Context. Galaxies migrate along filaments of the cosmic web from small groups to clusters, which creates the appearance that the evolution of their properties speeds up as environments get denser. Aims. We introduce the MUSE gAlaxy Groups in COSMOS (MAGIC) survey, which was built to study the impact of environment on galaxy evolution down to low stellar masses over the last 8 Gyr. Methods. The MAGIC survey consists of 17 Multi-Unit Spectrocopic Exporer (MUSE) fields targeting 14 massive, known structures at intermediate redshift (0.3 &lt; z &lt; 0.8) in the COSMOS area, with a total on-source exposure of 67 h. We securely measured the redshifts for 1419 sources and identified 76 galaxy pairs and 67 groups of at least three members using a friends-of-friends algorithm. The environment of galaxies is quantified from group properties, as well as from global and local density estimators. Results. The MAGIC survey has increased the number of objects with a secure spectroscopic redshift over its footprint by a factor of about 5 compared to previous extensive spectroscopic campaigns on the COSMOS field. Most of the new redshifts have apparent magnitudes in the z++ band zapp++ &gt; 21.5. The spectroscopic redshift completeness is high: in the redshift range of [O II] emitters (0.25 ≤ z &lt; 1.5), where most of the groups are found, it globally reaches a maximum of 80% down to zapp++ = 25.9, and locally decreases from ∼100% to ∼50% in magnitude bins from zapp++ = 23−24 to zapp++ = 25.5. We find that the fraction of quiescent galaxies increases with local density and with the time spent in groups. A morphological dichotomy is also found between bulge-dominated quiescent and disk-dominated star-forming galaxies. As environment gets denser, the peak of the stellar mass distribution shifts towards M⋆ &gt; 1010 M⊙, and the fraction of galaxies with M⋆ &lt; 109 M⊙ decreases significantly, even for star-forming galaxies. We also highlight peculiar features such as close groups, extended nebulae, and a gravitational arc. Conclusions. Our results suggest that galaxies are preprocessed in groups of increasing mass before entering rich groups and clusters. We publicly release two catalogs containing the properties of galaxies and groups, respectively.

Environmental dependence of AGN activity and star formation in galaxy clusters from Magneticum simulations

Context. The environment inside and on the outskirts of galaxy clusters has a profound impact on the star formation rate and active galactic nucleus (AGN) activity in cluster galaxies. While the overall star formation and AGN suppression in the inner cluster regions has been thoroughly studied in the past, recent X-ray studies also indicate that conditions on the cluster outskirts may promote AGN activity. Aims. We investigate how the environment and the properties of host galaxies impact the levels of AGN activity and star formation in galaxy clusters. We aim to identify significant trends in different galaxy populations and suggest possible explanations. Methods. We studied galaxies with stellar mass log M*(M⊙) &gt; 10.15 in galaxy clusters with mass M500 &gt; 1013 M⊙ extracted from box2b (640 comoving Mpc h−1) of the Magneticum Pathfinder suite of cosmological hydrodynamical simulations at redshifts 0.25 and 0.90. We examined the influence of stellar mass, distance to the nearest neighbouring galaxy, cluster-centric radius, substructure membership, and large-scale surroundings on the fraction of galaxies hosting an AGN, star formation rate, and the ratio between star-forming and quiescent galaxies. Results. We find that in low-mass galaxies, AGN activity and star formation are similarly affected by the environment and decline towards the cluster centre. In massive galaxies, the impact is different; star-formation level increases in the inner regions and peaks between 0.5 and 1 R500 with a rapid decline in the centre, whereas AGN activity declines in the inner regions and rapidly rises below R500 towards the centre. We suggest that this increase is a result of the larger black hole masses relative to stellar masses in the cluster centre. After disentangling the contributions of neighbouring cluster regions, we find an excess of AGN activity in massive galaxies on the cluster outskirts (∼3 R500). We also find that the local density, substructure membership, and stellar mass strongly influence star formation and AGN activity but verify that they cannot fully account for the observed radial trends.

At the end of cosmic noon: Short gas depletion times in unobscured quasars at z ∼ 1

Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas–rich mergers to gas–depleted, quenched galaxies. Studies of this population, however, find a high incidence of far–infrared–luminous sources–suggesting significant dust-obscured star formation–but direct observations of the cold molecular gas fuelling this star formation are still necessary. We present a NOEMA study of CO(2–1) emission, tracing the cold molecular gas, in ten lensed z = 1 − 1.5 unobscured QSOs. We detected CO(2–1) in seven of our targets, four of which also show continuum emission (λrest = 1.3 mm). After subtracting the foreground galaxy contribution to the photometry, spectral energy distribution fitting yielded stellar masses of 109 − 11 M⊙, with star formation rates of 25−160 M⊙ yr−1 for the host galaxies. These QSOs have lower LCO′ than star–forming galaxies with the same LIR, and show depletion times spanning a large range (50−900 Myr), but with a median of just 90(αCO/4) Myr. We find molecular gas masses in the range ≤2−40 × 109(αCO/4) M⊙, which suggest gas fractions above ∼50% for most of the targets. Despite the presence of an unobscured QSO, the host galaxies are able to retain significant amounts of cold gas. However, with a median depletion time of ∼90 Myr, the intense burst of star formation taking place in these targets will quickly deplete their molecular gas reservoirs in the absence of gas replenishment, resulting in a quiescent host galaxy. The non–detected QSOs are three of the four radio–loud QSOs in the sample, and their properties indicate that they are likely already transitioning into quiescence. Recent cosmological simulations tend to overestimate the depletion times expected for these z ∼ 1 QSO–host galaxies, which is likely linked to their difficulty producing starbursts across the general high-redshift galaxy population.

Faint mm NIKA2 dusty star-forming galaxies: Finding the high-redshift population

Aims. High-redshift dusty star-forming galaxies (DSFGs) are proposed to be the progenitors of massive quiescent galaxies arising at cosmic noon, providing a crucial insight into the formation, assembly, and early quenching of massive galaxies in the early Universe. However, their high redshift combined with high dust obscuration adds significant difficulties to their redshift measurement, which is mandatory for detailed studies of their physical properties. Blind mm spectral scans are the most unbiased way in prinicple for obtaining accurate spectroscopic redshifts for these sources, but identifying faint molecular and atomic lines within limited telescope time for faint DSFGs is also difficult with these scans. Methods. We developed a new framework to constrain the source redshift. The method jointly accounts for the detection and/or nondetection of spectral lines and the prior information from the photometric redshift and total infrared luminosity from spectral energy distribution analysis. The method uses the estimated total infrared luminosity to predict the line fluxes at given redshifts and generates model spectra. The redshift-dependent spectral models were then compared with the observed spectra to determine the redshift. Results. We applied this joint redshift analysis method to four high-z dusty star-forming galaxy candidates selected from the NIKA2 observations of the HLSJ091828.6+514223 (HLS) field that were further observed by NOEMA with blind spectral scans. These sources only have Herschel SPIRE photometry as ancillary data. They were selected because SPIRE counterparts are faint or entirely lacking and thus favor to select the highest-redshift candidates. The method finds a spectroscopic redshift of 4 in the five NOEMA-counterpart detected sources, with z &gt; 3. Based on these measurements, we derived the CO and [CI] lines and mm continuum fluxes from the NOEMA data and studied the properties of their interstellar medium and star formation. We find cold dust temperatures in some of the HLS sources compared to the general population of submm galaxies, which might be related to the bias introduced by the SPIRE-dropout selection. All sources except for one have a short gas-depletion time of a few hundred million years, which is typical of high-z submm galaxies. The only exception shows a longer gas-depletion time of up to a few billion years. This is comparable to the gas-depletion times of main-sequence galaxies at the same redshift. Furthermore, we identify a possible overdensity of dusty star-forming galaxies at z = 5.2 that is traced by two sources in our sample, as well as a lensed galaxy HLSJ091828.6+514223. Conclusions. We demonstrate that our method when applied to mm-selected DSFGs is able to determine the redshift accurately. This accuracy with only multiple emission lines with a low signal-to-noise ratio shows promising potential for the blind redshift search in large samples of high-z DSFGs, even in the absence of optical to near infrared photometric redshifts.

On the Existence, Rareness, and Uniqueness of Quenched H i-rich Galaxies in the Local Universe

Using data from ALFALFA, xGASS, H i-MaNGA, and the Sloan Digital Sky Survey (SDSS), we identify a sample of 47 “red but H i-rich” (RR) galaxies with near-UV (NUV) − r > 5 and unusually high H i-to-stellar mass ratios. We compare the optical properties and local environments between the RR galaxies and a control sample of “red and H i-normal” (RN) galaxies that are matched in stellar mass and color. The two samples are similar in the optical properties typical of massive red (quenched) galaxies in the local Universe. The RR sample tends to be associated with slightly lower-density environments and has lower clustering amplitudes and smaller neighbor counts at scales from several hundred kiloparsecs to a few megaparsecs. The results are consistent with the RR galaxies being preferentially located at the center of low-mass halos, with a median halo mass ∼1012 h −1 M ⊙ compared to ∼1012.5 h −1 M ⊙ for the RN sample. This result is confirmed by the SDSS group catalog, which reveals a central fraction of 89% for the RR sample, compared to ∼60% for the RN sample. If assumed to follow the H i size–mass relation of normal galaxies, the RR galaxies have an average H i-to-optical radius ratio of R HI/R 90 ∼ 4, four times the average ratio for the RN sample. We compare our RR sample with similar samples in previous studies, and quantify the population of RR galaxies using the SDSS complete sample. We conclude that the RR galaxies form a unique but rare population, accounting for only a small fraction of the massive quiescent galaxy population. We discuss the formation scenarios of the RR galaxies.

Cosmic Vine: A z = 3.44 large-scale structure hosting massive quiescent galaxies

We report the discovery of a large-scale structure at z = 3.44 revealed by JWST data in the Extended Groth Strip (EGS) field. This structure, called the Cosmic Vine, consists of 20 galaxies with spectroscopic redshifts at 3.43 &lt; z &lt; 3.45 and six galaxy overdensities (4 − 7σ) with consistent photometric redshifts, making up a vine-like structure extending over a ∼4 × 0.2 pMpc2 area. The two most massive galaxies (M* ≈ 1010.9 M⊙) of the Cosmic Vine are found to be quiescent with bulge-dominated morphologies (B/T &gt; 70%). Comparisons with simulations suggest that the Cosmic Vine would form a cluster with halo mass Mhalo &gt; 1014 M⊙ at z = 0, and the two massive galaxies are likely forming the brightest cluster galaxies (BCGs). The results unambiguously reveal that massive quiescent galaxies can form in growing large-scale structures at z &gt; 3, thus disfavoring the environmental quenching mechanisms that require a virialized cluster core. Instead, as suggested by the interacting and bulge-dominated morphologies, the two galaxies are likely quenched by merger-triggered starburst or active galactic nucleus (AGN) feedback before falling into a cluster core. Moreover, we found that the observed specific star formation rates of massive quiescent galaxies in z &gt; 3 dense environments are one to two orders of magnitude lower than that of the BCGs in the TNG300 simulation. This discrepancy potentially poses a challenge to the models of massive cluster galaxy formation. Future studies comparing a large sample with dedicated cluster simulations are required to solve the problem.

A Massive Quiescent Galaxy in a Group Environment at z = 4.53

We report on the spectroscopic confirmation of a massive quiescent galaxy at z spec = 4.53 in the COSMOS field. The object was first identified as a galaxy with suppressed star formation at z phot ∼ 4.65 from the COSMOS2020 catalog. The follow-up spectroscopy with Keck/MOSFIRE in the K band reveals faint [O ii] emission and the Balmer break, indicative of evolved stellar populations. We fit the spectral energy distribution using photometry and a spectrum to infer physical properties. The obtained stellar mass is high (M * ∼ 1010.8 M ⊙) and the current star formation rate is more than 1 dex below that of main-sequence galaxies at z = 4.5. Its star formation history suggests that this galaxy experienced rapid quenching from z ∼ 5. The galaxy is among the youngest quiescent galaxies confirmed so far at z spec > 3 with z form ∼ 5.2 (200 Myr ago), which is the epoch when 50% of the total stellar mass was formed. A unique aspect of the galaxy is that it is in an extremely dense region; there are four massive star-forming galaxies at 4.4 < z phot < 4.7 located within 150 physical kpc from the galaxy. Interestingly, three of them have virial radii that strongly overlap with that of the central quiescent galaxy (∼70 kpc), suggesting that the overdensity region is likely the highest-redshift candidate of a dense group with a spectroscopically confirmed quiescent galaxy at the center. The group provides us with a unique opportunity to gain insights into the role of the group environment in quenching at z ∼ 5, which corresponds to the formation epoch of massive elliptical galaxies in the local Universe.

PROVABGS: The Probabilistic Stellar Mass Function of the BGS One-percent Survey

We present the probabilistic stellar mass function (pSMF) of galaxies in the DESI Bright Galaxy Survey (BGS), observed during the One-percent Survey. The One-percent Survey was one of DESI’s survey validation programs conducted from 2021 April to May, before the start of the main survey. It used the same target selection and similar observing strategy as the main survey and successfully observed the spectra and redshifts of 143,017 galaxies in the r < 19.5 magnitude-limited BGS Bright sample and 95,499 galaxies in the fainter surface-brightness- and color-selected BGS Faint sample over z < 0.6. We derive pSMFs from posteriors of stellar mass, M *, inferred from DESI photometry and spectroscopy using the Hahn et al. PRObabilistic Value-Added BGS (PROVABGS) Bayesian spectral energy distribution modeling framework. We use a hierarchical population inference framework that statistically and rigorously propagates the M * uncertainties. Furthermore, we include correction weights that account for the selection effects and incompleteness of the BGS observations. We present the redshift evolution of the pSMF in BGS, as well as the pSMFs of star-forming and quiescent galaxies classified using average specific star formation rates from PROVABGS. Overall, the pSMFs show good agreement with previous stellar mass function measurements in the literature. Our pSMFs showcase the potential and statistical power of BGS, which in its main survey will observe >100 × more galaxies. Moreover, we present the statistical framework for subsequent population statistics measurements using BGS, which will characterize the global galaxy population and scaling relations at low redshifts with unprecedented precision.

The Emergence of a Brightest Cluster Galaxy in a Protocluster Core at z = 2.24

We report the detection of a pair of massive quiescent galaxies likely in the process of merging at the center of the spectroscopically confirmed, extremely massive protocluster BOSS1244 at z = 2.24 ± 0.02. These galaxies, BOSS1244-QG1 and BOSS1244-QG2, were detected with Hubble Space Telescope grism slitless spectroscopic observations. These two quiescent galaxies are among the brightest member galaxies, with z = 2.223–2.255 in BOSS1244, and reside at redshifts z = 2.244 and z = 2.242, with a half-light radius of 6.76 ± 0.50 kpc and 2.72 ± 0.16 kpc, respectively. BOSS1244-QG1 and BOSS1244-QG2 are separated by a projected distance of about 70 physical kpc, implying that the two galaxies likely merge to form a massive brightest cluster galaxy (BCG) with size and mass similar to the most massive BCGs in the local Universe. We thus infer that BCG formation through dry major mergers may happen earlier than the full assembly of a cluster core, which broadens our previous understanding of the coevolution of mature galaxy clusters and BCGs in the nearby Universe. Moreover, we find a strong density–star formation relation over a scale of ∼18 comoving Mpc in BOSS1244, i.e., star formation activity decreases as density increases, implying that the quenching of star formation in BCGs and their progenitors is likely governed by environment-related processes before the virialization of the cluster core.

A massive galaxy that formed its stars at z ≈ 11.

The formation of galaxies by gradual hierarchical co-assembly of baryons and cold dark matter halos is a fundamental paradigm underpinning modern astrophysics1,2 and predicts a strong decline in the number of massive galaxies at early cosmic times3-5. Extremely massive quiescent galaxies (stellar masses of more than 1011 M⊙) have now been observed as early as 1-2 billion years after the Big Bang6-13. These galaxies are extremely constraining on theoretical models, as they had formed 300-500 Myr earlier, and only some models can form massive galaxies this early12,14. Here we report on the spectroscopic observations with the JWST of a massive quiescent galaxy ZF-UDS-7329 at redshift 3.205 ± 0.005. It has eluded deep ground-based spectroscopy8, it is significantly redder than is typical and its spectrum reveals features typical of much older stellar populations. Detailed modelling shows that its stellar population formed around 1.5 billion years earlier in time (z ≈ 11) at an epoch when dark matter halos of sufficient hosting mass had not yet assembled in the standard scenario4,5. This observation may indicate the presence of undetected populations of early galaxies and the possibility of significant gaps in our understanding of early stellar populations, galaxy formation and the nature of dark matter.

A population of faint, old, and massive quiescent galaxies at 3

Here we present a sample of 12 massive quiescent galaxy candidates at z∼3-4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z\\sim 3-4$$\\end{document} observed with the James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec). These galaxies were pre-selected from the Hubble Space Telescope imaging and 10 of our sources were unable to be spectroscopically confirmed by ground based spectroscopy. By combining spectroscopic data from NIRSpec with multi-wavelength imaging data from the JWST Near Infrared Camera (NIRCam), we analyse their stellar populations and their formation histories. We find that all of our galaxies classify as quiescent based on the reconstruction of their star formation histories but show a variety of quenching timescales and ages. All our galaxies are massive (∼0.1-1.2×1011\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 0.1-1.2\ imes 10^{11}$$\\end{document} M⊙\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_\\odot$$\\end{document}), with masses comparable to massive galaxies in the local Universe. We find that the oldest galaxy in our sample formed ∼1.0×1011\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 1.0\ imes 10^{11}$$\\end{document} M⊙\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_\\odot$$\\end{document} of mass within the first few hundred million years of the Universe and has been quenched for more than a billion years by the time of observation at z∼3.2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z\\sim 3.2$$\\end{document} (∼2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 2$$\\end{document} billion years after the Big Bang). Our results point to very early formation of massive galaxies requiring a high conversion rate of baryons to stars in the early Universe.