Search

We present the measurements of the small-scale clustering for the emission line galaxy (ELG) sample from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) in the Sloan Digital Sky Survey IV (SDSS-IV). We use conditional abundance matching method to interpret the clustering measurements from $0.34h^{-1}\textrm{Mpc}$ to $70h^{-1}\textrm{Mpc}$. In order to account for the correlation between properties of emission line galaxies and their environment, we add a secondary connection between star formation rate of ELGs and halo accretion rate. Three parameters are introduced to model the ELG [OII] luminosity and to mimic the target selection of eBOSS ELGs. The parameters in our models are optimized using Markov Chain Monte Carlo (MCMC) method. We find that by conditionally matching star formation rate of galaxies and the halo accretion rate, we are able to reproduce the eBOSS ELG small scale clustering within 1$\sigma$ error level. Our best fit model shows that the eBOSS ELG sample only consists of $\sim 12\%$ of all star-forming galaxies, and the satellite fraction of eBOSS ELG sample is 19.3\%. We show that the effect of assembly bias is $\sim20\%$ on the two-point correlation function and $\sim5\%$ on the void probability function at scale of $r\sim 20 h^{-1}\rm Mpc$., Comment: MNRAS, Volume 519, Issue 3

Over the next five years, the Dark Energy Spectroscopic Instrument (DESI) will use 10 spectrographs with 5000 fibers on the 4m Mayall Telescope at Kitt Peak National Observatory to conduct the first Stage-IV dark energy galaxy survey. At $z < 0.6$, the DESI Bright Galaxy Survey (BGS) will produce the most detailed map of the Universe during the dark energy dominated epoch with redshifts of >10 million galaxies over 14,000 deg$^2$. In this work, we present and validate the final BGS target selection and survey design. From the Legacy Surveys, BGS will target a $r < 19.5$ magnitude-limited sample (BGS Bright); a fainter $19.5 < r < 20.175$ sample, color-selected to have high redshift efficiency (BGS Faint); and a smaller low-z quasar sample. BGS will observe these targets using exposure times, scaled to achieve uniform completeness, and visit each point on the footprint three times. We use observations from the Survey Validation programs conducted prior to the main survey along with realistic simulations to show that BGS can complete its strategy and make optimal use of `bright' time. We demonstrate that BGS targets have stellar contamination <1% and that their densities do not depend strongly on imaging properties. We also confirm that BGS Bright will achieve >80% fiber assignment efficiency. Finally, we show that BGS Bright and Faint will achieve >95% redshift success rates with no significant dependence on observing conditions. BGS meets the requirements for an extensive range of scientific applications. BGS will yield the most precise Baryon Acoustic Oscillations and Redshift-Space Distortions measurements at $z < 0.4$. It also presents opportunities to exploit new methods that require highly complete and dense galaxy samples (e.g. N-point statistics, multi-tracers). BGS further provides a powerful tool to study galaxy populations and the relations between galaxies and dark matter., Comment: AJ, submitted, 34 pages, 22 figures, one of a suite of 8 papers detailing targeting for DESI

We analyze spatially resolved and co-added SDSS-IV MaNGA spectra with signal-to-noise ~100 from 2200 passive central galaxies (z~0.05) to understand how central galaxy assembly depends on stellar mass (M*) and halo mass (Mh). We control for systematic errors in Mh by employing a new group catalog from Tinker (2020a,b) and the widely-used Yang et al. (2007) catalog. At fixed M*, the strength of several stellar absorption features varies systematically with Mh. Completely model-free, this is one of the first indications that the stellar populations of centrals with identical M* are affected by the properties of their host halos. To interpret these variations, we applied full spectral fitting with the code alf. At fixed M*, centrals in more massive halos are older, show lower [Fe/H], and have higher [Mg/Fe] with 3.5 sigma confidence. We conclude that halos not only dictate how much M* galaxies assemble, but also modulate their chemical enrichment histories. Turning to our analysis at fixed Mh, high-M* centrals are older, show lower [Fe/H], and have higher [Mg/Fe] for Mh>10^{12}Msun/h with confidence > 4 sigma. While massive passive galaxies are thought to form early and rapidly, our results are among the first to distinguish these trends at fixed Mh. They suggest that high-M* centrals experienced unique early formation histories, either through enhanced collapse and gas fueling, or because their halos were early-forming and highly concentrated, a possible signal of galaxy assembly bias., Comment: Accepted for publication in ApJ. 28 pages and 12 figures

We analyze clustering measurements of BOSS galaxies using a simulation-based emulator of two-point statistics. We focus on the monopole and quadrupole of the redshift-space correlation function, and the projected correlation function, at scales of $0.1\sim60~h^{-1}$Mpc. Although our simulations are based on $w$CDM with general relativity (GR), we include a scaling parameter of the halo velocity field, $\gamma_f$, defined as the amplitude of the halo velocity field relative to the GR prediction. We divide the BOSS data into three redshift bins. After marginalizing over other cosmological parameters, galaxy bias parameters, and the velocity scaling parameter, we find $f\sigma_{8}(z=0.25) = 0.413\pm0.031$, $f\sigma_{8}(z=0.4) = 0.470\pm0.026$ and $f\sigma_{8}(z=0.55) = 0.396\pm0.022$. Compared with Planck observations using a flat $\Lambda$CDM model, our results are lower by $1.9\sigma$, $0.3\sigma$ and $3.4\sigma$ respectively. These results are consistent with other recent simulation-based results at non-linear scales, including weak lensing measurements of BOSS LOWZ galaxies, two-point clustering of eBOSS LRGs, and an independent clustering analysis of BOSS LOWZ. All these results are generally consistent with a combination of $\gamma_f^{1/2}\sigma_8\approx 0.75$. We note, however, that the BOSS data is well fit assuming GR, i.e. $\gamma_f=1$. We cannot rule out an unknown systematic error in the galaxy bias model at non-linear scales, but near-future data and modeling will enhance our understanding of the galaxy--halo connection, and provide a strong test of new physics beyond the standard model., Comment: 30 pages, 23 figures and 4 tables; ApJ accepted, updated to the match the accepted version

The PRObabilistic Value-Added Bright Galaxy Survey (PROVABGS) catalog will provide measurements of galaxy properties, such as stellar mass ($M_*$), star formation rate (${\rm SFR}$), stellar metallicity ($Z_{\rm MW}$), and stellar age ($t_{\rm age, MW}$), for >10 million galaxies of the DESI Bright Galaxy Survey. Full posterior distributions of the galaxy properties will be inferred using state-of-the-art Bayesian spectral energy distribution (SED) modeling of DESI spectroscopy and Legacy Surveys photometry. In this work, we present the SED model, Bayesian inference framework, and methodology of PROVABGS. Furthermore, we apply the PROVABGS SED modeling on realistic synthetic DESI spectra and photometry, constructed using the L-GALAXIES semi-analytic model. We compare the inferred galaxy properties to the true galaxy properties of the simulation using a hierarchical Bayesian framework to quantify accuracy and precision. Overall, we accurately infer the true $M_*$, ${\rm SFR}$, $Z_{\rm MW}$, and $t_{\rm age, MW}$ of the simulated galaxies. However, the priors on galaxy properties induced by the SED model have a significant impact on the posteriors. They impose a ${\rm SFR}{>}10^{-1} M_\odot/{\rm yr}$ lower bound on ${\rm SFR}$, a ${\sim}0.3$ dex bias on $\log Z_{\rm MW}$ for galaxies with low spectral signal-to-noise, and $t_{\rm age, MW} < 8\,{\rm Gyr}$ upper bound on stellar age. This work also demonstrates that a joint analysis of spectra and photometry significantly improves the constraints on galaxy properties over photometry alone and is necessary to mitigate the impact of the priors. With the methodology presented and validated in this work, PROVABGS will maximize information extracted from DESI observations and provide a probabilistic value-added galaxy catalog that will extend current galaxy studies to new regimes and unlock cutting-edge probabilistic analyses., Comment: 40 pages, 17 figures, submitted to ApJ, the PROVABGS SED modeling pipeline is publicly available at https://github.com/changhoonhahn/provabgs

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 present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitude relations, $(g-r) - M_r$ and $(FUV-NUV)-M_r$, of the simulations to a $M_r < -20$ and UV complete SDSS galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed color-magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuation-slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low $A_V$ since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations., Comment: 26 pages, 15 figures

We apply a new galaxy group finder to the Main Galaxy Sample of the SDSS. This algorithm introduces new freedom to assign halos to galaxies that is self-calibrated by comparing the catalog to complementary data. These include galaxy clustering data and measurements of the total satellite luminosity from deep imaging data. We present constraints on the galaxy-halo connection for star-forming and quiescent populations. The results of the self-calibrated group catalog differ in several key ways from previous group catalogs and halo occupation analyses. The transition halo mass scale, where half of halos contain quiescent central galaxies, is at M_h~10^12.4 Msol/h, significantly higher than other constraints. Additionally, the width of the transition from predominantly star-forming halos to quiescent halos occurs over a narrower range in halo mass. Quiescent central galaxies in low-mass halos are significantly more massive than star-forming centrals at the same halo mass, but this difference reverses above the transition halo mass. We find that the scatter in logM* at fixed M_h is ~0.2 dex for massive halos, in agreement with previous estimates, but rises sharply at lower halo masses. The halo masses assigned by the group catalog are in good agreement with weak lensing estimates for star-forming and quiescent central galaxies. We discuss possible improvements to the algorithm made clear by this first application to data. The group catalog is made publicly available., Comment: 27 pages, 23 figures, submitted to ApJ. Group catalog available at https://www.galaxygroupfinder.net/

We use the technique of total satellite luminosity, L_sat, to probe the dark matter halos around active galactic nuclei in the SDSS Main Galaxy Sample. Our results focus on galaxies and AGN that are the central galaxy of their halo. Our two AGN samples are constructed from optical emission-line diagnostics and from WISE infrared colors. Both optically-selected and WISE-selected AGN have L_sat values twice as high as non-active galaxy samples when controlling for stellar mass and mean stellar age. This implies that the halos are twice as massive, but we cannot rule out that the increase in L_sat is due to these AGN residing in younger halos at the same mass. When only controlling for host galaxy stellar mass, WISE-selected AGN also have higher L_sat values than optical AGN at the factor of two level, consistent with previous results comparing the clustering of obscured and unobscured AGN. However, controlling for stellar age in the two populations of host galaxies removes half of this difference, attenuating the statistical significance of the difference. We perform permutation tests to quantify the different in the halo populations of each sample. The difference in star formation properties does not fully explain the difference in the two AGN populations, however. Although AGN luminosity correlates with mean stellar age, the difference in stellar age between the WISE and optical samples cannot be fully explained by differences in their AGN luminosity distributions., Comment: 10 pages, 7 figures, submitted to MNRAS

We describe an extension of the halo-based galaxy group-finding algorithm. We add freedom to the algorithm in order to more accurately determine which galaxies are central and which are satellites, and to provide unbiased estimates of halo masses. We focus on determination of the galaxy-halo relations for star-forming and quiescent galaxies. The added freedom in the group-finding algorithm is self-calibrated using observations of color-dependent galaxy clustering, as well as measurements of the total satellite luminosity in deep imaging data around stacked samples of spectroscopic central galaxies, L_sat. We test this approach on a series of mocks that vary the galaxy-halo connection, including one mock constructed from UniverseMachine results. Our self-calibrated algorithm shows marked improvement over previous methods in estimating the color-dependent satellite fraction of galaxies. It reduces the error in log M_halo for central galaxies by over a factor of two, to <~0.2 dex. Through the L_sat data, it can quantify differences in the luminosity-to-halo mass relations for star-forming and quiescent galaxies, even for groups with only one spectroscopic member. Whereas previous algorithms cannot constrain the scatter in L_gal at fixed M_halo, the self-calibration technique can provide a robust lower limit to this scatter., Comment: 20 pages, 13 figures, submitted to ApJ

The completed eBOSS catalogues contain redshifts of 344080 QSOs over 0.8

We present 2000 mock galaxy catalogs for the analysis of baryon acoustic oscillations in the Emission Line Galaxy (ELG) sample of the Extended Baryon Oscillation Spectroscopic Survey Data Release 16 (eBOSS DR16). Each mock catalog has a number density of $6.7 \times 10^{-4} h^3 \rm Mpc^{-3}$, covering a redshift range from 0.6 to 1.1. The mocks are calibrated to small-scale eBOSS ELG clustering measurements at scales of around 10 $h^{-1}$Mpc. The mock catalogs are generated using a combination of GaLAxy Mocks (GLAM) simulations and the Quick Particle-Mesh (QPM) method. GLAM simulations are used to generate the density field, which is then assigned dark matter halos using the QPM method. Halos are populated with galaxies using a halo occupation distribution (HOD). The resulting mocks match the survey geometry and selection function of the data, and have slightly higher number density which allows room for systematic analysis. The large-scale clustering of mocks at the baryon acoustic oscillation (BAO) scale is consistent with data and we present the correlation matrix of the mocks., Comment: 12 pages, 11 figures, submitted to MNRAS

We produce 1000 realizations of synthetic clustering catalogues for each type of the tracers used for the baryon acoustic oscillation and redshift space distortion analysis of the Sloan Digital Sky Surveys-IV extended Baryon Oscillation Spectroscopic Survey final data release (eBOSS DR16), covering the redshift range from 0.6 to 2.2, to provide reliable estimates of covariance matrices and test the robustness of the analysis pipeline with respect to observational systematics. By extending the Zel'dovich approximation density field with an effective tracer bias model calibrated with the clustering measurements from the observational data, we accurately reproduce the two- and three-point clustering statistics of the eBOSS DR16 tracers, including their cross-correlations in redshift space with very low computational costs. In addition, we include the gravitational evolution of structures and sample selection biases at different redshifts, as well as various photometric and spectroscopic systematic effects. The agreements on the auto-clustering statistics between the data and mocks are generally within 1 $\sigma$ variances inferred from the mocks, for scales down to a few $h^{-1}\,{\rm Mpc}$ in configuration space, and up to $0.3\,h\,{\rm Mpc}^{-1}$ in Fourier space. For the cross correlations between different tracers, the same level of consistency presents in configuration space, while there are only discrepancies in Fourier space for scales above $0.15\,h\,{\rm Mpc}^{-1}$. The accurate reproduction of the data clustering statistics permits reliable covariances for multi-tracer analysis., Comment: 25 pages, 22 figures; updated to match the version accepted for publication at MNRAS

We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Ly$\alpha$ forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter distances and Hubble distances relative to the sound horizon, $r_d$, from eight different samples and six measurements of the growth rate parameter, $f\sigma_8$, from redshift-space distortions (RSD). This composite sample is the most constraining of its kind and allows us to perform a comprehensive assessment of the cosmological model after two decades of dedicated spectroscopic observation. We show that the BAO data alone are able to rule out dark-energy-free models at more than eight standard deviations in an extension to the flat, $\Lambda$CDM model that allows for curvature. When combined with Planck Cosmic Microwave Background (CMB) measurements of temperature and polarization the BAO data provide nearly an order of magnitude improvement on curvature constraints. The RSD measurements indicate a growth rate that is consistent with predictions from Planck primary data and with General Relativity. When combining the results of SDSS BAO and RSD with external data, all multiple-parameter extensions remain consistent with a $\Lambda$CDM model. Regardless of cosmological model, the precision on $\Omega_\Lambda$, $H_0$, and $\sigma_8$, remains at roughly 1\%, showing changes of less than 0.6\% in the central values between models. The inverse distance ladder measurement under a o$w_0w_a$CDM yields $H_0= 68.20 \pm 0.81 \, \rm km\, s^{-1} Mpc^{-1}$, remaining in tension with several direct determination methods. (abridged), Comment: A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://www.sdss.org/science/final-bao-and-rsd-measurements/ . The full cosmological interpretation of these measurements can be found here: https://www.sdss.org/science/cosmology-results-from-eboss/

We present large-scale structure catalogs from the completed extended Baryon Oscillation Spectroscopic Survey (eBOSS). Derived from Sloan Digital Sky Survey (SDSS) -IV Data Release 16 (DR16), these catalogs provide the data samples, corrected for observational systematics, and random positions sampling the survey selection function. Combined, they allow large-scale clustering measurements suitable for testing cosmological models. We describe the methods used to create these catalogs for the eBOSS DR16 Luminous Red Galaxy (LRG) and Quasar samples. The quasar catalog contains 343,708 redshifts with $0.8 < z < 2.2$ over 4,808\,deg$^2$. We combine 174,816 eBOSS LRG redshifts over 4,242\,deg$^2$ in the redshift interval $0.6 < z < 1.0$ with SDSS-III BOSS LRGs in the same redshift range to produce a combined sample of 377,458 galaxy redshifts distributed over 9,493\,deg$^2$. Improved algorithms for estimating redshifts allow that 98 per cent of LRG observations result in a successful redshift, with less than one per cent catastrophic failures ($\Delta z > 1000$ ${\rm km~s}^{-1}$). For quasars, these rates are 95 and 2 per cent (with $\Delta z > 3000$ ${\rm km~s}^{-1}$). We apply corrections for trends between the number densities of our samples and the properties of the imaging and spectroscopic data. For example, the quasar catalog obtains a $\chi^2$/DoF$= 776/10$ for a null test against imaging depth before corrections and a $\chi^2$/DoF$=6/8$ after. The catalogs, combined with careful consideration of the details of their construction found here-in, allow companion papers to present cosmological results with negligible impact from observational systematic uncertainties., Comment: Matches version accepted by MNRAS, very minor changes. A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://www.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/

We analyse the clustering of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 16 luminous red galaxy sample (DR16 eBOSS LRG) in combination with the high redshift tail of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey Data Release 12 (DR12 BOSS CMASS). We measure the redshift space distortions (RSD) and also extract the longitudinal and transverse baryonic acoustic oscillation (BAO) scale from the anisotropic power spectrum signal inferred from 377,458 galaxies between redshifts 0.6 and 1.0, with effective redshift of $z_{\rm eff}=0.698$ and effective comoving volume of $2.72\,{\rm Gpc}^3$. After applying reconstruction we measure the BAO scale and infer $D_H(z_{\rm eff})/r_{\rm drag} = 19.30\pm 0.56$ and $D_M(z_{\rm eff})/r_{\rm drag} =17.86 \pm 0.37$. When we perform a redshift space distortions analysis on the pre-reconstructed catalogue on the monopole, quadrupole and hexadecapole we find, $D_H(z_{\rm eff})/r_{\rm drag} = 20.18\pm 0.78$, $D_M(z_{\rm eff})/r_{\rm drag} =17.49 \pm 0.52$ and $f\sigma_8(z_{\rm eff})=0.454\pm0.046$. We combine both sets of results along with the measurements in configuration space of \cite{LRG_corr} and report the following consensus values: $D_H(z_{\rm eff})/r_{\rm drag} = 19.77\pm 0.47$, $D_M(z_{\rm eff})/r_{\rm drag} = 17.65\pm 0.30$ and $f\sigma_8(z_{\rm eff})=0.473\pm 0.044$, which are in full agreement with the standard $\Lambda$CDM and GR predictions. These results represent the most precise measurements within the redshift range $0.6\leq z \leq 1.0$ and are the culmination of more than 8 years of SDSS observations., Comment: 41 pages, 20 figures; Minor updates to match the published version in MNRAS. A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://www.sdss.org/science/final-bao-and-rsd-measurements/ . The full cosmological interpretation of these measurements can be found here: https://www.sdss.org/science/cosmology-results-from-eboss/

We present the Emission Line Galaxy (ELG) sample of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) from the Sloan Digital Sky Survey IV Data Release 16 (DR16). After describing the observations and redshift measurement for the 269,243 observed ELG spectra over 1170 deg$^2$, we present the large-scale structure catalogues, which are used for the cosmological analysis. These catalogues contain 173,736 reliable spectroscopic redshifts between 0.6 and 1.1, along with the associated random catalogues quantifying the extent of observations, and the appropriate weights to correct for non-cosmological fluctuations. We perform a spherically averaged baryon acoustic oscillations (BAO) measurement in configuration space, with density field reconstruction: the data 2-point correlation function shows a feature consistent with that of the BAO, providing a 3.2-percent measurement of the spherically averaged BAO distance $D_V(z_{\rm eff})/r_{\rm drag} = 18.23\pm 0.58$ at the effective redshift $z_{\rm eff}=0.845$., Comment: Submitted to MNRAS. A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://www.sdss.org/science/final-bao-and-rsd-measurements/. The full cosmological interpretation of these measurements can be found here: https://www.sdss.org/science/cosmology-results-from-eboss/

We demonstrate how the total luminosity in satellite galaxies is a powerful probe of dark matter halos around central galaxies. The method cross-correlates central galaxies in spectroscopic galaxy samples with fainter galaxies detected in photometric surveys. After background subtraction, the excess galaxies around the central galaxies represent faint satellite galaxies within the dark matter halo. Using abundance matching models, we show that the the total galaxy luminosity, L_sat, scales linearly with host halo mass, making L_sat an excellent proxy for M_h. L_sat is also sensitive to the formation time of the halo, as younger halos have more substructure at fixed M_h. We demonstrate that probes of galaxy large-scale environment can break this degeneracy. Although this is an indirect probe of the halo, it can yield a high-S/N measurement for galaxies expected to occupy halos at $<10^{12}$ M_sol, where other methods suffer from larger errors. In this paper we focus on observational and theoretical systematics in the L_sat method. We test the robustness of our method of finding central galaxies and our methods of estimating the number of background galaxies. We implement this method on central galaxies in SDSS data, with satellites identified in faint imaging from the DESI Legacy Imaging Surveys. We find excellent agreement between our theoretical predictions and the observational measurements. Finally, we compare our L_sat measurements to weak lensing estimates of M_h for red and blue subsamples. In the stellar mass range where the measurements overlap, we find consistent results, where red galaxies live in larger halos. However, the L_sat approach allows us to probe significantly lower mass galaxies. At these masses, the L_sat values are equivalent. This example shows the potential of L_sat as a probe of dark halos. (Abridged), Comment: 21 pages, submitted to MNRAS

The total luminosity of satellite galaxies around a central galaxy, L$_{sat}$, is a powerful metric for probing dark matter halos. In this paper we use data from the Sloan Digital Sky Survey and DESI Legacy Imaging Surveys to explore the relationship between L$_{sat}$ and various observable galaxy properties for a sample of 117,966 central galaxies out to $z = 0.15$. At fixed stellar mass, every galaxy property we explore shows a correlation with L$_{sat}$. This implies that dark matter halos play a possibly significant role in determining these secondary galaxy properties. We quantify these correlations by computing the mutual information between L$_{sat}$ and secondary properties and explore how this mutual information varies as a function of stellar mass and when separating the sample into star-forming and quiescent central galaxies. We find that absolute r-band magnitude correlates more strongly with L$_{sat}$ than stellar mass across all galaxy populations; and that effective radius, velocity dispersion, and S\'ersic index do so as well for star-forming and quiescent galaxies. L$_{sat}$ is sensitive to both the mass of the host halo as well as the halo formation history, with younger halos having higher L$_{sat}$. L$_{sat}$ by itself cannot distinguish between these two effects, but measurements of galaxy large-scale environment can break this degeneracy. For star-forming central galaxies, we find that r$_{\rm eff}$, $\sigma_v$, and S\'ersic index all correlate with large-scale density, implying that these halo age plays a role in determining these properties. For quiescent galaxies, we find that all secondary properties are independent of environment, implying that correlations with L$_{sat}$ are driven only by halo mass. These results are a significant step forward in quantifying the full extent of the galaxy-halo connection, and present a new test of galaxy formation models., Comment: 21 pages, 16 figures. To be submitted to MNRAS

We develop and implement two new methods for constraining the scatter in the relationship between galaxies and dark matter halos. These new techniques are sensitive to the scatter at low halo masses, making them complementary to previous constraints that are dependent on clustering amplitudes or rich galaxy groups, both of which are only sensitive to more massive halos. In both of our methods, we use a galaxy group finder to locate central galaxies in the SDSS main galaxy sample. Our first technique uses the small-scale cross-correlation of central galaxies with all lower mass galaxies. This quantity is sensitive to the satellite fraction of low-mass galaxies, which is in turn driven by the scatter between halos and galaxies. The second technique uses the kurtosis of the distribution of line-of-sight velocities between central galaxies and neighboring galaxies. This quantity is sensitive to the distribution of halo masses that contain the central galaxies at fixed stellar mass. Theoretical models are constructed using peak halo circular velocity, $V_{\rm peak}$, as our property to connect galaxies to halos. The cross-correlation technique yields a constraint of $\sigma[ M_\ast|V_{\rm peak}]=0.27\pm 0.05$ dex, corresponding to a scatter in $\log M_\ast$ at fixed $M_h$ of $\sigma[ M_\ast|M_h]=0.38\pm 0.06$ dex at $M_h=10^{11.8}$ Msun. The kurtosis technique yields $\sigma[ M_\ast|V_{\rm peak}]=0.30\pm0.03$, corresponding to $\sigma[ M_\ast|M_h]=0.34\pm 0.04$ at $M_h=10^{12.2}$ Msun. The values of $\sigma[ M_\ast|M_h]$ are significantly larger than the constraints at higher masses, in agreement with the results of hydrodynamic simulations. This increase is only partly due to the scatter between $V_{\rm peak}$ and $M_h$, and it represents an increase of nearly a factor of two relative to the values inferred from clustering and group studies at high masses., Comment: 14 pages, 13 figures, submitted to MNRAS

We present constraints on the timescale of star formation variability and the correlation between star formation and host halo accretion histories of star-forming (SF) central galaxies from the scatter of the stellar-to-halo mass relation (SHMR). SF galaxies are found to have a tight relationship between their star formation rates and stellar masses on the so-called "star-forming sequence" (SFS), which characterizes their SF histories and $M_*$ growths. Meanwhile, observed SHMR constraints connect $M_*$ growth to halo accretion history. Combining these observed trends with a cosmological $N$-body simulation, we present flexible models that track the SF, $M_*$, and halo accretion histories of SF central galaxies at $z<1$ while reproducing the observed stellar mass function and SFS in SDSS DR7. We find that the scatter in SHMR at $M_h=10^{12}M_\odot$, $\sigma_{M_*|M_h=10^{12}M_\odot}$, is sensitive to the timescale of star formation variability, $t_{\rm duty}$, and the correlation coefficient, $r$, between SF and halo accretion histories: shorter $t_{\rm duty}$ and higher $r$ both result in tighter $\sigma_{M_*|M_h=10^{12}M_\odot}$. To reproduce a constant 0.2 dex scatter over $z=1-0$, our models require $t_{\rm duty}\leq1.5$ Gyr for $r=0.99$ or $r>0.6$ for $t_{\rm duty}=0.1$ Gyr. For $r\sim0.6$, as found in the literature, $t_{\rm duty}<0.2$ Gyr is necessary. Meanwhile, to reproduce the tightening of $\sigma_{M_*|M_h=10^{12}M_\odot}=0.35$ to 0.2 dex from $z=1-0$ in hydrodynamical simulations, our models require $t_{\rm duty}=0.1$ Gyr for $r>0.5$. Although, the lack of consensus on $\sigma_{M_*|M_h=10^{12}M_\odot}$ at $M_h=10^{12}M_\odot$ and at $z=1$ from observations and galaxy formation models remains the main bottleneck in precisely constraining $r$ and $t_{\rm duty}$, we demonstrate that SHMR can be used to constrain the SF and host halo accretion histories of SF central galaxies., Comment: 27 pages, 9 figures

Models of the spatial distribution of dark matter halos must achieve new levels of precision and accuracy in order to satisfy the requirements of upcoming experiments. In this work, we present a halo bias emulator for modeling the clustering of halos on large scales. It incorporates the cosmological dependence of the bias beyond the mapping of halo mass to peak height. The emulator makes substantial improvements in accuracy compared to the widely used Tinker et al. (2010) model. Halos in this work are defined using an overdensity criteria of 200 relative to the mean background density. Halo catalogs are produced for 40 N-body simulations as part of the Aemulus project at snapshots from z=3 to z=0. The emulator is trained over the mass range $6\times10^{12}-7\times10^{15}\ h^{-1}M_{\odot}$. Using an additional suite of 35 simulations, we determine that the precision of the emulator is redshift dependent, achieving sub-percent levels for a majority of the redshift range. Two additional simulation suites are used to test the ability of the emulator to extrapolate to higher and lower masses. Our high-resolution simulation suite is used to develop an extrapolation scheme in which the emulator asymptotes to the Tinker et al. (2010) model at low mass, achieving ~3% accuracy down to $10^{11}\ h^{-1}M_{\odot}$. Finally, we present a method to propagate emulator modeling uncertainty into an error budget. Our emulator is made publicly available at \url{https://github.com/AemulusProject/bias_emulator}., Comment: 14 pages, 7 figures, Submitted to ApJ

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models., Comment: 28 pages, 15 figures

Using the N-body simulations of the AEMULUS Project, we construct an emulator for the non-linear clustering of galaxies in real and redshift space. We construct our model of galaxy bias using the halo occupation framework, accounting for possible velocity bias. The model includes 15 parameters, including both cosmological and galaxy bias parameters. We demonstrate that our emulator achieves ~ 1% precision at the scales of interest, 0.1

The rapidly growing statistical precision of galaxy surveys has lead to a need for ever-more precise predictions of the observables used to constrain cosmological and galaxy formation models. The primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper we present one such suite of simulations, forming the basis for the AEMULUS Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h^-1 Gpc)^3 simulations with mass resolution and force softening of 3.51\times 10^10 (Omega_m / 0.3) ~ h^-1 M_sun and 20 ~ h^-1 kpc respectively in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) for halos with more than 500 (200) particles respectively and scales of r>200 ~ h^-1 kpc in real space or k ~ 3 h Mpc^-1 in harmonic space for z\le 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k ~ 4 h^-1 Mpc. We provide the halo catalogs and snapshots detailed in this work to the community at https://AemulusProject.github.io., Comment: 16 pages, 12 figures, 3 Tables Project website: https://aemulusproject.github.io/

Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a halo mass function emulator and demonstrate improved accuracy relative to state-of-the-art analytic models. In this work, mass is defined using an overdensity criteria of 200 relative to the mean background density. Our emulator is constructed from the AEMULUS simulations, a suite of 40 N-body simulations with snapshots from z=3 to z=0. These simulations cover the flat wCDM parameter space allowed by recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results, varying the parameters w, Omega_m, Omega_b, sigma_8, N_{eff}, n_s, and H_0. We validate our emulator using five realizations of seven different cosmologies, for a total of 35 test simulations. These test simulations were not used in constructing the emulator, and were run with fully independent initial conditions. We use our test simulations to characterize the modeling uncertainty of the emulator, and introduce a novel way of marginalizing over the associated systematic uncertainty. We confirm non-universality in our halo mass function emulator as a function of both cosmological parameters and redshift. Our emulator achieves better than 1% precision over much of the relevant parameter space, and we demonstrate that the systematic uncertainty in our emulator will remain a negligible source of error for cluster abundance studies through at least the LSST Year 1 data set., Comment: https://aemulusproject.github.io/

In our modern understanding of galaxy formation, every galaxy forms within a dark matter halo. The formation and growth of galaxies over time is connected to the growth of the halos in which they form. The advent of large galaxy surveys as well as high-resolution cosmological simulations has provided a new window into the statistical relationship between galaxies and halos and its evolution. Here we define this galaxy-halo connection as the multi-variate distribution of galaxy and halo properties that can be derived from observations and simulations. This connection provides a key test of physical galaxy formation models; it also plays an essential role in constraints of cosmological models using galaxy surveys and in elucidating the properties of dark matter using galaxies. We review techniques for inferring the galaxy-halo connection and the insights that have arisen from these approaches. Some things we have learned are that galaxy formation efficiency is a strong function of halo mass; at its peak in halos around a pivot halo mass of 10^12 Msun, less than 20% of the available baryons have turned into stars by the present day; the intrinsic scatter in galaxy stellar mass is small, less than 0.2 dex at a given halo mass above this pivot mass; below this pivot mass galaxy stellar mass is a strong function of halo mass; the majority of stars over cosmic time were formed in a narrow region around this pivot mass. We also highlight key open questions about how galaxies and halos are connected, including understanding the correlations with secondary properties and the connection of these properties to galaxy clustering., Comment: 56 pages, 13 figures. Invited review, submitted to Annual Reviews of Astronomy & Astrophysics, Vol 56

We present the clustering measurements of quasars in configuration space based on the Data Release 14 (DR14) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. This dataset includes 148,659 quasars spread over the redshift range $0.8\leq z \leq 2.2$ and spanning 2112.9 square degrees. We use the Convolution Lagrangian Perturbation Theory (CLPT) approach with a Gaussian Streaming (GS) model for the redshift space distortions of the correlation function and demonstrate its applicability for dark matter halos hosting eBOSS quasar tracers. At the effective redshift $z_{\rm eff} = 1.52$, we measure the linear growth rate of structure $f\sigma_{8}(z_{\rm eff})= 0.426 \pm 0.077$, the expansion rate $H(z_{\rm eff})= 159^{+12}_{-13}(r_{s}^{\rm fid}/r_s){\rm km.s}^{-1}.{\rm Mpc}^{-1}$, and the angular diameter distance $D_{A}(z_{\rm eff})=1850^{+90}_{-115}\,(r_s/r_{s}^{\rm fid}){\rm Mpc}$, where $r_{s}$ is the sound horizon at the end of the baryon drag epoch and $r_{s}^{\rm fid}$ is its value in the fiducial cosmology. The quoted errors include both systematic and statistical contributions. The results on the evolution of distances are consistent with the predictions of flat $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmology with Planck parameters, and the measurement of $f\sigma_{8}$ extends the validity of General Relativity (GR) to higher redshifts($z>1$) This paper is released with companion papers using the same sample. The results on the cosmological parameters of the studies are found to be in very good agreement, providing clear evidence of the complementarity and of the robustness of the first full-shape clustering measurements with the eBOSS DR14 quasar sample., Comment: 25 pages, 22 figures

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the EZ mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of $0.8

We analyse the clustering of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample (DR14Q). We measure the redshift space distortions using the power spectrum monopole, quadrupole and hexadecapole inferred from 148,659 quasars between redshifts 0.8 and 2.2 covering a total sky footprint of 2112.9 deg$^2$. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f\sigma_8$, and the Alcock-Paczynski dilation scales which allow constraints to be placed on the angular diameter distance $D_A(z)$ and the Hubble $H(z)$ parameter. At the effective redshift of $z_{\rm eff}=1.52$, $f\sigma_8(z_{\rm eff})=0.420\pm0.076$, $H(z_{\rm eff})=[162\pm 12]\, (r_s^{\rm fid}/r_s)\,{\rm km\, s}^{-1}{\rm Mpc}^{-1}$, and $D_A(z_{\rm eff})=[1.85\pm 0.11]\times10^3\,(r_s/r_s^{\rm fid})\,{\rm Mpc}$, where $r_s$ is the comoving sound horizon at the baryon drag epoch and the superscript `fid' stands for its fiducial value. The errors take into account the full error budget, including systematics and statistical contributions. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$CDM) cosmological model inferred from Planck measurements. Finally, we compare our measurements with other eBOSS companion papers and find excellent agreement, demonstrating the consistency and complementarity of the different methods used for analysing the data., Comment: 37 pages, 25 figures, 14 tables. Data and covariance will be made public after the acceptance of the paper. Submitted to MNRAS

We present a public data release of halo catalogs from a suite of 125 cosmological $N$-body simulations from the Abacus project. The simulations span 40 $w$CDM cosmologies centered on the Planck 2015 cosmology at two mass resolutions, $4\times 10^{10}\;h^{-1}M_\odot$ and $1\times 10^{10}\;h^{-1}M_\odot$, in $1.1\;h^{-1}\mathrm{Gpc}$ and $720\;h^{-1}\mathrm{Mpc}$ boxes, respectively. The boxes are phase-matched to suppress sample variance and isolate cosmology dependence. Additional volume is available via 16 boxes of fixed cosmology and varied phase; a few boxes of single-parameter excursions from Planck 2015 are also provided. Catalogs spanning $z=1.5$ to $0.1$ are available for friends-of-friends and Rockstar halo finders and include particle subsamples. All data products are available at https://lgarrison.github.io/AbacusCosmos, Comment: 13 pages, 9 figures, 3 tables. Additional figures added for mass resolution convergence tests, and additional redshifts added for existing tests. Matches ApJS accepted version

We measure how the properties of star-forming central galaxies correlate with large-scale environment, $\delta$, measured on $10$Mpc/h scales. We use group catalogs to isolate a robust sample of central galaxies with high purity and completeness. The properties we investigate are star formation rate (SFR), exponential disk scale length $R_{\rm exp}$, and Sersic index of the light profile, $n$. We find that, at all stellar masses, there is an inverse correlation between SFR and $\delta$, meaning that above-average star forming centrals live in underdense regions. For $n$ and $R_{\rm exp}$, there is no correlation with $\delta$ at $M_{\rm star}\lesssim 10^{10.5}$ $M_\odot$, but at higher masses there are positive correlations; a weak correlation with $R_{\rm exp}$ and a strong correlation with $n$. These data are evidence of assembly bias within the star-forming population. The results for SFR are consistent with a model in which SFR correlates with present-day halo accretion rate, $\dot{M}_h$. In this model, galaxies are assigned to halos using the abundance matching ansatz, which maps galaxy stellar mass onto halo mass. At fixed halo mass, SFR is assigned to galaxies using the same approach,but $\dot{M}_h$ is used to map onto SFR. The best-fit model requires some scatter in the $\dot{M}_h$-SFR relation. The $R_{\rm exp}$ and $n$ measurements are consistent with a model in which these quantities are correlated with the spin parameter of the halo, $\lambda$. Halo spin does not correlate with $\delta$ at low halo masses, but for higher mass halos, high-spin halos live in higher density environments at fixed $M_h$. Put together with the earlier installments of this series, these data demonstrate that quenching processes have limited correlation with halo formation history, but the growth of active galaxies, as well as other detailed properties, are influenced by the details of halo assembly., Comment: 13 pages, submitted to MNRAS

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 dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $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 $\Lambda$CDM best-fit 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 $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature., Comment: Accepted by MNRAS; BAO distance likelihood available in source files 'QSOv1.9fEZmock_BAOchi2.dat'; full set of data to be public eventually from SDSS website

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 ratio 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 redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts 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 AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter 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 July 2016., Comment: Published in Astronomical Journal

Using group catalogs from the SDSS DR7, we attempt to measure galactic conformity in the local universe. We measure the quenched fraction of neighbor galaxies around isolated primary galaxies, dividing the isolated sample into star-forming and quiescent objects. We restrict our measurements to scales $>1$ Mpc to probe the correlations between the formation histories of distinct halos. Over the stellar mass range $10^{9.7} \le M_\ast/M_\odot \le 10^{10.9}$, we find minimal statistical evidence for conformity. We further compare these data to predictions of the halo age-matching model, in which the oldest galaxies are associated with the oldest halos at fixed $M_\ast$. For models with strong correlations between halo and stellar age, the conformity signal is too large to be consistent with the data. For weaker implementations of age-matching, galactic conformity is not a sensitive diagnostic of halo assembly bias, and would not produce a detectable signal in SDSS data. We reproduce the results of Kauffmann et al 2013, in which the star formation rates of neighbor galaxies are significantly reduced around primary galaxies when the primaries are themselves low star formers. However, we find this result is mainly driven by contamination in the isolation criterion, when using our group catalog to remove the small fraction of satellite galaxies in the sample, the conformity signal largely goes away. Lastly, we show that small conformity signals, i.e., 2-5% differences in the quenched fractions of neighbor galaxies, can be produced by mechanisms other than halo assembly bias. For example, if passive galaxies occupy more massive halos than star forming galaxies of the same stellar mass, a conformity signal that is consistent with recent measurements from PRIMUS (Berti et al 2016) can be produced., Comment: 13 pages, 7 figures, submitted to MNRAS

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey., Comment: 27 pages, 3 figures and one table. A supplementary document is included. The BOSS DR12 BAO data used in the work can be downloaded from the SDSS website https://sdss3.org//science/boss_publications.php (check https://data.sdss.org/sas/dr12/boss/papers/clustering/Zhao_etal_2016_DR12COMBINED_tomoBAO_powspec.tar.gz for the tomographic BAO measurement presented in https://arxiv.org/abs/1607.03153)

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

We investigate the potential sources of theoretical systematics in the anisotropic Baryon Acoustic Oscillation (BAO) distance scale measurements from the clustering of galaxies in configuration space using the final Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS). We perform a detailed study of the impact on BAO measurements from choices in the methodology such as fiducial cosmology, clustering estimators, random catalogues, fitting templates, and covariance matrices. The theoretical systematic uncertainties in BAO parameters are found to be 0.002 in the isotropic dilation $\alpha$ and 0.003 in the quadrupolar dilation $\epsilon$. The leading source of systematic uncertainty is related to the reconstruction techniques. Theoretical uncertainties are sub-dominant compared with the statistical uncertainties for BOSS survey, accounting $0.2\sigma_{stat}$ for $\alpha$ and $0.25\sigma_{stat}$ for $\epsilon$ ($\sigma_{\alpha,stat} \sim$0.010 and $\sigma_{\epsilon,stat}\sim$ 0.012 respectively). We also present BAO-only distance scale constraints from the anisotropic analysis of the correlation function. Our constraints on the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$, including both statistical and theoretical systematic uncertainties, are 1.5\% and 2.8\% at $z_{\rm eff}=0.38$, 1.4\% and 2.4\% at $z_{\rm eff}=0.51$, and 1.7\% and 2.6\% at $z_{\rm eff}=0.61$. This paper is part of a set that analyzes the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are cross-checked with other BAO analysis in \citet{Acacia16}. The systematic error budget concerning the methodology on post-reconstruction BAO analysis presented here is used in \citet{Acacia16} to produce the final cosmological constraints from BOSS., Comment: 38 pages, 22 figures, accepted for publication in MNRAS

Central galaxies make up the majority of the galaxy population, including the majority of the quiescent population at $\mathcal{M}_* > 10^{10}\mathrm{M}_\odot$. Thus, the mechanism(s) responsible for quenching central galaxies plays a crucial role in galaxy evolution as whole. We combine a high resolution cosmological $N$-body simulation with observed evolutionary trends of the "star formation main sequence," quiescent fraction, and stellar mass function at $z < 1$ to construct a model that statistically tracks the star formation histories and quenching of central galaxies. Comparing this model to the distribution of central galaxy star formation rates in a group catalog of the SDSS Data Release 7, we constrain the timescales over which physical processes cease star formation in central galaxies. Over the stellar mass range $10^{9.5}$ to $10^{11} \mathrm{M}_\odot$ we infer quenching e-folding times that span $1.5$ to $0.5\; \mathrm{Gyr}$ with more massive central galaxies quenching faster. For $\mathcal{M}_* = 10^{10.5}\mathrm{M}_\odot$, this implies a total migration time of $\sim 4~\mathrm{Gyrs}$ from the star formation main sequence to quiescence. Compared to satellites, central galaxies take $\sim 2~\mathrm{Gyrs}$ longer to quench their star formation, suggesting that different mechanisms are responsible for quenching centrals versus satellites. Finally, the central galaxy quenching timescale we infer provides key constraints for proposed star formation quenching mechanisms. Our timescale is generally consistent with gas depletion timescales predicted by quenching through strangulation. However, the exact physical mechanism(s) responsible for this still remain unclear., Comment: 16 pages, 11 figures

Fiber-fed multi-object spectroscopic surveys, with their ability to collect an unprecedented number of redshifts, currently dominate large-scale structure studies. However, physical constraints limit these surveys from successfully collecting redshifts from galaxies too close to each other on the focal plane. This ultimately leads to significant systematic effects on galaxy clustering measurements. Using simulated mock catalogs, we demonstrate that fiber collisions have a significant impact on the power spectrum, $P(k)$, monopole and quadrupole that exceeds sample variance at scales smaller than $k\sim0.1~h/Mpc$. We present two methods to account for fiber collisions in the power spectrum. The first, statistically reconstructs the clustering of fiber collided galaxy pairs by modeling the distribution of the line-of-sight displacements between them. It also properly accounts for fiber collisions in the shot-noise correction term of the $P(k)$ estimator. Using this method, we recover the true $P(k)$ monopole of the mock catalogs with residuals of $<0.5\%$ at $k=0.3~h/Mpc$ and $<4\%$ at $k=0.83~h/Mpc$ -- a significant improvement over existing correction methods. The quadrupole, however, does not improve significantly. The second method models the effect of fiber collisions on the power spectrum as a convolution with a configuration space top-hat function that depends on the physical scale of fiber collisions. It directly computes theoretical predictions of the fiber-collided $P(k)$ multipoles and reduces the influence of smaller scales to a set of nuisance parameters. Using this method, we reliably model the effect of fiber collisions on the monopole and quadrupole down to the scale limits of theoretical predictions. The methods we present in this paper will allow us to robustly analyze galaxy power spectrum multipole measurements to much smaller scales than previously possible., Comment: 18 pages, 10 figures