We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of $M_\mathrm{halo} \approx 10^{10}\,$M$_\odot$ at $z=0$, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r$_{1/2}<500$ pc) have lower $\sigma_\star/V_\mathrm{max}$ ratios, reinforcing the idea that smaller dwarfs may reside in halos that are more massive than is naively expected. The majority of dwarfs simulated with self-interactions actually experience contraction of their inner density profiles with the addition of baryons relative to the cores produced in dark-matter-only runs, though the simulated dwarfs are always less centrally dense than in $\Lambda$CDM. The V$_{1/2}-$r$_{1/2}$ relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Spatially-resolved rotation curves in the central regions ($<400$ pc) of small dwarfs could provide a way to distinguish between CDM, WDM, and SIDM, however: at the masses probed in this simulation suite, cored density profiles in dwarfs with small r$_{1/2}$ values can only originate from dark matter self-interactions.

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.

We simulate the tidal disruption of a collisionless N-body globular star cluster in a total of 300 different orbits selected to have galactocentric radii between 10 and 30 kpc in four dark matter halos: (a) a spherical halo with no subhalos, (b) a spherical halo with subhalos, (c) a realistic halo with no subhalos, and (d) a realistic halo with subhalos. This allows us to isolate and study how the halo's (lack of) dynamical symmetry and substructures affect the dispersal of tidal debris. The realistic halos are constructed from the snapshot of the Via Lactea II simulation at redshift zero. We find that the overall halo's lack of dynamical symmetry disperses tidal debris to make the streams fluffier, consistent with previous studies of tidal debris of dwarf galaxies in larger orbits than ours in this study. On the other hand, subhalos in realistic potentials can locally enhance the densities along streams, making streams denser than their counterparts in smooth potentials. We show that many long and thin streams can survive in a realistic and lumpy halo for a Hubble time. This suggests that upcoming stellar surveys will likely uncover more thin streams which may contain density gaps that have been shown to be promising probes for dark matter substructures.

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.

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.

If the dark matter (DM) were composed of axions, then structure formation in the Universe would be suppressed below the axion Jeans scale. Using an analytic model for the halo mass function of a mixed DM model with axions and cold dark matter, combined with the abundance-matching technique, we construct the UV-luminosity function. Axions suppress high-$z$ galaxy formation and the UV-luminosity function is truncated at a faintest limiting magnitude. From the UV-luminosity function, we predict the reionization history of the universe and find that axion DM causes reionization to occur at lower redshift. We search for evidence of axions using the Hubble Ultra Deep Field UV-luminosity function in the redshift range $z=6$-$10$, and the optical depth to reionization, $\tau$, as measured from cosmic microwave background polarization. All probes we consider consistently exclude $m_a\lesssim 10^{-23}\text{ eV}$ from contributing more than half of the DM, with our strongest constraint ruling this model out at more than $8\sigma$ significance. In conservative models of reionization a dominant component of DM with $m_a=10^{-22}\text{ eV}$ is in $3\sigma$ tension with the measured value of $\tau$, putting pressure on an axion solution to the cusp-core problem. Tension is reduced to $2\sigma$ for the axion contributing only half of the DM. A future measurement of the UV-luminosity function in the range $z=10$-$13$ by JWST would provide further evidence for or against $m_a=10^{-22}\text{ eV}$. Probing still higher masses of $m_a=10^{-21}\text{ eV}$ will be possible using future measurements of the kinetic Sunyaev-Zel'dovich effect by Advanced ACTPol to constrain the time and duration of reionization.

We simulate tidal streams in the presence and absence of substructures inside the zero redshift snapshot of the Via Lactea II (VL-2) simulation. A halo finder is used to remove and isolate the subhalos found inside the high resolution dark matter halo of VL-2, and the potentials for both the main halo and all the subhalos are constructed individually using the self-consistent field (SCF) method. This allows us to make direct comparison of tidal streams between a smooth halo and a lumpy halo without assuming idealized profiles or triaxial fits. We simulate the kinematics of a star cluster starting with the same orbital position but two different velocities. Although these two orbits are only moderately eccentric and have similar apo- and pericentric distances, we find that the two streams have very different morphologies. We conclude that our model of the potential of VL-2 can provide insights about tidal streams that have not been explored by previous studies using idealized or axisymmetric models.

Observations of diffuse Galactic gamma ray emission (DGE) by the Fermi Large Area Telescope (LAT) allow a detailed study of cosmic rays and the interstellar medium. However, diffuse emission models of the inner Galaxy underpredict the Fermi-LAT data at energies above a few GeV and hint at possible non-astrophysical sources including dark matter (DM) annihilations or decays. We present a study of the possible emission components from DM using the high-resolution Via Lactea II N-body simulation of a Milky Way-sized DM halo. We generate full-sky maps of DM annihilation and decay signals that include modeling of the adiabatic contraction of the host density profile, Sommerfeld enhanced DM annihilations, $p$-wave annihilations, and decaying DM. We compare our results with the DGE models produced by the Fermi-LAT team over different sky regions, including the Galactic center, high Galactic latitudes, and the Galactic anti-center. This work provides possible templates to fit the observational data that includes the contribution of the subhalo population to DM gamma-ray emission, with the significance depending on the annihilation/decay channels and the Galactic regions being considered.

The observed population of the Milky Way satellite galaxies offer a unique testing ground for galaxy formation theory on small-scales. Our novel approach was to investigate the clustering of the known Milky Way satellite galaxies and to quantify the amount of substructure within their distribution using a two-point correlation function statistic in each of three spaces: configuration space, line-of-sight velocity space, and four-dimensional phase-space. These results were compared to those for three sets of subhaloes in the Via Lactea II Cold Dark Matter simulation defined to represent the luminous dwarfs. We found no evidence at a significance level above 2-sigma of substructure within the distribution of the Milky Way satellite galaxies in any of the three spaces. The "luminous" subhalo sets are more strongly clustered than are the Milky Way satellites in all three spaces and over a broader range of scales in four-dimensional phase-space. Each of the "luminous" subhalo sets are clustered as a result of substructure within their line-of-sight velocity space distributions at greater than 3-sigma significance, whereas the Milky Way satellite galaxies are randomly distributed in line-of-sight velocity space. While our comparison is with only one Cold Dark Matter simulation, the inconsistencies between the Milky Way satellite galaxies and the Via Lactea II subhalo sets for all clustering methods suggest a potential new 'small-scale' tension between Cold Dark Matter theory and the observed Milky Way satellites. Future work will obtain a more robust comparison between the observed Milky Way satellites and Cold Dark Matter theory by studying additional simulations.

We evaluate the ability of future data sets to discriminate among different quintessence dark energy models. This approach gives an alternative measure for assessing the impact of future experiments, as compared with the large body of literature that compares experiments in abstract parameter spaces and more recent work that evaluates the constraining power of experiments on individual parameter spaces of specific quintessence models. We use the Dark Energy Task Force (DETF) models of future data sets, and compare the discriminative power of experiments designated by the DETF as Stages 2, 3, and 4. Our work reveals a minimal increase in discriminating power when comparing Stage 3 to Stage 2, but a very striking increase in discriminating power when going to Stage 4. We also see evidence that even modest improvements over DETF Stage 4 could result in even more dramatic discriminating power among quintessence dark energy models. We develop and demonstrate the technique of using the independently measured modes of the equation of state as a common parameter space in which to compare the different quintessence models, and we argue that this technique is a powerful one. We use the PNGB, Exponential, Albrecht-Skordis, and Inverse Tracker (or Inverse Power Law) quintessence models for this work. One of our main results is that the goal of discriminating among these models sets a concrete measure on the capabilities of future dark energy experiments. Experiments have to be somewhat better than DETF Stage 4 simulated experiments to fully meet this goal.

Current cosmological observations indicate a preference for a cosmological constant that is drastically smaller than what can be explained by conventional particle physics. The Causal Entropic Principle (Bousso, {\it et al}.) provides an alternative approach to anthropic attempts to predict our observed value of the cosmological constant by calculating the entropy created within a causal diamond. We have extended this work to use the Causal Entropic Principle to predict the preferred curvature within the "multiverse". We have found that values larger than $\rho_k = 40\rho_m$ are disfavored by more than 99.99% and a peak value at $\rho_{\Lambda} = 7.9 \times 10^{-123}$ and $\rho_k =4.3 \rho_m$ for open universes. For universes that allow only positive curvature or both positive and negative curvature, we find a correlation between curvature and dark energy that leads to an extended region of preferred values. Our universe is found to be disfavored to an extent depending the priors on curvature. We also provide a comparison to previous anthropic constraints on open universes and discuss future directions for this work.

We report on the results of a Markov Chain Monte Carlo (MCMC) analysis of an inverse power law (IPL) quintessence model using the Dark Energy Task Force (DETF) simulated data sets as a representation of future dark energy experiments. We generate simulated data sets for a Lambda-CDM background cosmology as well as a case where the dark energy is provided by a specific IPL fiducial model and present our results in the form of likelihood contours generated by these two background cosmologies. We find that the relative constraining power of the various DETF data sets on the IPL model parameters is broadly equivalent to the DETF results for the w_{0}-w_{a} parameterization of dark energy. Finally, we gauge the power of DETF "Stage 4" data by demonstrating a specific IPL model which, if realized in the universe, would allow Stage 4 data to exclude a cosmological constant at better than the 3-sigma level.

We consider the effect of future dark energy experiments on ``Albrecht-Skordis'' (AS) models of scalar field dark energy using the Monte-Carlo Markov chain method. We deal with the issues of parameterization of these models, and have included spatial curvature as a parameter, finding it to be important. We use the Dark Energy Task Force (DETF) simulated data to represent future experiments and report our results in the form of likelihood contours in the chosen parameter space. Simulated data is produced for cases where the background cosmology has a cosmological constant, as well as cases where the dark energy is provided by the AS model. The latter helps us demonstrate the power of DETF Stage 4 data in the context of this specific model. Though the AS model can produce equations of state functions very different from what is possible with the $w_0-w_a$ parametrization used by the DETF, our results are consistent with those reported by the DETF.

We present an analysis of a scalar field model of dark energy with an exponential potential using the Dark Energy Task Force (DETF) simulated data models. Using Markov Chain Monte Carlo sampling techniques we examine the ability of each simulated data set to constrain the parameter space of the exponential potential for data sets based on a cosmological constant and a specific exponential scalar field model. We compare our results with the constraining power calculated by the DETF using their ``$w_0-w_a$'' parametrization of the dark energy. We find that respective increases in constraining power from one stage to the next produced by our analysis give results consistent with DETF results. To further investigate the potential impact of future experiments, we also generate simulated data for an exponential model background cosmology which can not be distinguished from a cosmological constant at DETF ``Stage 2'', and show that for this cosmology good DETF Stage 4 data would exclude a cosmological constant by better than 3$\sigma$.

We analyze the constraining power of future dark energy experiments for Pseudo-Nambu Goldstone Boson (PNGB) quintessence. Following the Dark Energy Task Force methodology, we forecast data for three experimental ``stages'': Stage 2 represents in-progress projects relevant to dark energy; Stage 3 refers to medium sized experiments; Stage 4 comprises larger projects. We determine the posterior probability distribution for the parameters of the PNGB model using Markov Chain Monte Carlo analysis. Utilizing data generated on a $\Lambda CDM$ cosmology, we find that the relative power of the different data stages on PNGB quintessence is roughly comparable to the DETF results for the $w_0-w_a$ parametrization of dark energy. We also generate data based on a PNGB cosmological model that is consistent with a $\Lambda CDM$ fiducial model at Stage 2. We find that Stage 4 data based on this PNGB fiducial model will rule out a cosmological constant by at least $3 \sigma$.