Impact of a non-Gaussian density field on Sunyaev-Zeldovich observables

J. Silk; S. Sadeh; Y. Rephaeli

REVIEW

Impact of a non-Gaussian density field on Sunyaev-Zeldovich observables

S. Sadeh, Y. Rephaeli, J. Silk

Astrophysics · Submitted 5 Mar 2006 · 2 references

Pith's verdict: Useful S-Z calculations for a chi^2 mass function, but the large differences vs. Gaussian are not cleanly attributable to non-Gaussianity because the two models also differ in spectral index and transfer function.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv IDastro-ph/0603115 (version 1)

Pith IDRLEXVMOT, Pith's permanent record of this paper and its review

CitationJ. Silk, S. Sadeh, Y. Rephaeli. "Impact of a non-Gaussian density field on Sunyaev-Zeldovich observables." arXiv:astro-ph/0603115 (2006). https://pith.science/paper/RLEXVMOT

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

Two things to know. The paper is the first to compute S-Z power spectrum, number counts, and angular correlation function for a chi^2 distributed density field mass function, extending Mathis et al. (2004) into S-Z observables. The calculations look internally consistent, and the A2PCF derivation, with the Bessel function integral done in closed form, is a nice piece of work. The sensitivity tests to temperature-mass slope and core radius are sensible.

The soft spot is the model comparison. The Gaussian model uses a scale-invariant spectrum (n=1) with an adiabatic transfer function and sigma8=0.9. The chi^2_1 model uses n=-1.8 with an isocurvature transfer function and sigma8=0.73, chosen to match local cluster abundance. So the two models differ in three ways at once: the PDF, the spectral index, and the transfer function. Since sigma(M,z) in equation (2) and the dark-matter correlation in equation (25) both depend directly on P(k), the factor-of-seven-to-twenty differences in S-Z power, and the large differences in counts and clustering, could be driven mostly by the power-spectrum shape rather than by non-Gaussianity. The abstract and Section 5 attribute the effect specifically to the mass function, which overstates what the calculation supports. A matched-P(k) comparison—same spectral index and transfer function, only the mass function functional form changed—would isolate the PDF contribution. Without that, the headline claim is not cleanly demonstrated.

The comparison with CBI/ACBAR is qualitative and lacks error treatment, but that is a minor concern because the paper's main claim is comparative. The sigma8 calibration is a reasonable normalization, not a damaging circularity.

The paper is a legitimate, clearly written comparative study with new computations. It deserves a serious referee, but a referee should push for a cleaner decomposition of the PDF effect versus the power-spectrum effect. I would not cite it for the headline claim; I might cite the A2PCF formalism. Bring it to reading group if you want to discuss confounding variables in model comparison.

Referee report

Significance

If the central claim were established by a controlled comparison, the paper would demonstrate that S-Z power spectra, number counts, and cluster angular clustering are powerful discriminators between Gaussian and positively skewed non-Gaussian density fields, with order-of-magnitude differences that upcoming surveys could test. The analytic framework is standard, the sensitivity tests to temperature-mass normalization and core radius are informative, and the paper makes explicit, falsifiable predictions. However, because the two models differ not only in the PDF of the density field but also in the spectral index, the transfer function, and the σ₈ normalization, the current calculations do not cleanly isolate the effect of the mass function. The reported significance of the claim is therefore not yet supported, although it is likely salvageable with additional matched-power-spectrum runs.

Major objections (3)

[Section 2, Eqs. (8)–(9) and paragraph following Eq. (9)] The two models are not varied solely through the mass function. The Gaussian model uses n = 1 with an adiabatic transfer function (Eq. 8) and σ₈ = 0.9, while the χ²_1 model uses n = −1.8 with an isocurvature transfer function (Eq. 9) and σ₈ = 0.73, the latter chosen to match the local cluster abundance. Since σ(M,z) in Eq. (2) and the mass function in Eq. (1) depend directly on P(k), the computed differences in C_l, number counts, and w(θ) could be driven largely by the change in spectral index and transfer function rather than by the non-Gaussian PDF. The abstract and Section 5 attribute the differences specifically to the mass function, but the calculations as presented do not isolate the PDF contribution. I recommend adding a controlled comparison in which the same P(k) (same n, same transfer function, and either matched σ₈ or matched low-z abundance) is used for both PDFs, and reporting how much of the effect survives.
[Section 3.1, Fig. 6] The comparison with the CBI and ACBAR data is only semi-quantitative, as the authors themselves acknowledge, but the text states that the χ²_1 model 'provides a reasonable match to the observed power excess at ℓ ≳ 2000.' The plotted S-Z power spectrum is only one component of the measured power at these multipoles; primary CMB anisotropy, point sources, and other foregrounds are not modeled. Without a treatment of these components or an error budget for the data points, the statement goes somewhat beyond what the figure can support. The wording should be tempered or the analysis extended to include at least a schematic foreground model.
[Section 5 (formation-time discussion, Fig. 11)] The abstract claims that the formation-time distribution of clusters is 'demonstrated to be sensitive to the underlying mass function,' but the paper presents only a preliminary calculation. The formation-time analysis uses the χ²_2 model rather than the χ²_1 model used for the main observables, adopts an arbitrary definition of formation time (f = 0.5 or 0.75 of the final mass), and does not show the actual PDFs of formation times — only ratios of effective temperatures and concentration parameters. The qualitative claim about earlier formation in the non-Gaussian model is plausible, but the current demonstration is not at the same standard as the S-Z power-spectrum and counts calculations.

Minor objections (5)

[Section 3.1, first two paragraphs] There is an internal inconsistency: the text first says that both the Gaussian and non-Gaussian models 'peak at multipoles higher than ℓ > 10,000,' then later states that the Gaussian case (a) 'peaks at multipole ℓ ∼ 10,000.' Please reconcile these statements.
[Section 3.2, paragraph listing g(x) values] The sentence giving g(x) = −4, 3.2, 6.7 should explicitly state the corresponding frequencies in the same order (143, 353, 545 GHz) rather than relying on the preceding sentence to convey the mapping.
[Section 4.1, Eq. (32)] The notation in the integrand is confusing: 'sink√(...)' should be written as sin(kR)/(kR) with R explicitly the square-root expression, to avoid apparent dimensional inconsistencies.
[Section 3.2, final paragraph] The phrase 'total count levels are significantly higher in the χ²_1 model, roughly by factors of ∼7.5 and ∼3' is followed a few lines later by 'up to an order of magnitude.' These are compatible, but the wording would be clearer if the factor for each case were stated once with the corresponding caveat.
[Section 5, third paragraph] The sentence 'It is the highest-mass clusters that are most susceptible to the high-mass tail of the mass function' seems to mean 'most sensitive to the shape of the high-mass tail'; consider rewording.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

A positively skewed density field changes Sunyaev-Zeldovich power, counts, and clustering by up to an order of magnitude.

The reading

This paper tries to establish that the three principal statistical Sunyaev-Zeldovich observables—the angular power spectrum, cumulative cluster number counts, and the angular correlation function—are highly sensitive to whether the cluster mass function comes from a Gaussian or a positively skewed non-Gaussian density field. It compares the standard Press-Schechter mass function with one derived from a chi-square-distributed field, with both normalized to the same present-day cluster abundance. The non-Gaussian model produces up to an order of magnitude more S-Z power and counts, with the extra signal coming from massive clusters at high redshift. If the comparison is sound, S-Z surveys can act as a statistical probe of non-Gaussian initial conditions, complementary to CMB-based constraints.

What carries the argument

The load-bearing object is the cluster mass function. For a Gaussian field the Press-Schechter multiplicity is $F(\mu)=\sqrt{2/\pi}\,e^{-\mu^2/2}\,(\mu/\sigma_R)$, while for the $\chi^2_m$ model it takes the form in equation (7), where $m$ counts CDM fields added in quadratures and $m=1$ gives maximal positive skewness. This $F(\mu)$ enters $n(M,z)=-F(\mu)\rho_b/(M\sigma)\,d\sigma/dM$, and the same $n(M,z)$ is fed into the integrals for the S-Z power spectrum, the number counts, and the angular correlation function. The power-spectrum integral uses $G_\ell(M,z)$, the angular transform of the cluster Comptonization profile; the count integral uses the flux above 30 mJy; the correlation-function calculation multiplies the dark-matter correlation by a mass-dependent bias. Cluster gas properties enter through assumed scalings of temperature, density, and core radius with mass and redshift; the formation-time result comes from applying the excursion-set formalism to the same mass function.

Extended reading notes

Core claim

The central claim is that the shape of the density field's probability distribution, not just its power spectrum, controls the high-mass end of the cluster mass function and, through it, all three S-Z observables. With the chi-square-distributed field normalized to the same local cluster abundance as the Gaussian field, the predicted S-Z power spectrum is higher by roughly an order of magnitude and peaks at higher multipoles because high-redshift massive clusters that are rare in the Gaussian case become common. Flux-limited counts above 30 mJy are larger by factors of about 7.5 for an evolving gas temperature and about 3 for no temperature evolution, and the angular correlation function becomes steeper, with higher correlation at arcminute separations and lower correlation at large separations. The paper further claims that cluster formation times are generally earlier in the non-Gaussian model, which would raise the average concentration parameters and intracluster gas temperatures of low-mass clusters.

Load-bearing premise

The comparison treats the density field's probability distribution as the only relevant difference, even though the Gaussian case also has spectral index n = 1 with an adiabatic transfer function while the chi-squared case has n = -1.8 with an isocurvature transfer function and a renormalized sigma8; those differences, rather than the PDF shape, could be producing the S-Z differences.

Editorial extensions

If this is right

The S-Z power spectrum in the non-Gaussian m = 1 model exceeds the Gaussian prediction by a factor that climbs from about 7 to about 20 across multipoles 1000-10000, and the peak shifts to higher multipoles.
Flux-limited cluster counts above 30 mJy are about 7.5 times higher when cluster temperatures evolve and about 3 times higher when they do not, with the non-Gaussian counts drawing much more from massive, distant clusters.
The angular correlation function of S-Z clusters is steeper in the non-Gaussian model: about 2-4 times higher at 1-arcminute separations but lower at large separations, giving a distinctive clustering signature.
Cluster formation times are generally earlier for the skewed field, which the paper argues raises average concentration parameters and intracluster gas temperatures for low-mass clusters.
A less skewed field (m = 2) lowers the S-Z power by a factor of about 1.5-1.8 and reduces high-redshift counts strongly, so the observables can also constrain the degree of skewness.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A cleaner control would recompute the S-Z observables for the chi-squared PDF while keeping the Gaussian model's spectral index and adiabatic transfer function; if the large differences persist, they are really attributable to PDF shape rather than to the differing input power spectra.
The earlier-formation prediction implies non-Gaussianity should also show up in cluster mass-concentration relations and intracluster gas entropy, so internal cluster properties offer an independent, non-CMB check.
Because the non-Gaussian signal is concentrated in high-redshift massive clusters, a survey that can resolve clusters by redshift could separate the mass-function effect from the temperature-scaling uncertainty that dominates the total power.

Original abstract

The main statistical properties of the Sunyaev-Zeldovich (S-Z) effect - the power spectrum, cluster number counts, and angular correlation function - are calculated and compared within the framework of two density fields which differ in their predictions of the cluster mass function at high redshifts. We do so for the usual Press and Schechter mass function, which is derived on the basis of a Gaussian density fluctuation field, and for a mass function based on a chi^2 distributed density field. These three S-Z observables are found to be very significantly dependent on the choice of the mass function. The different predictions of the Gaussian and non-Gaussian density fields are probed in detail by investigating the behaviour of the three S-Z observables in terms of cluster mass and redshift. The formation time distribution of clusters is also demonstrated to be sensitive to the underlying mass function. A semi-quantitative assessment is given of its impact on the concentration parameter and the temperature of intracluster gas.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,000 input tokens.

What the abstract itself states, as probabilities: claims first 0.02 · states limitation 0.14 · human participants 0.01 · controlled comparison 0.63 · randomized assignment 0.01 · states quantitative result 0.08 · compares with prior results 0.77 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

The main statistical properties of the Sunyaev-Zeldovich (S-Z) effect - the power spectrum, cluster number counts, and angular correlation function - are calculated and compared within the framework of two density fields which differ in their predictions of the cluster mass function at high redshifts. We do so for the usual Press and Schechter mass function, which is derived on the basis of a Gaussian density fluctuation field, and for a mass function based on a chi^2 distributed density field. These three S-Z observables are found to be very significantly dependent on the choice of the mass function. The different predictions of the Gaussian and non-Gaussian density fields are probed in detail by investigating the behaviour of the three S-Z observables in terms of cluster mass and redshift. The formation time distribution of clusters is also demonstrated to be sensitive to the underlying mass function. A semi-quantitative assessment is given of its impact on the concentration parameter and the temperature of intracluster gas.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	The main statistical properties of the Sunyaev-Zeldovich (S-Z) effect - the power spectrum, cluster number counts, and angular correlation function - are calculated and compared within the framework of two density fields which differ in their predictions of the cluster mass function at high redshifts.	contribution 0.98	0.93	0.14	asserted 0.99
s1	We do so for the usual Press and Schechter mass function, which is derived on the basis of a Gaussian density fluctuation field, and for a mass function based on a chi^2 distributed density field.	method 0.93	0.61	0.03	asserted 0.94
s2	These three S-Z observables are found to be very significantly dependent on the choice of the mass function.	result 1.00	0.94	0.76	asserted 1.00
s3	The different predictions of the Gaussian and non-Gaussian density fields are probed in detail by investigating the behaviour of the three S-Z observables in terms of cluster mass and redshift.	method 0.57	0.89	0.12	asserted 0.91
s4	The formation time distribution of clusters is also demonstrated to be sensitive to the underlying mass function.	result 0.94	0.90	0.08	asserted 1.00
s5	A semi-quantitative assessment is given of its impact on the concentration parameter and the temperature of intracluster gas.	result 0.62	0.87	0.07	asserted 0.72

What kind of paper the classifier reports this is: theoretical result 0.84 · evidence from simulation or computation 0.69.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures

**Figure 1.** Angular power spectrum of the S-Z effect. Shown are levels of power contributed by clusters lying at redshift range of 0 < z < 6, and at 6 redshift intervals. Left and right-hand panels correspond to the Gaussian and χ 2 1 models, respectively. Upper and lower panels relate to cases (a) and (b), respectively . from clusters lying in the 1015h −1M⊙ < M < 1016h −1M⊙ mass range. Whereas it is practically negligible at all multipoles in the Gaussian model, this is not the case in the non-Gaussian model; in fact, these massive clusters dominate the total power up to ℓ ∼ 1000 and ℓ ∼ 2000…

**Figure 2.** The same as in the previous figure, but here power spectrum levels are calculated for 4 mass intervals, spanning the range 1013h−1M⊙ < M < 1016h−1M⊙. 102 103 104 101 102 l PSχ 2 1 / PSGaussian mode a model b

**Figure 3.** A comparison of S-Z power spectrum levels in the χ 2 1 and Gaussian models for cases (a) and (b). c 2006 RAS, MNRAS 000, 1–17

**Figure 4.** The impact of the temperature-mass relation on the S-Z power spectrum, T ∝ (M/M15) ξ . Shown are 3 cases: ξ = 0.7(3), ξ = 0.6(6), ξ = 0.6. The arrangement of the panels is the same as in the previous figures. Gaussian and non-Gaussian models in slightly different ways; we have already seen that in the Gaussian model the largest contribution to the power spectrum comes from clusters with masses 1014h −1M⊙ < M < 1015h −1M⊙. That is, an outcome of increasing α is that most of the contributing clusters undergo cooling, and the level of power decreases. Obviously, lowering α has the…

**Figure 5.** The impact of cluster core scaling on the S-Z power spectrum. Each of the three values considered (rc = 0.1, 0.125, 0.15 M pc h−1 ) denote the core radius of a cluster with mass 1015M⊙ h−1 at redshift z = 0. It is useful to mention some recent studies of S-Z power spectrum calculations and compare with our results. Case (a) of the Gaussian model is the most adequate for comparison, since past work has focused on Gaussian density fields, and whether analytically or inferred from hydrodynamical simulations, the IC gas temperature has been assumed (or observed) to increase with…

**Figure 6.** Total angular power spectrum of the S-Z effect as calculated above, together with results from the CBI and ACBAR experiments. largely be attributed to uncertainties in the theoretical modeling of IC gas and the techniques implemented in hydrodynamical calculations. To complete this section, we combine, in a single plot, the total power spectrum levels of the S-Z effect in the Gaussian and χ 2 1 models, cases (a) and (b), with results from the CBI and ACBAR experiments. The analytic results were rescaled so as to correspond to an observation frequency of 31 GHz, corresponding to the…

**Figure 7.** Cumulative S-Z counts as a function of redshift; shown are contributions from 4 mass intervals. Note that the contribution from the lowest mass interval (1013 < M < 1014 solar masses) vanishes in all cases considered. the) S-Z spectral function - which is essentially independent of the gas temperature - is g(x) (eq. 19). At these frequencies, g(x) = −4, 3.2, 6.7, respectively, so our results for ν = 353 GHz can be easily scaled accordingly. Adopting the same value of the limiting flux (30mJy), we find that values of the ratios of number counts at ν = 353 GHz and ν = 545 GHz are…

**Figure 8.** Cumulative S-Z counts as a function of redshift. Number counts are shown for the three specified frequency channels of the HFI instrument on the Planck satellite. limit similar to ours, 30 mJy. This translates into all-sky counts of ∼ 2.5 ∼ 104 clusters, whereas we find ∼ 4 · 103 . Holder et al. (2000) report a prediction of ∼ 25 clusters per square degree, or a total of over 106 clusters. However, it is not clear what was the flux limit they assumed. In addition, the slightly higher normalization they adopted, σ8 = 1, may be partly responsible for the inferred high counts. Benson,…

**Figure 9.** A comparison between S-Z power levels (left-hand panel) and number counts (right-hand panel) for a χ 2m-distributed density field with m = 1 and m = 2 is shown. Calculations were carried out for case (a). m = 1 case are higher by a factor of ∼ 1.5 − 1.8 across the multipole range of 100 6 ℓ 6 10000. It is the highest-mass clusters that are most susceptible to the high-mass tail of the mass function, which is strongly dependent on the parameter m. In particular, this is most noticeable at high redshifts where the predictions of the mass function for the massive clusters in the m = 1…

**Figure 10.** The A2PCF calculated for the Gaussian and χ 2 1 models for cases (a) and (b). The left-hand panel describes the correlation levels, whereas the right-hand panel depicts the ratios between the χ 2 1 and Gaussian models. where the last expression was obtained by changing the order of integration. The integral over u can be performed using the following change of variables: x ≡ q dr dz · u 2 + (dA(z) · θ) 2 dA(z) · θ , (33) which results in an integral that can be put in a closed form 2 k(dr/dz) Z ∞ 1 sin [k x dA(z)θ] dx √ x 2 − 1 = π k(dr/dz) J0[kdA(z)θ], (34) where x was defined…

**Figure 11.** Ratios of effective temperatures are shown. The average temperature predicted in the χ 2 2 model is divided by the corresponding quantity in the Gaussian model for f = 0.5 and f = 0.75. Also shown is the ratio between the effective temperature in the non-Gaussian model and the temperature predicted by the temperature-mass relation provided in equation (13). approach, the relevant redshift is the observation redshift. Clearly, this picture does not address the merging process that gives rise to the clusters we observe today. In the merger picture the formation redshift, zf , may…

Figures reproduced from arXiv:astro-ph/0603115 by the authors.

References

Works cited (2)

[1]

6(6), ξ = 0. 6. The arrangement of the panels is the same as in the previous ﬁgures. Gaussian and non-Gaussian models in slightly diﬀerent ways ; we have already seen that in the Gaussian model the largest contribution to the power spectrum comes from clusters with masses 10 14h−1M⊙ <M < 1015h−1M⊙. That is, an outcome of increasing α is that most of the c...

View on Pith 2000
[2]

1, 0. 125, 0. 15 M pc h−1) denote the core radius of a cluster with mass 10 15M⊙ h−1 at redshift z = 0. It is useful to mention some recent studies of S-Z power spect rum calculations and compare with our results. Case (a) of the Gaussian model is the most adequate for comparison, si nce past work has focused on Gaussian density ﬁelds, and whether analyti...

View on Pith Pith review arXiv 2001

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.

[1] [1]

6(6), ξ = 0. 6. The arrangement of the panels is the same as in the previous ﬁgures. Gaussian and non-Gaussian models in slightly diﬀerent ways ; we have already seen that in the Gaussian model the largest contribution to the power spectrum comes from clusters with masses 10 14h−1M⊙ <M < 1015h−1M⊙. That is, an outcome of increasing α is that most of the c...

View on Pith 2000

[2] [2]

1, 0. 125, 0. 15 M pc h−1) denote the core radius of a cluster with mass 10 15M⊙ h−1 at redshift z = 0. It is useful to mention some recent studies of S-Z power spect rum calculations and compare with our results. Case (a) of the Gaussian model is the most adequate for comparison, si nce past work has focused on Gaussian density ﬁelds, and whether analyti...

View on Pith Pith review arXiv 2001