Search

We present two cases of Ly$\alpha$ changing-look (CL) quasars (J1306 and J1512) along with two additional candidates (J1511 and J1602), all discovered serendipitously at $z >2$ through the Dark Energy Spectroscopic Instrument (DESI) and the Sloan Digital Sky Survey (SDSS). It is the first time to capture CL events in Ly$\alpha$ at high redshift, which is crucial for understanding underlying mechanisms driving the CL phenomenon and the evolution of high-redshift quasars and galaxies. The variability of all four sources is confirmed by the significant change of amplitude in the $r$ band ($|r_{\rm DESI}-r_{\rm SDSS}| >0.5 \ \rm mag$). We find that the accretion rate in the dim state for these CL objects corresponds to a relatively low value ($\mathscr{\dot M} \approx 2\times10^{-3}$), which suggests that the inner region of the accretion disk might be in transition between the Advection Dominated Accretion Flow ($\mathscr{\dot M}<10^{-3}\sim 10^{-2}$) and the canonical accretion disk (optically thick, geometrically thin). However, unlike in C {\sc iv} CL quasars in which broad Ly$\alpha$ remained, the broad C {\sc iv} may still persist after a CL event occurs in Ly$\alpha$, making the physical origin of the CL and ionization mechanism event more puzzling and interesting.

We present new Galactic reddening maps of the high Galactic latitude sky using DESI imaging and spectroscopy. We directly measure the reddening of 2.6 million stars by comparing the observed stellar colors in $g-r$ and $r-z$ from DESI imaging with the synthetic colors derived from DESI spectra from the first two years of the survey. The reddening in the two colors is on average consistent with the \cite{fitzpatrick_correcting_1999} extinction curve with $R_\mathrm{V}=3.1$. We find that our reddening maps differ significantly from the commonly used \cite{schlegel_maps_1998} (SFD) reddening map (by up to 80 mmag in $E(B-V)$), and we attribute most of this difference to systematic errors in the SFD map. To validate the reddening map, we select a galaxy sample with extinction correction based on our reddening map, and this yields significantly better uniformity than the SFD extinction correction. Finally, we discuss the potential systematic errors in the DESI reddening measurements, including the photometric calibration errors that are the limiting factor on our accuracy. The $E(g-r)$ and $E(g-r)$ maps presented in this work, and for convenience their corresponding $E(B-V)$ maps with SFD calibration, are publicly available., Comment: Submitted to the Open Journal of Astrophysics. Associated data files: https://data.desi.lbl.gov/public/papers/mws/desi_dust/y2/v1/maps/

The Dark Energy Spectroscopic Instrument (DESI) Peculiar Velocity Survey aims to measure the peculiar velocities of early and late type galaxies within the DESI footprint using both the Fundamental Plane and Tully-Fisher relations. Direct measurements of peculiar velocities can significantly improve constraints on the growth rate of structure, reducing uncertainty by a factor of approximately 2.5 at redshift 0.1 compared to the DESI Bright Galaxy Survey's redshift space distortion measurements alone. We assess the quality of stellar velocity dispersion measurements from DESI spectroscopic data. These measurements, along with photometric data from the Legacy Survey, establish the Fundamental Plane relation and determine distances and peculiar velocities of early-type galaxies. During Survey Validation, we obtain spectra for 6698 unique early-type galaxies, up to a photometric redshift of 0.15. 64\% of observed galaxies (4267) have relative velocity dispersion errors below 10\%. This percentage increases to 75\% if we restrict our sample to galaxies with spectroscopic redshifts below 0.1. We use the measured central velocity dispersion, along with photometry from the DESI Legacy Imaging Surveys, to fit the Fundamental Plane parameters using a 3D Gaussian maximum likelihood algorithm that accounts for measurement uncertainties and selection cuts. In addition, we conduct zero-point calibration using the absolute distance measurements to the Coma cluster, leading to a value of the Hubble constant, $H_0 = 76.05 \pm 0.35$(statistical) $\pm 0.49$(systematic FP) $\pm 4.86$(statistical due to calibration) $\mathrm{km \ s^{-1} Mpc^{-1}}$. This $H_0$ value is within $2\sigma$ of Planck Cosmic Microwave Background results and within $1\sigma$, of other low redshift distance indicator-based measurements., Comment: 18 pages, 9 figures, 2 tables. Submitted for publication in MNRAS

We present the identification of changing-look active galactic nuclei (CL-AGNs) from the Dark Energy Spectroscopic Instrument First Data Release and Sloan Digital Sky Survey Data Release 16 at z \leq 0.9. To confirm the CL-AGNs, we utilize spectral flux calibration assessment via an [O\,{\sc iii}]-based calibration, pseudo-photometry examination, and visual inspection. This rigorous selection process allows us to compile a statistical catalog of 561 CL-AGNs, encompassing 527 $\rm H\beta$, 149$\rm H\alpha$, and 129 Mg II CL behaviors. In this sample, we find 1) a 283:278 ratio of turn-on to turn-off CL-AGNs. 2) the critical value for CL events is confirmed around Eddington ratio \sim 0.01. 3) a strong correlation between the change in the luminosity of the broad emission lines (BEL) and variation in the continuum luminosity, with Mg II and $\rm H\beta$ displaying similar responses during CL phases. 4) the Baldwin-Phillips-Terlevich diagram for CL-AGNs shows no statistically difference from the general AGN catalog. 5) five CL-AGNs are associated with asymmetrical mid-infrared flares, possibly linked to tidal disruption events. Given the large CL-AGNs and the stochastic sampling of spectra, we propose that some CL events are inherently due to typical AGN variability during low accretion rates, particularly for CL events of the singular BEL. Finally, we introduce a Peculiar CL phase, characterized by a gradual decline over decades in the light curve and the complete disappearance of entire BEL in faint spectra, indicative of a real transition in the accretion disk., Comment: Submitted to ApJS, comments welcome

Post-starburst galaxies (PSBs) are young quiescent galaxies that have recently experienced a rapid decrease in star formation, allowing us to probe the fast-quenching period of galaxy evolution. In this work, we obtained HST WFC3/F110W imaging to measure the sizes of 171 massive ($\mathrm{log(M_{*}/M_{\odot})\sim\,11)}$ spectroscopically identified PSBs at $1

We infer the growth of large scale structure over the redshift range $0.4\lesssim z \lesssim 1$ from the cross-correlation of spectroscopically calibrated Luminous Red Galaxies (LRGs) selected from the Dark Energy Spectroscopic Instrument (DESI) legacy imaging survey with CMB lensing maps reconstructed from the latest Planck and ACT data. We adopt a hybrid effective field theory (HEFT) model that robustly regulates the cosmological information obtainable from smaller scales, such that our cosmological constraints are reliably derived from the (predominantly) linear regime. We perform an extensive set of bandpower- and parameter-level systematics checks to ensure the robustness of our results and to characterize the uniformity of the LRG sample. We demonstrate that our results are stable to a wide range of modeling assumptions, finding excellent agreement with a linear theory analysis performed on a restricted range of scales. From a tomographic analysis of the four LRG photometric redshift bins we find that the rate of structure growth is consistent with $\Lambda$CDM with an overall amplitude that is $\simeq5-7\%$ lower than predicted by primary CMB measurements with modest $(\sim2\sigma)$ statistical significance. From the combined analysis of all four bins and their cross-correlations with Planck we obtain $S_8 = 0.765\pm0.023$, which is less discrepant with primary CMB measurements than previous DESI LRG cross Planck CMB lensing results. From the cross-correlation with ACT we obtain $S_8 = 0.790^{+0.024}_{-0.027}$, while when jointly analyzing Planck and ACT we find $S_8 = 0.775^{+0.019}_{-0.022}$ from our data alone and $\sigma_8 = 0.772^{+0.020}_{-0.023}$ with the addition of BAO data. These constraints are consistent with the latest Planck primary CMB analyses at the $\simeq 1.6-2.2\sigma$ level, and are in excellent agreement with galaxy lensing surveys., Comment: 60 pages, 26 figures, comments welcome

We present a high-significance cross-correlation of CMB lensing maps from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) with spectroscopically calibrated luminous red galaxies (LRGs) from the Dark Energy Spectroscopic Instrument (DESI). We detect this cross-correlation at a significance of 38$\sigma$; combining our measurement with the Planck Public Release 4 (PR4) lensing map, we detect the cross-correlation at 50$\sigma$. Fitting this jointly with the galaxy auto-correlation power spectrum to break the galaxy bias degeneracy with $\sigma_8$, we perform a tomographic analysis in four LRG redshift bins spanning $0.4 \le z \le 1.0$ to constrain the amplitude of matter density fluctuations through the parameter combination $S_8^\times = \sigma_8 \left(\Omega_m / 0.3\right)^{0.4}$. Prior to unblinding, we confirm with extragalactic simulations that foreground biases are negligible and carry out a comprehensive suite of null and consistency tests. Using a hybrid effective field theory (HEFT) model that allows scales as small as $k_{\rm max}=0.6$ $h/{\rm Mpc}$, we obtain a 3.3% constraint on $S_8^\times = \sigma_8 \left(\Omega_m / 0.3\right)^{0.4} = 0.792^{+0.024}_{-0.028}$ from ACT data, as well as constraints on $S_8^\times(z)$ that probe structure formation over cosmic time. Our result is consistent with the early-universe extrapolation from primary CMB anisotropies measured by Planck PR4 within 1.2$\sigma$. Jointly fitting ACT and Planck lensing cross-correlations we obtain a 2.7% constraint of $S_8^\times = 0.776^{+0.019}_{-0.021}$, which is consistent with the Planck early-universe extrapolation within 2.1$\sigma$, with the lowest redshift bin showing the largest difference in mean. The latter may motivate further CMB lensing tomography analyses at $z<0.6$ to assess the impact of potential systematics or the consistency of the $\Lambda$CDM model over cosmic time., Comment: Prepared for submission to JCAP (47 pages, 13 figures)

The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterising remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way's accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution $\Lambda$CDM magneto-hydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum., Comment: Accepted for publication in MNRAS, 31 pages, 27 figues, 7 tables. The mock catalogues are available at https://data.desi.lbl.gov/public/papers/mws/auridesi/v1

Dark Energy Spectroscopic Instrument (DESI) uses more than 2.4 million Emission Line Galaxies (ELGs) for 3D large-scale structure (LSS) analyses in its Data Release 1 (DR1). Such large statistics enable thorough research on systematic uncertainties. In this study, we focus on spectroscopic systematics of ELGs. The redshift success rate ($f_{\rm goodz}$) is the relative fraction of secure redshifts among all measurements. It depends on observing conditions, thus introduces non-cosmological variations to the LSS. We, therefore, develop the redshift failure weight ($w_{\rm zfail}$) and a per-fibre correction ($\eta_{\rm zfail}$) to mitigate these dependences. They have minor influences on the galaxy clustering. For ELGs with a secure redshift, there are two subtypes of systematics: 1) catastrophics (large) that only occur in a few samples; 2) redshift uncertainty (small) that exists for all samples. The catastrophics represent 0.26\% of the total DR1 ELGs, composed of the confusion between O\,\textsc{ii} and sky residuals, double objects, total catastrophics and others. We simulate the realistic 0.26\% catastrophics of DR1 ELGs, the hypothetical 1\% catastrophics, and the truncation of the contaminated $1.31

We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). We update the Obiwan pipeline, which had previously been developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition makes it possible to simulate the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing power of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying Large-Scale Structure systematic contamination (e.g. Cosmic-Infrared Background) in the Galactic extinction maps is a likely root cause. We additionally observe that the DESI LRGs sample exhibits a complex dependency on a combination of seeing, depth, and intrinsic galaxy brightness, which is not replicated by Obiwan, suggesting discrepancies between the current simulation settings and the actual observations. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRG sample., Comment: 46 pages, 26 figures

We present the MOST Hosts survey (Multi-Object Spectroscopy of Transient Hosts). The survey is planned to run throughout the five years of operation of the Dark Energy Spectroscopic Instrument (DESI) and will generate a spectroscopic catalog of the hosts of most transients observed to date, in particular all the supernovae observed by most public, untargeted, wide-field, optical surveys (PTF/iPTF, SDSS II, ZTF, DECAT, DESIRT). Scientific questions for which the MOST Hosts survey will be useful include Type Ia supernova cosmology, fundamental plane and peculiar velocity measurements, and the understanding of the correlations between transients and their host galaxy properties. Here, we present the first release of the MOST Hosts survey: 21,931 hosts of 20,235 transients. These numbers represent 36% of the final MOST Hosts sample, consisting of 60,212 potential host galaxies of 38,603 transients (a transient can be assigned multiple potential hosts). Of these galaxies, 40% do not appear in the DESI primary target list and therefore require a specific program like MOST Hosts. Of all the transients in the MOST Hosts list, only 26.7% have existing classifications, and so the survey will provide redshifts (and luminosities) for nearly 30,000 transients. A preliminary Hubble diagram and a transient luminosity-duration diagram are shown as examples of future potential uses of the MOST Hosts survey. The survey will also provide a training sample of spectroscopically observed transients for photometry-only classifiers, as we enter an era when most newly observed transients will lack spectroscopic classification. The MOST Hosts DESI survey data will be released through the Wiserep platform on a rolling cadence and updated to match the DESI releases. Dates of future releases and updates are available through the https://mosthosts.desi.lbl.gov website., Comment: Submitted to ApJS

We present strong constraints on the spacetime variation of the fine-structure constant $\alpha$ using the Dark Energy Spectroscopic Instrument (DESI). In this pilot work, we utilize $\sim110,000$ galaxies with strong and narrow O III $\lambda\lambda$4959,5007 emission lines to measure the relative variation $\Delta\alpha/\alpha$ in space and time. The O III doublet is arguably the best choice for this purpose owing to its wide wavelength separation between the two lines and its strong emission in many galaxies. Our galaxy sample spans a redshift range of $0

We present a simple, differentiable method for predicting emission line strengths from rest-frame optical continua using an empirically-determined mapping. Extensive work has been done to develop mock galaxy catalogues that include robust predictions for galaxy photometry, but reliably predicting the strengths of emission lines has remained challenging. Our new mapping is a simple neural network implemented using the JAX Python automatic differentiation library. It is trained on Dark Energy Spectroscopic Instrument Early Release data to predict the equivalent widths (EWs) of the eight brightest optical emission lines (including H$\alpha$, H$\beta$, [O II], and [O III]) from a galaxy's rest-frame optical continuum. The predicted EW distributions are consistent with the observed ones when noise is accounted for, and we find Spearman's rank correlation coefficient $\rho_s > 0.87$ between predictions and observations for most lines. Using a non-linear dimensionality reduction technique (UMAP), we show that this is true for galaxies across the full range of observed spectral energy distributions. In addition, we find that adding measurement uncertainties to the predicted line strengths is essential for reproducing the distribution of observed line-ratios in the BPT diagram. Our trained network can easily be incorporated into a differentiable stellar population synthesis pipeline without hindering differentiability or scalability with GPUs. A synthetic catalogue generated with such a pipeline can be used to characterise and account for biases in the spectroscopic training sets used for training and calibration of photo-$z$'s, improving the modelling of systematic incompleteness for the Rubin Observatory LSST and other surveys., Comment: 17 pages, 8 figures, 1 table. Published in MNRAS

Galaxy-galaxy lensing (GGL) and clustering measurements from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) dataset promise to yield unprecedented combined-probe tests of cosmology and the galaxy-halo connection. In such analyses, it is essential to identify and characterise all relevant statistical and systematic errors. In this paper, we forecast the covariances of DESI Y1 GGL+clustering measurements and characterise the systematic bias due to redshift evolution in the lens samples. Focusing on the projected clustering and galaxy-galaxy lensing correlations, we compute a Gaussian analytical covariance, using a suite of N-body and log-normal simulations to characterise the effect of the survey footprint. Using the DESI One Percent Survey data, we measure the evolution of galaxy bias parameters for the DESI Luminous Red Galaxy (LRG) and Bright Galaxy Survey (BGS) samples. We find mild evolution in the LRGs in 0.4 < z < 0.8, subdominant compared to the expected statistical errors. For BGS, we find less evolution effects for brighter absolute magnitude cuts, at the cost of reduced sample size. We find that with a fiducial redshift bin width delta z = 0.1, evolution effects on GGL is negligible across all scales, all fiducial selection cuts, all fiducial redshift bins, given DESI Y1 sample size. Galaxy clustering is more sensitive to evolution due to the bias squared scaling. Nevertheless the redshift evolution effect is insignificant for clustering above the 1-halo scale of 0.1Mpc/h. For studies that wish to reliably access smaller scales, additional treatment of redshift evolution is likely needed. This study serves as a reference for GGL and clustering studies using the DESI Y1 sample, Comment: 19 pages, 15 figures, submitted to MNRAS

High-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created by high-energy emission from active galactic nuclei and certain types of transients. In cases known as extreme coronal line emitting galaxies (ECLEs), these CLs are strong and fade away on a timescale of years. The most likely progenitors of these variable CLs are tidal disruption events (TDEs), which produce sufficient high-energy emission to create and sustain the CLs over these timescales. To test the possible connection between ECLEs and TDEs, we present the most complete variable ECLE rate calculation to date and compare the results to TDE rates from the literature. To achieve this, we search for ECLEs in the Sloan Digital Sky Survey (SDSS). We detect sufficiently strong CLs in 16 galaxies, more than doubling the number previously found in SDSS. We find that none of the nine new ECLEs evolve in a manner consistent with that of the five previously discovered variable ECLEs. Using this sample of five variable ECLEs, we calculate the galaxy-normalized rate of variable ECLEs in SDSS to be $R_\mathrm{G}=3.6~^{+2.6}_{-1.8}~(\mathrm{statistical})~^{+5.1}_{-0.0} (\mathrm{systematic})\times10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$. The mass-normalised rate is $R_\mathrm{M}=3.1~^{+2.3}_{-1.5}~(\mathrm{statistical})~^{+4.4}_{-0.0}~(\mathrm{systematic})\times10^{-17}~\mathrm{M_\odot^{-1}}~\mathrm{yr}^{-1}$ and the volumetric rate is $R_\mathrm{V}=7~^{+20}_{-5}~(\mathrm{statistical})~^{+10}_{-0.0}~(\mathrm{systematic})\times10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$. Our rates are one to two orders of magnitude lower than TDE rates from the literature, which suggests that only 10 to 40 per cent of all TDEs produce variable ECLEs. Additional uncertainties in the rates arising from the structure of the interstellar medium have yet to be included., Comment: 29 pages, 14 figures. Accepted by MNRAS

In this investigation, we leverage the combination of Dark Energy Spectroscopic Instrument Legacy imaging Surveys Data Release 9 (DESI LS DR9), Survey Validation 3 (SV3), and Year 1 (Y1) data sets to estimate the conditional luminosity and stellar mass functions (CLFs & 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 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 $10^{9}h^{-2}L_{\odot}$ (or $h^{-2}M_{\odot}$) evince a compelling concordance with the subhalo mass functions. After correcting the cosmic variance effect of our local Universe following arXiv:1809.00523, the faint end slopes of the LFs/SMFs turn out to be also in good agreement with the slope of the halo mass function., Comment: 28 pages, 13 figures, Accepted for publication in ApJ

We estimate the redshift-dependent, anisotropic clustering signal in DESI's 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 (ELGs). 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 caused by structure growth on scales of 10 - 80 Mpc/h by about 0.15% for low redshifts (0.4

We search for narrow-line optical emission from dark matter decay by stacking dark-sky spectra from the Dark Energy Spectroscopic Instrument (DESI) at the redshift of nearby galaxies from DESI's Bright Galaxy and Luminous Red Galaxy samples. Our search uses regions separated by 5 to 20 arcsecond from the centers of the galaxies, corresponding to an impact parameter of approximately $50\,\rm kpc$. No unidentified spectral line shows up in the search, and we place a line flux limit of $10^{-19}\,\rm{ergs}/\rm{s}/\rm{cm}^{2}/\rm{arcsec}^{2}$ on emissions in the optical band ($3000\lesssim\lambda\lesssim9000 \,\mathring{\rm A}$), which corresponds to $34$ in AB magnitude in a normal broadband detection. This detection limit suggests that the line surface brightness contributed from all dark matter along the line of sight is two orders of magnitude lower than the measured extragalactic background light (EBL), which rules out the possibility that narrow optical-line emission from dark matter decay is a major source of the EBL.

Emission line galaxies (ELGs) are now the preeminent tracers of large-scale structure at z>0.8 due to their high density and strong emission lines, which enable accurate redshift measurements. However, relatively little is known about ELG evolution and the ELG-halo connection, exposing us to potential modeling systematics in cosmology inference using these sources. In this paper, we propose a physical picture of ELGs and improve ELG-halo connection modeling using a variety of observations and simulated galaxy models. We investigate DESI-selected ELGs in COSMOS data, and infer that ELGs are rapidly star-forming galaxies with a large fraction exhibiting disturbed morphology, implying that many of them are likely to be merger-driven starbursts. We further postulate that the tidal interactions from mergers lead to correlated star formation in central-satellite ELG pairs, a phenomenon dubbed "conformity." We argue for the need to include conformity in the ELG-halo connection using galaxy models such as IllustrisTNG, and by combining observations such as the DESI ELG auto-correlation, ELG cross-correlation with Luminous Red Galaxies (LRGs), and ELG-cluster cross-correlation. We also explore the origin of conformity using the UniverseMachine model and elucidate the difference between conformity and the well-known galaxy assembly bias effect., Comment: 25 pages, 20 figures, submitted to MNRAS, comments welcome

We explore the galaxy-halo connection information that is available in low-redshift samples from the early data release of the Dark Energy Spectroscopic Instrument (DESI). We model the halo occupation distribution (HOD) from z=0.1-0.3 using Survey Validation 3 (SV3; a.k.a., the One-Percent Survey) data of the DESI Bright Galaxy Survey (BGS). In addition to more commonly used metrics, we incorporate counts-in-cylinders (CiC) measurements, which drastically tighten HOD constraints. Our analysis is aided by the Python package, galtab, which enables the rapid, precise prediction of CiC for any HOD model available in halotools. This methodology allows our Markov chains to converge with much fewer trial points, and enables even more drastic speedups due to its GPU portability. Our HOD fits constrain characteristic halo masses tightly and provide statistical evidence for assembly bias, especially at lower luminosity thresholds: the HOD of central galaxies in $z\sim0.15$ samples with limiting absolute magnitude $M_r < -20.0$ and $M_r < -20.5$ samples is positively correlated with halo concentration with a significance of 99.9% and 99.5%, respectively. Our models also favor positive central assembly bias for the brighter $M_r < -21.0$ sample at $z\sim0.25$ (94.8% significance), but there is no significant evidence for assembly bias with the same luminosity threshold at $z\sim0.15$. We provide our constraints for each threshold sample's characteristic halo masses, assembly bias, and other HOD parameters. These constraints are expected to be significantly tightened with future DESI data, which will span an area 100 times larger than that of SV3.

We present two galaxy samples, selected from DESI Legacy Imaging Surveys (LS) DR9, with approximately 20,000 square degrees of coverage and spectroscopic redshift distributions designed for cross-correlations such as with CMB lensing, galaxy lensing, and the Sunyaev-Zel'dovich effect. The first sample is identical to the DESI Luminous Red Galaxy (LRG) sample, and the second sample is an extended LRG sample with 2-3 times the DESI LRG density. We present the improved photometric redshifts, tomographic binning and their spectroscopic redshift distributions and imaging systematics weights, and magnification bias coefficients. The catalogs and related data products will be made publicly available. The cosmological constraints using this sample and Planck lensing maps are presented in a companion paper. We also make public the new set of general-purpose photometric redshifts trained using DESI spectroscopic redshifts, which are used in this work, for all galaxies in LS DR9., Comment: Matches the journal version. Associated data files: https://data.desi.lbl.gov/public/papers/c3/lrg_xcorr_2023/. General-purpose photo-z catalogs: https://www.legacysurvey.org/dr9/files/#photo-z-sweeps-9-1-photo-z-sweep-brickmin-brickmax-pz-fits

Galactic conformity is the phenomenon in which a galaxy of a certain physical property is correlated with its neighbors of the same property, implying a possible causal relationship. The observed auto correlations of emission line galaxies (ELGs) from the highly complete DESI One-Percent survey exhibit a strong clustering signal on small scales, providing clear evidence for the conformity effect of ELGs. Building upon the original subhalo abundance matching (SHAM) method developed by Gao et al. (2022, 2023), we propose a concise conformity model to improve the ELG-halo connection. In this model, the number of satellite ELGs is boosted by a factor of $\sim 5$ in the halos whose central galaxies are ELGs. We show that the mean ELG satellite number in such central halos is still smaller than 1, and the model does not significantly increase the overall satellite fraction. With this model, we can well recover the ELG auto correlations to the smallest scales explored with the current data (i.e. $r_{\mathrm{p}} > 0.03$ $\mathrm{Mpc}\,h^{-1}$ in real space and at $s > 0.3$ $\mathrm{Mpc}\,h^{-1}$ in redshift space), while the cross correlations between luminous red galaxies (LRGs) and ELGs are nearly unchanged. Although our SHAM model has only 8 parameters, we further verify that it can accurately describe the ELG clustering in the entire redshift range from $z = 0.8$ to $1.6$. We therefore expect that this method can be used to generate high-quality ELG lightcone mocks for DESI., Comment: 18 pages, 10 figures, accepted by ApJ

In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological $N$-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis. We also provide an optimized package for computing the power spectra and other two-point statistics of an arbitrary galaxy catalog as well as a pipeline for obtaining CV-reduced measurements on any of the AbacusSummit cubic box outputs. We make our scripts publicly available and report a speed improvement of $\sim$10 for a grid size of $N_{\rm mesh} = 256^3$ compared with \texttt{nbodykit}., Comment: 17 pages, 9 figures, public package (for power spectrum and control variates estimation)

Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a Halo Occupation Distribution (HOD) model upon the \textsc{FastPM} simulations, such that the resulting galaxy clustering reproduces high resolution $N$-body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements -- monopole and quadrupole -- to a precision required by the DESI Year 1 Emission Line Galaxy sample down to non-linear scales, i.e. $k<0.5\,h\mathrm{Mpc}$ or $s>10\,\mathrm{Mpc}/h$. Furthermore, we compute covariance matrices based on the resulting \textsc{FastPM} galaxy clustering -- monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate [e.g. power spectrum] on the covariance matrix of the Fourier counterpart [e.g. correlation function]. We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between two to ten per cent Consequently, using the current pipeline we can precisely reproduce the clustering of $N$-body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies., Comment: 21 pages, 18 figures, submitted to MNRAS, Zenodo link https://doi.org/10.5281/zenodo.8185822

We present new spectroscopic and photometric follow-up observations of the known sample of extreme coronal line emitting galaxies (ECLEs) identified in the Sloan Digital Sky Survey (SDSS). With these new data, observations of the ECLE sample now span a period of two decades following their initial SDSS detections. We confirm the nonrecurrence of the iron coronal line signatures in five of the seven objects, further supporting their identification as the transient light echoes of tidal disruption events (TDEs). Photometric observations of these objects in optical bands show little overall evolution. In contrast, mid-infrared (MIR) observations show ongoing long-term declines. The remaining two objects had been classified as active galactic nuclei (AGN) with unusually strong coronal lines rather than being TDE related, given the persistence of the coronal lines in earlier follow-up spectra. We confirm this classification, with our spectra continuing to show the presence of strong, unchanged coronal-line features and AGN-like MIR colours and behaviour. We have constructed spectral templates of both subtypes of ECLE to aid in distinguishing the likely origin of newly discovered ECLEs. We highlight the need for higher cadence, and more rapid, follow-up observations of such objects to better constrain their properties and evolution. We also discuss the relationships between ECLEs, TDEs, and other identified transients having significant MIR variability., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. Note the corrected caption of Figure 1 continued, which in this version correctly refers to 'SDSS J124' rather than the erroneous 'SDSS J1341' in the published version. 29 Pages, 14 Figures

We use angular clustering of luminous red galaxies from the Dark Energy Spectroscopic Instrument (DESI) imaging surveys to constrain the local primordial non-Gaussianity parameter $\fnl$. Our sample comprises over 12 million targets, covering 14,000 square degrees of the sky, with redshifts in the range $0.2< z < 1.35$. We identify Galactic extinction, survey depth, and astronomical seeing as the primary sources of systematic error, and employ linear regression and artificial neural networks to alleviate non-cosmological excess clustering on large scales. Our methods are tested against simulations with and without $\fnl$ and systematics, showing superior performance of the neural network treatment. The neural network with a set of nine imaging property maps passes our systematic null test criteria, and is chosen as the fiducial treatment. Assuming the universality relation, we find $\fnl = 34^{+24(+50)}_{-44(-73)}$ at 68\%(95\%) confidence. We apply a series of robustness tests (e.g., cuts on imaging, declination, or scales used) that show consistency in the obtained constraints. We study how the regression method biases the measured angular power-spectrum and degrades the $\fnl$ constraining power. The use of the nine maps more than doubles the uncertainty compared to using only the three primary maps in the regression. Our results thus motivate the development of more efficient methods that avoid over-correction, protect large-scale clustering information, and preserve constraining power. Additionally, our results encourage further studies of $\fnl$ with DESI spectroscopic samples, where the inclusion of 3D clustering modes should help separate imaging systematics and lessen the degradation in the $\fnl$ uncertainty., Comment: 21 pages, 17 figures, 7 tables (Appendix excluded). Published in MNRAS

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 (LRG) data collected during the initial two 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, nontrivial 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 datasets 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., Comment: 19 pages, 1 figure. Published in MNRAS. Code available at https://github.com/oliverphilcox/RascalC, table and figure data available at https://dx.doi.org/10.5281/zenodo.7750637

The One-Percent survey of the Dark Energy Spectroscopic Instrument collected ~ 270k emission line galaxies (ELGs) at 0.8 < z < 1.6. The high completeness of the sample allowed the clustering to be measured down to scales never probed before, 0.04 Mpc/h in rp for the projected 2-point correlation function (2PCF) and 0.17 Mpc/h in galaxy pair separation s for the 2PCF monopole and quadrupole. The most striking feature of the measurements is a strong signal at the smallest scales, below 0.2 Mpc/h in rp and 1 Mpc/h in s. We analyze these data in the halo occupation distribution framework. We consider different distributions for central galaxies, a standard power law for satellites with no condition on the presence of a central galaxy and explore several extensions of these models. For all considered models, the mean halo mass of the sample is found to be log10 ~ 11.9. We obtain a satellite mean occupation function which agrees with physically motivated ELG models only if we introduce central-satellite conformity, meaning that the satellite occupation is conditioned by the presence of central galaxies of the same type. To achieve in addition a good modeling of the clustering between 0.1 and 1 Mpc/h in rp, we allow for ELG positioning outside of the halo virial radius and find 0.5% of ELGs residing in the outskirts of halos. Furthermore, the satellite velocity dispersion inside halos is found to be ~ 30% larger than that of the halo dark matter particles. These are the main findings of our work. We investigate assembly bias as a function of halo concentration, local density or local density anisotropies and observe no significant change in our results. We split the data sample in two redshift bins and report no significant evolution with redshift. Lastly, changing the cosmology in the modeling impacts only slightly our results.

We present the probabilistic stellar mass function (pSMF) of galaxies in the DESI Bright Galaxy Survey (BGS), observed during the One-Percent Survey. The One-Percent Survey was one of DESI's survey validation programs conducted from April to May 2021, before the start of the main survey. It used the same target selection and similar observing strategy as the main survey and successfully observed the spectra and redshifts of 143,017 galaxies in the $r < 19.5$ magnitude-limited BGS Bright sample and 95,499 galaxies in the fainter surface brightness and color selected BGS Faint sample over $z < 0.6$. We derive pSMFs from posteriors of stellar mass, $M_*$, inferred from DESI photometry and spectroscopy using the Hahn et al. (2022a; arXiv:2202.01809) PRObabilistic Value-Added BGS (PROVABGS) Bayesian SED modeling framework. We use a hierarchical population inference framework that statistically and rigorously propagates the $M_*$ uncertainties. Furthermore, we include correction weights that account for the selection effects and incompleteness of the BGS observations. We present the redshift evolution of the pSMF in BGS as well as the pSMFs of star-forming and quiescent galaxies classified using average specific star formation rates from PROVABGS. Overall, the pSMFs show good agreement with previous stellar mass function measurements in the literature. Our pSMFs showcase the potential and statistical power of BGS, which in its main survey will observe >100$\times$ more galaxies. Moreover, we present the statistical framework for subsequent population statistics measurements using BGS, which will characterize the global galaxy population and scaling relations at low redshifts with unprecedented precision., Comment: 25 pages, 12 figures; data used to generate figures is available at https://doi.org/10.5281/zenodo.8018936; submitted to ApJ

In the current Dark Energy Spectroscopic Instrument (DESI) survey, emission line galaxies (ELGs) and luminous red galaxies (LRGs) are essential for mapping the dark matter distribution at $z \sim 1$. We measure the auto and cross correlation functions of ELGs and LRGs at $0.80.1$ $\mathrm{Mpc}\,h^{-1}$. We can also reproduce the auto correlations of ELGs at $r_{\mathrm{p}}>0.3$ $\mathrm{Mpc}\,h^{-1}$ ($s>1$ $\mathrm{Mpc}\,h^{-1}$) in real (redshift) space. Although our model has only seven parameters, we show that it can be extended to higher redshifts and reproduces the observed auto correlations of ELGs in the whole range of $0.8

We present the first comprehensive Halo Occupation Distribution (HOD) analysis of the DESI One-Percent survey Luminous Red Galaxy (LRG) and Quasi-Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32 Mpc/h in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of fsat = 11+-1%, a mean halo mass of log10 Mh = 13.40+0.02-0.02, and a linear bias of blin = 1.93+0.06-0.04. For LRGs in 0.6 < z < 0.8, we find fsat = 14+-1%, log10 Mh = 13.24+0.02-0.02, and blin = 2.08+0.03-0.03. For QSOs, we infer fsat = 3+8-2%, log10 Mh = 12.65+0.09-0.04, and blin = 2.63+0.37-0.26 in redshift range 0.8 < z < 2.1. Using these fits, we generate a large suite of high-fidelity galaxy mocks. We also study the redshift-evolution of the DESI LRG sample from z = 0.4 up to z = 1.1, revealing significant and interesting trends in mean halo mass, linear bias, and satellite fraction., Comment: Submitted to MNRAS, comments welcome

The joint analysis of different cosmological probes, such as galaxy clustering and weak lensing, can potentially yield invaluable insights into the nature of the primordial Universe, dark energy and dark matter. However, the development of high-fidelity theoretical models that cover a wide range of scales and redshifts is a necessary stepping-stone. Here, we present public high-resolution weak lensing maps on the light cone, generated using the $N$-body simulation suite AbacusSummit in the Born approximation, and accompanying weak lensing mock catalogues, tuned via fits to the Early Data Release small-scale clustering measurements of the Dark Energy Spectroscopic Instrument (DESI). Available in this release are maps of the cosmic shear, deflection angle and convergence fields at source redshifts ranging from $z = 0.15$ to 2.45 with $\Delta z = 0.05$ as well as CMB convergence maps ($z \approx 1090$) for each of the 25 ${\tt base}$-resolution simulations ($L_{\rm box} = 2000\,h^{-1}{\rm Mpc}$, $N_{\rm part} = 6912^3$) as well as for the two ${\tt huge}$ simulations ($L_{\rm box} = 7500\,h^{-1}{\rm Mpc}$, $N_{\rm part} = 8640^3$) at the fiducial AbacusSummit cosmology ($Planck$ 2018). The pixel resolution of each map is 0.21 arcmin, corresponding to a HEALPiX $N_{\rm side}$ of 16384. The sky coverage of the ${\tt base}$ simulations is an octant until $z \approx 0.8$ (decreasing to about 1800 deg$^2$ at $z \approx 2.4$), whereas the ${\tt huge}$ simulations offer full-sky coverage until $z \approx 2.2$. Mock lensing source catalogues are sampled matching the ensemble properties of the Kilo-Degree Survey, Dark Energy Survey, and Hyper-Suprime Cam weak lensing datasets. The produced mock catalogues are validated against theoretical predictions for various clustering and lensing statistics such as galaxy clustering multipoles, galaxy-shear and shear-shear, showing excellent agreement., Comment: 21 pages, 13 figures, data available at https://app.globus.org/file-manager?origin_id=3dd6566c-eed2-11ed-ba43-09d6a6f08166&path=%2F

The full-shape correlations of the Lyman alpha (Ly$\alpha$) forest contain a wealth of cosmological information through the Alcock-Paczy\'{n}ski effect. However, these measurements are challenging to model without robustly testing and verifying the theoretical framework used for analyzing them. Here, we leverage the accuracy and volume of the $N$-body simulation suite \textsc{AbacusSummit} to generate high-resolution Ly$\alpha$ skewers and quasi-stellar object (QSO) catalogs. One of the main goals of our mocks is to aid in the full-shape Ly$\alpha$ analysis planned by the Dark Energy Spectroscopic Instrument (DESI) team. We provide optical depth skewers for six of the fiducial cosmology base-resolution simulations ($L_{\rm box} = 2\,h^{-1}{\rm Gpc}$, $N = 6912^3$) at $z = 2.5$. We adopt a simple recipe based on the Fluctuating Gunn-Peterson Approximation (FGPA) for constructing these skewers from the matter density in an $N$-body simulation and calibrate it against the 1D and 3D Ly$\alpha$ power spectra extracted from the hydrodynamical simulation IllustrisTNG (TNG; $L_{\rm box} = 205\,h^{-1}{\rm Mpc}$, $N = 2500^3$). As an important application, we study the non-linear broadening of the baryon acoustic oscillation (BAO) peak and show the cross-correlation between DESI-like QSOs and our Ly$\alpha$ forest skewers. We find differences on small scales between the Kaiser approximation prediction and our mock measurements of the Ly$\alpha$$\times$QSO cross-correlation, which would be important to account for in upcoming analyses. The \textsc{AbacusSummit} Ly$\alpha$ forest mocks open up the possibility for improved modelling of cross correlations between Ly$\alpha$ and cosmic microwave background (CMB) lensing and Ly$\alpha$ and QSOs, and for forecasts of the 3-point Ly$\alpha$ correlation function. Our catalogues and skewers are publicly available on Globus., Comment: 17 pages, 8 figures, Globus link: https://app.globus.org/file-manager?origin_id=9ce29982-eed1-11ed-9bb4-c9bb788c490e&path=%2F

We present the first detection of the baryon acoustic oscillations (BAO) signal obtained using unblinded data collected during the initial two 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% completeness level, we report a ~5 sigma level BAO detection and the measurement of the BAO location at a precision of 1.7%. 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% completeness, we also detect the BAO feature at ~3 sigma significance with a 2.6% precision. These first BAO measurements represent an important milestone, acting as a quality control on the optimal performance of the complex robotically-actuated, fiber-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-percent precision with the completed 5-year 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., Comment: 17 pages, 10 figures, 4 tables. MNRAS, 525, 5406. Author accepted manuscript matching the published version

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a 5-year survey that aims to measure the redshifts of tens of millions of galaxies and quasars, including 8 million luminous red galaxies (LRGs) in the redshift range of $0.4

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 ($c_2=b_2/b_1$), and demonstrate why this is the case. For a $\sim$1500 sq. deg. galaxy survey we recover $c_2$ with an error of approximately $0.1$ for a wide range of $c_2$ 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 $c_2$ 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., Comment: 10 pages, 5 figures

We investigate the possibility to detect primordial non-Gaussianity by analysing the bulk of the probability distribution function (PDF) of late-time cosmic density fluctuations. For this purpose we devise a new method to predict the impact of general non-Gaussian initial conditions on the late-time density PDF. At redshift $z=1$ and for a smoothing scale of 30Mpc/$h$ our predictions agree with the high-resolution Quijote N-body simulations to $\sim 0.2\%$ precision. This is within cosmic variance of a $\sim 100(\mathrm{Gpc}/h)^3$ survey volume. When restricting to this 30Mpc/$h$ smoothing scale and to mildly non-linear densities ($\delta[30\mathrm{Mpc}/h] \in [-0.3, 0.4]$) and also marginalizing over potential ignorance of the amplitude of the non-linear power spectrum an analysis of the PDF for such a survey volume can still measure the amplitude of different primordial bispectrum shapes to an accuracy of \smash{$\Delta f_{\mathrm{NL}}^{\mathrm{loc}}=\pm 7.4\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{equi}}=\pm 22.0\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{ortho}}=\pm 46.0$} . When pushing to smaller scales and assuming a joint analysis of the PDF with smoothing radii of 30Mpc/$h$ and 15Mpc/$h$ ($\delta[15\mathrm{Mpc}/h] \in [-0.4, 0.5]$) this improves to \smash{$\Delta f_{\mathrm{NL}}^{\mathrm{loc}}=\pm 3.3\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{equi}}=\pm 11.0\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{ortho}}=\pm 17.0\ $} - even when marginalizing over the non-linear variances at both scales as two free parameters. Especially, such an analysis could simultaneously measure $f_{\mathrm{NL}}$ and the amplitude and slope of the non-linear power spectrum. However, at 15Mpc/$h$ our predictions are only accurate to $\lesssim 0.8\%$ for the considered density range. We discuss how this has to be improved in order to push to these small scales and make full use of upcoming surveys with a PDF-based analysis., Comment: changes to agree with journal version; 22 pages, 12 figures; public code available at https://github.com/OliverFHD/CosMomentum

We present the GEOMAX algorithm and its Python implementation for a two-step compression of bispectrum measurements. The first step groups bispectra by the geometric properties of their arguments; the second step then maximises the Fisher information with respect to a chosen set of model parameters in each group. The algorithm only requires the derivatives of the data vector with respect to the parameters and a small number of mock data, producing an effective, non-linear compression. By applying GEOMAX to bispectrum monopole measurements from BOSS DR12 CMASS redshift-space galaxy clustering data, we reduce the $68\%$ credible intervals for the inferred parameters $\left(b_1,b_2,f,\sigma_8\right)$ by $\left(50.4\%,56.1\%,33.2\%,38.3\%\right)$ with respect to standard MCMC on the full data vector. We run the analysis and comparison between compression methods over one hundred galaxy mocks to test the statistical significance of the improvements. On average GEOMAX performs $\sim15\%$ better than geometrical or maximal linear compression alone and is consistent with being lossless. Given its flexibility, the GEOMAX approach has the potential to optimally exploit three-point statistics of various cosmological probes like weak lensing or line-intensity maps from current and future cosmological data-sets such as DESI, Euclid, PFS and SKA., Comment: 17 pages, 9 figures, accepted version by MNRAS

We investigate various phenomenological schemes for the rapid generation of 3D mock galaxy catalogues with a given power spectrum and bispectrum. We apply the fast bispectrum estimator \MODALLSS{} to these mock galaxy catalogues and compare to $N$-body simulation data analysed with the halo-finder \texttt{ROCKSTAR} (our benchmark data). We propose an assembly bias model for populating parent halos with subhalos by using a joint lognormal-Gaussian probability distribution for the subhalo occupation number and the halo concentration. This prescription enabled us to recover the benchmark power spectrum from $N$-body simulations to within 1\% and the bispectrum to within 4\% across the entire range of scales of the simulation. A small further boost adding an extra galaxy to all parent halos above the mass threshold $M>2\times10^{14}\,h^{-1} M_\odot$ obtained a better than 1\% fit to both power spectrum and bispectrum in the range $K/3<1.1\,h\,\text{Mpc}^{-1}$, where $K=k_1+k_2+k_3$. This statistical model should be applicable to fast dark matter codes, allowing rapid generation of mock catalogues which simultaneously reproduce the halo power spectrum and bispectrum obtained from $N$-body simulations. We also investigate alternative schemes using the Halo Occupation Distribution (HOD) which depend only on halo mass, but these yield results deficient in both the power spectrum (2\%) and the bispectrum (>4\%) at $k,K/3 \approx 0.2\,h\,\text{Mpc}^{-1}$, with poor scaling for the latter. Efforts to match the power spectrum by modifying the standard four-parameter HOD model result in overboosting the bispectrum (with a 10\% excess). We also characterise the effect of changing the halo profile on the power spectrum and bispectrum.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov., Comment: 9-page APC White Paper submitted to the Astro2020 Decadal Survey. To be published in BAAS. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov

We present a novel method to compress galaxy clustering three-point statistics and apply it to redshift space galaxy bispectrum monopole measurements from BOSS DR12 CMASS data considering a $k$-space range of $0.03-0.12\,h/\mathrm{Mpc}$. The method consists in binning together bispectra evaluated at sets of wave-numbers forming closed triangles with similar geometrical properties: the area, the cosine of the largest angle and the ratio between the cosines of the remaining two angles. This enables us to increase the number of bispectrum measurements for example by a factor of $23$ over the standard binning (from 116 to 2734 triangles used), which is otherwise limited by the number of mock catalogues available to estimate the covariance matrix needed to derive parameter constraints. The $68\%$ credible intervals for the inferred parameters $\left(b_1,b_2,f,\sigma_8\right)$ are thus reduced by $\left(-39\%,-49\%,-29\%,-22\%\right)$, respectively. We find very good agreement with the posteriors recently obtained by alternative maximal compression methods. This new method does not require the a-priori computation of the data-vector covariance matrix and has the potential to be directly applicable to other three-point statistics (e.g. glaxy clustering, weak gravitational lensing, 21 cm emission line) measured from future surveys such as DESI, Euclid, PFS and SKA., Comment: 5 pages, 2 figures, Accepted by MNRAS: Letters

We investigate the impact of prior models on the upper bound of the sum of neutrino masses, $\sum m_{\nu}$. We use data from Large Scale Structure of galaxies, Cosmic Microwave Background, Type Ia SuperNovae, and Big Bang Nucleosynthesis. We probe physically motivated neutrino mass models (respecting oscillation experiment constraints) and compare them to constraints using standard cosmological approximations. The former give a consistent upper bound of $\sum m_{\nu} \lesssim 0.26$ eV ($95\%$ CI) and yields a strong competitive upper bound for the lightest neutrino mass species, $m_0^{\nu} < 0.086$ eV ($95\%$ CI). By contrast one of the approximations, which is somewhat inconsistent with oscillation experiments, yields an upper bound of $\sum m_{\nu} \lesssim 0.15$ eV ($95\%$ CI), which differs substantially from the former upper bound. We, therefore, argue that cosmological neutrino mass and hierarchy determination should be pursued using physically motivated models since approximations might lead to incorrect and nonphysical upper bounds., Comment: 5 pages, 2 figures

We constrain cosmological parameters by analysing the angular power spectra of the Baryon Oscillation Spectroscopic Survey DR12 galaxies, a spectroscopic follow-up of around 1.3 million SDSS galaxies over 9,376 deg$^2$ with an effective volume of $\sim 6.5$ (Gpc $h^{-1}$)$^3$ in the redshift range $0.15 \leq z < 0.80$. We split this sample into 13 tomographic bins ($\Delta z = 0.05$); angular power spectra were calculated using a Pseudo-$C_{\ell}$ estimator, and covariance matrices were estimated using log-normal simulated maps. Cosmological constraints obtained from these data were combined with constraints from Planck CMB experiment as well as the JLA supernovae compilation. Considering a $w$CDM cosmological model measured on scales up to $k_{max} = 0.07h$ Mpc$^{-1}$, we constrain a constant dark energy equation-of-state with a $\sim 4\%$ error at the 1-$\sigma$ level: $w_0 = -0.993^{+0.046}_{-0.043}$, together with $\Omega_m = 0.330\pm 0.012$, $\Omega_b = 0.0505 \pm 0.002$, $S_8 \equiv \sigma_8 \sqrt{\Omega_m/0.3} = 0.863 \pm 0.016$, and $h = 0.661 \pm 0.012$. For the same combination of datasets, but now considering a $\Lambda$CDM model with massive neutrinos and the same scale cut, we find: $\Omega_m = 0.328 \pm 0.009$, $\Omega_b = 0.05017^{+0.0009}_{-0.0008}$, $S_8 = 0.862 \pm 0.017$, and $h = 0.663^{+0.006}_{-0.007}$ and a 95\% credible interval (CI) upper limit of $\sum m_{\nu} < 0.14$ eV for a normal hierarchy. These results are competitive if not better than standard analyses with the same dataset, and demonstrate this should be a method of choice for future surveys, opening the door for their full exploitation in cross-correlations probes., Comment: 28 pages, 28 figures, 3 appendices

We apply two compression methods to the galaxy power spectrum monopole/quadrupole and bispectrum monopole measurements from the BOSS DR12 CMASS sample. Both methods reduce the dimension of the original data-vector to the number of cosmological parameters considered, using the Karhunen-Lo\`eve algorithm with an analytic covariance model. In the first case, we infer the posterior through MCMC sampling from the likelihood of the compressed data-vector (MC-KL). The second, faster option, works by first Gaussianising and then orthogonalising the parameter space before the compression; in this option (G-PCA) we only need to run a low-resolution preliminary MCMC sample for the Gaussianization to compute our posterior. Both compression methods accurately reproduce the posterior distributions obtained by standard MCMC sampling on the CMASS dataset for a $k$-space range of $0.03-0.12\,h/\mathrm{Mpc}$. The compression enables us to increase the number of bispectrum measurements by a factor of $\sim 23$ over the standard binning (from 116 to 2734 triangles used), which is otherwise limited by the number of mock catalogues available. This reduces the $68\%$ credible intervals for the parameters $\left(b_1,b_2,f,\sigma_8\right)$ by $\left(-24.8\%,-52.8\%,-26.4\%,-21\%\right)$, respectively. The best-fit values we obtain are $(b_1=2.31\pm0.17,b_2=0.77\pm0.19,$ $f(z_{\mathrm{CMASS}})=0.67\pm0.06,\sigma_8(z_{\mathrm{CMASS}})=0.51\pm0.03)$. Using these methods for future redshift surveys like DESI, Euclid and PFS will drastically reduce the number of simulations needed to compute accurate covariance matrices and will facilitate tighter constraints on cosmological parameters., Comment: 21 pages, 9 figures, Accepted by MNRAS after addressing the comments by the referee

We explore two methods of compressing the redshift space galaxy power spectrum and bispectrum with respect to a chosen set of cosmological parameters. Both methods involve reducing the dimension of the original data-vector ( e.g. 1000 elements ) to the number of cosmological parameters considered ( e.g. seven ) using the Karhunen-Lo\`eve algorithm. In the first case, we run MCMC sampling on the compressed data-vector in order to recover the one-dimensional (1D) and two-dimensional (2D) posterior distributions. The second option, approximately 2000 times faster, works by orthogonalising the parameter space through diagonalisation of the Fisher information matrix before the compression, obtaining the posterior distributions without the need of MCMC sampling. Using these methods for future spectroscopic redshift surveys like DESI, EUCLID and PFS would drastically reduce the number of simulations needed to compute accurate covariance matrices with minimal loss of constraining power. We consider a redshift bin of a DESI-like experiment. Using the power spectrum combined with the bispectrum as a data-vector, both compression methods on average recover the 68% credible regions to within 0.7% and 2% of those resulting from standard MCMC sampling respectively. These confidence intervals are also smaller than the ones obtained using only the power spectrum by (81%, 80%, 82%) respectively for the bias parameter b_1, the growth rate f and the scalar amplitude parameter A_s., Comment: 27 pages, 8 figures, 1 table, Accepted 2018 January 28. Received 2018 January 25; in original form 2017 September 11. Added clarifications in the text on the bias modelling and compression limits following referee's comments. Removed tetraspectrum term from the pk-bk cross covariance + correction in the appendix