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

A novel 3D technique to study the kinematics of lensed galaxies

We present a 3D Bayesian method to model the kinematics of strongly lensed galaxies from spatially-resolved emission-line observations. This technique enables us to simultaneously recover the lens-mass distribution and the source kinematics directly from the 3D data cube. We have tested this new method with simulated OSIRIS observations for nine star-forming lensed galaxies with different kinematic properties. The simulated rotation curves span a range of shapes which are prototypes of different morphological galaxy types, from dwarf to massive spiral galaxies. We have found that the median relative accuracy on the inferred lens and kinematic parameters are at the level of 1 and 2 per cent, respectively. We have also tested the robustness of the technique against different inclination angles, signal-to-noise ratios, the presence of warps or non-circular motions and we have found that the accuracy stays within a few per cent in most cases. This technique represents a significant step forward with respect to the methods used until now, as the lens parameters and the kinematics of the source are derived from the same 3D data. This enables us to study the possible degeneracies between the two and estimate the uncertainties on all model parameters consistently.Comment: Accepted for publication in MNRA

On the density profile of dark matter substructure in gravitational lens galaxies

We consider three extensions of the Navarro, Frenk and White (NFW) profile and investigate the intrinsic degeneracies among the density profile parameters on the gravitational lensing effect of satellite galaxies on highly magnified Einstein rings. In particular, we find that the gravitational imaging technique can be used to exclude specific regions of the considered parameter space, and therefore, models that predict a large number of satellites in those regions. By comparing the lensing degeneracy with the intrinsic density profile degeneracies, we show that theoretical predictions based on fits that are dominated by the density profile at larger radii may significantly over- or under-estimate the number of satellites that are detectable with gravitational lensing. Finally, using the previously reported detection of a satellite in the gravitational lens system JVAS B1938+666 as an example, we derive for this detected satellite values of r[subscript max] and v[subscript max] that are, for each considered profile, consistent within 1σ with the parameters found for the luminous dwarf satellites of the Milky Way and with a mass density slope γ < 1.6. We also find that the mass of the satellite within the Einstein radius as measured using gravitational lensing is stable against assumptions on the substructure profile. In the future thanks to the increased angular resolution of very long baseline interferometry at radio wavelengths and of the E-ELT in the optical we will be able to set tighter constraints on the number of allowed substructure profiles

Quantifying dwarf satellites through gravitational imaging: the case of SDSS J120602.09+514229.5

SDSS J120602.09+514229.5 is a gravitational lens system formed by a group of galaxies at redshift z=0.422 lensing a bright background galaxy at redshift z=2.001. The main peculiarity of this system is the presence of a luminous satellite near the Einstein radius, that slightly deforms the giant arc. This makes SDSS J120602.09+514229.5 the ideal system to test our grid-based Bayesian lens modelling method, designed to detect galactic satellites independently from their mass-to-light ratio, and to measure the mass of this dwarf galaxy despite its high redshift. Thanks to the pixelized source and potential reconstruction technique of Vegetti and Koopmans 2009a we are able to detect the luminous satellite as a local positive surface density correction to the overall smooth potential. Assuming a truncated Pseudo-Jaffe density profile, the satellite has a mass M=(2.75+-0.04)10^10 M_sun inside its tidal radius of r_t=0.68". We determine for the satellite a luminosity of L_B=(1.6+-0.8)10^9 L_sun, leading to a total mass-to-light ratio within the tidal radius of (M/L)_B=(17.2+-8.5) M_sun/L_sun. The central galaxy has a sub-isothermal density profile as in general is expected for group members. From the SDSS spectrum we derive for the central galaxy a velocity dispersion of sigma_kinem=380+-60 km/s within the SDSS aperture of diameter 3". The logarithmic density slope of gamma=1.7+0.25-0.30 (68% CL), derived from this measurement, is consistent within 1-sigma with the density slope of the dominant lens galaxy gamma~1.6, determined from the lens model. This paper shows how powerful pixelized lensing techniques are in detecting and constraining the properties of dwarf satellites at high redshift.Comment: Submitted to MNRAS; Abstract abridge

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