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 vel
ocity-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.8times 10^2,{rm M}_odot$ and a gravitational softening length of $34,{rm pc}$ at $z=0$. In this first study we focus on the regime of mostly isolated dwarf galaxies with halo masses $sim10^{10},{rm M}_odot$ 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, self-scattering 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,{rm kpc}$ is also larger by up to $sim 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.
We demonstrate that self-interacting dark matter models with interactions mediated by light particles can have significant deviations in the matter power-spectrum and detailed structure of galactic halos when compared to a standard cold dark matter s
cenario. While these deviations can take the form of suppression of small scale structure that are in some ways similar to that of warm dark matter, the self-interacting models have a much wider range of possible phenomenology. A long-range force in the dark matter can introduce multiple scales to the initial power spectrum, in the form of dark acoustic oscillations and an exponential cut-off in the power spectrum. Using simulations we show that the impact of these scales can remain observationally relevant up to the present day. Furthermore, the self-interaction can continue to modify the small-scale structure of the dark matter halos, reducing their central densities and creating a dark matter core. The resulting phenomenology is unique to this type of models.
The IceCube Neutrino Observatory has observed highly energetic neutrinos in excess of the expected atmospheric neutrino background. It is intriguing to consider the possibility that such events are probing physics beyond the standard model. In this c
ontext, $mathcal{O}$(PeV) dark matter particles decaying to neutrinos have been considered while dark matter annihilation has been dismissed invoking the unitarity bound as a limiting factor. However, the latter claim was done ignoring the contribution from dark matter substructure, which for PeV Cold Dark Matter would extend down to a free streaming mass of $mathcal{O}$($10^{-18}$M$_odot$). Since the unitarity bound is less stringent at low velocities, ($sigma_{rm ann}$v)$leq4pi/m_chi^2v$, then, it is possible that these cold and dense subhalos would contribute dominantly to a dark-matter-induced neutrino flux and easily account for the events observed by IceCube. A Sommerfeld-enhanced dark matter model can naturally support such scenario. Interestingly, the spatial distribution of the events shows features that would be expected in a dark matter interpretation. Although not conclusive, 9 of the 37 events appear to be clustered around a region near the Galactic Center while 6 others spatially coincide, within the reported angular errors, with 5 of 26 Milky Way satellites. However, a simple estimate of the probability of the latter occurring by chance is $sim35%$. More events are needed to statistically test this hypothesis. PeV dark matter particles are massive enough that their abundance as standard thermal relics would overclose the Universe. This issue can be solved in alternative scenarios, for instance if the decay of new massive unstable particles generates significant entropy reheating the Universe to a slightly lower temperature than the freeze-out temperature, $T_{rm RH} lesssim T_{rm f}sim4times10^4$~GeV.
We present a model for the structure of the particle phase space average density ($P^2SAD$) in galactic haloes, introduced recently as a novel measure of the clustering of dark matter. Our model is based on the stable clustering hypothesis in phase s
pace, the spherical collapse model, and tidal disruption of substructures, which is calibrated against the Aquarius simulations. Using this model, we can predict the behaviour of $P^2SAD$ in the numerically unresolved regime, down to the decoupling mass limit of generic WIMP models. This prediction can be used to estimate signals sensitive to the small scale structure of dark matter. For example, the dark matter annihilation rate can be estimated for arbitrary velocity-dependent cross sections in a convenient way using a limit of $P^2SAD$ to zero separation in physical space. We illustrate our method by computing the global and local subhalo annihilation boost to that of the smooth dark matter distribution in a Milky-Way-size halo. Two cases are considered, one where the cross section is velocity independent and one that approximates Sommerfeld-enhanced models. We find that the global boost is $sim10-30$, which is at the low end of current estimates (weakening expectations of large extragalactic signals), while the boost at the solar radius is below the percent level. We make our code to compute $P^2SAD$ publicly available, which can be used to estimate various observables that probe the nanostructure of dark matter haloes.
The mass aggregation and merger histories of present-day distinct haloes selected from the cosmological Millennium Simulations I and II are mapped into stellar mass aggregation and galaxy merger histories of central galaxies by using empirical stella
r-to-halo and stellar-to-gas mass relations. The growth of bulges driven by the galaxy mergers/interactions is calculated using dynamical prescriptions. The predicted bulge demographics at redshift z~0 is consistent with observations (Zavala+2012). Here we present the evolution of the morphological mix (traced by the bulge-to-total mass ratio, B/T) as a function of mass up to z=3. This mix remains qualitatively the same up to z~1: B/T<0.1 galaxies dominate at low masses, 0.1<B/T<0.45 at intermediate masses, and B/T>0.45 at large masses. At z>1, the fractions of disc-dominated and bulgeless galaxies increase strongly, and by z~2 the era of pure disc galaxies is reached. Bulge-dominated galaxies acquire such a morphology, and most of their mass, following a downsizing trend. Since our results are consistent with most of the recent observational studies of the morphological mix at different redshifts, a LCDM-based scenario of merger-driven bulge assembly does not seem to face critical issues. However, if the stellar-to-halo mass relation changes too little with redshift, then some tensions with observations appear.
We present a novel perspective on the clustering of dark matter in phase space by defining the particle phase space average density ($P^2SAD$) as a two-dimensional extension of the two-point correlation function averaged within a certain volume in ph
ase space. This statistics is a sensitive measure of small scale (sub-)structure of dark matter haloes. By analysing the structure of $P^2SAD$ in Milky-Way-size haloes using the Aquarius simulations, we find it to be nearly universal at small scales, i.e. small separations in phase space, where substructures dominate. This remarkable universality occurs across time and in regions of substantially different ambient densities (by nearly four orders of magnitude), with typical variations in $P^2SAD$ of a factor of a few. The maximum variations occur in regions where substructures have been strongly disrupted. The universality is also preserved across haloes of similar mass but diverse mass accretion histories and subhalo distributions. The universality is also broken at large scales, where the smooth dark matter distribution in the halo dominates. Although at small scales the structure of $P^2SAD$ is roughly described by a subhalo model, we argue that the simulation data is better fitted by a family of superellipse contours. This functional shape is inspired by a model that extends the stable clustering hypothesis into phase space. In a companion paper, we refine this model and show its advantages as a method to obtain predictions for non-gravitational signatures of dark matter.
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 curr
ent constraints. Given the seemingly stringent nature of the latter, it was thought for nearly a decade that SIDM would be viable only if the cross section for self-scattering was strongly velocity-dependent. However, it has recently been suggested that a constant cross section per unit mass of sigma_T/m~0.1cm^2/g 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 sigma_T/m~1cm^2/g 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.
Self-interacting dark matter offers an interesting alternative to collisionless dark matter because of its ability to preserve the large-scale success of the cold dark matter model, while seemingly solving its challenges on small scales. We present h
ere the first study of the expected dark matter detection signal taking into account different self-scattering models. We demonstrate that models with constant and velocity dependent cross sections, which are consistent with observational constraints, lead to distinct signatures in the velocity distribution, because non-thermalised features found in the cold dark matter distribution are thermalised through particle scattering. Depending on the model, self-interaction can lead to a 10% reduction of the recoil rates at high energies, corresponding to a minimum speed that can cause recoil larger than 300 km/s, compared to the cold dark matter case. At lower energies these differences are smaller than 5% for all models. The amplitude of the annual modulation signal can increase by up to 25%, and the day of maximum amplitude can shift by about two weeks with respect to the cold dark matter expectation. Furthermore, the exact day of phase reversal of the modulation signal can also differ by about a week between the different models. In general, models with velocity dependent cross sections peaking at the typical velocities of dwarf galaxies lead only to minor changes in the detection signals, whereas allowed constant cross section models lead to significant changes. We conclude that different self-interacting dark matter scenarios might be distinguished from each other through the details of direct detection signals. Furthermore, detailed constraints on the intrinsic properties of dark matter based on null detections, should take into account the possibility of self-scattering and the resulting effects on the detector signal.
The existence of an extended hot gaseous corona surrounding clusters, groups and massive galaxies is well established by observational evidence and predicted by current theories of galaxy formation. When a small galaxy collides with a larger one, the
ir coronae are the first to interact, producing disturbances that remove gas from the smaller system and settle it into the corona of the larger one. For a Milky-Way-size galaxy merging into a low-mass group, ram pressure stripping and the Kelvin-Helmholtz instability are the most relevant of these disturbances. We argue that the turbulence generated by the latter mixes the material of both coronae in the wake of the orbiting satellite creating a warm phase mixture with a cooling time a factor of several shorter than that of the ambient intragroup gas. We reach this conclusion using analytic estimates, as well as adiabatic and dissipative high resolution numerical simulations of a spherical corona subject to the ablation process of a constant velocity wind with uniform density and temperature. Although this is a preliminary analysis, our results are promising and we speculate that the mixture could potentially trigger in situ star formation and/or be accreted into the central galaxy as a cold gas flow resulting in a new mode of star formation in galaxy groups and clusters.
We use the Millennium I and II cosmological simulations to revisit the impact of mergers in the growth of bulges in central galaxies in the LCDM scenario. We seed galaxies within the growing CDM haloes using semi-empirical relations to assign stellar
and gaseous masses, and an analytic treatment to estimate the transfer of stellar mass to the bulge of the remnant after a galaxy merger. We find that this model roughly reproduces the observed correlation between the bulge-to-total (B/T) mass ratio and stellar mass in present-day central galaxies as well as their observed demographics, although low-mass B/T<0.1 (bulgeless) galaxies might be scarce relative to the observed abundance. In our merger-driven scenario, bulges have a composite population made of (i) stars acquired from infalling satellites, (ii) stars transferred from the primary disc due to merger-induced perturbations, and (iii) newly formed stars in starbursts triggered by mergers. We find that (i) and (ii) are the main channels of mass assembly, with the first being dominant for massive galaxies, creating large bulges with different stellar populations than those of the inner discs, while the second is dominant for intermediate/low-mass galaxies creating small bulges with similar stellar populations to the inner discs. We associate the dominion of the first (second) channel to classical (pseudo) bulges, and compare the predicted fractions to observations. We remark that our treatment does not include other mechanisms of bulge growth such as intrinsic secular disc instabilities or misaligned gas accretion. We find that the evolution of the stellar and gaseous contents of the satellite as it moves towards the central galaxy is a key ingredient in setting the morphology of the remnant, and that a good match to the observed bulge demographics occurs when this evolution proceeds closely to that of the central galaxy.
