Search

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)

We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We validate this approach on simulations reproducing the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) with and without fiber assignment. We show that we recover unbiased cosmological constraints using small angular scale truncated estimators from simulations with fiber assignment incompleteness, with respect to standard estimators from complete simulations. Additionally, we present an approach to remove the sensitivity of the fits to high $k$ modes in the theoretical power spectrum, by applying a transformation to the data vector and window matrix. We find that our method efficiently mitigates the effect of fiber assignment incompleteness in two-point correlation function and power spectrum measurements, at low computational cost and with little statistical loss., Comment: 36 pages, 23 figures, typos corrected, clarifications added

We present the technical details on how large-scale structure (LSS) catalogs are constructed from redshifts measured from spectra observed by the Dark Energy Spectroscopic Instrument (DESI). The LSS catalogs provide the information needed to determine the relative number density of DESI tracers as a function of redshift and celestial coordinates and, e.g., determine clustering statistics. We produce catalogs that are weighted subsamples of the observed data, each matched to a weighted `random' catalog that forms an unclustered sampling of the probability density that DESI could have observed those data at each location. Precise knowledge of the DESI observing history and associated hardware performance allows for a determination of the DESI footprint and the number of times DESI has covered it at sub-arcsecond level precision. This enables the completeness of any DESI sample to be modeled at this same resolution. The pipeline developed to create LSS catalogs has been designed to easily allow robustness tests and enable future improvements. We describe how it allows ongoing work improving the match between galaxy and random catalogs, such as including further information when assigning redshifts to randoms, accounting for fluctuations in target density, accounting for variation in the redshift success rate, and accommodating blinding schemes., Comment: Accepted (by JCAP) version of supporting publication of DESI 2024II: Sample definitions, characteristics, and two-point clustering statistics

In preparation for the next generation of galaxy redshift surveys, and in particular the year-one data release from the Dark Energy Spectroscopic Instrument (DESI), we investigate the consistency of a variety of effective field theory models that describe the galaxy-galaxy power spectra in redshift space into the quasi-linear regime using 1-loop perturbation theory. These models are employed in the pipelines \texttt{velocileptors}, \texttt{PyBird}, and \texttt{Folps$\nu$}. While these models have been validated independently, a detailed comparison with consistent choices has not been attempted. After briefly discussing the theoretical differences between the models we describe how to provide a more apples-to-apples comparison between them. We present the results of fitting mock spectra from the \texttt{AbacusSummit} suite of N-body simulations provided in three redshift bins to mimic the types of dark time tracers targeted by the DESI survey. We show that the theories behave similarly and give consistent constraints in both the forward-modeling and ShapeFit compressed fitting approaches. We additionally generate (noiseless) synthetic data from each pipeline to be fit by the others, varying the scale cuts in order to show that the models agree within the range of scales for which we expect 1-loop perturbation theory to be applicable. This work lays the foundation of Full-Shape analysis with DESI Y1 galaxy samples where in the tests we performed, we found no systematic error associated with the modeling of the galaxy redshift space power spectrum for this volume., Comment: 33 pages, 12 figures. Supporting publication of DESI 2024 V: Analysis of the full shape of two-point clustering statistics from galaxies and quasars

The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we focus on the robustness of the BAO analysis against the Halo Occupation Distribution (HOD) modeling for the Emission Line Galaxy (ELG) tracer. Based on a common dark matter simulation, our analysis relies on HOD mocks tuned to early DESI data, namely the One-Percent survey data. To build the mocks, we use several HOD models for the ELG tracer as well as extensions to the baseline HOD models. Among these extensions, we consider distinct recipes for galactic conformity and assembly bias. We perform two independent analyses in the Fourier space and in the configuration space. We recover the BAO signal from two-point measurements after performing reconstruction on our mocks. Additionally, we also apply the control variates technique to reduce sample variance noise. Our BAO analysis can recover the isotropic BAO parameter $\alpha_\text{iso}$ within 0.1\% and the Alcock Paczynski parameter $\alpha_\text{AP}$ within 0.3\%. Overall, we find that our systematic error due to the HOD dependence is below 0.17\%, with the Fourier space analysis being more robust against the HOD systematics. We conclude that our analysis pipeline is robust enough against the HOD systematics for the ELG tracer in the DESI 2024 BAO analysis., Comment: Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

In this paper, we present the estimation of systematics related to the halo occupation distribution (HOD) modeling in the baryon acoustic oscillations (BAO) distance measurement of the Dark Energy Spectroscopic Instrument (DESI) 2024 analysis. This paper focuses on the study of HOD systematics for luminous red galaxies (LRG). We consider three different HOD models for LRGs, including the base 5-parameter vanilla model and two extensions to it, that we refer to as baseline and extended models. The baseline model is described by the 5 vanilla HOD parameters, an incompleteness factor and a velocity bias parameter, whereas the extended one also includes a galaxy assembly bias and a satellite profile parameter. We utilize the 25 dark matter simulations available in the AbacusSummit simulation suite at $z=$ 0.8 and generate mock catalogs for our different HOD models. To test the impact of the HOD modeling in the position of the BAO peak, we run BAO fits for all these sets of simulations and compare the best-fit BAO-scaling parameters $\alpha_{\rm iso}$ and $\alpha_{\rm AP}$ between every pair of HOD models. We do this for both Fourier and configuration spaces independently, using post-reconstruction measurements. We find a 3.3$\sigma$ detection of HOD systematic for $\alpha_{\rm AP}$ in configuration space with an amplitude of 0.19%. For the other cases, we did not find a 3$\sigma$ detection, and we decided to compute a conservative estimation of the systematic using the ensemble of shifts between all pairs of HOD models. By doing this, we quote a systematic with an amplitude of 0.07% in $\alpha_{\rm iso}$ for both Fourier and configuration spaces; and of 0.09% in $\alpha_{\rm AP}$ for Fourier space., Comment: 36 pages, 9 figures. Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

We present an optimized way of producing the fast semi-analytical covariance matrices for the Legendre moments of the two-point correlation function, taking into account survey geometry and mimicking the non-Gaussian effects. We validate the approach on simulated (mock) catalogs for different galaxy types, representative of the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, used in 2024 analyses. We find only a few percent differences between the mock sample covariance matrix and our results, which can be expected given the approximate nature of the mocks, although we do identify discrepancies between the shot-noise properties of the DESI fiber assignment algorithm and the faster approximation used in the mocks. Importantly, we find a close agreement (<~ 5% relative differences) in the projected errorbars for distance scale parameters for the baryon acoustic oscillation measurements. This confirms our method as an attractive alternative to simulation-based covariance matrices, especially for non-standard models or galaxy sample selections, in particular, relevant to the broad current and future analyses of DESI data., Comment: Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars (arXiv:2404.03000). 29 pages, 4 figures. Prepared for submission to JCAP. Code available at https://github.com/oliverphilcox/RascalC and https://github.com/misharash/RascalC-scripts/tree/DESI2024. Data points from the plots available at https://zenodo.org/doi/10.5281/zenodo.10895161

We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$\alpha$ forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range $0.1-1$ and $w_a<0$. This preference is 2.6$\sigma$ for the DESI+CMB combination, and persists or grows when SN~Ia are added in, giving results discrepant with the $\Lambda$CDM model at the $2.5\sigma$, $3.5\sigma$ or $3.9\sigma$ levels for the addition of Pantheon+, Union3, or DES-SN5YR datasets respectively. For the flat $\Lambda$CDM model with the sum of neutrino mass $\sum m_\nu$ free, combining the DESI and CMB data yields an upper limit $\sum m_\nu < 0.072$ $(0.113)$ eV at 95% confidence for a $\sum m_\nu>0$ $(\sum m_\nu>0.059)$ eV prior. These neutrino-mass constraints are substantially relaxed in models beyond $\Lambda$CDM. [Abridged.], Comment: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers). Typos corrected and a new figure and discussion added to Appendix A

We present the measurement of Baryon Acoustic Oscillations (BAO) from the Lyman-$\alpha$ (Ly$\alpha$) forest of high-redshift quasars with the first-year dataset of the Dark Energy Spectroscopic Instrument (DESI). Our analysis uses over $420\,000$ Ly$\alpha$ forest spectra and their correlation with the spatial distribution of more than $700\,000$ quasars. An essential facet of this work is the development of a new analysis methodology on a blinded dataset. We conducted rigorous tests using synthetic data to ensure the reliability of our methodology and findings before unblinding. Additionally, we conducted multiple data splits to assess the consistency of the results and scrutinized various analysis approaches to confirm their robustness. For a given value of the sound horizon ($r_d$), we measure the expansion at $z_{\rm eff}=2.33$ with 2\% precision, $H(z_{\rm eff}) = (239.2 \pm 4.8) (147.09~{\rm Mpc} /r_d)$ km/s/Mpc. Similarly, we present a 2.4\% measurement of the transverse comoving distance to the same redshift, $D_M(z_{\rm eff}) = (5.84 \pm 0.14) (r_d/147.09~{\rm Mpc})$ Gpc. Together with other DESI BAO measurements at lower redshifts, these results are used in a companion paper to constrain cosmological parameters., Comment: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers). Minor changes in v4, version accepted for publication in JCAP

Baryon acoustic oscillations (BAO) provide a robust standard ruler to measure the expansion history of the Universe through galaxy clustering. Density-field reconstruction is now a widely adopted procedure for increasing the precision and accuracy of the BAO detection. With the goal of finding the optimal reconstruction settings to be used in the DESI 2024 galaxy BAO analysis, we assess the sensitivity of the post-reconstruction BAO constraints to different choices in our analysis configuration, performing tests on blinded data from the first year of DESI observations (DR1), as well as on mocks that mimic the expected clustering and selection properties of the DESI DR1 target samples. Overall, we find that BAO constraints remain robust against multiple aspects in the reconstruction process, including the choice of smoothing scale, treatment of redshift-space distortions, fiber assignment incompleteness, and parameterizations of the BAO model. We also present a series of tests that DESI followed in order to assess the maturity of the end-to-end galaxy BAO pipeline before the unblinding of the large-scale structure catalogs., Comment: Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

We present the DESI 2024 galaxy and quasar baryon acoustic oscillations (BAO) measurements using over 5.7 million unique galaxy and quasar redshifts in the range 0.1

In geophysics, volcanoes are particularly difficult to image because of the multi-scale heterogeneities of fluids and rocks that compose them and their complex non-linear dynamics. By exploiting seismic noise recorded by a sparse array of geophones, we are able to reveal the magmatic and hydrothermal plumbing system of La Soufri\`ere volcano in Guadeloupe. Spatio-temporal cross-correlation of seismic noise actually provides the impulse responses between virtual geophones located inside the volcano. The resulting reflection matrix can be exploited to numerically perform an auto-focus of seismic waves on any reflector of the underground. An unprecedented view on the volcano's inner structure is obtained at a half-wavelength resolution. This innovative observable provides fundamental information for the conceptual modeling and high-resolution monitoring of volcanoes., Comment: 32 pages, 9 figures

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.

The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes good-quality spectral information from 466,447 objects targeted as part of the Milky Way Survey, 428,758 as part of the Bright Galaxy Survey, 227,318 as part of the Luminous Red Galaxy sample, 437,664 as part of the Emission Line Galaxy sample, and 76,079 as part of the Quasar sample. In addition, the release includes spectral information from 137,148 objects that expand the scope beyond the primary samples as part of a series of secondary programs. Here, we describe the spectral data, data quality, data products, Large-Scale Structure science catalogs, access to the data, and references that provide relevant background to using these spectra., Comment: 43 pages, 7 figures, 17 tables, accepted for publication in the Astronomical Journal

The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg$^2$ using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg$^2$ program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval $z<1.1$, 0.39% over the redshift interval $1.1

The halo occupation distribution (HOD) framework is an empirical method to describe the connection between dark matter halos and galaxies, which is constrained by small scale clustering data. Efficient fitting procedures are required to scan the HOD parameter space. This paper describes such a method based on Gaussian Processes to iteratively build a surrogate model of the posterior of the likelihood surface from a reasonable amount of likelihood computations, typically two orders of magnitude less than standard Monte Carlo Markov chain algorithms. Errors in the likelihood computation due to stochastic HOD modelling are also accounted for in the method we propose. We report results of reproducibility, accuracy and stability tests of the method derived from simulation, taking as a test case star-forming emission line galaxies, which constitute the main tracer of the Dark Energy Spectroscopic Instrument and have so far a poorly constrained galaxy-halo connection from observational data.

Standard full-shape clustering analyses in Fourier space rely on a fixed power spectrum template, defined at the fiducial cosmology used to convert redshifts into distances, and compress the cosmological information into the Alcock-Paczynski parameters and the linear growth rate of structure. In this paper, we propose an analysis method that operates directly in the cosmology parameter space and varies the power spectrum template accordingly at each tested point. Predictions for the power spectrum multipoles from the TNS model are computed at different cosmologies in the framework of $\Lambda \rm{CDM}$. Applied to the final eBOSS QSO and LRG samples together with the low-z DR12 BOSS galaxy sample, our analysis results in a set of constraints on the cosmological parameters $\Omega_{\rm cdm}$, $H_0$, $\sigma_8$, $\Omega_{\rm b}$ and $n_s$. To reduce the number of computed models, we construct an iterative process to sample the likelihood surface, where each iteration consists of a Gaussian process regression. This method is validated with mocks from N-body simulations. From the combined analysis of the (e)BOSS data, we obtain the following constraints: $\sigma_8=0.877\pm 0.049$ and $\Omega_{\rm m}=0.304^{+0.016}_{-0.010}$ without any external prior. The eBOSS quasar sample alone shows a $3.1\sigma$ discrepancy compared to the Planck prediction., Comment: 13 pages, 7 figures

Galaxy redshift surveys are subject to incompleteness and inhomogeneous sampling due to the various constraints inherent to spectroscopic observations. This can introduce systematic errors on the summary statistics of interest, which need to be mitigated in cosmological analysis to achieve high accuracy. Standard practices involve applying weighting schemes based on completeness estimates across the survey footprint, possibly supplemented with additional weighting schemes accounting for density-dependent effects. In this work, we concentrate on pure angular systematics and describe an alternative approach consisting in analysing the galaxy two-point correlation function where angular modes are nulled. By construction, this procedure removes all possible known and unknown sources of angular observational systematics, but also part of the cosmological signal.We use a modified Landy-Szalay estimator for the two-point correlation function that relies on an additional random catalogue where angular positions are randomly drawn from the galaxy catalogue, and provide an analytical model to describe this modified statistic. We test the model by performing an analysis of the full anisotropic clustering in mock catalogues of luminous red and emission-line galaxies at 0.43 < z < 1.1. We find that the model fully accounts for the modified correlation function in redshift space, without introducing new nuisance parameters. The derived cosmological parameters from the analysis of baryon acoustic oscillations and redshift-space distortions display slightly larger statistical uncertainties, mostly for the growth rate of structure parameter fs8 that exhibits a 50% statistical error increase, but free from angular systematic error., Comment: 16 pages, submitted to MNRAS

We construct cosmic void catalogues with the DIVE void finder upon SDSS BOSS DR12 and eBOSS DR16 galaxy samples with BAO reconstruction applied, and perform a joint BAO analysis using different types of galaxies and the corresponding voids. The BAO peak is evident for the galaxy-galaxy, galaxy-void, and void-void correlation functions of all datasets, including the ones cross correlating luminous red galaxy and emission line galaxy samples. Two multi-tracer BAO fitting schemes are then tested, one combining the galaxy and void correlation functions with a weight applied to voids, and the other using a single BAO dilation parameter for all clustering measurements of different tracers. Both methods produce consistent results with mock catalogues, and on average ~10 per cent improvements of the BAO statistical uncertainties are observed for all samples, compared to the results from galaxies alone. By combining the clustering of galaxies and voids, the uncertainties of BAO measurements from the SDSS data are reduced by 5 to 15 per cent, yielding 0.9, 0.8, 1.1, 2.3, and 2.9 per cent constraints on the distance $D_{_{\rm V}}(z)$, at effective redshifts of 0.38, 0.51, 0.70, 0.77, and 0.85, respectively. When combined with BAO measurements from SDSS MGS, QSO, and Ly$\alpha$ samples, as well as the BBN results, we obtain $H_0 = 67.58 \pm 0.91\,{\rm km}\,{\rm s}^{-1}\,{\rm Mpc}^{-1}$, $\Omega_{\rm m} = 0.290 \pm 0.015$, and $\Omega_\Lambda h^2 = 0.3241 \pm 0.0079$ in the flat-$\Lambda$CDM framework, where the 1$\,\sigma$ uncertainties are around 6, 6, and 17 per cent smaller respectively, compared to constraints from the corresponding anisotropic BAO measurements without voids and LRG-ELG cross correlations., Comment: 33 pages, 30 figures, published on MNRAS

We measure the small-scale clustering of the Data Release 16 extended Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy sample, corrected for fibre-collisions using Pairwise Inverse Probability weights, which give unbiased clustering measurements on all scales. We fit to the monopole and quadrupole moments and to the projected correlation function over the separation range $7-60\,h^{-1}$Mpc with a model based on the Aemulus cosmological emulator to measure the growth rate of cosmic structure, parameterized by $f\sigma_8$. We obtain a measurement of $f\sigma_8(z=0.737)=0.408\pm0.038$, which is $1.4\sigma$ lower than the value expected from 2018 Planck data for a flat $\Lambda$CDM model, and is more consistent with recent weak-lensing measurements. The level of precision achieved is 1.7 times better than more standard measurements made using only the large-scale modes of the same sample. We also fit to the data using the full range of scales $0.1-60\,h^{-1}$Mpc modelled by the Aemulus cosmological emulator and find a $4.5\sigma$ tension in the amplitude of the halo velocity field with the Planck+$\Lambda$CDM model, driven by a mismatch on the non-linear scales. This may not be cosmological in origin, and could be due to a breakdown in the Halo Occupation Distribution model used in the emulator. Finally, we perform a robust analysis of possible sources of systematics, including the effects of redshift uncertainty and incompleteness due to target selection that were not included in previous analyses fitting to clustering measurements on small scales., Comment: 19 pages, 16 figures, published in MNRAS

We investigate the large-scale clustering of the final spectroscopic sample of quasars from the recently completed extended Baryon Oscillation Spectroscopic Survey (eBOSS). The sample contains $343708$ objects in the redshift range $0.8

Particularly in the presence of a hydrothermal system, many volcanoes output large quantities of heat through the transport of water from deep within the edifice to the surface. Thus, heat flux is a prime tool for evaluating volcanic activity and unrest. We review the volcanic unrest at La Soufri\`ere de Guadeloupe (French West Indies) using an airborne thermal camera survey, and in-situ measurements of temperature and flow rate through temperature probes, Pitot-tube and MultiGAS measurements. We deduce mass and heat fluxes for the fumarolic, ground and thermal spring outputs and follow these over a period spanning 2000--2020. Our results are compared with published data and we performed a retrospective analysis of the temporal variations in heat flux over this period using the literature data. We find that the heat emitted by the volcano is 36.5 +/- 7.9 MW, of which the fumarolic heat flux is dominant at 28.3 +/- 6.8 MW. Given a total heated area of 26780 m2, this equates to a heat flux density of 627 +/- 94 W/m2, which is amongst the highest established for worldwide volcanoes with hydrothermal systems, particularly for dome volcanoes. A major change at La Soufri\`ere de Guadeloupe, however, is the development of a widespread region of ground heating at the summit where heat output has increased from 0.2 +/- 0.1 MW in 2010 to 5.7 +/- 0.9 MW in 2020. This change is concurrent with accelerating unrest at the volcano, and the emergence of two new high-flux fumaroles in recent years. Our findings highlight the importance of continued and enhanced surveillance and research strategies at La Soufri\`ere de Guadeloupe, the results of which can be used to better understand hydrothermal volcanism the world over.

The forthcoming Dark Energy Spectroscopic Instrument (DESI) experiment plans to measure the effects of dark energy on the expansion of the Universe and create a $3$D map of the Universe using galaxies up to $z \sim 1.6$ and QSOs up to $z \sim 3.5$. In order to create this map, DESI will obtain spectroscopic redshifts of over $30$ million objects; among them, a majority are \oii emitting star-forming galaxies known as emission-line galaxies (ELGs). These ELG targets will be pre-selected by drawing a selection region on the $g - r$ vs. $r - z$ colour-colour plot, where high redshift ELGs form a separate locus from the lower redshift ELGs and interlopers. In this paper, we study the efficiency of three ELG target selection algorithms -- the final design report (FDR) cut based on the DEEP2 photometry, Number Density Modelling and Random Forest -- to determine how the combination of these three algorithms can be best used to yield a simple selection boundary that will be best suited to meet DESI's science goals. To do this, we selected $17$ small patches in the DESI footprint where we run the three target selection algorithms to pre-select ELGs based on their photometry. We observed the pre-selected ELGs using the MMT Binospec, which is similar in functionality to the DESI instrument, to obtain their spectroscopic redshifts and fluxes of $1054$ ELGs. By analysing the redshift and fluxing distribution of these galaxies, we find that although NDM performed the best, simple changes in the FDR definition would also yield sufficient performance., Comment: 16 pages, 9 figures, Accepted by the MNRAS on 2020 July 27

We test general relativity (GR) at the effective redshift $\bar{z} \sim 1.5$ by estimating the statistic $E_G$, a probe of gravity, on cosmological scales $19 - 190\,h^{-1}{\rm Mpc}$. This is the highest-redshift and largest-scale estimation of $E_G$ so far. We use the quasar sample with redshifts $0.8 < z < 2.2$ from Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16) as the large-scale structure (LSS) tracer, for which the angular power spectrum $C_\ell^{qq}$ and the redshift-space distortion (RSD) parameter $\beta$ are estimated. By cross correlating with the $\textit{Planck}$ 2018 cosmic microwave background (CMB) lensing map, we detect the angular cross-power spectrum $C_\ell^{\kappa q}$ signal at $12\,\sigma$ significance. Both jackknife resampling and simulations are used to estimate the covariance matrix (CM) of $E_G$ at $5$ bins covering different scales, with the later preferred for its better constraints on the covariances. We find $E_G$ estimates agree with the GR prediction at $1\,\sigma$ level over all these scales. With the CM estimated with $300$ simulations, we report a best-fit scale-averaged estimate of $E_G(\bar{z})=0.30\pm 0.05$, which is in line with the GR prediction $E_G^{\rm GR}(\bar{z})=0.33$ with $\textit{Planck}$ 2018 CMB+BAO matter density fraction $\Omega_{\rm m}=0.31$. The statistical errors of $E_G$ with future LSS surveys at similar redshifts will be reduced by an order of magnitude, which makes it possible to constrain modified gravity models., Comment: 16 pages, 14 figures; references added, matches version accepted by MNRAS

Accurate mock catalogues are essential for assessing systematics in the cosmological analysis of large galaxy surveys. Anisotropic two-point clustering measurements from the same simulation show some scatter for different lines of sight (LOS), but are on average equal, due to cosmic variance. This results in scatter in the measured cosmological parameters. We use the OuterRim N-body simulation halo catalogue to investigate this, considering the 3 simulation axes as LOS. The quadrupole of the 2-point statistics is particularly sensitive to changes in the LOS, with sub-percent level differences in the velocity distributions resulting in ~1.5$\sigma$ shifts on large scales. Averaging over multiple LOS can reduce the impact of cosmic variance. We derive an expression for the Gaussian cross-covariance between the power spectrum multipole measurements, for any two LOS, including shot noise, and the corresponding reduction in variance in the average measurement. Quadrupole measurements are anti-correlated, and for three orthogonal LOS, the variance on the average measurement is reduced by more than 1/3. We perform a Fisher analysis to predict the corresponding gain in precision on the cosmological parameter measurements, which we compare against a set of 300 extended Baryon Oscillation Spectroscopic Survey (eBOSS) emission line galaxy (ELG) EZmocks. The gain in $f\sigma_8$, which measures the growth of structure, is also better than 1/3. Averaging over multiple LOS in future mock challenges will allow the RSD models to be constrained with the same systematic error, with less than 3 times the CPU time., Comment: 14 pages, 10 figures, 1 table, updated to match the version accepted for publication in MNRAS

We measure the anisotropic clustering of the quasar sample from Data Release 16 (DR16) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (eBOSS). A sample of $343,708$ spectroscopically confirmed quasars between redshift $0.8

We measure the clustering of quasars of the final data release (DR16) of eBOSS. The sample contains $343\,708$ quasars between redshifts $0.8\leq z\leq2.2$ over $4699\,\mathrm{deg}^2$. We calculate the Legendre multipoles (0,2,4) of the anisotropic power spectrum and perform a BAO and a Full-Shape (FS) analysis at the effective redshift $z{\rm eff}=1.480$. The errors include systematic errors that amount to 1/3 of the statistical error. The systematic errors comprise a modelling part studied using a blind N-Body mock challenge and observational effects studied with approximate mocks to account for various types of redshift smearing and fibre collisions. For the BAO analysis, we measure the transverse comoving distance $D_{\rm M}(z_{\rm eff})/r_{\rm drag}=30.60\pm{0.90}$ and the Hubble distance $D_{\rm H}(z_{\rm eff})/r_{\rm drag}=13.34\pm{0.60}$. This agrees with the configuration space analysis, and the consensus yields: $D_{\rm M}(z_{\rm eff})/r_{\rm drag}=30.69\pm{0.80}$ and $D_{\rm H}(z_{\rm eff})/r_{\rm drag}=13.26\pm{0.55}$. In the FS analysis, we fit the power spectrum using a model based on Regularised Perturbation Theory, which includes Redshift Space Distortions and the Alcock-Paczynski effect. The results are $D_{\rm M}(z_{\rm eff})/r_{\rm drag}=30.68\pm{0.90}$ and $D_{\rm H}(z_{\rm eff})/r_{\rm drag}=13.52\pm{0.51}$ and we constrain the linear growth rate of structure $f(z_{\rm eff})\sigma_8(z_{\rm eff})=0.476\pm{0.047}$. Our results agree with the configuration space analysis. The consensus analysis of the eBOSS quasar sample yields: $D_{\rm M}(z_{\rm eff})/r_{\rm drag}=30.21\pm{0.79}$, $D_{\rm H}(z_{\rm eff})/r_{\rm drag}=3.23\pm{0.47}$ and $f(z_{\rm eff})\sigma_8(z_{\rm eff})=0.462\pm{0.045}$ and is consistent with a flat $\Lambda {\rm CDM}$ cosmological model using Planck results., Comment: 22 pages, 16 figures, 12 tables A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://sdss.org/science/final-bao-and-rsd-measurements/ . The full cosmological interpretation of these measurements can be found here: https://sdss.org/science/cosmology-results-from-eboss/