Galaxy Counts Research Articles

A framework for the next generation of stationary cosmological models

According to a new tired-light cosmological model, where [Formula: see text], the number density of galaxies has been nearly constant over the last 10[Formula: see text]Gyr, at least, meaning that, as far as galaxy counts are concerned, the Universe has been stationary. In this context, an analysis of the luminosity distances of quasars and supernovae Ia shows that the Universe is far from being as transparent as assumed nowadays, the photon lifetime along the line-of-sight being one third of the Hubble time. Such a “low” value could mean that there are huge amounts of grey dust in the inter galactic medium, that have so far escaped detection. It could also be a signature of “new physics,” namely, a consequence of the decay of photons into lighter particles. The tired-light model advocated in this study would be falsified if, for instance, the time-dilation of remote events was shown to have a general character, that is, if it were definitely observed for phenomena other than the light curves of supernovae Ia. On the other hand, further developments are needed in order to turn this model into an actual cosmology. In particular, the physical origin of the current stability of the Universe on galactic scales remains to be identified.

The Universe is Brighter in the Direction of Our Motion: Galaxy Counts and Fluxes are Consistent with the CMB Dipole

An observer moving with respect to the cosmic rest frame should observe a concentration and brightening of galaxies in the direction of motion and a spreading and dimming in the opposite direction. The velocity inferred from this dipole should match that of the cosmic microwave background (CMB) temperature dipole if galaxies are on average at rest with respect to the CMB rest frame. However, recent studies have claimed a many-fold enhancement of galaxy counts and flux in the direction of the solar motion compared to the CMB expectation, calling into question the standard cosmology. Here we show that the sky distribution and brightness of extragalactic radio sources are consistent with the CMB dipole in direction and velocity. We use the first epoch of the Very Large Array Sky Survey combined with the Rapid Australian Square Kilometer Array Pathfinder Continuum Survey to estimate the dipole via several different methods, and all show similar results. Typical fits find a km s−1 velocity dipole with apex in Galactic coordinates from source counts and km s−1 toward from radio fluxes. These are consistent with the CMB-solar velocity, 370 km s−1 toward (ℓ, b) = (264, 48), and show that galaxies are on average at rest with respect to the rest frame of the early universe, as predicted by the canonical cosmology.

Seeing-Sorted Large Binocular Camera U-band Imaging of the Extended Groth Strip

Abstract We present the results of deep ground-based U-band imaging of the Extended Groth Strip (EGS) field, and the creation of a best resolution mosaic and a best depth mosaic following the seeing-sorted stacking method of Ashcraft et al. The analysis uses 324 images taken from 2012 March to 2021 June using the Large Binocular Camera. The best resolution mosaic includes the 16.4% of images with FWHM ≤ 0.″9, while the best depth mosaic includes the 73.5% of images with FWHM ≤ 1.″6. Each image’s zero-point was matched to SDSS u, to account for transparency variations. Images with stellar ellipticity (1 − b/a) ≥ 0.08 were omitted. We find that the galaxy counts turn over at U AB ∼ 25.75 mag in the best resolution mosaic and at ∼26 mag in the best depth mosaic. The EGS field is the final UVCANDELS field analyzed using the seeing-sorted stacking method.

A Parameterized Model for Differential Galaxy Counts at Any Wavelength

Abstract Certain astronomy problems, such as the probabilistic cross-matching algorithms of Wilson & Naylor, require descriptions of galaxy counts as a function of brightness in many different filter sets, for which we often lack dedicated observations. With this motivation in mind, I present a simple model for calculating differential galaxy counts for an arbitrary bandpass or wavelength. The model uses a double Schechter function to describe the number density of galaxies as a function of absolute magnitude, including redshift-driven galaxy evolution, for both “blue” and “red” galaxy types. Collating previously derived values for the parameterization of the evolving Schechter function— M 0 * , ϕ 0 * , P, Q, and α—I fit the dependencies of these parameters as a function of wavelength, allowing for the tabulation of differential galaxy counts in a new filter for which observations are unavailable. Finally, I briefly compare the model with some differential galaxy counts obtained from the literature.

The local hole: a galaxy underdensity covering 90 per cent of sky to ≈200 Mpc

ABSTRACT We investigate the ‘Local Hole’, an anomalous underdensity in the local galaxy environment, by extending our previous galaxy K-band number-redshift and number-magnitude counts to ≈90 per cent of the sky. Our redshift samples are taken from the 2MASS Redshift Survey (2MRS) and the 2M++ catalogues, limited to K < 11.5. We find that both surveys are in good agreement, showing an $\approx 21\!-\!22{{\ \rm per\ cent}}$ underdensity at z < 0.075 when compared to our homogeneous counts model that assumes the same luminosity function (LF) and other parameters as in our earlier papers. Using the Two Micron All Sky Survey (2MASS) for n(K) galaxy counts, we measure an underdensity relative to this model of $20 \pm 2 {{\ \rm per\ cent}}$ at K < 11.5, which is consistent in both form and scale with the observed n(z) underdensity. To examine further the accuracy of the counts model, we compare its prediction for the fainter n(K) counts of the Galaxy and Mass Assembly (GAMA) survey. We further compare these data with a model assuming the parameters of a previous study where little evidence for the Local Hole was found. At 13 < K < 16, we find a significantly better fit for our galaxy counts model, arguing for our higher LF normalization. Although our implied underdensity of $\approx 20{{\ \rm per\ cent}}$ means local measurements of the Hubble Constant have been overestimated by ≈3 per cent, such a scale of underdensity is in tension with a global ΛCDM cosmology at an ≈3σ level.

The effect of selection – a tale of cluster mass measurement bias induced by correlation and projection

ABSTRACT Cosmology analyses using galaxy clusters by the Dark Energy Survey have recently uncovered an issue of a previously unknown selection effect affecting weak-lensing mass estimates. In this letter, we use the Illustris-TNG simulation to demonstrate that selecting on galaxy counts induces a selection effect because of projection and correlation between different observables. We compute the weak-lensing-like projected mass estimations of dark matter haloes and examine their projected subhalo counts. In the two-dimensional (2D) projected space, haloes that are measured as more massive than truth have higher subhalo counts. Thus, projection along the line of sight creates cluster observables that are correlated with cluster mass measurement deviations, which in turn creates a mass measurement bias when the clusters are selected by this correlated observable. We demonstrate that the bias is predicted in a forward model using the observable-mass measurement correlation.

Optimal machine-driven acquisition of future cosmological data

We present a set of maps classifying regions of the sky according to their information gain potential as quantified by Fisher information. These maps can guide the optimal retrieval of relevant physical information with targeted cosmological searches. Specifically, we calculated the response of observed cosmic structures to perturbative changes in the cosmological model and we charted their respective contributions to Fisher information. Our physical forward-modeling machinery transcends the limitations of contemporary analyses based on statistical summaries to yield detailed characterizations of individual 3D structures. We demonstrate this advantage using galaxy counts data and we showcase the potential of our approach by studying the information gain of the Coma cluster. We find that regions in the vicinity of the filaments and cluster core, where mass accretion ensues from gravitational infall, are the most informative with regard to our physical model of structure formation in the Universe. Hence, collecting data in those regions would be most optimal for testing our model predictions. The results presented in this work are the first of their kind to elucidate the inhomogeneous distribution of cosmological information in the Universe. This study paves a new way forward for the performance of efficient targeted searches for the fundamental physics of the Universe, where search strategies are progressively refined with new cosmological data sets within an active learning framework.

Galaxy velocity bias in cosmological simulations: towards per cent-level calibration

ABSTRACT Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, bv – defined as the ratio of galaxy and DM velocity dispersions – as a function of redshift, host halo mass, and galaxy stellar mass threshold ($M_{\rm \star , sat}$), for massive haloes ($M_{\rm 200c}\gt 10^{13.5} \, {\rm M}_\odot$) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of bv increasing with M200c and decreasing with redshift and $M_{\rm \star , sat}$. The ensemble-estimated theoretical uncertainty in bv is 2–3 per cent, but becomes percent-level when considering only the three highest resolution simulations. We update the mass–richness normalization for an SDSS redMaPPer cluster sample, and find our improved bv estimates reduce the normalization uncertainty from 22 to 8 per cent, demonstrating that dynamical mass estimation is competitive with weak lensing mass estimation. We discuss necessary steps for further improving this precision. Our estimates for $b_v(M_{\rm 200c}, M_{\rm \star , sat}, z)$ are made publicly available.

Finding High-redshift Galaxies with JWST

Abstract One of the primary goals for the upcoming James Webb Space Telescope is to observe the first galaxies. Predictions for planned and proposed surveys have typically focused on average galaxy counts, assuming a random distribution of galaxies across the observed field. The first and most-massive galaxies, however, are expected to be tightly clustered, an effect known as cosmic variance. We show that cosmic variance is likely to be the dominant contribution to uncertainty for high-redshift mass and luminosity functions, and that median high-redshift and high-mass galaxy counts for planned observations lie significantly below average counts. Several different strategies are considered for improving our understanding of the first galaxies, including adding depth, area, and independent pointings. Adding independent pointings is shown to be the most efficient both for discovering the single highest-redshift galaxy and also for constraining mass and luminosity functions.

The mass and galaxy distribution around SZ-selected clusters

ABSTRACT We present measurements of the radial profiles of the mass and galaxy number density around Sunyaev–Zel’dovich (SZ)-selected clusters using both weak lensing and galaxy counts. The clusters are selected from the Atacama Cosmology Telescope Data Release 5 and the galaxies from the Dark Energy Survey Year 3 data set. With signal-to-noise ratio of 62 (45) for galaxy (weak lensing) profiles over scales of about 0.2–20 h−1 Mpc, these are the highest precision measurements for SZ-selected clusters to date. Because SZ selection closely approximates mass selection, these measurements enable several tests of theoretical models of the mass and light distribution around clusters. Our main findings are: (1) The splashback feature is detected at a consistent location in both the mass and galaxy profiles and its location is consistent with predictions of cold dark matter N-body simulations. (2) The full mass profile is also consistent with the simulations. (3) The shapes of the galaxy and lensing profiles are remarkably similar for our sample over the entire range of scales, from well inside the cluster halo to the quasilinear regime. We measure the dependence of the profile shapes on the galaxy sample, redshift, and cluster mass. We extend the Diemer & Kravtsov model for the cluster profiles to the linear regime using perturbation theory and show that it provides a good match to the measured profiles. We also compare the measured profiles to predictions of the standard halo model and simulations that include hydrodynamics. Applications of these results to cluster mass estimation, cosmology, and astrophysics are discussed.

Obtaining nonlinear galaxy bias constraints from galaxy-lensing phase differences

ABSTRACT We demonstrate the utility and constraining power of a new statistic for investigating galaxy bias: the galaxy-lensing phase difference. The statistic consists in taking the differences of the phases of the harmonic wave modes between the weak lensing convergence field and the galaxy count field. We use dark matter simulations populated with galaxies up to redshift z = 1 to test the performance of this estimator. We find that phase differences are sensitive to the absolute value of the second-order bias (c2 = b2/b1) and demonstrate why this is the case. For a $1500$ deg2 galaxy survey we recover c2 with an error of approximately 0.1 for a wide range of c2 values; current constraints from redshift surveys have errors of 0.1–0.6 depending on redshift. This new statistic is therefore expected to provide constraints for c2 which are complementary and competitive with constraining power by the conventional estimators from the power spectrum and bispectrum. For the Dark Energy Survey (DES), we predict leading measurements of second-order bias.

Magnification bias in galaxy surveys with complex sample selection functions

ABSTRACT Gravitational lensing magnification modifies the observed spatial distribution of galaxies and can severely bias cosmological probes of large-scale structure if not accurately modelled. Standard approaches to modelling this magnification bias may not be applicable in practice as many galaxy samples have complex, often implicit, selection functions. We propose and test a procedure to quantify the magnification bias induced in clustering and galaxy–galaxy lensing (GGL) signals in galaxy samples subject to a selection function beyond a simple flux limit. The method employs realistic mock data to calibrate an effective luminosity function slope, αobs, from observed galaxy counts that can then be used with the standard formalism. We demonstrate this method for two galaxy samples derived from the Baryon Oscillation Spectroscopic Survey (BOSS) in the redshift ranges 0.2 < z ≤ 0.5 and 0.5 < z ≤ 0.75, complemented by mock data built from the MICE2 simulation. We obtain αobs = 1.93 ± 0.05 and αobs = 2.62 ± 0.28 for the two BOSS samples. For BOSS-like lenses, we forecast a contribution of the magnification bias to the GGL signal between the multipole moments, ℓ, of 100 and 4600 with a cumulative signal-to-noise ratio between 0.1 and 1.1 for sources from the Kilo-Degree Survey (KiDS), between 0.4 and 2.0 for sources from the Hyper Suprime-Cam survey (HSC), and between 0.3 and 2.8 for ESA Euclid-like source samples. These contributions are significant enough to require explicit modelling in future analyses of these and similar surveys. Our code is publicly available within the MagBEt module (https://github.com/mwiet/MAGBET).

On the importance of lensing for galaxy clustering in photometric and spectroscopic surveys

We study the importance of gravitational lensing in modelling the number counts of galaxies for the first time in spectroscopic surveys. We confirm previous results for photometric surveys, showing that lensing cannot be neglected in a survey like LSST since it would infer a significant shift of standard cosmological parameters. For a spectroscopic survey like SKA2, we find that neglecting lensing in the monopole, quadrupole and hexadecapole of the correlation function can also induce an important shift of cosmological parameters. For ΛCDM parameters, the shift is moderate, of the order of 0.6σ or less. However, for a model-independent analysis, that measures the growth rate of structure in each redshift bins, neglecting lensing introduces a shift of up to 2.3σ at high redshift. Since the growth rate is directly used to test the theory of gravity, such a strong shift would wrongly be interpreted as the breakdown of General Relativity. This shows the importance of including lensing in the analysis of future surveys. For a survey like DESI, we find on the other hand that lensing is not important, mainly due to the value of the magnification bias parameter of DESI, s(z), which strongly reduces the lensing contribution at high redshift. This result relies on our theoretical modelling of s(z) in DESI and should therefore be confirmed with measurements of s(z) in simulations.We also propose a way of improving the analysis of spectroscopic surveys, by including the cross-correlations between different redshift bins (which is neglected in spectroscopic surveys) from the spectroscopic survey or from a different photometric sample. We show that including the cross-correlations in the SKA2 analysis does not improve the constraints. On the other hand replacing the cross-correlations from SKA2 by cross-correlations measured with LSST improves the constraints by 10% to 20%.The marginal improvement is mainly due to the density correlations between nearby bins and, therefore, does notstrongly depend on our knowledge of the magnification bias. Interestingly, for standard cosmological parameter estimation, the photometric survey LSST in its 12 redshift bin configurationand the spectroscopic SKA2 survey are highly complementary, since they are affected by different degeneracies between parameters:LSST yields the tightest constraints on Ωcdm, h and n s, while SKA2 better constrains Ωbaryon, A s and the bias.

Blending and obscuration in weak-lensing magnification

ABSTRACT We test the impact of some systematic errors in weak-lensing (WL) magnification measurements with the COSMOS 30-band photo-z Survey flux limited to Iauto < 25.0 using correlations of both source galaxy counts and magnitudes. Systematic obscuration effects are measured by comparing counts and magnification correlations. We use the Advanced Camera for Survey-Hubble Space Telescope catalogues to identify potential blending objects (close pairs) and perform the magnification analyses with and without blended objects. We find that blending effects start to be important (∼0.04-mag obscuration) at angular scales smaller than 0.1 arcmin. Extinction and other systematic obscuration effects can be as large as 0.10 mag (U band) but are typically smaller than 0.02 mag depending on the band. After applying these corrections, we measure a 3.9σ magnification signal that is consistent for both counts and magnitudes. The corresponding projected mass profiles of galaxies at redshift z ≃ 0.6 (MI ≃ −21) is Σ = 25 ± 6 M⊙ h3 pc–2 at 0.1 Mpc h–1, consistent with NFW type profile with M200 ≃ 2 × 1012 M⊙h pc–2. Tangential shear and flux-size magnification over the same lenses show similar mass profiles. We conclude that magnification from counts and fluxes using photometric redshifts has the potential to provide complementary WL information in future wide-field surveys once we carefully take into account systematic effects, such as obscuration and blending.

Reproducing submillimetre galaxy number counts with cosmological hydrodynamic simulations

ABSTRACT Matching the number counts of high-z submillimetre-selected galaxies (SMGs) has been a long-standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the submm emission from galaxies in the simba cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850 μm-selected sources. We find good agreement with the shape of the integrated 850 μm luminosity function, and the normalization is within 0.25 dex at >3 mJy, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by simba’s good match to infrared measures of the star formation rate (SFR) function between z = 2 and 4 at high SFRs. Also important is the self-consistent on-the-fly dust model in simba, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a light-cone to investigate the effect of far-field blending, and find that 52 per cent of sources are blends of multiple components, which makes a small contribution to the normalization of the bright end of the number counts. We provide new fits to the 850 μm luminosity as a function of SFR and dust mass. Our results demonstrate that solutions to the discrepancy between submm counts in simulations and observations, such as a top-heavy initial mass function, are unnecessary, and that submillimetre-bright phases are a natural consequence of massive galaxy evolution.

New Horizons Observations of the Cosmic Optical Background

Abstract We used existing data from the New Horizons Long-range Reconnaissance Imager (LORRI) to measure the optical-band (0.4 ≲ λ ≲ 0.9 μm) sky brightness within seven high–Galactic latitude fields. The average raw level measured while New Horizons was 42–45 au from the Sun is 33.2 ± 0.5 nW m−2 sr−1. This is ∼10× as dark as the darkest sky accessible to the Hubble Space Telescope, highlighting the utility of New Horizons for detecting the cosmic optical background (COB). Isolating the COB contribution to the raw total required subtracting scattered light from bright stars and galaxies, faint stars below the photometric detection limit within the fields, and diffuse Milky Way light scattered by infrared cirrus. We removed newly identified residual zodiacal light from the IRIS 100 μm all-sky maps to generate two different estimates for the diffuse Galactic light. Using these yielded a highly significant detection of the COB in the range 15.9 ± 4.2 (1.8 stat., 3.7 sys.) nW m−2 sr−1 to 18.7 ± 3.8 (1.8 stat., 3.3 sys.) nW m−2 sr−1 at the LORRI pivot wavelength of 0.608 μm. Subtraction of the integrated light of galaxies fainter than the photometric detection limit from the total COB level left a diffuse flux component of unknown origin in the range 8.8 ± 4.9 (1.8 stat., 4.5 sys.) nW m−2 sr−1 to 11.9 ± 4.6 (1.8 stat., 4.2 sys.) nW m−2 sr−1. Explaining it with undetected galaxies requires the assumption that the galaxy count faint-end slope steepens markedly at V > 24 or that existing surveys are missing half the galaxies with V < 30.

Studies of extragalactic background light with TeV BL Lacertae objects

ABSTRACT Very high energy (VHE; E ≥ 100 GeV) gamma-rays from cosmological distances are attenuated by the extragalactic background light (EBL) in the infrared to ultraviolet bands. By contrasting measured versus intrinsic emission,we can derive the EBL photon density. However, we do not know the intrinsic spectra and the EBL separately, only their combined effect. Here we first present a flexible model-dependent optical depth method to study the EBL by fitting the emission spectra of TeV BL Lacertae objects (BL Lacs) via a one-zone leptonic synchrotron self-Compton model (SSC). We have little information about electron energy distributions (EEDs) in the jet, which is critically important to build spectral energy distributions (SEDs) in the SSC scenario. Based on current particle acceleration models, we use two types of EEDs to fit the observed spectra: a power-law log-parabola (PLLP) EED and a broken power-law (BPL) EED. We find that the upper limit of the EBL density is about 30 n W m−2 sr−1, which is similar to the published measurement. Furthermore, we propose an unprecedented method to test the radiation mechanisms involved in TeV objects, by simply comparing the reduced EBL density with the limit obtained by galaxy counts. We demonstrate that for some BL Lacs, at least, the one-zone SSC model should be reconsidered.

Searching for cross-correlation between stochastic gravitational-wave background and galaxy number counts

ABSTRACT Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) and Advanced Virgo have recently published the upper limit measurement of persistent directional stochastic gravitational-wave background (SGWB) based on data from their first and second observing runs. In this paper, we investigate whether a correlation exists between this maximal likelihood SGWB map and the electromagnetic (EM) tracers of matter structure in the Universe, such as galaxy number counts. The method we develop will improve the sensitivity of future searches for anisotropy in the SGWB and expand the use of SGWB anisotropy to probe the formation of structure in the Universe. In order to compute the cross-correlation, we used the spherical harmonic decomposition of SGWB in multiple frequency bands and converted them into pixel-based sky maps in healpix basis. For the EM part, we use the Sloan Digital Sky Survey alaxy catalogue and form healpix sky maps of galaxy number counts at the same angular resolution as the SGWB maps. We compute the pixel-based coherence between these SGWB and galaxy count maps. After evaluating our results in different SGWB frequency bands and in different galaxy redshift bins, we conclude that the coherence between the SGWB and galaxy number count maps is dominated by the null measurement noise in the SGWB maps, and therefore not statistically significant. We expect the results of this analysis to be significantly improved by using the more sensitive upcoming SGWB measurements based on the third observing run of Advanced LIGO and Advanced Virgo.

Measuring angular N -point correlations of binary black hole merger gravitational-wave events with hierarchical Bayesian inference

Advanced LIGO and Virgo have detected ten binary black hole mergers by the end of their second observing run. These mergers have already allowed constraints to be placed on the population distribution of black holes in the Universe, which will only improve with more detections and increasing sensitivity of the detectors. In this paper we develop techniques to measure the angular distribution of black hole mergers by measuring their statistical N-point correlations through hierarchical Bayesian inference. We apply it to the special case of two-point angular correlations using a Legendre polynomial basis on the sky. Building on the mixture model formalism introduced in Ref.[1] we show how one can measure two-point correlations with no threshold on significance, allowing us to target the ensemble of sub-threshold binary black hole mergers not resolvable with the current generation of ground based detectors. We also show how one can use these methods to correlate gravitational waves with other probes of large scale angular structure like galaxy counts, and validate both techniques through simulations.

STRIDES: Spectroscopic and photometric characterization of the environment and effects of mass along the line of sight to the gravitational lenses DES J0408–5354 and WGD 2038–4008

ABSTRACT In time-delay cosmography, three of the key ingredients are (1) determining the velocity dispersion of the lensing galaxy, (2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and (3) estimating the external convergence κext from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quad quasar systems, DES J0408–5354 and WGD 2038–4008 . We use the Gemini, Magellan, and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408–5354 and two groups in WGD 2038–4008 using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the ‘flexion shift’ criterion. We determine the probability distribution function of the external convergence κext for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408–5354 is located in a significantly underdense environment, leading to a tight (width $\sim 3{{\ \rm per\ cent}}$), negative-value κext distribution. On the other hand, WGD 2038–4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter κext of $\sim 1{{\ \rm per\ cent}}$, as long as no external shear constraints are imposed.