The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($Lambda$CDM). Given the ambiguity
in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $Lambda$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $Lambda$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $Lambda$CDM predictions. In fact, ~10% of $Lambda$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively.
We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body simulation of structure formation in the $Lambda{rm CDM}{}$ model. COCO follows an approximately spherical region of radius $sim 17.4h^{-1},{rm Mpc}$ embedded in a much l
arger periodic cube that is followed at lower resolution. The high resolution volume has a particle mass of $1.135times10^5h^{-1}{rm M}_{odot}$ (60 times higher than the Millennium-II simulation). COCO gives the dark matter halo mass function over eight orders of magnitude in halo mass; it forms $sim 60$ haloes of galactic size, each resolved with about 10 million particles. We confirm the power-law character of the subhalo mass function, $bar{N}(>mu)proptomu^{-s}$, down to a reduced subhalo mass $M_{sub}/M_{200}equivmu=10^{-6}$, with a best-fit power-law index, $s=0.94$, for hosts of mass $langle M_{200}rangle=10^{12}h^{-1}{rm M}_{odot}$. The concentration-mass relation of COCO haloes deviates from a single power law for masses $M_{200}<textrm{a few}times 10^{8}h^{-1}{rm M}_{odot}$, where it flattens, in agreement with results by Sanchez-Conde et al. The host mass invariance of the reduced maximum circular velocity function of subhaloes, $ uequiv V_{max}/V_{200}$, hinted at in previous simulations, is clearly demonstrated over five orders of magnitude in host mass. Similarly, we find that the average, normalised radial distribution of subhaloes is approximately universal (i.e. independent of subhalo mass), as previously suggested by the Aquarius simulations of individual haloes. Finally, we find that at fixed physical subhalo size, subhaloes in lower mass hosts typically have lower central densities than those in higher mass hosts.
We investigate the angular and kinematic distributions of satellite galaxies around a large sample of bright isolated primaries in the spectroscopic and photometric catalogues of the Sloan Digital Sky Survey (SDSS). We detect significant anisotropy i
n the spatial distribution of satellites. To test whether this anisotropy could be related to the rotating disks of satellites recently found by Ibata et al. in a sample of SDSS galaxies, we repeat and extend their analysis. Ibata et al. found an excess of satellites on opposite sides of their primaries having anticorrelated radial velocities. We find that this excess is sensitive to small changes in the sample selection criteria which can greatly reduce its significance. In addition, we find no evidence for correspondingly correlated velocities for satellites observed on the same side of their primaries, which would be expected for rotating disks of satellites. We conclude that the detection of coherent rotation in the satellite population in current observational samples is not robust. We compare our data to the $Lambda$CDM Millennium simulations populated with galaxies according to the semi-analytic model of Guo et al. We find excellent agreement with the spatial distribution of satellites in the SDSS data and the lack of a strong signal from coherent rotation.
We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $Lambda$CDM Universe
and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*<5times10^{10}M_odot$. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass.
We use the Evolution and assembly of galaxies and their environments (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_odot$ of whi
ch we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $rho_*(r) propto r^{-3}$. Stars dominate the mass density for $r<10~rm{kpc}$, and, as a result, the total mass density profiles are steeper than the Navarro-Frenk-White (NFW) profile, in remarkable agreement with observations. The dark matter halo itself closely follows the NFW form at all resolved radii ($rgtrsim3.0~rm{kpc}$). The EAGLE BCGs have similar surface brightness and line-of-sight velocity dispersion profiles as the BCGs in the sample of Newman et al., which have the most detailed measurements currently available. After subtracting the contribution of the stars to the central density, Newman et al. infer significantly shallower slopes than the NFW value, in contradiction with the EAGLE results. We discuss possible reasons for this discrepancy, and conclude that an inconsistency between the kinematical model adopted by Newman et al. for their BCGs, which assumes isotropic stellar orbits, and the kinematical structure of the EAGLE BCGs, in which the orbital stellar anisotropy varies with radius and tends to be radially biased, could explain at least part of the discrepancy.
Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating f
actors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution.
We use the distribution of maximum circular velocities, $V_{max}$, of satellites in the Milky Way (MW) to constrain the virial mass, $M_{200}$, of the Galactic halo under an assumed prior of a $Lambda$CDM universe. This is done by analysing the subha
lo populations of a large sample of halos found in the Millennium II cosmological simulation. The observation that the MW has at most three subhalos with $V_{max}ge30 km/s$ requires a halo mass $M_{200}le1.4times10^{12} M_odot$, while the existence of the Magellanic Clouds (assumed to have $V_{max}ge60 km/s$) requires $M_{200}ge1.0times10^{12} M_odot$. The first of these conditions is necessary to avoid the too-big-to-fail problem highlighted by Boylan-Kolchin et al., while the second stems from the observation that massive satellites like the Magellanic Clouds are rare. When combining both requirements, we find that the MW halo mass must lie in the range $0.25 le M_{200}/(10^{12} M_odot) le 1.4$ at $90%$ confidence. The gap in the abundance of Galactic satellites between $30 km/sle V_{max} le 60 km/s$ places our galaxy in the tail of the expected satellite distribution.
The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only
simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.
We study the effect of baryons on the abundance of structures and substructures in a Lambda-CDM cosmology, using a pair of high resolution cosmological simulations from the GIMIC project. Both simulations use identical initial conditions, but while o
ne contains only dark matter, the other also includes baryons. We find that gas pressure, reionisation, supernova feedback, stripping, and truncated accretion systematically reduce the total mass and the abundance of structures below ~10^12 solar masses compared to the pure dark matter simulation. Taking this into account and adopting an appropriate detection threshold lowers the abundance of observed galaxies with maximum circular velocities below 100 km/s, significantly reducing the reported discrepancy between Lambda-CDM and the measured HI velocity function of the ALFALFA survey. We also show that the stellar-to-total mass ratios of galaxies with stellar masses of ~10^5 - 10^7 solar masses inferred from abundance matching of the (sub)halo mass function to the observed galaxy mass function increase by a factor of ~2. In addition, we find that an important fraction of low-mass subhaloes are completely devoid of stars. Accounting for the presence of dark subhaloes below 10^10 solar masses further reduces the abundance of observable objects, and leads to an additional increase in the inferred stellar-to-total mass ratio by factors of 2 - 10 for galaxies in haloes of 10^9 - 10^10 solar masses. This largely reconciles the abundance matching results with the kinematics of individual dwarf galaxies in Lambda-CDM. We propose approximate corrections to the masses of objects derived from pure dark matter calculations to account for baryonic effects.
Recent studies suggest that only three of the twelve brightest satellites of the Milky Way (MW) inhabit dark matter halos with maximum circular velocity, V_max, exceeding 30km/s. This is in apparent contradiction with the LCDM simulations of the Aqua
rius Project, which suggest that MW-sized halos should have at least 8 subhalos with V_max>30km/s. The absence of luminous satellites in such massive subhalos is thus puzzling and may present a challenge to the LCDM paradigm. We note, however, that the number of massive subhalos depends sensitively on the (poorly-known) virial mass of the Milky Way, and that their scarcity makes estimates of their abundance from a small simulation set like Aquarius uncertain. We use the Millennium Simulation series and the invariance of the scaled subhalo velocity function (i.e., the number of subhalos as a function of u, the ratio of subhalo V_max to host halo virial velocity, V_200) to secure improved estimates of the abundance of rare massive subsystems. In the range 0.1< u<0.5, N_sub(> u) is approximately Poisson-distributed about an average given by <N_sub>=10.2x( u/0.15)^(-3.11). This is slightly lower than in Aquarius halos, but consistent with recent results from the Phoenix Project. The probability that a LCDM halo has 3 or fewer subhalos with V_max above some threshold value, V_th, is then straightforward to compute. It decreases steeply both with decreasing V_th and with increasing halo mass. For V_th=30km/s, ~40% of M_halo=10^12 M_sun halos pass the test; fewer than 5% do so for M_halo>= 2x10^12 M_sun; and the probability effectively vanishes for M_halo>= 3x 10^12 M_sun. Rather than a failure of LCDM, the absence of massive subhalos might simply indicate that the Milky Way is less massive than is commonly thought.
