The impact of baryonic physics on the subhalo mass function and implications for gravitational lensing

We investigate the impact of baryonic physics on the subhalo population by analyzing the results of two recent hydrodynamical simulations (EAGLE and Illustris), which have very similar configuration, but a different model of baryonic physics. We concentrate on haloes with a mass between $10^{12.5}$ and $10^{14}M_{\odot}h^{-1}$ and redshift between 0.2 and 0.5, comparing with observational results and subhalo detections in early-type galaxy lenses. We compare the number and the spatial distribution of subhaloes in the fully hydro runs and in their dark matter only counterparts, focusing on the differences between the two simulations. We find that the presence of baryons reduces the number of subhaloes, especially at the low mass end ($\leq 10^{10}M_{\odot}h^{-1}$), by different amounts depending on the model. The variations in the subhalo mass function are strongly dependent on those in the halo mass function, which is shifted by the effect of stellar and AGN feedback. Finally, we search for analogues of the observed lenses (SLACS) in the simulations, selecting them in velocity dispersion and dynamical properties. We use the selected galaxies to quantify detection expectations based on the subhalo populations in the different simulations, calculating the detection probability and the predicted values for the projected dark matter fraction in subhaloes $f_{DM}$ and the slope of the mass function $\alpha$. We compare these values with those derived from subhalo detections in observations and conclude that the dark-matter-only and hydro EAGLE runs are both compatible with observational results, while results from the hydro Illustris run do not lie within the errors.Comment: 15 pages, 11 figures, accepted for publication in MNRA

Ellipsoidal halo finders and implications for models of triaxial halo formation

We describe an algorithm for identifying ellipsoidal haloes in numerical simulations, and quantify how the resulting estimates of halo mass and shape differ with respect to spherical halo finders. Haloes become more prolate when fit with ellipsoids, the difference being most pronounced for the more aspherical objects. Although the ellipsoidal mass is systematically larger, this is less than 10% for most of the haloes. However, even this small difference in mass corresponds to a significant difference in shape. We quantify these effects also on the initial mass and deformation tensors, on which most models of triaxial collapse are based. By studying the properties of protohaloes in the initial conditions, we find that models in which protohaloes are identified in Lagrangian space by three positive eigenvalues of the deformation tensor are tenable only at the masses well-above $M_*$. The overdensity $\delta$ within almost any protohalo is larger than the critical value associated with spherical collapse (increasing as mass decreases); this is in good qualitative agreement with models which identify haloes requiring that collapse have occured along all three principal axes, each axis having turned around from the universal expansion at a different time. The distributions of initial values are in agreement with the simplest predictions associated with ellipsoidal collapse, assuming initially spherical protohaloes, collapsed around random positions which were sufficiently overdense. However, most protohaloes are not spherical and departures from sphericity increase as protohalo mass decreases. [Abridged]Comment: 18 pages, 17 figures, accepted for publication in MNRA

Accretion of satellites onto central galaxies in clusters: merger mass ratios and orbital parameters

We study the statistical properties of mergers between central and satellite galaxies in galaxy clusters in the redshift range $0<z<1$, using a sample of dark-matter only cosmological N-body simulations from Le SBARBINE dataset. Using a spherical overdensity algorithm to identify dark-matter haloes, we construct halo merger trees for different values of the over-density $\Delta_c$. While the virial overdensity definition allows us to probe the accretion of satellites at the cluster virial radius $r_{vir}$, higher overdensities probe satellite mergers in the central region of the cluster, down to $\approx 0.06 r_{vir}$, which can be considered a proxy for the accretion of satellite galaxies onto central galaxies. We find that the characteristic merger mass ratio increases for increasing values of $\Delta_c$: more than $60\%$ of the mass accreted by central galaxies since $z\approx 1$ comes from major mergers. The orbits of satellites accreting onto central galaxies tend to be more tangential and more bound than orbits of haloes accreting at the virial radius. The obtained distributions of merger mass ratios and orbital parameters are useful to model the evolution of the high-mass end of the galaxy scaling relations without resorting to hydrodynamic cosmological simulations.Comment: accepted by MNRAS (minor comments

Universality of dark matter haloes shape over six decades in mass: Insights from the Millennium XXL and SBARBINE simulations

For the last 30 years many observational and theoretical evidences have shown that galaxy clusters are not spherical objects, and that their shape is much better described by a triaxial geometry. With the advent of multi-wavelength data of increasing quality, triaxial investigations of galaxy clusters is gathering a growing interest from the community, especially in the time of "precision cosmology". In this work, we aim to provide the first statistically significant predictions in the unexplored mass range above 3x10^14 Mo/h, using haloes from two redshifts (z=0 and z=1) of the Millennium XXL simulation. The size of this cosmological dark matter only simulation (4.1 Gpc) allows the formation of a statistically significant number of massive cluster scale haloes (about 500 with M>2x10^15 Mo/h and 780000 with M>10^14 Mo/h). Besides, we aim to extend this investigation to lower masses in order to look for universal predictions across nearly six orders of magnitude in mass, from 10^10 to almost 10^16 Mo/h. For this purpose we use the SBARBINE simulations, allowing to model haloes of masses starting from 10^10 Mo/h. We use an elliptical overdensity method to select haloes and compute the shapes of the unimodal ones (approximately 50%), while we discard the unrelaxed. The minor to major and intermediate to major axis ratio are found to be well described by simple functional forms. For a given mass we can fully characterize the shape of a halo and give predictions about the distribution of axis ratios for a given cosmology and redshift. Moreover, these results are in some disagreement with the findings of Jing & Suto (2002) which are widely used in the community even though they have to be extrapolated far beyond their original mass range. This "recipe" is made available to the community in this paper and in a dedicated web page.Comment: 13 pages, 16 figure

Flux-ratio anomalies from discs and other baryonic structures in the Illustris simulation

The flux ratios in the multiple images of gravitationally lensed quasars can provide evidence for dark matter substructure in the halo of the lensing galaxy if the flux ratios differ from those predicted by a smooth model of the lensing galaxy mass distribution. However, it is also possible that baryonic structures in the lensing galaxy, such as edge-on discs, can produce flux-ratio anomalies. In this work, we present the first statistical analysis of flux-ratio anomalies due to baryons from a numerical simulation perspective. We select galaxies with various morphological types in the Illustris simulation and ray-trace through the simulated halos, which include baryons in the main lensing galaxies but exclude any substructures, in order to explore the pure baryonic effects. Our ray-tracing results show that the baryonic components can be a major contribution to the flux-ratio anomalies in lensed quasars and that edge-on disc lenses induce the strongest anomalies. We find that the baryonic components increase the probability of finding high flux-ratio anomalies in the early-type lenses by about 8% and by about 10 - 20% in the disc lenses. The baryonic effects also induce astrometric anomalies in 13% of the mock lenses. Our results indicate that the morphology of the lens galaxy becomes important in the analysis of flux-ratio anomalies when considering the effect of baryons, and that the presence of baryons may also partially explain the discrepancy between the observed (high) anomaly frequency and what is expected due to the presence of subhalos as predicted by the CDM simulations.Comment: 16 pages, 11 figures, accepted by MNRA

The interplay of Self-Interacting Dark Matter and baryons in shaping the halo evolution

We use high-resolution hydrodynamical simulation to test the difference of halo properties in cold dark matter (CDM) and a self-interacting dark matter (SIDM) scenario with a constant cross-section of $\sigma^\text{T}/m_{\chi}=1\;\text{cm}^{2}\text{g}^{-1}$. We find that the interplay between dark matter self-interaction and baryonic physics induces a complex evolution of the halo properties, which depends on the halo mass and morphological type, as well as on the halo mass accretion history. While high mass haloes, selected as analogues of early-type galaxies, show cored profiles in the SIDM run, systems of intermediate mass and with a significant disk component can develop a profile that is similar or cuspier than in CDM. The final properties of SIDM haloes - measured at z=0.2 - correlate with the halo concentration and formation time, suggesting that the differences between different systems are due to the fact that we are observing the impact self-interaction. We also search for signatures of self-interacting dark matter in the lensing signal of the main haloes and hints of potential differences in the distribution of Einstein radii, which suggests that future wide-field survey might be able to distinguish between CDM and SIDM models on this basis. Finally, we find that the subhalo abundances are not altered in the adopted SIDM model with respect to CDM.Comment: 11 pages, 11 figures, accepted for publication in MNRA

Constraining the mass density of free-floating black holes using razor-thin lensing arcs

Strong lensing of active galactic nuclei in the radio can result in razor-thin arcs, with a thickness of less than a milli-arcsecond, if observed at the resolution achievable with very long baseline interferometry (VLBI). Such razor-thin arcs provide a unique window on the coarseness of the matter distribution between source and observer. In this paper, we investigate to what extent such razor-thin arcs can constrain the number density and mass function of `free-floating' black holes, defined as black holes that do not, or no longer, reside at the centre of a galaxy. These can be either primordial in origin or arise as by-products of the evolution of super-massive black holes in galactic nuclei. When sufficiently close to the line of sight, free-floating black holes cause kink-like distortions in the arcs, which are detectable by eye in the VLBI images as long as the black hole mass exceeds $\sim 1000$ Solar masses. Using a crude estimate for the detectability of such distortions, we analytically compute constraints on the matter density of free-floating black holes resulting from null-detections of distortions along a realistic, fiducial arc, and find them to be comparable to those from quasar milli-lensing. We also use predictions from a large hydrodynamical simulation for the demographics of free-floating black holes that are not primordial in origin, and show that their predicted mass density is roughly four orders of magnitude below the constraints achievable with a single razor-thin arc.Comment: 17 pages, 13 figures, 1 table, comments welcom

Properties and observables of massive galaxies in self-interacting dark matter cosmologies

We use hydrodynamical cosmological simulations to test the differences between cold and self-interacting dark matter models (CDM and SIDM) in the mass range of massive galaxies ($10^{12}M_{\odot}h^{-1}<M<10^{13.5}M_{\odot}h^{-1}$). We consider two SIDM models: one with constant cross section $\sigma/m_{\chi}=1\mathrm{cm^2g^{-1}}$ and one where the cross section is velocity-dependent. We analyse the halo density profiles and concentrations, comparing the predictions of dark-matter-only and hydrodynamical simulations in all scenarios. We calculate the best-fit Einasto profiles and compare the resulting parameters with previous studies and define the best-fit concentration-mass relations. We find that the inclusion of baryons reduces the differences between different dark matter models with respect to the DM-only case. In SIDM hydro runs, deviations from the CDM density profiles weakly depend on mass: the most massive systems ($M>10^{13}M_{\odot}h^{-1}$) show cored profiles, while the least massive ones ($M<10^{12.5}M_{\odot}h^{-1}$) have cuspier profiles. Finally, we compare the predictions of our simulations to observational results, by looking at the dark matter fractions and the distribution of strong lensing Einstein radii. We find that in SIDM the DM-fractions decrease more rapidly with increasing stellar mass than in CDM, leading to lower fractions at $M_{*}>10^{11}M_{\odot}$, a distinctive signature of SIDM. At the same time, the distribution of Einstein radii, derived from both CDM and SIDM hydro runs, is comparable to observed samples of strong lenses with $M>10^{13}M_{\odot}h^{-1}$. We conclude that the interplay between self-interaction and baryons can greatly reduce the expected differences between CDM and SIDM models at this mass scale, and that techniques able to separate the dark and luminous mass in the inner regions of galaxies are needed to constrain self-interactions.Comment: 13 pages, 12 figure

SHARP -- VII. New constraints on the dark matter free-streaming properties and substructure abundance from gravitationally lensed quasars

We present an analysis of seven strongly gravitationally lensed quasars and the corresponding constraints on the properties of dark matter. Our results are derived by modelling the lensed image positions and flux-ratios using a combination of smooth macro models and a population of low-mass haloes within the mass range 10^6 to 10^9 Msun. Our lens models explicitly include higher-order complexity in the form of stellar discs and luminous satellites, as well as low-mass haloes located along the observed lines of sight for the first time. Assuming a Cold Dark Matter (CDM) cosmology, we infer an average total mass fraction in substructure of f_sub = 0.012^{+0.007}_{-0.004} (68 per cent confidence limits), which is in agreement with the predictions from CDM hydrodynamical simulations to within 1 sigma. This result is closer to the predictions than those from previous studies that did not include line-of-sight haloes. Under the assumption of a thermal relic dark matter model, we derive a lower limit on the particle relic mass of m th > 5.58 keV (95 per cent confidence limits), which is consistent with a value of m_th > 5.3 keV from the recent analysis of the Ly-alpha forest. We also identify two main sources of possible systematic errors and conclude that deeper investigations in the complex structure of lens galaxies as well as the size of the background sources should be a priority for this field.Comment: 14 pages, 7 figures, accepted for publication in MNRA