---
title: 'NICE: NOEMA Forming Cluster Survey'
url: https://www.emergentmind.com/topics/noema-forming-cluster-survey-nice
type: topic
---

# NICE: NOEMA Forming Cluster Survey

Searching arXiv for the NICE survey papers and related records.
Massive protoclusters at $z\sim1.5$–4 occupy the peak of the cosmic star formation history and are therefore central to the study of how massive galaxies in present-day clusters assemble. The Noema formIng Cluster survEy (NICE) is a NOEMA Large Program with a complementary ALMA program designed to build a homogeneous, statistically selected, spectroscopically confirmed sample of forming cluster cores and compact groups at high redshift, and to characterize their star formation, molecular gas, baryons, and dark matter through CO lines, millimeter continuum, and uniform ancillary data [2507.09971]. Across the survey literature, NICE is defined by a uniform candidate selection based on overdensities of red IRAC sources spatially coincident with red Herschel colors, followed by mm-wave spectroscopic confirmation and multi-method halo-mass inference [2407.02973].

## 1. Survey definition, scope, and scientific rationale

NICE is described as a two-facility millimeter spectroscopic program designed to build a homogeneous, statistically powerful census of massive groups and proto-clusters in the early Universe and to characterize their baryons and dark matter [2407.02973]. The Large Programs comprise 159 hr with NOEMA and 40 hr with ALMA. One survey description states that 48 overdensities were selected across five extragalactic legacy fields for NOEMA Band 1 spectroscopy, with an ALMA Cycle 8 program complementing southern targets [2310.15925]. A later survey paper states that the program targets 69 massive galaxy group candidates at $z>2$ in six deep fields—COSMOS, Lockman Hole, Elais-N1, Boötes, XMM-LSS, and ECDFS—over $46\,\mathrm{deg}^2$ [2407.02973]. This suggests that NICE evolved from an initial operational sample to a broader multi-field program.

The primary scientific motivation is the lack of a statistically and homogeneously selected and spectroscopically confirmed sample of massive protoclusters at $2\lesssim z\lesssim4$, where cold gas accretion and concentrated star formation are expected to be efficient [2310.15925]. Earlier proto-cluster work is described as heterogeneous in tracer choice, often based on radio galaxies, narrow-band emission-line overdensities, or Ly$\alpha$ tomography, and therefore prone to differences in selection function and evolutionary stage. NICE was designed to mitigate those limitations by adopting a single color-plus-overdensity selection and a common spectroscopic confirmation strategy [2507.09971].

Within COSMOS, NICE-COSMOS consists of eight confirmed protoclusters at $1.5<z<4$, analyzed within $r\leq1\,R_{\rm vir}$, with virial radii of 141–429 pkpc and 13–31 members per structure [2507.09971]. The COSMOS pilot paper reports that these eight systems are confirmed at $1.65\leq z\leq3.61$ and have best halo-mass estimates of $\log(M_{\rm h}/{\rm M_\odot})=12.8$–13.7 with uncertainty of 0.3 dex [2407.02973]. The Lockman Hole discovery paper adds an even higher-redshift compact group, LH-SBC3 at $z=3.95$, confirmed via CO(4–3) and [CI](1–0) in four massive galaxies [2310.15925].

## 2. Candidate selection and spectroscopic confirmation

NICE candidate selection combines an IRAC overdensity criterion with Herschel/SPIRE color cuts. The IRAC selection requires
$$
[3.6]-[4.5] > 0.1,\qquad 20<[4.5]<23,
$$
with a surface-density overdensity threshold $\Sigma_N>5\sigma$, where
$$
\Sigma_N = \frac{N}{\pi d_N^2},
$$
and $N=5$ or 10 nearest neighbors [2507.09971]. The SPIRE “350 $\mu$m peaker” cuts are
$$
S_{500}>30~{\rm mJy},\qquad S_{350}/S_{250}>1.06,\qquad S_{500}/S_{350}>0.72
$$
in the survey overview, with one COSMOS pilot description listing $S_{350}/S_{250}>1.07$ [2407.02973]. Both formulations appear in the NICE literature. The stated motivation is that these criteria select overdense, intensively star-forming massive environments at $z\approx2$–4 [2507.09971].

This dual selection is meant to identify dusty, star-forming structures while reducing projection contamination relative to submillimeter overdensities alone [2310.15925]. At the same time, the survey literature explicitly notes selection effects: the IRAC/SPIRE-based targeting emphasizes intensively star-forming protoclusters and may over-represent systems in active growth phases [2507.09971]. Severe blending can also suppress the apparent IRAC overdensity in compact cores, as discussed for HPC1001, which satisfied the red Herschel color cuts but was not an IRAC overdensity because of source blending at IRAC resolution [2407.02973].

Spectroscopic confirmation is carried out with NOEMA and ALMA through CO and [CI] line detections, with Subaru/FMOS H$\alpha$ used for the $z=1.65$ COSMOS system where the ALMA tunings did not cover strong mm lines [2407.02973]. NOEMA observations at 3 mm were designed with two frequency setups per target, covering CO(3–2) and CO(4–3) over $2\lesssim z\lesssim4$, with typical line sensitivity of 0.13 mJy beam$^{-1}$ over 500 km s$^{-1}$ at $\sim100$ GHz, continuum sensitivity of $\sim17\,\mu$Jy beam$^{-1}$, and $\sim4.1''$ synthesized beams [2407.02973]. ALMA Bands 4 and 5 used four tunings spanning 135–183 GHz, targeting CO(3–2), CO(4–3), CO(5–4), and [CI]$(^1P_1$–$^3P_0)$ depending on redshift, with line sensitivity of 0.13 mJy beam$^{-1}$ over 500 km s$^{-1}$ at $\sim150$ GHz, continuum $\sim16\,\mu$Jy beam$^{-1}$, and typical resolution $\sim0.7''$ [2407.02973].

A unified NOEMA+ALMA line-search pipeline extracted spectra directly in the uv domain with GILDAS uvfit, adopted a $\nu^{3.7}$ power law for the continuum, and searched for lines with an integrated-S/N algorithm [2407.02973]. In the COSMOS pilot sample, 22 significant lines—CO(3–2), CO(4–3), CO(5–4), and [CI](1–0)—were detected in 20 galaxies, with typical $P_{\rm chance}<0.001$ over the full band [2407.02973]. Photometric-redshift PDFs from COSMOS2020 were used to disambiguate single-line identifications, and 58% (95%) of spectroscopic solutions lie within the $1\sigma$ ($2\sigma$) $z_{\rm phot}$ PDF [2407.02973].

## 3. Observational datasets and measurement framework

NICE combines mm-wave spectroscopy with broad ancillary datasets. The survey overview lists COSMOS2020 photometry and LePhare-based stellar masses and SFRs, VLA 3 GHz, the super-deblended FIR catalog, Herschel/SPIRE, and JWST/COSMOS-Web imaging for two structures [2507.09971]. In crowded cores, integrated FIR/sub-mm fluxes are measured with a dedicated super-deblending procedure: one prior per group is placed at the SCUBA-2 850 $\mu$m peak, sources within $10''$ are excluded from the ancillary prior list, PSF fitting is performed on Herschel/PACS and SPIRE, SCUBA-2, and MeerKAT images, and $24\,\mu$m totals are built from the sum of member fluxes [2407.02973]. The integrated SEDs are then fit with STARDUST using the Magdis et al. (2012) dust templates at the group spectroscopic redshift, treating $>3\sigma$ points as detections and the rest as $3\sigma$ upper limits [2407.02973].

For member-galaxy star formation rates, the primary estimator is VLA 3 GHz when ${\rm S/N}>3$, using the infrared–radio correlation $q_{\rm IR}(z,M_\ast)$ from Delvecchio et al. (2021) and a spectral index $\alpha=-0.8$; the conversion to SFR follows Kennicutt (1998), adjusted via $L_{\rm IR}$ inferred from radio via IRRC [2507.09971]. For 3 GHz ${\rm S/N}<3$, SFR is derived from $24\,\mu$m fluxes using Wuyts et al. (2008) [2507.09971]. For the integrated quantity $\Sigma_{\rm SFR}/M_{\rm halo}$, NICE uses total SFR from FIR SED fits [2507.09971].

The survey provides explicit cross-checks between SFR estimators for radio detections with ${\rm S/N}>3$. The median ratios $(\mathrm{SFR}_{\rm other}/\mathrm{SFR}_{\rm radio})$ are 0.90 for Wuyts $24\,\mu$m, 1.81 for deblended FIR, 0.37 for COSMOS2020/LePhare, and 1.57 for cluster-average FIR SED scaled to members [2507.09971]. A radio-loud galaxy, COSMOS-SBC6 ID 839791, uses $\mathrm{SFR}_{24\mu{\rm m}}$ to avoid AGN radio contamination because its radio luminosity exceeds the AGN-dominance threshold from Wang 2024 [2507.09971].

Stellar masses are taken from COSMOS2020 using LePhare, with uncertainties as reported in COSMOS2020 and not re-derived in the NICE analysis [2507.09971]. Molecular gas masses are preferentially estimated from Rayleigh–Jeans dust continuum at observed 2–3 mm, corresponding to rest-frame $\sim850\,\mu$m, using the prescription calibrated in Hughes et al. (2017); visibilities are fit with Gaussians for resolved sources or point sources for unresolved sources [2507.09971]. A CO-based cross-check converts CO(3–2) or CO(4–3) to CO(1–0) using the ASPECS main-sequence CO SLED at $z\approx2$–2.7, with $\alpha_{\rm CO}$ from the mass–SFR–metallicity approach of Genzel et al. (2015) and the $\alpha_{\rm CO}$–metallicity scaling of Tacconi et al. (2018) [2507.09971]. The two gas-mass methods agree within uncertainties, with
$$
M_{\rm gas,RJ}\approx1.5\times M_{\rm gas,CO},
$$
a median log ratio $\Delta\log(M_{\rm gas,RJ}-M_{\rm gas,CO})\approx-0.18$ dex, and scatter $\approx0.3$ dex; typical $\alpha_{\rm CO}$ values are 3.52–4.78 $M_\odot\,({\rm K\,km\,s^{-1}\,pc^2})^{-1}$ [2507.09971].

The main derived quantities are defined explicitly:
$$
f_{\rm gas} = \frac{M_{\rm gas}}{M_{\rm gas}+M_\ast},\qquad
\mu_{\rm gas} = \frac{M_{\rm gas}}{M_\ast},\qquad
\tau_{\rm gas} = \frac{M_{\rm gas}}{\rm SFR}.
$$
Main-sequence and starburstiness are parameterized using Schreiber et al. (2015), with
$$
\Delta{\rm MS} = \log {\rm SFR} - \log {\rm SFR}_{\rm MS}(M_\ast,z),
$$
and
$$
R_{\rm SB} = \frac{{\rm SFR}}{{\rm SFR}_{\rm MS}(M_\ast,z)},
$$
where starbursts are defined by $R_{\rm SB}=4$ [2507.09971].

## 4. Halo masses, structural inference, and the COSMOS confirmed sample

NICE introduces six halo-mass estimators and adopts a consensus mass based on three of them [2407.02973]. These methods include a BCG stellar-to-halo mass estimate, two total-stellar-mass to halo-mass conversions, an overdensity-plus-galaxy-bias estimate, and two NFW-based fits to the radial stellar mass density [2407.02973]. The adopted cosmology in the COSMOS pilot is flat with $H_0=70~{\rm km\,s^{-1}\,Mpc^{-1}}$, $\Omega_m=0.27$, and $\Omega_\Lambda=0.73$ [2407.02973]. The NICE census paper uses $H_0=70~{\rm km\,s^{-1}\,Mpc^{-1}}$, $\Omega_m=0.3$, and $\Omega_\Lambda=0.7$, with a Chabrier (2003) IMF [2507.09971].

The overdensity-plus-bias method defines the projected overdensity as
$$
\delta = \frac{\Sigma_{\rm group}-\Sigma_{\rm field}}{\Sigma_{\rm field}},
$$
and uses
$$
M_{\rm h} = \rho \times \delta \times b \times \frac{4}{3}\pi [R_{200}(1+z)]^3,
$$
where $\rho$ is the mean matter density at $z$ and $b$ is the linear bias calibrated from Tinker et al. (2010); because $b$ depends on $M_{\rm h}$, an iterative solution is adopted [2407.02973]. The NFW-based methods fit the projected stellar mass surface density $\Sigma M_\ast(R)$ in annuli from $\sim10$ pkpc to $\sim2$ pMpc, centered on the FIR-peak-weighted barycenter of spectroscopically confirmed dusty members [2407.02973].

A central result is that the radial stellar mass density of all eight COSMOS structures is consistent with an NFW profile, supporting the interpretation that they are collapsed structures hosted by a single dark matter halo [2407.02973]. The profile is written as
$$
\rho(r)=\frac{\rho_s}{(r/r_s)(1+r/r_s)^2}.
$$
Compared to Ludlow et al. (2016), the best-fit concentrations show a scatter of $\sim0.28$ dex and a mean offset of $\sim+0.4$ dex, implying higher concentration on average [2407.02973]. This suggests early assembly and relatively compact halos for several systems.

Across the halo-mass estimators $M_{\rm h}(2a)$, $M_{\rm h}(3)$, and $M_{\rm h}(4)$, the inter-method scatter is 0.2–0.3 dex with mean offsets $\lesssim0.4$ dex. NICE therefore adopts the median of these three as the “best” halo mass with a conservative $\pm0.3$ dex uncertainty [2407.02973]. The eight COSMOS systems span $\log(M_{\rm h}/M_\odot)=12.8$–13.7 with $R_{\rm vir}\approx140$–430 pkpc; examples include HPC1001 $(13.3\pm0.3;\ R_{\rm vir}=181\pm42~{\rm pkpc})$, COS-SBCX1 $(13.7\pm0.3;\ 331\pm70~{\rm pkpc})$, and COS-SBC4 $(13.7\pm0.3;\ 429\pm90~{\rm pkpc})$ [2407.02973].

The COSMOS pilot further derives baryonic accretion rates using the Goerdt et al. (2010) scaling,
$$
{\rm BAR}\simeq137\left(\frac{M_{\rm h}}{10^{12}M_\odot}\right)^{1.15}
\left[\frac{1+z}{1+3}\right]^{2.25} M_\odot~{\rm yr}^{-1},
$$
yielding ${\rm BAR}=(1.2$–$8.4)\times10^3\,M_\odot\,{\rm yr}^{-1}$ for the COSMOS groups [2407.02973]. A simple two-regime model relates SFR and BAR through the theoretical cold-stream threshold $M_{\rm stream}$, and the combined NICE+literature sample yields the empirical relation
$$
{\rm SFR/BAR}=10^{-0.46\pm0.22}\left(M_{\rm stream}/M_{\rm h}\right)^{0.71\pm0.16},
$$
with a scatter of $\sim0.40$ dex [2407.02973]. This provides an explicit observational link between gas supply and star formation in dense environments.

## 5. Star formation and cold gas properties in protocluster members

The main NICE-COSMOS census reports a steep increase in star formation rates per halo mass with redshift in intensively star-forming protoclusters [2507.09971]. Specifically, NICE-COSMOS, together with other NICE prototypes, lies approximately 1–2 dex above field values at $z>2$ when compared to field evolution proportional to $(1+z)^{3-4}$ [2507.09971]. The inferred scaling for (proto)clusters is consistent with
$$
\Sigma_{\rm SFR}/M_{\rm halo}\propto(1+z)^{5.4},
$$
extending trends observed at lower redshift to $z>2$, with the data showing an approximately 1 dex increase from $z\approx2$ to $z\approx4$ [2507.09971].

NICE argues that this enhancement is not driven by a higher starburst fraction. Member galaxies generally follow the star-forming main sequence, and the survey reports a low starburst fraction indicated by the absence of a systematic elevation in the $R_{\rm SB}$ distribution relative to the field [2507.09971]. Instead, the dominant effect is the concentration of massive, gas-rich star-forming galaxies in protocluster cores [2507.09971]. The radial mass segregation is strong: the median stellar mass within $0.1\,R_{\rm vir}$ is approximately $10^{10.6}\,M_\odot$, compared with approximately $10^{9.5}\,M_\odot$ in the outskirts, a difference of more than 1 dex [2507.09971]. Among 31 massive galaxies with $M_\ast>3\times10^{10}\,M_\odot$, 21 lie within $0.3\,R_{\rm vir}$, and all but one are detected in CO [2507.09971].

Detected protocluster galaxies have median gas properties
$$
\mu_{\rm gas,NICE}\approx1.32,\qquad \tau_{\rm gas,NICE}\approx0.78~{\rm Gyr}
$$
[2507.09971]. The most massive protocluster galaxies show enhanced gas reservoirs relative to field galaxies. For $M_\ast>10^{10.8}\,M_\odot$, also described as the most massive protocluster galaxies with $M_\ast>8\times10^{10}\,M_\odot$, $\mu_{\rm gas}$ is typically approximately $2\times$ the field level at $z\approx2.5$, declining to approximately $1.2\times$ at $z\approx4$, after normalization to A3COSMOS field scaling relations that account for $M_\ast$, $z$, and $R_{\rm SB}$ [2507.09971]. These galaxies remain on the main sequence rather than in the starburst regime, implying sustained growth via large gas reservoirs rather than brief starburst episodes [2507.09971].

Less massive members are detected mainly when they are strongly starbursting, and their $\mu_{\rm gas}$ and $\tau_{\rm gas}$ are comparable to field galaxies at the same $M_\ast$, $z$, and $R_{\rm SB}$ [2507.09971]. The survey sensitivity is relevant here: the $3\sigma$ continuum sensitivity corresponds to $M_{\rm gas}\approx10^{11}\,M_\odot$ at $z\approx2.5$, which biases detections toward massive members and lower-mass starbursts [2507.09971]. The mean radiation field $\langle U\rangle$ derived from integrated protocluster FIR SEDs follows main-sequence evolution, with some systems slightly lower, consistent with moderately longer gas depletion times inferred for members [2507.09971].

A notable clarification in the NICE census concerns a notation issue: the abstract threshold “$M_{\rm halo}>8\times10^{10}\,M_\odot$” is explicitly stated to refer to galaxy stellar mass and should be read as $M_\ast>8\times10^{10}\,M_\odot$ [2507.09971]. This addresses a potential misconception, since halo masses of the protoclusters themselves are much larger.

## 6. Early discoveries, environmental interpretation, and limitations

The first published NICE discovery, LH-SBC3 in the Lockman Hole at $z=3.95$, established the feasibility of the survey strategy for compact, IR-luminous structures [2310.15925]. Four compact members within approximately 180 pkpc were confirmed in CO(4–3), two also in [CI](1–0), with an average redshift $z_{\rm av}=3.949$ and velocity offsets of $-54$, $-348$, $+543$, and $-141$ km s$^{-1}$ [2310.15925]. The core has an estimated halo mass of $\sim10^{13}\,M_\odot$ and total SFR of $\sim4000\,M_\odot\,{\rm yr}^{-1}$ [2310.15925]. One member hosts a radio-loud AGN with
$$
L_{1.4\,{\rm GHz,rest}} = 3.0\times10^{25}~{\rm W\,Hz}^{-1},
$$
above the common radio-loud threshold [2310.15925].

Using $r_{41}=0.41\pm0.07$ to convert CO(4–3) to CO(1–0) and adopting a ULIRG-like $\alpha_{\rm CO}=0.8\,M_\odot\,({\rm K\,km\,s^{-1}\,pc^2})^{-1}$, the paper derives molecular gas masses of approximately $1.25\times10^{11}$, $4.2\times10^{10}$, $1.4\times10^{10}$, and $1.8\times10^{10}\,M_\odot$ for the four confirmed members [2310.15925]. Corresponding depletion times are approximately 0.10, 0.022, 0.037, and 0.025 Gyr, identifying the system as a compact starbursting core [2310.15925]. The halo-integrated quantity $\Sigma{\rm SFR}/M_{\rm halo}$ is approximately $(2$–$10)\times10^{-10}\,{\rm yr}^{-1}$, about 2 dex above clusters at $z\sim1$–2 [2310.15925].

The later NICE-COSMOS census indicates that such extreme systems should not be taken as representative of the entire survey population. In COSMOS, the member galaxies generally remain on the field main sequence, and the core enhancement in $\Sigma_{\rm SFR}/M_{\rm halo}$ arises primarily from the concentration of massive, gas-rich main-sequence galaxies rather than from a global excess of starbursts [2507.09971]. This establishes an environmental picture in which the dominant mechanism is not universally short-lived burst activity, but sustained star formation supported by substantial gas reservoirs in the most massive core galaxies [2507.09971]. A plausible implication is that NICE samples multiple evolutionary phases, from compact starbursting cores at $z\sim4$ to gas-rich, main-sequence-dominated protocluster cores at $z\sim2$–3.

The survey literature is explicit about limitations. Candidate selection favors dusty, intensely star-forming systems and may bias the sample toward active phases [2507.09971]. Membership at low stellar mass is based on photometric redshifts and is therefore subject to contamination; the high-mass end is more secure because the most massive and high-SFR members are often CO-detected and spectroscopically confirmed [2507.09971]. Small-number statistics remain important: eight protoclusters in COSMOS are described as substantial but still limited [2507.09971]. Gas masses also carry systematic uncertainties because the dust-continuum method depends on Rayleigh–Jeans calibration and dust properties, while the CO-based method depends on excitation corrections and $\alpha_{\rm CO}$–metallicity scaling; the two methods agree within $\approx0.3$ dex, but absolute $M_{\rm gas}$ uncertainties remain [2507.09971]. Halo masses inherit uncertainties from the use of multiple proxies, and the survey notes that underestimation of $M_{\rm halo}$ at high redshift could bias $\Sigma_{\rm SFR}/M_{\rm halo}$ high [2507.09971].

## 7. Evolutionary significance and place in the literature

NICE places its results in the context of previous studies of cluster and protocluster evolution. The steep evolution of $\Sigma_{\rm SFR}/M_{\rm halo}$, extending to $z>2$ and consistent with $(1+z)^{5.4}$, aligns with extrapolations of lower-redshift trends discussed by Webb, Alberts, and others [2507.09971]. The lack of a global starburst excess is presented as consistent with several H$\alpha$-emitter studies of young protoclusters, including Spiderweb-like systems [2507.09971]. At the same time, NICE differs from some extreme $z\gtrsim4$ starbursting cores such as SPT2349-56 and DRC: by $z\approx2.5$, the most massive NICE core galaxies are gas-richer than field counterparts, while at $z\approx4$ the enhancement is milder, approximately $1.2\times$ [2507.09971]. This suggests environmental effects and selection may evolve with redshift.

The COSMOS pilot compares the confirmed structures with simulations and concludes that all eight are consistent with being progenitors of present-day clusters with $M_{\rm h}(z=0)>10^{14}\,M_\odot$ [2407.02973]. Their most massive member galaxies have stellar masses consistent with brightest-cluster-galaxy progenitors in the TNG300 simulation [2407.02973]. Combined with the NFW-like stellar-mass profiles and the concentration of massive, gas-rich galaxies toward the core, this supports the interpretation that NICE is observing collapsed group-scale halos in late stages of assembly, in which proto-BCGs are rapidly growing [2407.02973].

The NICE census further argues that the concentration of gas-rich, massive galaxies with $M_\ast\gtrsim10^{11}\,M_\odot$ in protocluster cores at $z\approx2$–3, together with $\tau_{\rm gas}\approx0.8$ Gyr, suggests sustained in-situ growth and early assembly of central massive ellipticals in present-day clusters [2507.09971]. The elevated $\Sigma_{\rm SFR}/M_{\rm halo}$ is interpreted as an early, efficient build-up phase in cluster cores, consistent with a “downsizing” scenario [2507.09971]. This suggests that, in the environments selected by NICE, the early formation of massive cluster galaxies is driven by persistent gas-rich star formation within dense cores rather than exclusively by brief starburst episodes or immediate quenching.

Recommended future work in the NICE literature includes increasing spectroscopic completeness at lower stellar mass, obtaining more uniform $\alpha_{\rm CO}$ constraints through low-$J$ CO and metallicity measurements, conducting deeper mm-continuum observations below $M_{\rm gas}\sim10^{11}\,M_\odot$ at $z\approx2$–3, extending mosaics beyond $R_{\rm vir}$ to map environmental gradients and inflows, and comparing NICE systematically with differently selected protocluster samples and simulations of gas accretion and quenching timescales [2507.09971]. These directions follow directly from the survey’s central premise: a homogeneous, CO/continuum-anchored census can turn high-redshift cluster formation from a collection of heterogeneous case studies into a comparative empirical framework.

Source: https://www.emergentmind.com/topics/noema-forming-cluster-survey-nice