We present N‐body simulations of a new class of self‐interacting dark matter models, which do not violate any astrophysical constraints due to a non‐power‐law velocity dependence of the transfer cross‐section which is motivated by a Yukawa‐like new gauge boson interaction. Specifically, we focus on the formation of a Milky‐Way‐like dark matter halo taken from the Aquarius project and resimulate it for a couple of representative cases in the allowed parameter space of this new model. We find that for these cases, the main halo only develops a small core (∼1 kpc) followed by a density profile identical to that of the standard cold dark matter scenario outside of that radius. Neither the subhalo mass function nor the radial number density of subhaloes is altered in these models but there is a significant change in the inner density structure of subhaloes resulting in the formation of a large density core. As a consequence, the inner circular velocity profiles of the most massive subhaloes differ significantly from the cold dark matter predictions and we demonstrate that they are compatible with the observational data of the brightest Milky Way dwarf spheroidals (dSphs) in such a velocity‐dependent self‐interacting dark matter scenario. Specifically, and contrary to the cold dark matter case, there are no subhaloes that are more concentrated than what is inferred from the kinematics of the Milky Way dSphs. We conclude that these models offer an interesting alternative to the cold dark matter model that can reduce the recently reported tension between the brightest Milky Way satellites and the dense subhaloes found in cold dark matter simulations.

Self-Interacting Dark Matter is an attractive alternative to the Cold Dark Matter paradigm only if it is able to substantially reduce the central densities of dwarf-size haloes while keeping the densities and shapes of cluster-size haloes within current constraints. Given the seemingly stringent nature of the latter, it was thought for nearly a decade that Self-Interacting Dark Matter would be viable only if the cross section for self-scattering was strongly velocitydependent. However, it has recently been suggested that a constant cross section per unit mass of σ T /m ∼ 0.1 cm 2 g −1 is sufficient to accomplish the desired effect. We explicitly investigate this claim using high resolution cosmological simulations of a Milky-Way size halo and find that, similarly to the Cold Dark Matter case, such cross section produces a population of massive subhaloes that is inconsistent with the kinematics of the classical dwarf spheroidals, in particular with the inferred slopes of the mass profiles of Fornax and Sculptor. This problem is resolved if σ T /m ∼ 1 cm 2 g −1 at the dwarf spheroidal scales. Since this value is likely inconsistent with the halo shapes of several clusters, our results leave only a small window open for a velocity-independent Self-Interacting Dark Matter model to work as a distinct alternative to Cold Dark Matter.

We present a theoretical analysis of some unexplored aspects of relaxed Bose-Einstein condensate dark matter (BECDM) haloes. This type of ultralight bosonic scalar field dark matter is a viable alternative to the standard cold dark matter (CDM) paradigm, as it makes the same large-scale predictions as CDM and potentially overcomes CDM's small-scale problems via a galaxy-scale de Broglie wavelength. We simulate BECDM halo formation through mergers, evolved under the Schrödinger-Poisson equations. The formed haloes consist of a soliton core supported against gravitational collapse by the quantum pressure tensor and an asymptotic r −3 NFW-like profile. We find a fundamental relation of the core=to-halo mass with the dimensionless invariant Ξ ≡ |E|/M 3 /(Gm/ ) 2 or M c /M 2.6Ξ 1/3 , linking the soliton to global halo properties. For r ≥ 3.5 r c core radii, we find equipartition between potential, classical kinetic, and quantum gradient energies. The haloes also exhibit a conspicuous turbulent behavior driven by the continuous reconnection of vortex lines due to wave interference. We analyse the turbulence 1D velocity power spectrum and find a k −1.1 power-law. This suggests the vorticity in BECDM haloes is homogeneous, similar to thermally-driven counterflow BEC systems from condensed matter physics, in contrast to a k −5/3 Kolmogorov power-law seen in mechanically-driven quantum systems. The mode where the power spectrum peaks is approximately the soliton width, implying the soliton-sized granules carry most of the turbulent energy in BECDM haloes.

We present the first simulations within an effective theory of structure formation (ETHOS), which includes the effect of interactions between dark matter and dark radiation on the linear initial power spectrum and dark matter self-interactions during non-linear structure formation. We simulate a Milky Way-like halo in four different dark matter models and the cold dark matter case. Our highest resolution simulation has a particle mass of 2.8 × 10 4 M and a softening length of 72.4 pc. We demonstrate that all alternative models have only a negligible impact on large scale structure formation. On galactic scales, however, the models significantly affect the structure and abundance of subhaloes due to the combined effects of small scale primordial damping in the power spectrum and late time self-interactions. We derive an analytic mapping from the primordial damping scale in the power spectrum to the cutoff scale in the halo mass function and the kinetic decoupling temperature. We demonstrate that certain models within this extended effective framework that can alleviate the too-big-to-fail and missing satellite problems simultaneously, and possibly the core-cusp problem. The primordial power spectrum cutoff of our models naturally creates a diversity in the circular velocity profiles, which is larger than that found for cold dark matter simulations. We show that the parameter space of models can be constrained by contrasting model predictions to astrophysical observations. For example, some models may be challenged by the missing satellite problem if baryonic processes were to be included and even over-solve the too-big-to-fail problem; thus ruling them out.

We formulate an effective theory of structure formation (ETHOS) that enables cosmological structure formation to be computed in almost any microphysical model of dark matter physics. This framework maps the detailed microphysical theories of particle dark matter interactions into the physical effective parameters that shape the linear matter power spectrum and the self-interaction transfer cross section of non-relativistic dark matter. These are the input to structure formation simulations, which follow the evolution of the cosmological and galactic dark matter distributions. Models with similar effective parameters in ETHOS but with different dark particle physics would nevertheless result in similar dark matter distributions. We present a general method to map an ultraviolet complete or effective field theory of low energy dark matter physics into parameters that affect the linear matter power spectrum and carry out this mapping for several representative particle models. We further propose a simple but useful choice for characterizing the dark matter self-interaction transfer cross section that parametrizes self-scattering in structure formation simulations. Taken together, these effective parameters in ETHOS allow the classification of dark matter theories according to their structure formation properties rather than their intrinsic particle properties, paving the way for future simulations to span the space of viable dark matter physics relevant for structure formation.

Using constrained simulations of the local Universe for generic cold dark matter and for 1 keV warm dark matter, we investigate the difference in the abundance of dark matter halos in the local environment. We find that the mass function within 20 h −1 Mpc of the Local Group is ∼ 2 times larger than the universal mass function in the 10 9 − 10 13 h −1 M ⊙ mass range. Imposing the field of view of the on-going HI blind survey ALFALFA in our simulations, we predict that the velocity function in the Virgo-direction region exceeds the universal velocity function by a factor of 3. Furthermore, employing a scheme to translate the halo velocity function into a galaxy velocity function, we compare the simulation results with a sample of galaxies from the early catalog release of ALFALFA. We find that our simulations are able to reproduce the velocity function in the 80 − 300 km s −1 velocity range, having a value ∼ 10 times larger than the universal velocity function in the Virgo-direction region. In the low velocity regime, 35 − 80 km s −1 , the warm dark matter simulation reproduces the observed flattening of the velocity function. On the contrary, the simulation with cold dark matter predicts a steep rise in the velocity function towards lower velocities; for V max = 35 km s −1 , it forecasts ∼ 10 times more sources than the ones observed. If confirmed by the complete ALFALFA survey, our results indicate a potential problem for the cold dark matter paradigm or for the conventional assumptions about energetic feedback in dwarf galaxies.

Dark matter annihilation or decay could have a significant impact on the ionization and thermal history of the universe. In this paper, we study the potential contribution of dark matter annihilation (s-wave-or p-wave-dominated) or decay to cosmic reionization, via the production of electrons, positrons and photons. We map out the possible perturbations to the ionization and thermal histories of the universe due to dark matter processes, over a broad range of velocity-averaged annihilation cross-sections/decay lifetimes and dark matter masses. We have employed recent numerical studies of the efficiency with which annihilation/decay products induce heating and ionization in the intergalactic medium, and in this work extended them down to a redshift of 1 + z = 4 for two different reionization scenarios. We also improve on earlier studies by using the results of detailed structure formation models of dark matter haloes and subhaloes that are consistent with up-to-date N -body simulations, with estimates on the uncertainties that originate from the smallest scales. We find that for dark matter models that are consistent with experimental constraints, a contribution of more than 10% to the ionization fraction at reionization is disallowed for all annihilation scenarios. Such a contribution is possible only for decays into electron/positron pairs, for light dark matter with mass mχ 100 MeV, and a decay lifetime τχ ∼ 10 24 − 10 25 s.

We present the first cosmological simulations of dwarf galaxies, which include dark matter self-interactions and baryons. We study two dwarf galaxies within cold dark matter, and four different elastic self-interacting scenarios with constant and velocity-dependent cross sections, motivated by a new force in the hidden dark matter sector. Our highest resolution simulation has a baryonic mass resolution of 1.8 × 10 2 M and a gravitational softening length of 34 pc at z = 0. In this first study we focus on the regime of mostly isolated dwarf galaxies with halo masses ∼ 10 10 M where dark matter dynamically dominates even at sub-kpc scales. We find that while the global properties of galaxies of this scale are minimally affected by allowed self-interactions, their internal structures change significantly if the cross section is large enough within the inner sub-kpc region. In these dark-matter-dominated systems, selfscattering ties the shape of the stellar distribution to that of the dark matter distribution. In particular, we find that the stellar core radius is closely related to the dark matter core radius generated by self-interactions. Dark matter collisions lead to dwarf galaxies with larger stellar cores and smaller stellar central densities compared to the cold dark matter case. The central metallicity within 1 kpc is also larger by up to ∼ 15% in the former case. We conclude that the mass distribution, and characteristics of the central stars in dwarf galaxies can potentially be used to probe the self-interacting nature of dark matter.