---
title: 'CL J1226.9+3332: High‑z Galaxy Cluster'
url: https://www.emergentmind.com/topics/cl-j1226-9-3332
type: topic
---

# CL J1226.9+3332: High‑z Galaxy Cluster

CL J1226.9+3332, also designated RX J1226.9+3332, is a massive galaxy cluster at $z=0.89$ or $z=0.888$ in the cited analyses. It is the highest-redshift object in the NIKA2 Sunyaev-Zel'dovich Large Program and has been observed at millimeter, X-ray, lensing, and optical wavelengths. Across these studies, the cluster is described as fairly regular on large scales but with a disturbed core, and it has served as a proof of concept for multi-wavelength studies of thermodynamics, hydrostatic mass reconstruction, hydrostatic-to-lensing bias, and galaxy evolution at high redshift [1410.2808; 2209.07460; 2111.01685; 0910.0511].

## 1. Astrophysical identification and observational context

CL J1226.9+3332 sits at $z=0.89$, with one multi-probe study quoting $z=0.888$, and is explicitly identified as the highest-redshift object in the NIKA2 SZ Large Program. For the cosmology adopted in the multi-probe hydrostatic-mass analysis, $1'$ on the sky corresponds to $\simeq 466\ \mathrm{kpc}$ at $z=0.89$, and the cluster’s characteristic scale is typically reported as $R_{500}\simeq 900$–$1000\ \mathrm{kpc}$ [2209.07460].

The system has been investigated with NIKA, NIKA2, MUSTANG, Bolocam, NOEMA, Chandra, XMM-Newton, Planck, HST/CLASH, Gemini/GMOS, and HST/ACS. This breadth of coverage has made it a recurrent target for cross-calibration of tSZ, X-ray, and lensing mass estimators, as well as for optical studies of galaxy structure and the Fundamental Plane [1410.2808; 2209.07460; 0910.0511].

A persistent observational theme is the coexistence of large-scale regularity with central disturbance. Large-scale X-ray surface brightness is described as fairly regular or relaxed, whereas higher-resolution tSZ, X-ray temperature, radio, and lensing data reveal a south-west substructure and other signatures of a past merger [1410.2808; 2209.07460].

## 2. Millimeter-wave and tSZ observations

High-resolution tSZ work on CL J1226.9+3332 began with NIKA on the IRAM 30-m telescope. In the full NIKA analysis, the camera operated at 150 and 260 GHz with Gaussian beams of $18.2''$ and $12.0''$, a field of view of $\simeq 1.9\ \mathrm{arcmin}$, and an effective integration time of $7.8$ hours on source acquired with on-the-fly raster scans. Calibration was performed on Uranus, with overall uncertainty $7\%$ at 150 GHz and $12\%$ at 260 GHz. The 260 GHz map showed a single bright point source, PS$_{(260)}$, located $30''$ south-east of the X-ray center, with best-fit fluxes $6.8 \pm 0.7(\mathrm{stat}) \pm 1.0(\mathrm{cal})\ \mathrm{mJy}$ at 260 GHz and $1.9 \pm 0.2(\mathrm{stat}) \pm 0.1(\mathrm{cal})\ \mathrm{mJy}$ at 150 GHz. After subtraction of PS$_{(260)}$, the 150 GHz map revealed a smooth tSZ decrement peaking at $18\sigma$ per beam [1410.2808].

The earlier NIKA proceeding emphasized instrumental performance and feasibility. It reported simultaneous 150/260 GHz operation with on-sky beam sizes $\theta_1=18''$ and $\theta_2=12''$, typical per-detector sensitivities of order $\mathrm{NEFD}_1\approx 14\ \mathrm{mJy}\cdot \mathrm{s}^{1/2}$ at 150 GHz and $\mathrm{NEFD}_2\approx 40\ \mathrm{mJy}\cdot \mathrm{s}^{1/2}$ at 260 GHz, cross-linked on-the-fly raster scans, and total on-source time $\simeq 4\ \mathrm{h}$ at 150 GHz and $\simeq 3\ \mathrm{h}$ at 260 GHz. That dataset yielded a peak tSZ signal-to-noise of $\simeq 15$, with final tSZ maps on $2''\times 2''$ pixels; the 150 GHz map contained the cluster decrement, while the 260 GHz map was used to trace and subtract foreground contaminants [1409.1137].

Later interferometric imaging with NOEMA extended the angular-resolution frontier. The NOEMA study used ten antennas in the most compact D configuration, with $19.5\ \mathrm{h}$ on source, dual sidebands covering 70–78 GHz and 86–94 GHz, and a synthesized beam of $4.8''\times 3.6''$ at position angle $54^\circ$. After source subtraction, the final $1\sigma$ continuum sensitivity was $5.6\ \mu\mathrm{Jy}\ \mathrm{beam}^{-1}$. Dirty-map inspection and spectral-cube line search identified three contaminating point sources: PS9, identified with the brightest cluster galaxy at $r\simeq 2''$ from the tSZ peak, and PS260–N and PS260–S, two submillimetre-galaxy line emitters at $\sim 27''$ south-east. These were modeled with `uv_fit` in GILDAS/MAPPING and removed directly in the uv-plane before construction of a source-clean continuum dataset [2509.14048].

## 3. Pressure reconstruction and thermodynamic structure

The tSZ analyses adopt the standard Compton-$y$ description
$$
y(\theta)=\frac{\sigma_T}{m_e c^2}\int P_e(\ell)\,d\ell,
$$
with observed temperature change
$$
\Delta T(\nu,\theta)=T_{\rm CMB}\,f(\nu)\,y(\theta),
$$
and, to first order neglecting relativistic corrections,
$$
f(x)=x\,\coth\frac{x}{2}-4,\qquad x\equiv \frac{h\nu}{k_B T_{\rm CMB}}.
$$
For NIKA’s bands, the proceeding quoted $f(150\ \mathrm{GHz})\approx -7.6$ and $f(260\ \mathrm{GHz})\approx +1.8$, so the 150 GHz channel measures the decrement while the 260 GHz channel primarily diagnoses dusty or radio contamination [1409.1137].

Radial pressure reconstruction has been performed in both parametric and non-parametric forms. The NIKA analysis modeled the three-dimensional electron pressure with a generalized NFW form,
$$
P_e(r)=P_0\Bigl(\frac r{r_s}\Bigr)^{-c}\Bigl[1+\bigl(r/r_s\bigr)^a\Bigr]^{-(b-c)/a},
$$
while the Chandra density profile was fitted with a simplified Vikhlinin model. In the baseline PPC pressure model, the best-fit normalization and scale radius were reported as $P_0=(18.4 \pm 1.9)\ \mathrm{keV\ cm^{-3}}$ and $r_s=(512 \pm 48)\ \mathrm{kpc}$. A simultaneous MCMC fit of the tSZ plus point-source model to the 150 GHz map, together with a Planck-derived integrated signal in a $15'$ aperture,
$$
Y_{15'}=(0.94\pm0.36)\times10^{-3}\ \mathrm{arcmin}^2,
$$
yielded full three-dimensional profiles of pressure, density, temperature, entropy, and mass [1410.2808].

That reconstruction implied a core temperature $T_e(0)\approx 15\ \mathrm{keV}$, declining to $\approx 5\ \mathrm{keV}$ at $\sim 1.5\ \mathrm{Mpc}$, and an entropy law $K(r)\propto r^\alpha$ with flattening below $\sim 100\ \mathrm{kpc}$ and $K_0 \gtrsim 100\ \mathrm{keV\ cm^2}$. In that interpretation, the core was already thermodynamically disturbed despite the cluster’s relaxed large-scale appearance [1410.2808].

NOEMA enabled a core-to-outskirts pressure reconstruction from $\sim 15$ to $\sim 1500\ \mathrm{kpc}$. In the non-parametric model, the pressure bins were fitted at $r=[15,60,150,270]\ \mathrm{kpc}$, giving $0.163^{+0.153}_{-0.023}$, $0.210^{+0.050}_{-0.019}$, $0.090^{+0.035}_{-0.005}$, and $0.019^{+0.025}_{-0.004}\ \mathrm{keV\ cm^{-3}}$, respectively. A joint gNFW fit to NOEMA, IRAM 30 m/NIKA2, MUSTANG, and Bolocam yielded a continuous $P(r)$ with 16–84% and 2.5–97.5% credible regions covering $0.01$–$1.5\,R_{500}$ [2509.14048].

The cluster also became a worked example for joint X-ray–SZ forward modeling in JoXSZ. In that framework, the Abel-projected pressure profile is convolved with the beam and transfer function, X-ray emissivity is translated to count rate with XSPEC, missing data are omitted from the likelihood, and calibration or background systematics are represented by nuisance parameters. Applied to Chandra plus NIKA data, the full joint fit found $T_{\rm SZ}(r\sim 200\ \mathrm{kpc})\simeq 12.5^{+1.5}_{-1.2}\ \mathrm{keV}$ and $T_{\rm SZ}(r_{500})\simeq 7.0^{+1.0}_{-0.8}\ \mathrm{keV}$, while the posterior on the temperature-ratio parameter implied
$$
\frac{T_X}{T_{\rm SZ}}\approx 0.79^{+0.20}_{-0.18},
$$
with a Savage–Dickey Bayes factor $B\sim 3.1$, characterized as weak, inconclusive evidence for a true offset [2005.12315].

## 4. Mass scale, hydrostatic reconstruction, and lensing comparison

Hydrostatic analyses use the standard equilibrium relation
$$
M_{\rm HSE}(<r)=-\,\frac{r^2}{G\,\mu\,m_p\,n_e(r)}\,\frac{dP_e(r)}{dr},
$$
with
$$
M_{500}=\frac{4\pi}{3}\,500\,\rho_c(z)\,R_{500}^3.
$$
For CL J1226.9+3332, published values of $M_{500}$ span a relatively narrow range when tSZ and X-ray information are combined, but the multi-probe literature repeatedly notes that slight differences in the shape of the mass profile near $R_{500}$ can substantially shift the integrated mass [2209.07460].

| Analysis | Reported result | Label |
|---|---:|---|
| NIKA + Chandra + Planck, baseline PPC | $R_{500}=930^{+50}_{-43}\ \mathrm{kpc}$; $M_{500}=(5.96^{+1.02}_{-0.79})\times10^{14}\ M_\odot$ | tSZ/X-ray/Planck |
| JoXSZ full joint fit, HE assumed | $r_{500}=910^{+89}_{-72}\ \mathrm{kpc}$; $M_{500}=(5.57^{+1.81}_{-1.23})\times10^{14}\ M_\odot$ | Joint forward model |
| SZ + X, gNFW approach | $M_{500}^{\rm HSE}=(6.25^{+0.59}_{-0.60})\times10^{14}\ M_\odot$ | Multi-probe |
| SZ + X, NFW backward approach | $M_{500}^{\rm HSE}=(7.30^{+1.17}_{-1.17})\times10^{14}\ M_\odot$ | Multi-probe |
| X-ray-only HSE | $M_{500}^{\rm HSE}=(4.83^{+0.98}_{-0.96})\times10^{14}\ M_\odot$ | Multi-probe |
| Lensing | $M_{500}^{\rm lens}=(7.35^{+0.53}_{-0.51})\times10^{14}\ M_\odot$ | CLASH-based |
| Combined variants, proof of concept | $M_{500}^{\rm HSE}=(7.65 \pm 1.03)\times10^{14}\ M_\odot$; $M_{500}^{\rm lens}=(7.35 \pm 0.65)\times10^{14}\ M_\odot$ | NIKA2 LPSZ proof of concept |

The corresponding gas fractions are also consistent at the $\sim 0.13$–$0.15$ level. The NIKA+Chandra+Planck study obtained
$$
f_{\rm gas}(R_{500})=0.146^{+0.041}_{-0.030},
$$
whereas the JoXSZ full joint fit found
$$
M_{\rm gas}(<r_{500})\simeq 7.2^{+2.3}_{-1.8}\times10^{13}\ M_\odot,\qquad
f_{\rm gas}(r_{500})\simeq 0.13\pm0.02.
$$
An earlier NIKA proceeding, using a self-similar $Y_{500}$–$M_{500}$ scaling relation, estimated $M_{500}\simeq 6\times10^{14}\ M_\odot$ with statistical $\pm 15\%$ and $f_{\rm gas}\approx 0.12\ (\pm 0.02)$ [1410.2808; 2005.12315; 1409.1137].

The hydrostatic-to-lensing bias has been treated as a central systematic question. In the exhaustive multi-probe analysis, the observational estimates were $1-b=0.86^{+0.10}_{-0.11}$ for the gNFW HSE reconstruction, $1-b=1.01^{+0.17}_{-0.19}$ for the NFW backward reconstruction, and $1-b=0.66^{+0.14}_{-0.14}$ for the X-ray-only HSE reconstruction. The range $1-b\approx 0.7$–$1$ was presented as the net effect of data-set choice and mass-profile modeling. The earlier proof-of-concept analysis instead summarized the comparison as
$$
b=-0.04\pm0.20,
$$
that is, consistent with zero within 20%. Both studies emphasized that model dependence dominates instrumental systematics in the final mass budget, while transfer-function choices in the NIKA2 analysis were found to have $<10\%$ impact on $P_e(r)$ and a $\sim 5$–$10\%$ effect on the final mass [2209.07460; 2111.01685].

## 5. Morphology, core state, and dynamical interpretation

The 2014 NIKA map showed a cluster that was azimuthally smooth on scales $\gtrsim 1'$ but elongated to the south-west on $\sim 20''$–$30''$ scales. That elongation was associated with three previously reported features: an XMM/Chandra hotter south-west region, a MUSTANG 90 GHz ridge $\sim 20''$ long and $10''$ south-west of the center, and an HST lensing subclump $\sim 40''$ south-west. The interpretation given there was a minor merger that had stripped or heated core gas [1410.2808].

Subsequent multi-probe work reinforced the disturbed-core picture. XMM-Newton and Chandra reveal a hot, south-west substructure $\sim 40''$ from the X-ray peak, and optical plus LOFAR data show a secondary mass clump and a giant radio halo, presented as signatures of a past merger. At the same time, the X-ray surface brightness remained fairly regular on larger scales, preserving the recurrent description of a globally regular cluster with central dynamical complexity [2209.07460].

The NOEMA core mapping sharpened this interpretation. It provided the first tSZ-based map of the core, $r\lesssim 40\ \mathrm{kpc}$, in a $z=0.89$ cluster and reported that the pressure profile is flat in the core of the cluster. The best-fit gNFW inner slope was $\gamma \simeq 0.7$, corresponding to a shallow pressure gradient inside $\sim 0.1\,R_{500}$. Residual imaginary visibilities pointed to departures from circular symmetry, including sloshing and substructures, and the combined evidence from weak-lensing mass subclumps, X-ray temperature asymmetry, and diffuse radio emission was taken to confirm that CL J1226.9+3332 is dynamically disturbed. The same study described the center as having a “puffed-up” pressure distribution, consistent with ICM sloshing or an early-stage merger, and suggested possible interplay between intracluster medium and circumgalactic medium around the brightest cluster galaxy [2509.14048].

An earlier NIKA proceeding had compared the pressure profile to the low-$z$ “universal” gNFW template of Arnaud et al. (2010) and reported that CL J1226.9+3332 appears marginally more compact in the core, with its location on the $E(z)^{-2/3}D_A^2Y_{500}$–$M_{500}$ plane consistent, within $1\sigma$, with the local scaling relation. A plausible implication is that increasingly resolved millimeter data shifted the emphasis from global scaling consistency toward explicit core disturbance and non-circularity [1409.1137].

## 6. Galaxy population, structural parameters, and the Fundamental Plane

Optical work on RX J1226.9+3332 places the cluster in a second, complementary context: galaxy structural evolution. Gemini/GMOS spectroscopy targeted $\sim 40$–$50$ objects in the field, and the extended sample comprised 87 objects. HST/ACS imaging used four WFC pointings in F606W and F814W, each with $8\times 500\ \mathrm{s}$ exposures, covering about $3'\times 3'$ per pointing at $0.049''$ pixel scale. Structural parameters were measured for 112 galaxies, of which 54 were spectroscopically confirmed cluster members. Targets were chosen independently of morphology to avoid progenitor bias [0910.0511].

Surface-brightness modeling employed GALFIT with both de Vaucouleurs and Sérsic profiles and a $9''$ drizzled Tiny Tim PSF. For the 54 confirmed members, the half-light radius ranged from $0.03''$ to $1.2''$, corresponding to $0.2$–$9\ \mathrm{kpc}$ for $H_0=70$, $\Omega_m=0.3$, $\Omega_\Lambda=0.7$, with median $\approx 0.30''$ and $1\sigma$ scatter $\approx 0.15''$. The effective surface brightness ranged from $19.0$ to $23.5\ \mathrm{mag\ arcsec^{-2}}$, with median $\approx 21.2\ \mathrm{mag\ arcsec^{-2}}$ and scatter $\approx 0.8\ \mathrm{mag}$. Sérsic indices spanned $0.5$ to $8.0$, with mean $\approx 3.9$, median $\approx 4.0$, rms $\approx 0.8$, and fraction with $n>2$ of about $85\%$. Simulations showed that forcing an $r^{1/4}$ fit on a true Sérsic $n<4$ galaxy gives $\Delta \mathrm{mag}\approx -0.1\pm0.3\ \mathrm{mag}$ and $\Delta \log R_e \approx +0.1\pm0.2\ \mathrm{dex}$, but the Fundamental-Plane combination $\log R_e-\beta \log \langle I\rangle_e$ was affected only at the level $\Delta \mathrm{FP}_{\rm par}\approx -0.02\pm0.10$ [0910.0511].

For the early-type members, the Fundamental Plane was reported as
$$
\log R_e=(1.57\pm0.23)\,\log\sigma-(0.82\pm0.11)\,\log\langle I\rangle_e+(0.20\pm0.02).
$$
The plane remained tight, with rms in $\log R_e$ of $\lesssim 0.11$ dex, but its slope differed from the local Coma relation. The quoted interpretation was “downsizing”: lower-mass galaxies experienced more recent luminosity evolution than the most massive systems. The zero-point offset corresponded to $\Delta \log(M/L_B)\approx -0.54\pm0.03$ at $z=0.89$, consistent with passive evolution of an old stellar population formed at $z\gtrsim 2$. In that sense, the cluster simultaneously hosts a mature early-type galaxy population and an intracluster medium with a disturbed core, illustrating that stellar-population maturity and ICM dynamical complexity need not coincide [0910.0511].

Source: https://www.emergentmind.com/topics/cl-j1226-9-3332