Search

In this investigation, we leverage the combination of the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys Data Release 9, Survey Validation 3, and Year 1 data sets to estimate the conditional luminosity functions and conditional stellar mass functions (CLFs and CSMFs) of galaxies across various halo mass bins and redshift ranges. To support our analysis, we utilize a realistic DESI mock galaxy redshift survey (MGRS) generated from a high-resolution Jiutian simulation. An extended halo-based group finder is applied to both MGRS catalogs and DESI observation. By comparing the r- and z-band luminosity functions (LFs) and stellar mass functions (SMFs) derived using both photometric and spectroscopic data, we quantified the impact of photometric redshift (photo-z) errors on the galaxy LFs and SMFs, especially in the low-redshift bin at the low-luminosity/mass end. By conducting prior evaluations of the group finder using MGRS, we successfully obtain a set of CLF and CSMF measurements from observational data. We find that at low redshift, the faint-end slopes of CLFs and CSMFs below ∼109 h −2 L ⊙ (or h −2 M ⊙) evince a compelling concordance with the subhalo mass functions. After correcting the cosmic variance effect of our local Universe following Chen et al., the faint-end slopes of the LFs/SMFs turn out to also be in good agreement with the slope of the halo mass function.

We present the first result in exploring the gaseous halo and galaxy correlation using the Dark Energy Spectroscopic Instrument survey validation data in the Cosmic Evolution Survey (COSMOS) and Hyper Suprime-Cam field. We obtain multiphase gaseous halo properties in the circumgalactic medium by using 115 quasar spectra (signal-to-noise ratio > 3). We detect Mg ii absorption at redshift 0.6 < z < 2.5, C iv absorption at 1.6 < z < 3.6, and H i absorption associated with the Mg ii and C iv. By crossmatching the COSMOS2020 catalog, we identify the Mg ii and C iv host galaxies in 10 quasar fields at 0.9< z < 3.1. We find that within the impact parameter of 250 kpc, a tight correlation is seen between the strong Mg ii equivalent width and the host galaxy star formation rate. The covering fraction f c of the strong Mg ii selected galaxies, which is the ratio of the absorbing galaxy in a certain galaxy population, shows significant evolution in the main-sequence galaxies and marginal evolution in all the galaxy populations within 250 kpc at 0.9 < z < 2.2. The f c increase in the main-sequence galaxies likely suggests the coevolution of strong Mg ii absorbing gas and the main-sequence galaxies at the cosmic noon. Furthermore, Mg ii and C iv absorbing gas is detected out of the galaxy virial radius, tentatively indicating the feedback produced by the star formation and/or the environmental effects.

We perform SubHalo Abundance Matching (SHAM) studies on UNIT simulations with {σ, Vceil, vsmear}-SHAM and {σ, Vceil, fsat}-SHAM. They are designed to reproduce the clustering on 5-30 of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasi-stellar objects (QSOs) at 0.4 < z < 3.5 from DESI (Dark Energy Spectroscopic Instrument) One Percent Survey. Vceil is the incompleteness of the massive host (sub)haloes and is the key to the generalized SHAM. vsmear models the clustering effect of redshift uncertainties, providing measurements consistent with those from repeat observations. A free satellite fraction fsat is necessary to reproduce the clustering of ELGs. We find ELGs present a more complex galaxy-halo mass relation than LRGs reflected in their weak constraints on σ. LRGs, QSOs, and ELGs show increasing Vceil values, corresponding to the massive galaxy incompleteness of LRGs, the quenched star formation of ELGs and the quenched black hole accretion of QSOs. For LRGs, a Gaussian vsmear presents a better profile for subsamples at redshift bins than a Lorentzian profile used for other tracers. The impact of the statistical redshift uncertainty on ELG clustering is negligible. The best-fitting satellite fraction for DESI ELGs is around 4 per cent, lower than previous estimations for ELGs. The mean halo mass log10((Mvir)) in for LRGs, ELGs, and QSOs are 13.16 ± 0.01, 11.90 ± 0.06, and 12.66 ± 0.45, respectively. Our generalized SHAM algorithms facilitate the production of multitracer galaxy mocks for cosmological tests.

The observations from the Dark Energy Spectroscopic Instrument (DESI) will significantly increase the numbers of known extremely metal-poor stars by a factor of ∼10, improving the sample statistics to study the early chemical evolution of the Milky Way and the nature of the first stars. In this paper we report follow-up observations with high signal-to-noise ratio of nine metal-poor stars identified during the DESI commissioning with the Optical System for Imaging and Low-Resolution Integrated Spectroscopy (OSIRIS) instrument on the 10.4 m Gran Telescopio Canarias. The analysis of the data using a well-vetted methodology confirms the quality of the DESI spectra and the performance of the pipelines developed for the data reduction and analysis of DESI data.

We present the first detection of the baryon acoustic oscillations (BAOs) signal obtained using unblinded data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). From a selected sample of 261 291 luminous red galaxies spanning the redshift interval 0.4 < z < 1.1 and covering 1651 square degrees with a 57.9 per cent completeness level, we report a ∼5σ level BAO detection and the measurement of the BAO location at a precision of 1.7 per cent. Using a bright galaxy sample of 109 523 galaxies in the redshift range 0.1 < z < 0.5, over 3677 square degrees with a 50.0 per cent completeness, we also detect the BAO feature at ∼3σ significance with a 2.6 per cent precision. These first BAO measurements represent an important milestone, acting as a quality control on the optimal performance of the complex robotically actuated, fibre-fed DESI spectrograph, as well as an early validation of the DESI spectroscopic pipeline and data management system. Based on these first promising results, we forecast that DESI is on target to achieve a high-significance BAO detection at sub-per cent precision with the completed 5-yr survey data, meeting the top-level science requirements on BAO measurements. This exquisite level of precision will set new standards in cosmology and confirm DESI as the most competitive BAO experiment for the remainder of this decade.

We present findings of the detection of Magnesium II (Mg ii, λ = 2796, 2803 Å) absorbers from the early data release of the Dark Energy Spectroscopic Instrument (DESI). DESI is projected to obtain spectroscopy of approximately 3 million quasars (QSOs), of which over 99% are anticipated to be at redshifts greater than z > 0.3, such that DESI would be able to observe an associated or intervening Mg ii absorber illuminated by the background QSO. We have developed an autonomous supplementary spectral pipeline that detects these systems through an initial line-fitting process and then confirms the line properties using a Markov Chain Monte Carlo sampler. Based upon a visual inspection of the resulting systems, we estimate that this sample has a purity greater than 99%. We have also investigated the completeness of our sample in regard to both the signal-to-noise properties of the input spectra and the rest-frame equivalent width (W 0) of the absorber systems. From a parent catalog containing 83,207 quasars, we detect a total of 23,921 Mg ii absorption systems following a series of quality cuts. Extrapolating from this occurrence rate of 28.8% implies a catalog at the completion of the five-year DESI survey that will contain over eight hundred thousand Mg ii absorbers. The cataloging of these systems will enable significant further research because they carry information regarding circumgalactic medium environments, the distribution of intervening galaxies, and the growth of metallicity across the redshift range 0.3 ≤ z < 2.5.

We describe the target selection and characteristics of the DESI Peculiar Velocity Survey, the largest survey of peculiar velocities (PVs) using both the fundamental plane (FP) and the Tully-Fisher (TF) relationship planned to date. We detail how we identify suitable early-type galaxies (ETGs) for the FP and suitable late-type galaxies (LTGs) for the TF relation using the photometric data provided by the DESI Legacy Imaging Survey DR9. Subsequently, we provide targets for 373 533 ETGs and 118 637 LTGs within the Dark Energy Spectroscopic Instrument (DESI) 5-yr footprint. We validate these photometric selections using existing morphological classifications. Furthermore, we demonstrate using survey validation data that DESI is able to measure the spectroscopic properties to sufficient precision to obtain PVs for our targets. Based on realistic DESI fibre assignment simulations and spectroscopic success rates, we predict the final DESI PV Survey will obtain ∼133 000 FP-based and ∼53 000 TF-based PV measurements over an area of 14 000 deg2. We forecast the ability of using these data to measure the clustering of galaxy positions and PVs from the combined DESI PV and Bright Galaxy Surveys (BGS), which allows for cancellation of cosmic variance at low redshifts. With these forecasts, we anticipate a 4 per cent statistical measurement on the growth rate of structure at z < 0.15. This is over two times better than achievable with redshifts from the BGS alone. The combined DESI PV and BGS will enable the most precise tests to date of the time and scale dependence of large-scale structure growth at z < 0.15.

The shear measurement from the Dark Energy Camera Legacy Survey (DECaLS) provides an excellent opportunity for galaxy-galaxy lensing study with the Dark Energy Spectroscopic Instrument (DESI) galaxies, given the large (∼9000 deg2) sky overlap. We explore this potential by combining the DESI 1 per cent survey and DECaLS Data Release 8 (DR8). With ∼106 deg2 sky overlap, we achieve significant detection of galaxy-galaxy lensing for Bright Galaxy Survey (BGS) and luminous red galaxy (LRG) as lenses. Scaled to the full BGS sample, we expect the statistical errors to improve from to a promising level of at. This brings stronger requirements for future systematics control. To fully realize such potential, we need to control the residual multiplicative shear bias |m| < 0.006 and the bias in the mean redshift |Δz| < 0.008, requiring the introduced bias in the measurement is

We present an extended validation of semi-analytical, semi-empirical covariance matrices for the two-point correlation function (2PCF) on simulated catalogs representative of luminous red galaxies (LRGs) data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). We run the pipeline on multiple effective Zel'dovich (EZ) mock galaxy catalogs with the corresponding cuts applied and compare the results with the mock sample covariance to assess the accuracy and its fluctuations. We propose an extension of the previously developed formalism for catalogs processed with standard reconstruction algorithms. We consider methods for comparing covariance matrices in detail, highlighting their interpretation and statistical properties caused by sample variance, in particular, non-trivial expectation values of certain metrics even when the external covariance estimate is perfect. With improved mocks and validation techniques, we confirm a good agreement between our predictions and sample covariance. This allows one to generate covariance matrices for comparable data sets without the need to create numerous mock galaxy catalogs with matching clustering, only requiring 2PCF measurements from the data itself. The code used in this paper is publicly available at https://github.com/oliverphilcox/RascalC.

We measure the tidal alignment of the major axes of luminous red galaxies (LRGs) from the Legacy Imaging Survey and use it to infer the artificial redshift-space distortion signature that will arise from an orientation-dependent, surface-brightness selection in the Dark Energy Spectroscopic Instrument (DESI) survey. Using photometric redshifts to downweight the shape–density correlations due to weak lensing, we measure the intrinsic tidal alignment of LRGs. Separately, we estimate the net polarization of LRG orientations from DESI’s fibre-magnitude target selection to be of order 10-2 along the line of sight. Using these measurements and a linear tidal model, we forecast a 0.5 per cent fractional decrease on the quadrupole of the two-point correlation function for projected separations of 40–80 h-1 Mpc. We also use a halo catalogue from the ABACUSSUMMIT cosmological simulation suite to reproduce this false quadrupole.

We utilize 17,000 bright luminous red galaxies (LRGs) from the novel Dark Energy Spectroscopic Instrument Survey Validation spectroscopic sample, leveraging its deep (2.5 hr galaxy-1 exposure time) spectra to characterize the contribution of recently quenched galaxies to the massive galaxy population at 0.4 < z < 1.3. We use Prospector to infer nonparametric star formation histories and identify a significant population of recently quenched galaxies that have joined the quiescent population within the past 1 Gyr. The highest-redshift subset (277 at z > 1) of our sample of recently quenched galaxies represents the largest spectroscopic sample of post-starburst galaxies at that epoch. At 0.4 < z < 0.8, we measure the number density of quiescent LRGs, finding that recently quenched galaxies constitute a growing fraction of the massive galaxy population with increasing look-back time. Finally, we quantify the importance of this population among massive ( log(M⋆/M⊙) > 11.2) LRGs by measuring the fraction of stellar mass each galaxy formed in the gigayear before observation, f 1 Gyr. Although galaxies with f 1 Gyr > 0.1 are rare at z ∼0.4 ( 20.5% of the population), by z ∼0.8, they constitute 3% of massive galaxies. Relaxing this threshold, we find that galaxies with f 1 Gyr > 5% constitute 10% of the massive galaxy population at z ∼0.8. We also identify a small but significant sample of galaxies at z = 1.1-1.3 that formed with f 1 Gyr > 50%, implying that they may be analogs to high-redshift quiescent galaxies that formed on similar timescales. Future analysis of this unprecedented sample promises to illuminate the physical mechanisms that drive the quenching of massive galaxies after cosmic noon.

The PRObabilistic Value-added Bright Galaxy Survey (PROVABGS) catalog will provide measurements of galaxy properties, such as stellar mass (M *), star formation rate (SFR), stellar metallicity (Z), and stellar age (t age), for >10 million galaxies of the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey. Full posterior distributions of the galaxy properties will be inferred using state-of-the-art Bayesian spectral energy distribution (SED) modeling of DESI spectroscopy and Legacy Surveys photometry. In this work, we present the SED model, the neural emulator for the model, and the Bayesian inference framework of PROVABGS. Furthermore, we apply the PROVABGS SED modeling on realistic synthetic DESI spectra and photometry, constructed using the L-Galaxies semi-analytic model. We compare the inferred galaxy properties to the true values of the simulation using a hierarchical Bayesian framework to quantify accuracy and precision. Overall, we accurately infer the true M *, SFR, Z, and t age of the simulated galaxies. However, the priors on galaxy properties induced by the SED model have a significant impact on the posteriors, which we characterize in detail. This work also demonstrates that a joint analysis of spectra and photometry significantly improves the constraints on galaxy properties over photometry alone and is necessary to mitigate the impact of the priors. With the methodology presented and validated in this work, PROVABGS will maximize information extracted from DESI observations and extend current galaxy studies to new regimes and unlock cutting-edge probabilistic analyses. https://github.com/changhoonhahn/provabgs/

The Dark Energy Spectroscopic Instrument (DESI) survey will measure large-scale structures using quasars as direct tracers of dark matter in the redshift range 0.9 < z < 2.1 and using Lyα forests in quasar spectra at z > 2.1. We present several methods to select candidate quasars for DESI, using input photometric imaging in three optical bands (g, r, z) from the DESI Legacy Imaging Surveys and two infrared bands (W1, W2) from the Wide-field Infrared Survey Explorer. These methods were extensively tested during the Survey Validation of DESI. In this paper, we report on the results obtained with the different methods and present the selection we optimized for the DESI main survey. The final quasar target selection is based on a random forest algorithm and selects quasars in the magnitude range of 16.5 < r < 23. Visual selection of ultra-deep observations indicates that the main selection consists of 71% quasars, 16% galaxies, 6% stars, and 7% inconclusive spectra. Using the spectra based on this selection, we build an automated quasar catalog that achieves a fraction of true QSOs higher than 99% for a nominal effective exposure time of ∼1000 s. With a 310 deg−2 target density, the main selection allows DESI to select more than 200 deg−2 quasars (including 60 deg−2 quasars with z > 2.1), exceeding the project requirements by 20%. The redshift distribution of the selected quasars is in excellent agreement with quasar luminosity function predictions.

Using synthetic Lyman-α forests from the Dark Energy Spectroscopic Instrument (DESI) survey, we present a study of the impact of errors in the estimation of quasar redshift on the Lyman-α correlation functions. Estimates of quasar redshift have large uncertainties of a few hundred km s−1 due to the broadness of the emission lines and the intrinsic shifts from other emission lines. We inject Gaussian random redshift errors into the mock quasar catalogues, and measure the auto-correlation and the Lyman-α-quasar cross-correlation functions. We find a smearing of the BAO feature in the radial direction, but changes in the peak position are negligible. However, we see a significant unphysical correlation for small separations transverse to the line of sight which increases with the amplitude of the redshift errors. We interpret this contamination as a result of the broadening of emission lines in the measured mean continuum, caused by quasar redshift errors, combined with the unrealistically strong clustering of the simulated quasars on small scales.

We investigate using three-point statistics in constraining the galaxy-halo connection. We show that for some galaxy samples, the constraints on the halo occupation distribution parameters are dominated by the three-point function signal (over its two-point counterpart). We demonstrate this on mock catalogues corresponding to the Luminous red galaxies (LRGs), Emission-line galaxies (ELGs), and quasars (QSOs) targeted by the Dark Energy Spectroscopic Instrument (DESI) Survey. The projected three-point function for triangle sides less up to 20 h-1 Mpc measured from a cubic Gpc of data can constrain the characteristic minimum mass of the LRGs with a preci sion of 0.46 per cent. For comparison, similar constraints from the projected two-point function are 1.55 per cent. The improvements for the ELGs and QSOs targets are more modest. In the case of the QSOs, it is caused by the high shot-noise of the sample, and in the case of the ELGs, it is caused by the range of halo masses of the host haloes. The most time-consuming part of our pipeline is the measurement of the three-point functions. We adopt a tabulation method, proposed in earlier works for the two-point function, to significantly reduce the required compute time for the three-point analysis.

Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogues to validate data analysis pipelines and identify potential systematics. We compare three N-body simulation codes, abacus, gadget-2, and swift, to investigate the regimes in which their results agree. We run N-body simulations at three different mass resolutions, 6.25 × 108, 2.11 × 109, and 5.00 × 109 h-1 M⊙, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the Dark Energy Spectroscopic Instrument (DESI) statistical error for s > 20 h-1 Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of 2.1 × 109 h-1 M⊙ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than 20 h-1 Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1 per cent for k ≤ 10 h Mpc-1. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi-N-body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy-halo relation will be presented in a follow-up study.

Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches ∼ 20 h-3, Gpc3. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. AbacusSummit is a suite of high-resolution dark-matter-only simulations designed for this purpose, with 200,-3}, Gpc3 (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-N-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. k < 0.3h Mpc-1 for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in AbacusSummit suite to extend the effective volume ∼20 times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

We have updated and applied a convolutional neural network (CNN) machine-learning model to discover and characterize damped Lyα systems (DLAs) based on Dark Energy Spectroscopic Instrument (DESI) mock spectra. We have optimized the training process and constructed a CNN model that yields a DLA classification accuracy above 99% for spectra that have signal-to-noise ratios (S/N) above 5 per pixel. The classification accuracy is the rate of correct classifications. This accuracy remains above 97% for lower S/N ≈1 spectra. This CNN model provides estimations for redshift and H i column density with standard deviations of 0.002 and 0.17 dex for spectra with S/N above 3 pixel-1. Also, this DLA finder is able to identify overlapping DLAs and sub-DLAs. Further, the impact of different DLA catalogs on the measurement of baryon acoustic oscillations (BAO) is investigated. The cosmological fitting parameter result for BAO has less than 0.61% difference compared to analysis of the mock results with perfect knowledge of DLAs. This difference is lower than the statistical error for the first year estimated from the mock spectra: above 1.7%. We also compared the performances of the CNN and Gaussian Process (GP) models. Our improved CNN model has moderately 14% higher purity and 7% higher completeness than an older version of the GP code, for S/N > 3. Both codes provide good DLA redshift estimates, but the GP produces a better column density estimate by 24% less standard deviation. A credible DLA catalog for the DESI main survey can be provided by combining these two algorithms.

The quasar target selection for the upcoming survey of the Dark Energy Spectroscopic Instrument (DESI) will be fixed for the next 5 yr. The aim of this work is to validate the quasar selection by studying the impact of imaging systematics as well as stellar and galactic contaminants, and to develop a procedure to mitigate them. Density fluctuations of quasar targets are found to be related to photometric properties such as seeing and depth of the Data Release 9 of the DESI Legacy Imaging Surveys. To model this complex relation, we explore machine learning algorithms (random forest and multilayer perceptron) as an alternative to the standard linear regression. Splitting the footprint of the Legacy Imaging Surveys into three regions according to photometric properties, we perform an independent analysis in each region, validating our method using extended Baryon Oscillation Spectroscopic Survey (eBOSS) EZ-mocks. The mitigation procedure is tested by comparing the angular correlation of the corrected target selection on each photometric region to the angular correlation function obtained using quasars from the Sloan Digital Sky Survey (SDSS) Data Release 16. With our procedure, we recover a similar level of correlation between DESI quasar targets and SDSS quasars in two-thirds of the total footprint and we show that the excess of correlation in the remaining area is due to a stellar contamination that should be removed with DESI spectroscopic data. We derive the Limber parameters in our three imaging regions and compare them to previous measurements from SDSS and the 2dF QSO Redshift Survey.

We measure the 1D Ly α power spectrum P1D from Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ), The Spectral Quasar Absorption Database (SQUAD), and XQ-100 quasars using the optimal quadratic estimator. We combine KODIAQ and SQUAD at the spectrum level, but perfo a separate XQ-100 estimation to control its large resolution corrections in check. Our final analysis measures P1D at scales k < 0.1 s km-1 between redshifts $z$ = 2.0-4.6 using 538 quasars. This sample provides the largest number of high-resolution, high-S/N observations; and combined with the power of optimal estimator it provides exceptional precision at small scales. These small-scale modes (k 0.02 s km-1), unavailable in Sloan Digital Sky Survey and Dark Energy Spectroscopic Instrument analyses, are sensitive to the theal state and reionization history of the intergalactic medium, as well as the nature of dark matter. As an example, a simple Fisher forecast analysis estimates that our results can improve small-scale cut-off sensitivity by more than a factor of 2.

We characterize the selection cuts and clustering properties of a magnitude-limited sample of bright galaxies that is part of the Bright Galaxy Survey (BGS) of the Dark Energy Spectroscopic Instrument (DESI) using the ninth data release of the Legacy Imaging Surveys (DR9). We describe changes in the DR9 selection compared to the DR8 one and we also compare the DR9 selection in three distinct regions: BASS/MzLS in the north Galactic Cap (NGC), DECaLS in the NGC, and DECaLS in the south Galactic Cap (SGC). We investigate the systematics associated with the selection and assess its completeness by matching the BGS targets with the Galaxy and Mass Assembly (GAMA) survey. We measure the angular clustering for the overall bright sample (rmag ≤ 19.5) and as function of apparent magnitude and colour. This enables to deteine the clustering strength r0 and slope γby fitting a power-law model that can be used to generate accurate mock catalogues for this tracer. We use a counts-in-cells technique to explore higher order statistics and cross-correlations with external spectroscopic data sets in order to check the evolution of the clustering with redshift and the redshift distribution of the BGS targets using clustering redshifts. While this work validates the properties of the BGS bright targets, the final target selection pipeline and clustering properties of the entire DESI BGS will be fully characterized and validated with the spectroscopic data of Survey Validation.

We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging data set, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The photometric LRG sample in this study contains 2.7 million objects over the redshift range of 0.4 < z < 0.9 over 5655 deg2. We have developed new photometric redshift (photo-z) estimates using the Legacy Survey DECam and WISE photometry, with σNMAD = 0.02 precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise ratio while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a five-parameter HOD model fit using the MultiDark halo catalogue, we find that there is little evolution in HOD parameters except at the highest redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of Markov chain Monte Carlo fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogues, and present a new technique for finding best-fitting parameters in this situation. Accompanying this paper, we have released the Photometric Redshifts for the Legacy Surveys catalogue of photo-z's obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 data set. This catalogue provides accurate photometric redshifts for objects with z < 21 over more than 16 000 deg2 of sky.

We evaluate the impact of imaging systematics on the clustering of luminous red galaxies (LRG), emission-line galaxies (ELG), and quasars (QSO) targeted for the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. Using Data Release 7 of the DECam Legacy Survey, we study the effects of astrophysical foregrounds, stellar contamination, differences between north galactic cap and south galactic cap measurements, and variations in imaging depth, stellar density, galactic extinction, seeing, airmass, sky brightness, and exposure time before presenting survey masks and weights to mitigate these effects. With our sanitized samples in hand, we conduct a preliminary analysis of the clustering amplitude and evolution of the DESI main targets. From measurements of the angular correlation functions, we determine power law fits r0 = 7.78 ± 0.26 h−1Mpc, γ = 1.98 ± 0.02 for LRGs and r0 = 5.45 ± 0.1 h−1Mpc, γ = 1.54 ± 0.01 for ELGs. Additionally, from the angular power spectra, we measure the linear biases and model the scale-dependent biases in the weakly non-linear regime. Both sets of clustering measurements show good agreement with survey requirements for LRGs and ELGs, attesting that these samples will enable DESI to achieve precise cosmological constraints. We also present clustering as a function of magnitude, use cross-correlations with external spectroscopy to infer dN/dz and measure clustering as a function of luminosity, and probe higher order clustering statistics through counts-in-cells moments.

Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg2 overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-yr long multicolour (griz) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to r∼23.5 mag) powered by less massive black holes (with masses ≳ 108.5M⊙) at high redshift for the first time. We find five candidates with significant (at >99.74 per cent single-frequency significance in at least two bands with a global p-value of ∼7 × 10-4 -3 × 10-3 accounting for the look-elsewhere effect) periodicity with observed periods of ∼3-5 yr (i.e. 1-2 yr in rest frame) having ∼4-6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of ∼0.8+0.5-0.3 per cent or ∼1.1+0.7-0.5 quasar per deg2. Our detection rate is 4-80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. We discuss implications on the future direct detection of low-frequency gravitational waves. Continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins.

We report on SPT-CLJ2011-5228, a giant system of arcs created by a cluster at z = 1.06. The arc system is notable for the presence of a bright central image. The source is a Lyman break galaxy at z s = 2.39 and the mass enclosed within the Einstein ring of radius 14 arcsec is . We perform a full reconstruction of the light profile of the lensed images to precisely infer the parameters of the mass distribution. The brightness of the central image demands that the central total density profile of the lens be shallow. By fitting the dark matter as a generalized Navarro-Frenk-White profile - with a free parameter for the inner density slope - we find that the break radius is kpc, and that the inner density falls with radius to the power -0.38 ±0.04 at 68% confidence. Such a shallow profile is in strong tension with our understanding of relaxed cold dark matter halos; dark matter-only simulations predict that the inner density should fall as . The tension can be alleviated if this cluster is in fact a merger; a two-halo model can also reconstruct the data, with both clumps (density varying as and ) much more consistent with predictions from dark matter-only simulations. At the resolution of our Dark Energy Survey imaging, we are unable to choose between these two models, but we make predictions for forthcoming Hubble Space Telescope imaging that will decisively distinguish between them.

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey and the Large Synoptic Survey Telescope, which will discover SNe by the thousands. Spectroscopic resources are limited, and so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey.