Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all halo mass, including mass in subhaloes, in studies of halo density profiles. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host. This difference makes comparisons between theoretical and observed results difficult. In this work we investigate density profiles of the smooth components of host haloes by excluding mass contained within subhaloes. We find that the density profiles of the smooth halo component (without subhaloes) differs substantially from the conventional halo density profile. Smooth profiles decline more rapidly at large radii and are not well characterised by the standard NFW profile. We also find that concentrations derived from smooth density profiles exhibit less scatter at fixed mass and a weaker mass dependence than standard concentrations. Both smooth and standard halo profiles can be described by a generalised Einasto profile, an Einasto profile with a modified central slope, with smaller residuals than either an NFW or Einasto profile. These results hold for both Milky Way-mass and cluster-mass haloes. This new characterisation of smooth halo profiles can be useful for many analyses, such as lensing and dark matter annihilation, in which the smooth and clumpy components of a halo should be accounted for separately.

We investigate analytically whether in a close encounter with a star, a Dark Matter particle can be accelerated above the escape speed of a Globular Cluster and be ejected. We find that this mechanism is not sufficient to eject a massive, extended Dark Matter halo by the present time. Combined with observations of isolated Globular Clusters that may not have had their halos tidally stripped, these results cast doubt on the scenario in which Globular Clusters formed in Dark Matter halos.

The concentration parameter is a key characteristic of a dark matter halo that conveniently connects the halo's present-day structure with its assembly history. Using 'Dark Sky', a suite of cosmological $N$-body simulations, we investigate how halo concentration evolves with time and emerges from the mass assembly history. We also explore the origin of the scatter in the relation between concentration and assembly history. We show that the evolution of halo concentration has two primary modes: (1) smooth increase due to pseudo-evolution; and (2) intense responses to physical merger events. Merger events induce lasting and substantial changes in halo structures, and we observe a universal response in the concentration parameter. We argue that merger events are a major contributor to the uncertainty in halo concentration at fixed halo mass and formation time. In fact, even haloes that are typically classified as having quiescent formation histories experience multiple minor mergers. These minor mergers drive small deviations from pseudo-evolution, which cause fluctuations in the concentration parameters and result in effectively irreducible scatter in the relation between concentration and assembly history. Hence, caution should be taken when using present-day halo concentration parameter as a proxy for the halo assembly history, especially if the recent merger history is unknown.

We study the non-adiabatic cosmological production of ultra light dark matter (ULDM) under a minimal set of assumptions: a free ultra light real scalar as a spectator field in its Bunch-Davies vacuum state during inflation and instantaneous reheating into a radiation dominated era. For (ULDM) fields minimally coupled to gravity, non-adiabatic particle production yields a \emph{distribution function} peaked at \emph{low} comoving momentum $\mathcal{N}_k \propto 1/k^3$. The infrared behavior is a remnant of the infrared enhancement of light minimally coupled fields during inflation. We obtain the full energy momentum tensor, show explicity its equivalence with the fluid-kinetic one in the adiabatic regime, and extract the abundance, equation of state and free streaming length (cutoff in the matter power spectrum). Taking the upper bound on the scale of inflation from Planck, the (UDLM) saturates the dark matter abundance for $m \simeq 1.5\,\times 10^{-5}\mathrm{eV}$ with an equation of state parameter $w \simeq 10^{-14}$ and a free streaming length $\lambda_{fs} \simeq 70\,\mathrm{pc}$. Thus this cosmologically produced (ULDM) yields a \emph{cold} dark matter particle. We argue that the abundance from non-adiabatic production yields a \emph{lower bound} on generic (ULDM) and axion-like particles that must be included in any assessment of (ULDM) as a dark matter candidate.

We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging dataset, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The LRG sample contains 2.7 million objects over the redshift range of $0.4 < z < 0.9$ over 5655 sq. degrees. We have developed new photometric redshift (photo-$z$) estimates using the Legacy Survey DECam and WISE photometry, with $\sigma_{\mathrm{NMAD}} = 0.02$ precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a 5-parameter HOD model fit using the MultiDark halo catalog, we find that there is little evolution in HOD parameters except at the highest-redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of MCMC fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogs, and present a new technique for finding best-fit parameters in this situation. Accompanying this paper we have released the PRLS (Photometric Redshifts for the Legacy Surveys) catalog of photo-$z$'s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 dataset. This catalog provides accurate photometric redshifts for objects with $z < 21$ over more than 16,000 square degrees of sky.

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.

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.

We present a Bayesian hierarchical inference formalism (Basilisk) to constrain the galaxy-halo connection using satellite kinematics. Unlike traditional methods, Basilisk does not resort to stacking the kinematics of satellite galaxies in bins of central luminosity, and does not make use of summary statistics, such as satellite velocity dispersion. Rather, Basilisk leaves the data in its raw form and computes the corresponding likelihood. In addition, Basilisk can be applied to flux-limited, rather than volume-limited samples, greatly enhancing the quantity and dynamic range of the data. And finally, Basilisk is the only available method that simultaneously solves for halo mass and orbital anisotropy of the satellite galaxies, while properly accounting for scatter in the galaxy-halo connection. Basilisk uses the conditional luminosity function to model halo occupation statistics, and assumes that satellite galaxies are a relaxed tracer population of the host halo's potential with kinematics that obey the spherical Jeans equation. We test and validate Basilisk using mocks of varying complexity, and demonstrate that it yields unbiased constraints on the galaxy-halo connection and at a precision that rivals galaxy-galaxy lensing. In particular, Basilisk accurately recovers the full PDF of the relation between halo mass and central galaxy luminosity, and simultaneously constrains the orbital anisotropy of the satellite galaxies. Basilisk's inference is not affected by potential velocity bias of the central galaxies, or by slight errors in the inferred, radial profile of satellite galaxies that arise as a consequence of interlopers and sample impurity.

On small scales there have been a number of claims of discrepancies between the standard Cold Dark Matter (CDM) model and observations. The 'missing satellites problem' infamously describes the overabundance of subhalos from CDM simulations compared to the number of satellites observed in the Milky Way. A variety of solutions to this discrepancy have been proposed; however, the impact of the specific properties of the Milky Way halo relative to the typical halo of its mass have yet to be explored. Motivated by recent studies that identified ways in which the Milky Way is atypical (e.g., Licquia et al. 2015), we investigate how the properties of dark matter halos with mass comparable to our Galaxy's --- including concentration, spin, shape, and scale factor of the last major merger --- correlate with the subhalo abundance. Using zoom-in simulations of Milky Way-like halos, we build two models of subhalo abundance as functions of host halo properties and conclude that the Milky Way should be expected to have 22%-44% fewer subhalos with low maximum rotation velocities ($V_{\rm max}^{\rm sat} \sim 10$kms$^{-1}$) at the 95% confidence level and up to 72% fewer than average subhalos with high rotation velocities ($V_{\rm max}^{\rm sat} \gtrsim 30$kms$^{-1}$, comparable to the Magellanic Clouds) than would be expected for a typical halo of the Milky Way's mass. Concentration is the most informative single parameter for predicting subhalo abundance. Our results imply that models tuned to explain the missing satellites problem assuming typical subhalo abundances for our Galaxy will be over-correcting.

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/ .

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.

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.

We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter haloes, as well as the radial profile of satellite galaxies. Confirming previous results, we find that red central galaxies live in more massive haloes than blue galaxies at fixed luminosity. Additionally, our results suggest that satellite galaxies have a radial profile less centrally concentrated than dark matter but not as cored as resolved subhaloes in dark matter-only simulations. Compared to previous works using satellite kinematics by More et al., we find much more competitive constraints on the galaxy-halo connection, on par with those derived from a combination of galaxy clustering and galaxy-galaxy lensing. We compare our results on the galaxy-halo connection to other studies using galaxy clustering and group catalogues, showing very good agreement between these different techniques. We discuss future applications of satellite kinematics in the context of constraining cosmology and the relationship between galaxies and dark matter haloes.

The kinematics of satellite galaxies moving in a dark matter halo are a direct probe of the underlying gravitational potential. Thus, the phase-space distributions of satellites represent a powerful tool to determine the galaxy-halo connection from observations. By stacking the signal of a large number of satellite galaxies this potential can be unlocked even for haloes hosting a few satellites on average. In this work, we test the impact of various modelling assumptions on constraints derived from analysing satellite phase-space distributions in the non-linear, 1-halo regime. We discuss their potential to explain the discrepancy between average halo masses derived from satellite kinematics and gravitational lensing previously reported. Furthermore, we develop an updated, more robust analysis to extract constraints on the galaxy-halo relation from satellite properties in spectroscopic galaxy surveys such as the SDSS. We test the accuracy of this approach using a large number of realistic mock catalogues. Furthermore, we find that constraints derived from such an analysis are complementary and competitive with respect to the commonly used galaxy clustering and galaxy-galaxy lensing observables.

We study scalar particle decay during the radiation and matter dominated epochs of a standard cosmological model. An adiabatic approximation is introduced that is valid for degrees of freedom with typical wavelengths much smaller than the particle horizon ($\propto$~Hubble radius) at a given time. We implement a non-perturbative method that includes the cosmological expansion and obtain a cosmological Fermi's Golden Rule that enables one to compute the decay law of a parent particle of mass $m_1$, along with the build up of the population of daughter particles of mass $m_2$. The survival probability of the decaying particle is $P(t)=e^{-\widetilde{\Gamma}_k(t)\,t}$ with $\widetilde{\Gamma}_k(t)$ being an \emph{effective momentum and time dependent decay rate}. It features a transition time scale $t_{nr}$ between the relativistic and non-relativistic regimes and for $k \neq 0$ is always smaller than the analogous rate in Minkowski spacetime, as a consequence of (local) time dilation and the cosmological redshift. For $t \ll t_{nr}$ the decay law is a "stretched exponential" $P(t) = e^{-(t/t^*)^{3/2}}$, whereas for the non-relativistic stage with $t \gg t_{nr}$, we find $P(t) = e^{-\Gamma_0 t}\,(t/t_{nr})^{\Gamma_0\,t_{nr}/2}$. The Hubble time scale $\propto 1/H(t)$ introduces an energy uncertainty $\Delta E \sim H(t)$ which relaxes the constraints of kinematic thresholds. This opens new decay channels into heavier particles for $2\pi E_k(t) H(t) \gg 4m^2_2-m^2_1$, with $E_k(t)$ the (local) comoving energy of the decaying particle. As the expansion proceeds this channel closes and the usual two particle thresholds restrict the decay kinematics.

Secondary halo bias, commonly known as 'assembly bias,' is the dependence of halo clustering on a halo property other than mass. This prediction of the Lambda-Cold Dark Matter cosmology is essential to modelling the galaxy distribution to high precision and interpreting clustering measurements. As the name suggests, different manifestations of secondary halo bias have been thought to originate from halo assembly histories. We show conclusively that this is incorrect for cluster-size haloes. We present an up-to-date summary of secondary halo biases of high-mass haloes due to various halo properties including concentration, spin, several proxies of assembly history, and subhalo properties. While concentration, spin, and the abundance and radial distribution of subhaloes exhibit significant secondary biases, properties that directly quantify halo assembly history do not. In fact, the entire assembly histories of haloes in pairs are nearly identical to those of isolated haloes. In general, a global correlation between two halo properties does not predict whether or not these two properties exhibit similar secondary biases. For example, assembly history and concentration (or subhalo abundance) are correlated for both paired and isolated haloes, but follow slightly different conditional distributions in these two cases. This results in a secondary halo bias due to concentration (or subhalo abundance), despite the lack of assembly bias in the strict sense for cluster-size haloes. Due to this complexity, caution must be exercised in using any one halo property as a proxy to study the secondary bias due to another property.

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.

Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $\Delta$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $\Delta$. For conventional halo definitions ($\Delta \sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $\Delta \sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $\Delta$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $\Delta \gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass.

We revisit SN1987A constraints on light, hidden sector gauge bosons ("dark photons") that are coupled to the standard model through kinetic mixing with the photon. These constraints are realized because excessive bremsstrahlung radiation of the dark photon can lead to rapid cooling of the SN1987A progenitor core, in contradiction to the observed neutrinos from that event. The models we consider are of interest as phenomenological models of strongly self-interacting dark matter. We clarify several possible ambiguities in the literature and identify errors in prior analyses. We find constraints on the dark photon mixing parameter that are in rough agreement with the early estimates of Dent et al., but only because significant errors in their analyses fortuitously canceled. Our constraints are in good agreement with subsequent analyses by Rrapaj & Reddy and Hardy & Lasenby. We estimate the dark photon bremsstrahlung rate using one-pion exchange (OPE), while Rrapaj & Reddy use a soft radiation approximation (SRA) to exploit measured nuclear scattering cross sections. We find that the differences between mixing parameter constraints obtained through the OPE approximation or the SRA approximation are roughly a factor of $\sim 2-3$. Hardy & Laseby include plasma effects in their calculations finding significantly weaker constraints on dark photon mixing for dark photon masses below $\sim 10\, \mathrm{MeV}$. We do not consider plasma effects. Lastly, we point out that the properties of the SN1987A progenitor core remain somewhat uncertain and that this uncertainty alone causes uncertainty of at least a factor of $\sim 2-3$ in the excluded values of the dark photon mixing parameter. Further refinement of these estimates is unwarranted until either the interior of the SN1987A progenitor is more well understood or additional, large, and heretofore neglected effects, are identified.

We compute the weak lensing convergence power spectrum, $C^{\kappa\kappa}(\theta)$, in a dust-filled universe using fully non-linear general relativistic simulations. The spectrum is then compared to more standard, approximate calculations by computing the Bardeen (Newtonian) potentials in linearized gravity and utilizing the Born approximation. We find corrections to the angular power spectrum amplitude of order ten percent at very large angular scales, $\ell ~ 2-3$, and percent-level corrections at intermediate angular scales of $\ell ~ 20-30$.

Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxies produced in these models compared to those based on $V_{\rm peak}$. Despite success in predicting clustering, a simple $V_{\rm peak}$ based SHAM model results in predictions for galaxy growth that are at odds with observations. We evaluate three possible remedies that could "save" mass-based SHAM: (1) SHAM models require a significant population of "orphan" galaxies as a result of artificial disruption/merging of sub-haloes in modern high resolution dark matter simulations; (2) satellites must grow significantly after their accretion; and (3) stellar mass is significantly affected by halo assembly history. No solution is entirely satisfactory. However, regardless of the particulars, we show that popular SHAM models based on $M_{\rm peak}$ cannot be complete physical models as presented. Either $V_{\rm peak}$ truly is a better predictor of stellar mass at $z\sim 0$ and it remains to be seen how the correlation between stellar mass and $V_{\rm peak}$ comes about, or SHAM models are missing vital component(s) that significantly affect galaxy clustering.

A strong correlation has been measured between the observed centripetal accelerations in galaxies and the accelerations implied by the baryonic components of galaxies. This empirical radial acceleration relation must be accounted for in any viable model of galaxy formation. We measure and compare the radial accelerations contributed by baryons and by dark matter in disk galaxies in the MassiveBlack-II hydrodynamic galaxy formation simulation. The sample of 1594 galaxies spans three orders of magnitude in luminosity and four in surface brightness, comparable to the observed sample from the Spitzer Photometry & Accurate Rotation Curves (SPARC) dataset used by McGaugh et al. (2016). We find that radial accelerations contributed by baryonic matter only and by total matter are highly correlated, with only small scatter around their mean or median relation, despite the wide ranges of galaxy luminosity and surface brightness. We further find that the radial acceleration relation in this simulation differs from that of the SPARC sample, and can be described by a simple power law in the acceleration range we are probing.

Traditional cosmological inference using Type Ia supernovae (SNeIa) have used stretch- and color-corrected fits of SN Ia light curves and assumed a resulting fiducial mean and symmetric intrinsic dispersion for the resulting relative luminosity. As systematics become the main contributors to the error budget, it has become imperative to expand supernova cosmology analyses to include a more general likelihood to model systematics to remove biases with losses in precision. To illustrate an example likelihood analysis, we use a simple model of two populations with a relative luminosity shift, independent intrinsic dispersions, and linear redshift evolution of the relative fraction of each population. Treating observationally viable two-population mock data using a one-population model results in an inferred dark energy equation of state parameter $w$ that is biased by roughly 2 times its statistical error for a sample of N $ \gtrsim$ 2500 SNeIa. Modeling the two-population data with a two-population model removes this bias at a cost of an approximately $\sim20\%$ increase in the statistical constraint on $w$. These significant biases can be realized even if the support for two underlying SNeIa populations, in the form of model selection criteria, is inconclusive. With the current observationally-estimated difference in the two proposed populations, a sample of N $ \gtrsim$ 10,000 SNeIa is necessary to yield conclusive evidence of two populations.

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.

We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold dark matter (CDM), warm dark matter (WDM), and decaying dark matter (DDM) simulations for our models of structure formation. The subhalo candidates in CDM typically have smaller mass and higher concentrations at z = 0 than the corresponding candidates in WDM and DDM. We examine the formation histories of the ~ 100 Fornax candidate subhalos identified in CDM simulations and, using approximate luminosity-mass relationships for subhalos, we find two of these subhalos that are consistent with both the Fornax luminosity and kinematics. These two subhalos have a peak mass over ten times larger than their z = 0 mass. We suggest that in CDM the dark matter halo hosting Fornax must have been severely stripped of mass and that it had an infall time into the Milky Way of ~ 9 Gyr ago. In WDM, we find that candidate subhalos consistent with the properties of Fornax have a similar infall time and a similar degree of mass loss, while in DDM we find a later infall time of ~ 3 - 4 Gyr ago and significantly less mass loss. We discuss these results in the context of the Fornax star formation history, and show that these predicted subhalo infall times can be linked to different star formation quenching mechanisms. This emphasizes the links between the properties of the dark matter and the mechanisms that drive galaxy evolution.

White Dwarfs (WD) capture Dark Matter (DM) as they orbit within their host halos. These captured particles may subsequently annihilate, heating the stellar core and preventing the WD from cooling. The potential wells of WDs are considerably deeper and core temperatures significantly cooler than those of main sequence stars. Consequently, DM evaporation is less important in WDs and DM with masses $M_{\chi} \gtrsim 100\, \kev$ and annihilation cross-sections orders of magnitude below the canonical thermal cross-section ($\sigmav \gtrsim 10^{-46}\, \cm^3$/s) can significantly alter WD cooling in particular astrophysical environments. We consider WDs in globular clusters (GCs) and dwarf galaxies. If the parameters of the DM particle are known, then the temperature of the coolest WD in a GC can be used to constrain the DM density of the cluster's halo (potentially even ruling out the presence of a halo if the inferred density is of order the ambient Galactic density). Recently several direct detection experiments have seen signals whose origins might be due to low mass DM. In this paper, we show that if these claims from CRESST, DAMA, CDMS-Si, and CoGeNT could be interpreted as DM, then observations of NGC 6397 limit the fraction of DM in that cluster to be $f_{\mathrm{DM}} \lesssim 10^{-3}$. This would be an improvement over existing constraints of 3 orders of magnitude and clearly rule out a significant DM halo for this cluster...

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.

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 examine the importance of baryonic feedback effects on the matter power spectrum on small scales, and the implications for the precise measurement of neutrino masses through gravitational weak lensing. Planned large galaxy surveys such as the Large Synoptic Sky Telescope (LSST) and Euclid are expected to measure the sum of neutrino masses to extremely high precision, sufficient to detect non-zero neutrino masses even in the minimal mass normal hierarchy. We show that weak lensing of galaxies while being a very good probe of neutrino masses, is extremely sensitive to baryonic feedback processes. We use publicly available results from the Overwhelmingly Large Simulations (OWLS) project to investigate the effects of active galactic nuclei feedback, the nature of the stellar initial mass function, and gas cooling rates, on the measured weak lensing shear power spectrum. Using the Fisher matrix formalism and priors from CMB+BAO data, we show that when one does not account for feedback, the measured neutrino mass may be substantially larger or smaller than the true mass, depending on the dominant feedback mechanism, with the mass error |\Delta m_nu| often exceeding the mass m_nu itself. We also consider gravitational lensing of the cosmic microwave background (CMB) and show that it is not sensitive to baryonic feedback on scales l < 2000, although CMB experiments that aim for sensitivities sigma(m_nu) < 0.02 eV will need to include baryonic effects in modeling the CMB lensing potential. A combination of CMB lensing and galaxy lensing can help break the degeneracy between neutrino masses and baryonic feedback processes. We conclude that future large galaxy lensing surveys such as LSST and Euclid can only measure neutrino masses accurately if the matter power spectrum can be measured to similar accuracy.

We present a set of N-body simulations of a class of models in which an unstable dark matter particle decays into a stable non-interacting dark matter particle, with decay lifetime comparable to the Hubble time. We study the effects of the kinematic recoil velocity received by the stable dark matter on the structures of dark matter halos ranging from galaxy-cluster to Milky Way mass scales. For Milky Way-mass halos, we use high-resolution, zoom-in simulations to explore the effects of decays on Galactic substructure. In general, halos with circular velocities comparable to the magnitude of kick velocity are most strongly affected by decays. We show that decaying dark matter models with lifetimes comparable to Hubble time and recoil speeds about 20-40 km/s can significantly reduce both the abundance of Galactic subhalos and the internal densities of the subhalos. We also compare subhalo circular velocity profiles with observational constraints on the Milky Way dwarf satellite galaxies. Interestingly, we find that decaying dark matter models that do not violate current astrophysical constraints, can significantly mitigate both the well-documented "missing satellites problem" and the more recent "too big to fail problem" associated with the abundances and densities of Local Group dwarf satellite galaxies. A relatively unique feature of late decaying dark matter models is that they predict significant evolution of halos as a function of time. This is an important consideration because at high redshifts, prior to decays, decaying models exhibit the same sequence of structure formation as cold dark matter. We conclude that models of decaying dark matter make predictions that are relevant for the interpretation of observations of small galaxies in the Local Group and can be tested or constrained by the kinematics of Local Group dwarf galaxies as well as by forthcoming large-scale surveys.

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.

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.

We have devised a method to select galaxies that are isolated in their dark matter halo (N=1 systems) and galaxies that reside in a group of exactly two (N=2 systems). Our N=2 systems are widely-separated (up to $\sim$\,200\,$h^{-1}$\,kpc), where close galaxy-galaxy interactions are not dominant. We apply our selection criteria to two volume-limited samples of galaxies from SDSS DR6 with $M_{r}-5 \log_{10} h \leq$ -19 and -20 to study the effects of the environment of very sparse groups on galaxy colour. For satellite galaxies in a group of two, we find a red excess attributed to star formation quenching of 0.15\,$\pm$\,0.01 and 0.14\,$\pm$\,0.01 for the -19 and -20 samples, respectively, relative to isolated galaxies of the same stellar mass. Assuming N=1 systems are the progenitors of N=2 systems, an immediate-rapid star formation quenching scenario is inconsistent with these observations. A delayed-then-rapid star formation quenching scenario with a delay time of 3.3 and 3.7\,Gyr for the -19 and -20 samples, respectively, yields a red excess prediction in agreement with the observations. The observations also reveal that central galaxies in a group of two have a slight blue excess of 0.06\,$\pm$\,0.02 and 0.02\,$\pm$\,0.01 for the -19 and -20 samples, respectively, relative to N=1 populations of the same stellar mass. Our results demonstrate that even the environment of very sparse groups of luminous galaxies influence galaxy evolution and in-depth studies of these simple systems are an essential step towards understanding galaxy evolution in general.

We present an analysis of high-resolution N-body simulations of decaying dark matter cosmologies focusing on the statistical properties of the transmitted Lyman-alpha forest flux in the high-redshift intergalactic medium. In this type of model a dark matter particle decays into a slightly less massive stable dark matter daughter particle and a comparably light particle. The small mass splitting provides a non-relativistic kick velocity V_k to the daughter particle resulting in free-streaming and subsequent damping of small-scale density fluctuations. Current Lyman-alpha forest power spectrum measurements probe comoving scales up to ~ 2-3 h^-1 Mpc at redshifts z ~ 2-4, providing one of the most robust ways to probe cosmological density fluctuations on relatively small scales. The suppression of structure growth due to the free-streaming of dark matter daughter particles also has a significant impact on the neutral hydrogen cloud distribution, which traces the underlying dark matter distribution well at high redshift. We exploit Lyman-alpha forest power spectrum measurements to constrain the amount of free-streaming of dark matter in such models and thereby place limits on decaying dark matter based only on the dynamics of cosmological perturbations without any assumptions about the interactions of the decay products. We find that SDSS 1D Lyman-alpha forest power spectrum data place a lifetime-dependent upper limit V_k < 30-70 km/s for decay lifetimes < 10 Gyr. This is the most stringent model-independent bound on invisible dark matter decays with small mass splittings. For large mass splittings (large V_k), Lyman-alpha forest data restrict the dark matter lifetime to Gamma^-1 > 40 Gyr. Forthcoming BOSS data should be able to provide more stringent constraints on exotic dark matter, mainly because the larger BOSS quasar spectrum sample will significantly reduce statistical errors.

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 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 introduce a new technique that uses galaxy clustering to constrain how satellite galaxies lose stellar mass and contribute to the diffuse "intrahalo light" (IHL). We implement two models that relate satellite galaxy stellar mass loss to the detailed knowledge of subhalo dark matter mass loss. Model 1 assumes that the fractional stellar mass loss of a galaxy is proportional to the fractional amount of dark matter mass loss of its subhalo. Model 2 accounts for a delay in the time that stellar mass is lost since the galaxy resides deep in the potential well of the subhalo which may experience dark matter mass loss for some time before the galaxy is affected. We use these models to predict the stellar masses of a population of galaxies and use abundance matching to predict the clustering of several r-band luminosity threshold samples from the Sloan Digital Sky Survey. Abundance matching assuming no stellar mass loss (akin to abundance matching at the time of subhalo infall) over-estimates the correlation function on small scales (<~ 1 Mpc), while allowing too much stellar mass loss leads to an under-estimate. For each sample, we are thus able to constrain the amount of stellar mass loss required to match the observed clustering. We find that less luminous satellite galaxies experience more efficient stellar mass loss than luminous satellites. From these models, we can infer the amount of stellar mass that is deposited into the IHL. We find that both of our model predictions for the mean amount of IHL as a function of halo mass are consistent with current observational measurements. However, our two models predict a different amount of scatter in the IHL from halo to halo, with Model 2 being favored by observations. This demonstrates that a comparison to IHL measurements provides independent verification of our stellar mass loss models. (Abridged)

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.

Recent kinematical constraints on the internal densities of the Milky Way's dwarf satellites have revealed a discrepancy with the subhalo populations of simulated Galaxy-scale halos in the standard CDM model of hierarchical structure formation. This has been dubbed the "too big to fail" problem, with reference to the improbability of large and invisible companions existing in the Galactic environment. In this paper, we argue that both the Milky Way observations and simulated subhalos are consistent with the predictions of the standard model for structure formation. Specifically, we show that there is significant variation in the properties of subhalos among distinct host halos of fixed mass and suggest that this can reasonably account for the deficit of dense satellites in the Milky Way. We exploit well-tested analytic techniques to predict the properties in a large sample of distinct host halos with a variety of masses spanning the range expected of the Galactic halo. The analytic model produces subhalo populations consistent with both Via Lactea II and Aquarius, and our results suggest that natural variation in subhalo properties suffices to explain the discrepancy between Milky Way satellite kinematics and these numerical simulations. At least ~10% of Milky Way-sized halos host subhalo populations for which there is no "too big to fail" problem, even when the host halo mass is as large as M_host = 10^12.2 h^-1 M_sun. Follow-up studies consisting of high-resolution simulations of a large number of Milky Way-sized hosts are necessary to confirm our predictions. In the absence of such efforts, the "too big to fail" problem does not appear to be a significant challenge to the standard model of hierarchical formation. [abridged]

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.

Simulations of cluster formation have demonstrated that condensation of baryons into central galaxies during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder, even at radii as large as half of the virial radius. However, such simulations generally predict stellar fractions within cluster virial radii that are ~2 to 3 times larger than the stellar masses deduced from observations. In this work we compare ellipticity profiles of clusters simulated with and without baryonic cooling to the cluster ellipticity profiles derived from Chandra and ROSAT observations in an effort to constrain the fraction of gas that cools and condenses into the central galaxies within clusters. We find that the observed ellipticity profiles are fairly constant with radius, with an average ellipticity of 0.18 +/- 0.05. The observed ellipticity profiles are in good agreement with the predictions of non-radiative simulations. On the other hand, the ellipticity profiles of the clusters in simulations that include radiative cooling, star formation, and supernova feedback (but no AGN feedback) deviate significantly from the observed ellipticity profiles at all radii. The simulations with cooling overpredict (underpredict) ellipticity in the inner (outer) regions of galaxy clusters. By comparing the simulations with and without cooling, we show that the cooling of gas via cooling flows in the central regions of simulated clusters causes the gas distribution to be more oblate in the central regions, but makes the outer gas distribution more spherical. We find that late-time gas cooling and star formation are responsible for the significantly oblate gas distributions in cluster cores, but the gas shapes outside of cluster cores are set primarily by baryon dissipation at high redshift z > 2.

We analyze self-consistent N-body simulations of the Milky Way disk and the ongoing disruption of the Sagittarius dwarf satellite to study the effect of Sagittarius tidal debris on dark matter detection experiments. In agreement with significant previous work, we reiterate that the standard halo model is insufficient to describe the non-Maxwellian velocity distribution of the Milky Way halo in our equilibrium halo-only and halo/galaxy models, and offer suggestions for correcting for this discrepancy. More importantly, we emphasize that the dark matter component of the leading tidal arm of the Sagittarius dwarf is significantly more extended than the stellar component of the arm, since the dark matter and stellar streams are not necessarily coaxial and may be offset by several kpc at the point at which they impact the Galactic disk. This suggests that the dark matter component of the Sagittarius debris is likely to have a non-negligible influence on dark matter detection experiments even when the stellar debris is centered several kpc from the solar neighborhood. Relative to models without an infalling Sagittarius dwarf, the Sagittarius dark matter debris in our models induces an energy-dependent enhancement of direct search event rates of as much as ~20 - 45%, an energy-dependent reduction in the amplitude of the annual modulation of the event rate by as much as a factor of two, a shift in the phase of the annual modulation by as much as ~20 days, and a shift in the recoil energy at which the modulation reverses phase. These influences of Sagittarius are of general interest in the interpretation of dark matter searches, but may be particularly important in the case of relatively light (m_X < 20 GeV) dark matter because the Sagittarius stream impacts the solar system at high speed compared to the primary halo dark matter.

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.

In this paper we explore the effect of decaying dark matter (DDM) on large-scale structure and possible constraints from galaxy imaging surveys. DDM models have been studied, in part, as a way to address apparent discrepancies between the predictions of standard cold dark matter models and observations of galactic structure. Our study is aimed at developing independent constraints on these models. In such models, DDM decays into a less massive, stable dark matter (SDM) particle and a significantly lighter particle. The small mass splitting between the parent DDM and the daughter SDM provides the SDM with a recoil or "kick" velocity vk, inducing a free-streaming suppression of matter fluctuations. This suppression may be probed via weak lensing power spectra measured by a number of forthcoming imaging surveys that aim primarily to constrain dark energy. Using scales on which linear perturbation theory alone is valid (multipoles < 300), surveys like Euclid or LSST can be sensitive to vk > 90 km/s for lifetimes ~ 1-5 Gyr. To estimate more aggressive constraints, we model nonlinear corrections to lensing power using a simple halo evolution model that is in good agreement with numerical simulations. In our most ambitious forecasts, using multipoles < 3000, we find that imaging surveys can be sensitive to vk ~ 10 km/s for lifetimes < 10 Gyr. Lensing will provide a particularly interesting complement to existing constraints in that they will probe the long lifetime regime far better than contemporary techniques. A caveat to these ambitious forecasts is that the evolution of perturbations on nonlinear scales will need to be well calibrated by numerical simulations before they can be realized. This work motivates the pursuit of such a numerical simulation campaign to constrain dark matter with cosmological weak lensing.

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 perform a pixel-by-pixel analysis of 467 galaxies in the GOODS-VIMOS survey to study systematic effects in extracting properties of stellar populations (age, dust, metallicity and SFR) from pixel colors using the pixel-z method. The systematics studied include the effect of the input stellar population synthesis model, passband limitations and differences between individual SED fits to pixels and global SED-fitting to a galaxy's colors. We find that with optical-only colors, the systematic errors due to differences among the models are well constrained. The largest impact on the age and SFR e-folding time estimates in the pixels arises from differences between the Maraston models and the Bruzual&Charlot models, when optical colors are used. This results in systematic differences larger than the 2{\sigma} uncertainties in over 10 percent of all pixels in the galaxy sample. The effect of restricting the available passbands is more severe. In 26 percent of pixels in the full sample, passband limitations result in systematic biases in the age estimates which are larger than the 2{\sigma} uncertainties. Systematic effects from model differences are reexamined using Near-IR colors for a subsample of 46 galaxies in the GOODS-NICMOS survey. For z > 1, the observed optical/NIR colors span the rest frame UV-optical SED, and the use of different models does not significantly bias the estimates of the stellar population parameters compared to using optical-only colors. We then illustrate how pixel-z can be applied robustly to make detailed studies of substructure in high redshift galaxies such as (a) radial gradients of age, SFR, sSFR and dust and (b) the distribution of these properties within subcomponents such as spiral arms and clumps. Finally, we show preliminary results from the CANDELS survey illustrating how the new HST/WFC3 data can be exploited to probe substructure in z~1-3 galaxies.

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.

We explore the halo-to-halo variation of dark matter substructure in galaxy-sized dark matter halos, focusing on its implications for strongly gravitational lensed systems. We find that the median value for projected substructure mass fractions within projected radii of 3% of the host halo virial radius is approximately fsub ~ 0.25%, but that the variance is large with a 95-percentile range of 0 <= fsub <= 1%. We quantify possible effects of substructure on quadruply-imaged lens systems using the cusp relation and the simple statistic, Rcusp. We estimate that the probability of obtaining the large values of the Rcusp which have been observed from substructure effects is roughly ~10^-3 to ~10^-2. We consider a variety of possible correlations between host halo properties and substructure properties in order to probe possible sample biases. In particular, low-concentration host dark matter halos have more large substructures and give rise to large values of Rcusp more often. However, there is no known observational bias that would drive observed quadruply-imaged quasars to be produced by low-concentration lens halos. Finally, we show that the substructure mass fraction is a relatively reliable predictor of the value of Rcusp.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We model the evolution of galaxy clustering through cosmic time to investigate the nature of the power-law shape of xi(r), the galaxy two-point correlation function. Departures from a power law are mainly governed by galaxy pair counts on small scales, subject to non-linear dynamics. We assume that galaxies reside within dark matter halos and subhalos and use a semi-analytic substructure evolution model to study subhalo populations within host halos. We find that tidal mass loss and, to a lesser degree, dynamical friction deplete the number of subhalos within larger host halos over time by ~90%, just the right amount for achieving a power-law xi(r) at z = 0. We find that xi(r) breaks from a power law at high masses, implying that only galaxies of luminosities <= Lstar should exhibit power-law clustering. We also demonstrate that xi(r) evolves from being far from a power law at high redshift, toward a near power law at z = 0 and deviates in the future. This is mainly caused by the evolving competition between the accretion and destruction rates of subhalos, which happen to strike just the right balance at z~0. We show that key ingredients determining the shape of xi(r) are the fraction of galaxies that are satellites, the relative difference in mass between the halos of isolated galaxies and halos that contain a single satellite on average, and the rareness of halos that host galaxies. These pieces are intertwined and we find no simple, universal rule for which a power-law xi(r) will occur. However, we do show that the physics responsible for setting the galaxy content of halos do not care about the conditions needed to achieve a power law xi(r) and these conditions are met only in a narrow mass and redshift range. We conclude that the power-law nature of xi(r) for Lstar and fainter galaxy samples at low redshift is a cosmic coincidence. (Abridged)

We present analysis of the three-dimensional shape of intracluster gas in clusters formed in cosmological simulations of the Lambda-CDM cosmology and compare it to the shape of dark matter distribution and the shape of the overall isopotential surfaces. We find that in simulations with radiative cooling, star formation and stellar feedback (CSF), intracluster gas outside the cluster core is more spherical compared to non-radiative (NR) simulations, while in the core the gas in the CSF runs is more triaxial and has a distinctly oblate shape. The latter reflects the ongoing cooling of gas, which settles into a thick oblate ellipsoid as it loses thermal energy. The shape of the gas in the inner regions of clusters can therefore be a useful diagnostic of gas cooling. We find that gas traces the shape of the underlying potential rather well outside the core, as expected in hydrostatic equilibrium. At smaller radii, however, the gas and potential shapes differ significantly. In the CSF runs, the difference reflects the fact that gas is partly rotationally supported. Interestingly, we find that in NR simulations the difference between gas and potential shape at small radii is due to random gas motions, which make the gas distribution more spherical than the equipotential surfaces. Finally, we use mock Chandra X-ray maps to show that the differences in shapes observed in three-dimensional distribution of gas are discernible in the ellipticity of X-ray isophotes. Contrasting the ellipticities measured in simulated clusters against observations can therefore constrain the amount of cooling of the intracluster medium and the presence of random gas motions in cluster cores.

Upcoming imaging surveys such as the Large Synoptic Survey Telescope will repeatedly scan large areas of sky and have the potential to yield million-supernova catalogs. Type Ia supernovae are excellent standard candles and will provide distance measures that suffice to detect mean pairwise velocities of their host galaxies. We show that when combining these distance measures with photometric redshifts for either the supernovae or their host galaxies, the mean pairwise velocities of the host galaxies will provide a dark energy probe which is competitive with other widely discussed methods. Adding information from this test to type Ia supernova photometric luminosity distances from the same experiment, plus the cosmic microwave background power spectrum from the Planck satellite, improves the Dark Energy Task Force Figure of Merit by a factor of 1.8. Pairwise velocity measurements require no additional observational effort beyond that required to perform the traditional supernova luminosity distance test, but may provide complementary constraints on dark energy parameters and the nature of gravity. Incorporating additional spectroscopic redshift follow-up observations could provide important dark energy constraints from pairwise velocities alone. Mean pairwise velocities are much less sensitive to systematic redshift errors than the luminosity distance test or weak lensing techniques, and also are only mildly affected by systematic evolution of supernova luminosity.

The nature of the dark matter remains a mystery. The possibility of an unstable dark matter particle decaying to invisible daughter particles has been explored many times in the past few decades. Meanwhile, weak gravitational lensing shear has gained a lot of attention as a probe of dark energy. Weak lensing is a useful tool for constraining the stability of the dark matter. In the coming decade a number of large, galaxy imaging surveys will be undertaken and will measure the statistics of cosmological weak lensing with unprecedented precision. Weak lensing statistics are sensitive to unstable dark matter in at least two ways. Dark matter decays alter the matter power spectrum and change the angular diameter distance-redshift relation. We show how measurements of weak lensing shear correlations may provide the most restrictive, model-independent constraints on the lifetime of unstable dark matter. Our results rely on assumptions regarding nonlinear evolution of density fluctuations in scenarios of unstable dark matter and one of our aims is to stimulate interest in theoretical work on nonlinear structure growth in unstable dark matter models.

Environmental statistics provide a necessary means of comparing the properties of galaxies in different environments and a vital test of models of galaxy formation within the prevailing, hierarchical cosmological model. We explore counts-in-cylinders, a common statistic defined as the number of companions of a particular galaxy found within a given projected radius and redshift interval. Galaxy distributions with the same two-point correlation functions do not necessarily have the same companion count distributions. We use this statistic to examine the environments of galaxies in the Sloan Digital Sky Survey, Data Release 4. We also make preliminary comparisons to four models for the spatial distributions of galaxies, based on N-body simulations, and data from SDSS DR4 to study the utility of the counts-in-cylinders statistic. There is a very large scatter between the number of companions a galaxy has and the mass of its parent dark matter halo and the halo occupation, limiting the utility of this statistic for certain kinds of environmental studies. We also show that prevalent, empirical models of galaxy clustering that match observed two- and three-point clustering statistics well fail to reproduce some aspects of the observed distribution of counts-in-cylinders on 1, 3 and 6-Mpc/h scales. All models that we explore underpredict the fraction of galaxies with few or no companions in 3 and 6-Mpc/h cylinders. Roughly 7% of galaxies in the real universe are significantly more isolated within a 6 Mpc/h cylinder than the galaxies in any of the models we use. Simple, phenomenological models that map galaxies to dark matter halos fail to reproduce high-order clustering statistics in low-density environments.

We calculate the probability distribution function (PDF) of the expected annihilation luminosities of dark matter subhalos as a function of subhalo mass and distance from the Galactic center using a semi-analytical model of halo evolution. We find that the PDF of luminosities is relatively broad, exhibiting a spread of as much as an order of magnitude at fixed subhalo mass and halo-centric distance. The luminosity PDF allows for simple construction of mock samples of gamma-ray luminous subhalos and assessment of the variance in among predicted gamma-ray signals from dark matter annihilation. Other applications include quantifying the variance among the expected luminosities of dwarf spheroidal galaxies, assessing the level at which dark matter annihilation can be a contaminant in the expected gamma-ray signal from other astrophysical sources, as well as estimating the level at which nearby subhalos can contribute to the antimatter flux.

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.

We study the impact of theoretical uncertainty in the dark matter halo mass function and halo bias on dark energy constraints from imminent galaxy cluster surveys. We find that for an optical cluster survey like the Dark Energy Survey, the accuracy required on the predicted halo mass function to make it an insignificant source of error on dark energy parameters is ~ 1%. The analogous requirement on the predicted halo bias is less stringent (~ 5%), particularly if the observable-mass distribution can be well constrained by other means. These requirements depend upon survey area but are relatively insensitive to survey depth. The most stringent requirements are likely to come from a survey over a significant fraction of the sky that aims to observe clusters down to relatively low mass, Mth ~ 10^13.7 Msun/h; for such a survey, the mass function and halo bias must be predicted to accuracies of ~ 0.5% and ~ 1%, respectively. These accuracies represent a limit on the practical need to calibrate ever more accurate halo mass and bias functions. We find that improving predictions for the mass function in the low-redshift and low-mass regimes is the most effective way to improve dark energy constraints.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

A potential flux of high-energy neutrinos from the annihilation of dark matter particles trapped within the Sun has been exploited to place indirect limits on particle dark matter. In most models, the dark matter interacts weakly, but the possibility of a dark matter particle with a large cross section for elastic scattering with itself has been proposed in several contexts. I study the consequences of such dark matter self-interactions for the high-energy neutrino flux from annihilation within the Sun. The self-interaction may allow dark matter in the halo to be captured within the Sun by scattering off of previously-captured dark matter within the Sun. This effect is not negligible in acceptable and accessible regions of parameter space. Enhancements in the predicted high-energy neutrino flux from the Sun of tens to hundreds of percent can be realized in broad regions of parameter space. Enhancements as large as factors of several hundred may be realized in extreme regions of the viable parameter space. Large enhancements require the dark matter annihilation cross section to be relatively small, sigma*v <~ 10^-27 cm^3/s. This phenomenology is interesting. First, self-capture is negligible for the Earth, so dark matter self-interactions break the correspondence between the solar and terrestrial neutrino signals. Likewise, the correspondence between indirect and direct detection limits on scattering cross sections on nuclei is broken by the self-interaction. These broken correspondences may evince strong dark matter self-interactions. In some cases, self-capture can lead to observable indirect signals in regions of parameter space where limits on dark matter-nucleon cross sections from direct detection would indicate that no such signal should be observable.

(Abridged) We perform dissipationless N-body simulations to elucidate the dynamical response of thin disks to bombardment by cold dark matter (CDM) substructure. Our method combines (1) cosmological simulations of the formation of Milky Way (MW)-sized CDM halos to derive the properties of substructure and (2) controlled numerical experiments of consecutive subhalo impacts onto an initially-thin, fully-formed MW type disk galaxy. The present study is the first to account for the evolution of satellite populations over cosmic time in such an investigation of disk structure. We find that accretions of massive subhalos onto the central regions of host halos, where the galactic disks reside, since z~1 should be common. One host halo accretion history is used to initialize the controlled simulations of satellite-disk encounters. We show that these accretion events severely perturb the thin galactic disk and produce a wealth of distinctive dynamical signatures on its structure and kinematics. These include (1) considerable thickening and heating at all radii, with the disk thickness and velocity ellipsoid nearly doubling at the solar radius; (2) prominent flaring associated with an increase in disk thickness greater than a factor of 4 in the disk outskirts; (3) surface density excesses at large radii, beyond ~5 disk scale lengths, resembling those of observed antitruncated disks; (4) lopsidedness at levels similar to those measured in observational samples of disk galaxies; and (5) substantial tilting. The interaction with the most massive subhalo drives the disk response while subsequent bombardment is much less efficient at disturbing the disk. We conclude that substructure-disk encounters of the kind expected in the LCDM paradigm play a significant role in setting the structure of disk galaxies and driving galaxy evolution.

The formation and evolution of Black Holes inevitably affects the distribution of dark and baryonic matter in the neighborhood of the Black Hole. These effects may be particularly relevant around Supermassive and Intermediate Mass Black Holes (IMBHs), the formation of which can lead to large Dark Matter overdensities, called {\em spikes} and {\em mini-spikes} respectively. Despite being larger and more dense, spikes evolve at the very centers of galactic halos, in regions where numerous dynamical effects tend to destroy them. Mini-spikes may be more likely to survive, and they have been proposed as worthwhile targets for indirect Dark Matter searches. We review here the formation scenarios and the prospects for detection of mini-spikes, and we present new estimates for the abundances of mini-spikes to illustrate the sensitivity of such predictions to cosmological parameters and uncertainties regarding the astrophysics of Black Hole formation at high redshift. We also connect the IMBHs scenario to the recent measurements of cosmic-ray electron and positron spectra by the PAMELA, ATIC, H.E.S.S., and Fermi collaborations.

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.

The detection of byproducts from particle annihilations in galactic halos would provide important information about the nature of the dark matter. Observational evidence for a local excess of high-energy positrons has motivated recent models with an additional interaction between dark matter particles that can result in a Sommerfeld enhancement to the cross section for annihilation. In such models, the cross section becomes velocity-dependent and may enhance the dark matter annihilation rate in the solar neighborhood relative to the rate in the early universe sufficiently to source observed fluxes of high-energy positrons. We demonstrate that, for particle interaction cross sections that increase with decreasing velocity, the kinematical structures of dark matter halos with interior density profiles shallower than isothermal, such as Navarro-Frenk-White or Einasto halos, may induce a further enhancement owing to the position-dependent velocity distribution. We provide specific examples for the increase in the annihilation rate with a cross section enhanced by the Sommerfeld effect. In dark matter halos like that of the Milky Way and Local Group dwarf galaxies, the effective cross section at the halo center can be significantly larger than its local value. The additional enhancement owing to halo kinematics depends upon the parameters of any model, but is a prediction of certain models aimed at explaining measured positron fluxes and can exceed an order of magnitude.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

Our heuristic understanding of the abundance of dark matter halos centers around the concept of a density threshold, or "barrier", for gravitational collapse. If one adopts the ansatz that regions of the linearly evolved density field smoothed on mass scale M with an overdensity that exceeds the barrier will undergo gravitational collapse into halos of mass M, the corresponding abundance of such halos can be estimated simply as a fraction of the mass density satisfying the collapse criterion divided by the mass M. The key ingredient of this ansatz is therefore the functional form of the collapse barrier as a function of mass M or, equivalently, of the variance sigma^2(M). Several such barriers based on the spherical, Zel'dovich, and ellipsoidal collapse models have been extensively discussed. Using large scale cosmological simulations, we show that the relation between the linear overdensity and the mass variance for regions that collapse to form halos by the present epoch resembles expectations from dynamical models of ellipsoidal collapse. However, we also show that using such a collapse barrier with the excursion set ansatz predicts a halo mass function inconsistent with that measured directly in cosmological simulations. This inconsistency demonstrates a failure of the excursion set ansatz as a physical model for halo collapse. We discuss implications of our results for understanding the collapse epoch for halos as a function of mass, and avenues for improving consistency between analytical models for the collapse epoch and the results of cosmological simulations.

We study the utility of a large sample of type Ia supernovae that might be observed in an imaging survey that rapidly scans a large fraction of the sky for constraining dark energy. We consider information from the traditional luminosity distance test as well as the spread in SNeIa fluxes at fixed redshift induced by gravitational lensing. We include a treatment of photometric redshift uncertainties in our analysis. Our primary result is that the information contained in the mean distance moduli of SNeIa and the dispersion among SNeIa distance moduli complement each other, breaking a degeneracy between the present dark energy equation of state and its time variation without the need for a high-redshift supernova sample. To address photometric redshift uncertainties, we present dark energy constraints as a function of the size of an external set of spectroscopically-observed SNeIa that may be used for redshift calibration, nspec. We find that an imaging survey can constrain the dark energy equation of state at the epoch where it is best constrained with a 1-sigma error of sigma(wpiv)~0.03-0.09$, depending upon various assumptions. In addition, the marginal improvement in the error sigma(wpiv) from an increase in the spectroscopic calibration sample drops once nspec ~ 10^3. This result is important because it is of the order of the size of calibration samples likely to be compiled in the coming decade and because, for samples of this size, the spectroscopic and imaging surveys individually place comparable constraints on the dark energy equation of state. In all cases, it is best to calibrate photometric redshifts with a set of spectroscopically-observed SNeIa with relatively more objects at high redshift than the parent sample of imaging SNeIa.

We make predictions for the metallicity of diffuse stellar components in systems ranging from small spiral galaxies to rich galaxy clusters. We extend the formalism of Purcell et al. (2007), in which diffuse stellar mass is produced via galaxy disruption, and we convolve this result with the observed mass-metallicity relation for galaxies in order to analyze the chemical abundance of intrahalo light (IHL) in host halos with virial mass 10^10.5 M_sun < M_host < 10^15 M_sun. We predict a steep rise of roughly two dex in IHL metallicity from the scales of small to large spiral galaxies. In terms of the total dynamical mass M_host of the host systems under consideration, we predict diffuse light metallicities ranging from Z_IHL < -2.5 for M_host ~ 10^11 M_sun, to Z_IHL ~ -1.0 for M_host ~ 10^12 M_sun. In larger systems, we predict a gradual flattening of this trend with Z_IHL ~ -0.4 for M_host ~ 10^13 M_sun, increasing to Z_IHL ~ 0.1 for M_host ~ 10^15 M_sun. This behavior is coincident with a narrowing of the intrahalo metallicity distribution as host mass increases. The observable distinction in surface brightness between old, metal-poor IHL stars and more metal-rich, dynamically-younger tidal streams is of crucial importance when estimating the chemical abundance of an intrahalo population with multiple origins.

We perform a set of high-resolution, dissipationless N-body simulations to investigate the influence of cold dark matter (CDM) substructure on the dynamical evolution of thin galactic disks. Our method combines cosmological simulations of galaxy-sized CDM halos to derive the properties of substructure populations and controlled numerical experiments of consecutive subhalo impacts onto initially-thin, fully-formed disk galaxies. We demonstrate that close encounters between massive subhalos and galactic disks since z~1 should be common occurrences in LCDM models. In contrast, extremely few satellites in present-day CDM halos are likely to have a significant impact on the disk structure. One typical host halo merger history is used to seed controlled N-body experiments of subhalo-disk encounters. As a result of these accretion events, the disk thickens considerably at all radii with the disk scale height increasing in excess of a factor of 2 in the solar neighborhood. We show that interactions with the subhalo population produce a wealth of distinctive morphological signatures in the disk stars including: conspicuous flares; bars; low-lived, ring-like features in the outskirts; and low-density, filamentary structures above the disk plane. We compare a resulting dynamically-cold, ring-like feature in our simulations to the Monoceros ring stellar structure in the MW. The comparison shows quantitative agreement in both spatial distribution and kinematics, suggesting that such observed complex stellar components may arise naturally as disk stars are excited by encounters with subhalos. These findings highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution.

(Abridged) We conduct a series of high-resolution, dissipationless N-body simulations to investigate the cumulative effect of substructure mergers onto thin disk galaxies in the context of the LCDM paradigm of structure formation. Our simulation campaign is based on a hybrid approach. Substructure properties are culled directly from cosmological simulations of galaxy-sized cold dark matter (CDM) halos. In contrast to what can be inferred from statistics of the present-day substructure populations, accretions of massive subhalos onto the central regions of host halos, where the galactic disk resides, since z~1 should be common occurrences. One host halo merger history is subsequently used to seed controlled numerical experiments of repeated satellite impacts on an initially-thin Milky Way-type disk galaxy. We show that these accretion events produce several distinctive observational signatures in the stellar disk including: a ring-like feature in the outskirts; a significant flare; a central bar; and faint filamentary structures that (spuriously) resemble tidal streams. The final distribution of disk stars exhibits a complex vertical structure that is well-described by a standard ``thin-thick'' disk decomposition. We conclude that satellite-disk encounters of the kind expected in LCDM models can induce morphological features in galactic disks that are similar to those being discovered in the Milky Way, M31, and in other disk galaxies. These results highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution. Upcoming galactic structure surveys and astrometric satellites may be able to distinguish between competing cosmological models by testing whether the detailed structure of galactic disks is as excited as predicted by the CDM paradigm.

We study the merger histories of galaxy dark matter halos using a high resolution LCDM N-body simulation. Our merger trees follow ~17,000 halos with masses M_0 = (10^11--10^13) Msun at z=0 and track accretion events involving objects as small as m = 10^10 Msun. We find that mass assembly is remarkably self-similar in m/M_0, and dominated by mergers that are ~10% of the final halo mass. While very large mergers, m > 0.4 M_0, are quite rare, sizeable accretion events, m ~ 0.1 M_0, are common. Over the last 10 Gyr, an overwhelming majority (~95%) of Milky Way-sized halos with M_0 = 10^12 Msun have accreted at least one object with greater total mass than the Milky Way disk (m > 5x10^10 Msun), and approximately 70% have accreted an object with more than twice that mass (m > 10^11 Msun). Our results raise serious concerns about the survival of thin-disk dominated galaxies within the current paradigm for galaxy formation in a CDM universe. In order to achieve a ~70% disk-dominated fraction in Milky Way-sized CDM halos, mergers involving m ~ 2x10^11 Msun objects must not destroy disks. Considering that most thick disks and bulges contain old stellar populations, the situation is even more restrictive: these mergers must not heat disks or drive gas into their centers to create young bulges.

Numerical studies indicate that uncertainties in the treatment of baryonic physics can affect predictions for shear power spectra at a level that is significant for forthcoming surveys such as DES, SNAP, and LSST. Correspondingly, we show that baryonic effects can significantly bias dark energy parameter measurements. Eliminating such biases by neglecting information in multipoles beyond several hundred leads to weaker parameter constraints by a factor of approximately 2 to 3 compared with using information out to multipoles of several thousand. Fortunately, the same numerical studies that explore the influence of baryons indicate that they primarily affect power spectra by altering halo structure through the relation between halo mass and mean effective halo concentration. We explore the ability of future weak lensing surveys to constrain both the internal structures of halos and the properties of the dark energy simultaneously as a first step toward self calibrating for the physics of baryons. This greatly reduces parameter biases and no parameter constraint is degraded by more than 40% in the case of LSST or 30% in the cases of SNAP or DES. Modest prior knowledge of the halo concentration relation greatly improves even these forecasts. Additionally, we find that these surveys can constrain effective halo concentrations near m~10^14 Msun/h and z~0.2 to better than 10% with shear power spectra alone. These results suggest that inferring dark energy parameters with measurements of shear power spectra can be made robust to baryonic effects and may simultaneously be competitive with other methods to inform models of galaxy formation. (Abridged)

The lensing cross section of triaxial halos depends on the relative orientation between a halo's principal axes and its line of sight. Consequently, a lensing subsample of randomly oriented halos is not, in general,randomly oriented. Using an isothermal mass model for the lensing galaxies and their host halos, we show that the lensing subsample of halos that produces doubles is preferentially aligned along the lines of sight, whereas halos that produce quads tend to be projected along their middle axes. These preferred orientations result in different projected ellipticity distributions for quad, doubles, and random galaxies. We show that ~300 lens systems must be discovered to detect this effect at the 95% confidence level. We also investigate the importance of halo shape for predicting the quad-to-double ratio and find that the latter depends quite sensitively on the distribution of the short-to-long axis ratio, but is otherwise nearly independent of halo shape. Finally, we estimate the impact of the preferred orientation of lensing galaxies on their projected substructure mass fraction, and find that the observed alignment between the substructure distribution and the mass distribution of halos result in a negligible bias.

Galaxy pairs provide a potentially powerful means of studying triggered star formation from galaxy interactions. We use a large cosmological N-body simulation coupled with a well-tested semi-analytic substructure model to demonstrate that the majority of galaxies in close pairs reside within cluster or group-size halos and therefore represent a biased population, poorly suited for direct comparison to ``field'' galaxies. Thus, the frequent observation that some types of galaxies in pairs have redder colors than ``field'' galaxies is primarily a selection effect. We select galaxy pairs that are isolated in their dark matter halos with respect to other massive subhalos (N=2 halos) and a control sample of isolated galaxies (N=1 halos) for comparison. We then apply these selection criteria to a volume-limited subset of the 2dF Galaxy Redshift Survey with M_Bj <= -19 and obtain the first clean measure of the typical fraction of galaxies affected by triggered star formation and the average elevation in the star formation rate. We find that 24% (30.5%) of these L^\star and sub-L^{\star} galaxies in isolated 50 (30) kpc/h pairs exhibit star formation that is boosted by a factor of >~ 5 above their average past value, while only 10% of isolated galaxies in the control sample show this level of enhancement. Thus, 14% (20 %) of the galaxies in these close pairs show clear triggered star formation. The isolation criteria we develop provide a means to constrain star formation and feedback prescriptions in hydrodynamic simulations and a very general method of understanding the importance of triggered star formation in a cosmological context. (Abridged.)

The discovery of luminous quasars at redshift z ~ 6 indicates the presence of supermassive black holes (SMBHs) of mass ~10^9 Msun when the Universe was less than one billion years old. This finding presents several challenges for theoretical models. Here, we present the first multi-scale simulations that, together with a self-regulated model for the SMBH growth, produce a luminous quasar at z ~ 6.5 in the LCDM paradigm. We follow the hierarchical assembly history of the most massive halo in a ~ 3 Gpc^3 volume, and find that this halo of ~ 8x 10^{12} Msun forming at z ~ 6.5 after several major mergers is able to reproduce a number of observed properties of SDSS J1148+5251, the most distant quasar detected at z =6.42 (Fan et al. 2003). Moreover, the SMBHs grow through gas accretion below the Eddington limit in a self-regulated manner owing to feedback. We find that the progenitors experience significant star formation (up to 10^4 Msun/yr) preceding the major quasar phase such that the stellar mass of the quasar host reaches 10^{12} Msun at z ~ 6.5, consistent with observations of significant metal enrichment in SDSS J1148+5251. Our results provide a viable formation mechanism for z ~ 6 quasars in the standard LCDM cosmology, and demonstrate a common, merger-driven origin for the rarest quasars and the fundamental SMBH-host correlation in a hierarchical Universe.(Abridged)

We make predictions for diffuse stellar mass fractions in dark matter halos from the scales of small spiral galaxies to those of large galaxy clusters. We use an extensively-tested analytic model for subhalo infall and evolution and empirical constraints from galaxy survey data to set the stellar mass in each accreted subhalo to model diffuse light. We add stellar mass to the diffuse light as subhalos become disrupted due to interactions within their host halos. We predict that the stellar mass fraction in diffuse, intrahalo light should rise on average from ~0.5% to approximately 20% from small galaxy halos to poor groups. The trend with mass flattens considerably beyond the group scale, increasing weakly from a fraction of ~20% in poor galaxy clusters (~10^14 M_sun) to roughly ~30% in massive clusters (~10^15 M_sun). The mass-dependent diffuse light fraction is governed primarily by the empirical fact that the mass-to-light ratio in galaxy halos must vary as a function of halo mass. Galaxy halos have little diffuse light because they accrete most of their mass in small subhalos that themselves have high mass-to-light ratios; stellar halos around galaxies are built primarily from disrupted dwarf-irregular-type galaxies with M*~10^8.5 M_sun. The diffuse light in group and cluster halos is built from satellite galaxies that form stars efficiently and have correspondingly low mass-to-light ratios; intracluster light is dominated by material liberated from massive galaxies with M*~10^11 M_sun. Our results are consistent with existing observations spanning the galaxy, group, and cluster scale; however, they can be tested more rigorously in future deep surveys for faint diffuse light.

(Abridged) We present a semi-analytic model for Cold Dark Matter halo substructure that can be used as a framework for studying galaxy formation and an ingredient in halo models of galaxy clustering. We perform a comprehensive comparison of the model to the results of a suite of high-resolution cosmological simulations. The comparisons reveal that subhalo statistics, such as velocity and mass functions, radial distributions, and the halo occupation distributions, agree well over three orders of magnitude in host halo mass at various redshifts. Both in the simulations and in our model, the radial distributions of subhalos are significantly shallower than that of the dark matter density. The abundance of subhalos in a host is set by competition between destruction and new accretion. High mass halos and halos at high redshift tend to host more subhalos because the subhalos have, on average, been accreted more recently. Similarly, at fixed mass, halos that formed more recently host more subhalos. Observed "fossil groups" may represent an extreme tail of this correlation. We find a related correlation between host halo concentration and subhalo abundance at fixed host mass, Nsat ~ c^-a, where a changes with redshift and host-to-subhalo mass ratio. Lastly, we use our model to populate host halos in one of our high-resolution cosmological simulations, replacing the subhalos resolved in the simulation with subhalos computed according to our model. We show that the resulting 2pt correlation function of such a "hybrid" halo ensemble is indistinguishable from that measured directly in the simulation. This supports one of the key tenets of the standard halo model -- the assumption that the halo occupation distribution is statistically independent of host halo environment.

We study the importance of baryonic physics on predictions of the matter power spectrum as it is relevant for forthcoming weak lensing surveys. We quantify the impact of baryonic physics using a set of three cosmological numerical simulations. Each simulation has the same initial density field, but models a different set of physical processes. The first simulation evolves the density field using gravity alone, the second includes non-radiative gasdynamics, and the third includes radiative heating and cooling of baryons, star formation, and supernova feedback. We find that baryonic processes alter predictions for the matter power spectrum significantly relative to models that include only gravitational interactions. Our results imply that future weak lensing experiments such as LSST and SNAP will be very sensitive to the poorly-understood physics governing the nonlinear evolution of the baryonic component of the universe. The net effect is significantly larger in the case of the model with cooling and star formation, in which case our results imply that contemporary surveys such as the CFHT Wide survey may also be sensitive to baryonic processes. In particular, this effect could be important for forecasts of the constraining power of future surveys if information from scales larger than l ~ 1000 is included in the analysis. We find that deviations are caused primarily by the rearrangement of matter within individual dark matter halos relative to the gravity-only case, rather than a large-scale rearrangement of matter. Consequently, we propose a simple model, based on the phenomenological halo model of dark matter clustering, for baryonic effects that can be used to aid in the interpretation of forthcoming weak lensing data.

I review the excursion set theory (EST) of dark matter halo formation and clustering. I recount the Press-Schechter argument for the mass function of bound objects and review the derivation of the Press-Schechter mass function in EST. The EST formalism is powerful and can be applied to numerous problems. I review the EST of halo bias and the properties of void regions. I spend considerable time reviewing halo growth in the EST. This section culminates with descriptions of two Monte Carlo methods for generating halo mass accretion histories. In the final section, I emphasize that the standard EST approach is the result of several simplifying assumptions. Dropping these assumptions can lead to more faithful predictions and a more versatile formalism. One such assumption is the constant height of the barrier for nonlinear collapse. I review implementations of the excursion set approach with arbitrary barrier shapes. An application of this is the now well-known improvement to standard EST that follows from the ellipsoidal-collapse barrier. Additionally, I emphasize that the statement that halo accretion histories are independent of halo environments is a simplifying assumption, rather than a prediction of the theory. I review the method for constructing correlated random walks of the density field in more general cases. I construct a simple toy model with correlated walks and I show that excursion set theory makes a qualitatively simple and general prediction for the relation between halo accretion histories and halo environments: regions of high density preferentially contain late-forming halos and conversely for regions of low density. I conclude with a brief discussion of this prediction in the context of recent numerical studies of the environmental dependence of halo properties. (Abridged)

Standard LCDM predicts that the major merger rate of galaxy-size dark matter halos rises rapidly with redshift. The average number of close companions per galaxy, Nc, is often used to infer the galaxy merger rate, however, recent observational studies suggest that Nc evolves very little with redshift. Here we use a "hybrid" N- body simulation plus analytic substructure model to predict Nc directly. We identify dark matter subhalos with galaxies and show that the observed lack of close pair count evolution arises because the high merger rate per halo at early times is counteracted by a decrease in the number of halos massive enough to host a galaxy pair. We compare our results to data compiled from the DEEP2, SSRS2, and the UZC redshift surveys. Observed pair counts match our predictions if we assume a monotonic mapping between galaxy luminosity and the maximum circular velocity that each subhalo had when it was first accreted onto its host halo. This suggests that satellite galaxies are significantly more resilient to destruction than are dissipationless dark matter subhalos. We argue that while Nc does not provide a direct measure of the halo merger rate, it offers a powerful means to constrain the Halo Occupation Distribution and the spatial distribution of galaxies within halos. Interpreted in this way, close pair counts provide a useful test of galaxy formation processes on < 100 kpc scales.

We explore the angular distribution of two samples of satellite galaxies orbiting isolated hosts extracted from the Sloan Digital Sky Survey Data Release 4. We find a clear alignment of the satellites along the major axis of their hosts when restricting the analysis to red hosts. The anisotropy is most pronounced for red satellites of red hosts. We find that the distribution of the satellites about blue, isolated hosts is consistent with isotropy. We show that under the assumption that the true, underlying distribution of satellites of blue hosts exhibits the same anisotropy as the satellites of red hosts, the sample of blue hosts is too small to measure this anisotropy at a statistically-significant level. The anisotropy that we detect for satellites about red primaries is independent of the projected radius. In particular, it is evident at large projected distances from primaries (300 < r_p < 500 kpc).

We investigate the dependence of dark matter halo clustering on halo formation time, density profile concentration, and subhalo occupation number, using high-resolution numerical simulations of a LCDM cosmology. We confirm results that halo clustering is a function of halo formation time, and that this trend depends on halo mass. For the first time, we show unequivocally that halo clustering is a function of halo concentration and show that the dependence of halo bias on concentration, mass, and redshift can be accurately parameterized in a simple way: b(c,M|z) = b(M|z) b(c|M/M*). The scaling between bias and concentration changes sign with the value of M/M*: high concentration (early forming) objects cluster more strongly for M <~ M* while low concentration (late forming) objects cluster more strongly for rare high-mass halos, M >~ M*. We show the first explicit demonstration that host dark halo clustering depends on the halo occupation number (of dark matter subhalos) and discuss implications for halo model calculations of dark matter power spectra and galaxy clustering statistics. The effect of these halo properties on clustering is strongest for early-forming dwarf-mass halos, which are significantly more clustered than typical halos of their mass. Our results suggest that isolated low-mass galaxies (e.g. low surface-brightness dwarfs) should have more slowly-rising rotation curves than their clustered counterparts, and may have consequences for the dearth of dwarf galaxies in voids. They also imply that self calibrating richness-selected cluster samples with their clustering properties might overestimate cluster masses and bias cosmological parameter estimation.

Precise proper motion measurements (sigma_mu ~ 10 mkas/yr) of the recently discovered hyper-velocity star (HVS) SDSS J090745.0+024507 would yield significant constraints on the axis ratios and orientation of a triaxial model for the Galactic halo. Triaxiality of dark matter halos is predicted by Cold Dark Matter models of galaxy formation and may be used to probe the nature of dark matter. However, unless the distance to this star is determined to better than 10%, these constraints suffer from one-dimensional degeneracies, which we quantify. We show how proper motion measurements of several HVSs could simultaneously resolve the distance degeneracies of all such stars and produce a detailed picture of the triaxial halo. Additional HVSs may be found from radial velocity surveys or from parallax/proper-motion data derived from GAIA. High-precision proper-motion measurements of these stars using the Space Interferometry Mission (SIM PlanetQuest) would substantially tighten the constraints they yield on the Galactic potential.

We present a comprehensive series of dissipationless N-body simulations to investigate the evolution of density distribution in equal-mass mergers between dark matter (DM) halos and multicomponent galaxies. The DM halo models are constructed with various asymptotic power-law indices ranging from steep cusps to core-like profiles and the structural properties of the galaxy models are motivated by the LCDM paradigm of structure formation. The adopted force resolution allows robust density profile estimates in the inner ~1% of the virial radii of the simulated systems. We demonstrate that the central slopes and overall shapes of the remnant density profiles are virtually identical to those of the initial systems suggesting that the remnants retain a remarkable memory of the density structure of their progenitors, despite the relaxation that accompanies merger activity. We also find that halo concentrations remain approximately constant through hierarchical merging involving identical systems and show that remnants contain significant fractions of their bound mass well beyond their formal virial radii. These conclusions hold for a wide variety of initial asymptotic density slopes, orbital energies, and encounter configurations, including sequences of consecutive merger events, simultaneous mergers of severals ystems, and mergers of halos with embedded cold baryonic components in the form of disks, spheroids, or both. As an immediate consequence, the net effect of gas cooling, which contracts and steepens the inner density profiles of DM halos, should be preserved through a period of dissipationless major merging. Our results imply that the characteristic universal shape of DM density profiles may be set early in the evolution of halos.

We study the statistical properties of magnification perturbations by substructures in strong lensed systems using linear perturbation theory and an analytical substructure model including tidal truncation and a continuous substructure mass spectrum. We demonstrate that magnification perturbations are dominated by perturbers found within a tidal radius of an image, and that sizable magnification perturbations may arise from small, coherent contributions from several substructures within the lens halo. We find that the root-mean-square (rms) fluctuation of the magnification perturbation is 10% to 20% and both the average and rms perturbations are sensitive to the mass spectrum and density profile of the perturbers. Interestingly, we find that relative to a smooth model of the same mass, the average magnification in clumpy models is lower (higher) than that in smooth models for positive (negative) parity images. This is opposite from what is observed if one assumes that the image magnification predicted by the best-fit smooth model of a lens is a good proxy for what the observed magnification would have been if substructures were absent. While it is possible for this discrepancy to be resolved via nonlinear perturbers, we argue that a more likely explanation is that the assumption that the best-fit lens model is a good proxy for the magnification in the absence of substructure is not correct. We conclude that a better theoretical understanding of the predicted statistical properties of magnification perturbations by CDM substructure is needed in order to affirm that CDM substructures have been unambiguously detected.

We study the gravitational wave background produced from the formation and assembly of supermassive black holes within the cosmological paradigm of hierarchical structure formation. In particular, we focus on a supermassive black hole formation scenario in which the present-day population of supermassive black holes is built from high-mass seed black holes and we compute the concomitant spectrum of gravitational radiation produced by mergers of the seed black holes. We find that this scenario predicts a large, gravitational wave background that should be resolved into individual sources with space interferometers such as the proposed Laser Interferometric Space Antenna (LISA). The number of inspiral, merger and ringdown events above a signal to noise ratio of 5 that result from massive black hole seeds is of order 10^3. This prediction is robust and insensitive to several of the details of the model. We conclude that an interferometer such as LISA will be able to effectively rule out or confirm a class of models where supermassive black holes grow from high-mass seed black holes and may be able to place strong limits on the role of mergers as a channel for supermassive black hole growth. Supermassive black hole seeds likely form in the earliest proto-galactic structures at high redshift and the masses of supermassive black holes are known to be strongly correlated with the potentials of the spheroids in which they reside, therefore these results imply that space interferometers can be used as a powerful probe of the physics of galaxy formation and proto-galaxy formation at very high redshift.

I review recent results on the spatial distribution of substructure in CDM halos. I show that the spatial distribution of subhalos is anisotropic and generally prolate with a long axis that is closely aligned with the long axis of the total mass distribution of the host halo. I show that the correlation between the subhalo distribution and the long axis of the host halo is strong both in dissipationless and dissipational gasdynamical simulations. More massive subhalos tend to be more strongly clustered along the major axis of the host halo reflecting filamentary accretion. The anisotropy of subhalos has potential implications for the interpretation of several observations in the Local Group and beyond. For example, I show that while the mean projected mass fraction in substructure in the central regions of CDM halos is fsub \approx 0.4%, fsub is a strong function of projection angle and is ~ 5 to 6 times higher for projections that are nearly collinear with the major axis of the host halo.

We study the prospects for detecting gamma-rays from Dark Matter (DM) annihilations in enhancements of the DM density (mini-spikes) around intermediate-mass black holes with masses in the range $10^2 \lsim M / \msun \lsim 10^6$. Focusing on two different IMBH formation scenarios, we show that, for typical values of mass and cross section of common DM candidates, mini-spikes, produced by the adiabatic growth of DM around pregalactic IMBHs, would be bright sources of gamma-rays, which could be easily detected with large field-of-view gamma-ray experiments such as GLAST, and further studied with smaller field-of-view, larger-area experiments like Air Cherenkov Telescopes CANGAROO, HESS, MAGIC and VERITAS. The detection of many gamma-ray sources not associated with a luminous component of the Local Group, and with identical cut-offs in their energy spectra at the mass of the DM particle, would provide a potential smoking-gun signature of DM annihilations and shed new light on the nature of intermediate and supermassive Black Holes.

Observational evidence indicates a mismatch between the shapes of collisionless dark matter (DM) halos and those of observed systems. Using hydrodynamical cosmological simulations we investigate the effect of baryonic dissipation on halo shapes. We show that dissipational simulations produce significantly rounder halos than those formed in equivalent dissipationless simulations. Gas cooling causes an average increase in halo principal axis ratios of ~ 0.2-0.4 in the inner regions and a systematic shift that persists out to the virial radius, alleviating any tension between theory and observations. Although the magnitude of the effect may be overestimated due to overcooling, cluster formation simulations designed to reproduce the observed fraction of cold baryons still produce substantially rounder halos. Subhalos also exhibit a trend of increased axis ratios in dissipational simulations. Moreover, we demonstrate that subhalos are generally rounder than corresponding field halos even in dissipationless simulations. Lastly, we analyze a series of binary, equal-mass merger simulations of disk galaxies. Collisionless mergers reveal a strong correlation between DM halo shape and stellar remnant morphology. In dissipational mergers, the combination of strong gas inflows and star formation leads to an increase of the DM axis ratios in the remnant. All of these results highlight the vital role of baryonic processes in comparing theory with observations and warn against over-interpreting discrepancies with collisionless simulations on small scales.

We analyze the angular distribution and the orbital rotation directions of a sample of carefully-selected satellite galaxies extracted from the Sloan Digital Sky Survey (SDSS). We also study these statistics in an N-body simulation of cosmological structure formation set within the LCDM paradigm under various assumptions for the orientations of disk angular momenta. Assuming that the angular momenta of the disks are aligned with the angular momenta of the inner regions of their host dark matter halos, we find that the fraction of simulated satellite halos that exhibit prograde motion is 0.55-0.60, with larger satellites more likely to be prograde. In our observational sample, approximately 60% of the satellites exhibit prograde motion, a result that is broadly consistent with the simulated sample. Contrary to several recent studies, our observational sample of satellite galaxies show no evidence for being anisotropically distributed about their primary disks. Again, this result is broadly consistent with our simulated sample of satellites under the assumption that disk and halo angular momenta are aligned. However, the small size of our observational sample does not yet allow us to distinguish between various assumptions regarding the orientations of disks in their halos. Finally, we assessed the importance of contamination by interlopers on the measured prograde and retrograde statistics.

We present a study of the spatial distribution of subhalos in galactic dark matter halos using dissipationless cosmological simulations of the concordance LCDM model. We find that subhalos are distributed anisotropically and are preferentially located along the major axes of the triaxial mass distributions of their hosts. The Kolmogorov-Smirnov probability for drawing our simulated subhalo sample from an isotropic distribution is P_KS \simeq 1.5 \times 10^{-4}. An isotropic distribution of subhalos is thus not the correct null hypothesis for testing the CDM paradigm. The nearly planar distribution of observed Milky Way (MW) satellites is marginally consistent (probability \simeq 0.02) with being drawn randomly from the subhalo distribution in our simulations. Furthermore, if we select the subhalos likely to be luminous, we find a distribution that is consistent with the observed MW satellites. In fact, we show that subsamples of the subhalo population with a centrally-concentrated radial distribution, similar to that of the MW dwarfs, typically exhibit a comparable degree of planarity. We explore the origin of the observed subhalo anisotropy and conclude that it is likely due to (1) preferential accretion of subhalos along filaments, often closely aligned with the major axis of the host halo, and (2) evolution of satellite orbits within the prolate, triaxial potentials typical of CDM halos. Agreement between predictions and observations requires the major axis of the outer dark matter halo of the Milky Way to be nearly perpendicular to the disk. We discuss possible observational tests of such disk-halo alignment with current large galaxy surveys.

In this proceeding, we briefly review three recent results. First, we show that halos formed in simulations with gas cooling are significantly rounder than halos formed in dissipationless $N$-body simulations. The increase in principle axis ratios is $\delta (c/a) ~ 0.2 - 0.4$ in the inner halo and remains significant at large radii. Second, we discuss the CDM substructure crisis and demonstrate the sensitivity of the crisis to the spectrum of primordial density fluctuations on small scales. Third, we assess the ability of experiments like VERITAS and GLAST to detect $\gamma$-rays from neutralino dark matter annihilation in dark subhalos about the MW.

We investigate the effect of dark energy on the density profiles of dark matter haloes with a suite of cosmological N-body simulations and use our results to test analytic models. We consider constant equation of state models, and allow both w>-1 and w<-1. Using five simulations with w ranging from -1.5 to -0.5, and with more than ~1600 well-resolved haloes each, we show that the halo concentration model of Bullock et al. (2001) accurately predicts the median concentrations of haloes over the range of w, halo masses, and redshifts that we are capable of probing. We find that the Bullock et al. (2001) model works best when halo masses and concentrations are defined relative to an outer radius set by a cosmology-dependent virial overdensity. For a fixed power spectrum normalization and fixed-mass haloes, larger values of w lead to higher concentrations and higher halo central densities, both because collapse occurs earlier and because haloes have higher virial densities. While precise predictions of halo densities are quite sensitive to various uncertainties, we make broad comparisons to galaxy rotation curve data. At fixed power spectrum normalization (fixed sigma_8), w>-1 quintessence models seem to exacerbate the central density problem relative to the standard w=-1 model. Meanwhile w<-1 models help to reduce the apparent discrepancy. We confirm that the Jenkins et al. (2001) halo mass function provides an excellent approximation to the abundance of haloes in our simulations and extend its region of validity to include models with w<-1.

We analyze the effect of dissipation on the shapes of dark matter (DM) halos using high-resolution cosmological gasdynamics simulations of clusters and galaxies in the LCDM cosmology. We find that halos formed in simulations with gas cooling are significantly more spherical than corresponding halos formed in adiabatic simulations. Gas cooling results in an average increase of the principle axis ratios of halos by ~ 0.2-0.4 in the inner regions. The systematic difference decreases slowly with radius but persists almost to the virial radius. We argue that the differences in simulations with and without cooling arise both during periods of quiescent evolution, when gas cools and condenses toward the center, and during major mergers. We perform a series of high-resolution N-body simulations to study the shapes of remnants in major mergers of DM halos and halos with embedded stellar disks. In the DM halo-only mergers, the shape of the remnants depends only on the orbital angular momentum of the encounter and not on the internal structure of the halos. However, significant shape changes in the DM distribution may result if stellar disks are included. In this case the shape of the DM halos is correlated with the morphology of the stellar remnants.

We estimate the probability of detecting gamma-rays from the annihilation of neutralino dark matter in dense, central regions of Milky Way substructure. We characterize Galactic substructure statistically based on Monte Carlo realizations of the formation of a Milky Way-like halo using a semi-analytic method that has been calibrated against the results of high-resolution N-body simulations. We find that it may be possible for the upcoming experiments GLAST and VERITAS, working in concert, to detect gamma-rays from dark matter substructure if the neutralino is relatively light ($\Mchi \lsim 100 \GeV$), while for $\Mchi \gsim 500 \GeV$ such a detection would be unlikely. We perform most of our calculations within the framework of the standard $\Lambda$CDM cosmological model; however, we also investigate the robustness of our results to various assumptions and find that the probability of detection is sensitive to poorly-constrained input parameters, particularly those that characterize the primordial power spectrum. Specifically, the best-fitting power spectrum of the WMAP team, with a running spectral index, predicts roughly a factor of fifty fewer detectable subhalos compared to the standard $\LCDM$ cosmological model with scale-invariant power spectrum. We conclude that the lack of a detected gamma-ray signal gives very little information about the supersymmetric parameter space due to uncertainties associated with both the properties of substructure and cosmological parameters.

We estimate the probability of detecting gamma-rays from the annihilation of neutralino dark matter in the substructure of the Milky Way. We characterize substructure statistically based on Monte Carlo realizations of the formation of a Milky Way-like halo using semi-analytic methods that have been calibrated against N-body simulations. We find that it may be possible for the upcoming GLAST and VERITAS experiments, working in concert, to detect gamma-rays from dark matter substructure if $\Mchi \lsim 100 \GeV$, while for $\Mchi \gsim 500 \GeV$ such a detection seems unlikely. We investigate the effects of the underlying cosmological model and find that the probability of detection is sensitive to the primordial power spectrum of density fluctuations on small (galactic and sub-galactic) scales. We conclude that the lack of such a detection reveals little about the supersymmetric parameter space due to the uncertainties associated with the properties of substructure and cosmological parameters.

(ABRIDGED) We present a semi-analytic model to explore merger histories, destruction rates, and survival probabilities of substructure in dark matter halos and use it to study the substructure populations of galaxy-sized halos as a function of the power spectrum. We successfully reproduce the subhalo velocity function and radial distribution seen in N-body simulations for standard LCDM. We explore the implications of spectra with normalizations and tilts spanning sigma_8 = 0.65-1 and n = 0.8-1. We also study a running index (RI) model with dn/dlnk=-0.03, as discussed in the first year WMAP report, and several WDM models with masses m_W = 0.75, 1.5, 3.0 keV. The substructure mass fraction is relatively insensitive to the tilt and overall normalization of the power spectrum. All CDM-type models yield projected substructure mass fractions that are consistent with, but on the low side of, estimates from strong lens systems: f = 0.4-1.5% (64 percentile) in systems M_sub < 10^9 Msun. Truncated models produce significantly smaller fractions and are disfavored by lensing results. We compare our predicted subhalo velocity functions to the dwarf satellite population of the Milky Way. Assuming isotropic velocity dispersions, we find the standard n=1 model overpredicts the number of MW satellites as expected. Models with less small-scale power are more successful because there are fewer subhalos of a given circular velocity and the mapping between observed velocity dispersion and halo circular velocity is markedly altered. The RI model, or a fixed tilt with sigma_8=0.75, can account for the MW dwarfs without the need for differential feedback; however, these comparisons depend sensitively on the assumption of isotropic velocities in satellite galaxies.

In this proceeding, we present the results of a semi-analytic study of CDM substructure as a function of the primordial power spectrum. We apply our method to several tilted models in the LCDM framework with n=0.85-1.1, sigma_8=0.65-1.2 when COBE normalized. We also study a more extreme, warm dark matter-like spectrum that is sharply truncated below a scale of 10^10 h^-1 Msun. We show that the mass fraction of halo substructure is not a strong function of spectral slope, so it likely will be difficult to constrain tilt using flux ratios of gravitationally lensed quasars. On the positive side, all of our CDM-type models yield projected mass fractions in good agreement with strong lensing estimates: f \sim 1.5% at M \sim 10^8 Msun. The truncated model produces a significantly smaller fraction, f \lsim 0.3%, suggesting that warm dark matter-like spectra may be distinguished from CDM spectra using lensing. We also discuss the issue of dwarf satellite abundances, with emphasis on the cosmological dependence of the map between the observed central velocity dispersion of Milky Way satellites and the maximum circular velocities of their host halos. In agreement with earlier work, we find that standard LCDM over-predicts the estimated count of Milky Way satellites at fixed Vmax by an order of magnitude, but tilted models do better because subhalos are less concentrated. Interestingly, under the assumption that dwarfs have isotropic velocity dispersion tensors, models with significantly tilted spectra (n \lsim 0.85, sigma_8 \lsim 0.7) may under-predict the number of large Milky Way satellites with Vmax \gsim 40 km/s.

A problem with high central densities in dark halos has arisen in the context of LCDM cosmologies with scale-invariant initial power spectra. Although n=1 is often justified by appealing to the inflation scenario, inflationary models with mild deviations from scale-invariance are not uncommon and models with significant running of the spectral index are plausible. Even mild deviations from scale-invariance can be important because halo collapse times and densities depend on the relative amount of small-scale power. We choose several popular models of inflation and work out the ramifications for galaxy central densities. For each model, we calculate its COBE-normalized power spectrum and deduce the implied halo densities using a semi-analytic method calibrated against N-body simulations. We compare our predictions to a sample of dark matter-dominated galaxies using a non-parametric measure of the density. While standard n=1, LCDM halos are overdense by a factor of 6, several of our example inflation+CDM models predict halo densities well within the range preferred by observations. We also show how the presence of massive (0.5 eV) neutrinos may help to alleviate the central density problem even with n=1. We conclude that galaxy central densities may not be as problematic for the CDM paradigm as is sometimes assumed: rather than telling us something about the nature of the dark matter, galaxy rotation curves may be telling us something about inflation and/or neutrinos. An important test of this idea will be an eventual consensus on the value of sigma_8, the rms overdensity on the scale 8 h^-1 Mpc. Our successful models have values of sigma_8 approximately 0.75, which is within the range of recent determinations. Finally, models with n>1 (or sigma_8 > 1) are highly disfavored.

Although the currently favored cold dark matter plus cosmological constant model for structure formation assumes an n=1 scale-invariant initial power spectrum, most inflation models produce at least mild deviations from n=1. Because the lever arm from the CMB normalization to galaxy scales is long, even a small ``tilt'' can have important implications for galactic observations. Here we calculate the COBE-normalized power spectra for several well-motivated models of inflation and compute implications for the substructure content and central densities of galaxy halos. Using an analytic model, normalized against N-body simulations, we show that while halos in the standard (n=1) model are overdense by a factor of ~6 compared to observations, several of our example inflation+LCDM models predict halo densities well within the range of observations, which prefer models with n ~ 0.85. We go on to use a semi-analytic model (also normalized against N-body simulations) to follow the merger histories of galaxy-sized halos and track the orbital decay, disruption, and evolution of the merging substructure. Models with n ~0.85 predict a factor of ~3 fewer subhalos at a fixed circular velocity than the standard $n = 1$ case. Although this level of reduction does not resolve the ``dwarf satellite problem'', it does imply that the level of feedback required to match the observed number of dwarfs is sensitive to the initial power spectrum. Finally, the fraction of galaxy-halo mass that is bound up in substructure is consistent with limits imposed by multiply imaged quasars for all models considered: f_sub > 0.01 even for an effective tilt of n ~0.8.We conclude that, at their current level, lensing constraints of this kind do not provide strong limits on the primordial power spectrum.

We discuss bounds on the cosmological relativistic energy density as a function of redshift, reviewing the big bang nucleosynthesis and cosmic microwave background bounds, updating bounds from large scale structure, and introducing a new bound from the magnitude-redshift relation of Type Ia supernovae. We conclude that the standard and well-motivated assumption that relativistic energy is negligible during recent epochs is not necessitated by extant data. We then demonstrate the utility of these bounds by constraining the mass and lifetime of a hypothetical massive big bang relic particle.

Primordial nucleosynthesis provides a probe of the universal abundance of baryons when the universe was only a few minutes old. Recent observations of anisotropy in the cosmic microwave background (CMB) probe the baryon abundance when the universe was several hundred thousand years old. Observations of type Ia supernovae and clusters of galaxies in the very recent past, when the universe is several billion years old and older, provide a complementary measure of the baryon density in excellent agreement with the early universe values. The general agreement among the three measurements represents an impressive confirmation of the standard model of cosmology. However, there is a hint that the CMB observations may not be in perfect agreement with those from big bang nucleosynthesis (BBN). If this "tension" between BBN and the CMB persists, the standard model of cosmology may need to be modified. Here, in a contribution dedicated to Silvia Torres-Peimbert and Manuel Peimbert, we describe how an asymmetry between neutrinos and antineutrinos ("neutrino degeneracy") has the potential for resolving this possible conflict between BBN and the CMB.

Data from Type Ia supernovae, along with X-ray cluster estimates of the universal baryon fraction and Big Bang Nucleosynthesis (BBN) determinations of the baryon-to-photon ratio, are used to provide estimates of several global cosmological parameters at epochs near zero redshift. We show that our estimate of the present baryon density is in remarkably good agreement with that inferred from BBN at high redshift, provided the primordial abundance of deuterium is relatively low and the Universe is flat. We also compare these estimates to the baryon density at z = 1100 as inferred from the CMB angular power spectrum.