I show that model dependencies in the analysis by Dessert, Rodd & Safdi (2020) relax their claimed constraint by a factor of ~20. After including conservative model choices, the derived limits are comparable to or slightly better than limits from previous searches. Further model tests and expansion of the data energy may enhance or relax sensitivity of the methodology.

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 assess Galactic Dark Matter (DM) sensitivities to photons from annihilation and decay using the spatial and kinematic information determined by state-of-the-art simulations in the Latte suite of Feedback In Realistic Environments (FIRE-2). For kinematic information, we study the energy shift pattern of DM narrow emission lines predicted in FIRE-2 and discuss its potential as DM-signal diagnosis, showing for the first time the power of symmetric observations around $l=0^{\circ}$. We find that the exposures needed to resolve the line separation of DM to gas by XRISM at $5\sigma$ to be large, $\gtrsim 4$ Ms, while exposures are smaller for Athena ($\lesssim 50$ ks) and Lynx ($\lesssim 100$ ks). We find that large field-of-view exposures remain the most sensitive methods for detection of DM annihilation or decay by the luminosity of signals in the field of view dominating velocity information. The $\sim$4 sr view of the Galactic Center region by the Wide Field Monitor (WFM) aboard the eXTP mission will be highly sensitive to DM signals, with a prospect of $\sim 10^5$ to $10^6$ events from the 3.5 keV line in a 100 ks exposure, with the range dependent on photon acceptance in WFM's field of view. We also investigate detailed all-sky luminosity maps for both DM annihilation and decay signals - evaluating the signal-to-noise for a DM detection with realistic X-ray and gamma-ray backgrounds - as a guideline for what could be a forthcoming era of DM astronomy.

Gamma-ray observations have revealed strong variability in blazar luminosities in the gamma-ray band over time scales as short as minutes. We show, for the first time, that the correlation of the spectrum with intensity is consistent with the behavior of the luminosity variation of blazar SEDs along a blazar sequence for low synchrotron peak blazars. We show that the observational signatures of variability with flux are consistent with wakefield acceleration of electrons initiated by instabilities in the blazar accretion disk. This mechanism reproduces the observed time variations as short as 100 seconds. The wakefield mechanism also predicts a reduction of the electron spectral index with increased gamma-ray luminosity, which could be detected in higher energy observations well above the inverse Compton peak.

We discuss numerous mechanisms for production of sterile neutrinos, which can account for all or a fraction of dark matter, and which can range from warm to effectively cold dark matter, depending on the cosmological scenario. We investigate production by Higgs boson decay, $(B-L)$ gauge boson production at high temperature, as well as production via resonant and nonresonant neutrino oscillations. We calculate the effects on structure formation in these models, some for the first time. If two populations of sterile neutrinos, one warm and one cold, were produced by different mechanisms, or if sterile neutrinos account for only a fraction of dark matter, while the remainder is some other cold dark matter particle, the resulting multi-component dark matter may alleviate some problems in galaxy formation. We examine the X-ray constraints and the candidate signal at 3.5 keV. Finally, we also show that the $\sigma_8$ problem can be a signature of fractional dark matter in the form of sterile neutrinos in several mechanisms.

Planned CMB Stage IV experiments have the potential to measure the effective number of relativistic degrees of freedom in the early Universe, $N_\text{eff}$, with percent-level accuracy. This probes new thermalized light particles and also constrains possible new-physics interactions of Dirac neutrinos. Many Dirac-neutrino models that aim to address the Dirac stability, the smallness of neutrino masses or the matter--anti-matter asymmetry of our Universe endow the right-handed chirality partners $\nu_R$ with additional interactions that can thermalize them. Unless the reheating temperature of our Universe was low, this leads to testable deviations in $N_\text{eff}$. We discuss well-motivated models for $\nu_R$ interactions such as gauged $U(1)_{B-L}$ and the neutrinophilic two-Higgs-doublet model, and compare the sensitivity of SPT-3G, Simons Observatory, and CMB-S4 to other experiments, in particular the LHC.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

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 study the impact of a warm dark matter (WDM) cosmology on dwarf galaxy formation through a suite of cosmological hydrodynamical zoom-in simulations of $M_{\rm halo} \approx10^{10}\,M_{\odot}$ dark matter halos as part of the Feedback in Realistic Environments (FIRE) project. A main focus of this paper is to evaluate the combined effects of dark matter physics and stellar feedback on the well-known small-scale issues found in cold dark matter (CDM) models. We find that the $z=0$ stellar mass of a galaxy is strongly correlated with the central density of its host dark matter halo at the time of formation, $z_{\rm f}$, in both CDM and WDM models. WDM halos follow the same $M_{\star}(z=0)-V_{\rm max}(z_{\rm f})$ relation as in CDM, but they form later, are less centrally dense, and therefore contain galaxies that are less massive than their CDM counterparts. As a result, the impact of baryonic effects on the central gravitational potential is typically diminished relative to CDM. However, the combination of delayed formation in WDM and energy input from stellar feedback results in dark matter profiles with lower overall densities. The WDM galaxies studied here have a wider diversity of star formation histories (SFHs) than the same systems simulated in CDM, and the two lowest $M_{\star}$ WDM galaxies form all of their stars at late times. The discovery of young ultra-faint dwarf galaxies with no ancient star formation -- which do not exist in our CDM simulations -- would therefore provide evidence in support of WDM.

Sterile neutrinos are natural extensions to the standard model of particle physics in neutrino mass generation mechanisms. If they are relatively light, less than approximately 10 keV, they can alter cosmology significantly, from the early Universe to the matter and radiation energy density today. Here, we review the cosmological role such light sterile neutrinos can play from the early Universe, including production of keV-scale sterile neutrinos as dark matter candidates, and dynamics of light eV-scale sterile neutrinos during the weakly-coupled active neutrino era. We review proposed signatures of light sterile neutrinos in cosmic microwave background and large scale structure data. We also discuss keV-scale sterile neutrino dark matter decay signatures in X-ray observations, including recent candidate $\sim$3.5 keV X-ray line detections consistent with the decay of a $\sim$7 keV sterile neutrino dark matter particle.

The Milky Way's Galactic Center harbors a gamma-ray excess that is a candidate signal of annihilating dark matter. Dwarf galaxies remain predominantly dark in their expected commensurate emission. In this work we quantify the degree of consistency between these two observations through a joint likelihood analysis. In doing so we incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma-Ray Space Telescope. The preferred range of annihilation rates and masses expands when including these unknowns. Even so, using two recent determinations of the Milky Way halo's local density leave the GCE preferred region of single-channel dark matter annihilation models to be in strong tension with annihilation searches in combined dwarf galaxy analyses. A third, higher Milky Way density determination, alleviates this tension. Our joint likelihood analysis allows us to quantify this inconsistency. We provide a set of tools for testing dark matter annihilation models' consistency within this combined dataset. As an example, we test a representative inverse Compton sourced self-interacting dark matter model, which is consistent with both the GCE and dwarfs.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

We test the sensitivity of neutrino parameter constraints from combinations of CMB and LSS data sets to the assumed form of the primordial power spectrum (PPS) using Bayesian model selection. Significantly, none of the tested combinations, including recent high-precision local measurements of $\mathrm{H}_0$ and cluster abundances, indicate a signal for massive neutrinos or extra relativistic degrees of freedom. For PPS models with a large, but fixed number of degrees of freedom, neutrino parameter constraints do not change significantly if the location of any features in the PPS are allowed to vary, although neutrino constraints are more sensitive to PPS features if they are known a priori to exist at fixed intervals in $\log k$. Although there is no support for a non-standard neutrino sector from constraints on both neutrino mass and relativistic energy density, we see surprisingly strong evidence for features in the PPS when it is constrained with data from Planck 2015, SZ cluster counts, and recent high-precision local measurements of $\mathrm{H}_0$. Conversely combining Planck with matter power spectrum and BAO measurements yields a much weaker constraint. Given that this result is sensitive to the choice of data this tension between SZ cluster counts, Planck and $\mathrm{H}_0$ measurements is likely an indication of unmodeled systematic bias that mimics PPS features, rather than new physics in the PPS or neutrino sector.

We incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma Ray Space Telescope. The range of particle annihilation rate and masses expand when including these unknowns. However, two of the most precise empirical determinations of the Milky Way halo's local density and density profile leave the signal region to be in considerable tension with dark matter annihilation searches from combined dwarf galaxy analyses for single-channel dark matter annihilation models. The GCE and dwarf tension can be alleviated if: one, the halo is very highly concentrated or strongly contracted; two, the dark matter annihilation signal differentiates between dwarfs and the GC; or, three, local stellar density measures are found to be significantly lower, like that from recent stellar counts, increasing the local dark matter density.

The anomalous 3.55 keV X-ray line recently detected towards a number of massive dark matter objects may be interpreted as the radiative decays of 7.1 keV mass sterile neutrino dark matter. Depending on its parameters, the sterile neutrino can range from cold to warm dark matter with small-scale suppression that differs in form from commonly-adopted thermal warm dark matter. Here, we numerically investigate the subhalo properties for 7.1 keV sterile neutrino dark matter produced via the resonant Shi-Fuller mechanism. Using accurate matter power spectra, we run cosmological zoom-in simulations of a Milky Way-sized halo and explore the abundance of massive subhalos, their radial distributions, and their internal structure. We also simulate the halo with thermal 2.0 keV warm dark matter for comparison and discuss quantitative differences. We find that the resonantly produced sterile neutrino model for the 3.55 keV line provides a good description of structures in the Local Group, including the number of satellite dwarf galaxies and their radial distribution, and largely mitigates the too-big-to-fail problem. Future searches for satellite galaxies by deep surveys, such as the Dark Energy Survey, Large Synoptic Survey Telescope, and Wide Field Infrared Survey Telescope, will be a strong direct test of warm dark matter scenarios.

Sterile neutrinos comprise an entire class of dark matter models that, depending on their production mechanism, can be hot, warm, or cold dark matter. We simulate the Local Group and representative volumes of the Universe in a variety of sterile neutrino models, all of which are consistent with the possible existence of a radiative decay line at ~3.5 keV. We compare models of production via resonances in the presence of a lepton asymmetry (suggested by Shi & Fuller 1999) to "thermal" models. We find that properties in the highly nonlinear regime - e.g., counts of satellites and internal properties of halos and subhalos - are insensitive to the precise fall-off in power with wavenumber, indicating that nonlinear evolution essentially washes away differences in the initial (linear) matter power spectrum. In the quasi-linear regime at higher redshifts, however, quantitative differences in the 3D matter power spectra remain, raising the possibility that such models can be tested with future observations of the Lyman-alpha forest. While many of the sterile neutrino models largely eliminate multiple small-scale issues within the Cold Dark Matter (CDM) paradigm, we show that these models may be ruled out in the near future via discoveries of additional dwarf satellites in the Local Group.

We perform a detailed study of the weak interactions of standard model neutrinos with the primordial plasma and their effect on the resonant production of sterile neutrino dark matter. Motivated by issues in cosmological structure formation on small scales, and reported X-ray signals that could be due to sterile neutrino decay, we consider $7$ keV-scale sterile neutrinos. Oscillation-driven production of such sterile neutrinos occurs at temperatures $T \gtrsim 100$ MeV, where we study two significant effects of weakly charged species in the primordial plasma: (1) the redistribution of an input lepton asymmetry; (2) the opacity for active neutrinos. We calculate the redistribution analytically above and below the quark-hadron transition, and match with lattice QCD calculations through the transition. We estimate opacities due to tree level processes involving leptons and quarks above the quark-hadron transition, and the most important mesons below the transition. We report final sterile neutrino dark matter phase space densities that are significantly influenced by these effects, and yet relatively robust to remaining uncertainties in the nature of the quark-hadron transition. We also provide transfer functions for cosmological density fluctuations with cutoffs at $k \simeq 10 \ h \ {\rm Mpc}^{-1}$, that are relevant to galactic structure formation.

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.

The High Altitude Water Cherenkov (HAWC) observatory is a wide field-of-view detector sensitive to gamma rays of 100 GeV to a few hundred TeV. Located in central Mexico at 19 degrees North latitude and 4100 m above sea level, HAWC will observe gamma rays and cosmic rays with an array of water Cherenkov detectors. The full HAWC array is scheduled to be operational in Spring 2015. In this paper, we study the HAWC sensitivity to the gamma-ray signatures of high-mass (multi- TeV) dark matter annihilation. The HAWC observatory will be sensitive to diverse searches for dark matter annihilation, including annihilation from extended dark matter sources, the diffuse gamma-ray emission from dark matter annihilation, and gamma-ray emission from non-luminous dark matter subhalos. Here we consider the HAWC sensitivity to a subset of these sources, including dwarf galaxies, the M31 galaxy, the Virgo cluster, and the Galactic center. We simulate the HAWC response to gamma rays from these sources in several well-motivated dark matter annihilation channels. If no gamma-ray excess is observed, we show the limits HAWC can place on the dark matter cross-section from these sources. In particular, in the case of dark matter annihilation into gauge bosons, HAWC will be able to detect a narrow range of dark matter masses to cross-sections below thermal. HAWC should also be sensitive to non-thermal cross-sections for masses up to nearly 1000 TeV. The constraints placed by HAWC on the dark matter cross-section from known sources should be competitive with current limits in the mass range where HAWC has similar sensitivity. HAWC can additionally explore higher dark matter masses than are currently constrained.

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.

Using the temperature data from Planck we search for departures from a power-law primordial power spectrum, employing Bayesian model-selection and posterior probabilities. We parametrize the spectrum with $n$ knots located at arbitrary values of $\log{k}$, with both linear and cubic splines. This formulation recovers both slow modulations and sharp transitions in the primordial spectrum. The power spectrum is well-fit by a featureless, power-law at wavenumbers $k>10^{-3} \, \mathrm{Mpc}^{-1}$. A modulated primordial spectrum yields a better fit relative to $\Lambda$CDM at large scales, but there is no strong evidence for a departure from a power-law spectrum. Moreover, using simulated maps we show that a local feature at $k \sim 10^{-3} \, \mathrm{Mpc}^{-1}$ can mimic the suppression of large-scale power. With multi-knot spectra we see only small changes in the posterior distributions for the other free parameters in the standard $\Lambda$CDM universe. Lastly, we investigate whether the hemispherical power asymmetry is explained by independent features in the primordial power spectrum in each ecliptic hemisphere, but find no significant differences between them.

An inflationary gravitational wave background consistent with BICEP2 is difficult to reconcile with a simple power-law spectrum of primordial scalar perturbations. Tensor modes contribute to the temperature anisotropies at multipoles with $l\lesssim 100$, and this effect --- together with a prior on the form of the scalar perturbations --- was the source of previous bounds on the tensor-to-scalar ratio. We compute Bayesian evidence for combined fits to BICEP2 and Planck for three nontrivial primordial spectra: a) a running spectral index, b) a cutoff at fixed wavenumber, and c) a spectrum described by a linear spline with a single internal knot. We find no evidence for a cutoff, weak evidence for a running index, and significant evidence for a "broken" spectrum. Taken at face-value, the BICEP2 results require two new inflationary parameters in order to describe both the broken scale invariance in the perturbation spectrum and the observed tensor-to-scalar ratio. Alternatively, this tension may be resolved by additional data and more detailed analyses.

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.

This is a report on the status and prospects of the quantification of neutrino properties through the cosmological neutrino background for the Cosmic Frontier of the Division of Particles and Fields Community Summer Study long-term planning exercise. Experiments planned and underway are prepared to study the cosmological neutrino background in detail via its influence on distance-redshift relations and the growth of structure. The program for the next decade described in this document, including upcoming spectroscopic galaxy surveys eBOSS and DESI and a new Stage-IV CMB polarization experiment CMB-S4, will achieve sigma(sum m_nu) = 16 meV and sigma(N_eff) = 0.020. Such a mass measurement will produce a high significance detection of non-zero sum m_nu, whose lower bound derived from atmospheric and solar neutrino oscillation data is about 58 meV. If neutrinos have a minimal normal mass hierarchy, this measurement will definitively rule out the inverted neutrino mass hierarchy, shedding light on one of the most puzzling aspects of the Standard Model of particle physics --- the origin of mass. This precise a measurement of N_eff will allow for high sensitivity to any light and dark degrees of freedom produced in the big bang and a precision test of the standard cosmological model prediction that N_eff = 3.046.

Fluctuations in the intensity and polarization of the cosmic microwave background (CMB) and the large-scale distribution of matter in the universe each contain clues about the nature of the earliest moments of time. The next generation of CMB and large-scale structure (LSS) experiments are poised to test the leading paradigm for these earliest moments---the theory of cosmic inflation---and to detect the imprints of the inflationary epoch, thereby dramatically increasing our understanding of fundamental physics and the early universe. A future CMB experiment with sufficient angular resolution and frequency coverage that surveys at least 1% of the sky to a depth of 1 uK-arcmin can deliver a constraint on the tensor-to-scalar ratio that will either result in a 5-sigma measurement of the energy scale of inflation or rule out all large-field inflation models, even in the presence of foregrounds and the gravitational lensing B-mode signal. LSS experiments, particularly spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will complement the CMB effort by improving current constraints on running of the spectral index by up to a factor of four, improving constraints on curvature by a factor of ten, and providing non-Gaussianity constraints that are competitive with the current CMB bounds.

We use $N$-body simulations to show that high-redshift galaxy counts provide an interesting constraint on the nature of dark matter, specifically Warm Dark Matter (WDM), owing to the lack of early structure formation these models. Our simulations include three WDM models with thermal-production masses of 0.8 keV, 1.3 keV, and 2.6 keV, as well as CDM. Assuming a relationship between dark halo mass and galaxy luminosity that is set by the observed luminosity function at bright magnitudes, we find that 0.8 keV WDM is disfavored by direct galaxy counts in the Hubble Ultra Deep Field at $>\!\!10\sigma$. Similarly, 1.3 keV WDM is statistically inconsistent at $2.2\sigma$. Future observations with JWST (and possibly HST via the Frontier Fields) could rule out $1.3$ keV WDM at high significance, and may be sensitive to WDM masses greater than 2.6 keV. We also examine the ability of galaxies in these WDM models to reionize the universe, and find that 0.8 keV and 1.3 keV WDM produce optical depths to the Cosmic Microwave Background (CMB) that are inconsistent at 68% C.L. with current Planck results, even with extremely high ionizing radiation escape fractions, and 2.6 keV WDM requires an optimistic escape fraction to yield an optical depth consistent with Planck data. Although CMB optical depth calculations are model dependent, we find a strong challenge for stellar processes alone to reionize the universe in a 0.8 keV and 1.3 keV WDM cosmology.

Sterile neutrinos produced through a resonant Shi-Fuller mechanism are arguably the simplest model for a dark matter interpretation origin of the recent unidentified X-ray line seen toward a number of objects harboring dark matter. Here, I calculate the exact parameters required in this mechanism to produce the signal. The suppression of small scale structure predicted by these models is consistent with Local Group and high-$z$ galaxy count constraints. Very significantly, the parameters necessary in these models to produce the full dark matter density fulfill previously determined requirements to successfully match the Milky Way Galaxy's total satellite abundance, the satellites' radial distribution and their mass density profile, or "too big to fail problem." I also discuss how further precision determinations of the detailed properties of the candidate sterile neutrino dark matter can probe the nature of the quark-hadron transition, which takes place during the dark matter production.

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.

The American Physical Society's Division of Particles and Fields initiated a long-term planning exercise over 2012-13, with the goal of developing the community's long term aspirations. The sub-group "Dark Energy and CMB" prepared a series of papers explaining and highlighting the physics that will be studied with large galaxy surveys and cosmic microwave background experiments. This paper summarizes the findings of the other papers, all of which have been submitted jointly to the arXiv.

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.

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 study the spectra of gamma ray signals that arise from dark matter annihilation in the universe. We focus on the large class of theories where the photon spectrum includes both continuum spectrum of gamma rays that arise from annihilation into Standard Model states at tree level, as well as monochromatic gamma rays arising from annihilation directly into two photons at the one loop level. In this class of theories we obtain lower bounds on the ratio of the strength of the gamma ray line relative to the gamma ray continuum as a function of the dark matter mass and spin. These limits are obtained from the unitarity relation between the tree level amplitude of the primary annihilation channel and the imaginary part of the loop level amplitude for annihilation directly into photons, with the primary decay products running in the loop. These results are exact in the limit that dark matter annihilation is exclusively to a single Standard Model species, occurs through the lowest partial wave and respects CP. Away from this limit the bounds are approximate. Our conclusions agree with the known results in the literature in the case of the Minimal Supersymmetric Standard Model (MSSM). We use the Fermi-LAT observations to translate these limits into upper bounds on the dark matter annihilation cross section into any specific Standard Model state.

We study the relation between the measured anisotropies in the extragalactic diffuse gamma-ray background (DGRB) and the DGRB spectral intensity, and their potential origin from the unresolved blazar population. Using a physical-evolution model for blazars with a luminosity dependent evolution (LDDE) and an observationally-determined luminosity-dependent blazar spectral energy distribution (SED), we find that blazars can account for the observed anisotropy of the DGRB consistent with their observed source-count distribution, but are in turn constrained in contributing significantly to the observed DGRB intensity. For the best-fit LDDE model accounting for the DGRB anisotropy and source-count distribution, blazars only contribute 5.7^{+2.1}_{-1.0}% (68% CL) of the DGRB intensity above 1 GeV. Requiring a higher fraction of the DGRB intensity contribution by blazars overproduces the DGRB anisotropy, and therefore blazars in the LDDE+SED-sequence model cannot simultaneously account for the DGRB intensity as well as anisotropy. We discuss the limitations of LDDE models. However, these models do not require the many unjustified and observationally-inconsistent simplifying assumptions---including a single power law for all blazar spectra and a simple broken power-law model for their source-count distribution---that are present in much previous work.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

We show the conclusions claimed in the manuscript arXiv:1202.5309v1 by Cuoco, Komatsu and Siegal-Gaskins (CKS) are not generally valid. The results in CKS are based on a number of simplifying assumptions regarding the source population below the detection threshold and the threshold flux itself, and do not apply to many physical models of the blazar population. Physical blazar population models that match the measured source counts above the observational threshold can account for 60% of the diffuse gamma-ray background intensity between 1-10 GeV, while the assumptions in CKS limit the intensity to <30%. The shortcomings of the model considered in CKS arise from an over-simplified blazar source model. A number of the simplifying assumptions are unjustified, including: first, the adoption of an assumed power-law source-count distribution, dN/dS, to arbitrary low source fluxes, which is not exhibited in physical models of the blazar population; and, second, the lack of blazar spectral information in calculating the anisotropy of unresolved gamma-ray blazar emission. We also show that the calculation of the unresolved blazars' anisotropy is very sensitive to the spectral distribution of the unresolved blazars through the adopted source resolution threshold value, and must be taken into account in an accurate anisotropy calculation.

We perform a detailed examination of current constraints on annihilating and decaying dark matter models from both prompt and inverse-Compton emission photons, including both model-dependent and model-independent bounds. We also show that the observed isotropic diffuse gamma-ray background (DGRB), which provides one of the most conservative constraints on models of annihilating weak-scale dark matter particles, may enhance its sensitivity by a factor of ~2 to 3 (95% C.L.) as the Fermi-LAT experiment resolves DGRB contributing blazar sources with five years of observation. For our forecasts, we employ the results of constraints to the luminosity-dependent density evolution plus blazar spectral energy distribution sequence model, which is constrained by the DGRB and blazar source count distribution function.

We examine the constraints on models of weakly interacting massive particle (WIMP) dark matter from the recent observations of the Galactic center by the High Energy Spectroscopic System (HESS) telescope. We analyze canonical WIMP annihilation into Standard Model particle final states, including b/b-bar, t/t-bar and W^+ W^-. The constraints on annihilation into b/b-bar is within an order of magnitude of the thermal cross section at ~3 TeV, while the tau^+/tau^- channel is within a factor of ~2 of thermal. We also study constraints on Sommerfeld-enhanced dark matter annihilation models, and find that the gamma-ray observational constraints here rule out all of the parameter space consistent with dark matter annihilation interpretations of PAMELA and the Fermi-LAT e^+/e^- spectrum, in specific classes of models, and strongly constrains these interpretations in other classes. The gamma-ray constraints we find are more constraining on these models, in many cases, than current relic density, cosmic microwave background, halo shape and naturalness constraints.

We examine the constraints on the luminosity-dependent density evolution model for the evolution of blazars given the observed spectrum of the diffuse gamma-ray background (DGRB), blazar source-count distribution, and the blazar spectral energy distribution sequence model, which relates the observed the blazar spectrum to its luminosity. We show that the DGRB observed by the Large Area Telescope (LAT) aboard the Fermi Gamma Ray Space Telescope can be produced entirely by gamma-ray emission from blazars and nonblazar active galactic nuclei, and that our blazar evolution model is consistent with and constrained by the spectrum of the DGRB and flux source-count distribution function of blazars observed by Fermi-LAT. Our results are consistent with previous work that used EGRET spectral data to forecast the Fermi-LAT DGRB. The model includes only three free parameters, and forecasts that >~ 95% of the flux from blazars will be resolved into point sources by Fermi-LAT with 5 years of observation, with a corresponding reduction of the flux in the DGRB by a factor of ~2 to 3 (95% confidence level), which has implications for the Fermi-LAT's sensitivity to dark matter annihilation photons.

We present a comprehensive analysis of the halo model of cosmological large to small-scale structure statistics in the case of warm dark matter (WDM) structure formation scenarios. We include the effects of WDM on the linear matter power spectrum, halo density profile, halo concentration relation, halo mass function, subhalo density profile, subhalo mass function and biasing of the smooth dark matter component. As expected, we find large differences at the smallest physical scales in the nonlinear matter power spectrum predicted in the halo model between WDM and cold dark matter even for reasonably high-scale WDM particle masses. We find that significant effects are contributed from the alteration of the halo density profile and concentration, as well as the halo mass function. We further find that the effects of WDM on the subhalo population are important but sub-dominant. Clustering effects of the biasing of the smooth component in WDM is not largely significant.

Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach.

I show that the spectrum and morphology of a recent Fermi-LAT observation of the Galaxy center are consistent with a millisecond pulsar population in the nuclear Central stellar cluster of the Milky Way. The Galaxy Center gamma-ray spectrum is consistent with the spectrum of four of eight globular clusters that have been detected in the gamma-ray. A dark matter annihilation interpretation cannot be ruled out, though no unique features exist that would require this conclusion.

We examine the constraints on final state radiation from Weakly Interacting Massive Particle (WIMP) dark matter candidates annihilating into various standard model final states, as imposed by the measurement of the isotropic diffuse gamma-ray background by the Large Area Telescope aboard the Fermi Gamma-Ray Space Telescope. The expected isotropic diffuse signal from dark matter annihilation has contributions from the local Milky Way (MW) as well as from extragalactic dark matter. The signal from the MW is very insensitive to the adopted dark matter profile of the halos, and dominates the signal from extragalactic halos, which is sensitive to the low mass cut-off of the halo mass function. We adopt a conservative model for both the low halo mass survival cut-off and the substructure boost factor of the Galactic and extragalactic components, and only consider the primary final state radiation. This provides robust constraints which reach the thermal production cross-section for low mass WIMPs annihilating into hadronic modes. We also reanalyze limits from HESS observations of the Galactic Ridge region using a conservative model for the dark matter halo profile. When combined with the HESS constraint, the isotropic diffuse spectrum rules out all interpretations of the PAMELA positron excess based on dark matter annihilation into two lepton final states. Annihilation into four leptons through new intermediate states, although constrained by the data, is not excluded.

The WMAP haze is an excess in microwave emission coming from the center of the Milky Way galaxy. In the case of synchrotron emission models of the haze, we present tests for the source of radiating high-energy electrons/positrons. We explore several models in the case of a pulsar population or dark matter annihilation as the source. These morphological signatures of these models are small behind the WMAP Galactic mask, but are testable and constrain the source models. We show that detailed measurements of the morphology may distinguish between the pulsar and dark matter interpretations as well as differentiate among different pulsar models and dark matter profile models individually. Specifically, we find that a zero central density Galactic pulsar population model is in tension with the observed WMAP haze. The Planck Observatory's greater sensitivity and expected smaller Galactic mask should potentially provide a robust signature of the WMAP haze as either a pulsar population or the dark matter.

We propose an X-ray mission called Xenia to search for decaying superweakly interacting Dark Matter particles (super-WIMP) with a mass in the keV range. The mission and its observation plan are capable of providing a major break through in our understanding of the nature of Dark Matter (DM). It will confirm, or reject, predictions of a number of particle physics models by increasing the sensitivity of the search for decaying DM by about two orders of magnitude through a wide-field imaging X-ray spectrometer in combination with a dedicated observation program. The proposed mission will provide unique limits on the mixing angle and mass of neutral leptons, right handed partners of neutrinos, which are important Dark Matter candidates. The existence of these particles is strongly motivated by observed neutrino flavor oscillations and the problem of baryon asymmetry of the Universe. In super-WIMP models, the details of the formation of the cosmic web are different from those of LambdaCDM. The proposed mission will, in addition to the search for decaying Dark Matter, provide crucial insight into the nature of DM by studying the structure of the "cosmic web". This will be done by searching for missing baryons in emission, and by using gamma-ray bursts as backlight to observe the warm-hot intergalactic media in absorption.

This paper describes the Seventh Data Release of the Sloan Digital Sky Survey (SDSS), marking the completion of the original goals of the SDSS and the end of the phase known as SDSS-II. It includes 11663 deg^2 of imaging data, with most of the roughly 2000 deg^2 increment over the previous data release lying in regions of low Galactic latitude. The catalog contains five-band photometry for 357 million distinct objects. The survey also includes repeat photometry over 250 deg^2 along the Celestial Equator in the Southern Galactic Cap. A coaddition of these data goes roughly two magnitudes fainter than the main survey. The spectroscopy is now complete over a contiguous area of 7500 deg^2 in the Northern Galactic Cap, closing the gap that was present in previous data releases. There are over 1.6 million spectra in total, including 930,000 galaxies, 120,000 quasars, and 460,000 stars. The data release includes improved stellar photometry at low Galactic latitude. The astrometry has all been recalibrated with the second version of the USNO CCD Astrograph Catalog (UCAC-2), reducing the rms statistical errors at the bright end to 45 milli-arcseconds per coordinate. A systematic error in bright galaxy photometr is less severe than previously reported for the majority of galaxies. Finally, we describe a series of improvements to the spectroscopic reductions, including better flat-fielding and improved wavelength calibration at the blue end, better processing of objects with extremely strong narrow emission lines, and an improved determination of stellar metallicities. (Abridged)

There exists considerable interest in dark matter candidates that can reduce cosmological structure at sub-galactic scales through a suppression of the power spectrum of primordial perturbations as well as have a primordial velocity distribution to produce cores in the smallest-scale dwarf halos. Light keV-mass-scale sterile neutrinos can be such a "warm" dark matter candidate and are still viable candidates in broad regions of their parameter space. We review the status of this candidate and the current constraints from its radiative decay in the X-ray and from structure formation. We also provide a forecast for the sensitivity of the International X-ray Observatory to this decay, which we show may have the ability to detect or exclude sterile neutrinos within its full parameter space of interest as a dark matter candidate.

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 possibility of light sterile neutrinos allows for the resonant production of lepton number in the early universe through matter-affected neutrino mixing. For a given a mixing of the active and sterile neutrino states it has been found that the lepton number generation process is chaotic and strongly oscillatory. We undertake a new study of this process' sensitivity to initial conditions through the quantum rate equations. We confirm the chaoticity of the process in this solution, and moreover find that the resultant lepton number and the sign of the asymmetry produces a fractal in the parameter space of mass, mixing angle and initial baryon number. This has implications for future searches for sterile neutrinos, where arbitrary high sensitivity could not be determinate in forecasting the lepton number of the universe.

We measure the mass function of dark matter halos in a large set of collisionless cosmological simulations of flat LCDM cosmology and investigate its evolution at z<~2. Halos are identified as isolated density peaks, and their masses are measured within a series of radii enclosing specific overdensities. We argue that these spherical overdensity masses are more directly linked to cluster observables than masses measured using the friends-of-friends algorithm (FOF), and are therefore preferable for accurate forecasts of halo abundances. Our simulation set allows us to calibrate the mass function at z=0 for virial masses in the range 10^{11} Msol/h < M < 10^{15} Msol/h, to <~ 5%. We derive fitting functions for the halo mass function in this mass range for a wide range of overdensities, both at z=0 and earlier epochs. In addition to these formulae, which improve on previous approximations by 10-20%, our main finding is that the mass function cannot be represented by a universal fitting function at this level of accuracy. The amplitude of the "universal" function decreases monotonically by ~20-50%, depending on the mass definition, from z=0 to 2.5. We also find evidence for redshift evolution in the overall shape of the mass function.

We present upper limits on line emission in the Cosmic X-ray background (CXB) that would be produced by decay of sterile neutrino dark matter. We employ the spectra of the unresolved component of the CXB in the Chandra Deep Fields North and South obtained with the Chandra CCD detector in the E=0.8-9 keV band. The expected decay flux comes from the dark matter on the lines of sight through the Milky Way galactic halo. Our constraints on the sterile neutrino decay rate are sensitive to the modeling of the Milky Way halo. The highest halo mass estimates provide a limit on the sterile neutrino mass of m_s<2.9 keV in the Dodelson-Widrow production model, while the lowest halo mass estimates provide the conservative limit of m_s<5.7 keV (2-sigma). We also discuss constraints from a short observation of the softer (E<1 keV) X-ray background with a rocket-borne calorimeter by McCammon and collaborators.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

We study primordial nucleosynthesis abundance yields for assumed ranges of cosmological lepton numbers, sterile neutrino mass-squared differences and active-sterile vacuum mixing angles. We fix the baryon-to-photon ratio at the value derived from the cosmic microwave background (CMB) data and then calculate the deviation of the 2H, 4He, and 7Li abundance yields from those expected in the zero lepton number(s), no-new-neutrino-physics case. We conclude that high precision (< 5% error) measurements of the primordial 2H abundance from, e.g., QSO absorption line observations coupled with high precision (< 1% error) baryon density measurements from the CMB could have the power to either: (1) reveal or rule out the existence of a light sterile neutrino if the sign of the cosmological lepton number is known; or (2) place strong constraints on lepton numbers, sterile neutrino mixing properties and resonance sweep physics. Similar conclusions would hold if the primordial 4He abundance could be determined to better than 10%.

We identify galaxy groups and clusters in volume-limited samples of the SDSS redshift survey, using a redshift-space friends-of-friends algorithm. We optimize the friends-of-friends linking lengths to recover galaxy systems that occupy the same dark matter halos, using a set of mock catalogs created by populating halos of N-body simulations with galaxies. Extensive tests with these mock catalogs show that no combination of perpendicular and line-of-sight linking lengths is able to yield groups and clusters that simultaneously recover the true halo multiplicity function, projected size distribution, and velocity dispersion. We adopt a linking length combination that yields, for galaxy groups with ten or more members: a group multiplicity function that is unbiased with respect to the true halo multiplicity function; an unbiased median relation between the multiplicities of groups and their associated halos; a spurious group fraction of less than ~1%; a halo completeness of more than ~97%; the correct projected size distribution as a function of multiplicity; and a velocity dispersion distribution that is ~20% too low at all multiplicities. These results hold over a range of mock catalogs that use different input recipes of populating halos with galaxies. We apply our group-finding algorithm to the SDSS data and obtain three group and cluster catalogs for three volume-limited samples that cover 3495.1 square degrees on the sky. We correct for incompleteness caused by fiber collisions and survey edges, and obtain measurements of the group multiplicity function, with errors calculated from realistic mock catalogs. These multiplicity function measurements provide a key constraint on the relation between galaxy populations and dark matter halos.

We present a comprehensive analysis of constraints on the sterile neutrino as a dark matter candidate. The minimal production scenario with a standard thermal history and negligible cosmological lepton number is in conflict with conservative radiative decay constraints from the cosmic X-ray background in combination with stringent small-scale structure limits from the Lyman-alpha forest. We show that entropy release through massive particle decay after production does not alleviate these constraints. We further show that radiative decay constraints from local group dwarf galaxies are subject to large uncertainties in the dark matter density profile of these systems. Within the strongest set of constraints, resonant production of cold sterile neutrino dark matter in non-zero lepton number cosmologies remains allowed.

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.

I consider constraints from observations on a cutoff scale in clustering due to free streaming of the dark matter in a warm dark matter cosmological model with a cosmological constant. The limits are derived in the framework of a sterile neutrino warm dark matter universe, but can be applied to gravitinos and other models with small scale suppression in the linear matter power spectrum. With freedom in all cosmological parameters including the free streaming scale of the sterile neutrino dark matter, limits are derived using observations of the fluctuations in the cosmic microwave background, the 3D clustering of galaxies and 1D clustering of gas in the Lyman-alpha (Ly-alpha) forest in the Sloan Digital Sky Survey (SDSS), as well as the Ly-alpha forest in high-resolution spectroscopic observations. In the most conservative case, using only the SDSS main-galaxy 3D power-spectrum shape, the limit is m_s > 0.11 keV; including the SDSS Ly-alpha forest, this limit improves to m_s > 1.7 keV. More stringent constraints may be placed from the inferred matter power spectrum from high-resolution Ly-alpha forest observations, which has significant systematic uncertainties; in this case, the limit improves to m_s > 3.0 keV (all at 95% CL).

Sterile neutrinos, fermions with no standard model couplings [SU(2) singlets], are predicted by most extensions of the standard model, and may be the dark matter. I describe the nonthermal production and linear perturbation evolution in the early universe of this dark matter candidate. I calculate production of sterile neutrino dark matter including effects of Friedmann dynamics dictated by the quark-hadron transition and particle population, the alteration of finite temperature effective mass of active neutrinos due to the presence of thermal leptons, and heating of the coupled species due to the disappearance of degrees of freedom in the plasma. These effects leave the sterile neutrinos with a non-trivial momentum distribution. I also calculate the evolution of sterile neutrino density perturbations in the early universe through the linear regime and provide a fitting function form for the transfer function describing the suppression of small scale fluctuations for this warm dark matter candidate. The results presented here differ quantitatively from previous work due to the inclusion here of the relevant physical effects during the production epoch.

We study post weak decoupling coherent active-sterile and active-active matter-enhanced neutrino flavor transformation in the early universe. We show that flavor conversion efficiency at Mikheyev-Smirnov-Wolfenstein resonances is likely to be high (adiabatic evolution) for relevant neutrino parameters and energies. However, we point out that these resonances cannot sweep smoothly and continuously with the expansion of the universe. We show how neutrino flavor conversion in this way can leave both the active and sterile neutrinos with non-thermal energy spectra, and how, in turn, these distorted energy spectra can affect the neutron-to-proton ratio, primordial nucleosynthesis, and cosmological mass/closure constraints on sterile neutrinos. We demonstrate that the existence of a light sterile neutrino which mixes with active neutrinos can change fundamentally the relationship between the cosmological lepton numbers and the primordial nucleosynthesis He-4 yield.

The power spectrum is traditionally parameterized by a truncated Taylor series: $ln P(k) = ln P_* + (n_*-1) ln(k/k_*) + {1/2} n'_* ln^2(k/k_*)$. It is reasonable to truncate the Taylor series if $|n'_* ln(k/k_*)| << |n_*-1|$, but it is not if $|n'_* ln(k/k_*)| \gtrsim |n_*-1|$. We argue that there is no good theoretical reason to prefer $|n'_*| << |n_*-1|$, and show that current observations are consistent with $|n'_* ln(k/k_*)| ~ |n_*-1|$ even for $|ln(k/k_*)| ~ 1$. Thus, there are regions of parameter space, which are both theoretically and observationally relevant, for which the traditional truncated Taylor series parameterization is inconsistent, and hence it can lead to incorrect parameter estimations. Motivated by this, we propose a simple extension of the traditional parameterization, which uses no extra parameters, but that, unlike the traditional approach, covers well motivated inflationary spectra with $|n'_*| ~ |n_*-1|$. Our parameterization therefore covers not only standard-slow-roll inflation models but also a much wider class of inflation models. We use this parameterization to perform a likelihood analysis for the cosmological parameters.

The predicted mass function of dark matter halos is essential in connecting observed galaxy cluster counts and models of galaxy clustering to the properties of the primordial density field. We determine the mass function in the concordance $\Lambda$CDM cosmology, as well as its uncertainty, using sixteen $1024^3$-particle nested-volume dark-matter simulations, spanning a mass range of over five orders of magnitude. Using the nested volumes and single-halo tests, we find and correct for a systematic error in the friends-of-friends halo-finding algorithm. We find a fitting form and full error covariance for the mass function that successfully describes the simulations' mass function and is well-behaved outside the simulations' resolutions. Estimated forecasts of uncertainty in cosmological parameters from future cluster count surveys have negligible contribution from remaining statistical uncertainties in the central cosmology multiplicity function. There exists a potentially non-negligible cosmological dependence (non-universality) of the halo multiplicity function.

Measurements of the linear power spectrum of galaxies have placed tight constraints on neutrino masses. We extend the framework of the halo model of cosmological nonlinear matter clustering to include the effect of massive neutrino infall into cold dark matter (CDM) halos. The magnitude of the effect of neutrino clustering for three degenerate mass neutrinos with m_nu=0.9 eV is of order ~1%, within the potential sensitivity of upcoming weak lensing surveys. In order to use these measurements to further constrain--or eventually detect--neutrino masses, accurate theoretical predictions of the nonlinear power spectrum in the presence of massive neutrinos will be needed, likely only possible through high-resolution multiple particle (neutrino, CDM and baryon) simulations.

We use the projected correlation function w_p(r_p) of a volume-limited subsample of the Sloan Digital Sky Survey (SDSS) main galaxy redshift catalogue to measure the halo occupation distribution (HOD) of the galaxies of the sample. Simultaneously, we allow the cosmology to vary within cosmological constraints imposed by cosmic microwave background experiments in a Lambda-CDM model. We find that combining w_p(r_p) for this sample alone with the observations by WMAP, ACBAR, CBI and VSA can provide one of the most precise techniques available to measure cosmological parameters. For a minimal flat six-parameter Lambda-CDM model with an HOD with three free parameters, we find Omega_m=0.278^{+0.027}_{-0.027}, sigma_8=0.812^{+0.028}_{-0.027}, and H_0=69.8^{+2.6}_{-2.6}km s^{-1} Mpc^{-1}; these errors are significantly smaller than from CMB alone and similar to those obtained by combining CMB with the large-scale galaxy power spectrum assuming scale-independent bias. The corresponding HOD parameters describing the minimum halo mass and the normalization and cut-off of the satellite mean occupation are M_min=(3.03^{+0.36}_{-0.36})x 10^{12} h^{-1} M_sun, M_1 = (4.58^{+0.60}_{-0.60})x 10^{13} h^{-1} M_sun, and kappa=4.44^{+0.51}_{-0.69}. When more parameters are added to the HOD model, the error bars on the HOD parameters increase because of degeneracies, but the error bars on the cosmological parameters do not increase greatly. Similar modeling for other galaxy samples could reduce the statistical errors on these results, while more thorough investigations of the cosmology dependence of nonlinear halo bias and halo mass functions are needed to eliminate remaining systematic uncertainties, which may be comparable to statistical uncertainties.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

The emerging structure of the neutrino mass matrix, when combined with the primordial element abundances, places the most stringent constraint on the flavor asymmetries in the cosmological neutrino background and therefore its energy density. I review the mechanism of synchronized neutrino oscillations in the an early universe with degenerate (asymmetric) neutrino and antineutrino densities and the implications of refined measurements of neutrino parameters.

The Sloan Digital Sky Survey has validated and made publicly available its First Data Release. This consists of 2099 square degrees of five-band (u, g, r, i, z) imaging data, 186,240 spectra of galaxies, quasars, stars and calibrating blank sky patches selected over 1360 square degrees of this area, and tables of measured parameters from these data. The imaging data go to a depth of r ~ 22.6 and are photometrically and astrometrically calibrated to 2% rms and 100 milli-arcsec rms per coordinate, respectively. The spectra cover the range 3800--9200 A, with a resolution of 1800--2100. Further characteristics of the data are described, as are the data products themselves.

Weak gravitational lensing of background galaxies by intervening matter directly probes the mass distribution in the universe. This distribution, and its evolution at late times, is sensitive to both the dark energy, a negative pressure energy density component, and neutrino mass. We examine the potential of lensing experiments to measure features of both simultaneously. Focusing on the radial information contained in a future deep 4000 square degree survey, we find that the expected (1-sigma) error on a neutrino mass is 0.1 eV, if the dark energy parameters are allowed to vary. The constraints on dark energy parameters are similarly restrictive, with errors on w of 0.09. Much of the restrictive power on the dark energy comes not from the evolution of the gravitational potential but rather from how distances vary as a function of redshift in different cosmologies.

Recent models invoking extra spacelike dimensions inhabited by (bulk) neutrinos are shown to have significant cosmological effects if the size of the largest extra dimension is R > 1 fm. We consider effects on cosmic microwave background anisotropies, big bang nucleosynthesis, deuterium and Li-6 photoproduction, diffuse photon backgrounds, and structure formation. The resulting constraints can be stronger than either bulk graviton overproduction constraints or laboratory constraints.

New results, namely the independent determination of the deuterium abundance in several quasar absorption systems, and the complementary determination of the cosmological baryon density by observations of anisotropies in the cosmic microwave background (CMB), allow for a reevaluation of the constraints on the relativistic particle content of the universe at primordial nucleosynthesis. Expressed in terms of the neutrino energy density, we find 1.7 < N_\nu < 3.5 (95% CL). In particular, we show that phenomenological four neutrino models including a sterile state (not participating in SU(2)_L x U(1)_Y interactions) unavoidably thermalize a fourth neutrino, and are highly disfavored in the standard minimal model of primordial nucleosynthesis, if the systematic uncertainty in the primordial helium abundance is small. We describe plausible extensions of the minimal model which evade this constraint.

We point out that the relic densities of singlet (sterile) neutrinos of interest in viable warm and cold dark matter scenarios, depend on the characteristics of the QCD transition in the early universe. In the most promising of these dark matter scenarios the production of the singlets occurs at or near the QCD transition. Since production of the singlets, their dilution, and the disappearance of weak scatterers occur simultaneously, we calculate these processes contemporaneously to obtain accurate predictions of relic sterile neutrino mass density. Therefore, a determination of the mass and superweak mixing of the singlet neutrino through, for example, its radiative decay, along with a determination of its contribution to the critical density, can provide insight into the finite-temperature QCD transition.

We assess a mechanism which can transform neutrino-antineutrino asymmetries between flavors in the early universe, and confirm that such transformation is unavoidable in the near bi-maximal framework emerging for the neutrino mixing matrix. We show that the process is a standard Mikheyev-Smirnov-Wolfenstein flavor transformation dictated by a synchronization of momentum states. We also show that flavor ``equilibration'' is a special feature of maximal mixing, and carefully examine new constraints placed on neutrino asymmetries. In particular, the big bang nucleosynthesis limit on electron neutrino degeneracy xi_e < 0.04 does not apply directly to all flavors, yet confirmation of the large-mixing-angle solution to the solar neutrino problem will eliminate the possibility of degenerate big bang nucleosynthesis.

We report the discovery of a highly luminous optical transient (OT), SDSS J124602.54+011318.8, associated with a galaxy at a redshift of 0.385. In this paper we consider the possibility that the OT may be a GRB afterglow. Three sets of images and two sets of spectra were obtained as part of the normal operations of the Sloan Digital Sky Survey (SDSS). In the first two image sets, observed two nights apart, the object appears as a point source at $r^{*}\approx 17$. The third image set, observed about 410 days later, shows an extended source which is more than 2.5 magnitudes fainter. The spectra were observed about 400 and 670 days after the first two image sets, and both show an apparently normal galaxy at a redshift of 0.385. Associating the OT with the galaxy, the absolute magnitude was $M_{r^*}=-24.8$, which is over 4 magnitudes brighter than the most luminous supernova ever measured. The spectral energy distributions of the galaxy-subtracted OT derived from the first two image sets are well-fit by single power-laws with indices of $\beta_{\nu}=-0.92$ and -1.29 respectively, similar to most GRB afterglows. Based upon the luminosity of the OT, non-detections in contemporaneous ROTSE-I images, and the change in spectral slope, the OT, if an afterglow, was likely discovered early during a ``plateau'' or slowly-fading phase. The discovery of a GRB afterglow at this stage of the SDSS is consistent with expectations, but only if the optical emission is much less strongly beamed than the gamma-rays. We emphasize that other explanations for the OT cannot be ruled out; a recent follow-up study by [galyam02] provides strong evidence that this source is in fact an unusual AGN.

We point out a serendipitous link between warm dark matter (WDM) models for structure formation on the one hand and the high sensitivity energy range (1-10 keV) for x-ray photon detection on the Chandra and XMM-Newton observatories on the other. This fortuitous match may provide either a direct detection of the dark matter or exclusion of many candidates. We estimate expected x-ray fluxes from field galaxies and clusters of galaxies if the dark matter halos of these objects are composed of WDM candidate particles with rest masses in the structure formation-preferred range (~1 keV to ~20 keV) and with small radiative decay branches. Existing observations lead us to conclude that for singlet neutrinos (possessing a very small mixing with active neutrinos) to be a viable WDM candidate they must have rest masses < 5 keV in the zero lepton number production mode. Future deeper observations may detect or exclude the entire parameter range for the zero lepton number case, perhaps restricting the viability of singlet neutrino WDM models to those where singlet production is driven by a significant lepton number. The Constellation X project has the capability to detect/exclude singlet neutrino WDM for lepton number values up to 10% of the photon number. We also consider diffuse x-ray background constraints on these scenarios. These same x-ray observations additionally may constrain parameters of active neutrino and gravitino WDM candidates.

Dark energy models which alter the relative scaling behavior of dark energy and matter could provide a natural solution to the cosmic coincidence problem - why the densities of dark energy and dark matter are comparable today. A generalized class of dark energy models is introduced which allows non-canonical scaling of the ratio of dark matter and dark energy with the Robertson-Walker scale factor a(t). Upcoming observations, such as a high redshift supernova survey, application of the Alcock-Paczynski test to quasar pairs, and cluster evolution, will strongly constrain the relative scaling of dark matter and dark energy as well as the equation of state of the dark energy. Thus, whether there actually is a coincidence problem, and the extent of cosmic coincidence in the universe's recent past can be answered observationally in the near future. Determining whether today is a special time in the history of the universe will be a SNAP.

We calculate the incoherent resonant and non-resonant scattering production of sterile neutrinos in the early universe. We find ranges of sterile neutrino masses, vacuum mixing angles, and initial lepton numbers which allow these species to constitute viable hot, warm, and cold dark matter (HDM, WDM, CDM) candidates which meet observational constraints. The constraints considered here include energy loss in core collapse supernovae, energy density limits at big bang nucleosynthesis, and those stemming from sterile neutrino decay: limits from observed cosmic microwave background anisotropies, diffuse extragalactic background radiation, and Li-6/D overproduction. Our calculations explicitly include matter effects, both effective mixing angle suppression and enhancement (MSW resonance), as well as quantum damping. We for the first time properly include all finite temperature effects, dilution resulting from the annihilation or disappearance of relativistic degrees of freedom, and the scattering-rate-enhancing effects of particle-antiparticle pairs (muons, tauons, quarks) at high temperature in the early universe.

We demonstrate a qualitatively new aspect of the dynamics of Gamma-Ray Burst (GRB) fireballs: the development of a substantial dispersion in the proton component in fireballs in which neutron decoupling occurs and is sufficiently pronounced. This effect depends sensitively on the neutron to proton ratio in the fireball, becoming more dramatic with increasing neutron excess. Simple physical arguments and transport calculations indicate that the dispersion in Lorentz factor of the protons can be of order the final mean Lorentz factor of the fireball. We show how plasma instabilities could play an important role in the evolution of the fireball and how they might ultimately govern the development of such a velocity dispersion in the proton component. The role of these instabilities in setting/diminishing a proton Lorentz factor dispersion represents a new and potentially important venue for the study of plasma instabilities. Significant dispersion in the proton velocities translates into fewer protons attaining the highest Lorentz factors. This is tantamount to a reduction in the total energy required to attain a given Lorentz factor for the highest energy protons. As well, a proton component dispersion can have consequences for the electromagnetic and neutrino signature of GRBs.

We point out that the baryon-loading problem in Gamma-Ray Burst (GRB) models can be amelioriated if a significant fraction of the baryons which inertially confine the fireball are converted to neutrons. A high neutron fraction in some circumstances can result in a reduced transfer of energy from relativistic light particles in the fireball to baryons. The energy needed to produce the required relativistic flow in the GRB is consequently reduced, in some cases by orders of magnitude. This could be relevant to GRB models because a high neutron-to-proton ratio has been calculated in neutron star-merger fireball environments. Significant neutron excess also could occur near compact objects with high neutrino fluxes.

We assess the effects on Big Bang Nucleosynthesis (BBN) of lepton number generation in the early universe resulting from the two-doublet four-neutrino mass/mixing scheme. It has been argued that this neutrino mass/mixing arrangement gives the most viable fit to the existing data. We study full 4 x 4 mixing matrices and show how possible symmetries in these can affect the BBN He-4 abundance yields. Though there is as yet no consensus on the reliability of BBN calculations with neutrino flavor mixing, we show that, in the case where the sign of the lepton number asymmetry is unpredictable, BBN considerations may pick out specific relationships between mixing angles. In particular, reconciling the observed light element abundances with a \bar\nu_\mu <-> \bar\nu_e oscillation interpretation of LSND would allow unique new constraints on the neutrino mixing angles in this model.

We examine in detail the changes in the production of primordial helium resulting from nonthermal neutrino momentum distributions produced by resonant transformation of electron-type neutrinos to steriles. These transformations, anti-nu_e -> anti-nu_s (nu_e->nu_s), amplify a positive (negative) lepton number asymmetry. We find that the resulting supression relative to a thermal distribution of low energy nu_e reduces n -> p conversion to a greater extent than does the enhancement of n -> p from an identical suppression of anti-nu_e. Thus, equal lepton-number asymmetries of opposite sign have unequal effects on the resulting helium yield in primordial nucleosynthesis.

This is a reply to the preprint "On Neutrino-Mixing-Generated Lepton Asymmetry and the Primordial Helium-4 Abundance" by M. V. Chizhov and D. P. Kirilova (hep-ph/9908525), which criticised our recent publication (X. Shi, G. M. Fuller and K. Abazajian Phys. Rev. D 60, 063002 (1999)). Here we point out factual errors in their description of what our paper says. We also show that their criticisms of our work have no merit.

It has been proposed that an asymmetry in the electron neutrino sector may be generated by resonant active-sterile neutrino transformations during Big Bang Nucleosynthesis (BBN). We calculate the change in the primordial $^4$He yield $Y$ resulting from this asymmetry, taking into account both the time evolution of the $\nu_e$ and $\bar\nu_e$ distribution function and the spectral distortions in these. We calculate this change in two schemes: (1) a lepton asymmetry directly generated by $\nu_e$ mixing with a lighter right-handed sterile neutrino $\nu_s$; and (2) a lepton asymmetry generated by a $\nu_\tau\leftrightarrow\nu_s$ or $\nu_\mu\leftrightarrow\nu_s$ transformation which is subsequently partially converted to an asymmetry in the $\nu_e\bar\nu_e$ sector by a matter-enhanced active-active neutrino transformation. In the first scheme, we find that the percentage change in $Y$ is between -1% and 9% (with the sign depending on the sign of the asymmetry), bounded by the Majorana mass limit $m_{\nu_e}\la 1$ eV. In the second scheme, the maximal percentage reduction in $Y$ is 2%, if the lepton number asymmetry in neutrinos is positive; Otherwise, the percentage increase in $Y$ is $\la 5%$ for $m^2_{\nu_\mu,\nu_\tau}-m^2_{\nu_s}\la 10^4$ eV. We conclude that the change in the primordial $^4$He yield induced by a neutrino-mixing-generated lepton number asymmetry can be substantial in the upward direction, but limited in the downward direction.

We investigate the effects of matter-enhanced (MSW) transformation of neutrinos in the Early Universe on the primordial helium yield (Y_p). We find that Y_p is affected much more by the MSW-induced alterations in the neutrino energy spectra than by the associated change in expansion rate. Specifically, the absence due to transformation of low energy electron neutrinos can significantly affect neutron-proton weak interconversion rates through the lifting of Fermi-blocking of neutron decay at low energies and through halting low-energy neutrino capture on neutrons. We find that the change of Y_p within a causal horizon is -0.005\le \delta Y_p \leq 0.013 for m_{\nu_{\mu,\tau}}^2 - m_{\nu_s}^2 \le 10^4 eV^2 in the case of \nu_{\mu,\tau}-\nu_s-\nu_e mixing, with the lower limit at m_{\nu_{\mu,\tau}}^2 - m_{\nu_s}^2 \approx 100 eV^2, and -0.002\le \delta Y_p \leq 0.020 for m_{\nu_e}^2 - m_{\nu_s}^2 \le 1 eV^2 in the case of \nu_e-\nu_s mixing.

The formation of supermassive black holes through the gravitational collapse of supermassive objects ($M \ga 10^4 M_\odot$) has been proposed as a source of cosmological gamma-ray bursts. The major advantage of this model is that such collapses are far more energetic than stellar-remnant mergers. The major drawback of this idea is the severe baryon loading problem in one-dimensional models. We can show that the observed log N - log P (number vs. peak flux) distribution for gamma-ray bursts in the BATSE database is not inconsistent with an identification of supermassive object collapse as the origin of the gamma-ray bursts. This conclusion is valid for a range of plausible cosmological and gamma-ray burst spectral parameters.

We investigate the effects of baryon/anti-baryon inhomogeneity on primordial nucleosynthesis. Recent work claims that electroweak baryogenesis could give rise to distinct regions of net baryon and anti-baryon number, which could survive until the nucleosynthesis epoch. We discuss neutron diffusion effects on nucleosynthesis yields in these models and speculate on the prospects for obtaining constraints.