Source author record

David H. Weinberg

David H. Weinberg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

149works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

149 published item(s)

preprint2026arXiv

Fisher Forecasts for Cosmological Yields from $3\!\times\!2$pt Analysis of the Roman Space Telescope High Latitude Imaging Survey

The High Latitude Imaging Survey (HLIS) of NASA's Nancy Grace Roman Space Telescope will provide powerful tests of cosmological models through sensitive measurements of cosmic shear, galaxy-galaxy lensing (GGL), and galaxy clustering. As part of the HLIS Project Infrastructure Team's Data Challenge 1 (DC1), we carry out Fisher forecasts of cosmological parameter constraints from combinations of these probes, focusing on inverse-variance figures of merit (FoMs) for the parameters $σ_8$ and $Ω_{\rm{m}}$, which scale the amplitude of weak lensing signals. We find good agreement between Fisher analysis and Markov chain Monte Carlo (MCMC) analysis of the DC1 baseline data vector, and we exploit the flexibility of Fisher analysis to investigate varied priors on cosmological parameters and on nuisance parameters describing unknown biases in photometric redshifts or shear measurements. Given the benchmark DC1 priors, the forecast constraints from GGL+clustering are substantially stronger than those from cosmic shear, with the combination of all three probes (``$3\!\times\!2$pt'') providing moderate further improvement. Adding tight external priors on the power spectrum shape parameters $n_{\rm{s}}$, $Ω_{\rm{b}}$, and $h_0$ can improve the $(σ_8, Ω_{\rm{m}})$ FoMs by factors of $1.2$--$3.5$. The smallest scale angular bins provide much more information than the largest scale bins, and the highest redshift tomographic bins provide more information than the lowest redshift bins. Factor-of-two changes in the priors on photo-$z$ and shear biases, relative to the benchmark values based on anticipated calibration accuracy, produce changes of $\lesssim 20\%$ in FoMs, implying robust cosmological performance if this demanding level of accuracy can be achieved.

preprint2022arXiv

A New Model For Including Galactic Winds in Simulations of Galaxy Formation I: Introducing the Physically Evolved Winds (PhEW) Model

The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is over-simplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds (PhEW), that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a sub-grid prescription to model galactic winds more robustly both physically and numerically.

preprint2022arXiv

A New Model For Including Galactic Winds in Simulations of Galaxy Formation II: Implementation of PhEW in Cosmological Simulations

Although galactic winds play a critical role in regulating galaxy formation, hydrodynamic cosmological simulations do not resolve the scales that govern the interaction between winds and the ambient circumgalactic medium (CGM). We implement the Physically Evolved Wind (PhEW) model of Huang et al. (2020) in the GIZMO hydrodynamics code and perform test cosmological simulations with different choices of model parameters and numerical resolution. PhEW adopts an explicit subgrid model that treats each wind particle as a collection of clouds that exchange mass, metals, and momentum with their surroundings and evaporate by conduction and hydrodynamic instabilities as calibrated on much higher resolution cloud scale simulations. In contrast to a conventional wind algorithm, we find that PhEW results are robust to numerical resolution and implementation details because the small scale interactions are defined by the model itself. Compared to conventional wind simulations with the same resolution, our PhEW simulations produce similar galaxy stellar mass functions at $z\geq 1$ but are in better agreement with low-redshift observations at $M_* < 10^{11}M_\odot$ because PhEW particles shed mass to the CGM before escaping low mass halos. PhEW radically alters the CGM metal distribution because PhEW particles disperse metals to the ambient medium as their clouds dissipate, producing a CGM metallicity distribution that is skewed but unimodal and is similar between cold and hot gas. While the temperature distributions and radial profiles of gaseous halos are similar in simulations with PhEW and conventional winds, these changes in metal distribution will affect their predicted UV/X-ray properties in absorption and emission.

preprint2022arXiv

Elucidating Galaxy Assembly Bias in SDSS

We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy sample using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as $\mathcal{Q}_\mathrm{cen}$ and $\mathcal{Q}_\mathrm{sat}$, respectively, to predict a suite of one- and two-point observables. In particular, our fiducial constraint employs the probability distribution of the galaxy number counts measured on $8\,\mathrm{Mpc}\,h^{-1}$ scales $N_8^g$ and the projected cross-correlation functions of quintiles of galaxies selected by $N_8^g$ with our entire galaxy sample. We perform extensive tests of the efficacy of our method by fitting the same observables to mock data using both constrained and non-constrained simulations. We discover that in many cases the level of cosmic variance between the two simulations can produce biased constraints that lead to an erroneous detection of galaxy assembly bias if the non-constrained simulation is used. When applying our method to the SDSS data, the ELUCID reconstruction effectively removes an otherwise strong degeneracy between cosmic variance and galaxy assembly bias in SDSS, enabling us to derive an accurate and stringent constraint on the latter. Our fiducial ELUCID constraint, for galaxies above a stellar mass threshold $M_*{=}10^{10.2}\,h^{-2}\,M_\odot$, is $\mathcal{Q}_\mathrm{cen}{=}{-}0.09\pm{0.05}$ and $\mathcal{Q}_\mathrm{sat}{=}0.09\pm{0.10}$, indicating no evidence for a significant~($>2σ$) galaxy assembly bias in the local Universe probed by SDSS. Finally, our method provides a promising path to the robust modelling of the galaxy-halo connection within future surveys like DESI and PFS.

preprint2022arXiv

How Many Elements Matter?

Some studies of stars' multi-element abundance distributions suggest at least 5-7 significant dimensions, but others show that many elemental abundances can be predicted to high accuracy from [Fe/H] and [Mg/Fe] (or [Fe/H] and age) alone. We show that both propositions can be, and are, simultaneously true. We adopt a machine learning technique known as normalizing flow to reconstruct the probability distribution of Milky Way disk stars in the space of 15 elemental abundances measured by APOGEE. Conditioning on Teff and log g minimizes the differential systematics. After further conditioning on [Fe/H] and [Mg/Fe], the residual scatter for most abundances is $σ_{[X/{\rm H}]} \lesssim 0.02$ dex, consistent with APOGEE's reported statistical uncertainties of $\sim$0.01-0.015 dex and intrinsic scatter of 0.01-0.02 dex. Despite the small scatter, residual abundances display clear correlations between elements, which we show are too large to be explained by measurement uncertainties or by the finite sampling noise. We must condition on at least seven elements to reduce correlations to a level consistent with observational uncertainties. Our results demonstrate that cross-element correlations are a much more sensitive probe of hidden structure than dispersion, and they can be measured precisely in a large sample even if star-by-star measurement noise is comparable to the intrinsic scatter. We conclude that many elements have an independent story to tell, even for the "mundane" disk stars and elements produced by core-collapse and Type Ia supernovae. The only way to learn these lessons is to measure the abundances directly, and not merely infer them.

preprint2022arXiv

Optical selection bias and projection effects in stacked galaxy cluster weak lensing

Cosmological constraints from current and upcoming galaxy cluster surveys are limited by the accuracy of cluster mass calibration. In particular, optically identified galaxy clusters are prone to selection effects that can bias the weak lensing mass calibration. We investigate the selection bias of the stacked cluster lensing signal associated with optically selected clusters, using clusters identified by the redMaPPer algorithm in the Buzzard simulations as a case study. We find that at a given cluster halo mass, the residuals of redMaPPer richness and weak lensing signal are positively correlated. As a result, for a given richness selection, the stacked lensing signal is biased high compared with what we would expect from the underlying halo mass probability distribution. The cluster lensing selection bias can thus lead to overestimated mean cluster mass and biased cosmology results. We show that the lensing selection bias exhibits a strong scale-dependence and is approximately 20 to 60 percent for $ΔΣ$ at large scales. This selection bias largely originates from spurious member galaxies within +/- 20 to 60 Mpc/h along the line of sight, highlighting the importance of quantifying projection effects associated with the broad redshift distribution of member galaxies in photometric cluster surveys. While our results qualitatively agree with those in the literature, accurate quantitative modelling of the selection bias is needed to achieve the goals of cluster lensing cosmology and will require synthetic catalogues covering a wide range of galaxy-halo connection models.

preprint2022arXiv

Primordial helium-3 redux: The helium isotope ratio of the Orion nebula

We report the first direct measurement of the helium isotope ratio, 3He/4He, outside of the Local Interstellar Cloud, as part of science verification observations with the upgraded CRyogenic InfraRed Echelle Spectrograph (CRIRES). Our determination of 3He/4He is based on metastable HeI* absorption along the line-of-sight towards Tet02 Ori A in the Orion Nebula. We measure a value 3He/4He=(1.77+/-0.13)x10^{-4}, which is just ~40 per cent above the primordial relative abundance of these isotopes, assuming the Standard Model of particle physics and cosmology, (3He/4He)_p = (1.257+/-0.017)x10^-4. We calculate a suite of galactic chemical evolution simulations to study the Galactic build up of these isotopes, using the yields from Limongi & Chieffi (2018) for stars in the mass range M=8-100 M_sun and Lagarde (2011,2012) for M=0.8-8 M_sun. We find that these simulations simultaneously reproduce the Orion and protosolar 3He/4He values if the calculations are initialized with a primordial ratio (3He/4He)_p=(1.043+/-0.089)x10^-4. Even though the quoted error does not include the model uncertainty, this determination agrees with the Standard Model value to within ~2sigma. We also use the present-day Galactic abundance of deuterium (D/H), helium (He/H), and 3He/4He to infer an empirical limit on the primordial 3He abundance, (3He/H)_p < (1.09+/-0.18)x10^-5, which also agrees with the Standard Model value. We point out that it is becoming increasingly difficult to explain the discrepant primordial 7Li/H abundance with non-standard physics, without breaking the remarkable simultaneous agreement of three primordial element ratios (D/H, 4He/H, and 3He/4He) with the Standard Model values.

preprint2021arXiv

APOGEE Chemical Abundance Patterns of the Massive Milky Way Satellites

The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey has obtained high-resolution spectra for thousands of red giant stars distributed among the massive satellite galaxies of the Milky Way (MW): the Large and Small Magellanic Clouds (LMC/SMC), the Sagittarius Dwarf (Sgr), Fornax (Fnx), and the now fully disrupted \emph{Gaia} Sausage/Enceladus (GSE) system. We present and analyze the APOGEE chemical abundance patterns of each galaxy to draw robust conclusions about their star formation histories, by quantifying the relative abundance trends of multiple elements (C, N, O, Mg, Al, Si, Ca, Fe, Ni, and Ce), as well as by fitting chemical evolution models to the [$α$/Fe]-[Fe/H] abundance plane for each galaxy. Results show that the chemical signatures of the starburst in the MCs observed by Nidever et al. in the $α$-element abundances extend to C+N, Al, and Ni, with the major burst in the SMC occurring some 3-4 Gyr before the burst in the LMC. We find that Sgr and Fnx also exhibit chemical abundance patterns suggestive of secondary star formation epochs, but these events were weaker and earlier ($\sim$~5-7 Gyr ago) than those observed in the MCs. There is no chemical evidence of a second starburst in GSE, but this galaxy shows the strongest initial star formation as compared to the other four galaxies. All dwarf galaxies had greater relative contributions of AGB stars to their enrichment than the MW. Comparing and contrasting these chemical patterns highlight the importance of galaxy environment on its chemical evolution.

preprint2021arXiv

Chemical Cartography with APOGEE: Mapping Disk Populations with a Two-Process Model and Residual Abundances

We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] vs. [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core collapse supernovae and a "delayed" component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals Δ[X/H] from this two-parameter fit. The rms residuals range from ~0.01-0.03 dex for the most precisely measured APOGEE abundances to ~0.1 dex for Na, V, and Ce. The correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow comparison of populations controlled for differences in metallicity and [α/Fe]. Relative to the main disk (R=3-13 kpc, |Z|<2 kpc), we find nearly identical abundance patterns in the outer disk (R=15-17 kpc), 0.05-0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4-1 dex) of multiple elements in ωCen. Residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements.

preprint2021arXiv

Exploiting Non-linear Scales in Galaxy-Galaxy Lensing and Galaxy Clustering: A Forecast for the Dark Energy Survey

The combination of galaxy-galaxy lensing (GGL) and galaxy clustering is a powerful probe of low redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we extend the N-body and halo occupation distribution (HOD) emulator method of arxiv:1907.06293 to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding parameters that describe central galaxy incompleteness, galaxy assembly bias, and a scale-independent multiplicative lensing bias $A_{lens}$. We use this emulator to forecast cosmological constraints attainable from the GGL surface density profile $ΔΣ(r_p)$ and the projected galaxy correlation function $w_{p,gg}(r_p)$ in the final (Year 6) DES data set over scales $r_p=0.3-30h^{-1}$ Mpc. For a $3\%$ prior on $A_{lens}$ we forecast precisions of $1.9\%$, $2.0\%$, and $1.9\%$ on $Ω_m$, $σ_8$, and $S_8 \equiv σ_8Ω_m^{0.5}$, marginalized over all halo occupation distribution (HOD) parameters as well as $A_{lens}$ and a point-mass contribution to $ΔΣ$. Adding scales $r_p=0.3-3h^{-1}$ Mpc improves the $S_8$ precision by a factor of $\sim1.6$ relative to a large scale ($3.0-30.0h^{-1}$ Mpc) analysis, equivalent to increasing the survey area by a factor of ${\sim}2.6$. Sharpening the $A_{lens}$ prior to $1\%$ further improves the $S_8$ precision by a factor of $1.7$ (to $1.1\%$), and it amplifies the gain from including non-linear scales. Our emulator achieves percent-level accuracy similar to the projected DES statistical uncertainties, demonstrating the feasibility of a fully non-linear analysis. Obtaining precise parameter constraints from multiple galaxy types and from measurements that span linear and non-linear clustering offers many opportunities for internal cross-checks, which can diagnose systematics and demonstrate the robustness of cosmological results.

preprint2021arXiv

Residual Abundances in GALAH DR3: Implications for Nucleosynthesis and Identification of Unique Stellar Populations

We investigate the [X/Mg] abundances of 16 elements for 82,910 Galactic disk stars from GALAH+ DR3. We fit the median trends of low-Ia and high-Ia populations with a two-process model, which describes stellar abundances in terms of a prompt core-collapse and delayed Type-Ia supernova component. For each sample star, we fit the amplitudes of these two components and compute the residual $Δ$[X/H] abundances from this two-parameter fit. We find RMS residuals $\lesssim 0.07$ dex for well-measured elements and correlated residuals among some elements (such as Ba, Y, and Zn) that indicate common enrichment sources. From a detailed investigation of stars with large residuals, we infer that roughly $40\%$ of the large deviations are physical and $60\%$ are caused by problematic data such as unflagged binarity, poor wavelength solutions, and poor telluric subtraction. As one example of a population with distinctive abundance patterns, we identify 15 stars that have 0.3-0.6 dex enhancements of Na but normal abundances of other elements from O to Ni and positive average residuals of Cu, Zn, Y, and Ba. We measure the median elemental residuals of 14 open clusters, finding systematic $\sim0.1-0.4$ dex enhancements of O, Ca, K, Y, and Ba and $\sim0.2$ dex depletion of Cu in young clusters. Finally, we present a restricted three-process model where we add an asymptotic giant branch star (AGB) component to better fit Ba and Y. With the addition of the third process, we identify a population of stars, preferentially young, that have much higher AGB enrichment than expected from their SNIa enrichment.

preprint2021arXiv

The High Latitude Spectroscopic Survey on the Nancy Grace Roman Space Telescope

The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0-1.93 $μ$m, $R=435-865$) and the 0.28 deg$^2$ wide field camera. We present a reference HLSS which maps 2000 deg$^2$ and achieves an emission line flux limit of 10$^{-16}$ erg/s/cm$^2$ at 6.5$σ$, requiring $\sim$0.6 yrs of observing time. We summarize the flowdown of the Roman science objectives to the science and technical requirements of the HLSS. We construct a mock redshift survey over the full HLSS volume by applying a semi-analytic galaxy formation model to a cosmological N-body simulation, and use this mock survey to create pixel-level simulations of 4 deg$^2$ of HLSS grism spectroscopy. We find that the reference HLSS would measure $\sim$ 10 million H$α$ galaxy redshifts that densely map large scale structure at $z=1-2$ and 2 million [OIII] galaxy redshifts that sparsely map structures at $z=2-3$. We forecast the performance of this survey for measurements of the cosmic expansion history with baryon acoustic oscillations and the growth of large scale structure with redshift space distortions. We also study possible deviations from the reference design, and find that a deep HLSS at $f_{\rm line}>7\times10^{-17}$erg/s/cm$^2$ over 4000 deg$^2$ (requiring $\sim$1.5 yrs of observing time) provides the most compelling stand-alone constraints on dark energy from Roman alone. This provides a useful reference for future optimizations. The reference survey, simulated data sets, and forecasts presented here will inform community decisions on the final scope and design of the Roman HLSS.

preprint2020arXiv

Inside Out and Upside-Down: The Roles of Gas Cooling and Dynamical Heating in Shaping the Stellar Age-Velocity Relation

Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relationship (AVR) of h277 is nearly identical to that of the analogous solar-neighborhood measurement in the Milky Way. A crucial element of this success is the simulation's dynamically cold multi-phase ISM, which allows young stars to form with a low velocity dispersion ($σ_{\mathrm{birth}}$$\sim 6 - 8\ \mathrm{km}\, \mathrm{s}^{-1}$) at late times. Older stars are born kinematically hotter (i.e., the disk settles over time in an "upside-down" formation scenario), and are subsequently heated after birth. The disk also grows "inside-out", and many of the older stars in the solar neighborhood at $z=0$ are present because of radial mixing. We demonstrate that the evolution of $σ_{\mathrm{birth}}$ in h277 can be explained by the same model used to describe the general decrease in velocity dispersion observed in disk galaxies from $z\sim 2-3$ to the present-day, in which the disk evolves in quasi-stable equilibrium and the ISM velocity dispersion decreases over time due to a decreasing gas fraction. Thus, our results tie together local observations of the Milky Way's AVR with observed kinematics of high $z$ disk galaxies.

preprint2020arXiv

The Impact of Wind Scalings on Stellar Growth and the Baryon Cycle in Cosmological Simulations

Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed $v_w$ and the scaling of the mass loading factor $η$ with velocity dispersion $σ$. We design our fiducial model to reproduce the measured wind properties at 25% of the virial radius from the Feedback In Realistic Environments (FIRE) simulations. We find that a strong dependence of $η\sim σ^5$ in low mass haloes with $σ< 106\ \mathrm{km\ s^{-1}}$ is required to match the faint end of the stellar mass functions at $z > 1$. The wind speed also has a major impact, with faster winds significantly reducing wind recycling and heating more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z=2$ and $z=0$ within a single model, suggesting that an additional feedback source such as AGN might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $\mathrm{Mpc/h}$, $2\times576^3$ simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $z\sim0-2$. In closing, the results from simulations of galaxy formation are much more sensitive to small changes in the feedback implementation than to the hydrodynamic technique.

preprint2020arXiv

The Robustness of Cosmological Hydrodynamic Simulation Predictions to Changes in Numerics and Cooling Physics

We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$α$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

preprint2019arXiv

Abundance ratios in GALAH DR2 and their implications for nucleosynthesis

Using a sample of 70 924 stars from the second data release of the GALAH optical spectroscopic survey, we construct median sequences of [X/Mg] vs. [Mg/H] for 21 elements, separating the high-$α$/``low-Ia'' and low-$α$/``high-Ia'' stellar populations through cuts in [Mg/Fe]. Previous work with the near-IR APOGEE survey has shown that such sequences are nearly independent of location in the Galactic disk, implying that they are determined by stellar nucleosynthesis yields with little sensitivity to other chemical evolution aspects. The separation between the two [X/Mg] sequences indicates the relative importance of prompt and delayed enrichment mechanisms, while the sequences' slopes indicate metallicity dependence of the yields. GALAH and APOGEE measurements agree for some of their common elements, but differ in sequence separation or metallicity trends for others. GALAH offers access to nine new elements. We infer that about $75\%$ of solar C comes from core collapse supernovae and $25\%$ from delayed mechanisms. We find core collapse fractions of $60-80\%$ for the Fe-peak elements Sc, Ti, Cu, and Zn, with strong metallicity dependence of the core collapse Cu yield. For the neutron capture elements Y, Ba, and La, we infer large delayed contributions with non-monotonic metallicity dependence. The separation of the [Eu/Mg] sequences implies that at least $\sim30\%$ of Eu enrichment is delayed with respect to star formation. We compare our results to predictions of several supernova and AGB yield models; C, Na, K, Mn, and Ca all show discrepancies with models that could make them useful diagnostics of nucleosynthesis physics.

preprint2019arXiv

Cosmology with galaxy-galaxy lensing on non-perturbative scales: Emulation method and application to BOSS LOWZ

We describe our nonlinear emulation (i.e., interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with $N$-body simulations of nonlinear structure formation, designed to accurately predict the projected clustering and galaxy-galaxy lensing signals from luminous red galaxies (LRGs) in the redshift range $0.16 < z < 0.36$ on comoving scales $0.6 < r_p < 30$ \hMpc. The interpolation accuracy is $\lesssim 1-2$ per cent across the entire physically plausible range of parameters for all scales considered. We correctly recover the true value of the cosmological parameter $S_8 = ({σ_8}/{0.8228}) ({Ω_{\text{m}}}/{0.3107})^{0.6}$ from mock measurements produced via subhalo abundance matching (SHAM)-based lightcones designed to approximately match the properties of the SDSS LOWZ galaxy sample. Applying our model to Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 14 (DR14) LOWZ galaxy clustering and galaxy-shear cross-correlation measurements made with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) imaging, we perform a prototype cosmological analysis marginalizing over $w$CDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of $S_8 = 0.847 \pm 0.037$, in $3.5σ$ tension with the Planck cosmological results of $1.00 \pm 0.02$. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy.

preprint2019arXiv

Probing black hole accretion tracks, scaling relations and radiative efficiencies from stacked X-ray active galactic nuclei

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

preprint2016arXiv

Galaxy Three-Point Correlation Functions and Halo/Subhalo Models

We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.

preprint2016arXiv

Modeling Lyman-α Forest Cross-Correlations with LyMAS

We use the Ly-$α$ Mass Association Scheme (LyMAS; Peirani et al. 2014) to predict cross-correlations at $z=2.5$ between dark matter halos and transmitted flux in the Ly-$α$ forest, and compare to cross-correlations measured for quasars and damped Ly-$α$ systems (DLAs) from the Baryon Oscillation Spectroscopic Survey (BOSS) by Font-Ribera et al. (2012, 2013). We calibrate LyMAS using Horizon-AGN hydrodynamical cosmological simulations of a $(100\ h^{-1}\ \mathrm{Mpc})^3$ comoving volume. We apply this calibration to a $(1\ h^{-1}\ \mathrm{Gpc})^3$ simulation realized with $2048^3$ dark matter particles. In the 100 $h^{-1}$ Mpc box, LyMAS reproduces the halo-flux correlations computed from the full hydrodynamic gas distribution very well. In the 1 $h^{-1}$ Gpc box, the amplitude of the large scale cross-correlation tracks the halo bias $b_h$ as expected. We provide empirical fitting functions that describe our numerical results. In the transverse separation bins used for the BOSS analyses, LyMAS cross-correlation predictions follow linear theory accurately down to small scales. Fitting the BOSS measurements requires inclusion of random velocity errors; we find best-fit RMS velocity errors of 399 km s$^{-1}$ and 252 km s$^{-1}$ for quasars and DLAs, respectively. We infer bias-weighted mean halo masses of $M_h/10^{12}\ h^{-1}M_\odot=2.19^{+0.16}_{-0.15}$ and $0.69^{+0.16}_{-0.14}$ for the host halos of quasars and DLAs, with $\sim 0.2$ dex systematic uncertainty associated with redshift evolution, IGM parameters, and selection of data fitting range.

preprint2016arXiv

Modelling galaxy clustering: halo occupation distribution versus subhalo matching

We model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $χ^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

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-$α$ 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.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

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= λ/Δλ$ 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.

preprint2016arXiv

The effects of assembly bias on cosmological inference from galaxy-galaxy lensing and galaxy clusters

The combination of galaxy-galaxy lensing (GGL) and galaxy clustering is a promising route to measuring the amplitude of matter clustering and testing modified gravity theories of cosmic acceleration. Halo occupation distribution (HOD) modeling can extend the approach down to nonlinear scales, but galaxy assembly bias could introduce systematic errors by causing the HOD to vary with large scale environment at fixed halo mass. We investigate this problem using the mock galaxy catalogs created by Hearin & Watson (2013, HW13), which exhibit significant assembly bias because galaxy luminosity is tied to halo peak circular velocity and galaxy colour is tied to halo formation time. The preferential placement of galaxies (especially red galaxies) in older halos affects the cutoff of the mean occupation function $\langle N_\text{cen}(M_\text{min}) \rangle$ for central galaxies, with halos in overdense regions more likely to host galaxies. The effect of assembly bias on the satellite galaxy HOD is minimal. We introduce an extended, environment dependent HOD (EDHOD) prescription to describe these results and fit galaxy correlation measurements. Crucially, we find that the galaxy-matter cross-correlation coefficient, $r_{gm} \equiv ξ_{gm} \cdot [ ξ_{mm} ξ_{gg} ]^{-1/2}$, is insensitive to assembly bias on scales $r \gtrsim 1 \; h^{-1}\text{Mpc}$, even though $ξ_{gm} $ and $ξ_{gg} $ are both affected individually. We can therefore recover the correct $ξ_{mm} $ from the HW13 galaxy-galaxy and galaxy-matter correlations using either a standard HOD or EDHOD fitting method. For $M_r \leq -19$ or $M_r \leq -20$ samples the recovery of $ξ_{mm}$ is accurate to 2% or better. For a sample of red $M_r \leq -20$ galaxies we achieve 2% recovery at $r \gtrsim 2\;h^{-1}\text{Mpc}$ with EDHOD modeling but lower accuracy at smaller scales or with a standard HOD fit.

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

preprint2015arXiv

Abundances, Stellar Parameters, and Spectra From the SDSS-III/APOGEE Survey

The SDSS-III/APOGEE survey operated from 2011-2014 using the APOGEE spectrograph, which collects high-resolution (R~22,500), near-IR (1.51-1.70 microns) spectra with a multiplexing (300 fiber-fed objects) capability. We describe the survey data products that are publicly available, which include catalogs with radial velocity, stellar parameters, and 15 elemental abundances for over 150,000 stars, as well as the more than 500,000 spectra from which these quantities are derived. Calibration relations for the stellar parameters (Teff, log g, [M/H], [alpha/M]) and abundances (C, N, O, Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, Ni) are presented and discussed. The internal scatter of the abundances within clusters indicates that abundance precision is generally between 0.05 and 0.09 dex across a broad temperature range; within more limited ranges and at high S/N, it is smaller for some elemental abundances. We assess the accuracy of the abundances using comparison of mean cluster metallicities with literature values, APOGEE observations of the solar spectrum and of Arcturus, comparison of individual star abundances with other measurements, and consideration of the locus of derived parameters and abundances of the entire sample, and find that it is challenging to determine the absolute abundance scale; external accuracy may be good to 0.1-0.2 dex. Uncertainties may be larger at cooler temperatures (Teff<4000K). Access to the public data release and data products is described, and some guidance for using the data products is provided.

preprint2015arXiv

ASPCAP: The Apogee Stellar Parameter and Chemical Abundances Pipeline

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) has built the largest moderately high-resolution (R=22, 500) spectroscopic map of the stars across the Milky Way, and including dust-obscured areas. The APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) is the software developed for the automated analysis of these spectra. ASPCAP determines atmospheric parameters and chemical abundances from observed spectra by comparing observed spectra to libraries of theoretical spectra, using chi-2 minimization in a multidimensional parameter space. The package consists of a fortran90 code that does the actual minimization, and a wrapper IDL code for book-keeping and data handling. This paper explains in detail the ASPCAP components and functionality, and presents results from a number of tests designed to check its performance. ASPCAP provides stellar effective temperatures, surface gravities, and metallicities precise to 2%, 0.1 dex, and 0.05 dex, respectively, for most APOGEE stars, which are predominantly giants. It also provides abundances for up to 15 chemical elements with various levels of precision, typically under 0.1 dex. The final data release (DR12) of the Sloan Digital Sky Survey III contains an APOGEE database of more than 150,000 stars. ASPCAP development continues in the SDSS-IV APOGEE-2 survey.

preprint2015arXiv

Baryon Cycling in the Low-Redshift Circumgalactic Medium: A Comparison of Simulations to the COS-Halos Survey

We analyze the low-redshift (z~0.2) circumgalactic medium by comparing absorption-line data from the COS-Halos Survey to absorption around a matched galaxy sample from two cosmological hydrodynamic simulations. The models include different prescriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), and constant winds (cw). We extract for comparison direct observables including equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are generally in good agreement with these observations for HI and certain low ionization metal lines, but show poorer agreement with higher ionization metal lines including SiIII and OVI that are well-observed by COS-Halos. These discrepancies suggest that both wind models predict too much cool, metal-enriched gas and not enough hot gas, and/or that the metals are not sufficiently well-mixed. This may reflect our model assumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a heuristic prescription to quench massive galaxies by super-heating its ISM gas, which we show yields sufficient low ionisation absorption to be broadly consistent with observations, but also substantial OVI absorption that is inconsistent with data, suggesting that gas around quenched galaxies in the real Universe does not cool. At impact parameters of <50 kpc, recycling winds dominate the absorption of low ions and even HI, while OVI almost always arises from metals ejected longer than 1 Gyr ago. The similarity between the wind models is surprising, since we show that they differ substantially in their predicted amount and phase distribution of halo gas. We show that this similarity owes mainly to our comparison here at fixed stellar mass rather than at fixed halo mass in our previous works, which suggests that CGM properties are more closely tied to the stellar mass of galaxies rather than halo mass.

preprint2015arXiv

Chemical Cartography with APOGEE: Metallicity Distribution Functions and the Chemical Structure of the Milky Way Disk

Using a sample of 69,919 red giants from the SDSS-III/APOGEE Data Release 12, we measure the distribution of stars in the [$α$/Fe] vs. [Fe/H] plane and the metallicity distribution functions (MDF) across an unprecedented volume of the Milky Way disk, with radius $3<R<15$ kpc and height $|z|<2$ kpc. Stars in the inner disk ($R<5$ kpc) lie along a single track in [$α$/Fe] vs. [Fe/H], starting with $α$-enhanced, metal-poor stars and ending at [$α$/Fe]$\sim0$ and [Fe/H]$\sim+0.4$. At larger radii we find two distinct sequences in [$α$/Fe] vs. [Fe/H] space, with a roughly solar-$α$ sequence that spans a decade in metallicity and a high-$α$ sequence that merges with the low-$α$ sequence at super-solar [Fe/H]. The location of the high-$α$ sequence is nearly constant across the disk, however there are very few high-$α$ stars at $R>11$ kpc. The peak of the midplane MDF shifts to lower metallicity at larger $R$, reflecting the Galactic metallicity gradient. Most strikingly, the shape of the midplane MDF changes systematically with radius, with a negatively skewed distribution at $3<R<7$ kpc, to a roughly Gaussian distribution at the solar annulus, to a positively skewed shape in the outer Galaxy. For stars with $|z|>1$ kpc or [$α$/Fe]$>0.18$, the MDF shows little dependence on $R$. The positive skewness of the outer disk MDF may be a signature of radial migration; we show that blurring of stellar populations by orbital eccentricities is not enough to explain the reversal of MDF shape but a simple model of radial migration can do so.

preprint2015arXiv

Clustering of intermediate redshift quasars using the final SDSS III-BOSS sample

We measure the two-point clustering of spectroscopically confirmed quasars from the final sample of the Baryon Oscillation Spectroscopic Survey (BOSS) on comoving scales of 4 < s < 22 Mpc/h. The sample covers 6950 deg^2 (~ 19 (Gpc/h)^3) and, over the redshift range 2.2 < z < 2.8, contains 55,826 homogeneously selected quasars, which is twice as many as in any similar work. We deduce b_Q = 3.54 +/- 0.10 ; the most precise measurement of quasar bias to date at these redshifts. This corresponds to a host halo mass of ~ 2 x 10^12 ~ M_sun/h with an implied quasar duty cycle of ~1 percent. The real-space projected correlation function is well-fit by a power law of index -2 and correlation length r0 = (8.12 +/- 0.22), Mpc/h over scales of 4 < rp < 25 ~ Mpc/h. To better study the evolution of quasar clustering at moderate redshift, we extend the redshift range of our study to z ~ 3.4 and measure the bias and correlation length of three subsamples over 2.2 < z < 3.4. We find no significant evolution of r0 or bias over this range, implying that the host halo mass of quasars decreases somewhat with increasing redshift. We find quasar clustering remains similar over a decade in luminosity, contradicting a scenario in which quasar luminosity is monotonically related to halo mass at z ~ 2.5. Our results are broadly consistent with previous BOSS measurements, but they yield more precise constraints based upon a larger and more uniform data set.

preprint2015arXiv

Cosmological implications of baryon acoustic oscillation (BAO) measurements

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $Ω_m=0.301 \pm 0.008$ and curvature $Ω_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $χ^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

preprint2015arXiv

Large-scale clustering of Lyman-alpha emission intensity from SDSS/BOSS

(Abridged) We detect the large-scale structure of Lya emission in the Universe at redshifts z=2-3.5 by measuring the cross-correlation of Lya surface brightness with quasars in SDSS/BOSS. We use a million spectra targeting Luminous Red Galaxies at z<0.8, after subtracting a best fit model galaxy spectrum from each one, as an estimate of the high-redshift Lya surface brightness. The quasar-Lya emission cross-correlation we detect has a shape consistent with a LambdaCDM model with Omega_M =0.30^+0.10-0.07. The predicted amplitude of this cross-correlation is proportional to the product of the mean Lya surface brightness, <mu_alpha>, the amplitude of mass fluctuations, and the quasar and Lya emission bias factors. Using known values, we infer <mu_alpha>(b_alpha/3) = (3.9 +/- 0.9) x 10^-21 erg/s cm^-2 A^-1 arcsec^-2, where b_alpha is the Lya emission bias factor. If the dominant sources of Lya emission are star forming galaxies, we infer rho_SFR = (0.28 +/- 0.07) (3/b_alpha) /yr/Mpc^3 at z=2-3.5. For b_alpha=3, this value is a factor of 21-35 above previous estimates from individually detected Lya emitters, although consistent with the total rho_SFR derived from dust-corrected, continuum UV surveys. 97% of the Lya emission in the Universe at these redshifts is therefore undetected in previous surveys of Lya emitters. Our measurement is much greater than seen from stacking analyses of faint halos surrounding previously detected Lya emitters, but we speculate that it arises from similar Lya halos surrounding all luminous star-forming galaxies. We also detect redshift space anisotropy of the quasar-Lya emission cross-correlation, finding evidence at the 3.0 sigma level that it is radially elongated, consistent with distortions caused by radiative-transfer effects (Zheng et al. (2011)). Our measurements represent the first application of the intensity mapping technique to optical observations.

preprint2015arXiv

Modelling The Redshift-Space Three-Point Correlation Function in SDSS-III

We present the measurements of the redshift-space three-point correlation function (3PCF) for z~0.5 luminous red galaxies of the CMASS sample in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. The 3PCF measurements are interpreted within the halo occupation distribution framework using high-resolution N-body simulations, and the model successfully reproduces the 3PCF on scales larger than 1Mpc/h. As with the case for the redshift-space two-point correlation functions, we find that the redshift-space 3PCF measurements also favour the inclusion of galaxy velocity bias in the model. In particular, the central galaxy in a halo is on average in motion with respect to the core of the halo. We discuss the potential of the small-scale 3PCF to tighten the constraints on the relation between galaxies and dark matter haloes and on the phase-space distribution of galaxies.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

preprint2014arXiv

A Budget and Accounting of Metals at z~0: Results from the COS-Halos Survey

We present a budget and accounting of metals in and around star-forming galaxies at $z\sim 0$. We combine empirically derived star formation histories with updated supernova and AGB yields and rates to estimate the total mass of metals produced by galaxies with present-day stellar mass of $10^{9.3}$--$10^{11.6} M_{\odot}$. On the accounting side of the ledger, we show that a surprisingly constant 20--25% mass fraction of produced metals remain in galaxies' stars, interstellar gas and interstellar dust, with little dependence of this fraction on the galaxy stellar mass (omitting those metals immediately locked up in remnants). Thus, the bulk of metals are outside of galaxies, produced in the progenitors of today's $L^*$ galaxies. The COS-Halos survey is uniquely able to measure the mass of metals in the circumgalactic medium (to impact parameters of $< 150$ kpc) of low-redshift $\sim L^*$ galaxies. Using these data, we map the distribution of CGM metals as traced by both the highly ionized OVI ion and a suite of low-ionization species; combined with constraints on circumgalactic dust and hotter X-ray emitting gas out to similar impact parameters, we show that $\sim 40$% of metals produced by $M_{\star}\sim 10^{10}M_{\odot}$ galaxies can be easily accounted for out to 150 kpc. With the current data, we cannot rule out a constant mass of metals within this fixed physical radius. This census provides a crucial boundary condition for the eventual fate of metals in galaxy evolution models.

preprint2014arXiv

A measurement of the Alcock-Paczynski effect using cosmic voids in the SDSS

We perform an Alcock-Paczynski test using stacked cosmic voids identified in the SDSS Data Release 7 main sample and Data Release 10 LOWZ and CMASS samples. We find ~1,500 voids out to redshift $0.6$ using a heavily modified and extended version of the watershed algorithm ZOBOV, which we call VIDE (Void IDentification and Examination). To assess the impact of peculiar velocities we use the mock void catalogs presented in Sutter et al. (2013). We find a constant uniform flattening of 14% along the line of sight when peculiar velocities are included. This flattening appears universal for all void sizes at all redshifts and for all tracer densities. We also use these mocks to identify an optimal stacking strategy. After correcting for systematic effects we find that our Alcock-Paczynski measurement leads to a preference of our best-fit value of $Ω_{\rm M}\sim 0.15$ over $Ω_{\rm M} = 1.0$ by a likelihood ratio of 10. Likewise, we find a factor of $4.5$ preference of the likelihood ratio for a $Λ$CDM $Ω_{\rm M} = 0.3$ model and a null measurement. Taken together, we find substantial evidence for the Alcock-Paczynski signal in our sample of cosmic voids. Our assessment using realistic mocks suggests that measurements with future SDSS releases and other surveys will provide tighter cosmological parameter constraints. The void-finding algorithm and catalogs used in this work will be made publicly available at http://www.cosmicvoids.net.

preprint2014arXiv

Baryon Acoustic Oscillations in the Lyα forest of BOSS DR11 quasars

We report a detection of the baryon acoustic oscillation (BAO) feature in the flux-correlation function of the Lyα forest of high-redshift quasars with a statistical significance of five standard deviations. The study uses 137,562 quasars in the redshift range $2.1\le z \le 3.5$ from the Data Release 11 (DR11) of the Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS-III. This sample contains three times the number of quasars used in previous studies. The measured position of the BAO peak determines the angular distance, $D_A(z=2.34)$ and expansion rate, $H(z=2.34)$, both on a scale set by the sound horizon at the drag epoch, $r_d$. We find $D_A/r_d=11.28\pm0.65(1σ)^{+2.8}_{-1.2}(2σ)$ and $D_H/r_d=9.18\pm0.28(1σ)\pm0.6(2σ)$ where $D_H=c/H$. The optimal combination, $\sim D_H^{0.7}D_A^{0.3}/r_d$ is determined with a precision of $\sim2\%$. For the value $r_d=147.4~{\rm Mpc}$, consistent with the CMB power spectrum measured by Planck, we find $D_A(z=2.34)=1662\pm96(1σ)~{\rm Mpc}$ and $H(z=2.34)=222\pm7(1σ)~{\rm km\,s^{-1}Mpc^{-1}}$. Tests with mock catalogs and variations of our analysis procedure have revealed no systematic uncertainties comparable to our statistical errors. Our results agree with the previously reported BAO measurement at the same redshift using the quasar-Lyα forest cross-correlation. The auto-correlation and cross-correlation approaches are complementary because of the quite different impact of redshift-space distortion on the two measurements. The combined constraints from the two correlation functions imply values of $D_A/r_d$ and $D_H/r_d$ that are, respectively, 7% low and 7% high compared to the predictions of a flat $Λ$CDM cosmological model with the best-fit Planck parameters. With our estimated statistical errors, the significance of this discrepancy is $\approx 2.5σ$.

preprint2014arXiv

IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function

The Ly$α$ forest transmission probability distribution function (PDF) is an established probe of the intergalactic medium (IGM) astrophysics, especially the temperature-density relationship of the IGM. We measure the transmission PDF from 3393 Baryon Oscillations Spectroscopic Survey (BOSS) quasars from SDSS Data Release 9, and compare with mock spectra that include careful modeling of the noise, continuum, and astrophysical uncertainties. The BOSS transmission PDFs, measured at $\langle z \rangle = [2.3,2.6,3.0]$, are compared with PDFs created from mock spectra drawn from a suite of hydrodynamical simulations that sample the IGM temperature-density relationship, $γ$, and temperature at mean-density, $T_0$, where $T(Δ) = T_0 Δ^{γ-1}$. We find that a significant population of partial Lyman-limit systems with a column-density distribution slope of $β_\mathrm{pLLS} \sim -2$ are required to explain the data at the low-transmission end of transmission PDF, while uncertainties in the mean Ly$α$ forest transmission affect the high-transmission end. After modelling the LLSs and marginalizing over mean-transmission uncertainties, we find that $γ=1.6$ best describes the data over our entire redshift range, although constraints on $T_0$ are affected by systematic uncertainties. Within our model framework, isothermal or inverted temperature-density relationships ($γ\leq 1$) are disfavored at a significance of over 4$σ$, although this could be somewhat weakened by cosmological and astrophysical uncertainties that we did not model.

preprint2014arXiv

LyMAS: Predicting Large-Scale Lyman-alpha Forest Statistics from the Dark Matter Density Field

[abridged] We describe LyMAS (Ly-alpha Mass Association Scheme), a method of predicting clustering statistics in the Ly-alpha forest on large scales from moderate resolution simulations of the dark matter distribution, with calibration from high-resolution hydrodynamic simulations of smaller volumes. We use the "Horizon MareNostrum" simulation, a 50 Mpc/h comoving volume evolved with the adaptive mesh hydrodynamic code RAMSES, to compute the conditional probability distribution P(F_s|delta_s) of the transmitted flux F_s, smoothed (1-dimensionally) over the spectral resolution scale, on the dark matter density contrast delta_s, smoothed (3-dimensionally) over a similar scale. In this study we adopt the spectral resolution of the SDSS-III BOSS at z=2.5, and we find optimal results for a dark matter smoothing length sigma=0.3 Mpc/h (comoving). In extended form, LyMAS exactly reproduces both the 1-dimensional power spectrum and 1-point flux distribution of the hydro simulation spectra. Applied to the MareNostrum dark matter field, LyMAS accurately predicts the 2-point conditional flux distribution and flux correlation function of the full hydro simulation for transverse sightline separations as small as 1 Mpc/h, including redshift-space distortion effects. It is substantially more accurate than a deterministic density-flux mapping ("Fluctuating Gunn-Peterson Approximation"), often used for large volume simulations of the forest. With the MareNostrum calibration, we apply LyMAS to 1024^3 N-body simulations of a 300 Mpc/h and 1.0 Gpc/h cube to produce large, publicly available catalogs of mock BOSS spectra that probe a large comoving volume. LyMAS will be a powerful tool for interpreting 3-d Ly-alpha forest data, thereby transforming measurements from BOSS and other massive quasar absorption surveys into constraints on dark energy, dark matter, space geometry, and IGM physics.

preprint2014arXiv

Parametrising Star Formation Histories

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

preprint2014arXiv

Probing the Circumgalactic Medium at High-Redshift Using Composite BOSS Spectra of Strong Lyman-alpha Forest Absorbers

We present composite spectra constructed from a sample of 242,150 Lyman-alpha (Lya) forest absorbers at redshifts 2.4<z<3.1 identified in quasar spectra from the Baryon Oscillation Spectroscopic Survey (BOSS) as part of Data Release 9 of the Sloan Digital Sky Survey III. We select forest absorbers by their flux in bins 138 km/s wide (approximately the size of the BOSS resolution element). We split these absorbers into five samples spanning the range of flux -0.05 < F<0.45. Tests on a smaller sample of high-resolution spectra show that our three strongest absorption bins would probe circumgalactic regions (projected separation < 300 proper kpc and |Delta v| < 300km/s) in about 60% of cases for very high signal-to-noise ratio. Within this subset, weakening Lya absorption is associated with decreasing purity of circumgalactic selection once BOSS noise is included. Our weaker two Lya absorption samples are dominated by the intergalactic medium. We present composite spectra of these samples and a catalogue of measured absorption features from HI and 13 metal ionization species, all of which we make available to the community. We compare measurements of seven Lyman series transitions in our composite spectra to single line models and obtain further constraints from their associated excess Lyman limit opacity. This analysis provides results consistent with column densities over the range 14.4 <~ Log (N_HI) <~ 16.45. We compare our measurements of metal absorption to a variety of simple single-line, single-phase models for a preliminary interpretation. Our results imply clumping on scales down to ~30 pc and near-solar metallicities in the circumgalactic samples, while high-ionization metal absorption consistent with typical IGM densities and metallicities is visible in all samples.

preprint2014arXiv

Sparse sampling, galaxy bias, and voids

To study the impact of sparsity and galaxy bias on void statistics, we use a single large-volume, high-resolution N-body simulation to compare voids in multiple levels of subsampled dark matter, halo populations, and mock galaxies from a Halo Occupation Distribution model tuned to different galaxy survey densities. We focus our comparison on three key observational statistics: number functions, ellipticity distributions, and radial density profiles. We use the hierarchical tree structure of voids to interpret the impacts of sampling density and galaxy bias, and theoretical and empirical functions to describe the statistics in all our sample populations. We are able to make simple adjustments to theoretical expectations to offer prescriptions for translating from analytics to the void properties measured in realistic observations. We find that sampling density has a much larger effect on void sizes than galaxy bias. At lower tracer density, small voids disappear and the remaining voids are larger, more spherical, and have slightly steeper profiles. When a proper lower mass threshold is chosen, voids in halo distributions largely mimic those found in galaxy populations, except for ellipticities, where galaxy bias leads to higher values. We use the void density profile of Hamaus et al. (2014) to show that voids follow a self-similar and universal trend, allowing simple translations between voids studied in dark matter and voids identified in galaxy surveys. We have added the mock void catalogs used in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

preprint2014arXiv

The APOGEE red-clump catalog: Precise distances, velocities, and high-resolution elemental abundances over a large area of the Milky Way's disk

The Sloan Digital Sky Survey III's Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution near-infrared spectroscopic survey covering all of the major components of the Galaxy, including the dust-obscured regions of the inner Milky Way disk and bulge. Here we present a sample of 10,341 likely red-clump stars (RC) from the first two years of APOGEE operations, selected based on their position in color-metallicity-surface-gravity-effective-temperature space using a new method calibrated using stellar-evolution models and high-quality asteroseismology data. The narrowness of the RC locus in color-metallicity-luminosity space allows us to assign distances to the stars with an accuracy of 5 to 10%. The sample extends to typical distances of about 3 kpc from the Sun, with some stars out to 8 kpc, and spans a volume of approximately 100 kpc^3 over 5 kpc <~ R <~ 14 kpc, |Z| <~ 2 kpc, and -15 deg <~ Galactocentric azimuth <~ 30 deg. The APOGEE red-clump (APOGEE-RC) catalog contains photometry from 2MASS, reddening estimates, distances, line-of-sight velocities, stellar parameters and elemental abundances determined from the high-resolution APOGEE spectra, and matches to major proper motion catalogs. We determine the survey selection function for this data set and discuss how the RC selection samples the underlying stellar populations. We use this sample to limit any azimuthal variations in the median metallicity within the ~45 degree-wide azimuthal region covered by the current sample to be <= 0.02 dex, which is more than an order of magnitude smaller than the radial metallicity gradient. This result constrains coherent non-axisymmetric flows within a few kpc from the Sun.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $Λ$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,σ$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: modeling of the luminosity and colour dependence in the Data Release 10

We investigate the luminosity and colour dependence of clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The halo occupation distribution framework is adopted to model the projected two-point correlation function measurements on small and intermediate scales (from $0.02$ to $60\,h^{-1}{\rm {Mpc}}$) and to interpret the observed trends and infer the connection of galaxies to dark matter halos. We find that luminous red galaxies reside in massive halos of mass $M{\sim}10^{13}$--$10^{14}\,h^{-1}{\rm M_\odot}$ and more luminous galaxies are more clustered and hosted by more massive halos. The strong small-scale clustering requires a fraction of these galaxies to be satellites in massive halos, with the fraction at the level of 5--8 per cent and decreasing with luminosity. The characteristic mass of a halo hosting on average one satellite galaxy above a luminosity threshold is about a factor $8.7$ larger than that of a halo hosting a central galaxy above the same threshold. At a fixed luminosity, progressively redder galaxies are more strongly clustered on small scales, which can be explained by having a larger fraction of these galaxies in the form of satellites in massive halos. Our clustering measurements on scales below $0.4\,h^{-1}{\rm {Mpc}}$ allow us to study the small-scale spatial distribution of satellites inside halos. While the clustering of luminosity-threshold samples can be well described by a Navarro-Frenk-White (NFW) profile, that of the reddest galaxies prefers a steeper or more concentrated profile. Finally, we also use galaxy samples of constant number density at different redshifts to study the evolution of luminous galaxies, and find the clustering to be consistent with passive evolution in the redshift range of $0.5 \lesssim z \lesssim 0.6$.

preprint2014arXiv

The COS-Dwarfs Survey: The Carbon Reservoir Around sub-L* Galaxies

We report new observations of circumgalactic gas from the COS-Dwarfs survey, a systematic investigation of the gaseous halos around 43 low-mass z $\leq$ 0.1 galaxies using background QSOs observed with the Cosmic Origins Spectrograph. From the projected 1D and 2D distribution of C IV absorption, we find that C IV absorption is detected out to ~ 0.5 R$_{vir}$ of the host galaxies. The C IV absorption strength falls off radially as a power law and beyond 0.5 R$_{vir}$, no C IV absorption is detected above our sensitivity limit of ~ 50-100 m$Å$. We find a tentative correlation between detected C IV absorption strength and star formation, paralleling the strong correlation seen in highly ionized oxygen for L~L* galaxies by the COS-Halos survey. The data imply a large carbon reservoir in the CGM of these galaxies, corresponding to a minimum carbon mass of $\gtrsim$ 1.2$\times 10^6$ $M_\odot$ out to ~ 110 kpc. This mass is comparable to the carbon mass in the ISM and more than the carbon mass currently in stars of these galaxies. The C IV absorption seen around these sub-L* galaxies can account for almost two-thirds of all $W_r$> 100 m$Å$ C IV absorption detected at low z. Comparing the C IV covering fraction with hydrodynamical simulations, we find that an energy-driven wind model is consistent with the observations whereas a wind model of constant velocity fails to reproduce the CGM or the galaxy properties.

preprint2014arXiv

The COS-Halos Survey: Physical Conditions and Baryonic Mass in the Low-Redshift Circumgalactic Medium

We analyze the physical conditions of the cool, photoionized (T $\sim 10^4$ K) circumgalactic medium (CGM) using the COS-Halos suite of gas column density measurements for 44 gaseous halos within 160 kpc of $L \sim L^*$ galaxies at $z \sim 0.2$. These data are well described by simple photoionization models, with the gas highly ionized (n$_{\rm HII}$/n$_{\rm H} \gtrsim 99\%$) by the extragalactic ultraviolet background (EUVB). Scaling by estimates for the virial radius, R$_{\rm vir}$, we show that the ionization state (tracked by the dimensionless ionization parameter, U) increases with distance from the host galaxy. The ionization parameters imply a decreasing volume density profile n$_{\rm H}$ = (10$^{-4.2 \pm 0.25}$)(R/R$_{\rm vir})^{-0.8\pm0.3}$. Our derived gas volume densities are several orders of magnitude lower than predictions from standard two-phase models with a cool medium in pressure equilibrium with a hot, coronal medium expected in virialized halos at this mass scale. Applying the ionization corrections to the HI column densities, we estimate a lower limit to the cool gas mass M$_{\rm CGM}^{\rm cool} > 6.5 \times 10^{10}$ M$_{\odot}$ for the volume within R $<$ R$_{\rm vir}$. Allowing for an additional warm-hot, OVI-traced phase, the CGM accounts for at least half of the baryons purported to be missing from dark matter halos at the 10$^{12}$ M$_{\odot}$ scale.

preprint2014arXiv

The dark matter of galaxy voids

How do observed voids relate to the underlying dark matter distribution? To examine the spatial distribution of dark matter contained within voids identified in galaxy surveys, we apply Halo Occupation Distribution models representing sparsely and densely sampled galaxy surveys to a high-resolution N-body simulation. We compare these galaxy voids to voids found in the halo distribution, low-resolution dark matter, and high-resolution dark matter. We find that voids at all scales in densely sampled surveys - and medium- to large-scale voids in sparse surveys - trace the same underdensities as dark matter, but they are larger in radius by ~20%, they have somewhat shallower density profiles, and they have centers offset by ~0.4Rv rms. However, in void-to-void comparison we find that shape estimators are less robust to sampling, and the largest voids in sparsely sampled surveys suffer fragmentation at their edges. We find that voids in galaxy surveys always correspond to underdensities in the dark matter, though the centers may be offset. When this offset is taken into account, we recover almost identical radial density profiles between galaxies and dark matter. All mock catalogs used in this work are available at http://www.cosmicvoids.net.

preprint2014arXiv

The Photon Underproduction Crisis

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

preprint2014arXiv

Tracing chemical evolution over the extent of the Milky Way's Disk with APOGEE Red Clump Stars

We employ the first two years of data from the near-infrared, high-resolution SDSS-III/APOGEE spectroscopic survey to investigate the distribution of metallicity and alpha-element abundances of stars over a large part of the Milky Way disk. Using a sample of ~10,000 kinematically-unbiased red-clump stars with ~5% distance accuracy as tracers, the [alpha/Fe] vs. [Fe/H] distribution of this sample exhibits a bimodality in [alpha/Fe] at intermediate metallicities, -0.9<[Fe/H]<-0.2, but at higher metallicities ([Fe/H]=+0.2) the two sequences smoothly merge. We investigate the effects of the APOGEE selection function and volume filling fraction and find that these have little qualitative impact on the alpha-element abundance patterns. The described abundance pattern is found throughout the range 5<R<11 kpc and 0<|Z|<2 kpc across the Galaxy. The [alpha/Fe] trend of the high-alpha sequence is surprisingly constant throughout the Galaxy, with little variation from region to region (~10%). Using simple galactic chemical evolution models we derive an average star formation efficiency (SFE) in the high-alpha sequence of ~4.5E-10 1/yr, which is quite close to the nearly-constant value found in molecular-gas-dominated regions of nearby spirals. This result suggests that the early evolution of the Milky Way disk was characterized by stars that shared a similar star formation history and were formed in a well-mixed, turbulent, and molecular-dominated ISM with a gas consumption timescale (1/SFE) of ~2 Gyr. Finally, while the two alpha-element sequences in the inner Galaxy can be explained by a single chemical evolutionary track this cannot hold in the outer Galaxy, requiring instead a mix of two or more populations with distinct enrichment histories.

preprint2014arXiv

Velocity Bias from the Small Scale Clustering of SDSS-III BOSS Galaxies

We present the measurements and modelling of the projected and redshift-space clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. For a volume-limited luminous red galaxy sample in the redshift range of $0.48<z<0.55$, we perform halo occupation distribution modelling of the small- and intermediate-scale ($0.1$--$60h^{-1}{\rm {Mpc}}$) projected and redshift-space two-point correlation functions, with an accurate model built on high resolution $N$-body simulations. To interpret the measured redshift-space distortions, the distribution of galaxy velocities must differ from that of the dark matter inside haloes of $\sim 10^{13}$--$10^{14}h^{-1}{\rm M_{\odot}}$, i.e. the data require the existence of galaxy velocity bias. Most notably, central galaxies on average are not at rest with respect to the core of their host haloes (defined by the inner 25% of particles around the halo potential minimum), but rather move around it with a 1D velocity dispersion of $0.22^{+0.03}_{-0.04}$ times that of the dark matter, implying a spatial offset from the centre at the level of $\lesssim$1% of the halo virial radius. The luminous satellite galaxies move more slowly than the dark matter, with velocities $0.86^{+0.08}_{-0.03}$ times those of the dark matter, which suggests that the velocity and spatial distributions of these satellites cannot both be unbiased. The constraints mainly arise from the Fingers-of-God effect at nonlinear scales and the smoothing to the Kaiser effect in the translinear regime; the robustness of the results is demonstrated by a variety of tests. In addition, no clear evidence is found for a strong luminosity dependence of the velocity bias. We discuss the implications of the existence of galaxy velocity bias for investigations of galaxy formation and cosmology.

preprint2014arXiv

Voids in the SDSS DR9: observations, simulations, and the impact of the survey mask

We present and study cosmic voids identified using the watershed void finder VIDE in the Sloan Digital Sky Survey Data Release 9, compare these voids to ones identified in mock catalogs, and assess the impact of the survey mask on void statistics such as number functions, ellipticity distributions, and radial density profiles. The nearly 1,000 identified voids span three nearly volume-limited samples from redshift z = 0.43 to 0.7. For comparison we use 98 of the publicly available 2LPT-based mock galaxy catalogs of Manera et al., and also generate our own mock catalogs by applying a Halo Occupation Distribution model to an N-body simulation. We find that the mask reduces the number density of voids at all scales by a factor of three and slightly skews the relative size distributions. This engenders an increase in the mean ellipticity by roughly 30%. However, we find that radial density profiles are largely robust to the effects of the mask. We see excellent agreement between the data and both mock catalogs, and find no tension between the observed void properties and the properties derived from ΛCDM simulations. We have added the void catalogs from both data and mock galaxy populations discussed in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

preprint2013arXiv

A response to arXiv:1310.2791: A self-consistent public catalogue of voids and superclusters in the SDSS Data Release 7 galaxy surveys

Recently, Nadathur & Hotchkiss (2013) submitted a paper discussing a new cosmic void catalog. This paper includes claims about the void catalog described in Sutter et al. (2012). In this note, we respond to those claims, clarify some discrepancies between the text of Sutter et al. (2012) and the most recent version of the catalog, and provide some comments on the differences between our catalog and that of Nadathur & Hotchkiss (2013). All updates and documentation for our catalog are available at http://www.cosmicvoids.net.

preprint2013arXiv

Baryon Acoustic Oscillations in the Ly-α forest of BOSS quasars

We report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted flux fraction in the \Lya forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range $2.1\le z \le 3.5$ from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean redshift $z=2.3$, we measure the monopole and quadrupole components of the correlation function for separations in the range $20\hMpc<r<200\hMpc$. A peak in the correlation function is seen at a separation equal to $(1.01\pm0.03)$ times the distance expected for the BAO peak within a concordance $Λ$CDM cosmology. This first detection of the BAO peak at high redshift, when the universe was strongly matter dominated, results in constraints on the angular diameter distance $\da$ and the expansion rate $H$ at $z=2.3$ that, combined with priors on $H_0$ and the baryon density, require the existence of dark energy. Combined with constraints derived from Cosmic Microwave Background (CMB) observations, this result implies $H(z=2.3)=(224\pm8){\rm km\,s^{-1}Mpc^{-1}}$, indicating that the time derivative of the cosmological scale parameter $\dot{a}=H(z=2.3)/(1+z)$ is significantly greater than that measured with BAO at $z\sim0.5$. This demonstrates that the expansion was decelerating in the range $0.7<z<2.3$, as expected from the matter domination during this epoch. Combined with measurements of $H_0$, one sees the pattern of deceleration followed by acceleration characteristic of a dark-energy dominated universe.

preprint2013arXiv

Cold dark matter: controversies on small scales

The cold dark matter (CDM) cosmological model has been remarkably successful in explaining cosmic structure over an enormous span of redshift, but it has faced persistent challenges from observations that probe the innermost regions of dark matter halos and the properties of the Milky Way's dwarf galaxy satellites. We review the current observational and theoretical status of these "small scale controversies." Cosmological simulations that incorporate only gravity and collisionless CDM predict halos with abundant substructure and central densities that are too high to match constraints from galaxy dynamics. The solution could lie in baryonic physics: recent numerical simulations and analytic models suggest that gravitational potential fluctuations tied to efficient supernova feedback can flatten the central cusps of halos in massive galaxies, and a combination of feedback and low star-formation efficiency could explain why most of the dark matter subhalos orbiting the Milky Way do not host visible galaxies. However, it is not clear that this solution can work in the lowest mass galaxies where discrepancies are observed. Alternatively, the small-scale conflicts could be evidence of more complex physics in the dark sector itself. For example, elastic scattering from strong dark matter self-interactions can alter predicted halo mass profiles, leading to good agreement with observations across a wide range of galaxy mass. Gravitational lensing and dynamical perturbations of tidal streams in the stellar halo provide evidence for an abundant population of low mass subhalos in accord with CDM predictions. These observational approaches will get more powerful over the next few years.

preprint2013arXiv

Dark Energy: A Short Review

The accelerating expansion of the universe is the most surprising cosmological discovery in many decades. In this short review, we briefly summarize theories for the origin of cosmic acceleration and the observational methods being used to test these theories. We then discuss the current observational state of the field, with constraints from the cosmic microwave background (CMB), baryon acoustic oscillations (BAO), Type Ia supernovae (SN), direct measurements of the Hubble constant ($H_0$), and measurements of galaxy and matter clustering. Assuming a flat universe and dark energy with a constant equation-of-state parameter $w = P/ρ$, the combination of Planck CMB temperature anisotropies, WMAP CMB polarization, the Union2.1 SN compilation, and a compilation of BAO measurements yields $w = -1.10^{+0.08}_{-0.07}$, consistent with a cosmological constant ($w=-1$). However, with these constraints the cosmological constant model predicts a value of $H_0$ that is lower than several of the leading recent estimates, and it predicts a parameter combination $σ_8(Ω_m)^{0.5}$ that is higher than many estimates from weak gravitational lensing, galaxy clusters, and redshift-space distortions. Individually these tensions are only significant at the ~$2σ$ level, but they arise in multiple data sets with independent statistics and distinct sources of systematic uncertainty. The tensions are stronger with Planck CMB data than they were with WMAP because of the smaller statistical errors and the higher central value of $Ω_m.$ With the improved measurements expected from the next generation of data sets, these tensions may diminish, or they may sharpen in a way that points towards a more complete physical understanding of cosmic acceleration.

preprint2013arXiv

Fitting Methods for Baryon Acoustic Oscillations in the Lyman-α Forest Fluctuations in BOSS Data Release 9

We describe fitting methods developed to analyze fluctuations in the Lyman-α forest and measure the parameters of baryon acoustic oscillations (BAO). We apply our methods to BOSS Data Release 9. Our method is based on models of the three-dimensional correlation function in physical coordinate space, and includes the effects of redshift-space distortions, anisotropic non-linear broadening, and broadband distortions. We allow for independent scale factors along and perpendicular to the line of sight to minimize the dependence on our assumed fiducial cosmology and to obtain separate measurements of the BAO angular and relative velocity scales. Our fitting software and the input files needed to reproduce our main BOSS Data Release 9 results are publicly available.

preprint2013arXiv

Hydrogen and Metal Line Absorption Around Low-Redshift Galaxies in Cosmological Hydrodynamic Simulations

We study the physical conditions of the circum-galactic medium (CGM) around z=0.25 galaxies as traced by HI and metal line absorption, using cosmological hydrodynamic simulations that include galactic outflows. Using lines of sight targeted at impact parameters from 10 kpc to 1 Mpc around galaxies with halo masses from 10^11-10^13 M_solar, we study the physical conditions and their variation with impact parameter b and line-of-sight velocity delta v in the CGM as traced by HI, MgII, SiIV, CIV, OVI, and NeVIII absorbers. All ions show a strong excess of absorption near galaxies compared to random lines of sight. The excess continues beyond 1 Mpc, reflecting the correlation of metal absorption with large-scale structure. Absorption is particularly enhanced within about v<300 km/sec and roughly 300 kpc of galaxies (with distances somewhat larger for the highest ion), approximately delineating the CGM; this range contains the majority of global metal absorption. Low ions like MgII and SiIV predominantly arise in denser gas closer to galaxies and drop more rapidly with b, while high ions OVI and NeVIII trace more diffusely distributed gas with a comparatively flat radial profile; CIV is intermediate. All ions predominantly trace T~10^4-4.5 K photo-ionised gas at all b, but when hot CGM gas is present (mostly in larger halos), we see strong collisionally-ionised OVI and NeVIII at b <= 100 kpc. Larger halo masses generally produce more absorption, though overall the trends are not as strong as that with impact parameter. These findings arise using our favoured outflow scalings as expected for momentum-driven winds; with no winds, the CGM gas remains mostly unenriched, while our outflow model with a constant velocity and mass loading factor produce hotter, more widely dispersed metals.

preprint2013arXiv

Inside Out and Upside Down: Tracing the Assembly of a Simulated Disk Galaxy Using Mono-Age Stellar Populations

We analyze the present-day structure and assembly history of a high resolution hydrodynamic simulation of the formation of a Milky Way (MW)-like disk galaxy, from the "Eris" simulation suite, dissecting it into cohorts of stars formed at different epochs of cosmic history. At z=0, stars with t_form < 2 Gyr mainly occupy the stellar spheroid, with the oldest (earliest forming) stars having more centrally concentrated profiles. The younger age cohorts populate disks of progressively longer radial scale length and shorter vertical scale height. At a given radius, the vertical density profiles and velocity dispersions of stars vary smoothly as a function of age, and the superposition of old, vertically-extended and young, vertically-compact cohorts gives rise to a double-exponential profile like that observed in the MW. Turning to formation history, we find that the trends of spatial structure and kinematics with stellar age are largely imprinted at birth, or immediately thereafter. Stars that form during the active merger phase at z>3 are quickly scattered into rounded, kinematically hot configurations. The oldest disk cohorts form in structures that are radially compact and relatively thick, while subsequent cohorts form in progressively larger, thinner, colder configurations from gas with increasing levels of rotational support. The disk thus forms "inside-out" in a radial sense and "upside-down" in a vertical sense. Secular heating and radial migration influence the final state of each age cohort, but the changes they produce are small compared to the trends established at formation. The predicted correlations of stellar age with spatial and kinematic structure are in good qualitative agreement with the correlations observed for mono-abundance stellar populations in the MW.

preprint2013arXiv

Measurement of Baryon Acoustic Oscillations in the Lyman-alpha Forest Fluctuations in BOSS Data Release 9

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) to detect and measure the position of the Baryonic Acoustic Oscillation (BAO) feature in the three-dimensional correlation function in the Lyman-alpha forest flux fluctuations at a redshift z=2.4. The feature is clearly detected at significance between 3 and 5 sigma (depending on the broadband model and method of error covariance matrix estimation) and is consistent with predictions of the standard LCDM model. We assess the biases in our method, stability of the error covariance matrix and possible systematic effects. We fit the resulting correlation function with several models that decouple the broadband and acoustic scale information. For an isotropic dilation factor, we measure 100x(alpha_iso-1) = -1.6 ^{+2.0+4.3+7.4}_{-2.0-4.1-6.8} (stat.) +/- 1.0 (syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits) with respect to the acoustic scale in the fiducial cosmological model (flat LCDM with Omega_m=0.27, h=0.7). When fitting separately for the radial and transversal dilation factors we find marginalised constraints 100x(alpha_par-1) = -1.3 ^{+3.5+7.6 +12.3}_{-3.3-6.7-10.2} (stat.) +/- 2.0 (syst.) and 100x(alpha_perp-1) = -2.2 ^{+7.4+17}_{-7.1-15} +/- 3.0 (syst.). The dilation factor measurements are significantly correlated with cross-correlation coefficient of ~ -0.55. Errors become significantly non-Gaussian for deviations over 3 standard deviations from best fit value. Because of the data cuts and analysis method, these measurements give tighter constraints than a previous BAO analysis of the BOSS DR9 Lyman-alpha forest sample, providing an important consistency test of the standard cosmological model in a new redshift regime.

preprint2013arXiv

Observational Probes of Cosmic Acceleration

The accelerating expansion of the universe is the most surprising cosmological discovery in many decades, implying that the universe is dominated by some form of "dark energy" with exotic physical properties, or that Einstein's theory of gravity breaks down on cosmological scales. The profound implications of cosmic acceleration have inspired ambitious experimental efforts to measure the history of expansion and growth of structure with percent-level precision or higher. We review in detail the four most well established methods for making such measurements: Type Ia supernovae, baryon acoustic oscillations (BAO), weak gravitational lensing, and galaxy clusters. We pay particular attention to the systematic uncertainties in these techniques and to strategies for controlling them at the level needed to exploit "Stage IV" dark energy facilities such as BigBOSS, LSST, Euclid, and WFIRST. We briefly review a number of other approaches including redshift-space distortions, the Alcock-Paczynski test, and direct measurements of H_0. We present extensive forecasts for constraints on the dark energy equation of state and parameterized deviations from GR, achievable with Stage III and Stage IV experimental programs that incorporate supernovae, BAO, weak lensing, and CMB data. We also show the level of precision required for other methods to provide constraints competitive with those of these fiducial programs. We emphasize the value of a balanced program that employs several of the most powerful methods in combination, both to cross-check systematic uncertainties and to take advantage of complementary information. Surveys to probe cosmic acceleration produce data sets with broad applications, and they continue the longstanding astronomical tradition of mapping the universe in ever greater detail over ever larger scales.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological constraints from the full shape of the clustering wedges

We explore the cosmological implications of the clustering wedges, xi_perp(s) and xi_para(s), of the CMASS Data Release 9 (DR9) sample of the Baryon Oscillation Spectroscopic Survey (BOSS). These clustering wedges are defined by averaging the full two-dimensional correlation function, xi(mu,s), over the ranges 0<mu<0.5 and 0.5<mu<1, respectively. These measurements allow us to constrain the parameter combinations D_A(z)/r_s(z_d)=9.03 +- 0.21 and cz/(r_s(z_d)H(z)) = 12.14 +- 0.43 at the mean redsfhit of the sample, z=0.57. We combine the information from the clustering wedges with recent measurements of CMB, BAO and type Ia supernovae to obtain constraints on the cosmological parameters of the standard LCDM model and a number of potential extensions. The information encoded in the clustering wedges is most useful when the dark energy equation of state is allowed to deviate from its standard LCDM value. The combination of all datasets shows no evidence of a deviation from a constant dark energy equation of state, in which case we find w_DE = -1.013 +- 0.064, in complete agreement with a cosmological constant. We explore potential deviations from general relativity by constraining the growth rate f(z)=d ln D(a)/ d ln a, in which case the combination of the CMASS clustering wedges with CMB data implies f(z=0.57)=0.719 +- 0.094, in accordance with the predictions of GR. Our results clearly illustrate the additional constraining power of anisotropic clustering measurements with respect to that of angle-averaged quantities.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Luminosity and Color Dependence and Redshift Evolution

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over ~3,300 sq. deg in the redshift range 0.43<z<0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r_0=7.93\pm0.06Mpc/h and an index of γ=1.85\pm0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Measuring D_A and H at z=0.57 from the Baryon Acoustic Peak in the Data Release 9 Spectroscopic Galaxy Sample

We present measurements of the angular diameter distance to and Hubble parameter at z=0.57 from the measurement of the baryon acoustic peak in the correlation of galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey. Our analysis is based on a sample from Data Release 9 of 264,283 galaxies over 3275 square degrees in the redshift range 0.43<z<0.70. We use two different methods to provide robust measurement of the acoustic peak position across and along the line of sight in order to measure the cosmological distance scale. We find D_A(0.57) = 1408 +/- 45 Mpc and H(0.57) = 92.9 +/- 7.8 km/s/Mpc for our fiducial value of the sound horizon. These results from the anisotropic fitting are fully consistent with the analysis of the spherically averaged acoustic peak position presented in Anderson et al, 2012. Our distance measurements are a close match to the predictions of the standard cosmological model featuring a cosmological constant and zero spatial curvature.

preprint2013arXiv

The Neutral Hydrogen Content of Galaxies in Cosmological Hydrodynamic Simulations

We examine the global HI properties of galaxies in quarter-billion particle cosmological simulations using Gadget-2, focusing on how galactic outflows impact HI content. We consider four outflow models, including a new one (ezw) motivated by recent interstellar medium simulations in which the wind speed and mass loading factor scale as expected for momentum-driven outflows for larger galaxies and energy-driven outflows for dwarfs (sigma<75 km/s). To obtain predicted HI masses, we employ a simple but effective local correction for particle self-shielding, and an observationally-constrained transition from neutral to molecular hydrogen. Our ezw simulation produces an HI mass function whose faint-end slope of -1.3 agrees well with observations from the ALFALFA survey; other models agree less well. Satellite galaxies have a bimodal distribution in HI fraction versus halo mass, with smaller satellites and/or those in larger halos more often being HI-deficient. At a given stellar mass, HI content correlates with star formation rate and inversely correlates with metallicity, as expected if driven by stochasticity in the accretion rate. To higher redshifts, massive HI galaxies disappear and the mass function steepens. The global cosmic HI density conspires to remain fairly constant from z~5-0, but the relative contribution from smaller galaxies increases with redshift.

preprint2013arXiv

The one-dimensional Ly-alpha forest power spectrum from BOSS

We have developed two independent methods to measure the one-dimensional power spectrum of the transmitted flux in the Lyman-$α$ forest. The first method is based on a Fourier transform, and the second on a maximum likelihood estimator. The two methods are independent and have different systematic uncertainties. The determination of the noise level in the data spectra was subject to a novel treatment, because of its significant impact on the derived power spectrum. We applied the two methods to 13,821 quasar spectra from SDSS-III/BOSS DR9 selected from a larger sample of over 60,000 spectra on the basis of their high quality, large signal-to-noise ratio, and good spectral resolution. The power spectra measured using either approach are in good agreement over all twelve redshift bins from $<z> = 2.2$ to $<z> = 4.4$, and scales from 0.001 $\rm(km/s)^{-1}$ to $0.02 \rm(km/s)^{-1}$. We determine the methodological and instrumental systematic uncertainties of our measurements. We provide a preliminary cosmological interpretation of our measurements using available hydrodynamical simulations. The improvement in precision over previously published results from SDSS is a factor 2--3 for constraints on relevant cosmological parameters. For a $Λ$CDM model and using a constraint on $H_0$ that encompasses measurements based on the local distance ladder and on CMB anisotropies, we infer $σ_8 =0.83\pm0.03$ and $n_s= 0.97\pm0.02$ based on \ion{H}{i} absorption in the range $2.1<z<3.7$.

preprint2013arXiv

The Redshift-Space Cluster-Galaxy Cross-Correlation Function: I. Modeling Galaxy Infall onto Millennium Simulation Clusters and SDSS Groups

The large scale infall of galaxies around massive clusters provides a potentially powerful diagnostic of structure growth, dark energy, and cosmological deviations from General Relativity. We develop and test a method to recover galaxy infall kinematics (GIK) from measurements of the redshift-space cluster-galaxy cross-correlation function ξ_{cg}(r_p,r_π). Using galaxy and halo samples from the Millennium simulation, we calibrate an analytic model of the galaxy kinematic profiles comprised of a virialized component with an isotropic Gaussian velocity distribution and an infall component described by a skewed 2D t-distribution with a characteristic infall velocity v_r and separate radial and tangential dispersions. We show that convolving the real-space cross-correlation function with this velocity distribution accurately predicts the redshift-space ξ_{cg}, and we show that measurements of ξ_{cg} can be inverted to recover the four distinct elements of the GIK profiles. These in turn provide diagnostics of cluster mass profiles, and we expect the characteristic infall velocity v_r(r) in particular to be insensitive to galaxy formation physics that can affect velocity dispersions within halos. As a proof of concept we measure ξ_{cg} for rich galaxy groups in the Sloan Digital Sky Survey and recover GIK profiles for groups in two bins of central galaxy stellar mass. The higher mass bin has a v_r(r) curve very similar to that of 10^{14} Msun halos in the Millennium simulation, and the recovered kinematics follow the expected trends with mass. GIK modeling of cluster-galaxy cross-correlations can be a valuable complement to stacked weak lensing analyses, allowing novel tests of modified gravity theories that seek to explain cosmic acceleration.

preprint2013arXiv

Tracing Inflows and Outflows with Absorption Lines in Circumgalactic Gas

We examine how HI and metal absorption lines within low-redshift galaxy halos trace the dynamical state of circumgalactic gas, using cosmological hydrodynamic simulations that include a well-vetted heuristic model for galactic outflows. We categorize inflowing, outflowing, and ambient gas based on its history and fate as tracked in our simulation. Following our earlier work showing that the ionisation level of absorbers was a primary factor in determining the physical conditions of absorbing gas, we show here that it is also a governing factor for its dynamical state. Low-ionisation metal absorbers (e.g. MgII) tend to arise in gas that will fall onto galaxies within several Gyr, while high-ionisation metal absorbers (e.g. OVI) generally trace material that was deposited by outflows many Gyr ago. Inflowing gas is dominated by enriched material that was previously ejected in an outflow, hence accretion at low redshifts is typically substantially enriched. Recycling wind material is preferentially found closer to galaxies, and is more dominant in lower-mass halos since high-mass halos have more hot gas that is able to support itself against infall. Low-mass halos also tend to re-eject more of their accreted material, owing to our outflow prescription that employs higher mass loading factors for lower-mass galaxies. Typical HI absorbers trace unenriched ambient material that is not participating in the baryon cycle, but stronger HI absorbers arise in cool, enriched inflowing gas. Instantaneous radial velocity measures of absorbers are generally poor at distinguishing between inflowing and outflowing gas, except in the case of very recent outflows. These results suggest that probing halo gas using a range of absorbers can provide detailed information about the amount and physical conditions of material that is participating in the baryon cycle.

preprint2012arXiv

A first application of the Alcock-Paczynski test to stacked cosmic voids

We report on the first application of the Alcock-Paczynski test to stacked voids in spectroscopic galaxy redshift surveys.We use voids from the Sutter et al. (2012) void catalog, which was derived from the Sloan Digital Sky Survey Data Release 7 main sample and luminous red galaxy catalogs. The construction of that void catalog removes potential shape measurement bias by using a modified version of the ZOBOV algorithm and by removing voids near survey boundaries and masks. We apply the shape-fitting procedure presented in Lavaux & Wandelt (2012) to ten void stacks out to redshift z=0.36. Combining these measurements, we determine the mean cosmologically induced "stretch" of voids in three redshift bins, with 1-sigma errors of 5-15%. The mean stretch is consistent with unity, providing no indication of a distortion induced by peculiar velocities. While the statistical errors are too large to detect the Alcock-Paczynski effect over our limited redshift range, this proof-of-concept analysis defines procedures that can be applied to larger spectroscopic galaxy surveys at higher redshifts to constrain dark energy using the expected statistical isotropy of structures that are minimally affected by uncertainties in galaxy velocity bias.

preprint2012arXiv

A public void catalog from the SDSS DR7 Galaxy Redshift Surveys based on the watershed transform

We produce the most comprehensive public void catalog to date using the Sloan Digital Sky Survey Data Release 7 main sample out to redshift z=0.2 and the luminous red galaxy sample out to z=0.44. Using a modified version of the parameter-free void finder ZOBOV, we fully take into account the presence of the survey boundary and masks. Our strategy for finding voids is thus appropriate for any survey configuration. We produce two distinct catalogs: a complete catalog including voids near any masks, which would be appropriate for void galaxy surveys, and a bias-free catalog of voids away from any masks, which is necessary for analyses that require a fair sampling of void shapes and alignments. Our discovered voids have effective radii from 5 to 135 h^-1 Mpc. We discuss basic catalog statistics such as number counts and redshift distributions and describe some additional data products derived from our catalog, such as radial density profiles and projected density maps. We find that radial profiles of stacked voids show a qualitatively similar behavior across nearly two decades of void radii and throughout the full redshift range.

preprint2012arXiv

Accretion-Driven Evolution of Black Holes: Eddington Ratios, Duty Cycles, and Active Galaxy Fractions

We develop semi-empirical models of the supermassive black hole and active galactic nucleus (AGN) populations, which incorporate the black hole growth implied by the observed AGN luminosity function assuming a radiative efficiency ε, and a distribution of Eddington ratios λ. By generalizing these continuity-equation models to allow a distribution P(λ|mbh,z) we are able to draw on constraints from observationally estimated P(λ) distributions and active galaxy fractions while accounting for the luminosity thresholds of observational samples. We consider models with a Gaussian distribution of log λ, and Gaussians augmented with a power-law tail to low λ. Within our framework, reproducing the high observed AGN fractions at low redshift requires a characteristic Eddington ratio λ_c that declines at late times, and matching observed Eddington ratio distributions requires P(λ) that broadens at low redshift. To reproduce the observed increase of AGN fraction with black hole or galaxy mass, we also require a λ_c that decreases with increasing black hole mass, reducing the AGN luminosity associated with the most massive black holes. Finally, achieving a good match to the high mass end of the local black hole mass function requires an increased radiative efficiency at high black hole mass. We discuss the potential impact of black hole mergers or a λ-dependent bolometric correction, and we compute evolutionary predictions for black hole and galaxy specific accretion rates. Despite the flexibility of our framework, no one model provides a good fit to all the data we consider.

preprint2012arXiv

Cosmological Constraints from the Large Scale Weak Lensing of SDSS MaxBCG Clusters

We derive constraints on the matter density \Om and the amplitude of matter clustering \sig8 from measurements of large scale weak lensing (projected separation R=5-30\hmpc) by clusters in the Sloan Digital Sky Survey MaxBCG catalog. The weak lensing signal is proportional to the product of \Om and the cluster-mass correlation function \xicm. With the relation between optical richness and cluster mass constrained by the observed cluster number counts, the predicted lensing signal increases with increasing \Om or \sig8, with mild additional dependence on the assumed scatter between richness and mass. The dependence of the signal on scale and richness partly breaks the degeneracies among these parameters. We incorporate external priors on the richness-mass scatter from comparisons to X-ray data and on the shape of the matter power spectrum from galaxy clustering, and we test our adopted model for \xicm against N-body simulations. Using a Bayesian approach with minimal restrictive priors, we find \sig8(\Om/0.325)^{0.501}=0.828 +/- 0.049, with marginalized constraints of \Om=0.325_{-0.067}^{+0.086} and \sig8=0.828_{-0.097}^{+0.111}, consistent with constraints from other MaxBCG studies that use weak lensing measurements on small scales (R<=2\hmpc). The (\Om,\sig8) constraint is consistent with and orthogonal to the one inferred from WMAP CMB data, reflecting agreement with the structure growth predicted by GR for an LCDM cosmological model. A joint constraint assuming LCDM yields \Om=0.298 +/- 0.020 and \sig8=0.831 +/- 0.020. Our cosmological parameter errors are dominated by the statistical uncertainties of the large scale weak lensing measurements, which should shrink sharply with current and future imaging surveys.

preprint2012arXiv

Not Dead Yet: Cool Circumgalactic Gas in the Halos of Early Type Galaxies

We report new observations of circumgalactic gas in the halos of early type galaxies obtained by the COS-Halos Survey with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. We find that detections of HI surrounding early type galaxies are typically as common and strong as around star-forming galaxies, implying that the total mass of circumgalactic material is comparable in the two populations. For early type galaxies, the covering fraction for HI absorption above 10^16 cm^2 is ~40-50% within ~150 kpc. Line widths and kinematics of the detected material show it to be cold (T ~< 10^5 K) in comparison to the virial temperature of the host halos. The implied masses of cool, photoionized CGM baryons may be up to 10^9 --- 10^11 Msun. Contrary to some theoretical expectations, strong halo HI absorbers do not disappear as part of the quenching of star-formation. Even passive galaxies retain significant reservoirs of halo baryons which could replenish the interstellar gas reservoir and eventually form stars. This halo gas may feed the diffuse and molecular gas that is frequently observed inside ETGs.

preprint2012arXiv

Principal Component Abundance Analysis of Microlensed Bulge Dwarf and Subgiant Stars

Elemental abundance patterns can provide vital clues to the formation and enrichment history of a stellar population. Here we present an investigation of the Galactic bulge, where we apply principal component abundance analysis (PCAA)---a principal component decomposition of relative abundances [X/Fe]---to a sample of 35 microlensed bulge dwarf and subgiant stars, characterizing their distribution in the 12-dimensional space defined by their measured elemental abundances. The first principal component PC1, which suffices to describe the abundance patterns of most stars in the sample, shows a strong contribution from alpha-elements, reflecting the relative contributions of Type II and Type Ia supernovae. The second principal component PC2 is characterized by a Na--Ni correlation, the likely product of metallicity-dependent Type II supernova yields. The distribution in PC1 is bimodal, showing that the bimodality previously found in the [Fe/H] values of these stars is robustly and independently recovered by looking at only their relative abundance patterns. The two metal-rich stars that are alpha-enhanced have outlier values of PC2 and PC3, respectively, further evidence that they have distinctive enrichment histories. Applying PCAA to a sample of local thin and thick disk dwarfs yields a nearly identical PC1; in PC1, the metal-rich and metal-poor bulge dwarfs track kinematically selected thin and thick disk dwarfs, respectively, suggesting broadly similar alpha-enrichment histories. However, the disk PC2 is dominated by a Y--Ba correlation, likely indicating a greater contribution of s-process enrichment from long-lived asymptotic giant branch stars that is absent from the bulge PC2 because of its rapid formation.

preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

preprint2012arXiv

The BOSS Lyman-alpha Forest Sample from SDSS Data Release 9

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

preprint2012arXiv

The Center is Everywhere

"The Center is Everywhere" is a sculpture by Josiah McElheny, currently (through October 14, 2012) on exhibit at the Institute of Contemporary Art, Boston. The sculpture is based on data from the Sloan Digital Sky Survey (SDSS), using hundreds of glass crystals and lamps suspended from brass rods to represent the three-dimensional structure mapped by the SDSS through one of its 2000+ spectroscopic plugplates. This article describes the scientific ideas behind this sculpture, emphasizing the principle of the statistical homogeneity of cosmic structure in the presence of local complexity. The title of the sculpture is inspired by the work of the French revolutionary Louis Auguste Blanqui, whose 1872 book "Eternity Through The Stars: An Astronomical Hypothesis" was the first to raise the spectre of the infinite replicas expected in an infinite, statistically homogeneous universe. Puzzles of infinities, probabilities, and replicas continue to haunt modern fiction and contemporary discussions of inflationary cosmology.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 9 Spectroscopic Galaxy Sample

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5σin both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7σ. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the large-scale two-point correlation function

We obtain constraints on cosmological parameters from the spherically averaged redshift-space correlation function of the CMASS Data Release 9 (DR9) sample of the Baryonic Oscillation Spectroscopic Survey (BOSS). We combine this information with additional data from recent CMB, SN and BAO measurements. Our results show no significant evidence of deviations from the standard flat-Lambda CDM model, whose basic parameters can be specified by Omega_m = 0.285 +- 0.009, 100 Omega_b = 4.59 +- 0.09, n_s = 0.96 +- 0.009, H_0 = 69.4 +- 0.8 km/s/Mpc and sigma_8 = 0.80 +- 0.02. The CMB+CMASS combination sets tight constraints on the curvature of the Universe, with Omega_k = -0.0043 +- 0.0049, and the tensor-to-scalar amplitude ratio, for which we find r < 0.16 at the 95 per cent confidence level (CL). These data show a clear signature of a deviation from scale-invariance also in the presence of tensor modes, with n_s <1 at the 99.7 per cent CL. We derive constraints on the fraction of massive neutrinos of f_nu < 0.049 (95 per cent CL), implying a limit of sum m_nu < 0.51 eV. We find no signature of a deviation from a cosmological constant from the combination of all datasets, with a constraint of w_DE = -1.033 +- 0.073 when this parameter is assumed time-independent, and no evidence of a departure from this value when it is allowed to evolve as w_DE(a) = w_0 + w_a (1 - a). The achieved accuracy on our cosmological constraints is a clear demonstration of the constraining power of current cosmological observations.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: measurements of the growth of structure and expansion rate at z=0.57 from anisotropic clustering

We analyze the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which consists of 264,283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3,275 square degrees. Both peculiar velocities and errors in the assumed redshift-distance relation ("Alcock-Paczynski effect") generate correlations between clustering amplitude and orientation with respect to the line-of-sight. Together with the sharp baryon acoustic oscillation (BAO) standard ruler, our measurements of the broadband shape of the monopole and quadrupole correlation functions simultaneously constrain the comoving angular diameter distance (2190 +/- 61 Mpc) to z=0.57, the Hubble expansion rate at z=0.57 (92.4 +/- 4.5 km/s/Mpc), and the growth rate of structure at that same redshift (d sigma8/d ln a = 0.43 +/- 0.069). Our analysis provides the best current direct determination of both DA and H in galaxy clustering data using this technique. If we further assume a LCDM expansion history, our growth constraint tightens to d sigma8/d ln a = 0.415 +/- 0.034. In combination with the cosmic microwave background, our measurements of DA, H, and growth all separately require dark energy at z > 0.57, and when combined imply Ω_Λ = 0.74 +/- 0.016, independent of the Universe's evolution at z<0.57. In our companion paper (Samushia et al. prep), we explore further cosmological implications of these observations.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: the low redshift sample

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.

preprint2012arXiv

The clustering of intermediate redshift quasars as measured by the Baryon Oscillation Spectroscopic Survey

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

preprint2012arXiv

The Milky Way's circular velocity curve between 4 and 14 kpc from APOGEE data

We measure the Milky Way's rotation curve over the Galactocentric range 4 kpc <~ R <~ 14 kpc from the first year of data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We model the line-of-sight velocities of 3,365 stars in fourteen fields with b = 0 deg between 30 deg < l < 210 deg out to distances of 10 kpc using an axisymmetric kinematical model that includes a correction for the asymmetric drift of the warm tracer population (σ_R ~ 35 km/s). We determine the local value of the circular velocity to be V_c(R_0) = 218 +/- 6 km/s and find that the rotation curve is approximately flat with a local derivative between -3.0 km/s/kpc and 0.4 km/s/kpc. We also measure the Sun's position and velocity in the Galactocentric rest frame, finding the distance to the Galactic center to be 8 kpc < R_0 < 9 kpc, radial velocity V_{R,sun} = -10 +/- 1 km/s, and rotational velocity V_{ϕ,sun} = 242^{+10}_{-3} km/s, in good agreement with local measurements of the Sun's radial velocity and with the observed proper motion of Sgr A*. We investigate various systematic uncertainties and find that these are limited to offsets at the percent level, ~2 km/s in V_c. Marginalizing over all the systematics that we consider, we find that V_c(R_0) < 235 km/s at >99% confidence. We find an offset between the Sun's rotational velocity and the local circular velocity of 26 +/- 3 km/s, which is larger than the locally-measured solar motion of 12 km/s. This larger offset reconciles our value for V_c with recent claims that V_c >~ 240 km/s. Combining our results with other data, we find that the Milky Way's dark-halo mass within the virial radius is ~8x10^{11} M_sun.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: The Quasar Luminosity Function from Data Release Nine

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

preprint2012arXiv

The Sloan Digital Sky Survey quasar catalog: ninth data release

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $Λ$CDM cosmology with H0 = 70 km/s/Mpc, $Ω_{\rm M}$ = 0.3, and $Ω_Λ$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

preprint2011arXiv

Cosmological Constraints from Galaxy Clustering and the Mass-to-Number Ratio of Galaxy Clusters

We place constraints on the average density (Omega_m) and clustering amplitude (sigma_8) of matter using a combination of two measurements from the Sloan Digital Sky Survey: the galaxy two-point correlation function, w_p, and the mass-to-galaxy-number ratio within galaxy clusters, M/N, analogous to cluster M/L ratios. Our w_p measurements are obtained from DR7 while the sample of clusters is the maxBCG sample, with cluster masses derived from weak gravitational lensing. We construct non-linear galaxy bias models using the Halo Occupation Distribution (HOD) to fit both w_p and M/N for different cosmological parameters. HOD models that match the same two-point clustering predict different numbers of galaxies in massive halos when Omega_m or sigma_8 is varied, thereby breaking the degeneracy between cosmology and bias. We demonstrate that this technique yields constraints that are consistent and competitive with current results from cluster abundance studies, even though this technique does not use abundance information. Using w_p and M/N alone, we find Omega_m^0.5*sigma_8=0.465+/-0.026, with individual constraints of Omega_m=0.29+/-0.03 and sigma_8=0.85+/-0.06. Combined with current CMB data, these constraints are Omega_m=0.290+/-0.016 and sigma_8=0.826+/-0.020. All errors are 1-sigma. The systematic uncertainties that the M/N technique are most sensitive to are the amplitude of the bias function of dark matter halos and the possibility of redshift evolution between the SDSS Main sample and the maxBCG sample. Our derived constraints are insensitive to the current level of uncertainties in the halo mass function and in the mass-richness relation of clusters and its scatter, making the M/N technique complementary to cluster abundances as a method for constraining cosmology with future galaxy surveys.

preprint2011arXiv

Galaxy Clustering in the Completed SDSS Redshift Survey: The Dependence on Color and Luminosity

We measure the luminosity and color dependence of galaxy clustering in the SDSS DR7 main galaxy sample, focusing on the projected correlation function w_p(r_p) of volume-limited samples. We interpret our measurements using halo occupation distribution (HOD) modeling assuming a Lambda-CDM cosmology. The amplitude of w_p(r_p) grows slowly with luminosity for L < L_* and increases sharply at higher luminosities, with bias factor b(>L)=1.06+0.23(L/L_*)^{1.12}. At fixed luminosity, redder galaxies have a stronger and steeper w_p(r_p), a trend that runs steadily from the bluest galaxies to the reddest galaxies. The individual luminosity trends for the red and blue galaxy populations are strikingly different. Blue galaxies show a slow but steady increase of w_p(r_p) with luminosity, at all scales. The large-scale clustering of red galaxies shows little luminosity dependence until a sharp increase at L > 4L_*, but the lowest luminosity red galaxies (0.04-0.25 L_*) show very strong clustering on scales r_p < 2 Mpc/h. Most of the observed trends can be naturally understood within the LCDM+HOD framework. The growth of w_p(r_p) with luminosity reflects an overall shift in the halo mass scale, in particular an increase in the minimum host halo mass Mmin. The mass at which a halo has, on average, one satellite galaxy brighter than L is M_1 ~ 17 Mmin(L) over most of the luminosity range. The growth and steepening of w_p(r_p) for redder galaxies reflects the increasing fraction of galaxies that are satellite systems in high mass halos instead of central systems in low mass halos, a trend that is especially marked at low luminosities. Our extensive measurements, provided in tabular form, will allow detailed tests of theoretical models of galaxy formation, a firm grounding of semi-empirical models of the galaxy population, and new cosmological tests.

preprint2011arXiv

Intergalactic Dust Extinction in Hydrodynamic Cosmological Simulations

Recently Menard et al. detected a subtle but systematic change in the mean color of quasars as a function of their projected separation from foreground galaxies, extending to comoving separations of ~10Mpc/h, which they interpret as a signature of reddening by intergalactic dust. We present theoretical models of this remarkable observation, using SPH cosmological simulations of a (50Mpc/h)^3 volume. Our primary model uses a simulation with galactic winds and assumes that dust traces the intergalactic metals. The predicted galaxy-dust correlation function is similar in form to the galaxy-mass correlation function, and reproducing the MSFR data requires a dust-to-metal mass ratio of 0.24, about half the value in the Galactic ISM. Roughly half of the reddening arises in dust that is more than 100Kpc/h from the nearest massive galaxy. We also examine a simulation with no galactic winds, which predicts a much smaller fraction of intergalactic metals (3% vs. 35%) and therefore requires an unphysical dust-to-metal ratio of 2.18 to reproduce the MSFR data. In both models, the signal is dominated by sightlines with E(g-i)=0.001-0.1. The no-wind simulation can be reconciled with the data if we also allow reddening to arise in galaxies up to several x 10^10 Msun. The wind model predicts a mean visual extinction of A_V ~0.0133 mag out to z=0.5, with a sightline-to-sightline dispersion similar to the mean, which could be significant for future supernova cosmology studies. Reproducing the MSFR results in these simulations requires that a large fraction of ISM dust survive its expulsion from galaxies and its residence in the intergalactic medium. Future observational studies that provide higher precision and measure the dependence on galaxy type and environment will allow detailed tests for models of enriched galactic outflows and the survival of IG dust.

preprint2011arXiv

Radial Mixing in Galactic Disks: The Effects of Disk Structure and Satellite Bombardment

We use a suite of numerical simulations to investigate the mechanisms and effects of radial migration of stars in disk galaxies like the Milky Way (MW). An isolated, collisionless stellar disk with a MW-like scale-height shows only the radial "blurring" expected from epicyclic orbits. Reducing the disk thickness or adding gas to the disk substantially increases the level of radial migration, induced by interaction with transient spiral arms and/or a central bar. We also examine collisionless disks subjected to gravitational perturbations from a cosmologically motivated satellite accretion history. In the perturbed disk that best reproduces the observed properties of the MW, 20% of stars that end up in the solar annulus 7 kpc < R < 9 kpc started at R < 6 kpc, and 7% started at R > 10 kpc. This level of migration would add considerable dispersion to the age-metallicity relation of solar neighborhood stars. In the isolated disk models, the probability of migration traces the disk's radial mass profile, but in perturbed disks migration occurs preferentially at large radii, where the disk is more weakly bound. The orbital dynamics of migrating particles are also different in isolated and perturbed disks: satellite perturbations drive particles to lower angular momentum for a given change in radius. Thus, satellite perturbations appear to be a distinct mechanism for inducing radial migration, which can operate in concert with migration induced by bars and spiral structure. We investigate correlations between changes in radius and changes in orbital circularity or vertical energy, identifying signatures that might be used to test models and distinguish radial migration mechanisms in chemo-dynamical surveys of the MW disk.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

Self-Similar Bumps and Wiggles: Isolating the Evolution of the BAO Peak with Power-law Initial Conditions

Motivated by cosmological surveys that demand accurate theoretical modeling of the baryon acoustic oscillation (BAO) feature in galaxy clustering, we analyze N-body simulations in which a BAO-like gaussian bump modulates the linear theory correlation function ξ_L(r)=(r_0/r)^{n+3} of an underlying self-similar model with initial power spectrum P(k)=A k^n. These simulations test physical and analytic descriptions of BAO evolution far beyond the range of most studies, since we consider a range of underlying power spectra (n=-0.5, -1, -1.5) and evolve simulations to large effective correlation amplitudes (equivalent to σ_8=4-12 for r_bao = 100 Mpc/h). In all cases, non-linear evolution flattens and broadens the BAO bump in ξ(r) while approximately preserving its area. This evolution resembles a "diffusion" process in which the bump width σ_bao is the quadrature sum of the linear theory width and a length proportional to the rms relative displacement Σ_pair(r_bao}) of particle pairs separated by r_bao. For n=-0.5 and n=-1, we find no detectable shift of the location of the BAO peak, but the peak in the n=-1.5 model shifts steadily to smaller scales, following r_peak}/r_bao = 1-1.08(r_0/r_bao)^{1.5}. The "SimpleRG" perturbation theory scheme and, to a lesser extent, standard 1-loop perturbation theory are fairly successful at explaining the non-linear evolution of the fourier power spectrum of our models. Analytic models also explain why the ξ(r) peak shifts much more for n=-1.5 than for n >= -1, though no ab initio model we have examined reproduces all of our numerical results. Simulations with L_box = 10 r_bao} and L_box = 20 r_bao yield consistent results for ξ(r) at the BAO scale, provided one corrects for the integral constraint imposed by the uniform density box.

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2011arXiv

The intergalactic medium over the last 10 billion years II: Metal-line absorption and physical conditions

We investigate the metallicity evolution and content of the intergalactic medium (IGM) and galactic halo gas from z=2->0 using 110-million particle cosmological hydrodynamic simulations. We focus on the detectability and physical properties of UV resonance metal-line absorbers observable with Hubble's Cosmic Origins Spectrograph (COS). We confirm that galactic superwind outflows are required to enrich the IGM to observed levels down to z=0 using three wind prescriptions contrasted to a no-wind simulation. Our favoured momentum-conserved wind prescription deposits metals closer to galaxies owing to its moderate energy input, while the more energetic constant wind model enriches the warm-hot IGM 6.4x more. Despite these significant differences, all wind models produce metal-line statistics within a factor of two of existing observations. This is because OVI, CIV, SiIV, and NeVIII absorbers primarily arise from T<10^5 K, photo-ionised gas that is enriched to similar levels in the three feedback schemes. OVI absorbers trace the diffuse phase with delta<100, which is enriched to ~1/50 Zsol at z=0, although the absorbers themselves usually exceed 0.3 Zsol and arise from inhomogeneously distributed, un-mixed winds. CIV and SiIV absorbers trace primarily T~10^4 K gas inside haloes. We predict COS will observe a population of NeVIII photo-ionised absorbers tracing T<10^5 K, delta~10 gas. MgX and SiXII are rarely detected in COS S/N=30 simulated sight lines although detections trace T=10^(6-7) K halo gas. In general, the IGM is enriched in an outside-in manner, where wind-blown metals released at higher redshift reach lower overdensities, resulting in higher ionisation species tracing lower-density, older metals. At z=0, the 90% of baryons outside of galaxies are enriched to 0.096 Zsol, but the 65% of unbound baryons in the IGM have 0.018 Zsol and contain only 4% of all metals.

preprint2011arXiv

The Large, Oxygen-Rich Halos of Star-Forming Galaxies Are A Major Reservoir of Galactic Metals

The circumgalactic medium (CGM) is fed by galaxy outflows and accretion of intergalactic gas, but its mass, heavy element enrichment, and relation to galaxy properties are poorly constrained by observations. In a survey of the outskirts of 42 galaxies with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope, we detected ubiquitous, large (150 kiloparsec) halos of ionized oxygen surrounding star-forming galaxies, but we find much less ionized oxygen around galaxies with little or no star formation. This ionized CGM contains a substantial mass of heavy elements and gas, perhaps far exceeding the reservoirs of gas in the galaxies themselves. It is a basic component of nearly all star-forming galaxies that is removed or transformed during the quenching of star formation and the transition to passive evolution.

preprint2011arXiv

The Lyman-alpha forest in three dimensions: measurements of large scale flux correlations from BOSS 1st-year data

Using a sample of approximately 14,000 z>2.1 quasars observed in the first year of the Baryon Oscillation Spectroscopic Survey (BOSS), we measure the three-dimensional correlation function of absorption in the Lyman-alpha forest. The angle-averaged correlation function of transmitted flux (F = exp(-tau)) is securely detected out to comoving separations of 60 Mpc/h, the first detection of flux correlations across widely separated sightlines. A quadrupole distortion of the redshift-space correlation function by peculiar velocities, the signature of the gravitational instability origin of structure in the Lyman-alpha forest, is also detected at high significance. We obtain a good fit to the data assuming linear theory redshift-space distortion and linear bias of the transmitted flux, relative to the matter fluctuations of a standard LCDM cosmological model (inflationary cold dark matter with a cosmological constant). At 95% confidence, we find a linear bias parameter 0.16<b<0.24 and redshift-distortion parameter 0.44<beta<1.20, at central redshift z=2.25, with a well constrained combination b(1+β)=0.336 +/- 0.012. The errors on beta are asymmetric, with beta=0 excluded at over 5 sigma confidence level. The value of beta is somewhat low compared to theoretical predictions, and our tests on synthetic data suggest that it is depressed (relative to expectations for the Lyman-alpha forest alone) by the presence of high column density systems and metal line absorption. These results set the stage for cosmological parameter determinations from three-dimensional structure in the Lyman-alpha forest, including anticipated constraints on dark energy from baryon acoustic oscillations.

preprint2011arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: Quasar Target Selection for Data Release Nine

The SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), a five-year spectroscopic survey of 10,000 deg^2, achieved first light in late 2009. One of the key goals of BOSS is to measure the signature of baryon acoustic oscillations in the distribution of Ly-alpha absorption from the spectra of a sample of ~150,000 z>2.2 quasars. Along with measuring the angular diameter distance at z\approx2.5, BOSS will provide the first direct measurement of the expansion rate of the Universe at z > 2. One of the biggest challenges in achieving this goal is an efficient target selection algorithm for quasars over 2.2 < z < 3.5, where their colors overlap those of stars. During the first year of the BOSS survey, quasar target selection methods were developed and tested to meet the requirement of delivering at least 15 quasars deg^-2 in this redshift range, out of 40 targets deg^-2. To achieve these surface densities, the magnitude limit of the quasar targets was set at g <= 22.0 or r<=21.85. While detection of the BAO signature in the Ly-alpha absorption in quasar spectra does not require a uniform target selection, many other astrophysical studies do. We therefore defined a uniformly-selected subsample of 20 targets deg^-2, for which the selection efficiency is just over 50%. This "CORE" subsample will be fixed for Years Two through Five of the survey. In this paper we describe the evolution and implementation of the BOSS quasar target selection algorithms during the first two years of BOSS operations. We analyze the spectra obtained during the first year. 11,263 new z>2.2 quasars were spectroscopically confirmed by BOSS. Our current algorithms select an average of 15 z > 2.2 quasars deg^-2 from 40 targets deg^-2 using single-epoch SDSS imaging. Multi-epoch optical data and data at other wavelengths can further improve the efficiency and completeness of BOSS quasar target selection. [Abridged]

preprint2010arXiv

A Search for Oxygen in the Low-Density Lyman-alpha Forest Using the Sloan Digital Sky Survey

We use 2167 Sloan Digital Sky Survey (SDSS) quasar spectra to search for low-density oxygen in the Intergalactic Medium (IGM). Oxygen absorption is detected on a pixel-by-pixel basis by its correlation with Lyman-alpha forest absorption. We have developed a novel Locally Calibrated Pixel (LCP) search method that uses adjacent regions of the spectrum to calibrate interlopers and spectral artifacts, which would otherwise limit the measurement of OVI absorption. Despite the challenges presented by searching for weak OVI within the Lyman-alpha forest in spectra of moderate resolution and signal-to-noise, we find a highly significant detection of absorption by oxygen at 2.7 < z < 3.2 (the null hypothesis has a chi^2=80 for 9 data points). We interpret our results using synthetic spectra generated from a lognormal density field assuming a mixed quasar-galaxy photoionizing background (Haardt & Madau 2001) and that it dominates the ionization fraction of detected OVI. The LCP search data can be fit by a constant metallicity model with [O/H] = -2.15_(-0.09)^(+0.07), but also by models in which low-density regions are unenriched and higher density regions have a higher metallicity. The density-dependent enrichment model by Aguirre et al. (2008) is also an acceptable fit. All our successful models have similar mass-weighted oxygen abundance, corresponding to [<O/H>_MW] = -2.45+-0.06. This result can be used to find the cosmic oxygen density in the Lyman-alpha forest, Omega_(Oxy, IGM) = 1.4(+-0.2)x10^(-6) = 3x10^(-4) Omega_b. This is the tightest constraint on the mass-weighted mean oxygen abundance and the cosmic oxygen density in the Lyman-alpha forest to date and indicates that it contains approximately 16% of metals produced by star formation (Bouché et al. 2008) up to z = 3.

preprint2010arXiv

Constraints on black hole duty cycles and the black hole-halo relation from SDSS quasar clustering

We use Shen et al.'s (2009) measurements of luminosity-dependent clustering in the SDSS Data Release 5 Quasar Catalog, at redshifts 0.4 < z < 2.5, to constrain the relation between quasar luminosity and host halo mass and to infer the duty cycle f_opt, the fraction of black holes that shine as optically luminous quasars at a given time. We assume a monotonic mean relation between quasar luminosity and host halo mass, with log-normal scatter Σ. For specified f_opt and Σ, matching the observed quasar space density determines the normalization of the luminosity-halo mass relation, from which we predict the clustering bias. The data show no change of bias between the faint and bright halves of the quasar sample but a modest increase in bias for the brightest 10%. At the mean redshift z=1.45 of the sample, the data can be well described either by models with small intrinsic scatter (Σ=0.1 dex) and a duty cycle f_opt=6*10^(-4) or by models with much larger duty cycles and larger values of the scatter. "Continuity equation" models of the black hole mass population imply f_opt > 2*10^(-3) in this range of masses and redshifts, and the combination of this constraint with the clustering measurements implies scatter Σ> 0.4 dex. These findings contrast with those inferred from the much stronger clustering of high-luminosity quasars at z~4, which require minimal scatter between luminosity and halo mass and duty cycles close to one.

preprint2010arXiv

Feedback and Recycled Wind Accretion: Assembling the z=0 Galaxy Mass Function

We analyse cosmological hydrodynamic simulations that include observationally-constrained prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z<1. This recycled wind accretion, or wind mode, provides a third physically distinct accretion channel in addition to the "hot" and "cold" modes emphasised in recent theoretical studies. Because of the interaction between outflows and gas in and around halos, the recycling timescale of wind material (t_rec) is shorter in higher-mass systems, which reside in denser gaseous environments. In these simulations, this differential recycling plays a central role in shaping the present-day galaxy stellar mass function (GSMF). If we remove all particles that were ever ejected in a wind, then the predicted GSMFs are much steeper than observed; galaxy masses are suppressed both by the direct removal of gas and by the hydrodynamic heating of their surroundings, which reduces subsequent infall. With wind recycling included, the simulation that incorporates our favoured momentum-driven wind scalings reproduces the observed GSMF for stellar masses 10^9 < M < 5x10^10 Msolar. At higher masses, wind recycling leads to excessive galaxy masses and excessive star formation rates relative to observations. In these massive systems, some quenching mechanism must suppress the re-accretion of gas ejected from star-forming galaxies. In short, as has long been anticipated, the form of the GSMF is governed by outflows; the unexpected twist here for our simulated winds is that it is not primarily the ejection of material but how the ejected material is re-accreted that governs the GSMF.

preprint2010arXiv

Galaxy Clustering Topology in the Sloan Digital Sky Survey Main Galaxy Sample: a Test for Galaxy Formation Models

We measure the topology of the main galaxy distribution using the Seventh Data Release of the Sloan Digital Sky Survey, examining the dependence of galaxy clustering topology on galaxy properties. The observational results are used to test galaxy formation models. A volume-limited sample defined by $M_r<-20.19$ enables us to measure the genus curve with amplitude of $G=378$ at $6h^{-1}$Mpc smoothing scale, with 4.8\% uncertainty including all systematics and cosmic variance. The clustering topology over the smoothing length interval from 6 to $10 h^{-1}$Mpc reveals a mild scale-dependence for the shift ($Δν$) and void abundance ($A_V$) parameters of the genus curve. We find substantial bias in the topology of galaxy clustering with respect to the predicted topology of the matter distribution, which varies with luminosity, morphology, color, and the smoothing scale of the density field. The distribution of relatively brighter galaxies shows a greater prevalence of isolated clusters and more percolated voids. Even though early (late)-type galaxies show topology similar to that of red (blue) galaxies, the morphology dependence of topology is not identical to the color dependence. In particular, the void abundance parameter $A_V$ depends on morphology more strongly than on color. We test five galaxy assignment schemes applied to cosmological N-body simulations of a $Λ$CDM universe to generate mock galaxies: the Halo-Galaxy one-to-one Correspondence model, the Halo Occupation Distribution model, and three implementations of Semi-Analytic Models (SAMs). None of the models reproduces all aspects of the observed clustering topology; the deviations vary from one model to another but include statistically significant discrepancies in the abundance of isolated voids or isolated clusters and the amplitude and overall shift of the genus curve. (Abridged)

preprint2010arXiv

MARVELS-1b: A Short-Period, Brown Dwarf Desert Candidate from the SDSS-III MARVELS Planet Search

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

preprint2010arXiv

On the radiative efficiencies, Eddington ratios, and duty cycles of luminous high-redshift quasars

We investigate the characteristic radiative efficiency ε, Eddington ratio λ, and duty cycle P_0 of high-redshift active galactic nuclei (AGN), drawing on measurements of the AGN luminosity function at z=3-6 and, especially, on recent measurements of quasar clustering at z=3-4.5 from the Sloan Digital Sky Survey. The free parameters of our models are ε, λ, and the normalization, scatter, and redshift evolution of the relation between black hole mass \mbh and halo virial velocity V_vir. We compute the luminosity function from the implied growth of the black hole mass function and the quasar correlation length from the bias of the host halos. We test our adopted formulae for the halo mass function and halo bias against measurements from the large N-body simulation developed by the MICE collaboration. The strong clustering of AGNs observed at z=3 and, especially, at z=4 implies that massive black holes reside in rare, massive dark matter halos. Reproducing the observed luminosity function then requires high efficiency εand/or low Eddington ratio λ, with a lower limit (based on 2σagreement with the measured z=4 correlation length) ε> 0.7λ/(1+0.7λ), implying ε> 0.17 for λ> 0.25. Successful models predict high duty cycles, P_0~0.2, 0.5, and 0.9 at z=3.1, 4.5 and 6, respectively, and they require that the fraction of halo baryons locked in the central black hole is much larger than the locally observed value. The rapid drop in the abundance of the massive and rare host halos at z>7 implies a proportionally rapid decline in the number density of luminous quasars, much stronger than simple extrapolations of the z=3-6 luminosity function would predict. (abridged)

preprint2010arXiv

Pressure Support vs. Thermal Broadening in the Lyman-alpha Forest I: Effects of the Equation of State on Longitudinal Structure

In the low density intergalactic medium (IGM) that gives rise to the Lyman-alpha forest, gas temperature and density are tightly correlated. The velocity scale of thermal broadening and the Hubble flow across the gas Jeans scale are of similar magnitude (Hlambda_J ~ sigma_th). To separate the effects of gas pressure support and thermal broadening on the Lya forest, we compare spectra extracted from two smoothed particle hydrodynamics (SPH) simulations evolved with different photoionization heating rates (and thus different Jeans scales), imposing different temperature-density relations on the evolved particle distributions. The turnover scales in the flux power spectrum and flux autocorrelation function are determined mainly by thermal broadening rather than pressure. However, the insensitivity to pressure arises partly from a cancellation effect with a sloped temperature-density relation (T ~ rho^{0.6} in our simulations): the high density peaks in the colder, lower pressure simulation are less smoothed by pressure support than in the hotter simulation, and it is this higher density gas that experiences the strongest thermal broadening. Changes in thermal broadening and pressure support have comparably important effects on the flux probability distribution (PDF), which responds directly to the gas overdensity distribution rather than the scale on which it is smooth. Tests on a lower resolution simulation show that our statistical results are converged even at this lower resolution. While thermal broadening generally dominates the longitudinal structure in the Lya forest, we show in Paper II that pressure support determines the transverse coherence of the forest observed towards close quasar pairs. [ABRIDGED]

preprint2010arXiv

Pressure Support vs. Thermal Broadening in the Lyman-alpha Forest II: Effects of the Equation of State on Transverse Structure

We examine the impact of gas pressure on the transverse coherence of high-redshift (2 <= z <= 4) Lyman-alpha forest absorption along neighboring lines of sight that probe the gas Jeans scale (projected separation Delta r <= 500 kpc/h comoving; angular separation Delta theta <= 30"). We compare predictions from two smoothed particle hydrodynamics (SPH) simulations that have different photoionization heating rates and thus different temperature-density relations in the intergalactic medium (IGM). We also compare spectra computed from the gas distributions to those computed from the pressureless dark matter. The coherence along neighboring sightlines is markedly higher for the hotter, higher pressure simulation, and lower for the dark matter spectra. We quantify this coherence using the flux cross-correlation function and the conditional distribution of flux decrements as a function of transverse and line-of-sight (velocity) separation. Sightlines separated by Delta theta <= 15" are ideal for probing this transverse coherence. Higher pressure decreases the redshift-space anisotropy of the flux correlation function, while higher thermal broadening increases the anisotropy. In contrast to the longitudinal (line-of-sight) structure of the Lya forest, the transverse structure on these scales is dominated by pressure effects rather than thermal broadening. With the rapid recent growth in the number of known close quasar pairs, paired line-of-sight observations offer a promising new route to probe the IGM temperature-density relation and test the unexpectedly high temperatures that have been inferred from single sightline analyses.

preprint2010arXiv

The clustering of massive galaxies at z~0.5 from the first semester of BOSS data

We calculate the real- and redshift-space clustering of massive galaxies at z~0.5 using the first semester of data by the Baryon Oscillation Spectroscopic Survey (BOSS). We study the correlation functions of a sample of 44,000 massive galaxies in the redshift range 0.4<z<0.7. We present a halo-occupation distribution modeling of the clustering results and discuss the implications for the manner in which massive galaxies at z~0.5 occupy dark matter halos. The majority of our galaxies are central galaxies living in halos of mass 10^{13}Msun/h, but 10% are satellites living in halos 10 times more massive. These results are broadly in agreement with earlier investigations of massive galaxies at z~0.5. The inferred large-scale bias (b~2) and relatively high number density (nbar=3e-4 h^3 Mpc^{-3}) imply that BOSS galaxies are excellent tracers of large-scale structure, suggesting BOSS will enable a wide range of investigations on the distance scale, the growth of large-scale structure, massive galaxy evolution and other topics.

preprint2010arXiv

The Composite Spectrum of Strong Lyman-alpha Forest Absorbers

We present a new method for probing the physical conditions and metal enrichment of the Intergalactic Medium: the composite spectrum of Ly-alpha forest absorbers. We apply this technique to a sample of 9480 Ly-alpha absorbers with redshift 2 < z < 3.5 identified in the spectra of 13,279 high-redshift quasars from the Sloan Digital Sky Survey (SDSS) Fifth Data Release (DR5). Absorbers are selected as local minima in the spectra with 2.4 < tau_Ly-alpha < 4.0; at SDSS resolution (~ 150km/s FWHM), these absorbers are blends of systems that are individually weaker. In the stacked spectra we detect seven Lyman-series lines and metal lines of O VI, N V, C IV, C III, Si IV, C II, Al II, Si II, Fe II, Mg II, and O I. Many of these lines have peak optical depths of < 0.02, but they are nonetheless detected at high statistical significance. Modeling the Lyman-series measurements implies that our selected systems have total H I column densities N_HI ~ 10^15.4cm-2. Assuming typical physical conditions rho / <rho> = 10, T = 10^4 - 10^4.5 K, and [Fe/H]= -2 yields reasonable agreement with the line strengths of high-ionization species, but it underpredicts the low-ionization species by two orders of magnitude or more. This discrepancy suggests that the low ionization lines arise in dense, cool, metal-rich clumps, present in some absorption systems.

preprint2010arXiv

The intergalactic medium over the last 10 billion years I: Lyman alpha absorption and physical conditions

The intergalactic medium (IGM) is the dominant reservoir of baryons at all cosmic epochs. We investigate the evolution of the IGM from z=2-0 in 48 Mpc/h, 110-million particle cosmological hydrodynamic simulations using three prescriptions for galactic outflows. We focus on the evolution of IGM physical properties, and how such properties are traced by Ly-alpha absorption as detectable using HST/COS. Our results broadly confirm the canonical picture that most Ly-alpha absorbers arise from highly ionized gas tracing filamentary large-scale structure. Growth of structure causes gas to move from the diffuse photoionized IGM into other cosmic phases, namely stars, cold and hot gas within galaxy halos, and the unbound and shock-heated warm-hot intergalactic medium (WHIM). By today, baryons are roughly equally divided between bound phases (35%), the diffuse IGM (41%), and the WHIM (24%). Here we (re)define the WHIM as gas with overdensities lower than that in halos and temperatures >10^5 K, in order to more closely align it with "missing baryons". When we tune our photoionizing background to match the observed evolution of the Ly-alpha mean flux decrement, we obtain a line count evolution that broadly agrees with available data. We predict a column density distribution slope of -1.70 for our favored momentum-driven wind model, in agreement with recent observations, and it becomes shallower with redshift. With improved statistics, the frequency of strong lines can be a valuable diagnostic of outflows, and our favored wind model matches existing data best among our models. The relationship between column density and physical density is fairly tight from z=2-0, and evolves as rho N_HI^0.74 10^(-0.37z) for diffuse absorbers. Linewidths only loosely reflect the temperature of the absorbing gas, which will hamper attempts to quantify the WHIM using broad Ly-alpha absorbers. [Abridged]

preprint2010arXiv

The Nature of Sub-millimetre Galaxies in Cosmological Hydrodynamic Simulations

We study the nature of rapidly star-forming galaxies at z=2 in cosmological hydrodynamic simulations, and compare their properties to observations of sub-millimetre galaxies (SMGs). We identify simulated SMGs as the most rapidly star-forming systems that match the observed number density of SMGs. In our models, SMGs are massive galaxies sitting at the centres of large potential wells, being fed by smooth infall and gas-rich satellites at rates comparable to their star formation rates (SFR). They are not typically undergoing major mergers that significantly boost their quiescent SFR, but they still often show complex gas morphologies and kinematics. Our simulated SMGs have stellar masses of log M*/Mo~11-11.7, SFRs of ~180-500 Mo/yr, a clustering length of 10 Mpc/h, and solar metallicities. The SFRs are lower than those inferred from far-IR data by a factor of 3, which we suggest may owe to one or more systematic effects in the SFR calibrations. SMGs at z=2 live in ~10^13 Mo halos, and by z=0 they mostly end up as brightest group galaxies in ~10^14 Mo halos. We predict that higher-M* SMGs should have on average lower specific SFRs, less disturbed morphologies, and higher clustering. We also predict that deeper far-IR surveys will smoothly join SMGs onto the massive end of the SFR-M* relationship defined by lower-mass z=2 galaxies. Overall, our simulated rapid star-formers provide as good a match to available SMG data as merger-based scenarios, offering an alternative scenario that emerges naturally from cosmological simulations.

preprint2010arXiv

The Sloan Digital Sky Survey Quasar Catalog V. Seventh Data Release

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

preprint2009arXiv

A quantitative explanation of the observed population of Milky Way satellite galaxies

We revisit the well known discrepancy between the observed number of Milky Way (MW) dwarf satellite companions and the predicted population of cold dark matter (CDM) sub-halos, in light of the dozen new low luminosity satellites found in SDSS imaging data and our recent calibration of the SDSS satellite detection efficiency, which implies a total population far larger than these dozen discoveries. We combine a dynamical model for the CDM sub-halo population with simple, physically motivated prescriptions for assigning stellar content to each sub-halo, then apply observational selection effects and compare to the current observational census. As expected, models in which the stellar mass is a constant fraction F(Omega_b/Omega_m) of the sub-halo mass M_sat at the time it becomes a satellite fail for any choice of F. However, previously advocated models that invoke suppression of gas accretion after reionization in halos with circular velocity v_c <~ 35 km/s can reproduce the observed satellite counts for -15 < M_V < 0, with F ~ 10^{-3}. Successful models also require strong suppression of star formation BEFORE reionization in halos with v_c <~ 10 km/s; models without pre-reionization suppression predict far too many satellites with -5 < M_V < 0. Our models also reproduce the observed stellar velocity dispersions ~ 5-10 km/s of the SDSS dwarfs given the observed sizes of their stellar distributions, and model satellites have M(<300 pc) ~ 10^7 M_sun as observed even though their present day total halo masses span more than two orders of magnitude. Our modeling shows that natural physical mechanisms acting within the CDM framework can quantitatively explain the properties of the MW satellite population as it is presently known, thus providing a convincing solution to the `missing satellite' problem.

preprint2009arXiv

Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

preprint2009arXiv

Binary Quasars at High Redshift II: Sub-Mpc Clustering at z ~ 3-4

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

preprint2009arXiv

Clustering of Low-Redshift (z <= 2.2) Quasars from the Sloan Digital Sky Survey

We present measurements of the quasar two-point correlation function, ξ_{Q}, over the redshift range z=0.3-2.2 based upon data from the SDSS. Using a homogeneous sample of 30,239 quasars with spectroscopic redshifts from the DR5 Quasar Catalogue, our study represents the largest sample used for this type of investigation to date. With this redshift range and an areal coverage of approx 4,000 deg^2, we sample over 25 h^-3 Gpc^3 (comoving) assuming the current LCDM cosmology. Over this redshift range, we find that the redshift-space correlation function, xi(s), is adequately fit by a single power-law, with s_{0}=5.95+/-0.45 h^-1 Mpc and γ_{s}=1.16+0.11-0.16 when fit over s=1-25 h^-1 Mpc. Using the projected correlation function we calculate the real-space correlation length, r_{0}=5.45+0.35-0.45 h^-1 Mpc and γ=1.90+0.04-0.03, over scales of rp=1-130 h^-1 Mpc. Dividing the sample into redshift slices, we find very little, if any, evidence for the evolution of quasar clustering, with the redshift-space correlation length staying roughly constant at s_{0} ~ 6-7 h^-1 Mpc at z<2.2 (and only increasing at redshifts greater than this). Comparing our clustering measurements to those reported for X-ray selected AGN at z=0.5-1, we find reasonable agreement in some cases but significantly lower correlation lengths in others. We find that the linear bias evolves from b~1.4 at z=0.5 to b~3 at z=2.2, with b(z=1.27)=2.06+/-0.03 for the full sample. We compare our data to analytical models and infer that quasars inhabit dark matter haloes of constant mass M ~2 x 10^12 h^-1 M_Sol from redshifts z~2.5 (the peak of quasar activity) to z~0. [ABRIDGED]

preprint2009arXiv

Cosmological Constraints from the Clustering of the Sloan Digital Sky Survey DR7 Luminous Red Galaxies

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find Ω_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and Ω_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find Ω_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_ν < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find Ω_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

preprint2009arXiv

Cosmological Constraints from the SDSS maxBCG Cluster Catalog

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog to simultaneously constrain cosmology and the richness--mass relation of the clusters. Assuming a flat \LambdaCDM cosmology, we find σ_8(Ω_m/0.25)^{0.41} = 0.832\pm 0.033 after marginalization over all systematics. In common with previous studies, our error budget is dominated by systematic uncertainties, the primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and uncertainty in the scatter of the richness--mass relation. Our constraints are fully consistent with the WMAP five-year data, and in a joint analysis we find σ_8=0.807\pm 0.020 and Ω_m=0.265\pm 0.016, an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster abundances. The remarkable consistency among these results demonstrates that cluster abundance constraints are not only tight but also robust, and highlight the power of optically-selected cluster samples to produce precision constraints on cosmological parameters.

preprint2009arXiv

Galaxies in a simulated $Λ$CDM Universe I: cold mode and hot cores

We study the formation of galaxies in a (50 Mpc/h)^3 cosmological simulation (2x288^3 particles), evolved using the entropy conserving SPH code Gadget-2. Most of the baryonic mass in galaxies of all masses is originally acquired through filamentary "cold mode" accretion of gas that was never shock heated to its halo virial temperature, confirming the key feature of our earlier results obtained with a different SPH code (Keres et al. 2005). Atmospheres of hot, virialized gas develop in halos above ~2.5e11 Msun, a transition mass that is nearly constant from z=3 to z=0. Cold accretion persists in halos above the transition mass, especially at z>=2. It dominates the growth of galaxies in low mass halos at all times, and it is the main driver of the cosmic star formation history. Satellite galaxies have accretion rates similar to central galaxies of the same baryonic mass at high redshifts, but they have less accretion than comparable central galaxies at low redshift. Relative to our earlier results, the Gadget-2 simulations predict much lower rates of "hot mode" accretion from the virialized gas component of massive halos. At z<=1, typical hot accretion rates in halos above 5e12 Msun are below 1 Msun/yr, even though our simulation does not include AGN heating or other forms of "preventive" feedback. The inner density profiles of hot gas in these halos are shallow, with long associated cooling times. The cooling recipes typically used in semi-analytic models can overestimate the accretion rates in these halos by orders of magnitude, so such models may overemphasize the role of preventive feedback in producing observed galaxy masses and colors. A fraction of the massive halos develop cuspy profiles and significant cooling rates between z=1 and z=0, a redshift trend similar to the observed trend in the frequency of cooling flow clusters.

preprint2009arXiv

Galaxies in a Simulated $Λ$CDM Universe II: Observable Properties and Constraints on Feedback

We compare the properties of galaxies that form in a cosmological simulation without strong feedback to observations at z=0. We confirm previous findings that models without strong feedback overproduce the observed galaxy baryonic mass function, especially at the low and high mass extremes. Through post-processing we investigate what kinds of feedback would be required to reproduce observed galaxy masses and star formation rates. To mimic an extreme form of "preventive" feedback (e.g., AGN radio mode) we remove all baryonic mass that was originally accreted via "hot mode" from shock-heated gas. This does not bring the high mass end of the galaxy mass function into agreement with observations because much of the stellar mass in these systems formed at high redshift from baryons that originally accreted via "cold mode" onto lower mass progenitors. An efficient "ejective" feedback mechanism, such as supernova driven winds, must reduce the masses of these progenitors. Feedback must also reduce the masses of lower mass z=0 galaxies, which assemble at lower redshifts and have much lower star formation rates. If we monotonically re-map galaxy masses to reproduce the observed mass function, but retain the simulation's predicted star formation rates, we obtain fairly good agreement with the observed sequence of star-forming galaxies but fail to recover the observed population of passive, low star formation rate galaxies. Suppressing all hot mode accretion improves agreement for high mass galaxies but worsens the agreement at intermediate masses. Reproducing these z=0 observations requires a feedback mechanism that dramatically suppresses star formation in a fraction of galaxies, increasing with mass, while leaving star formation rates of other galaxies essentially unchanged.

preprint2009arXiv

Halo Occupation Distribution Modeling of Clustering of Luminous Red Galaxies

We perform Halo Occupation Distribution (HOD) modeling to interpret small-scale and intermediate-scale clustering of 35,000 luminous early-type galaxies and their cross-correlation with a reference imaging sample of normal L* galaxies in the Sloan Digital Sky Survey. The modeling results show that most of these luminous red galaxies (LRGs) are central galaxies residing in massive halos of typical mass M ~ a few times 10^13 to 10^14 Msun/h, while a few percent of them have to be satellites within halos in order to produce the strong auto-correlations exhibited on smaller scales. The mean luminosity Lc of central LRGs increases with the host halo mass, with a rough scaling relation of Lc \propto M^0.5. The halo mass required to host on average one satellite LRG above a luminosity threshold is found to be about 10 times higher than that required to host a central LRG above the same threshold. We find that in massive halos the distribution of L* galaxies roughly follows that of the dark matter and their mean occupation number scales with halo mass as M^1.5. The HOD modeling results also allows for an intuitive understanding of the scale-dependent luminosity dependence of the cross-correlation between LRGs and L_* galaxies. Constraints on the LRG HOD provide tests to models of formation and evolution of massive galaxies, and they are also useful for cosmological parameter investigations. In one of the appendices, we provide LRG HOD parameters with dependence on cosmology inferred from modeling the two-point auto-correlation functions of LRGs.

preprint2009arXiv

Lyman-alpha Emission From Cosmic Structure I: Fluorescence

We present predictions for the fluorescent Lyman-alpha emission signature arising from photoionized, optically thick structures in Smoothed Particle Hydrodynamic (SPH) cosmological simulations of a Lambda-CDM universe using a Monte Carlo Lyman-alpha radiative transfer code. We calculate the expected Lyman-alpha image and 2-dimensional spectra for gas exposed to a uniform ultraviolet ionizing background as well as gas exposed additionally to the photoionizing radiation from a local quasar, after correcting for the self-shielding of hydrogen. As a test of our numerical methods and for application to current observations, we examine simplified analytic structures that are uniformly or anisotropically illuminated. We compare these results with recent observations. We discuss future observing campaigns on large telescopes and realistic strategies for detecting fluorescence owing to the ambient metagalactic ionization and in regions close to bright quasars. While it will take hundreds of hours on the current generation of telescopes to detect fluorescence caused by the ultraviolet background (UVB) alone, our calculations suggest that of order ten sources of quasar-induced fluorescent Lyman-alpha emission should be detectable after a 10 hour exposure in a 10 arcmin^2 field around a bright quasar. These observations will help probe the physical conditions in the densest regions of the intergalactic medium as well as the temporal light curves and isotropy of quasar radiation.

preprint2008arXiv

Extending Recovery of the Primordial Matter Power Spectrum

The shape of the primordial matter power spectrum Plin(k) encodes critical information on cosmological parameters. At large scales, the observable galaxy power spectrum Pobs(k) is expected to follow the shape of Plin(k), but on smaller scales the effects of nonlinearity and galaxy bias make the ratio Pobs(k)/Plin(k) scale-dependent. We develop a method that can extend the dynamic range of the Plin(k) recovery by incorporating constraints on the galaxy halo occupation distribution (HOD) from the projected galaxy correlation function wp. We devise an analytic model to calculate Pobs(k) in real-space and redshift-space. Once HOD parameters are determined by matching wp for a given cosmological model, galaxy bias is completely specified, and our analytic model predicts both the shape and normalization of Pobs(k). Applying our method to SDSS main galaxy samples, we find that the real-space Pobs(k) follows the shape of the nonlinear matter power spectrum at the 1-2% level up to k=0.2 h/Mpc. When we apply our method to SDSS LRG samples, the linear bias approximation is accurate to 5% at k<0.08 h/Mpc, but the scale-dependence of LRG bias prevents the use of linear theory at k>0.08 h/Mpc. Our HOD model prediction is in good agreement with the recent SDSS LRG Pobs(k) measurements at all measured scales (k<0.2 h/Mpc), naturally explaining the shape of Pobs(k). The "Q-model" prescription is a poor description of galaxy bias for the LRG samples, and it can lead to biased cosmological parameter estimates when measurements at k>0.1 h/Mpc are included in the analysis. We quantify the potential bias and constraints on cosmological parameters that arise from applying linear theory and Q-model fitting, and we demonstrate the utility of HOD modeling of future high precision measurements of Pobs(k) on quasi-linear scales.

preprint2008arXiv

Quasar Clustering from SDSS DR5: Dependences on Physical Properties

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2σ$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$.

preprint2008arXiv

Self-Consistent Models of the AGN and Black Hole Populations: Duty Cycles, Accretion Rates, and the Mean Radiative Efficiency

We construct evolutionary models of the populations of AGN and supermassive black holes, in which the black hole mass function grows at the rate implied by the observed luminosity function, given assumptions about the radiative efficiency and the Eddington ratio. We draw on a variety of recent X-ray and optical measurements to estimate the bolometric AGN luminosity function and compare to X-ray background data and the independent estimate of Hopkins et al. (2007) to assess remaining systematic uncertainties. The integrated AGN emissivity closely tracks the cosmic star formation history, suggesting that star formation and black hole growth are closely linked at all redshifts. Observational uncertainties in the local black hole mass function remain substantial, with estimates of the integrated black hole mass density ρ_BH spanning the range 3-5.5x10^5 Msun/Mpc^3. We find good agreement with estimates of the local mass function for a reference model where all active black holes have efficiency \eps=0.065 and L_bol/L_Edd~0.4. In this model, the duty cycle of 10^9 Msun black holes declines from 0.07 at z=3 to 0.004 at z=1 and 0.0001 at z=0. The decline is shallower for less massive black holes, a signature of "downsizing" evolution in which more massive black holes build their mass earlier. The predicted duty cycles and AGN clustering bias in this model are in reasonable accord with observational estimates. If the typical Eddington ratio declines at z<2, then the "downsizing" of black hole growth is less pronounced. Matching the integrated AGN emissivity to the local black hole mass density implies \eps=0.075 (ρ_BH/4.5x10^5 Msun/Mpc^3)^{-1} for our standard luminosity function estimate (25% higher for Hopkins et al.'s), lower than the values \eps=0.16-0.20 predicted by MHD simulations of disk accretion.

preprint2008arXiv

The Growth of Central and Satellite Galaxies in Cosmological Smoothed Particle Hydrodynamics Simulations

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to 7e9 M[Sun]. Most friends-of-friends halos in the simulation have a distinct central galaxy, typically two to five times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive halos form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5-1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Since z=1, 27% of central galaxies (above 3e10 M[Sun]) and 22% of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain "central" objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with the observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.

preprint2006arXiv

Baryon Dynamics, Dark Matter Substructure, and Galaxies

By comparing a collisionless cosmological N-body simulation (DM) to an SPH simulation with the same initial conditions, we investigate the correspondence between the dark matter subhalos produced by collisionless dynamics and the galaxies produced by dissipative gas dynamics in a dark matter background. When galaxies in the SPH simulation become satellites in larger groups, they retain local dark matter concentrations (SPH subhalos) whose mass is typically five times their baryonic mass. The more massive subhalos of the SPH simulation have corresponding subhalos of similar mass and position in the DM simulation; at lower masses, there is fairly good correspondence, but some DM subhalos are in different spatial positions and some suffer tidal stripping or disruption. The halo occupation statistics of DM subhalos -- the mean number of subhalos, pairs, and triples as a function of host halo mass -- are very similar to those of SPH subhalos and SPH galaxies. Gravity of the dissipative baryon component amplifies the density contrast of subhalos in the SPH simulation, making them more resistant to tidal disruption. Relative to SPH galaxies and SPH subhalos, the DM subhalo population is depleted in the densest regions of the most massive halos. The good agreement of halo occupation statistics between the DM subhalo and SPH galaxy populations leads to good agreement of their two-point correlation functions and higher order moments on large scales. The depletion of DM subhalos in dense regions depresses their clustering at R<1 Mpc/h. In these simulations, the "conversation" between dark matter and baryons is mostly one-way, with dark matter dynamics telling galaxies where to form and how to cluster, but the "back talk" of the baryons influences small scale clustering by enhancing the survival of substructure in the densest environments.

preprint2006arXiv

Black Hole Masses and Eddington Ratios at 0.3<z<4

We study the distribution of Eddington luminosity ratios, L_bol/L_edd, of active galactic nuclei (AGNs) discovered in the AGN and Galaxy Evolution Survey (AGES). We combine H-beta, MgII, and CIV line widths with continuum luminosities to estimate black hole (BH) masses in 407 AGNs, covering the redshift range z~0.3-4 and the bolometric luminosity range L_bol~10^45-10^47 erg/s. The sample consists of X-ray or mid-infrared (24 micron) point sources with optical magnitude R<=21.5 mag and optical emission line spectra characteristic of AGNs. For the range of luminosity and redshift probed by AGES, the distribution of estimated Eddington ratios is well described as log-normal with a peak at L_bol/L_edd ~ 1/4 and a dispersion of 0.3 dex. Since additional sources of scatter are minimal, this dispersion must account for contributions from the scatter between estimated and true BH mass and the scatter between estimated and true bolometric luminosity. Therefore, we conclude that: (1) neither of these sources of error can contribute more than \~0.3 dex rms; and (2) the true Eddington ratios of optically luminous AGNs are even more sharply peaked. Because the mass estimation errors must be smaller than ~0.3 dex, we can also investigate the distribution of Eddington ratios at fixed BH mass. We show for the first time that the distribution of Eddington ratios at fixed BH mass is peaked, and that the dearth of AGNs at a factor ~10 below Eddington is real and not an artifact of sample selection. These results provide strong evidence that supermassive BHs gain most of their mass while radiating close to the Eddington limit, and they suggest that the fueling rates in luminous AGNs are ultimately determined by BH self-regulation of the accretion flow rather than galactic scale dynamical disturbances.

preprint2006arXiv

Cosmic Voids and Galaxy Bias in the Halo Occupation Framework

(Abridged) We investigate the power of void statistics to constrain galaxy bias and the amplitude of dark matter fluctuations. We use the halo occupation distribution (HOD) framework to describe the relation between galaxies and dark matter. After choosing HOD parameters that reproduce the mean space density n_gal and projected correlation function w_p measured for galaxy samples with M_r<-19 and M_r<-21 from the Sloan Digital Sky Survey (SDSS), we predict the void probability function (VPF) and underdensity probability function (UPF) of these samples by populating the halos of a large, high-resolution N-body simulation. If we make the conventional assumption that the HOD is independent of large scale environment at fixed halo mass, then models constrained to match n_gal and w_p predict nearly identical void statistics, independent of the scatter between halo mass and central galaxy luminosity or uncertainties in HOD parameters. Models with sigma_8=0.7 and sigma_8=0.9 also predict very similar void statistics. However, the VPF and UPF are sensitive to environmental variations of the HOD in a regime where these variations have little impact on w_p. For example, doubling the minimum host halo mass in regions with large scale (5 Mpc/h) density contrast delta<-0.65 has a readily detectable impact on void probabilities of M_r<-19 galaxies, and a similar change for delta<-0.2 alters the void probabilities of M_r<-21 galaxies at a detectable level. The VPF and UPF provide complementary information about the onset and magnitude of density- dependence in the HOD. By detecting or ruling out HOD changes in low density regions, void statistics can reduce systematic uncertainties in the cosmological constraints derived from HOD modeling, and, more importantly, reveal connections between halo formation history and galaxy properties.

preprint2006arXiv

On the Luminosity Dependence of the Galaxy Pairwise Velocity Dispersion

(Abridged) We make predictions for the pairwise velocity dispersion (PVD) of galaxies with models that are constrained to match the projected correlation function and luminosity function of galaxies in the Two-Degree Field Galaxy Redshift Survey (2dFGRS). We use these data to constrain the halo occupation distribution (HOD), then calculate the PVD by populating the halos of a high resolution N-body simulation. We examine the luminosity and scale dependence of the predicted PVD. At r<1 Mpc/h and r>3 Mpc/h, we find that the PVD decreases with increasing galaxy luminosity. This result is mostly driven by the fraction of satellite galaxies f_sat, which is well-constrained by the correlation function. We find f_sat~25% for galaxies fainter than L_star, while for brighter galaxies the satellite fraction rapidly declines, creating the decrease in the PVD with luminosity. At r=1 Mpc/h, the PVD has no dependence on luminosity because satellite galaxies dominate the statistics for all objects. Recent measurements of the PVD in Fourier space using the "dispersion model" have reported a strong decline in PVD with increasing luminosity at k=1 h/Mpc. We test this method with our HOD models, finding that there is no consistent comparison between the PVD at a given k and the true dispersion at a given value of r. This results in a luminosity dependence in k-space that is stronger than in configuration space. The luminosity dependence of the HOD results in Fourier space are consistent with those measured at k=1 h/Mpc; thus the recent measurements of the PVD are fully explainable in the context of halo occupation models. The normalization of the PVD is lower than predicted by our fiducial model, and reproducing it requires a lower value of Omega_m (~0.2 instead of 0.3), a lower value of sigma_8 (~0.7 instead of 0.9), or strong velocity bias.

preprint2005arXiv

On the Mass-to-Light Ratio of Large Scale Structure

We examine the dependence of the mass-to-light (M/L) ratio of large-scale structure on cosmological parameters, in models that are constrained to match observations of the projected galaxy correlation function w(rp). For a sequence of cosmological models with a fixed P(k) shape and increasing normalization \sig8, we find parameters of the galaxy halo occupation distribution (HOD) that reproduce SDSS measurements of w(rp) as a function of luminosity. Using these HOD models we calculate mean M/L ratios as a function of halo mass and populate halos of N-body simulations to compute M/L in larger scale environments, including cluster infall regions. For all cosmological models, the M/L ratio in high mass halos or high density regions is approximately independent of halo mass or smoothing scale. However, the "plateau" value of M/L depends on \sig8 as well as Ω_m, and it represents the universal mass-to-light ratio <M/L> only for models in which the galaxy correlation function is approximately unbiased, i.e., with \sig8 ~ \sig8_gal. Our results for cluster mass halos follow the trend M/L = 577(Ω_m/0.3)(\sig8/0.9)^{1.7} h Msun/Lsun. Combined with Carlberg et al.'s (1996) mean M/L ratio of CNOC galaxy clusters, this relation implies (\sig8/0.9)(Ω_m/0.3)^{0.6} = 0.75 +/- 0.06. M/L ratios of clusters from the SDSS and CAIRNS surveys yield similar results. This constraint is inconsistent with parameter values Ω_m ~ 0.3, \sig8 ~ 0.9 favored by recent joint analyses of CMB measurements and other large-scale structure data. We discuss possible resolutions, none of which seems entirely satisfactory. Appendices present an improved formula for halo bias factors and an improved analytic technique for calculating the galaxy correlation function from a given cosmological model and HOD. (Abridged)

preprint2003arXiv

Accretion Driven Evolution of Quasars and Black Holes: Theoretical Models

We present a flexible framework for constructing physical models of quasar evolution that can incorporate a variety of observational constraints, such as multi-wavelength luminosity functions, estimated masses and accretion rates of active black holes, space densities of quasar hosts, and the local black hole mass function. We focus on the accretion rate distribution p(mdot|M,z), the probability that a black hole of mass M at redshift z accretes at a rate mdot in Eddington units. Given the radiative efficiency as a function of mdot, the quasar luminosity function (QLF) is determined by a convolution of p(mdot|M,z) with the black hole mass function n(M,z). In the absence of mergers p(mdot|M,z) also determines the full evolution of n(M,z), given a "boundary value" of n(M) at some redshift. Matching the observed decline of the QLF break luminosity at z<2 requires either a shift towards lower characteristic accretion rates or an evolving mass dependence of p(mdot) that preferentially shuts off accretion onto high mass black holes. These two scenarios make different predictions for the masses and accretion rates of active black holes. If the first mechanism dominates, then the QLF changes character between z=2 and z=0, shifting from a sequence of black hole mass towards a sequence of L/Ledd. We construct and compare five models that illustrate different assumptions about the quasar population: short and long lifetime models dominated by unobscured thin-disk accretion, a model with a high fraction of obscured quasars, and models in which mergers or ADAF accretion produce substantial black hole growth at low redshift. We discuss the observational advances that would be most valuable for distinguishing such models and for pinning down the physics that drives black hole and quasar evolution. (Abridged)

preprint2003arXiv

The Halo Occupation Distribution and the Physics of Galaxy Formation

The halo occupation distribution (HOD) describes the bias between galaxies and dark matter by specifying (a) the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class and (b) the relative distributions of galaxies and dark matter within halos. We calculate predicted HODs for a Lambda-CDM cosmological model using an SPH hydrodynamic simulation and a semi-analytic (SA) galaxy formation model. Although the two methods predict different galaxy mass functions, their HOD predictions agree remarkably well. For mass-selected samples, the mean occupation <N(M)> exhibits a sharp cutoff at low halo masses, a slowly rising plateau for <N>~1-2, and a more steeply rising high occupancy regime. At low <N>, the mean pair and triple counts are well below Poisson expectations, with important consequences for small scale behavior of 2- and 3-point correlation functions. The HOD depends strongly on galaxy age, with high mass halos populated mainly by old galaxies and low mass halos by young galaxies. The SPH simulation supports several simplifying assumptions about HOD bias: the most massive galaxy in a halo usually lies close to the center and moves near the halo's mean velocity; satellite galaxies have the same radial profile and velocity dispersion as the dark matter; and the mean occupation at fixed halo mass is independent of the halo's larger scale environment. By applying the SPH and SA HODs to a large volume N-body simulation, we show that both methods predict slight, observable departures from a power-law galaxy correlation function. The predicted HODs are closely tied to the underlying galaxy formation physics, they offer useful guidance to theoretical models of galaxy clustering, and they will be tested empirically by ongoing analyses of galaxy redshift surveys. (Shortened)

preprint2002arXiv

How Do Galaxies Get Their Gas?

Not the way one might have thought. In hydrodynamic simulations of galaxy formation, some gas follows the traditionally envisioned route, shock heating to the halo virial temperature before cooling to the much lower temperature of the neutral ISM. But most gas enters galaxies without ever heating close to the virial temperature, gaining thermal energy from weak shocks and adiabatic compression, and radiating it just as quickly. This ``cold mode'' accretion is channeled along filaments, while the conventional, ``hot mode'' accretion is quasi-spherical. Cold mode accretion dominates high redshift growth by a substantial factor, while at z<1 the overall accretion rate declines and hot mode accretion has greater relative importance. The decline of the cosmic star formation rate at low z is driven largely by geometry, as the typical cross section of filaments begins to exceed that of the galaxies at their intersections.

preprint2002arXiv

Lyman Alpha Absorber Correlations and the Bias of the Lyman Alpha Forest

Lyman alpha absorber correlations contain information about the underlying density distribution associated with a particular class of absorbers. As such, they provide an opportunity to independently measure the "bias" of the Lyman alpha forest, i.e. the relationship between HI column density and underlying dark matter density. In these proceedings we use hydrodynamic simulations to investigate whether the evolution of this bias is measurable from observable correlations. Unfortunately, the increasingly complex physics in the IGM at z<1 makes a direct measurement of the bias difficult. Nevertheless, current simulations do make predictions for HI absorber correlations that are in broad agreement with observations at both high and low redshift, thus reinforcing the bias evolution predictions given by these models.

preprint2002arXiv

Towards a Precise Measurement of Matter Clustering: Lyman-alpha Forest Data at Redshifts 2-4

We measure the filling factor, correlation function, and power spectrum of transmitted flux in a large sample of Lya forest spectra, comprised of 30 Keck HIRES spectra and 23 Keck LRIS spectra. We infer the linear matter power spectrum P(k) from the flux power spectrum P_F(k), using an improved version of the method of Croft et al. (1998) that accounts for the influence of z-space distortions, non- linearity, and thermal broadening on P_F(k). The evolution of the shape and amplitude of P(k) over the range z= 2-4 is consistent with gravitational instability, implying that non-gravitational fluctuations do not make a large contribution. Our fiducial measurement of P(k) comes from data with <z> = 2.72. It has amplitude Delta^2(k_p)=0.74^0.20_-0.16 at wavenumber k_p=0.03 (km/s)^-1 and is well described by a power-law of index -2.43 +/- 0.06 or by a CDM-like power spectrum with shape parameter Gamma'=1.3^+0.7_-0.5*10^-3 (km/s) at z=2.72. For Omega_m=0.4, Omega_Lam=0.6, the best-fit Gamma =0.16 (h^-1mpc)^-1, in good agreement with the 2dF Galaxy Redshift Survey, and the best-fit sigma_8=0.82 (Gamma/0.15)^-0.44. Matching the observed cluster mass function and our Delta^2(k_p) in spatially flat models requires Omega_m=0.38^+0.10_-0.08 + 2.2 (Gamma-0.15). Matching Delta^2(k_p) in COBE-normalized, flat CDM models with no tensor fluctuations requires Omega_m = (0.29 +/-0.04) n^-2.89 h_65^-1.9. The Lya forest complements other probes of P(k) by constraining a regime of redshift and lengthscale not accessible by other means, and the consistency of these inferred parameters with independent estimates provides further support for inflation, cold dark matter, and vacuum energy (abridged).

preprint2001arXiv

Simulations of Damped Lyman-Alpha and Lyman Limit Absorbers in Different Cosmologies: Implications for Structure Formation at High Redshift

We use hydrodynamic cosmological simulations to study damped Lyman-alpha (DLA) and Lyman limit (LL) absorption at redshifts z=2-4 in five variants of the cold dark matter scenario. Our standard simulations resolve the formation of dense concentrations of neutral gas in halos with circular velocity v_c roughly 140 km/s for Omega_m=1 and 90 km/s for Omega_m=0.4, at z=2; an additional LCDM simulation resolves halos down to v_c approximately 50 km/s at z=3. We find a clear relation between HI column density and projected distance to the center of the nearest galaxy, with DLA absorption usually confined to galactocentric radii less than 10-15 kpc and LL absorption arising out to projected separations of 30 kpc or more. Detailed examination provides evidence of non-equilibrium effects on absorption cross-section. If we consider only absorption in the halos resolved by our standard simulations, then all five models fall short of reproducing the observed abundance of DLA and LL systems at these redshifts. If we extrapolate to lower halo masses, we find all four models are consistent with the observed abundance of DLA systems if the the extrapolated behavior extends to circular velocities roughly 50-80 km/s, and they may produce too much absorption if the relation continues to 40 km/s. Our results suggest that LL absorption is closely akin to DLA absorption, arising in less massive halos or at larger galactocentric radii but not caused by processes acting on a radically different mass scale.

preprint2000arXiv

Baryons in the Warm-Hot Intergalactic Medium

Approximately 30-40% of all baryons in the present day universe reside in a warm-hot intergalactic medium (WHIM), with temperatures between 10^5<T<10^7 K. This is a generic prediction from six hydrodynamic simulations of currently favored structure formation models having a wide variety of numerical methods, input physics, volumes, and spatial resolutions. Most of these warm-hot baryons reside in diffuse large-scale structures with a median overdensity around 10-30, not in virialized objects such as galaxy groups or galactic halos. The evolution of the WHIM is primarily driven by shock heating from gravitational perturbations breaking on mildly nonlinear, non-equilibrium structures such as filaments. Supernova feedback energy and radiative cooling play lesser roles in its evolution. WHIM gas is consistent with observations of the 0.25 keV X-ray background without being significantly heated by non-gravitational processes because the emitting gas is very diffuse. Our results confirm and extend previous work by Cen & Ostriker and Dave' et al.

preprint2000arXiv

Biased Galaxy Formation And Measurements Of Beta

Measurements of the cosmological density parameter Omega using techniques that exploit the gravity-induced motions of galaxies constrain, in linear perturbation theory, the degenerate parameter combination beta = Omega^{0.6}/b, where the linear bias parameter b is the ratio of the fluctuation amplitudes of the galaxy and mass distributions. However, the relation between the galaxy and mass density fields depends on the complex physics of galaxy formation, and it can in general be non-linear, stochastic, and perhaps non-local. The one-parameter linear bias model is almost certainly oversimplified, which leads to the obvious question: What is the quantity beta that is actually measured by different techniques? To address this question, we estimate beta from galaxy distributions that are constructed by applying a variety of locally biased galaxy formation models to cosmological N-body simulations. We compare the values of beta estimated using three different techniques: a density-density comparison similar to the POTENT analysis, a velocity-velocity comparison similar to the VELMOD analysis, and an anisotropy analysis of the redshift-space power spectrum. In most cases, we find that beta estimated using all three methods is similar to the asymptotic value of Omega^{0.6}/b_{sigma}(R) at large R, where b_{sigma}(R) is the ratio of rms galaxy fluctuations to rms mass fluctuations on scale R. Thus, something close to the conventional interpretation of beta continues to hold even for complex bias models. Moreover, we find that beta estimates made using these three methods should, in principle, agree with each other. It is thus unlikely that non-linear or scale-dependent bias is responsible for the discrepancies that exist among current measurements of beta from different techniques.

preprint2000arXiv

Cooling Radiation and the Lyman-alpha Luminosity of Forming Galaxies

We examine the cooling radiation from forming galaxies in hydrodynamic simulations of the LCDM model (cold dark matter with a cosmological constant), focusing on the Ly-alpha line luminosities of high-redshift systems. Primordial composition gas condenses within dark matter potential wells, forming objects with masses and sizes comparable to the luminous regions of observed galaxies. As expected, the energy radiated in this process is comparable to the gravitational binding energy of the baryons, and the total cooling luminosity of the galaxy population peaks at z ~= 2. However, in contrast to the classical picture of gas cooling from the \sim 10^6 K virial temperature of a typical dark matter halo, we find that most of the cooling radiation is emitted by gas with T < 20,000 K. As a consequence, roughly 50% of this cooling radiation emerges in the Ly-alpha line. While a galaxy's cooling luminosity is usually smaller than the ionizing continuum luminosity of its young stars, the two are comparable in the most massive systems, and the cooling radiation is produced at larger radii, where the Ly-alpha photons are less likely to be extinguished by dust. We suggest, in particular, that cooling radiation could explain the two large (\sim 100 kpc), luminous (L_{Ly-alpha} \sim 10^{44} erg s^{-1}) ``blobs'' of Ly-alpha emission found in Steidel et al.'s (1999) narrow band survey of a z = 3 proto-cluster. Our simulations predict objects of the observed luminosity at about the right space density, and radiative transfer effects can account for the observed sizes and line widths. We discuss observable tests of this hypothesis for the nature of the Ly-alpha blobs, and we present predictions for the contribution of cooling radiation to the Ly-alpha luminosity function of galaxies as a function of redshift.

preprint2000arXiv

Topology from the Simulated Sloan Digital Sky Survey

We measure the topology (genus curve) of the galaxy distribution in a mock redshift catalog designed to resemble the upcoming Sloan Digital Sky Survey (SDSS). The catalog, drawn from a large N-body simulation of a Lambda-CDM cos- mological model, mimics the anticipated spectroscopic selection procedures of the SDSS in some detail. Sky maps, redshift slices, and 3-D contour maps of the mock survey reveal a rich and complex structure, including networks of voids and superclusters that resemble the patterns seen in the CfA redshift survey and the Las Campanas Redshift Survey (LCRS). The 3-D genus curve can be measured from the simulated catalog with superb precision; this curve has the general shape predicted for Gaussian, random phase initial conditions, but the error bars are small enough to demonstrate with high significance the subtle departures from this shape caused by non-linear gravitational evolution. These distortions have the form predicted by Matsubara's (1994) perturbative anal- ysis, but they are much smaller in amplitude. We also measure the 3-D genus curve of the radial peculiar velocity field measured by applying distance- indicator relations (with realistic errors) to the mock catalog. This genus curve is consistent with the Gaussian random phase prediction, though it is of relatively low precision because of the large smoothing length required to overcome noise in the measured velocity field. Finally, we measure the 2-D topology in redshift slices, similar to early slices from the SDSS and to slices already observed in the LCRS. The genus curves of these slices are consistent with the observed genus curves of the LCRS, providing further evidence in favor of the inflationary CDM model with Omega_M~0.4. The catalog is publicly available at http://www.astronomy.ohio-state.edu/~dhw/sdss.html.

preprint1998arXiv

Reconstruction Analysis of Galaxy Redshift Surveys: A Hybrid Reconstruction Method

In reconstruction analysis of galaxy redshift surveys, one works backwards from the observed galaxy distribution to the primordial density field in the same region, then evolves the primordial fluctuations forward in time with an N-body code. This incorporates assumptions about the cosmological parameters, the properties of primordial fluctuations, and the biasing relation between galaxies and mass. These can be tested by comparing the reconstruction to the observed galaxy distribution, and to peculiar velocity data. This paper presents a hybrid reconstruction method that combines the `Gaussianization'' technique of Weinberg(1992) with the dynamical schemes of Nusser & Dekel(1992) and Gramann(1993). We test the method on N-body simulations and on N-body mock catalogs that mimic the depth and geometry of the Point Source Catalog Redshift Survey and the Optical Redshift Survey. This method is more accurate than Gaussianization or dynamical reconstruction alone. Matching the observed morphology of clustering can limit the bias factor b, independent of Omega. Matching the cluster velocity dispersions and z-space distortions of the correlation function xi(s,mu) constrains the parameter beta=Omega^{0.6}/b. Relative to linear or quasi-linear approximations, a fully non-linear reconstruction makes more accurate predictions of xi(s,mu) for a given beta, thus reducing the systematic biases of beta measurements and offering further scope for breaking the degeneracy between Omega and b. It also circumvents the cosmic variance noise that limits conventional analyses of xi(s,mu). It can also improve the determination of Omega and b from joint analyses of redshift & peculiar velocity surveys as it predicts the fully non-linear peculiar velocity distribution at each point in z-space.

preprint1997arXiv

Large-scale structure in COBE-normalized cold dark matter cosmogonies

We study the clustering of the mass distribution in COBE-normalized open and flat CDM models using large N-body simulations. With an age of the universe of 14 Gyr (12 Gyr) for the flat (open) models and a baryon density fixed by nucleosynthesis constraints, the observed abundance of rich galaxy clusters leads to tight constraints on the density parameter; 0.25<Omega_0<0.4 (0.4<Omega_0<0.5) for flat (open) models. The upper limits on Omega_0 can be relaxed if one lowers the H_0 and increases the age of the universe, but h<0.25 is required to allow Omega_0=1. The constraints also change if one allows tilted primordial power spectra. An acceptable Omega_0=1 model with h=0.5 can be constructed, with a tilt of n=0.8. We compare the evolved mass correlation functions of the most promising models with those of APM galaxies. The flat models have steep correlation functions at small scales and require galaxies to be antibiased on scales r<8 Mpc/h. The open models require little or no antibias on small scales and a positive bias on large scales. The lack of a positive bias on small scales is difficult to reconcile with the virial analysis of galaxy clusters which, if Omega_0>0.2, implies that galaxies are overabundant in clusters relative to the field. The tilted Omega_0=1 model, on the other hand, does require that galaxies be positively biased on all scales. We also compute the topology of isodensity contours in these models, obtaining theoretical predictions that are less sensitive to galaxy bias.

preprint1995arXiv

Cosmological Simulations with TreeSPH

We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989) generalized to handle comoving coordinates and periodic boundary conditions, combines smoothed--particle hydrodynamics (SPH) with the hierarchical tree method for computing gravitational forces. The Lagrangian hydrodynamics approach and individual time steps for gas particles give the algorithm a large dynamic range, which is essential for studies of galaxy formation in a cosmological context. The code incorporates radiative cooling for an optically thin, primordial composition gas in ionization equilibrium with a user-specified ultraviolet background. We adopt a phenomenological prescription for star formation that gradually turns cold, dense, Jeans-unstable gas into collisionless stars, returning supernova feedback energy to the surrounding medium. In CDM simulations, some of the baryons that fall into dark matter potential wells dissipate their acquired thermal energy and condense into clumps with roughly galactic masses. The resulting galaxy population is insensitive to assumptions about star formation; we obtain similar baryonic mass functions and galaxy correlation functions from simulations with star formation and from simulations without star formation in which we identify galaxies directly from the cold, dense gas.

preprint1995arXiv

Hydrodynamic Simulations of Galaxy Formation. II. Photoionization and the Formation of Low Mass Galaxies

Photoionization by the high-redshift ultraviolet radiation background heats low density gas before it falls into dark matter potential wells, and it eliminates the neutral hydrogen and singly ionized helium that dominate cooling of primordial gas at temperatures of $10^4-10^5$K. We investigate the influence of photoionization on galaxy formation using high-resolution simulations with a 1-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code. We find that the presence of a photoionizing background suppresses the formation of galaxies with circular velocities $v_{circ} \simlt 30\,$\kms\ and substantially reduces the mass of cooled baryons in systems with circular velocities up to $v_{circ} \sim 50\,$\kms. Above $v_{circ} \sim 75\,$\kms, photoionization has no significant effect. Photoionization exerts its influence primarily by heating gas before collapse; the elimination of line cooling processes is less important. We discuss the implications of these results for hierarchical theories of galaxy formation.

preprint1995arXiv

Small Scale Structure and High Redshift HI

Cosmological simulations with gas dynamics suggest that the Lyman-alpha forest is produced mainly by "small scale structure" --- filaments and sheets that are the high redshift analog of today's galaxy superclusters. There is no sharp distinction between Lyman-alpha clouds and "Gunn-Peterson" absorption produced by the fluctuating IGM -- the Lyman-alpha forest {\it is} the Gunn-Peterson effect. Lyman limit and damped Lyman-alpha absorption arises in the radiatively cooled gas of forming galaxies. At $z~2-3$, most of the gas is in the photoionized, diffuse medium associated with the Lyman-alpha forest, but most of the {\it neutral} gas is in damped Lyman-alpha systems. We discuss generic evolution of cosmic gas in a hierarchical scenario of structure formation, with particular attention to the prospects for detecting 21cm emission from high redshift HI. A scaling argument based on the present-day cluster mass function suggests that objects with $M_{HI} >~ 5e11 h^{-1} \msun$ should be extremely rare at $z~3$, so detections with existing instruments will be difficult. An instrument like the proposed Square Kilometer Array could detect individual damped Lyman-alpha systems at high redshift, making it possible to map structure in the high redshift universe in much the same way that today's galaxy redshift surveys map the local large scale structure.

preprint1994arXiv

Constraints onOmega from the IRAS Redshift Surveys

We measure the anisotropy of the redshift-space power spectrum in the 1.2-Jy and QDOT redshift surveys of IRAS-selected galaxies. On large scales, this anisotropy is caused by coherent peculiar motions, and gravitational instability theory predicts a distortion of the power spectrum that depends only on the ratio $β\equiv f(Ω)/b \approx Ω^{0.6}/b$, where Omega is the cosmological density parameter and $b$ is the bias parameter. On small scales, the distortion is dominated by the random velocity dispersion in non-linear structures. We fit the observed anisotropy with an analytic model that incorporates two parameters, beta, and a small-scale velocity dispersion sigma_v. Tests on N-body simulations show that this model recovers beta quite accurately on the scales accessible to the existing IRAS redshift surveys. Applying our procedure to the 1.2-Jy and QDOT surveys, we find beta=0.52 +/- 0.13 and beta=0.54 +/- 0.3, respectively. These results imply Omega approximately 0.35 if galaxies trace mass, or a bias factor of about 2 if Omega=1.

preprint1993arXiv

Cosmic Voids and Biased Galaxy Formation

Using cosmological $N$-body simulations and the void probability function (VPF), we investigate the statistical properties of voids within a wide range of initially Gaussian models for the origin of large-scale structure. We pay particular attention to the ability of the VPF to diagnose `biased galaxy formation': the preferential formation of galaxies in regions of high background density and corresponding suppression of galaxy formation in regions of low background density. Biased models have systematically higher VPFs than unbiased models, but for a given biasing scheme the VPF is relatively insensitive to the value of the bias factor $b$. Thus, while the VPF can distinguish unbiased models from some biased models, it is probably not a useful way to constrain the bias factor; uncertainties in the appropriate choice of biasing prescription overwhelm the mild dependence on $b$. We compare the predictions of our models to the most extensive VPF observations published to date. These data do {\it not} require strong biasing; Gaussian models in which galaxies trace mass can reproduce the VPF data to within the errors expected from the current finite volume fluctuations. Models with the moderate biasing predicted by cosmological simulations with gas dynamics yield a slightly better match to the data. Models in which galaxy formation is strongly suppressed in low density regions produce an excess of large, empty voids.

preprint1993arXiv

Fourier Analysis of Redshift Space Distortions and the Determination of Omega

The peculiar velocities of galaxies distort the pattern of galaxy clustering in redshift space, making the redshift space power spectrum anisotropic. In the linear regime, the strength of this distortion depends only on the ratio $β\equiv f(Ω)/b \approx Ω^{0.6}/b$, where $Ω$ is the cosmological density parameter and $b$ is the bias parameter. We derive a linear theory estimator for $β$ based on the harmonic moments of the redshift space power spectrum. Using N-body simulations, we examine the impact of non-linear gravitational clustering on the power spectrum anisotropy and on our $β$-estimator. Non-linear effects can be important out to wavelengths $λ\sim 50$Mpc/h or larger; in most cases, they lower the quadrupole moment of the power spectrum and thereby depress the estimate of $β$ below the true value. With a sufficiently large redshift survey, the scaling of non-linear effects may allow separate determinations of $Ω$ and $b$. We describe a practical technique for measuring the anisotropy of the power spectrum from galaxy redshift surveys, and we test the technique on mock catalogues drawn from the N-body simulations. Preliminary application of our methods to the 1.2 Jy IRAS galaxy survey yields $β_{est} \sim 0.3-0.4 $ at wavelengths $λ\sim 30-40$Mpc/h . Non-linear effects remain important at these scales, so this estimate of $β$ is probably lower than the true value.