We present a Bayesian method to identify multiple (chemodynamic) stellar populations in dwarf spheroidal galaxies (dSphs) using velocity, metallicity, and positional stellar data without the assumption of spherical symmetry. We apply this method to a new Keck/DEIMOS spectroscopic survey of the Ursa Minor (UMi) dSph. We identify 892 likely members, making this the largest UMi sample with line-of-sight velocity and metallicity measurements. Our Bayesian method detects two distinct chemodynamic populations with high significance ($\ln{B}\sim33$). The metal-rich ($[{\rm Fe/H}]=-2.05\pm0.03$) population is kinematically colder (radial velocity dispersion of $\sigma_v=4.9\pm0.8 \, {\rm km \, s^{-1}}$) and more centrally concentrated than the metal-poor ($[{\rm Fe/H}]=-2.29\pm0.05$) and kinematically hotter population ($\sigma_v =11.5\pm0.9\, {\rm km \, s^{-1}}$). Furthermore, we apply the same analysis to an independent MMT/Hectochelle data set and confirm the existence of two chemodynamic populations in UMi. In both data sets, the metal-rich population is significantly flattened ($\epsilon=0.75\pm0.03$) and the metal-poor population is closer to spherical ($\epsilon=0.33_{-0.09}^{+0.12}$). Despite the presence of two populations, we are unable to robustly estimate the slope of the dynamical mass profile. We found hints for prolate rotation of order $\sim 2 \, {\rm km \, s^{-1}}$ in the MMT data set, but further observations are required to verify this. The flattened metal-rich population invalidates assumptions built into simple dynamical mass estimators, so we computed new astrophysical dark matter annihilation (J) and decay profiles based on the rounder, hotter metal-poor population and inferred $\log_{10}{(J(0.5^{\circ})/{\rm GeV^{2} \, cm^{-5}})}\approx19.1$ for the Keck data set. Our results paint a more complex picture of the evolution of Ursa Minor than previously discussed.

The extended excess towards the Galactic Center (GC) in gamma rays inferred from Fermi-LAT observations has been interpreted as being due to dark matter (DM) annihilation. Here, we perform new likelihood analyses of the GC and show that when including templates for the stellar galactic and nuclear bulges, the GC shows no significant detection of a DM annihilation template, even after generous variations in the Galactic diffuse emission models and a wide range of DM halo profiles. We include Galactic diffuse emission models with combinations of 3D inverse Compton maps, variations of interstellar gas maps, and a central source of electrons. For the DM profile, we include both spherical and ellipsoidal DM morphologies and a range of radial profiles from steep cusps to kiloparsec-sized cores, motivated in part by hydrodynamical simulations. Our derived upper limits on the dark matter annihilation flux place strong constraints on DM properties. In the case of the pure $b$-quark annihilation channel, our limits on the annihilation cross section are more stringent than those from the Milky Way dwarfs up to DM masses of $\sim$TeV, and rule out the thermal relic cross section up to $\sim$300 GeV. Better understanding of the DM profile, as well as the Fermi-LAT data at its highest energies, would further improve the sensitivity to DM properties.

We use a semianalytic approach that is calibrated to N-body simulations to study the evolution of self-interacting dark matter cores in galaxies. We demarcate the regime where the temporal evolution of the core density follows a well-defined track set by the initial halo parameters and the cross section. Along this track, the central density reaches a minimum value set by the initial halo density. Further evolution leads to an outward heat transfer, inducing gravothermal core collapse such that the core shrinks as its density increases. We show that the time scale for the core collapse is highly sensitive to the outer radial density profile. Satellite galaxies with significant mass loss due to tidal stripping should have larger central densities and significantly faster core collapse compared to isolated halos. Such a scenario could explain the dense and compact cores of dwarf galaxies in the Local Group like Tucana (isolated from the Milky Way), the classical Milky Way satellite Draco, and some of the ultrafaint satellites. If the ultimate fate of core collapse is black hole formation, then the accelerated time scale provides a new mechanism for creating intermediate mass black holes.

We explore the implications of a rapid appearance of dark energy between the redshifts ($z$) of one and two on the expansion rate and growth of perturbations. Using both Gaussian process regression and a parameteric model, we show that this is the preferred solution to the current set of low-redshift ($z<3$) distance measurements if $H_0=73~\rm km\,s^{-1}\,Mpc^{-1}$ to within 1\% and the high-redshift expansion history is unchanged from the $\Lambda$CDM inference by the Planck satellite. Dark energy was effectively non-existent around $z=2$, but its density is close to the $\Lambda$CDM model value today, with an equation of state greater than $-1$ at $z<0.5$. If sources of clustering other than matter are negligible, we show that this expansion history leads to slower growth of perturbations at $z<1$, compared to $\Lambda$CDM, that is measurable by upcoming surveys and can alleviate the $\sigma_8$ tension between the Planck CMB temperature and low-redshift probes of the large-scale structure.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM). The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the $\textit{Fermi}$ bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the ${\tt GALPROP}$ team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an $\textit{admixture formation}$ scenario, in which some of the relevant binaries are $\textit{primordial}$ and the rest are formed $\textit{dynamically}$, is preferred over a primordial-only formation scenario at $7.6\sigma$ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site.

We point out an anti-correlation between the central dark matter (DM) densities of the bright Milky Way dwarf spheroidal galaxies (dSphs) and their orbital pericenter distances inferred from Gaia data. The dSphs that have not come close to the Milky Way center (like Fornax, Carina and Sextans) are less dense in DM than those that have come closer (like Draco and Ursa Minor). The same anti-correlation cannot be inferred for the ultra-faint dSphs due to large scatter. Including ultra-faints, a trend that dSphs with more extended stellar distributions tend to have lower DM densities emerges. A fresh look at solutions to the too-big-to-fail problem is warranted in light of these observations.

The Maunakea Spectroscopic Explorer is a next-generation massively multiplexed spectroscopic facility currently under development in Hawaii. It is completely dedicated to large-scale spectroscopic surveys and will enable transformative science. In this white paper we summarize the science case and describe the current state of the project.

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

We present results from N-body simulations of self-interacting dark matter (SIDM) subhalos, which could host ultra-faint dwarf spheroidal galaxies, inside a Milky-Way-like main halo. We find that high-concentration subhalos are driven to gravothermal core collapse, while low-concentration subhalos develop large (kpc-sized) low-density cores, with both effects depending sensitively on the satellite's orbit and the self-interaction cross section over mass $\sigma/m$. The overall effect for $\sigma/m \gtrsim 3 \ \rm cm^2/g$ is to increase the range of inner densities, potentially explaining the observed diversity of Milky Way satellites, which include compact systems like Draco and Segue 1 that are dense in dark matter, and less dense, diffuse systems like Sextans and Crater II. We discuss possible ways of distinguishing SIDM models from collisionless dark matter models using the inferred dark matter densities and stellar sizes of the dwarf spheroidal galaxies.

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

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.

We perform high-resolution simulations of a MW-like galaxy in a self-interacting cold dark matter model with elastic cross section over mass of $1~\rm cm^2/g$ (SIDM) and compare to a model without self-interactions (CDM). We run our simulations with and without a time-dependent embedded potential to capture effects of the baryonic disk and bulge contributions. The CDM and SIDM simulations with the embedded baryonic potential exhibit remarkably similar host halo profiles, subhalo abundances and radial distributions within the virial radius. The SIDM host halo is denser in the center than the CDM host and has no discernible core, in sharp contrast to the case without the baryonic potential (core size $\sim 7 \, \rm kpc$). The most massive subhalos (with $V_{\mathrm{peak}}> 20 \, \rm km/s$) in our SIDM simulations, expected to host the classical satellite galaxies, have density profiles that are less dense than their CDM analogs at radii less than 500 pc but the deviation diminishes for less massive subhalos. With the baryonic potential included in the CDM and SIDM simulations, the most massive subhalos do not display the too-big-to-fail problem. However, the least dense among the massive subhalos in both these simulations tend to have the smallest pericenter values, a trend that is not apparent among the bright MW satellite galaxies.

We propose that the Ultra-Diffuse Galaxy (UDG) population represents a set of satellite galaxies born in $\sim10^{10}-10^{11}$ M$_\odot$ halos, similar to field dwarfs, which suffer a dramatic reduction in surface brightness due to tidal stripping and heating. This scenario is observationally motivated by the radial alignment of UDGs in Coma as well as the significant dependence of UDG abundance on cluster mass. As a test of this formation scenario, we apply a semi-analytic model describing the change in stellar mass and half-light radius of dwarf satellites, occupying either cored or cuspy halos, to cluster subhalos in the Illustris-dark simulation. Key to this model are results from simulations which indicate that galaxies in cored dark-matter halos expand significantly in response to tidal stripping and heating, whereas galaxies in cuspy halos experience limited size evolution. Our analysis indicates that a population of tidally-stripped dwarf galaxies, residing in cored halos (like those hosting low-surface brightness field dwarfs), is able to reproduce the observed sizes and stellar masses of UDGs in clusters remarkably well.

We test the ability of strong lensing data to constrain the size of a central core in the dark matter halos of galaxy clusters, using Abell 611 as a prototype. Using simulated data, we show that modeling a cluster halo with ellipticity in the gravitational potential can bias the inferred mass and concentration, which may bias the inferred central density when weak lensing or X-ray data are added. We also the highlight the possibility for spurious constraints on the core size if the radial density profile is different from the assumed model. These systematics can be ameliorated if central images are present in the data. Applying our methodology to Abell 611 and imposing a reasonable prior on the stellar mass-to-light ratio restricts the core size to be less than about 4 kpc, with a minimum reduced $\chi^2$ of 0.28 for 0."2 positional errors. Such small cores imply a constraint on the dark matter self-interaction cross section of the order of $0.1\ \mathrm{cm^2/g}$ at relative velocities of about 1500 km/s.

(Abridged) This is the Maunakea Spectroscopic Explorer 2018 book. It is intended as a concise reference guide to all aspects of the scientific and technical design of MSE, for the international astronomy and engineering communities, and related agencies. The current version is a status report of MSE's science goals and their practical implementation, following the System Conceptual Design Review, held in January 2018. MSE is a planned 10-m class, wide-field, optical and near-infrared facility, designed to enable transformative science, while filling a critical missing gap in the emerging international network of large-scale astronomical facilities. MSE is completely dedicated to multi-object spectroscopy of samples of between thousands and millions of astrophysical objects. It will lead the world in this arena, due to its unique design capabilities: it will boast a large (11.25 m) aperture and wide (1.52 sq. degree) field of view; it will have the capabilities to observe at a wide range of spectral resolutions, from R2500 to R40,000, with massive multiplexing (4332 spectra per exposure, with all spectral resolutions available at all times), and an on-target observing efficiency of more than 80%. MSE will unveil the composition and dynamics of the faint Universe and is designed to excel at precision studies of faint astrophysical phenomena. It will also provide critical follow-up for multi-wavelength imaging surveys, such as those of the Large Synoptic Survey Telescope, Gaia, Euclid, the Wide Field Infrared Survey Telescope, the Square Kilometre Array, and the Next Generation Very Large Array.

Galactic rotation curves exhibit diverse behavior in the inner regions, while obeying an organizing principle, i.e., they can be approximately described by a radial acceleration relation or the Modified Newtonian Dynamics phenomenology. We analyze the rotation curve data from the SPARC sample, and explicitly demonstrate that both the diversity and uniformity are naturally reproduced in a hierarchical structure formation model with the addition of dark matter self-interactions. The required concentrations of the dark matter halos are fully consistent with the concentration-mass relation predicted by the Planck cosmological model. The inferred stellar mass-to-light ($3.6 \mu m$) ratios scatter around $0.5 M_\odot/L_\odot$, as expected from population synthesis models, leading to a tight radial acceleration relation and baryonic Tully-Fisher relation. The inferred stellar-halo mass relation is consistent with the expectations from abundance matching. These results indicate that the inner dark matter halos of galaxies are thermalized due to the self-interactions of dark matter particles.

We model the expansion history of the Universe as a Gaussian Process and find constraints on the dark energy density and its low-redshift evolution using distances inferred from the Luminous Red Galaxy (LRG) and Lyman-alpha (Ly$\alpha$) datasets of the Baryon Oscillation Spectroscopic Survey, supernova data from the Joint Light-curve Analysis (JLA) sample, Cosmic Microwave Background (CMB) data from the Planck satellite, and local measurement of the Hubble parameter from the Hubble Space Telescope ($\mathsf H0$). Our analysis shows that the CMB, LRG, Ly$\alpha$, and JLA data are consistent with each other and with a $\Lambda$CDM cosmology, but the ${\mathsf H0}$ data is inconsistent at moderate significance. Including the presence of dark radiation does not alleviate the ${\mathsf H0}$ tension in our analysis. While some of these results have been noted previously, the strength here lies in that we do not assume a particular cosmological model. We calculate the growth of the gravitational potential in General Relativity corresponding to these general expansion histories and show that they are well-approximated by $\Omega_{\rm m}^{0.55}$ given the current precision. We assess the prospects for upcoming surveys to measure deviations from $\Lambda$CDM using this model-independent approach.

We consider a self-interacting dark matter model in which the massive dark photon mediating the self-interaction decays to light dark fermions to avoid over-closing the universe. We find that if the model is constrained to explain the dark matter halos inferred for spiral galaxies and galaxy clusters simultaneously, there is a strong indication that dark matter is produced asymmetrically in the early universe. It also implies the presence of dark radiation, late kinetic decoupling for dark matter, and a suppressed linear power spectrum due to dark acoustic damping. The Lyman-$\alpha$ forest power spectrum measurements put a strong upper limit on the damping scale and the model has little room to reduce the abundances of satellite galaxies. Future observations in the matter power spectrum and the CMB, in tandem with the impact of self-interactions in galactic halos, makes it possible to measure the gauge coupling and masses of the dark sector particles even when signals in conventional dark matter searches are absent.

In this paper, we conduct a search in the latest large-scale structure measurements for signatures of the dark matter-dark radiation interaction proposed by Buen-Abad et al. (2015). We show that prior claims of a detection of this interaction rely on a use of the SZ cluster mass function that ignores uncertainty in the mass-observable relationship. Including this uncertainty we find that the inferred level of interaction remains consistent with the data, but so does zero interaction; i.e., there is no longer a detection. We also point out that inference of the shape and amplitude of the matter power spectrum from Ly$\alpha$ forest measurements is highly inconsistent with the predictions of the $\Lambda$CDM model conditioned on Planck CMB temperature, polarization, and lensing power spectra, and that the dark matter-dark radiation model can restore that consistency. We also phenomenologically generalize the model of Buen-Abad et al. (2015) to allow for interaction rates with different scalings with temperature, and find that the original scaling is preferred by the data.

We present an empirical approach for interpreting gravitational wave signals of binary black hole mergers under the assumption that the underlying black hole population is sourced by remnants of stellar evolution. Using the observed relationship between galaxy mass and stellar metallicity, we predict the black hole count as a function of galaxy stellar mass. We show, for example, that a galaxy like the Milky Way should host millions of $\sim 30~M_\odot$ black holes and dwarf satellite galaxies like Draco should host $\sim 100$ such remnants, with weak dependence on the assumed IMF and stellar evolution model. Most low-mass black holes ($\sim10 M_\odot$) typically reside within massive galaxies ($M_\star \simeq 10^{11} M_\odot$) while massive black holes ($\sim 50~M_\odot$) typically reside within dwarf galaxies ($M_\odot \simeq 10^9 M_\odot$) today. If roughly $1\%$ of black holes are involved in a binary black hole merger, then the reported merger rate densities from Advanced LIGO can be accommodated for a range of merger timescales, and the detection of mergers with $> 50~M_\odot$ black holes should be expected within the next decade. Identifying the host galaxy population of the mergers provides a way to constrain both the binary neutron star or black hole formation efficiencies and the merger timescale distributions; these events would be primarily localized in dwarf galaxies if the merger timescale is short compared to the age of the universe and in massive galaxies otherwise. As more mergers are detected, the prospect of identifying the host galaxy population, either directly through the detection of electromagnetic counterparts of binary neutron star mergers or indirectly through the anisotropy of the events, will become a realistic possibility.

We present an analysis of the radial angular profile of the galacto-isotropic (GI) $\gamma$-ray flux--the statistically uniform flux in circular annuli about the Galactic center. Two different approaches are used to measure the GI flux profile in 85 months of Fermi-LAT data: the BDS statistic method which identifies spatial correlations, and a new Poisson ordered-pixel method which identifies non-Poisson contributions. Both methods produce similar GI flux profiles. The GI flux profile is well-described by an existing model of bremsstrahlung, $\pi^0$ production, inverse Compton scattering, and the isotropic background. Discrepancies with data in our full-sky model are not present in the GI component, and are therefore due to mis-modeling of the non-GI emission. Dark matter annihilation constraints based solely on the observed GI profile are close to the thermal WIMP cross section below 100 GeV, for fixed models of the dark matter density profile and astrophysical $\gamma$-ray foregrounds. Refined measurements of the GI profile are expected to improve these constraints by a factor of a few.

A few dark matter substructures have recently been detected in strong gravitational lenses though their perturbations of highly magnified images. We derive a characteristic scale for lensing perturbations and show that this is significantly larger than the perturber's Einstein radius. We show that the perturber's projected mass enclosed within this radius, scaled by the log-slope of the host galaxy's density profile, can be robustly inferred even if the inferred density profile and tidal radius of the perturber are biased. We demonstrate the validity of our analytic derivation by using several gravitational lens simulations where the tidal radii and the inner log-slopes of the density profile of the perturbing subhalo are allowed to vary. By modeling these simulated data we find that our mass estimator, which we call the effective subhalo lensing mass, is accurate to within about 10\% or smaller in each case, whereas the inferred total subhalo mass can potentially be biased by nearly an order of magnitude. We therefore recommend that the effective subhalo lensing mass be reported in future lensing reconstructions, as this will allow for a more accurate comparison with the results of dark matter simulations.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

We study the atomic physics and the astrophysical implications of a model in which the dark matter is the analog of hydrogen in a secluded sector. The self interactions between dark matter particles include both elastic scatterings as well as inelastic processes due to a hyperfine transition. The self-interaction cross sections are computed by numerically solving the coupled Schr\"{o}dinger equations for this system. We show that these self interactions exhibit the right velocity dependence to explain the low dark matter density cores seen in small galaxies while being consistent with all constraints from observations of clusters of galaxies. For a viable solution, the dark hydrogen mass has to be in 10--100 GeV range and the dark fine-structure constant has to be larger than 0.02. Precisely for this range of parameters, we show that significant cooling losses may occur due to inelastic excitations to the hyperfine state and subsequent decays, with implications for the evolution of low-mass halos and the early growth of supermassive black holes. Cooling from excitations to higher $n$ levels of dark hydrogen and subsequent decays is possible at the cluster scale, with a strong dependence on halo mass. Finally, we show that the minimum halo mass is in the range of $10^{3.5}$ to $10^7 M_\odot$ for the viable regions of parameter space, significantly larger than the typical predictions for weakly-interacting dark matter models. This pattern of observables in cosmological structure formation is unique to this model, making it possible to rule in or rule out hidden sector hydrogen as a viable dark matter model.

We perform a composite likelihood analysis of subdivided regions within the central $26^\circ\times20^\circ$ of the Milky Way, with the aim of characterizing the spectrum of the gamma-ray galactic center excess in regions of varying galactocentric distance. Outside of the innermost few degrees, we find that the radial profile of the excess is background-model dependent and poorly constrained. The spectrum of the excess emission is observed to extend upwards of 10 GeV outside $\sim5^\circ$ in radius, but cuts off steeply between 10--20 GeV only in the innermost few degrees. If interpreted as a real feature of the excess, this radial variation in the spectrum has important implications for both astrophysical and dark matter interpretations of the galactic center excess. Single-component dark matter annihilation models face challenges in reproducing this variation; on the other hand, a population of unresolved millisecond pulsars contributing both prompt and secondary inverse Compton emission may be able to explain the spectrum as well as its spatial dependency. We show that the expected differences in the photon-count distributions of a smooth dark matter annihilation signal and an unresolved point source population are an order of magnitude smaller than the fluctuations in residuals after fitting the data, which implies that mismodeling is an important systematic effect in point source analyses aimed at resolving the gamma-ray excess.

The rotation curves of spiral galaxies exhibit a diversity that has been difficult to understand in the cold dark matter (CDM) paradigm. We show that the self-interacting dark matter (SIDM) model provides excellent fits to the rotation curves of a sample of galaxies with asymptotic velocities in the 25 to 300 km/s range that exemplify the full range of diversity. We only assume the halo concentration-mass relation predicted by the CDM model and a fixed value of the self-interaction cross section.In dark matter dominated galaxies, thermalization due to self-interactions creates large cores and reduces dark matter densities. In contrast, thermalization leads to denser and smaller cores in more luminous galaxies, and naturally explains the flat rotation curves of the highly luminous galaxies. Our results demonstrate that the impact of the baryons on the SIDM halo profile and the scatter from the assembly history of halos as encoded in the concentration-mass relation can explain the diverse rotation curves of spiral galaxies.

We investigate the response of self-interacting dark matter (SIDM) halos to the growth of galaxy potentials using idealized simulations, each run in tandem with standard collisionless Cold Dark Matter (CDM). We find a greater diversity in the SIDM halo profiles compared to the CDM halo profiles. If the stellar gravitational potential strongly dominates in the central parts of a galaxy, then SIDM halos can be as dense as CDM halos on observable scales. For extreme cases with highly compact disks core collapse can occur, leading to SIDM halos that are denser and cuspier than their CDM counterparts. If the stellar potential is not dominant, then SIDM halos retain constant density cores with densities far below CDM predictions. When a disk potential is present, the inner SIDM halo becomes \em{more flattened} in the disk plane than the CDM halo. These results are in excellent quantitative agreement with the predictions of Kaplinghat et al. (2014). We also simulated a galaxy cluster halo with a central stellar distribution similar to the brightest central galaxy of the cluster A2667. A SIDM halo simulated with cross section over mass $\sigma/m = 0.1\ \mathrm{cm^2 g^{-1}}$ provides a good match to the measured dark matter density profile of A2667, while an adiabatically-contracted CDM halo is denser and cuspier. The cored profile of the same halo simulated with $\sigma/m = 0.5\ \mathrm{cm^2 g^{-1}}$ is not dense enough to match A2667. Our findings are in agreement with previous results that $\sigma/m \gtrsim 0.1\ \mathrm{cm^2 g^{-1}}$ is disfavored for dark matter collision velocities in excess of about 1500 km/s. More generally, the predictive cross-talk between baryonic potentials and SIDM density distributions offers new directions for constraining SIDM cross sections in massive galaxies where baryons are dynamically important.

Astrophysical observations spanning dwarf galaxies to galaxy clusters indicate that dark matter (DM) halos are less dense in their central regions compared to expectations from collisionless DM N-body simulations. Using detailed fits to DM halos of galaxies and clusters, we show that self-interacting DM (SIDM) may provide a consistent solution to the DM deficit problem across all scales, even though individual systems exhibit a wide diversity in halo properties. Since the characteristic velocity of DM particles varies across these systems, we are able to measure the self-interaction cross section as a function of kinetic energy and thereby deduce the SIDM particle physics model parameters. Our results prefer a mildly velocity-dependent cross section, from $\sigma/m \simeq 2\; {\rm cm^2/g}$ on galaxy scales to $\sigma/m \simeq 0.1\; {\rm cm^2/g}$ on cluster scales, consistent with the upper limits from merging clusters. Our results dramatically improve the constraints on SIDM models and may allow the masses of both DM and dark mediator particles to be measured even if the dark sector is completely hidden from the Standard Model, which we illustrate for the dark photon model.

We present a new extended gamma ray excess detected with the Fermi Satellite Large Area Telescope toward the Galactic Center that traces the morphology of infrared starlight emission. Combined with its measured spectrum, this new extended source is approximately consistent with inverse Compton emission from a high-energy electron-positron population with energies up to about 10 GeV. Previously detected emissions tracing the 20 cm radio, interpreted as bremsstrahlung radiation, and the Galactic Center Extended emission tracing a spherical distribution and peaking at 2 GeV, are also detected. We show that the inverse Compton and bremsstrahlung emissions are likely due to the same source of electrons and positrons. All three extended emissions may be explained within the framework of a model where the dark matter annihilates to leptons or a model with unresolved millisecond pulsars in the Galactic Center.

We consider a simple supersymmetric hidden sector: pure SU(N) gauge theory. Dark matter is made up of hidden glueballinos with mass $m_X$ and hidden glueballs with mass near the confinement scale $\Lambda$. For $m_X \sim 1\,\text{TeV}$ and $\Lambda \sim 100\,\text{MeV}$, the glueballinos freeze out with the correct relic density and self-interact through glueball exchange to resolve small-scale structure puzzles. An immediate consequence is that the glueballino spectrum has a hyperfine splitting of order $\Lambda^2 / m_X \sim 10\,\text{keV}$. We show that the radiative decays of the excited state can explain the observed 3.5 keV X-ray line signal from clusters of galaxies, Andromeda, and the Milky Way.

There is strong evidence in favor of the idea that dark matter is self interacting, with the cross section-to-mass ratio $\sigma / m \sim 1\,\mathrm{cm^2/g} \sim 1\,\mathrm{barn/GeV}$. We show that viable models of dark matter with this large cross section are straightforwardly realized with non-Abelian hidden sectors. In the simplest of such models, the hidden sector is a pure gauge theory, and the dark matter is composed of hidden glueballs with a mass around $100\,\mathrm{MeV}$. Alternatively, the hidden sector may be a supersymmetric pure gauge theory with a $\sim 10\,\mathrm{TeV}$ gluino thermal relic. In this case, the dark matter is largely composed of glueballinos that strongly self interact through the exchange of light glueballs. We present a unified framework that realizes both of these possibilities in anomaly-mediated supersymmetry breaking, where, depending on a few model parameters, the dark matter may be composed of hidden glueballinos, hidden glueballs, or a mixture of the two. These models provide simple examples of multicomponent dark matter, have interesting implications for particle physics and cosmology, and include cases where a subdominant component of dark matter may be extremely strongly self interacting, with interesting astrophysical consequences.

In this report we summarize the many dark matter searches currently being pursued through four complementary approaches: direct detection, indirect detection, collider experiments, and astrophysical probes. The essential features of broad classes of experiments are described, each with their own strengths and weaknesses. The complementarity of the different dark matter searches is discussed qualitatively and illustrated quantitatively in two simple theoretical frameworks. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program drawing from all four approaches.

We point out three correlated predictions of the axion monodromy inflation model: large amplitude of gravitational waves, suppression of power on horizon scales and on scales relevant for the formation of dwarf galaxies. While these predictions are likely generic to models with oscillations in the inflaton potential, the axion monodromy model naturally accommodates the required running spectral index through Planck-scale corrections to the inflaton potential. Applying this model to a combined data set of Planck, ACT, SPT, and WMAP low-$\ell$ polarization cosmic microwave background (CMB) data, we find a best-fit tensor-to-scalar ratio $r_{0.05} = 0.07^{+0.05}_{-0.04}$ due to gravitational waves, which may have been observed by the BICEP2 experiment. Despite the contribution of gravitational waves, the total power on large scales (CMB power spectrum at low multipoles) is lower than the standard $\Lambda$CDM cosmology with a power-law spectrum of initial perturbations and no gravitational waves, thus mitigating some of the tension on large scales. There is also a reduction in the matter power spectrum of 20-30\% at scales corresponding to $k = 10~{\rm Mpc}^{-1}$, which are relevant for dwarf galaxy formation. This will alleviate some of the unsolved small-scale structure problems in the standard $\Lambda$CDM cosmology. The inferred matter power spectrum is also found to be consistent with recent Lyman-$\alpha$ forest data, which is in tension with the Planck-favored $\Lambda$CDM model with power-law primordial power spectrum.

Observations by the Fermi-LAT telescope have uncovered a significant $\gamma$-ray excess toward the Milky Way Galactic Center. There has been no detection of a similar signal in the direction of the Milky Way dwarf spheroidal galaxies. Additionally, astronomical observations indicate that dwarf galaxies and other faint galaxies are less dense than predicted by the simplest cold dark matter models. We show that a self-interacting dark matter model with a particle mass of roughly 50 GeV annihilating to the mediator responsible for the strong self-interaction can simultaneously explain all three observations. The mediator is necessarily unstable and its mass must be below about 100 MeV in order to lower densities in faint galaxies. If the mediator decays to electron-positron pairs with a cross section on the order of the thermal relic value, then we find that these pairs can up-scatter the interstellar radiation field and produce the observed $\gamma$-ray excess. We show that this model is compatible with all current constraints and highlight detectable signatures unique to self-interacting dark matter models.

The BICEP2 results, when interpreted as a gravitational wave signal and combined with other CMB data, suggest a roll-off in power towards small scales in the primordial matter power spectrum. Among the simplest possibilities is a running of the spectral index. Here we show that the preferred level of running alleviates small-scale issues within the $\Lambda$CDM model, more so even than viable WDM models. We use cosmological zoom-in simulations of a Milky Way-size halo along with full-box simulations to compare predictions among four separate cosmologies: a BICEP2-inspired running index model ($\alpha_s$ = -0.024), two fixed-tilt $\Lambda$CDM models motivated by Planck, and a 2.6 keV thermal WDM model. We find that the running BICEP2 model reduces the central densities of large dwarf-size halos ($V_\mathrm{max}$ ~ 30 - 80 km s$^{-1}$) and alleviates the too-big-to-fail problem significantly compared to our adopted Planck and WDM cases. Further, the BICEP2 model suppresses the count of small subhalos by ~50% relative to Planck models, and yields a significantly lower "boost" factor for dark matter annihilation signals. Our findings highlight the need to understand the shape of the primordial power spectrum in order to correctly interpret small-scale data.

We construct empirical models of the diffuse gamma-ray background toward the Galactic Center. Including all known point sources and a template of emission associated with interactions of cosmic rays with molecular gas, we show that the extended emission observed previously in the Fermi Large Area Telescope data toward the Galactic Center is detected at high significance for all permutations of the diffuse model components. However, we find that the fluxes and spectra of the sources in our model change significantly depending on the background model. In particular, the spectrum of the central Sgr A$^\ast$ source is less steep than in previous works and the recovered spectrum of the extended emission has large systematic uncertainties, especially at lower energies. If the extended emission is interpreted to be due to dark matter annihilation, we find annihilation into pure $b$-quark and $\tau$-lepton channels to be statistically equivalent goodness of fits. In the case of the pure $b$-quark channel, we find a dark matter mass of $39.4\left(^{+3.7}_{-2.9}\rm\ stat.\right)\left(\pm 7.9\rm\ sys.\right)\rm\ GeV$, while a pure $\tau^{+} \tau^{-}$-channel case has an estimated dark matter mass of $9.43\left(^{+0.63}_{-0.52}\rm\ stat.\right)(\pm 1.2\rm\ sys.)\ GeV$. Alternatively, if the extended emission is interpreted to be astrophysical in origin such as due to unresolved millisecond pulsars, we obtain strong bounds on dark matter annihilation, although systematic uncertainties due to the dependence on the background models are significant.

Self-interacting dark matter (SIDM) models have been proposed to solve the small-scale issues with the collisionless cold dark matter (CDM) paradigm. We derive equilibrium solutions in these SIDM models for the dark matter halo density profile including the gravitational potential of both baryons and dark matter. Self-interactions drive dark matter to be isothermal and this ties the core sizes and shapes of dark matter halos to the spatial distribution of the stars, a radical departure from previous expectations and from CDM predictions. Compared to predictions of SIDM-only simulations, the core sizes are smaller and the core densities are higher, with the largest effects in baryon-dominated galaxies. As an example, we find a core size around 0.5 kpc for dark matter in the Milky Way, more than an order of magnitude smaller than the core size from SIDM-only simulations, which has important implications for indirect searches of SIDM candidates.

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 4, on the Cosmic Frontier, discusses the program of research relevant to cosmology and the early universe. This area includes the study of dark matter and the search for its particle nature, the study of dark energy and inflation, and cosmic probes of fundamental symmetries.

We measured velocities of 862 individual red giant stars in seven isolated dwarf galaxies in the Local Group: NGC 6822, IC 1613, VV 124 (UGC 4879), the Pegasus dwarf irregular galaxy (DDO 216), Leo A, Cetus, and Aquarius (DDO 210). We also computed velocity dispersions, taking into account the measurement uncertainties on individual stars. None of the isolated galaxies is denser than the densest Local Group satellite galaxy. Furthermore, the isolated dwarf galaxies have no obvious distinction in the velocity dispersion--half-light radius plane from the satellite galaxies of the Milky Way and M31. The similarity of the isolated and satellite galaxies' dynamics and structural parameters imposes limitations on environmental solutions to the too-big-to-fail problem, wherein there are fewer dense dwarf satellite galaxies than would be expected from cold dark matter simulations. This data set also has many other applications for dwarf galaxy evolution, including the transformation of dwarf irregular into dwarf spheroidal galaxies. We intend to explore these issues in future work.

We present 3D kinematic observations of stars within the central 0.5 pc of the Milky Way nuclear star cluster using adaptive optics imaging and spectroscopy from the Keck telescopes. Recent observations have shown that the cluster has a shallower surface density profile than expected for a dynamically relaxed cusp, leading to important implications for its formation and evolution. However, the true three dimensional profile of the cluster is unknown due to the difficulty in de-projecting the stellar number counts. Here, we use spherical Jeans modeling of individual proper motions and radial velocities to constrain for the first time, the de-projected spatial density profile, cluster velocity anisotropy, black hole mass ($M_\mathrm{BH}$), and distance to the Galactic center ($R_0$) simultaneously. We find that the inner stellar density profile of the late-type stars, $\rho(r)\propto r^{-\gamma}$ to have a power law slope $\gamma=0.05_{-0.60}^{+0.29}$, much more shallow than the frequently assumed Bahcall $\&$ Wolf slope of $\gamma=7/4$. The measured slope will significantly affect dynamical predictions involving the cluster, such as the dynamical friction time scale. The cluster core must be larger than 0.5 pc, which disfavors some scenarios for its origin. Our measurement of $M_\mathrm{BH}=5.76_{-1.26}^{+1.76}\times10^6$ $M_\odot$ and $R_0=8.92_{-0.55}^{+0.58}$ kpc is consistent with that derived from stellar orbits within 1$^{\prime\prime}$ of Sgr A*. When combined with the orbit of S0-2, the uncertainty on $R_0$ is reduced by 30% ($8.46_{-0.38}^{+0.42}$ kpc). We suggest that the MW NSC can be used in the future in combination with stellar orbits to significantly improve constraints on $R_0$.

We show that the canonical oscillation-based (non-resonant) production of sterile neutrino dark matter is inconsistent at $>99$% confidence with observations of galaxies in the Local Group. We set lower limits on the non-resonant sterile neutrino mass of $2.5$ keV (equivalent to $0.7$ keV thermal mass) using phase-space densities derived for dwarf satellite galaxies of the Milky Way, as well as limits of $8.8$ keV (equivalent to $1.8$ keV thermal mass) based on subhalo counts of $N$-body simulations of M 31 analogues. Combined with improved upper mass limits derived from significantly deeper X-ray data of M 31 with full consideration for background variations, we show that there remains little room for non-resonant production if sterile neutrinos are to explain $100$% of the dark matter abundance. Resonant and non-oscillation sterile neutrino production remain viable mechanisms for generating sufficient dark matter sterile neutrinos.

In this Report we discuss the four complementary searches for the identity of dark matter: direct detection experiments that look for dark matter interacting in the lab, indirect detection experiments that connect lab signals to dark matter in our own and other galaxies, collider experiments that elucidate the particle properties of dark matter, and astrophysical probes sensitive to non-gravitational interactions of dark matter. The complementarity among the different dark matter searches is discussed qualitatively and illustrated quantitatively in several theoretical scenarios. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program based on all four of those approaches.

Dark matter self-interactions have important implications for the distributions of dark matter in the Universe, from dwarf galaxies to galaxy clusters. We present benchmark models that illustrate characteristic features of dark matter that is self-interacting through a new light mediator. These models have self-interactions large enough to change dark matter densities in the centers of galaxies in accord with observations, while remaining compatible with large-scale structure data and all astrophysical observations such as halo shapes and the Bullet Cluster. These observations favor a mediator mass in the 10 - 100 MeV range and large regions of this parameter space are accessible to direct detection experiments like LUX, SuperCDMS, and XENON1T.

Super-weakly interacting massive particles produced in the late decays of weakly interacting massive particles (WIMPs) are generic in large regions of supersymmetric parameter space and other frameworks for physics beyond the standard model. If their masses are similar to that of the decaying WIMP, then they could naturally account for all of the cosmological dark matter abundance. Their astrophysical consequences and collider signatures are distinct and different from WIMP candidates. In particular, they could modify Big Bang Nucleosynthesis, distort the Cosmic Microwave Background, reduce galactic substructure and lower central densities of low-mass galaxies.

We show the existence of a statistically significant, robust detection of a gamma-ray source in the Milky Way Galactic Center that is consistent with a spatially extended signal using about 4 years of Fermi-LAT data. The gamma-ray flux is consistent with annihilation of dark matter particles with a thermal annihilation cross-section if the spatial distribution of dark matter particles is similar to the predictions of dark matter only simulations. We find statistically significant detections of an extended source with gamma-ray spectrum that is consistent with dark matter particle masses of approximately 10 GeV to 1 TeV annihilating to b/b-bar quarks, and masses approximately 10 GeV to 30 GeV annihilating to tau+ tau- leptons. However, a part of the allowed region in this interpretation is in conflict with constraints from Fermi observations of the Milky Way satellites. The biggest improvement over the fit including just the point sources is obtained for a 30 GeV dark matter particle annihilating to b/b-bar quarks. The gamma-ray intensity and spectrum are also well fit with emission from a millisecond pulsar (MSP) population following a density profile like that of low-mass X-ray binaries observed in M31. The greatest goodness-of-fit of the extended emission is with spectra consistent with known astrophysical sources like MSPs in globular clusters or cosmic ray bremsstrahlung on molecular gas. Therefore, we conclude that the bulk of the emission is likely from an unresolved or spatially extended astrophysical source. However, the interesting possibility of all or part of the extended emission being from dark matter annihilation cannot be excluded at present.

Segue 2, discovered by Belokurov et al. (2009), is a galaxy with a luminosity of only 900 L_sun. We present Keck/DEIMOS spectroscopy of 25 members of Segue 2--a threefold increase in spectroscopic sample size. The velocity dispersion is too small to be measured with our data. The upper limit with 90% (95%) confidence is sigma_v < 2.2 (2.6) km/s, the most stringent limit for any galaxy. The corresponding limit on the mass within the 3-D half-light radius (46 pc) is M_1/2 < 1.5 (2.1) x 10^5 M_sun. Segue 2 is the least massive galaxy known. We identify Segue 2 as a galaxy rather than a star cluster based the wide dispersion in [Fe/H] (from -2.85 to -1.33) among the member stars. The stars' [alpha/Fe] ratios decline with increasing [Fe/H], indicating that Segue 2 retained Type Ia supernova ejecta despite its presently small mass and that star formation lasted for at least 100 Myr. The mean metallicity, <[Fe/H]> = -2.22 +/- 0.13 (about the same as the Ursa Minor galaxy, 330 times more luminous than Segue 2), is higher than expected from the luminosity-metallicity relation defined by more luminous dwarf galaxy satellites of the Milky Way. Segue 2 may be the barest remnant of a tidally stripped, Ursa Minor-sized galaxy. If so, it is the best example of an ultra-faint dwarf galaxy that came to be ultra-faint through tidal stripping. Alternatively, Segue 2 could have been born in a very low-mass dark matter subhalo (v_max < 10 km/s), below the atomic hydrogen cooling limit.

We find that the viability of a cosmological model that incorporates 2 sterile neutrinos with masses around 1 eV each, as favored by global neutrino oscillation analyses including short baseline results, is significantly dependent on the choice of datasets included in the analysis and the ability to control the systematic uncertainties associated with these datasets. Our analysis includes a variety of cosmological probes including the cosmic microwave background (WMAP7+SPT), Hubble constant (HST), galaxy power spectrum (SDSS-DR7), and supernova distances (SDSS and Union2 compilations). In the joint observational analysis, our sterile neutrino model is equally favored as a LCDM model when using the MLCS light curve fitter for the supernova measurements, and strongly disfavored by the data at \Delta\chi^2 ~ 18 when using the SALT2 fitter. When excluding the supernova measurements, the sterile neutrino model is disfavored by the other datasets at \Delta\chi^2 ~ 12, and at best becomes mildly disfavored at \Delta\chi^2 ~ 3 when allowing for curvature, evolving dark energy, additional relativistic species, running of the spectral index, and freedom in the primordial helium abundance. No single additional parameter accounts for most of this effect. Therefore, if laboratory experiments continue to favor a scenario with roughly eV mass sterile neutrinos, and if this becomes decisively disfavored by cosmology, then a more exotic cosmological model than explored here may become necessary.

We present a method for identifying localized secondary populations in stellar velocity data using Bayesian statistical techniques. We apply this method to the dwarf spheroidal galaxy Ursa Minor and find two secondary objects in this satellite of the Milky Way. One object is kinematically cold with a velocity dispersion of $4.25 \pm 0.75\ \kms$ and centered at $(9.1\arcmin \pm 1.5, 7.2\arcmin \pm 1.2)$ in relative RA and DEC with respect to the center of Ursa Minor. The second object has a large velocity offset of $-12.8^{+1.75}_{-1.5}\ \kms$ compared to Ursa Minor and centered at $(-14.0\arcmin^{+2.4}_{-5.8}, -2.5\arcmin^{+0.4}_{-1.0})$. The kinematically cold object has been found before using a smaller data set but the prediction that this cold object has a velocity dispersion larger than $2.0\ \kms$ at 95% C.L. differs from previous work. We use two and three component models along with the information criteria and Bayesian evidence model selection methods to argue that Ursa Minor has one or two localized secondary populations. The significant probability for a large velocity dispersion in each secondary object raises the intriguing possibility that each has its own dark matter halo, that is, it is a satellite of a satellite of the Milky Way.

We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass \sigma/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consistently from the Boltzmann equation and that reproduces analytic expectations in controlled numerical experiments. We find that well-resolved SIDM halos have constant-density cores, with significantly lower central densities than their CDM counterparts. In contrast, the subhalo content of SIDM halos is only modestly reduced compared to CDM, with the suppression greatest for large hosts and small halo-centric distances. Moreover, the large-scale clustering and halo circular velocity functions in SIDM are effectively identical to CDM, meaning that all of the large-scale successes of CDM are equally well matched by SIDM. From our largest cross section runs we are able to extract scaling relations for core sizes and central densities over a range of halo sizes and find a strong correlation between the core radius of an SIDM halo and the NFW scale radius of its CDM counterpart. We construct a simple analytic model, based on CDM scaling relations, that captures all aspects of the scaling relations for SIDM halos. Our results show that halo core densities in \sigma/m = 1 cm^2/g models are too low to match observations of galaxy clusters, low surface brightness spirals (LSBs), and dwarf spheroidal galaxies. However, SIDM with \sigma/m ~ 0.1 cm^2/g appears capable of reproducing reported core sizes and central densities of dwarfs, LSBs, and galaxy clusters without the need for velocity dependence. (abridged)

If dark matter has a large self-interaction scattering cross section, then interactions among dark-matter particles will drive galaxy and cluster halos to become spherical in their centers. Work in the past has used this effect to rule out velocity-independent, elastic cross sections larger than sigma/m ~ 0.02 cm^2/g based on comparisons to the shapes of galaxy cluster lensing potentials and X-ray isophotes. In this paper, we use cosmological simulations to show that these constraints were off by more than an order of magnitude because (a) they did not properly account for the fact that the observed ellipticity gets contributions from the triaxial mass distribution outside the core set by scatterings, (b) the scatter in axis ratios is large and (c) the core region retains more of its triaxial nature than estimated before. Including these effects properly shows that the same observations now allow dark matter self-interaction cross sections at least as large as sigma/m = 0.1 cm^2/g. We show that constraints on self-interacting dark matter from strong-lensing clusters are likely to improve significantly in the near future, but possibly more via central densities and core sizes than halo shapes.

We forecast the expected cosmological constraints from a combination of probes of both the universal expansion rate and matter perturbation growth, in the form of weak lensing tomography, galaxy tomography, supernovae, and the cosmic microwave background incorporating all cross-correlations between the observables for an extensive cosmological parameter set. We allow for non-zero curvature and parameterize our ignorance of the early universe by allowing for a non-negligible fraction of dark energy (DE) at high redshifts. We find that early DE density can be constrained to 0.2% of the critical density of the universe with Planck combined with a ground-based LSST-like survey, while curvature can be constrained to 0.06%. However, these additional degrees of freedom degrade our ability to measure late-time dark energy and the sum of neutrino masses. We find that the combination of cosmological probes can break degeneracies and constrain the sum of neutrino masses to 0.04 eV, present DE density also to 0.2% of the critical density, and the equation of state to 0.01 - roughly a factor of two degradation in the constraints overall compared to the case without allowing for early DE. The constraints for a space-based mission are similar. Even a modest 1% dark energy fraction of the critical density at high redshift, if not accounted for in future analyses, biases the cosmological parameters by up to 2 sigma. Our analysis suggests that throwing out nonlinear scales (multipoles > 1000) may not result in significant degradation in future parameter measurements when multiple cosmological probes are combined. We find that including cross-correlations between the different probes can result in improved constraints by up to a factor of 2 for the sum of neutrino masses and early dark energy density.

We use the Aquarius simulations to show that the most massive subhalos in galaxy-mass dark matter halos in LCDM are grossly inconsistent with the dynamics of the brightest Milky Way dwarf spheroidal galaxies. While the best-fitting hosts of the dwarf spheroidals all have 12 < Vmax < 25 km/s, LCDM simulations predict at least ten subhalos with Vmax > 25 km/s. These subhalos are also among the most massive at earlier times, and significantly exceed the UV suppression mass back to z ~ 10. No LCDM-based model of the satellite population of the Milky Way explains this result. The problem lies in the satellites' densities: it is straightforward to match the observed Milky Way luminosity function, but doing so requires the dwarf spheroidals to have dark matter halos that are a factor of ~5 more massive than is observed. Independent of the difficulty in explaining the absence of these dense, massive subhalos, there is a basic tension between the derived properties of the bright Milky Way dwarf spheroidals and LCDM expectations. The inferred infall masses of these galaxies are all approximately equal and are much lower than standard LCDM predictions for systems with their luminosities. Consequently, their implied star formation efficiencies span over two orders of magnitude, from 0.2% to 20% of baryons converted into stars, in stark contrast with expectations gleaned from more massive galaxies. We explore possible solutions to these problems within the context of LCDM and find them to be unconvincing. In particular, we use controlled simulations to demonstrate that the small stellar masses of the bright dwarf spheroidals make supernova feedback an unlikely explanation for their low inferred densities.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

Dwarf spheroidal galaxies have a large mass to light ratio and low astrophysical background, and are therefore considered one of the most promising targets for dark matter searches in the gamma-ray band. By applying a joint likelihood analysis, the power of resultant limits in case of no detection can be enhanced and robust constraints on the dark matter parameter space can be obtained. We present results from a combined analysis of 10 dwarf spheroidal galaxies using Fermi-LAT data. Different annihilation channels have been analyzed and uncertainties from astrophysical properties have been taken into account.

Primordial non-Gaussianity is a crucial test of inflationary cosmology. We consider the impact of non-Gaussianity on the ionization power spectrum from 21 cm emission during the epoch of reionization. We focus on the power spectrum on large scales at redshifts of 7 to 8 and explore the expected constraint on the local non-Gaussianity parameter f_NL for current and next-generation 21 cm experiments. We show that experiments such as SKA and MWA could measure f_NL values of order 10. This can be improved by an order of magnitude with a fast-Fourier transform telescope like Omniscope.

We introduce a comprehensive analysis of multi-epoch stellar line-of-sight velocities to determine the intrinsic velocity dispersion of the ultrafaint satellites of the Milky Way. Our method includes a simultaneous Bayesian analysis of both membership probabilities and the contribution of binary orbital motion to the observed velocity dispersion within a 14-parameter likelihood. We apply our method to the Segue 1 dwarf galaxy and conclude that Segue 1 is a dark-matter-dominated galaxy at high probability with an intrinsic velocity dispersion of 3.7^{+1.4}_{-1.1} km/sec. The dark matter halo required to produce this dispersion must have an average density of 2.5^{+4.1}_{-1.9} solar mass/pc^3 within a sphere that encloses half the galaxy's stellar luminosity. This is the highest measured density of dark matter in the Local Group. Our results show that a significant fraction of the stars in Segue 1 may be binaries with the most probable mean period close to 10 years, but also consistent with the 180 year mean period seen in the solar vicinity at about 1 sigma. Despite this binary population, the possibility that Segue 1 is a bound star cluster with the observed velocity dispersion arising from the orbital motion of binary stars is disfavored by the multi-epoch stellar velocity data at greater than 99% C.L. Finally, our treatment yields a projected (two-dimensional) half-light radius for the stellar profile of Segue 1 of 28^{+5}_{-4} pc, in excellent agreement with photometric measurements.

We show that dissipationless LCDM simulations predict that the majority of the most massive subhaloes of the Milky Way are too dense to host any of its bright satellites (L_V > 10^5 L_sun). These dark subhaloes have circular velocities at infall of 30-70 km/s and infall masses of [0.2-4] x 10^10 M_sun. Unless the Milky Way is a statistical anomaly, this implies that galaxy formation becomes effectively stochastic at these masses. This is in marked contrast to the well-established monotonic relation between galaxy luminosity and halo circular velocity (or halo mass) for more massive haloes. We show that at least two (and typically four) of these massive dark subhaloes are expected to produce a larger dark matter annihilation flux than Draco. It may be possible to circumvent these conclusions if baryonic feedback in dwarf satellites or different dark matter physics can reduce the central densities of massive subhaloes by order unity on a scale of 0.3 - 1 kpc.

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

The annihilation cross section of thermal relic dark matter determines both its relic density and indirect detection signals. We determine how large indirect signals may be in scenarios with Sommerfeld-enhanced annihilation, subject to the constraint that the dark matter has the correct relic density. This work refines our previous analysis through detailed treatments of resonant Sommerfeld enhancement and the effect of Sommerfeld enhancement on freeze out. Sommerfeld enhancements raise many interesting issues in the freeze out calculation, and we find that the cutoff of resonant enhancement, the equilibration of force carriers, the temperature of kinetic decoupling, and the efficiency of self-interactions for preserving thermal velocity distributions all play a role. These effects may have striking consequences; for example, for resonantly-enhanced Sommerfeld annihilation, dark matter freezes out but may then chemically recouple, implying highly suppressed indirect signals, in contrast to naive expectations. In the minimal scenario with standard astrophysical assumptions, and tuning all parameters to maximize the signal, we find that, for force-carrier mass m_phi = 250 MeV and dark matter masses m_X = 0.1, 0.3, and 1 TeV, the maximal Sommerfeld enhancement factors are S_eff = 7, 30, and 90, respectively. Such boosts are too small to explain both the PAMELA and Fermi excesses. Non-minimal models may require smaller boosts, but the bounds on S_eff could also be more stringent, and dedicated freeze out analyses are required. For concreteness, we focus on 4 mu final states, but we also discuss 4 e and other modes, deviations from standard astrophysical assumptions and non-minimal particle physics models, and we outline the steps required to determine if such considerations may lead to a self-consistent explanation of the PAMELA or Fermi excesses.

We derive an accurate mass estimator for dispersion-supported stellar systems and demonstrate its validity by analyzing resolved line-of-sight velocity data for globular clusters, dwarf galaxies, and elliptical galaxies. Specifically, by manipulating the spherical Jeans equation we show that the dynamical mass enclosed within the 3D deprojected half-light radius r_1/2 can be determined with only mild assumptions about the spatial variation of the stellar velocity dispersion anisotropy. We find M_1/2 = 3 \sigma_los^2 r_1/2 / G ~ 4 \sigma_los^2 R_eff / G, where \sigma_los^2 is the luminosity-weighted square of the line-of-sight velocity dispersion and R_eff is the 2D projected half-light radius. While deceptively familiar in form, this formula is not the virial theorem, which cannot be used to determine accurate masses unless the radial profile of the total mass is known a priori. We utilize this finding to show that all of the Milky Way dwarf spheroidal galaxies (MW dSphs) are consistent with having formed within a halo of mass approximately 3 x 10^9 M_sun in Lambda CDM cosmology. The faintest MW dSphs seem to have formed in dark matter halos that are at least as massive as those of the brightest MW dSphs, despite the almost five orders of magnitude spread in luminosity. We expand our analysis to the full range of observed dispersion-supported stellar systems and examine their I-band mass-to-light ratios (M/L). The M/L vs. M_1/2 relation for dispersion-supported galaxies follows a U-shape, with a broad minimum near M/L ~ 3 that spans dwarf elliptical galaxies to normal ellipticals, a steep rise to M/L ~ 3,200 for ultra-faint dSphs, and a more shallow rise to M/L ~ 800 for galaxy cluster spheroids.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We predict that there is a population of low-luminosity dwarf galaxies orbiting within the halo of the Milky Way that have surface brightnesses low enough to have escaped detection in star-count surveys. The overall count of stealth galaxies is sensitive to the presence (or lack) of a low-mass threshold in galaxy formation. These systems have luminosities and stellar velocity dispersions that are similar to those of known ultrafaint dwarf galaxies but they have more extended stellar distributions (half light radii greater than about 100 pc) because they inhabit dark subhalos that are slightly less massive than their higher surface brightness counterparts. As a result, the typical peak surface brightness is fainter than 30 mag per square arcsec. One implication is that the inferred common mass scale for Milky Way dwarfs may be an artifact of selection bias. If there is no sharp threshold in galaxy formation at low halo mass, then ultrafaint galaxies like Segue 1 represent the high-mass, early forming tail of a much larger population of objects that could number in the hundreds and have typical peak circular velocities of about 8 km/s and masses within 300 pc of about 5 million solar masses. Alternatively, if we impose a low-mass threshold in galaxy formation in order to explain the unexpectedly high densities of the ultrafaint dwarfs, then we expect only a handful of stealth galaxies in the halo of the Milky Way. A complete census of these objects will require deeper sky surveys, 30m-class follow-up telescopes, and more refined methods to identify extended, self-bound groupings of stars in the halo.

Dark matter with Sommerfeld-enhanced annihilation has been proposed to explain observed cosmic ray positron excesses in the 10 GeV to TeV energy range. We show that the required enhancement implies thermal relic densities that are too small to be all of dark matter. We also show that the dark matter is sufficiently self-interacting that observations of elliptical galactic dark matter halos exclude large Sommerfeld enhancement for light force carriers. Resonant Sommerfeld enhancement does not modify these conclusions, and the astrophysical boosts required to resolve these discrepancies are disfavored, especially when significant self-interactions suppress halo substructure.

Warm dark matter (WDM) and self-interacting dark matter (SIDM) are often motivated by the inferred cores in the dark matter halos of low surface brightness (LSB) galaxies. We test thermal WDM, non-thermal WDM, and SIDM using high-resolution rotation curves of nine LSB galaxies. We fit these dark matter models to the data and determine the halo core radii and central densities. While the minimum core size in WDM models is predicted to decrease with halo mass, we find that the inferred core radii increase with halo mass and also cannot be explained with a single value of the primordial phase space density. Moreover, if the core size is set by WDM particle properties, then even the smallest cores we infer would require primordial phase space density values that are orders of magnitude smaller than lower limits obtained from the Lyman alpha forest power spectra. We also find that the dark matter halo core densities vary by a factor of about 30 from system to system while showing no systematic trend with the maximum rotation velocity of the galaxy. This strongly argues against the core size being directly set by large self-interactions (scattering or annihilation) of dark matter. We therefore conclude that the inferred cores do not provide motivation to prefer WDM or SIDM over other dark matter models.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

We show that measured velocity dispersions of dwarf spheroidal galaxies from about 4 to 10 km/s are unlikely to be inflated by more than 20% due to the orbital motion of binary stars, and demonstrate that the intrinsic velocity dispersions can be determined to within a few percent accuracy using two-epoch observations with 1-2 years as the optimal time interval. The crucial observable is the threshold fraction--the fraction of stars that show velocity changes larger than a given threshold between measurements. The threshold fraction is tightly correlated with the dispersion introduced by binaries, independent of the underlying binary fraction and distribution of orbital parameters. We outline a simple procedure to correct the velocity dispersion to within a few percent accuracy by using the threshold fraction and provide fitting functions for this method. We also develop a methodology for constraining properties of binary populations from both single- and two-epoch velocity measurements by including the binary velocity distribution in a Bayesian analysis.

The WMAP haze is an excess in the 22 to 93 GHz frequency bands of WMAP extending about 10 degrees from the galactic center. We show that synchrotron emission from electron-positron pairs injected into the interstellar medium by the galactic population of pulsars with energies in the 1 to 100 GeV range can explain the frequency spectrum of the WMAP haze and the drop in the average haze power with latitude. The same spectrum of high energy electron-positron pairs from pulsars, which gives rise to the haze, may also generate the observed excesses in AMS, HEAT and PAMELA. We discuss the spatial morphology of the pulsar synchrotron signal and its deviation from spherical symmetry, which may provide an avenue to determine the pulsar contribution to the haze.

Massive satellite accretions onto early galactic disks can lead to the deposition of dark matter in disk-like configurations that co-rotate with the galaxy. This phenomenon has potentially dramatic consequences for dark matter detection experiments. We utilize focused, high-resolution simulations of accretion events onto disks designed to be Galaxy analogues, and compare the resultant disks to the morphological and kinematic properties of the Milky Way's thick disk in order to bracket the range of co-rotating accreted dark matter. We find that the Milky Way's merger history must have been unusually quiescent compared to median LCDM expectations and therefore its dark disk must be relatively small: the fraction of accreted dark disk material near the Sun is about 20% of the host halo density or smaller and the co-rotating dark matter fraction near the Sun, defined as particles moving with a rotational velocity lag less than 50 km/s, is enhanced by about 30% or less compared to a standard halo model. Such a dark disk could contribute dominantly to the low energy (of order keV for a dark matter particle with mass 100 GeV) nuclear recoil event rate of direct dectection experiments, but it will not change the likelihood of detection significantly. These dark disks provide testable predictions of weakly-interacting massive particle dark matter models and should be considered in detailed comparisons to experimental data. Our findings suggest that the dark disk of the Milky Way may provide a detectable signal for indirect detection experiments, contributing up to about 25% of the dark matter self-annihilation signal in the direction of the center of the Galaxy, lending the signal a noticeably oblate morphology.

Can dark matter be stabilized by charge conservation, just as the electron is in the standard model? We examine the possibility that dark matter is hidden, that is, neutral under all standard model gauge interactions, but charged under an exact U(1) gauge symmetry of the hidden sector. Such candidates are predicted in WIMPless models, supersymmetric models in which hidden dark matter has the desired thermal relic density for a wide range of masses. Hidden charged dark matter has many novel properties not shared by neutral dark matter: (1) bound state formation and Sommerfeld-enhanced annihilation after chemical freeze out may reduce its relic density, (2) similar effects greatly enhance dark matter annihilation in protohalos at redshifts of z ~ 30, (3) Compton scattering off hidden photons delays kinetic decoupling, suppressing small scale structure, and (4) Rutherford scattering makes such dark matter self-interacting and collisional, potentially impacting properties of the Bullet Cluster and the observed morphology of galactic halos. We analyze all of these effects in a WIMPless model in which the hidden sector is a simplified version of the minimal supersymmetric standard model and the dark matter is a hidden sector stau. We find that charged hidden dark matter is viable and consistent with the correct relic density for reasonable model parameters and dark matter masses in the range 1 GeV < m_X < 10 TeV. At the same time, in the preferred range of parameters, this model predicts cores in the dark matter halos of small galaxies and other halo properties that may be within the reach of future observations. These models therefore provide a viable and well-motivated framework for collisional dark matter with Sommerfeld enhancement, with novel implications for astrophysics and dark matter searches.

We present a general methodology for determining the gamma-ray flux from annihilation of dark matter particles in Milky Way satellite galaxies, focusing on two promising satellites as examples: Segue 1 and Draco. We use the SuperBayeS code to explore the best-fitting regions of the Constrained Minimal Supersymmetric Standard Model (CMSSM) parameter space, and an independent MCMC analysis of the dark matter halo properties of the satellites using published radial velocities. We present a formalism for determining the boost from halo substructure in these galaxies and show that its value depends strongly on the extrapolation of the concentration-mass (c(M)) relation for CDM subhalos down to the minimum possible mass. We show that the preferred region for this minimum halo mass within the CMSSM with neutralino dark matter is ~10^-9-10^-6 solar masses. For the boost model where the observed power-law c(M) relation is extrapolated down to the minimum halo mass we find average boosts of about 20, while the Bullock et al (2001) c(M) model results in boosts of order unity. We estimate that for the power-law c(M) boost model and photon energies greater than a GeV, the Fermi space-telescope has about 20% chance of detecting a dark matter annihilation signal from Draco with signal-to-noise greater than 3 after about 5 years of observation.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Over the past five years, searches in Sloan Digital Sky Survey data have more than doubled the number of known dwarf satellite galaxies of the Milky Way, and have revealed a population of ultra-faint galaxies with luminosities smaller than typical globular clusters, L ~ 1000 Lsun. These systems are the faintest, most dark matter dominated, and most metal poor galaxies in the universe. Completeness corrections suggest that we are poised on the edge of a vast discovery space in galaxy phenomenology, with hundreds more of these extreme galaxies to be discovered as future instruments hunt for the low-luminosity threshold of galaxy formation. Dark matter dominated dwarfs of this kind probe the small-scale power-spectrum, provide the most stringent limits on the phase-space packing of dark matter, and offer a particularly useful target for dark matter indirect detection experiments. Full use of dwarfs as dark matter laboratories will require synergy between deep, large-area photometric searches; spectroscopic and astrometric follow-up with next-generation optical telescopes; and subsequent observations with gamma-ray telescopes for dark matter indirect detection.

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.

The hierarchical theory of galaxy formation rests on the idea that smaller galactic structures merge to form the galaxies that we see today. The past decade has provided remarkable observational support for this scenario, driven in part by advances in spectroscopic instrumentation. Multi-object spectroscopy enabled the discovery of kinematically cold substructures around the Milky Way and M31 that are likely the debris of disrupting satellites. Improvements in high-resolution spectroscopy have produced key evidence that the abundance patterns of the Milky Way halo and its dwarf satellites can be explained by Galactic chemical evolution models based on hierarchical assembly. These breakthroughs have depended almost entirely on observations of nearby stars in the Milky Way and luminous red giant stars in M31 and Local Group dwarf satellites. In the next decade, extremely large telescopes will allow observations far down the luminosity function in the known dwarf galaxies, and they will enable observations of individual stars far out in the Galactic halo. The chemical abundance census now available for the Milky Way will become possible for our nearest neighbor, M31. Velocity dispersion measurements now available in M31 will become possible for systems beyond the Local Group such as Sculptor and M81 Group galaxies. Detailed studies of a greater number of individual stars in a greater number of spiral galaxies and their satellites will test hierarchical assembly in new ways because dynamical and chemical evolution models predict different outcomes for halos of different masses in different environments.

Gravitational lensing of the cosmic microwave background by large-scale structure in the late universe is both a source of cosmological information and a potential contaminant of primordial gravity waves. Because lensing imprints growth of structure in the late universe on the CMB, measurements of CMB lensing will constrain parameters to which the CMB would not otherwise be sensitive, such as neutrino mass. If the instrumental noise is sufficiently small (<~ 5 uK-arcmin), the gravitational lensing contribution to the large-scale B-mode will be the limiting source of contamination when constraining a stochastic background of gravity waves in the early universe, one of the most exciting prospects for future CMB polarization experiments. High-sensitivity measurements of small-scale B-modes can reduce this contamination through a lens reconstruction technique that separates the lensing and primordial contributions to the B-mode on large scales. A fundamental design decision for a future CMB polarization experiment such as CMBpol is whether to have coarse angular resolution so that only the large-scale B-mode (and the large-scale E-mode from reionization) is measured, or high resolution to additionally measure CMB lensing. The purpose of this white paper is to evaluate the science case for CMB lensing in polarization: constraints on cosmological parameters, increased sensitivity to the gravity wave B-mode via lens reconstruction, expected level of contamination from non-CMB foregrounds, and required control of beam systematics.

The Milky Way has at least twenty-three known satellite galaxies that shine with luminosities ranging from about a thousand to a billion times that of the Sun. Half of these galaxies were discovered in the past few years in the Sloan Digital Sky Survey, and they are among the least luminous galaxies in the known Universe. A determination of the mass of these galaxies provides a test of galaxy formation at the smallest scales and probes the nature of the dark matter that dominates the mass density of the Universe. Here we use new measurements of the velocities of the stars in these galaxies to show that they are consistent with them having a common mass of about 10^7 M_\odot within their central 300 parsecs. This result demonstrates that the faintest of the Milky Way satellites are the most dark matter-dominated galaxies known, and could be a hint of a new scale in galaxy formation or a characteristic scale for the clustering of dark matter.

Proper motions for several Milky Way dwarf spheroidal (dSph) galaxies have been determined using both ground and space-based imaging. These measurements require long baselines and repeat observations and typical errors are of order ten milli-arcseconds per century. In this paper, we utilize the effect of "perspective rotation" to show that systematic proper motion of some dSphs can be determined to a similar precision using only stellar line-of-sight velocities. We show that including the effects of small intrinsic rotation in dSphs increases the proper motion errors by about a factor of two. We provide error projections for future data sets, and show that proposed thirty meter class telescopes will measure the proper motion of a few dSphs with milli-arcsecond per century precision.

We study the strong gravitational lens statistics of triaxial cold dark matter (CDM) halos occupied by central early-type galaxies. We calculate the image separation distribution for double, cusp and quad configurations. The ratios of image multiplicities at large separations are consistent with the triaxial NFW model, and at small separations are consistent with the singular isothermal ellipsoid (SIE) model. At all separations, the total lensing probability is enhanced by adiabatic contraction. If no adiabatic contraction is assumed, naked cusp configurations become dominant at approximately 2.5'', which is inconsistent with the data. We also show that at small-to-moderate separations, the image multiplicities depend sensitively on the alignment of the shapes of the luminous and dark matter projected density profiles. In constrast to other properties that affect these ratios, the degree of alignment does not have a significant effect on the total lensing probability. These correlations may therefore be constrained by comparing the theoretical image separation distribution to a sufficiently large lens sample from future wide and deep sky surveys such as Pan-Starrs, LSST and JDEM. Understanding the correlations in the shapes of galaxies and their dark matter halo is important for future weak lensing surveys.

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

We use kinematic data from three new, nearby, extremely low-luminosity Milky Way dwarf galaxies (Ursa Major II, Willman 1, and Coma Berenices) to constrain the properties of their dark matter halos, and from these make predictions for the gamma-ray flux from annihilation of dark matter particles in these halos. We show that these 10^3 solar luminosity dwarfs are the most dark matter dominated galaxies in the Universe, with total masses within 100 pc in excess of 10^6 solar masses. Coupled with their relative proximity, their large masses imply that they should have mean gamma-ray fluxes comparable to or greater than any other known satellite galaxy of the Milky Way. Our results are robust to both variations of the inner slope of the density profile and the effect of tidal interactions. The fluxes could be boosted by up to two orders of magnitude if we include the density enhancements caused by surviving dark matter substructure.

We discuss an instability in a large class of models where dark energy is coupled to matter. In these models the mass of the scalar field is much larger than the expansion rate of the Universe. We find models in which this instability is absent, and show that these models generically predict an apparent equation of state for dark energy smaller than -1, i.e., superacceleration. These models have no acausal behavior or ghosts.

We examine the prospects for detecting gamma-rays from dark matter annihilation in the six most promising dwarf spheroidal (dSph) satellite galaxies of the Milky Way. We use recently-measured velocity dispersion profiles to provide a systematic investigation of the dark matter mass distribution of each galaxy, and show that the uncertainty in the gamma-ray flux from mass modeling is less than a factor of ~ 5 for each dSph if we assume a smooth NFW profile. We show that Ursa Minor and Draco are the most promising dSphs for gamma-ray detection with GLAST and other planned observatories. For each dSph, we investigate the flux enhancement resulting from halo substructure, and show that the enhancement factor relative to a smooth halo flux cannot be greater than about 100. This enhancement depends very weakly on the lower mass cut-off scale of the substructure mass function. While the amplitude of the expected flux from each dSph depends sensitively on the dark matter model, we show that the flux ratios between the six Sphs are known to within a factor of about 10. The flux ratios are also relatively insensitive to the current theoretical range of cold dark matter halo central slopes and substructure fractions.

WIMPs with electroweak scale masses (neutralinos, etc.) remain in kinetic equilibrium with other particle species until temperatures approximately in the range of 10 MeV to 1 GeV, leading to the formation of dark matter substructure with masses as small as $10^{-4} M_{\odot}$ to $10^{-12} M_{\odot}$. However, if dark matter consists of particles with MeV scale masses, as motivated by the observation of 511 keV emission from the Galactic Bulge, such particles are naturally expected to remain in kinetic equilibrium with the cosmic neutrino background until considerably later times. This would lead to a strong suppression of small scale structure with masses below about $10^7 M_{\odot}$ to $10^4 M_{\odot}$. This cutoff scale has important implications for present and future searches for faint Local Group satellite galaxies and for the missing satellites problem.

We show that dark matter emerging from late decays (z < 1000) produces a linear power spectrum identical to that of Cold Dark Matter (CDM) on all observationally relavant scales (> 0.1 Mpc), and simultaneously generates observable constant-density cores in small dark matter halos. We refer to this class of models as meta-Cold Dark Matter (mCDM), because it is born with non-relativistic velocities from the decays of cold thermal relics. The constant-density cores are a result of the low phase-space density of mCDM at birth. Warm dark matter cannot produce similar size phase-space limited cores without saturating the Ly-alpha power spectrum bounds. Dark matter dominated galaxy rotation curves and stellar velocity dispersion profiles may provide the best means to discriminate between mCDM and CDM. mCDM candidates are motivated by the particle spectrum of supersymmetric and extra dimensional extensions to the standard model of particle physics.

B-modes of Cosmic Microwave Background (CMB) polarization can be created by a primordial gravitational wave background. If this background was created by Inflation, then the amplitude of the polarization signal is proportional the energy density of the universe during inflation. The primordial signal will be contaminated by polarized foregrounds including dust and synchrotron emission within the galaxy. In light of the WMAP polarization maps, we consider the ability of several hypothetical CMB polarization experiments to separate primordial CMB B-mode signal from galactic foregrounds. We also study the optimization of a CMB experiment with a fixed number of detectors in the focal plane to determine how the detectors should be distributed in different frequency bands to minimize foreground confusion. We show that the optimal configuration requires observations in at least 5 channels spread over the frequency range between 30 GHz and 500 GHz with substantial coverage around 150 GHz. If a low-resolution space experiment using 1000 detectors to reach a noise level of about 1000 nK$^2$ concentrates on roughly 66% of the sky with the least foreground contamination the minimum detectable level of the tensor-to-scalar ratio would be about 0.002 at the 99% confidence level for an optical depth of 0.1.

We show that measurements of stellar proper motions in dwarf spheroidal galaxies provide a powerful probe of the nature of dark matter. Allowing for general dark matter density profiles and stellar velocity anisotropy profiles, we show that the log-slope of the dark matter profile at about twice the stellar core (King) radius can be measured to within \pm 0.2 when the proper motions of 200 stars are added to standard line-of-sight velocity dispersion data. This measurement of the log-slope provides a test of Cold and Warm Dark Matter theories at a sensitivity not possible with line-of-sight velocity dispersion measurements alone. The upcoming SIM PlanetQuest will have the sensitivity to obtain the required number of proper motions in Milky Way dwarf spheroidal galaxies.

Charged particles (X) decaying after primordial nucleosynthesis are constrained by the requirement that their decay products should not change the light element abundances drastically. If the decaying particle is negatively charged (X-) then it will bind to the nuclei. We consider the effects of the decay of X when bound to Helium-4 and show that this will modify the Lithium abundances.

Gravitational lensing distorts the cosmic microwave background (CMB) temperature and polarization fields and encodes valuable information on distances and growth rates at intermediate redshifts into the lensed power spectra. The non-Gaussian bandpower covariance induced by the lenses is negligible to l=2000 for all but the B polarization field where it increases the net variance by up to a factor of 10 and favors an observing strategy with 3 times more area than if it were Gaussian. To quantify the cosmological information, we introduce two lensing observables, characterizing nearly all of the information, which simplify the study of non-Gaussian impact, parameter degeneracies, dark energy models, and complementarity with other cosmological probes. Information on the intermediate redshift parameters rapidly becomes limited by constraints on the cold dark matter density and initial amplitude of fluctuations as observations improve. Extraction of this information requires deep polarization measurements on only 5-10% of the sky, and can improve Planck lensing constraints by a factor of ~2-3 on any one of the parameters w_0, w_a, Omega_K, sum(m_nu) with the others fixed. Sensitivity to the curvature and neutrino mass are the highest due to the high redshift weight of CMB lensing but degeneracies between the parameters must be broken externally.

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.

The search for the first illuminated astronomical sources in the universe is at the edge of the cosmic frontier. Promising techniques for discovering the first objects and their effects span the electromagnetic spectrum and include gravitational waves. We summarize a workshop on discovering and understanding these sources which was held in May 2005 through the Center for Cosmology at the University of California, Irvine.

Two leading dark matter candidates from supersymmetry and other theories of physics beyond the standard model are WIMPs and weak scale gravitinos. If the lightest stable particle is a gravitino, then a WIMP will decay into it with a natural lifetime of order a month ~ M_{pl}^2/M_{weak}^3. We show that if the bulk of dark matter today came from decays of neutral particles with lifetimes of order a year or smaller, then it could lead to a reduction in the amount of small scale substructure, less concentrated halos and constant density cores in the smallest mass halos. Such beneficial effects may therefore be realized naturally, as discussed by Cembranos, Feng, Rajaraman, and Takayama, in the case of supersymmetry.

We propose a method to probe the phenomenological nature of dark energy which makes no assumptions about the evolution of its energy density. We exemplify this method with a test for a super-acceleration phase of the universe i.e., a phase when the dark energy density grows as the universe expands. We show how such a phase can be detected by combining SNIa (SNAP-like) and CMB (Planck) data without making any assumptions about the evolution of the dark energy equation of state, or about the value of the matter density parameter.

In 2002, Ray Davis and Masatoshi Koshiba were awarded the Nobel Prize in Physics ``for pioneering contributions to astrophysics, in particular for the detection of cosmic neutrinos.'' However, while astronomy has undergone a revolution in understanding by synthesizing data taken at many wavelengths, the universe has only barely been glimpsed in neutrinos, just the Sun and the nearby SN 1987A. An entire universe awaits, and since neutrinos can probe astrophysical objects at densities, energies, and distances that are otherwise inaccessible, the results are expected to be particularly exciting. Similarly, the revolution in quantitative cosmology has heightened the need for very precise tests that depend on the effects of neutrinos, and prominent among them is the search for the effects of neutrino mass, since neutrinos are a small but known component of the dark matter. In this report, we highlight some of the key opportunties for progress in neutrino astrophysics and cosmology, and the implications for other areas of physics.

We calculate the Supernova Relic Neutrino (SRN) background flux for the KamLAND and Super-Kamiokande (SK) detectors, motivated by the reduction in background at SK and new results for the star formation history (e.g., from the Sloan Digital Sky Survey (SDSS)). Our best estimate for the flux at SK is slightly below, but very close to the current SK upper limit. The SK upper limit is already inconsistent with a range of star formation histories allowed by the SDSS data. We estimate that the SRN background should be detected (at 1-sigma) at SK with a total of about 9 years (including the existing 4 years) of data. While KamLAND is a much smaller detector compared to SK, it profits from being practically background-free and from its sensitivity to the lower energy supernova neutrinos. KamLAND could make a 1-sigma detection of the SRN with a total of about 5 years of data. Given the small expected SRN event rate, we also consider the detection of the SRN in a modified SK detector with a lower threshold and reduced background where the time to detection can be reduced by a factor of 10 relative to the existing SK estimate.

The B modes generated by the lensing of CMB polarization are a primary target for the upcoming generation of experiments and can potentially constrain quantities such as the neutrino mass and dark energy equation of state. The net sample variance on the small scale B modes out to l=2000 exceeds Gaussian expectations by a factor of 10 reflecting the variance of the larger scale lenses that generate them. It manifests itself as highly correlated band powers with correlation coefficients approaching 70% for wide bands of Delta l/l \~0.25. It will double the total variance for experiments that achieve a sensitivity of approximately 4 uK-arcmin and a beam of several arcminutes or better. This non-Gaussianity must be taken into account in the analysis of experiments that go beyond first detection.

Distortions of Cosmic Microwave Background (CMB) temperature and polarization maps caused by gravitational lensing, observable with high angular resolution and high sensitivity, can be used to measure the neutrino mass. Assuming two massless species and one with mass m_nu we forecast sigma(m_nu) = 0.15 eV from the Planck satellite and sigma(m_nu)=0.04 eV from observations with twice the angular resolution and about 20 times the sensitivity. A detection is likely at this higher sensitivity since the observation of atmospheric neutrino oscillations require mass-squared differences greater than about (0.04 eV)^2.

We show that deviations of the quantum state of the inflaton from the thermal vacuum of inflation may leave an imprint in the CMB anisotropies. The quantum dynamics of the inflaton in such a state produces corrections to the inflationary fluctuations, which may be observable. Because these effects originate from IR physics below the Planck scale, they will dominate over any trans-Planckian imprints in any theory which obeys decoupling. Inflation sweeps away these initial deviations and forces its quantum state closer to the thermal vacuum. We view this as the quantum version of the cosmic no-hair theorem. Such imprints in the CMB may be a useful, independent test of the duration of inflation, or of significant features in the inflaton potential about 60 e-folds before inflation ended, instead of an unlikely discovery of the signatures of quantum gravity. The absence of any such substructure would suggest that inflation lasted uninterrupted much longer than ${\cal O}(100)$ e-folds.

With WMAP putting the phenomenological standard model of cosmology on a strong footing, one can look forward to mining the cosmic microwave background (CMB) for fundamental physics with higher sensitivity and on smaller scales. Future CMB observations have the potential to measure absolute neutrino masses, test for cosmic acceleration independent of supernova Ia observations, probe for the presence of dark energy at redshifts of 2 and larger, illuminate the end of the dark ages, measure the scale--dependence of the primordial power spectrum and detect gravitational waves generated by inflation.

The normal parameters are a non--linear transformation of the cosmological parameters whose likelihood function is very well--approximated by a normal distribution. This transformation serves as an extreme form of data compression allowing for practically instantaneous calculation of the likelihood of any given model, as long as the model is in the parameter space originally considered. The compression makes all the information about cosmological parameter constraints from a given set of experiments available in a useable manner. Here we explicitly define the normal parameters that work for the current CMB data, and give their mean and covariance matrix which best fit the likelihood function calculated by the Monte Carlo Markov Chain method. Along with standard parameter estimation results, we propose that future CMB parameter analyses define normal parameters and quote their mean and covariance matrix.

The Davis Meeting on Cosmic Inflation marked an exciting milestone on the road to precision cosmology. This is the index page for the proceedings of the conference. Individual proceedings contributions, when they appear on this archive, are linked from this page.

In light of the recent inference of a high optical depth to Thomson scattering, tau, from the WMAP data we investigate the effects of extended periods of partial ionization and ask if the value of tau inferred by assuming a single sharp transition is an unbiased estimate. We construct and consider several representative ionization models and evaluate their signatures in the CMB. If tau is estimated with a single sharp transition we show that there can be a significant bias in the derived value (and therefore a bias in sigma8 as well). For WMAP noise levels the bias in tau is smaller than the statistical uncertainty, but for Planck or a cosmic variance limited experiment the tau bias could be much larger than the statistical uncertainties. This bias can be reduced in the ionization models we consider by fitting a slightly more complicated ionization history, such as a two-step ionization process. Assuming this two-step process we find the Planck satellite can simultaneously determine the initial redshift of reionization to +-2 and tau to +-0.01 Uncertainty about the ionization history appears to provide a limit of about 0.005 on how well tau can be estimated from CMB polarization data, much better than expected from WMAP but significantly worse than expected from cosmic-variance limits.

The recent discovery of a Gunn--Peterson (GP) trough in the spectrum of the redshift 6.28 SDSS quasar has raised the tantalizing possibility that we have detected the reionization of the universe. However, a neutral fraction (of hydrogen) as small as 0.1% is sufficient to cause the GP trough, hence its detection alone cannot rule out reionization at a much earlier epoch. The Cosmic Microwave Background (CMB) polarization anisotropy offers an alternative way to explore the dark age of the universe. We show that for most models constrained by the current CMB data and by the discovery of a GP trough (showing that reionization occurred at z > 6.3), MAP can detect the reionization signature in the polarization power spectrum. The expected 1-sigma error on the measurement of the electron optical depth is around 0.03 with a weak dependence on the value of that optical depth. Such a constraint on the optical depth will allow MAP to achieve a 1-sigma error on the amplitude of the primordial power spectrum of 6%. MAP with two years (Planck with one year) of observation can distinguish a model with 50% (6%) partial ionization between redshifts of 6.3 and 20 from a model in which hydrogen was completely neutral at redshifts greater than 6.3. Planck will be able to distinguish between different reionization histories even when they imply the same optical depth to electron scattering for the CMB photons.

We have developed a fast method for predicting the angular power spectrum, C_l, of the cosmic microwave background given cosmological parameters and a primordial power spectrum of perturbations. After pre--computing the radiation temperature and gravitational potential transfer functions over a small sub--space of the total model parameter space, the rest of the model space (six or more cosmological parameters and arbitrarily many primordial power spectrum parameters) is reached via rapid analytic and semi--analytic approximations which are highly accurate on all angular scales for which linear perturbation theory applies. A single power spectrum can be calculated in ~ 1 second on a desktop computer. We discuss applications to cosmological parameter estimation.

Milgrom noticed the remarkable fact that the gravitational effect of dark matter in galaxies only becomes important where accelerations are less than about 10^{-8} cm s^{-2} ~ cH_0. This forms the basis for his Modified Newtonian Dynamics (MOND), an alternative to particle dark matter. However, any successful theory of galactic dynamics must account for Milgrom's Law. We show how Milgrom's Law comes about in the Cold Dark Matter (CDM) theory of structure formation.

In conventional general relativity, the expansion rate H of a Robertson-Walker universe is related to the energy density by the Friedmann equation. Aside from the present day, the only epoch at which we can constrain the expansion history in a model-independent way is during Big-Bang Nucleosynthesis (BBN). We consider a simple two-parameter characterization of the behavior of H during BBN and derive constraints on this parameter space, finding that the allowed region of parameter space is essentially one-dimensional. We also study the effects of a large neutrino asymmetry within this framework. Our results provide a simple way to compare an alternative cosmology to the observational requirement of matching the primordial abundances of the light elements.

Recent cosmological observations strongly suggest that the universe is dominated by an unknown form of energy with negative pressure. Why is this dark energy density of order the critical density today? We propose that the dark energy has periodically dominated in the past so that its preponderance today is natural. We illustrate this paradigm with a model potential and show that its predictions are consistent with all observations.

Big-bang Nucleosynthesis (BBN) and Cosmic Microwave Background (CMB) anisotropy measurements give independent, accurate measurements of the baryon density and can test the framework of the standard cosmology. Early CMB data are consistent with the longstanding conclusion from BBN that baryons constitute a small fraction of matter in the Universe, but may indicate a slightly higher value for the baryon density. We clarify precisely what the two methods determine, and point out that differing values for the baryon density can indicate either an inconsistency or physics beyond the standard models of cosmology and particle physics. We discuss other signatures of the new physics in CMB anisotropy.

Structure formation with cold dark matter (CDM) predicts halos with a central density cusp, which are observationally disfavored. If CDM particles have an annihilation cross section sigma*v ~ 10**(-29) (m/GeV) cm**2, then annihilations will soften the cusps. We discuss plausible scenarios for avoiding the early Universe annihilation catastrophe that could result from such a large cross section. The predicted scaling of core density with halo mass depends upon the velocity dependence of sigma*v, and s-wave annihilation leads to a core density nearly independent of halo mass, which seems consistent with observations.

An upper bound to the supernova relic neutrino background from all past Type II supernovae is obtained using observations of the Universal metal enrichment history. We show that an unambiguous detection of these relic neutrinos by the Super-Kamiokande detector is unlikely. We also analyze the event rate in the Sudbury Neutrino Observatory (where coincident neutrons from anti-nu_e + D --> n + n + e+ might enhance background rejection), and arrive at the same conclusion. If the relic neutrino flux should be observed to exceed our upper bound and if the observations of the metal enrichment history (for z<1) are not in considerable error, then either the Type II supernova rate does not track the metal enrichment history or some mechanism may be responsible for transforming anti-nu_{mu,tau} --> anti-nu_e.

We have recently considered cosmologies in which the Universal scale factor varies as a power of the age of the Universe and concluded that they cannot satisfy the observational constraints on the present age, the magnitude-redshift relation for SN Ia, and the primordial element (D, He3, He4, and Li7) abundances. This claim has been challenged in a proposal that suggested a high baryon density model (Omega_B*h*h = 0.3) with an expansion factor varing linearly with time could be consistent with the observed abundance of primoridal helium-4, while satisfying the age and magnitude-redshift constraints. In this paper we further explore primordial nucleosynthesis in generic power-law cosmologies, including the linear case, concluding that models selected to satisfy the other observational constraints are incapable of accounting for all the light element abundances.

A massive neutrino which decays after recombination (t > 10^{13} sec) into relativistic decay products produces an enhanced integrated Sachs-Wolfe effect, allowing constraints to be placed on such neutrinos from present cosmic microwave background anisotropy data. Previous treatments of this problem have approximated the decay products as an additional component of the neutrino background. This approach violates energy-momentum conservation, and we show that it leads to serious errors for some neutrino masses and lifetimes. We redo this calculation more accurately, by correctly incorporating the spatial distribution of the decay products. For low neutrino masses and long lifetimes, we obtain a much smaller distortion in the CMB fluctuation spectrum than have previous treatments. We combine these new results with a recent set of CMB data to exclude the mass and lifetime range m_h > 100 eV, \tau > 10^{12} sec. Masses as low as 30 eV are excluded for a narrower range in lifetime.

Any time variation in the fine-structure constant alters the ionization history of the universe and therefore changes the pattern of cosmic microwave background fluctuations. We calculate the changes in the spectrum of these fluctuations as a function of the change in the fine-structure constant, and we find that these changes are dominated by the change in the redshift of recombination due to the shift in the binding energy of hydrogen. We estimate the accuracy with which the next generation of cosmic microwave background experiments might constrain any variation in the fine-structure constant at z of roughly 1000. We find that such experiments could potentially be sensitive to $|\Delta \alpha/\alpha| \sim 10^{-2} - 10^{-3}$.

In a class of models designed to solve the cosmological constant problem by coupling scalar or tensor classical fields to the space-time curvature, the universal scale factor grows as a power law in the age, $a \propto t^\alpha$, regardless of the matter content or cosmological epoch. We investigate constraints on such "power-law cosmologies" from the present age of the Universe, the magnitude-redshift relation, and from primordial nucleosynthesis. Constraints from the current age of the Universe and from the high-redshift supernovae data require "large" $\alpha$ ($\approx 1$), while consistency with the inferred primordial abundances of deuterium and helium-4 forces $\alpha$ to lie in a very narrow range around a lower value ($\approx 0.55$). Inconsistency between these independent cosmological constraints suggests that such power-law cosmologies are not viable.