Non-Gaussian shape discrimination with spectroscopic galaxy surveys

We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A bacusS ummitN-body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

We examine the clustering properties of low-power radio galaxies at redshift 0.4 < z < 0.8, using data from the 2SLAQ Luminous Red Galaxy (LRG) survey, and find that radio-detected LRGs (with typical optical luminosities of 3–5L* and 1.4 GHz radio powers in the range 1024–1026W Hz−1) are significantly more clustered than a matched population of radio-quiet (≲1024W Hz−1) LRGs with the same distribution in optical luminosity and colour. The measured scalelength of the two-point cross-correlation function between the full LRG sample and the radio-detected LRGs is 9.57 ± 0.50 h−1Mpc, compared to 8.47 ± 0.27 h−1Mpc for the matched sample of radio-quiet LRGs; while the implied scalelength of the auto-correlation function, r0, is 12.3 ± 1.2 h−1Mpc and 9.02 ± 0.52 h−1Mpc for the radio-detected and radio-quiet samples, respectively. We further interpret our clustering measurements in the halo model framework and demonstrate that the radio-detected LRGs have typical halo masses of 10.1 ± 1.4 × 1013h−1M⊙ and bias of 2.96 ± 0.17, compared to 6.44 ± 0.32 × 1013h−1M⊙ and 2.49 ± 0.02 for the radio-quiet sample. A model in which the radio-detected LRGs are almost all central galaxies within haloes provides the best fit to the measured clustering, and we estimate that at least 30 per cent of all 2SLAQ LRGs with the same clustering amplitude as the radio-detected LRGs are currently radio loud. Our results imply that radio-detected galaxies in the 2SLAQ LRG sample typically occupy more massive haloes than other LRGS of the same optical luminosity, so the probability of finding a radio-loud active galactic nucleus (AGN) in a massive galaxy at z∼ 0.55 is influenced by the halo mass and/or cluster environment in addition to the well-known dependence on optical luminosity. If we model the radio-detected fraction of LRGs, Frad, as a function of halo mass M, then the data are well-fitted by a power law of the form Frad∝M0.65±0.23. The observed relationship between radio emission and clustering strength could plausibly arise either through a higher fuelling rate of gas on to the central black holes of galaxies in the most massive haloes (producing more powerful radio jets) or through the presence of a denser IGM (which would provide a more efficient working surface for the jets, thus boosting their observed radio luminosity). Further work is needed to determine which of these effects is dominant.

Context. Empirical models of galaxy-halo connection such as the halo occupation distribution (HOD) model have been widely used over the past decades to intensively test perturbative models on quasi-linear scales. However, these models fail to reproduce the galaxygalaxy lensing signal on non-linear scales, over-predicting the observed signal by up to 40%. Aims. With ongoing Stage-IV galaxy surveys such as DESI and Euclid that will measure cosmological parameters at sub-percent precision, it is now crucial to precisely model the galaxy-halo connection in order to accurately estimate the theoretical uncertainties of perturbative models. Methods. This paper compares a standard HOD (based on halo mass only) to an extended HOD that incorporates as additional features galaxy assembly bias and local environmental dependencies on halo occupation. These models were calibrated against the observed clustering and galaxy-galaxy lensing signal of eBOSS luminous red galaxies and emission line galaxies in the range 0.6 &lt; z &lt; 1.1. We performed a combined clustering-lensing cosmological analysis on the simulated galaxy samples of both HODs to quantify the systematic budget of perturbative models. Results. By considering not only the mass of the dark matter halos but also these secondary properties, the extended HOD offers a more comprehensive understanding of the connection between galaxies and their surroundings. In particular, we found that the luminous red galaxies preferentially occupy denser and more anisotropic environments. Our results highlight the importance of considering environmental factors in empirical models with an extended HOD that reproduces the observed signal within 20% on scales below 10 h−1 Mpc. Our cosmological analysis reveals that our perturbative model yields similar constraints regardless of the galaxy population, with a better goodness of fit for the extended HOD. These results suggest that the extended HOD should be used to quantify modelling systematics. This extended framework should also prove useful for forward modelling techniques.

The Dark Energy Spectroscopic Instrument (DESI) Survey has obtained a set of spectroscopic measurements of galaxies to validate the final survey design and target selections. To assist in these tasks, we visually inspect DESI spectra of approximately 2500 bright galaxies, 3500 luminous red galaxies (LRGs), and 10,000 emission-line galaxies (ELGs) to obtain robust redshift identifications. We then utilize the visually inspected redshift information to characterize the performance of the DESI operation. Based on the visual inspection (VI) catalogs, our results show that the final survey design yields samples of bright galaxies, LRGs, and ELGs with purity greater than 99%. Moreover, we demonstrate that the precision of the redshift measurements is approximately 10 km s−1 for bright galaxies and ELGs and approximately 40 km s−1 for LRGs. The average redshift accuracy is within 10 km s−1 for the three types of galaxies. The VI process also helps improve the quality of the DESI data by identifying spurious spectral features introduced by the pipeline. Finally, we show examples of unexpected real astronomical objects, such as Lyα emitters and strong lensing candidates, identified by VI. These results demonstrate the importance and utility of visually inspecting data from incoming and upcoming surveys, especially during their early operation phases.

ABSTRACTWe present a void clustering analysis in configuration-space using the completed Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) DR16 samples. These samples consist of Luminous Red Galaxies (LRGs) combined with the high-redshift tail of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) DR12 CMASS galaxies (called as LRG+CMASS sample), Emission Line Galaxies (ELGs), and quasars (QSOs). We build void catalogues from the three eBOSS DR16 samples using a ZOBOV-based algorithm, providing 2814 voids, 1801 voids, and 4347 voids in the LRG+CMASS, ELG, and QSO samples, respectively, spanning the redshift range 0.6 &amp;lt; z &amp;lt; 2.2. We measure the redshift space distortions around voids using the anisotropic void-galaxy cross-correlation function and we extract the distortion parameter β. We test the methodology on realistic simulations before applying it to the data, and we investigate all our systematic errors on these mocks. We find βLRG(z = 0.74) = 0.415 ± 0.087, βELG(z = 0.85) = 0.665 ± 0.125 and βQSO(z = 1.48) = 0.313 ± 0.134, for the LRG+CMASS, ELG, and QSO sample, respectively. The quoted errors include systematic and statistical contributions. In order to convert our measurements in terms of the growth rate fσ8, we use consensus values of linear bias from the eBOSS DR16 companion papers, resulting in the following constraints: fσ8(z = 0.74) = 0.50 ± 0.11, fσ8(z = 0.85) = 0.52 ± 0.10, and fσ8(z = 1.48) = 0.30 ± 0.13. Our measurements are consistent with other measurements from eBOSS DR16 using conventional clustering techniques.

In the standard model of non-linear structure formation, a cosmic web of dark-matter dominated filaments connects dark matter halos. In this paper, we stack the weak lensing signal of an ensemble of filaments between groups and clusters of galaxies. Specifically, we detect the weak lensing signal, using CFHTLenS galaxy ellipticities, from stacked filaments between SDSS-III/BOSS luminous red galaxies (LRGs). As a control, we compare the physical LRG pairs with projected LRG pairs that are more widely separated in redshift space. We detect the excess filament mass density in the projected pairs at the $5\sigma$ level, finding a mass of $(1.6 \pm 0.3) \times 10^{13} M_{\odot}$ for a stacked filament region 7.1 $h^{-1}$ Mpc long and 2.5 $h^{-1}$ Mpc wide. This filament signal is compared with a model based on the three-point galaxy-galaxy-convergence correlation function, as developed in Clampitt, Jain & Takada (2014), yielding reasonable agreement.

We estimate the redshift-dependent, anisotropic clustering signal in the Dark Energy Spectroscopic Instrument (DESI) Year 1 Survey created by tidal alignments of Luminous Red Galaxies (LRGs) and a selection-induced galaxy orientation bias. To this end, we measured the correlation between LRG shapes and the tidal field with DESI’s Year 1 redshifts, as traced by LRGs and Emission-Line Galaxies. We also estimate the galaxy orientation bias of LRGs caused by DESI’s aperture-based selection, and find it to increase by a factor of seven between redshifts 0.4−1.1 due to redder, fainter galaxies falling closer to DESI’s imaging selection cuts. These effects combine to dampen measurements of the quadrupole of the correlation function (ξ2) caused by structure growth on scales of 10–80 h−1 Mpc by about 0.15 per cent for low redshifts (0.4 &amp;lt; z &amp;lt; 0.6) and 0.8 per cent for high (0.8 &amp;lt; z &amp;lt; 1.1), a significant fraction of DESI’s error budget. We provide estimates of the ξ2 signal created by intrinsic alignments that can be used to correct this effect, which is necessary to meet DESI’s forecasted precision on measuring the growth rate of structure. While imaging quality varies across DESI’s footprint, we find no significant difference in this effect between imaging regions in the Legacy Imaging Survey.

The next generation of spectroscopic surveys will have a wealth of photometric data available for use in target selection. Selecting the best targets is likely to be one of the most important hurdles in making these spectroscopic campaigns as successful as possible. Our ability to measure dark energy depends strongly on the types of targets that we are able to select with a given photometric data set. We show in this paper that we will be able to successfully select the targets needed for the next generation of spectroscopic surveys. We also investigate the details of this selection, including optimisation of instrument design and survey strategy in order to measure dark energy. We use color-color selection as well as neural networks to select the best possible emission line galaxies and luminous red galaxies for a cosmological survey. Using the Fisher matrix formalism we forecast the efficiency of each target selection scenario. We show how the dark energy figures of merit change in each target selection regime as a function of target type, survey time, survey density and other survey parameters. We outline the optimal target selection scenarios and survey strategy choices which will be available to the next generation of spectroscopic surveys.

Weak gravitational lensing convergence peaks, the local maxima in weak lensing convergence maps, have been shown to contain valuable cosmological information complementary to commonly used two-point statistics. To exploit the full power of weak lensing for cosmology, we must model baryonic feedback processes because these reshape the matter distribution on non-linear and mildly non-linear scales. We study the impact of baryonic physics on the number density of weak lensing peaks using the FLAMINGO cosmological hydrodynamical simulation suite. We generate ray-traced full-sky convergence maps mimicking the characteristics of a Stage IV weak lensing survey. We compare the number densities of peaks in simulations that have been calibrated to reproduce the observed galaxy mass function and cluster gas fraction or to match a shifted version of these, and that use either thermally driven or jet active galactic nucleus feedback. We show that the differences induced by realistic baryonic feedback prescriptions (typically 5–30 per cent for κ = 0.1–0.4) are smaller than those induced by reasonable variations in cosmological parameters (20–60 per cent for κ = 0.1–0.4) but must be modelled carefully to obtain unbiased results. The reasons behind these differences can be understood by considering the impact of feedback on halo masses, or by considering the impact of different cosmological parameters on the halo mass function. Our analysis demonstrates that, for the range of models we investigated, the baryonic suppression is insensitive to changes in cosmology up to κ ≈ 0.4 and that the higher κ regime is dominated by Poisson noise and cosmic variance.

In Montero-Dorta et al., we show that luminous red galaxies (LRGs) from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) at z ∼ 0.55 can be divided into two groups based on their star formation histories. So-called fast-growing LRGs assemble 80 per cent of their stellar mass at z ∼ 5, whereas slow-growing LRGs reach the same evolutionary state at z ∼ 1.5. We further demonstrate that these two subpopulations present significantly different clustering properties on scales of ∼1−30 Mpc. Here, we measure the mean halo mass of each subsample using the galaxy–galaxy lensing technique, in the ${\sim }190\deg ^2$ overlap of the LRG catalogue and the CS82 and CFHTLenS shear catalogues. We show that fast- and slow-growing LRGs have similar lensing profiles, which implies that they live in haloes of similar mass: $\log (M_{\rm halo}^{\rm fast}/h^{-1}\mathrm{M}_{{\odot }}) = 12.85^{+0.16}_{-0.26}$ and $\log (M_{\rm halo}^{\rm slow}/h^{-1}\mathrm{M}_{{\odot }}) =12.92^{+0.16}_{-0.22}$. This result, combined with the clustering difference, suggests the existence of galaxy assembly bias, although the effect is too subtle to be definitively proven, given the errors on our current weak-lensing measurement. We show that this can soon be achieved with upcoming surveys like DES.

We use matter power spectra from cosmological hydrodynamic simulations to quantify the effect of baryon physics on the weak gravitational lensing shear signal. The simulations consider a number of processes, such as radiative cooling, star formation, supernovae and feedback from active galactic nuclei (AGN). Van Daalen et al. (2011) used the same simulations to show that baryon physics, in particular the strong feedback that is required to solve the overcooling problem, modifies the matter power spectrum on scales relevant for cosmological weak lensing studies. As a result, the use of power spectra from dark matter simulations can lead to significant biases in the inferred cosmological parameters. We show that the typical biases are much larger than the precision with which future missions aim to constrain the dark energy equation of state, w_0. For instance, the simulation with AGN feedback, which reproduces X-ray and optical properties of groups of galaxies, gives rise to a ~40% bias in w_0. We demonstrate that the modification of the power spectrum is dominated by groups and clusters of galaxies, the effect of which can be modelled. We consider an approach based on the popular halo model and show that simple modifications can capture the main features of baryonic feedback. Despite its simplicity, we find that our model, when calibrated on the simulations, is able to reduce the bias in w_0 to a level comparable to the size of the statistical uncertainties for a Euclid-like mission. While observations of the gas and stellar fractions as a function of halo mass can be used to calibrate the model, hydrodynamic simulations will likely still be needed to extend the observed scaling relations down to halo masses of 10 ^12 M_sun/h.

Understanding the links between the activity of supermassive black holes (SMBHs) at the centres of galaxies and their host dark matter haloes is a key question in modern astrophysics. The final data release of the SDSS-IV eBOSS provides the largest contemporary spectroscopic sample of galaxies and quasi-stellar objects (QSOs). Using this sample and covering the redshift interval z = 0.7–1.1, we have measured the clustering properties of the eBOSS QSOs, emission-line galaxies (ELGs), and luminous red galaxies (LRGs). We have also measured the fraction of QSOs as a function of the overdensity defined by the galaxy population. Using these measurements, we investigate how QSOs populate and sample the galaxy population, and how the host dark-matter haloes of QSOs sample the underlying halo distribution. We find that the probability of a galaxy hosting a QSO is independent of the host dark matter halo mass of the galaxy. We also find that about 60 per cent of eBOSS QSOs are hosted by LRGs and about 20–40 per cent of QSOs are hosted by satellite galaxies. We find a slight preference for QSOs to populate satellite galaxies over central galaxies. This is connected to the host halo mass distribution of different types of galaxies. Based on our analysis, QSOs should be hosted by a very broad distribution of haloes, and their occurrence should be modulated only by the efficiency of galaxy formation processes.

We apply a new model for the spherically averaged correlation function at large pair separations to the measurement of the clustering of luminous red galaxies (LRGs) made from the SDSS by Cabre and Gaztanaga(2009). Our model takes into account the form of the BAO peak and the large scale shape of the correlation function. We perform a Monte Carlo Markov chain analysis for different combinations of datasets and for different parameter sets. When used in combination with a compilation of the latest CMB measurements, the LRG clustering and the latest supernovae results give constraints on cosmological parameters which are comparable and in remarkably good agreement, resolving the tension reported in some studies. The best fitting model in the context of a flat, Lambda-CDM cosmology is specified by Omega_m=0.261+-0.013, Omega_b=0.044+-0.001, n_s=0.96+-0.01, H_0=71.6+-1.2 km/s/Mpc and sigma_8=0.80+-0.02. If we allow the time-independent dark energy equation of state parameter to vary, we find results consistent with a cosmological constant at the 5% level using all data sets: w_DE=-0.97+-0.05. The large scale structure measurements by themselves can constrain the dark energy equation of state parameter to w_DE=-1.05+-0.15, independently of CMB or supernovae data. We do not find convincing evidence for an evolving equation of state. We provide a set of "extended distance priors" that contain the most relevant information from the CMB power spectrum and the shape of the LRG correlation function which can be used to constrain dark energy models and spatial curvature. Our model should provide an accurate description of the clustering even in much larger, forthcoming surveys, such as those planned with NASA's JDEM or ESA's Euclid mission.

We present initial results from the Dark Energy Spectroscopic Instrument (DESI) complete calibration of the colour–redshift relation (DC3R2) secondary target survey. Our analysis uses 230 k galaxies that overlap with KiDS-VIKING ugriZYJHKs photometry to calibrate the colour–redshift relation and to inform photometric redshift (photo-z) inference methods of future weak lensing surveys. Together with emission line galaxies (ELGs), luminous red galaxies (LRGs), and the Bright Galaxy Survey (BGS) that provide samples of complementary colour, the DC3R2 targets help DESI to span 56 per cent of the colour space visible to Euclid and LSST with high confidence spectroscopic redshifts. The effects of spectroscopic completeness and quality are explored, as well as systematic uncertainties introduced with the use of common Self-Organizing Maps trained on different photometry than the analysis sample. We further examine the dependence of redshift on magnitude at fixed colour, important for the use of bright galaxy spectra to calibrate redshifts in a fainter photometric galaxy sample. We find that noise in the KiDS-VIKING photometry introduces a dominant, apparent magnitude dependence of redshift at fixed colour, which indicates a need for carefully chosen deep drilling fields, and survey simulation to model this effect for future weak lensing surveys.

Environmental effects have an important influence on cluster galaxies, but studies at very faint magnitudes (R>21) are almost exclusively based on imaging. We present here a very deep spectroscopic survey of galaxies on the line of sight to Coma, based on redshifts obtained with VLT/VIMOS for 715 galaxies in the unprecedented magnitude range 21<R<23 (absolute magnitude -14 to -12). We confirm the substructures previously identified in Coma by Adami et al. (2005a), and identify three new ones. We detect many groups behind Coma: a large structure at z~0.5, the SDSS Great Wall, and a large and very young structure at z~0.054. These structures account for the mass maps derived from a recent weak lensing analysis by Gavazzi et al. (2009). The orbits of dwarf galaxies are probably anisotropic and radial, and could originate from field galaxies radially falling into the cluster. Spectral characteristics of Coma dwarf galaxies show that red or absorption line galaxies have larger stellar masses and are older than blue or emission line galaxies. R<22 galaxies show less prominent absorption lines than R>22 galaxies. This trend is less clear for field galaxies, suggesting that part of the faint Coma galaxies could have been recently injected from the field following the NGC 4911 group infall. We present a list of five Ultra Compact Dwarf galaxy candidates. We also globally confirm spectroscopically our previous results on the galaxy luminosity functions and find that dwarf galaxies follow a red sequence similar to that drawn by bright galaxies. Dwarf galaxies are very abundant in Coma, and are partly field galaxies that have fallen onto the cluster along filaments.