Note on $ r $-central Lah numbers and $ r $-central Lah-Bell numbers

<abstract><p>The $ r $-Lah numbers generalize the Lah numbers to the $ r $-Stirling numbers in the same sense. The Stirling numbers and the central factorial numbers are one of the important tools in enumerative combinatorics. The $ r $-Lah number counts the number of partitions of a set with $ n+r $ elements into $ k+r $ ordered blocks such that $ r $ distinguished elements have to be in distinct ordered blocks. In this paper, the $ r $-central Lah numbers and the $ r $-central Lah-Bell numbers ($ r\in \mathbb{N} $) are introduced parallel to the $ r $-extended central factorial numbers of the second kind and $ r $-extended central Bell polynomials. In addition, some identities related to these numbers including the generating functions, explicit formulas, binomial convolutions are derived. Moreover, the $ r $-central Lah numbers and the $ r $-central Lah-Bell numbers are shown to be represented by Riemann integral, respectively.</p></abstract>

Covariant decomposition of the non-linear galaxy number counts and their monopole

We present a fully non-linear and relativistically covariant expression for the observed galaxy density contrast. Building on a null tetrad tailored to the cosmological observer's past light cone, we find a decomposition of the non-linear galaxy over-density into manifestly gauge-invariant quantities, each of which has a clear physical interpretation as a cosmological observable. This ensures that the monopole of the galaxy over-density field (the mean galaxy density as a function of observed redshift) is properly accounted for. We anticipate that this decomposition will be useful for future work on non-linearities in galaxy number counts, for example, deriving the relativistic expression for the galaxy bispectrum. We then specialise our results to conformal Newtonian gauge, with a Hubble parameter either defined globally or measured locally, illustrating the significance of the different contributions to the observed monopole of the galaxy density.

Cosmological simulations of number counts

In this paper we present for the first time the angular power spectra C ℓ(z,z') for number counts from relativistic N-body simulations. We use the relativistic N-body code gevolution with its exact integration of lightlike geodesics which include all relativistic scalar contributions to the number counts. We compare our non-perturbative numerical results with the results from class using the hmcode approximation for the non-linear matter power spectrum. We find that this simple description is excellent for both, the density and the convergence. On the other hand, the current implementation of redshift-space distortions in Boltzmann codes is not accurate. We also find that the largest contribution to the unequal-redshift power spectra is the cross-correlation of the density and the lensing contribution to the number counts, especially for redshift bins that are far apart. Correlating the number counts with the convergence map we find that the signal is dominated by the lensing-lensing term when the convergence field redshift is not higher than the number counts one, while it is dominated by the density-lensing term in the opposite case. In the present study, the issue of galaxy bias is deliberately left aside by considering only unbiased samples of matter particles from the simulations.

Virial Halo Mass Function in the Planck Cosmology

We study halo mass functions with high-resolution N-body simulations under a ΛCDM cosmology. Our simulations adopt the cosmological model that is consistent with recent measurements of the cosmic microwave backgrounds with the Planck satellite. We calibrate the halo mass functions for 108.5 ≲ M vir/(h −1 M ⊙) ≲ 1015.0–0.45 z , where M vir is the virial spherical-overdensity mass and redshift z ranges from 0 to 7. The halo mass function in our simulations can be fitted by a four-parameter model over a wide range of halo masses and redshifts, while we require some redshift evolution of the fitting parameters. Our new fitting formula of the mass function has a 5%-level precision, except for the highest masses at z ≤ 7. Our model predicts that the analytic prediction in Sheth & Tormen would overestimate the halo abundance at z = 6 with M vir = 108.5–10 h −1 M ⊙ by 20%–30%. Our calibrated halo mass function provides a baseline model to constrain warm dark matter (WDM) by high-z galaxy number counts. We compare a cumulative luminosity function of galaxies at z = 6 with the total halo abundance based on our model and a recently proposed WDM correction. We find that WDM with its mass lighter than 2.71 keV is incompatible with the observed galaxy number density at a 2σ confidence level.

Calibration of bias and scatter involved in cluster mass measurements using optical weak gravitational lensing

Cosmological inference from cluster number counts is systematically limited by the accuracy of the mass calibration, i.e. the empirical determination of the mapping between cluster selection observables and halo mass. In this work we demonstrate a method to quantitatively determine the bias and uncertainties in weak-lensing mass calibration. To this end, we extract a library of projected matter density profiles from hydrodynamical simulations. Accounting for shear bias and noise, photometric redshift uncertainties, mis-centering, cluster member contamination, cluster morphological diversity, and line-of-sight projections, we produce a library of shear profiles. Fitting a one-parameter model to these profiles, we extract the so-called weak lensing mass MWL. Relating the weak-lensing mass to the halo mass from gravity-only simulations with the same initial conditions as the hydrodynamical simulations allows us to estimate the impact of hydrodynamical effects on cluster number counts experiments. Creating new shear libraries for ∼1000 different realizations of the systematics, provides a distribution of the parameters of the weak-lensing to halo mass relation, reflecting their systematic uncertainty. This result can be used as a prior for cosmological inference. We also discuss the impact of the inner fitting radius on the accuracy, and determine the outer fitting radius necessary to exclude the signal from neighboring structures. Our method is currently being applied to different Stage III lensing surveys, and can easily be extended to Stage IV lensing surveys.

Understanding matched filters for precision cosmology

Matched filters are routinely used in cosmology in order to detect galaxy clusters from mm observations through their thermal Sunyaev–Zeldovich (tSZ) signature. In addition, they naturally provide an observable, the detection signal-to-noise or significance, which can be used as a mass proxy in number counts analyses of tSZ-selected cluster samples. In this work, we show that this observable is, in general, non-Gaussian, and that it suffers from a positive bias, which we refer to as optimization bias. Both aspects arise from the fact that the signal-to-noise is constructed through an optimization operation on noisy data, and hold even if the cluster signal is modelled perfectly well, no foregrounds are present, and the noise is Gaussian. After reviewing the general mathematical formalism underlying matched filters, we study the statistics of the signal-to-noise with a set Monte Carlo mock observations, finding it to be well-described by a unit-variance Gaussian for signal-to-noise values of 6 and above, and quantify the magnitude of the optimization bias, for which we give an approximate expression that may be used in practice. We also consider the impact of the bias on the cluster number counts of Planck and the Simons Observatory (SO), finding it to be negligible for the former and potentially significant for the latter.

Cluster assembly times as a cosmological test

The abundance of galaxy clusters in the low-redshift universe provides an important cosmological test, constraining a product of the initial amplitude of fluctuations and the amount by which they have grown since early times. The degeneracy of the test with respect to these two factors remains a limitation of abundance studies. Clusters will have different mean assembly times, however, depending on the relative importance of initial fluctuation amplitude and subsequent growth. Thus, structural probes of cluster age such as concentration, shape or substructure may provide a new cosmological test that breaks the main degeneracy in number counts. We review analytic predictions for how mean assembly time should depend on cosmological parameters, and test these predictions using cosmological simulations. Given the overall sensitivity expected, we estimate the cosmological parameter constraints that could be derived from the cluster catalogues of forthcoming surveys such as Euclid, the Nancy Grace Roman Space Telescope, eROSITA, or CMB-S4. We show that by considering the structural properties of their cluster samples, such surveys could easily achieve errors of Δσ8 = 0.01 or better.