Search

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 arcsec from the centers of the galaxies, corresponding to an impact parameter of approximately 50 kpc. No unidentified spectral line shows up in the search, and we place a line flux limit of 10-19 ergs/s/cm2/arcsec2 on emissions in the wavelength range of 2000-9000A∘. This places the tightest constraints yet on the two-photon decay of dark matter in the mass range of 5 to 12 eV, with a particle lifetime exceeding 3×1025 s. This detection limit also implies that the line surface brightness contributed from all dark matter along the line of sight is at least 2 orders of magnitude lower than the measured extragalactic background light (EBL), ruling out the possibility that narrow optical-line emission from dark matter decay is a major source of the EBL.

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, N-body simulations, such as the abacuspng set presented in this study, consisting of 9 boxes, each of size Lbox=2 Gpc/h and particle mass of 1.01×1010M⊙/h, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory (PT) and the "separate universe"method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, bφ. As a science application, we investigate the link between halo assembly bias and bφ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between bφ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the bφ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasistellar object (QSO) samples and comment on the effect on fNL constraints for the allowed galaxy-halo models (note that σ[fNL]∝σ[bφ]bφ). We find that the error on fNL is 15, 8, 7 for the LRGs at z=0.5 and z=0.8 and QSOs at z=1.4, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

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

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

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

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

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a five-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 0.4 < z ≲ 1.0. Here we present the selection of the DESI LRG sample and assess its spectroscopic performance using data from Survey Validation (SV) and the first two months of the Main Survey. The DESI LRG sample, selected using g, r, z, and W1 photometry from the DESI Legacy Imaging Surveys, is highly robust against imaging systematics. The sample has a target density of 605 deg−2 and a comoving number density of 5 × 10−4 h3 Mpc−3 in 0.4 < z < 0.8; this is a significantly higher density than previous LRG surveys (such as SDSS, BOSS, and eBOSS) while also extending to z ∼ 1. After applying a bright star veto mask developed for the sample, 98.9% of the observed LRG targets yield confident redshifts (with a catastrophic failure rate of 0.2% in the confident redshifts), and only 0.5% of the LRG targets are stellar contamination. The LRG redshift efficiency varies with source brightness and effective exposure time, and we present a simple model that accurately characterizes this dependence. In the appendices, we describe the extended LRG samples observed during SV.

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

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

A system of 5,020 robotic fiber positioners was installed in 2019 on theMayall Telescope, at Kitt Peak National Observatory. The robots automaticallyre-target their optical fibers every 10 - 20 minutes, each to a precision ofseveral microns, with a reconfiguration time less than 2 minutes. Over the nextfive years, they will enable the newly-constructed Dark Energy SpectroscopicInstrument (DESI) to measure the spectra of 35 million galaxies and quasars.DESI will produce the largest 3D map of the universe to date and measure theexpansion history of the cosmos. In addition to the 5,020 robotic positionersand optical fibers, DESI's Focal Plane System includes 6 guide cameras, 4wavefront cameras, 123 fiducial point sources, and a metrology camera mountedat the primary mirror. The system also includes associated structural, thermal,and electrical systems. In all, it contains over 675,000 individual parts. Wediscuss the design, construction, quality control, and integration of all thesecomponents. We include a summary of the key requirements, the review andacceptance process, on-sky validations of requirements, and lessons learned forfuture multi-object, fiber-fed spectrographs.

This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 survey that publicly releases infrared spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the subsurvey Time Domain Spectroscopic Survey data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey subsurvey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated value-added catalogs. This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper, Local Volume Mapper, and Black Hole Mapper surveys.

We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical haonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $ {O}(N_math{g}2)$ for Ng particles, or $ {O}(N_math{FFT}log N_math{FFT})$ when using a Fast Fourier Transfo with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical haonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in $100$ CPU-hours, including the corrections necessary for non-unifo survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

We search for the baryon acoustic oscillations in the projected cross-correlation function binned into transverse comoving radius between the SDSS-IV DR16 eBOSS quasars and a dense photometric sample of galaxies selected from the DESI Legacy Imaging Surveys. We estimate the density of the photometric sample of galaxies in this redshift range to be about 2900 deg−2, which is deeper than the official DESI emission line galaxy selection, and the density of the spectroscopic sample is about 20 deg−2. In order to mitigate the systematics related to the use of different imaging surveys close to the detection limit, we use a neural network approach that accounts for complex dependences between the imaging attributes and the observed galaxy density. We find that we are limited by the depth of the imaging surveys that affects the density and purity of the photometric sample and its overlap in redshift with the quasar sample, which thus affects the performance of the method. When cross-correlating the photometric galaxies with quasars in the range 0.6 ≤ z ≤ 1.2, the cross-correlation function can provide better constraints on the comoving angular distance DM (6 per cent precision) compared to the constraint on the spherically averaged distance DV (9 per cent precision) obtained from the autocorrelation. Although not yet competitive, this technique will benefit from the arrival of deeper photometric data from upcoming surveys that will enable it to go beyond the current limitations we have identified in this work.

We present the steps taken to produce a reliable and complete input galaxy catalogue for the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) using the photometric Legacy Survey DR8 DECam. We analyse some of the main issues faced in the selection of targets for the DESI BGS, such as star-galaxy separation, contamination by fragmented stars and bright galaxies. Our pipeline utilizes a new way to select BGS galaxies using Gaia photometry and we implement geometrical and photometric masks that reduce the number of spurious objects. The resulting catalogue is cross-matched with the Galaxy And Mass Assembly (GAMA) survey to assess the completeness of the galaxy catalogue and the performance of the target selection. We also validate the clustering of the sources in our BGS catalogue by comparing with mock catalogues and the Sloan Digital Sky Survey (SDSS) data. Finally, the robustness of the BGS selection criteria is assessed by quantifying the dependence of the target galaxy density on imaging and other properties. The largest systematic correlation we find is a 7 per cent suppression of the target density in regions of high stellar density.

Abstract: The DESI Milky Way Survey (MWS) will observe ≥8 million stars between 16

Abstract: The Dark Energy Spectroscopic Instrument (DESI) will execute a nearly magnitude-limited survey of low redshift galaxies (0.05 ≤ z ≤ 0.4, median z ≈ 0.2). Clustering analyses of this Bright Galaxy Survey (BGS) will yield the most precise measurements to date of baryon acoustic oscillations and redshift-space distortions at low redshift. DESI BGS will comprise two target classes: (i) BRIGHT (r

Abstract: DESI will precisely constrain cosmic expansion and the growth of structure by collecting ∼35 million redshifts across ∼80% of cosmic history and one third of the sky to study Baryon Acoustic Oscillations (BAO) and Redshift Space Distortions (RSD). We present a preliminary target selection for an Emission Line Galaxy (ELG) sample, which will comprise about half of all DESI tracers. The selection consists of a g-band magnitude cut and a (g − r) versus (r − z) color box, which we validate using HSC/PDR2 photometric redshifts and DEEP2 spectroscopy. The ELG target density should be ∼2400 deg−2, with ∼65% of ELG redshifts reliably within a redshift range of 0.6

Abstract: The DESI survey will observe more than 8 million candidate luminous red galaxies (LRGs) in the redshift range 0.3

Abstract: The DESI survey will measure large-scale structure using quasars as direct tracers of dark matter in the redshift range 0.9  2.1. We present two methods to select candidate quasars for DESI based on imaging in three optical (g, r, z) and two infrared (W1, W2) bands. The first method uses traditional color cuts and the second utilizes a machine-learning algorithm.

The DESI Legacy Imaging Surveys (http://legacysurvey.org/) are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ≈14,000 deg2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12, and 22 μm) observed by the Wide-field Infrared Survey Explorer satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

Despite two decades of tremendous experimental and theoretical progress, theriddle of the accelerated expansion of the Universe remains to be solved. Onthe experimental side, our understanding of the possibilities and limitationsof the major dark energy probes has evolved; here we summarize the major probesand their crucial challenges. On the theoretical side, the taxonomy ofexplanations for the accelerated expansion rate is better understood, providingclear guidance to the relevant observables. We argue that: i) improvingstatistical precision and systematic control by taking more data, supportingresearch efforts to address crucial challenges for each probe, usingcomplementary methods, and relying on cross-correlations is well motivated; ii)blinding of analyses is difficult but ever more important; iii) studies of darkenergy and modified gravity are related; and iv) it is crucial that R&D for avibrant dark energy program in the 2030s be started now by supporting studiesand technical R&D that will allow embryonic proposals to mature. Understandingdark energy, arguably the biggest unsolved mystery in both fundamental particlephysics and cosmology, will remain one of the focal points of cosmology in theforthcoming decade.

The expansion of the Universe is understood to have accelerated during twoepochs: in its very first moments during a period of Inflation and much morerecently, at $z < 1$, when Dark Energy is hypothesized to drive cosmicacceleration. The undiscovered mechanisms behind these two epochs representsome of the most important open problems in fundamental physics. The largecosmological volume at $2 < z < 5$, together with the ability to efficientlytarget high-$z$ galaxies with known techniques, enables large gains in thestudy of Inflation and Dark Energy. A future spectroscopic survey can test theGaussianity of the initial conditions up to a factor of ~50 better than ourcurrent bounds, crossing the crucial theoretical threshold of$\sigma(f_{NL}^{\rm local})$ of order unity that separates single field andmulti-field models. Simultaneously, it can measure the fraction of Dark Energyat the percent level up to $z = 5$, thus serving as an unprecedented test ofthe standard model and opening up a tremendous discovery space.

Extreme emission line galaxies (EELGs) exhibit large equivalent widths (EW) in their rest-optical emission lines ([O  iii ] |$\lambda 5007$| or H  |$\alpha$| rest-frame EW |$\gt 750$|  Å) which can be tied to a recent upturn in star formation rate (SFR), due to the sensitivity of the nebular line emission and the rest-optical continuum to young (⁠|$\lt 10$|  Myr) and evolved stellar populations, respectively. By studying a sample of 85 star-forming galaxies (SFGs), spanning the redshift and magnitude interval |$3 \lt z\lt 9.5$| and |$-16\gt $| M |$\rm _{UV}\gt -21$|⁠ , in the JWST Advanced Deep Extragalactic Survey (JADES) with NIRSpec/prism spectroscopy, we determine that SFGs initiate an EELG phase when entering a significant burst of star formation, with the highest EWs observed in EELGs with the youngest luminosity-weighted ages (⁠|$\lt 5$| Myr) and the highest burst intensity (those with the greatest excess between their current and long-term average SFR). We spectroscopically confirm that a greater proportion of SFGs are in an EELG phase at high redshift in our UV-selected sample (⁠|$61\pm 4~{{\ \rm per\ cent}}$| in our |$z\gt 5.7$| high-redshift bin, compared to |$23^{+4}_{-1}\%$| in our lowest redshift bin |$3\lt z\lt 4.1$|⁠) due to the combined evolution of metallicity, ionization parameter, and star formation histories with redshift. We report that the EELGs within our sample exhibit a higher average ionization efficiency (⁠|$\log _{10}(\xi _\mathrm{ion}^\mathrm{HII}/{\rm erg^{-1}Hz}) =25.5\pm 0.2$|⁠) than the non-EELGs. High-redshift EELGs therefore comprise a population of efficient ionizing photon producers. Additionally, we report that 53 per cent (9/17) of EELGs at |$z\gt 5.7$| have observed Ly  |$\alpha$| emission, potentially lying within large ionized regions. The high detection rate of Ly  |$\alpha$| emitters in our EELG selection suggests that the physical conditions associated with entering an EELG phase also promote the escape of Ly  |$\alpha$| photons. [ABSTRACT FROM AUTHOR]

We describe new ultradeep James Webb Space Telescope (JWST) NIRSpec PRISM and grating spectra for the galaxies JADES-GS-z11-0 ( z spec = 11.122 − 0.003 + 0.005 ) and JADES-GS-z13-0 ( z spec = 13.20 − 0.04 + 0.03 ), the most distant spectroscopically confirmed galaxy discovered in the first year of JWST observations. The extraordinary depth of these observations (75 hr and 56 hr, respectively) provides a unique opportunity to explore the redshifts, stellar properties, UV magnitudes, and slopes for these two sources. For JADES-GS-z11-0, we find evidence for multiple emission lines, including [O ii ] λ λ 3726, 3729 and [Ne iii ] λ 3869, resulting in a spectroscopic redshift we determine with 94% confidence. We present stringent upper limits on the emission-line fluxes and line equivalent widths for JADES-GS-z13-0. At this spectroscopic redshift, the Ly α break in JADES-GS-z11-0 can be fit with a damped Ly α absorber with log (N HI / cm − 2) = 22.42 − 0.120 + 0.093 . These results demonstrate how neutral hydrogen fraction and Lyman-damping wings may impact the recovery of spectroscopic redshifts for sources like these, providing insight into the overprediction of the photometric redshifts seen for distant galaxies observed with JWST. In addition, we analyze updated NIRCam photometry to calculate the morphological properties of these resolved sources, and find a secondary source 0.″3 south of JADES-GS-z11-0 at a similar photometric redshift, hinting at how galaxies grow through interactions in the early Universe. [ABSTRACT FROM AUTHOR]

We present an algorithm enabling computation of the anisotropic redshift-space galaxy threepoint correlation function (3PCF) scaling as N2, with N the number of galaxies. Our previous work showed how to compute the isotropic 3PCF with this scaling by expanding the radially binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form multipole moments. The N2 scaling occurred because this approach never explicitly required the relative angle between a galaxy pair about the primary galaxy. Here, we generalize this work, demonstrating that in the presence of azimuthally symmetric anisotropy produced by redshift-space distortions (RSD), the 3PCF can be described by two triangle side lengths, two independent total angular momenta, and a spin. This basis for the anisotropic 3PCF allows its computation with negligible additional work over the isotropic 3PCF. We also present the covariance matrix of the anisotropic 3PCF measured in this basis. Our algorithm tracks the full 5D redshift-space 3PCF, uses an accurate line of sight to each triplet, is exact in angle, and easily handles edge correction. It will enable use of the anisotropic large-scale 3PCF as a probe of RSD in current and upcoming large-scale redshift surveys.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since 2014 July. This paper describes the second data release from this phase, and the 14th from SDSS overall (making this Data Release Fourteen or DR14). This release makes the data taken by SDSS-IV in its first two years of operation (2014-2016 July) public. Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey; the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data-driven machine-learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from the SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS web site (www.sdss.org) has been updated for this release and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020 and will be followed by SDSS-V.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147 000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts 0.8 < z < 2.2 and measure their spherically averaged clustering in both configuration and Fourier space. Our observational data set and the 1400 simulated realizations of the data set allow us to detect a preference for BAO that is greater than 2.8σ. We determine the spherically averaged BAO distance to z = 1.52 to 3.8 per cent precision: DV (z = 1.52) = 3843 ± 147 (rd/rd,fid) Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flatΛCDMbest-fitting cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find ΩΛ > 0 at 6.6s significance when testing a ΛCDM model with free curvature.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in 2014 July. It pursues three core programs: the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2), Mapping Nearby Galaxies at APO (MaNGA), and the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). As well as its core program, eBOSS contains two major subprograms: the Time Domain Spectroscopic Survey (TDSS) and the SPectroscopic IDentification of ERosita Sources (SPIDERS). This paper describes the first data release from SDSS-IV, Data Release 13 (DR13). DR13 makes publicly available the first 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA. It includes new observations from eBOSS, completing the Sloan Extended QUasar, Emission-line galaxy, Luminous red galaxy Survey (SEQUELS), which also targeted variability-selected objects and X-ray-selected objects. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification, and new reductions of the SDSS-III APOGEE-1 data, improving stellar parameters for dwarf stars and cooler stars. DR13 provides more robust and precise photometric calibrations. Value-added target catalogs relevant for eBOSS, TDSS, and SPIDERS and an updated red-clump catalog for APOGEE are also available. This paper describes the location and format of the data and provides references to important technical papers. The SDSS web site, http://www.sdss.org, provides links to the data, tutorials, examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ∼6 yr operations of SDSS-IV.

We present a measurement of baryon acoustic oscillations (BAO) in the cross-correlation of quasars with the Lyα-forest flux transmission at a mean redshift of z = 2.40. The measurement uses the complete Sloan Digital Sky Survey (SDSS-III) data sample: 168 889 forests and 234 367 quasars from the SDSS data release DR12. In addition to the statistical improvement on our previous study using DR11, we have implemented numerous improvements at the analysis level enabling a more accurate measurement of this cross-correlation. We have also developed the first simulations of the cross-correlation that allow us to test different aspects of our data analysis and to search for potential systematic errors in the determination of the BAO peak position. We measure the two ratios DH(z = 2.40) /rd = 9.01 ± 0.36 and DM(z = 2.40) /rd = 35.7 ± 1.7, where the errors include marginalization over the non-linear velocity of quasars and the cross-correlation of metals and quasars, among other effects. These results are within 1.8σ of the prediction of the flat-ΛCDM model describing the observed cosmic microwave background anisotropies. We combine this study with the Lyα-forest auto-correlation function, yielding DH(z = 2.40) /rd = 8.94 ± 0.22 and DM(z = 2.40) /rd = 36.6 ± 1.2, within 2.3σ of the same flat-ΛCDM model.

The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to billions of galaxies. We present Galactos, a high-performance implementation of a novel, O(N2) algorithm that uses a load-balanced k-d tree and spherical harmonic expansions to compute the anisotropic 3PCF. Our implementation is optimized for the Intel Xeon Phi architecture, exploiting SIMD parallelism, instruction and thread concurrency, and significant L1 and L2 cache reuse, reaching 39% of peak performance on a single node. Galactos scales to the full Cori system, achieving 9.8 PF (peak) and 5.06 PF (sustained) across 9636 nodes, making the 3PCF easily computable for all galaxies in the observable universe.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg2 and volume of 18.7 Gpc3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51 and 0.61. We measure the angular diameter distance DM and Hubble parameterHfrom the baryon acoustic oscillation (BAO) method, in combinationwith a cosmic microwave background prior on the sound horizon scale, after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DMH from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by fσ8(z), from redshift-space distortions (RSD). We combine individual measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method; in particular, the AP measurement from sub-BAO scales sharpens constraints from post-reconstruction BAOs by breaking degeneracy between DM and H. Combined with Planck 2016 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature ωK = 0.0003 ± 0.0026 and a dark energy equation-of-state parameter ω =-1.01 ± 0.06, in strong affirmation of the spatially flat cold dark matter (CDM) model with a cosmological constant (ΛCDM). Our RSD measurements of fσ8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3 ± 1.0 kms-1 Mpc-1 even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8 ± 1.2 kms-1 Mpc-1. Assuming flat ΛCDM, we find ωm = 0.310 ± 0.005 and H0 = 67.6 ± 0.5 kms-1 Mpc-1, and we find a 95 per cent upper limit of 0.16 eV c-2 on the neutrino mass sum.

We present a configuration-space model of the large-scale galaxy 3-point correlation function (3PCF) based on leading-order perturbation theory and including redshift-space distortions (RSD). This model should be useful in extracting distance-scale information from the 3PCF via the baryon acoustic oscillation method. We include the first redshift-space treatment of biasing by the baryon-dark matter relative velocity. Overall, on large scales the effect of RSD is primarily a renormalization of the 3PCF that is roughly independent of both physical scale and triangle opening angle; for our adoptedm and bias values, the rescaling is a factor of ∼1.8. We also present an efficient scheme for computing 3PCF predictions from our model, important for allowing fast exploration of the space of cosmological parameters in future analyses.

We present a measurement of the linear growth rate of structure, f, from the Sloan Digital Sky Survey III (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) using convolution Lagrangian perturbation theory (CLPT) with Gaussian streaming redshift space distortions (GSRSD) to model the two-point statistics of BOSS galaxies in DR12. The BOSS-DR12 data set includes 1198 006 massive galaxies spread over the redshift range 0.2 < z < 0.75. These galaxy samples are categorized in three redshift bins. Using CLPT-GSRSD in our analysis of the combined sample of the three redshift bins, we report measurements of fσ8 for the three redshift bins. We find fσ8 = 0.430 ± 0.054 at zeff = 0.38, fσ8 = 0.452 ± 0.057 at zeff = 0.51 and fσ8 = 0.457 ± 0.052 at zeff = 0.61. Our results are consistent with the predictions of Planck Λ cold dark matter-general relativity. Our constraints on the growth rates of structure in the Universe at different redshifts serve as a useful probe, which can help distinguish between a model of the Universe based on dark energy and models based on modified theories of gravity. This paper is part of a set that analyses the final galaxy clustering data set from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al., to produce the final cosmological constraints from BOSS.

We present the large-scale three-point correlation function (3PCF) of the Sloan Digital Sky Survey DR12 Constant stellar Mass (CMASS) sample of 777 202 Luminous Red Galaxies, the largest-ever sample used for a 3PCF or bispectrum measurement. We make the first high-significance (4.5σ) detection of baryon acoustic oscillations (BAO) in the 3PCF. Using these acoustic features in the 3PCF as a standard ruler, we measure the distance to z = 0.57 to 1.7 per cent precision (statistical plus systematic). We find DV =2024±29 Mpc (stat)±20 Mpc (sys) for our fiducial cosmology (consistent with Planck 2015) and bias model. This measurement extends the use of the BAO technique from the twopoint correlation function (2PCF) and power spectrum to the 3PCF and opens an avenue for deriving additional cosmological distance information from future large-scale structure redshift surveys such as DESI. Our measured distance scale from the 3PCF is fairly independent from that derived from the pre-reconstruction 2PCF and is equivalent to increasing the length of BOSS by roughly 10 per cent; reconstruction appears to lower the independence of the distance measurements. Fitting a model including tidal tensor bias yields a moderatesignificance (2.6σ) detection of this bias with a value in agreement with the prediction from local Lagrangian biasing.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratios in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially resolved spectroscopy for thousands of nearby galaxies (median z ∼ 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between z ~ 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGNs and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5 m Sloan Foundation Telescope at the Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5 m du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in 2016 July.

We have used flux-transmission correlations in Lyα forests to measure the imprint of baryon acoustic oscillations (BAO). The study uses spectra of 157 783 quasars in the redshift range 2:1 ≤ z ≤ 3:5 from the Sloan Digital Sky Survey (SDSS) data release 12 (DR12). Besides the statistical improvements on our previous studies using SDSS DR9 and DR11, we have implemented numerous improvements in the analysis procedure, allowing us to construct a physical model of the correlation function and to investigate potential systematic errors in the determination of the BAO peak position. The Hubble distance, DH = c=H(z), relative to the sound horizon is DH(z = 2:33)=rd = 9:07 ± 0:31. The best-determined combination of comoving angular-diameter distance, DM, and the Hubble distance is found to be D0:7 H D0:3 M =rd = 13:94±0:35. This value is 1:028±0:026 times the prediction of the flat-ΛCDM model consistent with the cosmic microwave background (CMB) anisotropy spectrum. The errors include marginalization over the effects of unidentified high-density absorption systems and fluctuations in ultraviolet ionizing radiation. Independently of the CMB measurements, the combination of our results and other BAO observations determine the open-ΛCDM density parameters to be ΩM = 0:296 ± 0:029, ΩΛ = 0:699 ± 0:100 and Ωk = -0:002 ± 0:119.

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

Spectroscopy with the James Webb Space Telescope has opened the possibility of identifying moderate-luminosity active galactic nuclei (AGNs) in the early Universe, at and beyond the epoch of re-ionisation, complementing previous surveys of much more luminous (and much rarer) quasars. We present 12 new AGNs at 4 < z < 7 in the JADES survey (in addition to the previously identified AGN in GN-z11 at z = 10.6) revealed through the detection of a broad-line region (BLR) seen in the Balmer emission lines. The depth of JADES, together with the use of three different spectral resolutions, enables us to probe a lower-mass regime relative to previous studies. In a few cases, we find evidence for two broad components of Hα, which suggests that these could be candidate merging black holes (BHs), although a complex BLR geometry cannot be excluded. The inferred BH masses range from 8 × 107 M⊙ down to 4 × 105 M⊙, interestingly probing the regime expected for direct collapse BHs. The inferred AGN bolometric luminosities (∼1044 − 1045 erg/s) imply accretion rates that are < 0.5 times the Eddington rate in most cases. However, small BHs, with MBH ∼ 106 M⊙, tend to accrete at Eddington or super-Eddington rates. These BHs at z ∼ 4–11 are over-massive relative to their host galaxies' stellar masses when compared to the local MBH − Mstar relation, even approaching MBH ∼ Mstar, as was expected from heavy BH seeds and/or super-Eddington accretion scenarios. However, we find that these early BHs tend to be more consistent with the local relation between MBH and velocity dispersion, as well as between MBH and dynamical mass, suggesting that these are more fundamental and universal relations. On the classical, optical narrow-line excitation-diagnostic diagrams, these AGNs are located in the region that is locally occupied by star-forming galaxies, implying that they would be missed by the standard classification techniques if they did not display broad lines. Their location on the diagram is consistent with what is expected for AGNs hosted in metal-poor galaxies (Z ∼ 0.1 − 0.2 Z⊙). The fraction of broad-line AGNs with LAGN > 1044 erg/s among galaxies in the redshift range of 4 < z < 6 is about 10%, suggesting that the contribution of AGNs and their hosts to the re-ionisation of the Universe is > 10%. [ABSTRACT FROM AUTHOR]

We analyse the baryon acoustic oscillation (BAO) signal of the final Baryon Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is performed in the Fourier space, using the power spectrum monopole and quadrupole. The data set includes 1198 006 galaxies over the redshift range 0.2 < z < 0.75.We divide this data set into three (overlapping) redshift bins with the effective redshifts zeff = 0.38, 0.51 and 0.61. We demonstrate the reliability of our analysis pipeline using N-body simulations as well as ~1000 MultiDark- Patchy mock catalogues that mimic the BOSS-DR12 target selection. We apply density field reconstruction to enhance the BAO signal-to-noise ratio. By including the power spectrum quadrupole we can separate the line of sight and angular modes, which allows us to constrain the angular diameter distance DA(z) and the Hubble parameter H(z) separately. We obtain two independent 1.6 and 1.5 per cent constraints on DA(z) and 2.9 and 2.3 per cent constraints on H(z) for the low (zeff = 0.38) and high (zeff = 0.61) redshift bin, respectively. We obtain two independent 1 and 0.9 per cent constraints on the angular averaged distance DV(z), when ignoring the Alcock-Paczynski effect. The detection significance of the BAO signal is of the order of 8σ (post-reconstruction) for each of the three redshift bins. Our results are in good agreement with the Planck prediction within Λ cold dark matter. This paper is part of a set that analyses the final galaxy clustering data set from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. to produce the final cosmological constraints from BOSS.

We test for galactic conformity at 0.2 < z < 1.0 to a projected distance of 5Mpc using spectroscopic redshifts from the PRism MUlti-object Survey (PRIMUS). Our sample consists of ∼60,000 galaxies in five separate fields covering a total of ∼5.5 square degrees, which allows us to account for cosmic variance. We identify star-forming and quiescent "isolated primary" (i.e., central) galaxies using isolation criteria and cuts in specific star formation rate. We match the redshift and stellar mass distributions of these samples to control for correlations between quiescent fraction and redshift and stellar mass. We detect a significant (σ) one-halo conformity signal, or an excess of star-forming neighbors around star-forming central galaxies, of ∼5% on scales of 0-1 Mpc and a 2.5σ two-halo signal of ∼1% on scales of 1-3 Mpc. These signals are weaker than those detected in the Sloan Digital Sky Survey and are consistent with galactic conformity being the result of large-scale tidal fields and reflecting assembly bias. We also measure the star-forming fraction of central galaxies at fixed stellar mass as a function of large-scale environment and find that central galaxies are more likely to be quenched in overdense environments, independent of stellar mass. However, we find that environment does not affect the star formation efficiency of central galaxies, as long as they are forming stars. We test for redshift and stellar mass dependence of the conformity signal within our sample and show that large volumes and multiple fields are required at intermediate redshift to adequately account for cosmic variance.