The James Webb Space Telescope (JWST) is expected to observe galaxies at $z>10$ that are presently inaccessible. Here, we use a self-consistent empirical model, the UniverseMachine, to generate mock galaxy catalogues and lightcones over the redshift range $z=0-15$. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy-halo relationships, and galaxy-galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at $z<10$. We predict that Cycle 1 JWST surveys will very likely detect galaxies with $M_*>10^7 M_\odot$ and/or $M_{1500}<-17$ out to at least $z\sim 13.5$. Number density uncertainties at $z>12$ expand dramatically, so efforts to detect $z>12$ galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits ($M_\ast < 10^7M_\odot$ or $M_\mathrm{1500}>-17$). For reionization models, extrapolating luminosity functions with a constant faint-end slope from $M_{1500}=-17$ down to $M_{1500}=-12$ gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to $z\sim 12$. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and lightcones for common fields are available online.

The shape of dark matter haloes plays a critical role in constraining cosmology with upcoming large-scale structure surveys. In this paper, we study the correlations between the triaxial shapes and formation histories in dark matter haloes in the MultiDark Planck 2 N-body cosmological simulation. We find that halo ellipticity is strongly correlated with halo properties that serve as proxies of halo formation history, such as halo concentration and the peak-centroid offset. In particular, the correlation between halo ellipticity and halo concentration is nearly independent of the halo density peak height. We present a simple model for the correlation between halo ellipticity and concentration using conditional abundance matching, and provide fitting formulae for the multi-dimensional distributions of triaxial halo shape as a function of halo peak height. We apply our halo shape model to gauge the effects of halo ellipticity and orientation bias on the excess surface mass density profiles in cluster-size halos. Our model should be useful for exploring the impact of triaxial halo shape on cosmological constraints in upcoming weak lensing surveys of galaxy clusters.

The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA's Euclid, and NASA's Wide Field Infrared Survey Telescope (WFIRST). The Co-Chairs of TACS, Katrin Heitmann and Alina Kiessling, invited community scientists from the USA and Europe who are each subject matter experts and are also members of one or more of the surveys to contribute. The following report represents the input from TACS that was delivered to the Agencies in December 2018.

The spectral energy distribution of a galaxy emerges from the complex interplay of many physical ingredients, including its star formation history (SFH), metallicity evolution, and dust properties. Using Galaxpy, a new stellar population synthesis model, and SFHs predicted by the empirical model UniverseMachine and the cosmological hydrodynamical simulation IllustrisTNG, we isolate the influence of SFH on broad-band colours at z=0. By carrying out a principal component analysis, we show that physically-motivated SFH variations modify broad-band colours along a single direction in colour space: the SFH-direction. We find that the projection of a galaxy's present-day colours onto the SFH-direction is almost completely regulated by the fraction of stellar mass that the galaxy formed over the last billion years. Together with cosmic downsizing, this results in galaxies becoming redder as their host halo mass increases. We additionally study the change in galaxy colours due to variations in metallicity, dust attenuation, and nebular emission lines, finding that these properties vary broad-band colours along distinct directions in colour space relative to the SFH-direction. Finally, we show that the colours of low-redshift SDSS galaxies span an ellipsoid with significant extent along two independent dimensions, and that the SFH-direction is well-aligned with the major axis of this ellipsoid. Our analysis supports the conclusion that variations in star formation history are the dominant influence on present-day galaxy colours, and that the nature of this influence is strikingly simple.

As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-analytic models as physically-motivated galaxy libraries. GalSampler populates a new, larger-volume halo catalog with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modeling techniques, GalSampler inherits statistical accuracy from the empirical model and physically-motivated complexity from the baseline library. We have recently used GalSampler to produce the cosmoDC2 extragalactic catalog made for the LSST DESC Data Challenge 2. Using cosmoDC2 as a guiding example, we outline how GalSampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), WFIRST, and Euclid.

Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $\Delta \Sigma(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+\Delta \Sigma(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $\Delta \Sigma(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection.

Extracting accurate cosmological information from galaxy-galaxy and galaxy-matter correlation functions on non-linear scales ($\lesssim 10 h^{-1} \mathrm{Mpc}$) requires cosmological simulations. Additionally, one has to marginalise over several nuisance parameters of the galaxy-halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $\mathcal{O}(10^6)$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realisations in parameter space. Afterwards, this emulator is used as a surrogate model in an MCMC analysis. Here, we demonstrate a different method called Cosmological Evidence Modelling (CEM). First, for each simulation, we calculate the Bayesian evidence marginalised over the galaxy-halo connection by repeatedly populating the simulation with galaxies. We show that this Bayesian evidence is directly related to the posterior probability of cosmological parameters. Finally, we build a physically motivated model for how the evidence depends on cosmological parameters as sampled by the simulations. We demonstrate the feasibility of CEM by using simulations from the Aemulus simulation suite and forecasting cosmological constraints from BOSS CMASS measurements of redshift-space distortions. Our analysis includes an exploration of how galaxy assembly bias affects cosmological inference. Overall, CEM has several potential advantages over the more common approach of emulating summary statistics, including the ability to easily marginalise over highly complex models of the galaxy-halo connection and greater accuracy, thereby reducing the number of simulations required.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

We present a method to flexibly and self-consistently determine individual galaxies' star formation rates (SFRs) from their host haloes' potential well depths, assembly histories, and redshifts. The method is constrained by galaxies' observed stellar mass functions, SFRs (specific and cosmic), quenched fractions, UV luminosity functions, UV-SM relations, IRX-UV relations, auto- and cross-correlation functions (including quenched and star-forming subsamples), and quenching dependence on environment; each observable is reproduced over the full redshift range available, up to 0<z<10. Key findings include: galaxy assembly correlates strongly with halo assembly; quenching at z>1 correlates strongly with halo mass; quenched fractions at fixed halo mass decrease with increasing redshift; massive quenched galaxies reside in higher-mass haloes than star-forming galaxies at fixed galaxy mass; star-forming and quenched galaxies' star formation histories at fixed mass differ most at z<0.5; satellites have large scatter in quenching timescales after infall, and have modestly higher quenched fractions than central galaxies; Planck cosmologies result in up to 0.3 dex lower stellar mass-halo mass ratios at early times; and, nonetheless, stellar mass-halo mass ratios rise at z>5. Also presented are revised stellar mass-halo mass relations for all, quenched, star-forming, central, and satellite galaxies; the dependence of star formation histories on halo mass, stellar mass, and galaxy SSFR; quenched fractions and quenching timescale distributions for satellites; and predictions for higher-redshift galaxy correlation functions and weak lensing surface densities. The public data release (DR1) includes the massively parallel (>10^5 cores) implementation (the UniverseMachine), the newly compiled and remeasured observational data, derived galaxy formation constraints, and mock catalogues including lightcones.

Previous work has demonstrated that at a given stellar mass, quiescent galaxies are more strongly clustered than star-forming galaxies. The contribution to this signal from central, as opposed to satellite, galaxies is not known, which has strong implications for galaxy evolution models. To investigate the contribution from central galaxies, here we present measurements of the clustering of isolated primary (IP) galaxies, used as a proxy for central galaxies, at ${0.2<z<0.9}$ with data from the PRIMUS galaxy redshift survey. Using a sample of spectroscopic redshifts for $\sim60,000$ galaxies with ${M_*\gtrsim10^9~M_\odot}$ covering 5 deg$^2$ on the sky, we define IP galaxies using isolation cuts in spatial proximity and stellar mass of nearby galaxies. We find that at fixed stellar mass, quiescent IP galaxies are more strongly clustered than star-forming IP galaxies at $z\sim0.35$ $(10\sigma)$. Using mock galaxy catalogs based on recent halo occupation models of Behroozi et al. (2018) and designed to replicate the parameters of the PRIMUS survey dataset, we find that these clustering differences are due in part to quiescent central galaxies being more strongly clustered than star-forming central galaxies. This is consistent with either distinct stellar-to-halo mass relations for quiescent and star-forming central galaxies, and/or central galaxy assembly bias. We additionally use mock catalogs to assess the dependence of both incompleteness and satellite galaxy contamination in the IP galaxy samples on redshift, galaxy type, and stellar mass, and demonstrate how isolation criteria yield biased subsamples of central galaxies via {\it environmental incompleteness}, or the preferential exclusion of central galaxies in overdense environments.

We use the UniverseMachine to analyze the source of scatter between the central galaxy mass, the total stellar mass in the halo, and the dark matter halo mass. We also propose a new halo mass estimator, the cen+N mass: the sum of the stellar mass of the central and the N most massive satellites. We show that, when real space positions are perfectly known, the cen+N mass has scatter competitive with that of richness-based estimators. However, in redshift space, the cen+N mass suffers less from projection effects in the UniverseMachine model. The cen+N mass is therefore a viable low scatter halo mass estimator, and should be considered an important tool to constrain cosmology with upcoming spectroscopic data from DESI. We analyze the scatter in stellar mass at fixed halo mass and show that the total stellar mass in a halo is uncorrelated with secondary halo properties, but that the central stellar mass is a function of both halo mass and halo age. This is because central galaxies in older halos have had more time to grow via accretion. If the UniverseMachine model is correct, accurate galaxy-halo modeling of mass selected samples therefore needs to consider halo age in addition to mass.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Using deep images from the Hyper Suprime-Cam (HSC) survey and taking advantage of its unprecedented weak lensing capabilities, we reveal a remarkably tight connection between the stellar mass distribution of massive central galaxies and their host dark matter halo mass. Massive galaxies with more extended stellar mass distributions tend to live in more massive dark matter haloes. We explain this connection with a phenomenological model that assumes, (1) a tight relation between the halo mass and the total stellar content in the halo, (2) that the fraction of in-situ and ex-situ mass at $r<10$ kpc depends on halo mass. This model provides an excellent description of the stellar mass functions (SMF) of total stellar mass ($M_{\star}^{\rm Max}$) and stellar mass within inner 10 kpc ($M_{\star}^{10}$) and also reproduces the HSC weak lensing signals of massive galaxies with different stellar mass distributions. The best-fit model shows that halo mass varies significantly at fixed total stellar mass (as much as 0.4 dex) with a clear dependence on $M_{\star}^{10}$. Our two-parameter $M_{\star}^{\rm Max}$-$M_{\star}^{10}$ description provides a more accurate picture of the galaxy-halo connection at the high-mass end than the simple stellar-halo mass relation (SHMR) and opens a new window to connect the assembly history of halos with those of central galaxies. The model also predicts that the ex-situ component dominates the mass profiles of galaxies at $r< 10$ kpc for $\log M_{\star} \ge 11.7$). The code used for this paper is available online: https://github.com/dr-guangtou/asap

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

We empirically constrain how galaxy size relates to halo virial radius using new measurements of the size- and stellar mass-dependent clustering of galaxies in the Sloan Digital Sky Survey. We find that small galaxies cluster much more strongly than large galaxies of the same stellar mass. The magnitude of this clustering difference increases on small scales, and decreases with increasing stellar mass. Using Halotools to forward model the observations, we test an empirical model in which present-day galaxy size is proportional to the size of the virial radius at the time the halo reached its maximum mass. This simple model reproduces the observed size-dependence of galaxy clustering in striking detail. The success of this model provides strong support for the conclusion that satellite galaxies have smaller sizes relative to central galaxies of the same halo mass. Our findings indicate that satellite size is set prior to the time of infall, and that a remarkably simple, linear size--virial radius relation emerges from the complex physics regulating galaxy size. We make quantitative predictions for future measurements of galaxy-galaxy lensing, including dependence upon size, scale, and stellar mass, and provide a scaling relation of the ratio of mean sizes of satellites and central galaxies as a function of their halo mass that can be used to calibrate hydrodynamical simulations and semi-analytic models.

Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight velocity and spatial offsets between brightest galaxies and their neighbours. We compare those to detailed mock catalogues, constructed from high-resolution, dark-matter-only $N$-body simulations, in which it is assumed that satellite galaxies trace dark matter subhaloes. This allows us to place constraints on the fraction $f_{\rm BNC}$ of haloes in which the brightest galaxy is not the central. Compared to previous studies we explicitly take into account the unrelaxed state of the host haloes, velocity offsets of halo cores and correlations between $f_{\rm BNC}$ and the satellite occupation. We find that $f_{\rm BNC}$ strongly decreases with the luminosity of the brightest galaxy and increases with the mass of the host halo. Overall, in the halo mass range $10^{13} - 10^{14.5} h^{-1} M_\odot$ we find $f_{\rm BNC} \sim 30\%$, in good agreement with a previous study by Skibba et al. We discuss the implications of these findings for studies inferring the galaxy--halo connection from satellite kinematics, models of the conditional luminosity function and galaxy formation in general.

We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python.

Standard approaches to Bayesian parameter inference in large scale structure assume a Gaussian functional form (chi-squared form) for the likelihood. This assumption, in detail, cannot be correct. Likelihood free inferences such as Approximate Bayesian Computation (ABC) relax these restrictions and make inference possible without making any assumptions on the likelihood. Instead ABC relies on a forward generative model of the data and a metric for measuring the distance between the model and data. In this work, we demonstrate that ABC is feasible for LSS parameter inference by using it to constrain parameters of the halo occupation distribution (HOD) model for populating dark matter halos with galaxies. Using specific implementation of ABC supplemented with Population Monte Carlo importance sampling, a generative forward model using HOD, and a distance metric based on galaxy number density, two-point correlation function, and galaxy group multiplicity function, we constrain the HOD parameters of mock observation generated from selected "true" HOD parameters. The parameter constraints we obtain from ABC are consistent with the "true" HOD parameters, demonstrating that ABC can be reliably used for parameter inference in LSS. Furthermore, we compare our ABC constraints to constraints we obtain using a pseudo-likelihood function of Gaussian form with MCMC and find consistent HOD parameter constraints. Ultimately our results suggest that ABC can and should be applied in parameter inference for LSS analyses.

We constrain the newly-introduced decorated Halo Occupation Distribution (HOD) model using SDSS DR7 measurements of projected galaxy clustering or r-band luminosity threshold samples. The decorated HOD is a model for the galaxy-halo connection that augments the HOD by allowing for the possibility of galaxy assembly bias: galaxy luminosity may be correlated with dark matter halo properties besides mass, Mvir. We demonstrate that it is not possible to rule out galaxy assembly bias using DR7 measurements of galaxy clustering alone. Moreover, galaxy samples with Mr < -20 and Mr < -20.5 favor strong central galaxy assembly bias. These samples prefer scenarios in which high-concentration are more likely to host a central galaxy relative to low-concentration halos of the same mass. We exclude zero assembly bias with high significance for these samples. Satellite galaxy assembly bias is significant for the faintest sample, Mr < -19. We find no evidence for assembly bias in the Mr < -21 sample. Assembly bias should be accounted for in galaxy clustering analyses or attempts to exploit galaxy clustering to constrain cosmology. In addition to presenting the first constraints on HOD models that accommodate assembly bias, our analysis includes several improvements over previous analyses of these data. Therefore, our inferences supersede previously-published results even in the case of a standard HOD analysis.

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.

The connection between galaxies and dark matter halos is often inferred from data using probabilistic models, such as the Halo Occupation Distribution (HOD). Conventional HOD formulations assume that only halo mass governs the galaxy-halo connection. Violations of this assumption, known as galaxy assembly bias, threaten the HOD program. We introduce decorated HODs, a new, flexible class of models designed to account for assembly bias. Decorated HODs minimally expand the parameter space and maximize the independence between traditional and novel HOD parameters. We use decorated HODs to quantify the influence of assembly bias on clustering and lensing statistics. For SDSS-like samples, the impact of assembly bias on galaxy clustering can be as large as a factor of two on r ~ 200 kpc scales and ~15% in the linear regime. Assembly bias can either enhance or diminish clustering on large scales, but generally increases clustering on scales r <~ 1 Mpc. We performed our calculations with Halotools, an open-source, community-driven python package for studying the galaxy-halo connection (http://halotools.readthedocs.org). We conclude by describing the use of decorated HODs to treat assembly bias in otherwise conventional likelihood analyses.

The two-point clustering of dark matter halos is influenced by halo properties besides mass, a phenomenon referred to as halo assembly bias. Using the depth of the gravitational potential well, $V_{\rm max}$, as our secondary halo property, in this paper we present the first study of the scale-dependence assembly bias. In the large-scale linear regime, $r\geq10h^{-1}{\rm Mpc},$ our findings are in keeping with previous results. In particular, at the low-mass end ($M_{\rm vir}<M_{\rm coll}\approx10^{12.5}{\rm M}_{\odot}$), halos with high-$V_{\rm max}$ show stronger large-scale clustering relative to halos with low-$V_{\rm max}$ of the same mass, this trend weakens and reverses for $M_{\rm vir}\geq M_{\rm coll}.$ In the nonlinear regime, assembly bias in low-mass halos exhibits a pronounced scale-dependent "bump" at $500h^{-1}{\rm kpc}-5h^{-1}{\rm Mpc},$ a new result. This feature weakens and eventually vanishes for halos of higher mass. We show that this scale-dependent signature can primarily be attributed to a special subpopulation of ejected halos, defined as present-day host halos that were previously members of a higher-mass halo at some point in their past history. A corollary of our results is that galaxy clustering on scales of $r\sim1-2h^{-1}{\rm Mpc}$ can be impacted by up to $\sim15\%$ by the choice of the halo property used in the halo model, even for stellar mass-limited samples.

Currently-proposed galaxy quenching mechanisms predict very different behaviours during major halo mergers, ranging from significant quenching enhancement (e.g., clump-induced gravitational heating models) to significant star formation enhancement (e.g., gas starvation models). To test real galaxies' behaviour, we present an observational galaxy pair method for selecting galaxies whose host haloes are preferentially undergoing major mergers. Applying the method to central L* (10^10 Msun < M_* < 10^10.5 Msun) galaxies in the Sloan Digital Sky Survey (SDSS) at z<0.06, we find that major halo mergers can at most modestly reduce the star-forming fraction, from 59% to 47%. Consistent with past research, however, mergers accompany enhanced specific star formation rates for star-forming L* centrals: ~10% when a paired galaxy is within 200 kpc (approximately the host halo's virial radius), climbing to ~70% when a paired galaxy is within 30 kpc. No evidence is seen for even extremely close pairs (<30 kpc separation) rejuvenating star formation in quenched galaxies. For galaxy formation models, our results suggest: (1) quenching in L* galaxies likely begins due to decoupling of the galaxy from existing hot and cold gas reservoirs, rather than a lack of available gas or gravitational heating from infalling clumps, (2) state-of-the-art semi-analytic models currently over-predict the effect of major halo mergers on quenching, and (3) major halo mergers can trigger enhanced star formation in non-quenched central galaxies.

We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most colour-dependent statistics are recovered with reasonable accuracy. However, it is also clear that certain systematic errors arise as a consequence of correlated errors in group membership determination, central/satellite designation, and halo mass assignment. We introduce a new statistic, the halo transition probability (HTP), which captures the combined impact of all these errors. As a rule of thumb, errors tend to equalize the properties of distinct galaxy populations (i.e. red vs. blue galaxies or centrals vs. satellites), and to result in inferred occupation statistics that are more accurate for red galaxies than for blue galaxies. A statistic that is particularly poorly recovered from the group catalogues is the red fraction of central galaxies as function of halo mass. Group finders do a good job in recovering galactic conformity, but also have a tendency to introduce weak conformity when none is present. We conclude that proper inference of colour-dependent statistics from group catalogues is best achieved using forward modelling (i.e., running group finders over mock data), or by implementing a correction scheme based on the HTP, as long as the latter is not too strongly model-dependent.

Correlations between the star formation rates (SFRs) of nearby galaxies (so-called galactic conformity) have been observed for projected separations up to 4 Mpc, an effect not predicted by current semi-analytic models. We investigate correlations between the mass accretion rates (dMvir/dt) of nearby halos as a potential physical origin for this effect. We find that pairs of host halos "know about" each others' assembly histories even when their present-day separation is greater than thirty times the virial radius of either halo. These distances are far too large for direct interaction between the halos to explain the correlation in their dMvir/dt. Instead, halo pairs at these distances reside in the same large-scale tidal environment, which regulates dMvir/dt for both halos. Larger halos are less affected by external forces, which naturally gives rise to a mass dependence of the halo conformity signal. SDSS measurements of galactic conformity exhibit a qualitatively similar dependence on stellar mass, including how the signal varies with distance. Based on the expectation that halo accretion and galaxy SFR are correlated, we predict the scale-, mass- and redshift-dependence of large-scale galactic conformity, finding that the signal should drop to undetectable levels by z > 1. These predictions are testable with current surveys to z ~ 1; confirmation would establish a strong correlation between dark matter halo accretion rate and central galaxy SFR.

Among the most stringent constraints on the dark matter annihilation cross section are those derived from observations of dwarf galaxies by the Fermi Gamma-Ray Space Telescope. As current (e.g., Dark Energy Survey, DES) and future (Large Synoptic Survey Telescope, LSST) optical imaging surveys discover more of the Milky Way's ultra-faint satellite galaxies, they may increase Fermi's sensitivity to dark matter annihilations. In this study, we use a semi-analytic model of the Milky Way's satellite population to predict the characteristics of the dwarfs likely to be discovered by DES and LSST, and project how these discoveries will impact Fermi's sensitivity to dark matter. While we find that modest improvements are likely, the dwarf galaxies discovered by DES and LSST are unlikely to increase Fermi's sensitivity by more than a factor of ~2-4. However, this outlook may be conservative, given that our model underpredicts the number of ultra-faint galaxies with large potential annihilation signals actually discovered in the Sloan Digital Sky Survey. Our simulation-based approach focusing on the Milky Way satellite population demographics complements existing empirically-based estimates.

It is common practice for methods that use galaxy clustering to constrain the galaxy-halo relationship, such as the halo occupation distribution (HOD) and/or conditional luminosity function (CLF), to assume that halo mass alone suffices to determine a halo's resident galaxy population. Yet, the clustering strength of cold dark matter halos depends upon halo properties in addition to mass, such as formation time, an effect referred to as assembly bias. If galaxy characteristics are correlated with any of these auxiliary halo properties, the basic assumption of HOD/CLF methods is violated. We estimate the potential for assembly bias to induce systematic errors in inferred halo occupation statistics. We use halo abundance matching and age matching to construct fiducial mock galaxy catalogs that exhibit assembly bias as well as additional mock catalogs with identical HODs, but with assembly bias removed. We fit a parameterized HOD to the projected two-point clustering of mock galaxies in each catalog to assess the systematic errors induced by reasonable levels of assembly bias. In the absence of assembly bias, the inferred HODs generally describe the true underlying HODs well, validating the basic methodology. However, in all of the cases with assembly bias, the inferred HODs have systematic errors that are statistically significant. In most cases, these systematic errors cannot be identified using void statistics as auxiliary observables. We conclude that the galaxy-halo relationship inferred from galaxy clustering should be subject to a non-negligible systematic error induced by assembly bias. Our work suggests that efforts to model and/or constrain assembly bias should be high priorities as it is a threatening source of systematic error in galaxy evolution studies as well as the precision cosmology program.

The redshift-space distortion (RSD) of galaxies surrounding massive clusters is emerging as a promising testbed for theories of modified gravity. Conventional applications of this method rely upon the assumption that the velocity field in the cluster environment is uniquely determined by the cluster mass profile. Yet, real dark matter halos in N-body simulations are known to violate the assumption that virial mass determines the configuration space distribution, an effect known as assembly bias. In this Letter, I show that assembly bias in simulated dark matter halos also manifests in velocity space. In the 1-10 Mpc environment surrounding a cluster, high-concentration "tracer" halos exhibit a 10-20% larger pairwise-velocity dispersion profile relative to low-concentration tracer halos of the same mass. This difference is comparable to the size of the RSD signal predicted by f(R) models designed to account for the cosmic acceleration. I use the age matching technique to study how color-selection effects may influence the cluster RSD signal, finding a ~10% effect due to redder satellites preferentially occupying higher mass halos, and a ~5% effect due to assembly-biased colors of centrals. In order to use cluster RSD measurements to robustly constrain modified gravity, we likely will need to develop empirical galaxy formation models more sophisticated than any in the current literature.

In this paper, we test the age matching hypothesis that the star formation rate (SFR) of a galaxy of fixed stellar mass is determined by its dark matter halo formation history, and as such, that more quiescent galaxies reside in older halos. This simple model has been remarkably successful at predicting color-based galaxy statistics at low redshift as measured in the Sloan Digital Sky Survey (SDSS). To further test this method with observations, we present new SDSS measurements of the galaxy two-point correlation function and galaxy-galaxy lensing as a function of stellar mass and SFR, separated into quenched and star-forming galaxy samples. We find that our age matching model is in excellent agreement with these new measurements. We also employ a galaxy group finder and show that our model is able to predict: (1) the relative SFRs of central and satellite galaxies, (2) the SFR-dependence of the radial distribution of satellite galaxy populations within galaxy groups, rich groups, and clusters and their surrounding larger scale environments, and (3) the interesting feature that the satellite quenched fraction as a function of projected radial distance from the central galaxy exhibits an ~ r^-.15 slope, independent of environment. The accurate prediction for the spatial distribution of satellites is intriguing given the fact that we do not explicitly model satellite-specific processes after infall, and that in our model the virial radius does not mark a special transition region in the evolution of a satellite, contrary to most galaxy evolution models. The success of the model suggests that present-day galaxy SFR is strongly correlated with halo mass assembly history.

We present a detailed study of how dark matter haloes assemble their mass and grow their (central) potential well. We characterize these via their mass accretion histories (MAHs) and potential well growth histories (PWGHs), which we extract from the Bolshoi simulation and from semi-analytical merger trees supplemented with a method to compute the maximum circular velocity, Vmax, of progenitor haloes. The results of both methods are in excellent agreement, both in terms of the average and the scatter. We show that the MAH and PWGH are tightly correlated, and that growth of the central potential precedes the assembly of mass; the maximum circular velocity is already half the present day value by the time the halo has accreted only 2 percent of its final mass. Finally, we demonstrate that MAHs have a universal form, which we use to develop a new and improved universal model that can be used to compute the average or median MAH and PWGH for a halo of any mass in any LCDM cosmology, without having to run a numerical simulation or a set of halo merger trees.

Satellites within simulated massive clusters are significantly spatially correlated with each other, even when those satellites are not gravitationally bound to each other. This correlation is produced by satellites that entered their hosts relatively recently, and is undetectable for satellites that have resided in their hosts for multiple dynamical timescales. Therefore, a measurement of clustering statistics of cluster satellites may be used to determine the typical accretion redshifts of those satellites into their observed hosts. We argue that such measurements may be used to determine the fraction of satellite galaxies that were quenched by their current hosts, thereby discriminating among models for quenching of star formation in satellite galaxies.

Systematic uncertainties that have been subdominant in past large-scale structure (LSS) surveys are likely to exceed statistical uncertainties of current and future LSS data sets, potentially limiting the extraction of cosmological information. Here we present a general framework (PCA marginalization) to consistently incorporate systematic effects into a likelihood analysis. This technique naturally accounts for degeneracies between nuisance parameters and can substantially reduce the dimension of the parameter space that needs to be sampled. As a practical application, we apply PCA marginalization to account for baryonic physics as an uncertainty in cosmic shear tomography. Specifically, we use CosmoLike to run simulated likelihood analyses on three independent sets of numerical simulations, each covering a wide range of baryonic scenarios differing in cooling, star formation, and feedback mechanisms. We simulate a Stage III (Dark Energy Survey) and Stage IV (Large Synoptic Survey Telescope/Euclid) survey and find a substantial bias in cosmological constraints if baryonic physics is not accounted for. We then show that PCA marginalization (employing at most 3 to 4 nuisance parameters) removes this bias. Our study demonstrates that it is possible to obtain robust, precise constraints on the dark energy equation of state even in the presence of large levels of systematic uncertainty in astrophysical processes. We conclude that the PCA marginalization technique is a powerful, general tool for addressing many of the challenges facing the precision cosmology program.

Quenched central galaxies tend to reside in a preferentially quenched large-scale environment, a phenomenon that has been dubbed galactic conformity. Remarkably, this tendency persists out to scales far larger than the virial radius of the halo hosting the central. Therefore, conformity manifestly violates the widely adopted assumption that the dark matter halo mass Mvir exclusively governs galaxy occupation statistics. This paper is the first in a series studying the implications of the observed conformity signal for the galaxy-dark matter connection. We show that recent measurements of conformity on scales R ~1-5 Mpc imply that central galaxy quenching statistics cannot be correctly predicted with the knowledge of Mvir alone. We also demonstrate that ejected (or `backsplash') satellites cannot give rise to the signal. We then invoke the age matching model, which is predicated on the co-evolution of galaxies and halos. We find that this model produces a strong signal, and that central galaxies are solely responsible. We conclude that large-scale `2-halo' conformity represents a smoking gun of central galaxy assembly bias, and indicates that contemporary models of satellite quenching have systematically over-estimated the influence of post-infall processes.

The age matching model has recently been shown to predict correctly the luminosity L and g-r color of galaxies residing within dark matter halos. The central tenet of the model is intuitive: older halos tend to host galaxies with older stellar populations. In this paper, we demonstrate that age matching also correctly predicts the g-r color trends exhibited in a wide variety of statistics of the galaxy distribution for stellar mass M* threshold samples. In particular, we present new measurements of the galaxy two-point correlation function and the galaxy-galaxy lensing signal as a function of M* and g-r color from the Sloan Digital Sky Survey, and show that age matching exhibits remarkable agreement with these and other statistics of low-redshift galaxies. In so doing, we also demonstrate good agreement between the galaxy-galaxy lensing observed by SDSS and the signal predicted by abundance matching, a new success of this model. We describe how age matching is a specific example of a larger class of Conditional Abundance Matching models (CAM), a theoretical framework we introduce here for the first time. CAM provides a general formalism to study correlations at fixed mass between any galaxy property and any halo property. The striking success of our simple implementation of CAM provides compelling evidence that this technique has the potential to describe the same set of data as alternative models, but with a dramatic reduction in the required number of parameters. CAM achieves this reduction by exploiting the capability of contemporary N-body simulations to determine dark matter halo properties other than mass alone, which distinguishes our model from conventional approaches to the galaxy-halo connection.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

We present age distribution matching, a theoretical formalism for predicting how galaxies of luminosity L and color C occupy dark matter halos. Our model supposes that there are just two fundamental properties of a halo that determine the color and brightness of the galaxy it hosts: the maximum circular velocity Vmax, and the redshift z_starve that correlates with the epoch at which the star formation in the galaxy ceases. The halo property z_starve is intended to encompass physical characteristics of halo mass assembly that may deprive the galaxy of its cold gas supply and, ultimately, quench its star formation. The new, defining feature of the model is that, at fixed luminosity, galaxy color is in monotonic correspondence with z_starve, with the larger values of z_starve being assigned redder colors. We populate an N- body simulation with a mock galaxy catalog based on age distribution matching, and show that the resulting mock galaxy distribution accurately describes a variety of galaxy statistics. Our model suggests that halo and galaxy assembly are indeed correlated. We make publicly available our low-redshift, SDSS M_r <-19 mock galaxy catalog, and main progenitor histories of all z=0 halos, at http://logrus.uchicago.edu/~aphearin

We construct mock catalogs of galaxy groups using subhalo abundance matching (SHAM) and undertake several new tests of the SHAM prescription for the galaxy-dark matter connection. All SHAM models we studied exhibit significant tension with galaxy groups observed in the Sloan Digital Sky Survey (SDSS). The SHAM prediction for the field galaxy luminosity function (LF) is systematically too dim, and the group galaxy LF systematically too bright, regardless of the details of the SHAM prescription. SHAM models connecting r-band luminosity, Mr, to Vacc, the maximum circular velocity of a subhalo at the time of accretion onto the host, faithfully reproduce galaxy group abundance as a function of richness, g(N). However, SHAM models connecting Mr with Vpeak, the peak value of Vmax over the entire merger history of the halo, over-predict galaxy group abundance. Our results suggest that no SHAM model can simultaneously reproduce the observed g(N) and two-point projected galaxy clustering. Nevertheless, we also report a new success of SHAM: an accurate prediction for Phi(m12), the abundance of galaxy groups as a function of magnitude gap m12, defined as the difference between the r-band absolute magnitude of the two brightest group members. We show that it may be possible to use joint measurements of g(N) and Phi(m12) to tightly constrain the details of the SHAM implementation. Additionally, we show that the hypothesis that the luminosity gap is constructed via random draws from a universal LF provides a poor description of the data, contradicting recent claims in the literature. Finally, we test a common assumption of the Conditional Luminosity Function (CLF) formalism, that the satellite LF need only be conditioned by the brightness of the central galaxy. We find this assumption to be well-supported by the observed Phi(m12).

We investigate the potential to improve optical tracers of cluster mass by exploiting measurements of the magnitude gap, m12, defined as the difference between the r-band absolute magnitude of the two brightest cluster members. We find that in a mock sample of galaxy groups and clusters constructed from the Bolshoi simulation, the scatter about the mass-richness relation decreases by 15-20% when magnitude gap information is included. A similar trend is evident in a volume-limited, spectroscopic sample of galaxy groups observed in the Sloan Digital Sky Survey (SDSS). We find that SDSS groups with small magnitude gaps are richer than large-gap groups at fixed values of the one-dimensional velocity dispersion among group members sigma_v, which we use as a mass proxy. We demonstrate explicitly that m12 contains information about cluster mass that supplements the information provided by group richness and the luminosity of the brightest cluster galaxy, L_bcg. In so doing, we show that the luminosities of the members of a group with richness N are inconsistent with the distribution of luminosities that results from N random draws from the global galaxy luminosity function. As the cosmological constraining power of galaxy clusters is limited by the precision in cluster mass determination, our findings suggest a new way to improve the cosmological constraints derived from galaxy clusters.

One of the most pernicious theoretical systematics facing upcoming gravitational lensing surveys is the uncertainty introduced by the effects of baryons on the power spectrum of the convergence field. One method that has been proposed to account for these effects is to allow several additional parameters (that characterize dark matter halos) to vary and to fit lensing data to these halo parameters concurrently with the standard set of cosmological parameters. We test this method. In particular, we use this technique to model convergence power spectrum predictions from a set of cosmological simulations. We estimate biases in dark energy equation of state parameters that would be incurred if one were to fit the spectra predicted by the simulations either with no model for baryons, or with the proposed method. We show that neglecting baryonic effect leads to biases in dark energy parameters that are several times the statistical errors for a survey like the Dark Energy Survey. The proposed method to correct for baryonic effects renders the residual biases in dark energy equation of state parameters smaller than the statistical errors. These results suggest that this mitigation method may be applied to analyze convergence spectra from a survey like the Dark Energy Survey. For significantly larger surveys, such as will be carried out by the Large Synoptic Survey Telescope, the biases introduced by baryonic effects are much more significant. We show that this mitigation technique significantly reduces the biases for such larger surveys, but that a more effective mitigation strategy will need to be developed in order ensure that the residual biases in these surveys fall below the statistical errors.

Forthcoming projects such as DES, LSST, WFIRST, and Euclid aim to measure weak lensing shear correlations with unprecedented precision, constraining the dark energy equation of state at the percent level. Reliance on photometrically-determined redshifts constitutes a major source of uncertainty for these surveys. Additionally, interpreting the weak lensing signal requires a detailed understanding of the nonlinear physics of gravitational collapse. We present a new analysis of the stringent calibration requirements for weak lensing analyses of future imaging surveys that addresses both photo-z uncertainty and errors in the calibration of the matter power spectrum. We find that when photo-z uncertainty is taken into account the requirements on the level of precision in the prediction for the matter power spectrum are more stringent than previously thought. Including degree-scale galaxy clustering statistics in a joint analysis with weak lensing not only strengthens the survey's constraining power by ~20%, but can also have a profound impact on the calibration demands, decreasing the degradation in dark energy constraints with matter power spectrum uncertainty by a factor of 2-5. Similarly, using galaxy clustering information significantly relaxes the demands on photo-z calibration. We compare these calibration requirements to the contemporary state-of-the-art in photometric redshift estimation and predictions of the power spectrum and suggest strategies to utilize forthcoming data optimally.

The cosmic microwave background (CMB) temperature distribution measured by the Wilkinson Microwave Anisotropy Probe (WMAP) exhibits anomalously low correlation at large angles. Quantifying the degree to which this feature in the temperature data is in conflict with standard Lambda-CDM cosmology is somewhat ambiguous because of the a posteriori nature of the observation. One physical mechanism that has been proposed as a possible explanation for the deficit in the large-angle temperature correlations is a suppression of primordial power on ~Gpc scales. To distinguish whether the anomaly is a signal of new physics, such as suppressed primordial power, it would be invaluable to perform experimental tests of the authenticity of this signal in data sets which are independent of the WMAP temperature measurements or even other CMB measurements. We explore the possibility of testing models of power suppression with large-scale structure observations, and compare the ability of planned photometric and spectroscopic surveys to constrain the power spectrum. Of the surveys planned for the next decade, a spectroscopic redshift survey such as BigBOSS will have a greater number of radial modes available for study, but we find that this advantage is outweighed by the greater surface density of high-redshift sources that will be observed by photometric surveys such as LSST or Euclid. We also find that the ability to constrain primordial power suppression is insensitive to the precision of the calibration of photometric redshifts. We conclude that very-wide-area imaging surveys have the potential to probe viable models for the missing power but that it will be difficult to use such surveys to conclusively rule out primordial power suppression as the mechanism behind the observed anomaly.

We study energy transport by asymmetric dark matter in the interiors of very low-mass stars and brown dwarfs. Our motivation is to explore astrophysical signatures of asymmetric dark matter, which otherwise may not be amenable to conventional indirect dark matter searches. In viable models, the additional cooling of very-low mass stellar cores can alter stellar properties. Asymmetric dark matter with mass 4 < Mx/GeV < 10 and a spin-dependent (spin-independent) cross sections of sigma \sim 10^{-37} cm^2 (sigma \sim 10^{-40} cm^2) can increase the minimum mass of main sequence hydrogen burning, partly determining whether or not the object is a star at all. Similar dark matter candidates reduce the luminosities of low-mass stars and accelerate the cooling of brown dwarfs. Such light dark matter is of particular interest given results from the DAMA, CoGeNT, and CRESST dark matter searches. We discuss possibilities for observing dark matter effects in stars in the solar neighborhood, globular clusters, and, of particular promise, local dwarf galaxies, among other environments, as well as exploiting these effects to constrain dark matter properties.

A goal of forthcoming imaging surveys is to use weak gravitational lensing shear measurements to constrain dark energy. We quantify the importance of uncalibrated photometric redshift outliers to the dark energy goals of forthcoming imaging surveys in a manner that does not assume any particular photometric redshift technique or template. In so doing, we provide an approximate blueprint for computing the influence of specific outlier populations on dark energy constraints. We find that outliers whose photo-z distributions are tightly localized about a significantly biased redshift must be controlled to a per-galaxy rate of <~ a few times 10^-3 to insure that systematic errors on dark energy parameters are rendered negligible. In the complementary limit, a subset of imaged galaxies with uncalibrated photometric redshifts distributed over a broad range must be limited to fewer than a per-galaxy error rate of <~ a few times 10^-4. Additionally, we explore the relative importance of calibrating the photo-z's of a core set of relatively well-understood galaxies as compared to the need to identify potential catastrophic photo-z outliers. We discuss the degradation of the statistical constraints on dark energy parameters induced by excising source galaxies at high- and low-photometric redshifts, concluding that removing galaxies with z_phot >~ 2.4 and z_phot <~ 0.3 may mitigate damaging catastrophic redshift outliers at a relatively small (~ 20%) cost in statistical error. In an appendix, we show that forecasts for the degradation in dark energy parameter constraints due to uncertain photometric redshifts depend sensitively on the treatment of the nonlinear matter power spectrum. Previous work using PD96 may have overestimated the photo-z calibration requirements of future surveys.

Forthcoming projects such as the DES, a JDEM, and LSST, aim to measure weak lensing shear correlations with unprecedented accuracy. Weak lensing observables are sensitive to both the distance-redshift relation and the growth of structure in the Universe. If the cause of accelerated cosmic expansion is dark energy within general relativity (GR), both cosmic distances and structure growth are governed by the properties of dark energy. Consequently, one may use lensing to check for this consistency and test GR. After reviewing the phenomenology of such tests, we address one major challenge to such a program. The evolution of the baryonic component of the Universe is highly uncertain and can influence lensing observables, manifesting as modified structure growth for a fixed cosmic distance scale. Using two proposed methods, we show that one could be led to reject the null hypothesis of GR when it is the true theory if this uncertainty in baryonic processes is neglected. Recent simulations suggest that we can correct for baryonic effects using a parametrized model in which the halo mass-concentration relation is modified. The correction renders biases small compared to statistical uncertainties. We study the ability of future weak lensing surveys to constrain the internal structures of halos and test the null hypothesis of GR simultaneously. Compared to nulling information from small-scales to mitigate sensitivity to baryonic physics, this internal calibration program should provide limits on deviations from GR that are several times more constraining. Specifically, we find that limits on general relativity in the case of internal calibration are degraded by only ~30% or less compared to the case of perfect knowledge of nonlinear structure.