Mock galaxy catalogs using the quick particle mesh method

Sophisticated analysis of modern large-scale structure surveys requires mock catalogs. Mock catalogs are used to optimize survey design, test reduction and analysis pipelines, make theoretical predictions for basic observables and propagate errors through complex analysis chains. We present a new method, which we call "quick particle mesh", for generating many large-volume, approximate mock catalogs at low computational cost. The method is based on using rapid, low-resolution particle mesh simulations that accurately reproduce the large-scale dark matter density field. Particles are sampled from the density field based on their local density such that they have N-point statistics nearly equivalent to the halos resolved in high-resolution simulations, creating a set of mock halos that can be populated using halo occupation methods to create galaxy mocks for a variety of possible target classes.Comment: 13 pages, 16 figures. Matches version accepted by MNRAS. Code available at http://github.com/mockFactor

A 2.5% measurement of the growth rate from small-scale redshift space clustering of SDSS-III CMASS galaxies

We perform the first fit to the anisotropic clustering of SDSS-III CMASS DR10 galaxies on scales of ~ 0.8 - 32 Mpc/h. A standard halo occupation distribution model evaluated near the best fit Planck LCDM cosmology provides a good fit to the observed anisotropic clustering, and implies a normalization for the peculiar velocity field of M ~ 2 x 10^13 Msun/h halos of f*sigma8(z=0.57) = 0.450 +/- 0.011. Since this constraint includes both quasi-linear and non-linear scales, it should severely constrain modified gravity models that enhance pairwise infall velocities on these scales. Though model dependent, our measurement represents a factor of 2.5 improvement in precision over the analysis of DR11 on large scales, f*sigma8(z=0.57) = 0.447 +/- 0.028, and is the tightest single constraint on the growth rate of cosmic structure to date. Our measurement is consistent with the Planck LCDM prediction of 0.480 +/- 0.010 at the ~1.9 sigma level. Assuming a halo mass function evaluated at the best fit Planck cosmology, we also find that 10% of CMASS galaxies are satellites in halos of mass M ~ 6 x 10^13 Msun/h. While none of our tests and model generalizations indicate systematic errors due to an insufficiently detailed model of the galaxy-halo connection, the precision of these first results warrant further investigation into the modeling uncertainties and degeneracies with cosmological parameters.Comment: 24 pages, 20 figures, submitted to MNRAS. v2 is 27 pages, 23 figures, accepted by MNRA

What Determines the Incidence and Extent of MgII Absorbing Gas Around Galaxies?

We study the connections between on-going star formation, galaxy mass, and extended halo gas, in order to distinguish between starburst-driven outflows and infalling clouds that produce the majority of observed MgII absorbers at large galactic radii (>~ 10 h^{-1} kpc) and to gain insights into halo gas contents around galaxies. We present new measurements of total stellar mass (M_star), H-alpha emission line strength (EW(H-alpha)), and specific star formation rate (sSFR) for the 94 galaxies published in H.-W. Chen et al. (2010). We find that the extent of MgII absorbing gas, R_MgII, scales with M_star and sSFR, following R_MgII \propto M_star^{0.28}\times sSFR^{0.11}. The strong dependence of R_MgII on M_star is most naturally explained, if more massive galaxies possess more extended halos of cool gas and the observed MgII absorbers arise in infalling clouds which will subsequently fuel star formation in the galaxies. The additional scaling relation of R_MgII with sSFR can be understood either as accounting for extra gas supplies due to starburst outflows or as correcting for suppressed cool gas content in high-mass halos. The latter is motivated by the well-known sSFR--M_star} inverse correlation in field galaxies. Our analysis shows that a joint study of galaxies and MgII absorbers along common sightlines provides an empirical characterization of halo gaseous radius versus halo mass. A comparison study of R_MgII around red- and blue-sequence galaxies may provide the first empirical constraint for resolving the physical origin of the observed sSFR--M_star} relation in galaxies.Comment: 6 pages, 3 figures; ApJL in pres

Tests of redshift-space distortions models in configuration space for the analysis of the BOSS final data release

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for observing the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. In preparation for analysis of redshift-space distortions from the Baryon Oscillation Spectroscopic Survey (BOSS) final data release we compare a number of analytic and phenomenological `streaming' models, specified in configuration space, to mock catalogs derived in different ways from several N-body simulations. The galaxies in each mock catalog have properties similar to those of the higher redshift galaxies measured by BOSS but differ in the details of how small-scale velocities and halo occupancy are determined. We find that all of the analytic models fit the simulations over a limited range of scales while failing at small scales. We discuss which models are most robust and on which scales they return reliable estimates of the rate of growth of structure: we find that models based on some form of resummation can fit our N-body data for BOSS-like galaxies above $30\,h^{-1}$Mpc well enough to return unbiased parameter estimates.Comment: 12 pages, 11 figures, matches version accepted by MNRA

Halo Occupation Distribution Modeling of Green Valley Galaxies

We present a clustering analysis of near ultraviolet (NUV) - optical color selected luminosity bin samples of green valley galaxies. These galaxy samples are constructed by matching the Sloan Digital Sky Survey Data Release 7 with the latest Galaxy Evolution Explorer source catalog which provides NUV photometry. We present cross-correlation function measurements and determine the halo occupation distribution of these transitional galaxies using a new multiple tracer analysis technique. We extend the halo-occupation formalism to model the cross-correlation function between a galaxy sample of interest and multiple tracer populations simultaneously. This method can be applied to commonly used luminosity threshold samples as well as to color and luminosity bin selected galaxy samples, and improves the accuracy of clustering analyses for sparse galaxy populations. We confirm the previously observed trend that red galaxies reside in more massive halos and are more likely to be satellite galaxies than average galaxies of similar luminosity. While the change in central galaxy host mass as a function of color is only weakly constrained, the satellite fraction and characteristic halo masses of green satellite galaxies are found to be intermediate between those of blue and red satellite galaxies.Comment: matches MNRAS accepted version; minor revisions, results unchange

Mass Function Predictions Beyond LCDM

The mass distribution of halos, as specified by the halo mass function, is a key input for several cosmological probes. The sizes of $N$-body simulations are now such that, for the most part, results need no longer be statistics-limited, but are still subject to various systematic uncertainties. We investigate and discuss some of the reasons for these differences. Quantifying error sources and compensating for them as appropriate, we carry out a high-statistics study of dark matter halos from 67 $N$-body simulations to investigate the mass function and its evolution for a reference $\Lambda$CDM cosmology and for a set of $w$CDM cosmologies. For the reference $\Lambda$CDM cosmology (close to WMAP5), we quantify the breaking of universality in the form of the mass function as a function of redshift, finding an evolution of as much as 10% away from the universal form between redshifts $z=0$ and $z=2$. For cosmologies very close to this reference we provide a fitting formula to our results for the (evolving) $\Lambda$CDM mass function over a mass range of $6\cdot 10^{11}-3\cdot 10^{15}$ M$_{\odot}$ to an estimated accuracy of about 2%. The set of $w$CDM cosmologies is taken from the Coyote Universe simulation suite. The mass functions from this suite (which includes a $\Lambda$CDM cosmology and others with $w\simeq-1$) are described by the fitting formula for the reference $\Lambda$CDM case at an accuracy level of 10%, but with clear systematic deviations. We argue that, as a consequence, fitting formulae based on a universal form for the mass function may have limited utility in high precision cosmological applications.Comment: 19 pages; 18 figures; accepted for publication in the Ap

Reductions in the dietary niche of southern sea otters (Enhydra lutris nereis) from the Holocene to the Anthropocene.

The sea otter (Enhydra lutris) is a marine mammal hunted to near extinction during the 1800s. Despite their well-known modern importance as a keystone species, we know little about historical sea otter ecology. Here, we characterize the ecological niche of ancient southern sea otters (E. lutris nereis) using δ13C analysis and δ15N analysis of bones recovered from archaeological sites spanning ~7,000 to 350 years before present (N = 112 individuals) at five regions along the coast of California. These data are compared with previously published data on modern animals (N = 165) and potential modern prey items. In addition, we analyze the δ15N of individual amino acids for 23 individuals to test for differences in sea otter trophic ecology through time. After correcting for tissue-specific and temporal isotopic effects, we employ nonparametric statistics and Bayesian niche models to quantify differences among ancient and modern animals. We find ancient otters occupied a larger isotopic niche than nearly all modern localities; likely reflecting broader habitat and prey use in prefur trade populations. In addition, ancient sea otters at the most southerly sites occupied an isotopic niche that was more than twice as large as ancient otters from northerly regions. This likely reflects greater invertebrate prey diversity in southern California relative to northern California. Thus, we suggest the potential dietary niche of sea otters in southern California could be larger than in central and northern California. At two sites, Año Nuevo and Monterey Bay, ancient otters had significantly higher δ15N values than modern populations. Amino acid δ15N data indicated this resulted from shifting baseline isotope values, rather than a change in sea otter trophic ecology. Our results help in better understanding the contemporary ecological role of sea otters and exemplify the strength of combing zooarchaeological and biological information to provide baseline data for conservation efforts

Connecting massive galaxies to dark matter halos in BOSS - I. Is galaxy color a stochastic process in high-mass halos?

We use subhalo abundance matching (SHAM) to model the stellar mass function (SMF) and clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) "CMASS" sample at $z\sim0.5$. We introduce a novel method which accounts for the stellar mass incompleteness of CMASS as a function of redshift, and produce CMASS mock catalogs which include selection effects, reproduce the overall SMF, the projected two-point correlation function $w_{\rm p}$, the CMASS $dn/dz$, and are made publicly available. We study the effects of assembly bias above collapse mass in the context of "age matching" and show that these effects are markedly different compared to the ones explored by Hearin et al. (2013) at lower stellar masses. We construct two models, one in which galaxy color is stochastic ("AbM" model) as well as a model which contains assembly bias effects ("AgM" model). By confronting the redshift dependent clustering of CMASS with the predictions from our model, we argue that that galaxy colors are not a stochastic process in high-mass halos. Our results suggest that the colors of galaxies in high-mass halos are determined by other halo properties besides halo peak velocity and that assembly bias effects play an important role in determining the clustering properties of this sample.Comment: 22 pages. Appendix. B added. Matches the version accepted by MNRAS. Mock galaxy catalog and HOD table are available at http://www.massivegalaxies.co

A Simple Model for Quasar Demographics

We present a simple model for the relationship between quasars, galaxies, and dark matter halos from 0.5<z<6. In the model, black hole (BH) mass is linearly related to galaxy mass, and galaxies are connected to dark matter halos via empirically constrained relations. A simple "scattered" light bulb model for quasars is adopted, wherein BHs shine at a fixed fraction of the Eddington luminosity during accretion episodes, and Eddington ratios are drawn from a lognormal distribution that is redshift-independent. This model has two free, physically meaningful parameters at each redshift: the normalization of the Mbh-Mgal relation and the quasar duty cycle; these parameters are fit to the observed quasar luminosity function (LF) over the interval 0.5<z<6. This simple model provides an excellent fit to the LF at all epochs, and also successfully predicts the observed projected two-point correlation of quasars from 0.5<z<2.5. It is significant that a single quasar duty cycle at each redshift is capable of reproducing the extant observations. The data are therefore consistent with a scenario wherein quasars are equally likely to exist in galaxies, and therefore dark matter halos, over a wide range in masses. The knee in the quasar LF is a reflection of the knee in the stellar mass-halo mass relation. Future constraints on the quasar LF and quasar clustering at high redshift will provide strong constraints on the model. In the model, the autocorrelation function of quasars becomes a strong function of luminosity only at the very highest luminosities, and will be difficult to observe because such quasars are so rare. Cross-correlation techniques may provide useful constraints on the bias of such rare objects.Comment: 12 pages, 12 figures, ApJ accepte

Towards an accurate model of the redshift space clustering of halos in the quasilinear regime

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for measuring the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. Galaxies occupy dark matter halos, whose redshift space clustering has a complex dependence on bias that cannot be inferred from the behavior of matter. We identify two distinct corrections on quasilinear scales (~ 30-80 Mpc/h): the non-linear mapping between real and redshift space positions, and the non-linear suppression of power in the velocity divergence field. We model the first non-perturbatively using the scale-dependent Gaussian streaming model, which we show is accurate at the <0.5 (2) per cent level in transforming real space clustering and velocity statistics into redshift space on scales s>10 (s>25) Mpc/h for the monopole (quadrupole) halo correlation functions. We use perturbation theory to predict the real space pairwise halo velocity statistics. Our fully analytic model is accurate at the 2 per cent level only on scales s > 40 Mpc/h. Recent models that neglect the corrections from the bispectrum and higher order terms from the non-linear real-to-redshift space mapping will not have the accuracy required for current and future observational analyses. Finally, we note that our simulation results confirm the essential but non-trivial assumption that on large scales, the bias inferred from real space clustering of halos is the same one that determines their pairwise infall velocity amplitude at the per cent level.Comment: 15 pages, 12 figures, submitted to MNRA