Galaxy Number Counts Research Articles

ASW^2DF: Census of the obscured star formation in a galaxy cluster in formation at z=2.2

We report the results of the deep and wide Atacama Large Millimeter/submillimeter Array (ALMA) 1.2\,mm mapping of the Spiderweb protocluster at $z=2.16$. The observations were divided into six contiguous fields covering a survey area of 19.3\,arcmin$^2$. With sim 13h of on-source time, the final maps in the six fields reach the 1sigma rms noise in a range of $40.3-57.1\ at a spatial resolution of $0 By using different source extraction codes and careful visual inspection, we detected 47 ALMA sources at a significance higher than 4sigma . We constructed the differential and cumulative number counts down to $ after the correction for purity and completeness obtained from Monte Carlo simulations. The ALMA 1.2\,mm number counts of dusty star-forming galaxies (DSFGs) in the Spiderweb protocluster are overall two times that of general fields, with some regions showing even higher overdensities (more than a factor of three). This is consistent with the results from previous studies over a larger scale using single-dish instruments. Comparison of the spatial distributions between different populations indicates that our ALMA sources are likely drawn from the same distribution as CO(1-0) emitters from the COALAS large program but are distinct from that of Halpha emitters. The cosmic star formation rate density of the ALMA sources is consistent with previous results (e.g., LABOCA 870\,mu m observations) after accounting for the difference in volume. We show that molecular gas masses estimates from dust measurements are not consistent with the ones derived from CO(1-0) and thus have to be taken with caution. The multiplicity fraction of single-dish DSFGs is higher than that of the field. Moreover, two extreme concentrations of ALMA sources were found on the outskirts of the Spiderweb protocluster, with an excess of more than 12 times that of the general fields. These results indicate that the ALMA-detected DSFGs are supplied through gas accretion along filaments and are triggered by intense star formation by accretion shocks before falling into the cluster center. The identified two galaxy groups are likely falling into the protocluster center and will trigger new merger events eventually, as indicated in simulations.

Detecting relativistic Doppler by multi-tracing a single galaxy population

New data from ongoing galaxy surveys, such as the Euclid satellite and the Dark Energy Spectroscopic Instrument (DESI), are expected to unveil physics on the largest scales of our universe. Dramatically affected by cosmic variance, these scales are of interest to large-scale structure studies as they exhibit relevant corrections due to general relativity (GR) in the n-point statistics of cosmological random fields. We focus on the relativistic, sample-dependent Doppler contribution to the observed clustering of galaxies, whose detection will further confirm the validity of GR in cosmological regimes. Sample- and scale-dependent, the Doppler term is more likely to be detected via cross-correlation measurements, where it acts as an imaginary correction to the power spectrum of fluctuations in galaxy number counts. We present a method allowing us to exploit multi-tracer benefits from a single data set, by subdividing a galaxy population into two sub-samples, according to galaxies’ luminosity/magnitude. To overcome cosmic variance we rely on a multi-tracer approach, and to maximise the detectability of the relativistic Doppler contribution in the data, we optimise sample selection. As a result, we find the optimal split and forecast the relativistic Doppler detection significance for both a DESI-like Bright Galaxy Sample and a Euclid-like Hα galaxy population.

Neural network reconstruction of density and velocity fields from the 2MASS Redshift Survey

Aims. Our aim is to reconstruct the 3D matter density and peculiar velocity fields in the local Universe up to a distance of 200 h−1 Mpc from the Two-Micron All-Sky Redshift Survey (2MRS) using a neural network (NN). Methods. We employed an NN with a U-net autoencoder architecture and a weighted mean squared error loss function trained separately to output either the density or velocity field for a given input grid of galaxy number counts. The NN was trained on mocks derived from the Quijote N-body simulations, incorporating redshift-space distortions (RSDs), galaxy bias, and selection effects closely mimicking the characteristics of 2MRS. The trained NN was benchmarked against a standard Wiener filter (WF) on a validation set of mocks before applying it to 2MRS. Results. The NN reconstructions effectively approximate the mean posterior estimate of the true density and velocity fields conditioned on the observations. They consistently outperform the WF in terms of reconstruction accuracy and effectively capture the nonlinear relation between velocity and density. The NN-reconstructed bulk flow of the total survey volume exhibits a significant correlation with the true mock bulk flow, demonstrating that the NN is sensitive to information on “super-survey” scales encoded in the RSDs. When applied to 2MRS, the NN successfully recovers the main known clusters, some of which are partially in the Zone of Avoidance. The reconstructed bulk flows in spheres of different radii less than 100 h−1 Mpc are in good agreement with a previous 2MRS analysis that required an additional external bulk flow component inferred from directly observed peculiar velocities. The NN-reconstructed peculiar velocity of the Local Group closely matches the observed Cosmic Microwave Background dipole in amplitude and Galactic latitude, and only deviates by 18° in longitude. The NN-reconstructed fields are publicly available.

5–25 μm Galaxy Number Counts from Deep JWST Data

Galaxy number counts probe the evolution of galaxies over cosmic time and serve as a valuable comparison point to theoretical models of galaxy formation. We present new galaxy number counts in eight photometric bands between 5 and 25 μm from the Systematic Mid-infrared Instrument Legacy Extragalactic Survey and the JWST Advanced Deep Extragalactic Survey deep MIRI parallel, extending to unprecedented depth. By combining our new MIRI counts with existing data from Spitzer and AKARI, we achieve counts across 3–5 orders of magnitude in flux in all MIRI bands. Our counts diverge from predictions from recent semianalytical models of galaxy formation, likely due to their treatment of mid-IR aromatic features. Finally, we integrate our combined JWST−Spitzer counts at 8 and 24 μm to measure the cosmic infrared background (CIB) light at these wavelengths; our measured CIB fluxes are consistent with those from previous mid-IR surveys but larger than predicted by models based on TeV blazar data.

A3COSMOS and A3GOODSS: Continuum source catalogues and multi-band number counts

Context. Galaxy submillimetre number counts are a fundamental measurement in our understanding of galaxy evolution models. Most early measurements are obtained via single-dish telescopes with substantial source confusion, whereas recent interferometric observations are limited to small areas. Aims. We used a large database of ALMA continuum observations to accurately measure galaxy number counts in multiple (sub)millimetre bands, thus bridging the flux density range between single-dish surveys and deep interferometric studies. Methods. We continued the Automated Mining of the ALMA Archive in the COSMOS Field project (A3COSMOS) and extended it with observations from the GOODS-South field (A3GOODSS). The database consists of ∼4000 pipeline-processed continuum images from the public ALMA archive, yielding 2050 unique detected sources, including sources with and without a known optical counterpart. To infer galaxy number counts, we constructed a method to reduce the observational bias inherent to targeted pointings that dominate the database. This method comprises a combination of image selection, masking, and source weighting. The effective area was calculated by accounting for inhomogeneous wavelengths, sensitivities, and resolutions and for the spatial overlap between images. We tested and calibrated our method with simulations. Results. We derived the number counts in a consistent and homogeneous way in four different ALMA bands covering a relatively large area. The results are consistent with number counts retrieved from the literature within the uncertainties. In Band 7, at the depth of the inferred number counts, ∼40% of the cosmic infrared background is resolved into discrete sources. This fraction, however, decreases with increasing wavelength, reaching ∼4% in Band 3. Finally, we used the number counts to test models of dusty galaxy evolution, and find a good agreement within the uncertainties. Conclusions. By continuing the A3COSMOS and A3GOODSS archival effort, we obtained the deepest archive-based (sub)millimetre number counts measured to date over such a wide area. This database proves to be a valuable resource that, thanks to its substantial size, can be used for statistical analyses after having applied certain conservative restrictions.

Impact of property covariance on cluster weak lensing scaling relations

ABSTRACT We present an investigation into a hitherto unexplored systematic that affects the accuracy of galaxy cluster mass estimates with weak gravitational lensing. Specifically, we study the covariance between the weak lensing signal, ΔΣ, and the ‘true’ cluster galaxy number count, Ngal, as measured within a spherical volume that is void of projection effects. By quantifying the impact of this covariance on mass calibration, this work reveals a significant source of systematic uncertainty. Using the MDPL2 simulation with galaxies traced by the SAGE semi-analytic model, we measure the intrinsic property covariance between these observables within the three-dimensional vicinity of the cluster, spanning a range of dynamical mass and redshift values relevant for optical cluster surveys. Our results reveal a negative covariance at small radial scales (R ≲ R200c) and a null covariance at large scales (R ≳ R200c) across most mass and redshift bins. We also find that this covariance results in a $2{\!-\!}3~{{\ \rm per\ cent}}$ bias in the halo mass estimates in most bins. Furthermore, by modelling Ngal and ΔΣ as multi-(log)-linear equations of secondary halo properties, we provide a quantitative explanation for the physical origin of the negative covariance at small scales. Specifically, we demonstrate that the Ngal–ΔΣ covariance can be explained by the secondary properties of haloes that probe their formation history. We attribute the difference between our results and the positive bias seen in other works with (mock)-cluster finders to projection effects. These findings highlight the importance of accounting for the covariance between observables in cluster mass estimation, which is crucial for obtaining accurate constraints on cosmological parameters.

Inferring the dark matter splashback radius from cluster gas and observable profiles in the FLAMINGO simulations

ABSTRACT The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev–Zel’dovich (SZ) observations of the hot intracluster gas. Here, we use the FLAMINGO simulations to investigate whether an extremum gradient in a similar position to the splashback radius is predicted to occur in the cluster gas profiles. We find that the minimum in the gradient of the stacked 3D gas density and pressure profiles, and the maximum in the gradient of the entropy profile, broadly align with the splashback feature though there are significant differences. While the dark matter splashback radius varies with specific mass accretion rate, in agreement with previous work, the radial position of the deepest minimum in the log-slope of the gas density is more sensitive to halo mass. In addition, we show that a similar minimum is also present in projected 2D pseudo-observable profiles: emission measure (X-ray), Compton-y (SZ), and surface mass density (weak lensing). We find that the latter traces the dark matter results reasonably well albeit the minimum occurs at a slightly smaller radius. While results for the gas profiles are largely insensitive to accretion rate and various observable proxies for dynamical state, they do depend on the strength of the feedback processes.

Perturbations of cosmological redshift drift

In this paper we calculate the linear perturbations of the cosmological redshift drift. We show explicitly that our expressions are gauge-invariant and compute the power spectrum of the redshift drift perturbations and its correlations with galaxy number counts within linear perturbation theory. Our findings show that the perturbations are small, and that the peculiar velocity and acceleration terms are dominating and cannot be neglected when modeling the full perturbative expression for the redshift drift. We also find that the cross-correlations with galaxy number count fluctuations might increase the detectability of the effect and can help to separate the perturbative effects from the background cosmological redshift drift signal.

An enhanced abundance of bright galaxies in protocluster candidates at z ∼ 3–5

ABSTRACT We present a protocluster search covering z ∼ 3 to z ∼ 5 based on the combination of the Hyper SuprimeCam Subaru Strategic Programme and the CFHT Large Area U-band Deep Survey. We identify about 30 protocluster candidates per unit redshift over the $\sim 25\, \mathrm{deg^2}$ area of the Deep/UltraDeep layer. Protocluster candidates are selected as regions with a significantly enhanced surface density of dropout galaxies. With this large sample, we characterize the properties of their individual member galaxies. We compare the number counts of dropout galaxies in protocluster candidates with that of coeval field galaxies. Rest-frame ultraviolet (UV) bright galaxies are overabundant in protocluster candidates, a trend seen across the full redshift range studied. We do not find evidence for their spatial distribution within protocluster candidates to be distinct from their fainter counterparts, nor for their UV colour to be different from that of field galaxies with the same brightness. Cosmological simulations predict this bright-end excess, with the main cause being a richer population of massive galaxies, with only a minor contribution from an enhancement in star formation activity (and therefore UV emission) at fixed mass. U-to-K SED modelling of our observed samples supports this interpretation. This environmental differentiation in number counts is already in place at z ∼ 5, with no significant redshift dependence over the range in lookback times probed . These observational results and model predictions suggest that the cosmic clock is ahead in high-density environments.

NIKA2 Cosmological Legacy Survey

Context. Finding and characterizing the heavily obscured galaxies with extreme star formation up to very high redshift is key for constraining the formation of the most massive galaxies in the early Universe. It has been shown that these obscured galaxies are major contributors to the accumulation of stellar mass to z ~ 4. At higher redshift, and despite recent progress, the contribution of dust-obscured galaxies remains poorly known. Aims. Deep surveys in the millimeter domain are necessary in order to probe the dust-obscured galaxies at high redshift. We conducted a large observing program at 1.2 and 2 mm with the NIKA2 camera installed on the IRAM 30m telescope. This NIKA2 Cosmological Legacy Survey (N2CLS) covers two emblematic fields: GOODS-N and COSMOS. We introduce the N2CLS survey and present new 1.2 and 2 mm number counts measurements based on the tiered N2CLS observations (from October 2017 to May 2021) covering 1169 arcmin2. Methods. After a careful data reduction and source extraction, we develop an end-to-end simulation that combines an input sky model with the instrument noise and data reduction pipeline artifacts. This simulation is used to compute the sample purity, flux boosting, pipeline transfer function, completeness, and effective area of the survey (taking into account the non-homogeneous sky coverage). For the input sky model, we used the 117 square degree SIDES simulations, which include galaxy clustering. Our formalism allows us to correct the source number counts to obtain galaxy number counts, the difference between the two being due to resolution effects caused by the blending of several galaxies inside the large beam of single-dish instruments. Results. The N2CLS-May2021 survey is already the deepest and largest ever made at 1.2 and 2 mm. It reaches an average 1σ- noise level of 0.17 and 0.048 mJy on GOODS-N over 159 arcmin2, and 0.46 and 0.14 mJy on COSMOS over 1010 arcmin2, at 1.2 and 2 mm, respectively. For a purity threshold of 80%, we detect 120 and 67 sources in GOODS-N and 195 and 76 sources in COSMOS at 1.2 and 2 mm, respectively. At 1.2 mm, the number counts measurement probes consistently 1.5 orders of magnitude in flux density, covering the full flux density range from previous single-dish surveys and going a factor of 2 deeper into the sub-mJy regime. Our measurement connects the bright single-dish to the deep interferometric number counts. At 2 mm, our measurement matches the depth of the deepest interferometric number counts and extends a factor of 2 above the brightest constraints. After correcting for resolution effects, our results reconcile the single-dish and interferometric number counts, which can be further accurately compared with model predictions. Conclusions. While the observation in GOODS-N have already reached the target depth, we expect the final N2CLS survey to be 1.5 times deeper for COSMOS. Thanks to its volume-complete flux selection, the final N2CLS sample will be an ideal reference for conducting a full characterization of dust-obscured galaxies at high redshift.

Hierarchical Inference of the Lensing Convergence from Photometric Catalogs with Bayesian Graph Neural Networks

We present a Bayesian graph neural network (BGNN) that can estimate the weak lensing convergence (κ) from photometric measurements of galaxies along a given line of sight (LOS). The method is of particular interest in strong gravitational time-delay cosmography (TDC), where characterizing the “external convergence” (κ ext) from the lens environment and LOS is necessary for precise Hubble constant (H 0) inference. Starting from a large-scale simulation with a κ resolution of ∼1′, we introduce fluctuations on galaxy–galaxy lensing scales of ∼1″ and extract random sight lines to train our BGNN. We then evaluate the model on test sets with varying degrees of overlap with the training distribution. For each test set of 1000 sight lines, the BGNN infers the individual κ posteriors, which we combine in a hierarchical Bayesian model to yield constraints on the hyperparameters governing the population. For a test field well sampled by the training set, the BGNN recovers the population mean of κ precisely and without bias (within the 2σ credible interval), resulting in a contribution to the H 0 error budget well under 1%. In the tails of the training set with sparse samples, the BGNN, which can ingest all available information about each sight line, extracts a stronger κ signal compared to a simplified version of the traditional method based on matching galaxy number counts, which is limited by sample variance. Our hierarchical inference pipeline using BGNNs promises to improve the κ ext characterization for precision TDC. The code is available as a public Python package, Node to Joy ⏬.

Dependency of high-mass satellite galaxy abundance on cosmology in Magneticum simulations

Context. Observational studies carried out to calibrate the masses of galaxy clusters often use mass–richness relations to interpret galaxy number counts. Aims. Here, we aim to study the impact of the richness–mass relation modelled with cosmological parameters on mock mass calibrations. Methods. We build a Gaussian process regression emulator of high-mass satellite abundance normalisation and log-slope based on cosmological parameters Ωm, Ωb, σ8, h0, and redshift z. We train our emulator using Magneticum hydrodynamic simulations that span different cosmologies for a given set of feedback scheme parameters. Results. We find that the normalisation depends, albeit weakly, on cosmological parameters, especially on Ωm and Ωb, and that their inclusion in mock observations increases the constraining power of these latter by 10%. On the other hand, the log-slope is ≈1 in every setup, and the emulator does not predict it with significant accuracy. We also show that satellite abundance cosmology dependency differs between full-physics simulations, dark-matter only, and non-radiative simulations. Conclusions. Mass-calibration studies would benefit from modelling of the mass–richness relations with cosmological parameters, especially if the satellite abundance cosmology dependency.

Faint Galaxy Number Counts in the Durham and SDSS Catalogues

Galaxy number counts in the K-, H-, I-, R-, B- and U-bands from the Durham Extragalactic Astronomy and Cosmology catalogue could be well-fitted over their whole range using luminosity function (LF) parameters derived from the SDSS at the bright region and required only modest luminosity evolution with the steepening of the LF slope (α), except for a sudden steep increase in the B-band and a less steep increase in the U-band at faint magnitudes that required a starburst evolutionary model to account for the excess faint number counts. A cosmological model treating Hubble expansion as an Einstein curvature required less correction at faint magnitudes than a standard ΛCDM model, without requiring dark matter or dark energy. Data from DR17 of the SDSS in the g, i, r, u and z bands over two areas of the sky centred on the North Galactic Cap (NGC) and above the South Galactic Cap (SGC), with areas of 5954 and 859 sq. deg., respectively, and a combined count of 622,121 galaxies, were used to construct bright galaxy number counts and galaxy redshift/density plots within the limits of redshift ≤0.4 and mag ≤20. Their comparative densities confirmed an extensive void in the Southern sky with a deficit of 26% out to a redshift z ≤ 0.15. Although not included in the number count data set because of its incompleteness at fainter magnitudes, extending the SDSS redshift-number count survey to fainter and more distant galaxies with redshift ≤ 1.20 showed a secondary peak in the number counts with many QSOs, bright X-ray and radio sources, and evolving irregular galaxies with rapid star formation rates. This sub-population at redshifts of 0.45–0.65 may account for the excess counts observed in the B-band. Recent observations from the HST and James Webb Space Telescope (JWST) have also begun to reveal a high density of massive galaxies at high redshifts (z>7) with high UV and X-ray emissions, and future observations by the JWST may reveal the assembly of galaxies in the early universe going back to the first light in the universe.

The N5K Challenge: Non-Limber Integration for LSST Cosmology

The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angular power spectrum of galaxy number counts and cosmic shear data without invoking the so-called 'Limber approximation', in the context of the Rubin Observatory Legacy Survey of Space and Time (LSST). We quantify the performance, in terms of accuracy and speed, of three non-Limber implementations: ${\tt FKEM (CosmoLike)}$, ${\tt Levin}$, and ${\tt matter}$, themselves based on different integration schemes and approximations. We find that in the challenge's fiducial 3x2pt LSST Year 10 scenario, ${\tt FKEM (CosmoLike)}$ produces the fastest run time within the required accuracy by a considerable margin, positioning it favourably for use in Bayesian parameter inference. This method, however, requires further development and testing to extend its use to certain analysis scenarios, particularly those involving a scale-dependent growth rate. For this and other reasons discussed herein, alternative approaches such as ${\tt matter}$ and ${\tt Levin}$ may be necessary for a full exploration of parameter space. We also find that the usual first-order Limber approximation is insufficiently accurate for LSST Year 10 3x2pt analysis on $\ell=200-1000$, whereas invoking the second-order Limber approximation on these scales (with a full non-Limber method at smaller $\ell$) does suffice.

Velocity reconstruction with the cosmic microwave background and galaxy surveys

The kinetic Sunyaev Zel'dovich (kSZ) and moving lens effects, secondary contributions to the cosmic microwave background (CMB), carry significant cosmological information due to their dependence on the large-scale peculiar velocity field. Previous work identified a promising means of extracting this cosmological information using a set of quadratic estimators for the radial and transverse components of the velocity field. These estimators are based on the statistically anisotropic components of the cross-correlation between the CMB and a tracer of large scale structure, such as a galaxy redshift survey. In this work, we assess the challenges to the program of velocity reconstruction posed by various foregrounds and systematics in the CMB and galaxy surveys, as well as biases in the quadratic estimators. To do so, we further develop the quadratic estimator formalism and implement a numerical code for computing properly correlated spectra for all the components of the CMB (primary/secondary blackbody components and foregrounds) and a photometric redshift survey, with associated redshift errors, to allow for accurate forecasting. We create a simulation framework for generating realizations of properly correlated CMB maps and redshift binned galaxy number counts, assuming the underlying fields are Gaussian, and use this to validate a velocity reconstruction pipeline and assess map-based systematics such as masking. We highlight the most significant challenges for velocity reconstruction, which include biases associated with: modelling errors, characterization of redshift errors, and coarse graining of cosmological fields on our past light cone. Despite these challenges, the outlook for velocity reconstruction is quite optimistic, and we use our reconstruction pipeline to confirm that these techniques will be feasible with near-term CMB experiments and photometric galaxy redshift surveys.

A simultaneous solution to the Hubble tension and observed bulk flow within 250 h−1 Mpc

ABSTRACT The Λ cold dark matter (ΛCDM) standard cosmological model is in severe tension with several cosmological observations. Foremost is the Hubble tension, which exceeds 5σ confidence. Galaxy number counts show the Keenan–Barger–Cowie (KBC) supervoid, a significant underdensity out to 300 Mpc that cannot be reconciled with ΛCDM cosmology. Haslbauer et al. previously showed that a high local Hubble constant arises naturally due to gravitationally driven outflows from the observed KBC supervoid. The main prediction of this model is that peculiar velocities are typically much larger than expected in the ΛCDM framework. This agrees with the recent discovery by Watkins et al. that galaxies in the CosmicFlows-4 catalogue have significantly faster bulk flows than expected in the ΛCDM model on scales of $100-250 \, h^{-1}$ Mpc. The rising bulk flow curve is unexpected in standard cosmology, causing 4.8σ tension at $200 \, h^{-1}$ Mpc. In this work, we determine what the semi-analytic void model of Haslbauer et al. predicts for the bulk flows on these scales. We find qualitative agreement with the observations, especially if our vantage point is chosen to match the observed bulk flow on a scale of $50 \, h^{-1}$ Mpc. This represents a highly non-trivial success of a previously published model that was not constrained by bulk flow measurements, but which was shown to solve the Hubble tension and explain the KBC void consistently with the peculiar velocity of the Local Group. Our results suggest that several cosmological tensions can be simultaneously resolved if structure grows more efficiently than in the ΛCDM paradigm on scales of tens to hundreds of Mpc.

Galaxy number-count dipole and superhorizon fluctuations

In view of the growing tension between the dipole anisotropy of number counts of cosmologically distant sources and of the cosmic microwave background (CMB), we investigate the number count dipole induced by primordial perturbations with wavelength comparable to or exceeding the Hubble radius today. First, we find that neither adiabatic nor isocurvature superhorizon modes can generate an intrinsic number count dipole. However a superhorizon isocurvature mode does induce a relative velocity between the CMB and the (dark) matter rest frames and thereby affects the CMB dipole. We revisit the possibility that it has an intrinsic component due to such a mode, thus enabling consistency with the galaxy number count dipole if the latter is actually kinematic in origin. Although this scenario is not particularly natural, there are possible links with other anomalies and it predicts a concommitant galaxy number count quadrupole which may be measurable in future surveys. We also investigate the number count dipole induced by modes smaller than the Hubble radius, finding that subject to CMB constraints this is too small to reconcile the dipole tension.

Galaxy source counts at 7.7, 10, and 15 μm with the James Webb Space Telescope

ABSTRACT We present mid-infrared (IR) galaxy number counts based on the Early Release Observations obtained by the James Webb Space Telescope (JWST) at 7.7-, 10-, and 15-μm (F770W, F1000W, and F1500W, respectively) bands of the Mid-Infrared Instrument (MIRI). Due to the superior sensitivity of JWST, the 80-per cent completeness limits reach 0.32, 0.79, and 2.0 μJy in F770W, F1000W, and F1500W filters, respectively, i.e. ∼100 times deeper than previous space IR telescopes such as Spitzer or AKARI. The number counts reach much deeper than the broad bump around 0.05∼0.5 mJy due to polycyclic aromatic hydrocarbon (PAH) emissions. An extrapolation towards fainter flux from the evolutionary models in the literature agrees amazingly well with the new data, where the extrapolated faint-end of IR luminosity functions combined with the cosmic star-formation history to higher redshifts can reproduce the deeper number counts by JWST. Our understanding of the faint IR sources has been confirmed by the observed data due to the superb sensitivity of JWST.

Elucidating galaxy assembly bias in SDSS

We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy sample using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as $\mathcal{Q}_\mathrm{cen}$ and $\mathcal{Q}_\mathrm{sat}$, respectively, to predict a suite of one- and two-point observables. In particular, our fiducial constraint employs the probability distribution of the galaxy number counts measured on $8\,\mathrm{Mpc}\,h^{-1}$ scales $N_8^g$ and the projected cross-correlation functions of quintiles of galaxies selected by $N_8^g$ with our entire galaxy sample. We perform extensive tests of the efficacy of our method by fitting the same observables to mock data using both constrained and non-constrained simulations. We discover that in many cases the level of cosmic variance between the two simulations can produce biased constraints that lead to an erroneous detection of galaxy assembly bias if the non-constrained simulation is used. When applying our method to the SDSS data, the ELUCID reconstruction effectively removes an otherwise strong degeneracy between cosmic variance and galaxy assembly bias in SDSS, enabling us to derive an accurate and stringent constraint on the latter. Our fiducial ELUCID constraint, for galaxies above a stellar mass threshold $M_*{=}10^{10.2}\,h^{-2}\,M_\odot$, is $\mathcal{Q}_\mathrm{cen}{=}{-}0.09\pm{0.05}$ and $\mathcal{Q}_\mathrm{sat}{=}0.09\pm{0.10}$, indicating no evidence for a significant~($>2\sigma$) galaxy assembly bias in the local Universe probed by SDSS. Finally, our method provides a promising path to the robust modelling of the galaxy-halo connection within future surveys like DESI and PFS.

Redshift weighted galaxy number counts

In this paper we introduce the 'redshift fluctuation' as a gauge-invariant cosmological observable and give its fully relativistic expression at first order in cosmological perturbation theory. We show that this corresponds effectively to number counts with a radial window function with vanishing mean which therefore resolve smaller scale radial modes than standard number counts. In a detailed analysis of the angular power spectrum of this new variable, we study the relevance of different relativistic contributions, and how it differs from the conventional observable galaxy number count fluctuations. In order to investigate its utility for future spectroscopic surveys, we perform Fisher forecasts for a Euclid-like and an SKAII-like configuration, as examples. Particular focus is placed on the dependence of the results on the size of the redshift bins and on the cutoff in ℓ adopted in the analysis.