Source author record

S. Vegetti

S. Vegetti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2026arXiv

A possible challenge for Cold and Warm Dark Matter

Measuring the density profile and mass concentration of dark-matter haloes is a key test of the standard cold dark matter paradigm. Such objects are dark and thus challenging to characterise, but they can be studied via gravitational lensing. Recently, a million-solar-mass object was discovered superposed on an extended and extremely thin gravitational arc. Here we report on extensive tests of various assumptions for the mass density profile and redshift of this object. We find models that best describe the data have two components: an unresolved point-mass of radius $\leq10$ pc centred on an extended mass distribution with an almost constant surface density out to a truncation radius of 139 pc. These properties do not resemble any known astronomical object. However, if the object is dark matter-dominated, its structure is incompatible with cold dark matter models, but may be compatible with a self-interacting dark matter halo where the central region has collapsed to form a black hole. This detection could thus carry substantial implications for our current understanding of dark matter.

preprint2020arXiv

A dynamically cold disk galaxy in the early Universe

The extreme astrophysical processes and conditions that characterize the early Universe are expected to result in young galaxies that are dynamically different from those observed today. This is because the strong effects associated with galaxy mergers and supernova explosions would lead to most young star-forming galaxies being dynamically hot, chaotic and strongly unstable. Here we report the presence of a dynamically cold, but highly star-forming, rotating disk in a galaxy at redshift ($z$) 4.2, when the Universe was just 1.4 billion years old. Galaxy SPT-S J041839-4751.9 is strongly gravitationally lensed by a foreground galaxy at $z = 0.263$, and it is a typical dusty starburst, with global star-forming and dust properties that are in agreement with current numerical simulations and observations of its galaxy population. Interferometric imaging at a spatial resolution of about 60 pc reveals a ratio of rotational-to-random motions of $V/σ= 9.7\pm 0.4$, which is at least four times larger than expected from any galaxy evolution model at this epoch, but similar to the ratios of spiral galaxies in the local Universe. We derive a rotation curve with the typical shape of nearby massive spiral galaxies, which demonstrates that at least some young galaxies are dynamically akin to those observed in the local Universe, and only weakly affected by extreme physical processes.

preprint2020arXiv

Full of Orions: a 200-pc mapping of the interstellar medium in the redshift-3 lensed dusty star-forming galaxy SDP.81

We present a sub-kpc resolved study of the interstellar medium properties in SDP.81, a z=3.042 strongly gravitationally lensed dusty star-forming galaxy, based on high-resolution, multi-band ALMA observations of the FIR continuum, CO ladder and the [CII] line. Using a visibility-plane lens modelling code, we achieve a median source-plane resolution of ~200 pc. We use photon-dominated region (PDR) models to infer the physical conditions - far-UV field strength, density, and PDR surface temperature - of the star-forming gas on 200-pc scales, finding a FUV field strength of ~10^3-10^4 G0, gas density of ~10^5 cm^-3 and cloud surface temperatures up to 1500 K, similar to those in the Orion Trapezium region. The [CII] emission is significantly more extended than that FIR continuum: ~50 per cent of [CII] emission arises outside the FIR-bright region. The resolved [CII]/FIR ratio varies by almost 2 dex across the source, down to ~2x10^-4 in the star-forming clumps. The observed [CII]/FIR deficit trend is consistent with thermal saturation of the C+ fine-structure level occupancy at high gas temperatures. We make the source-plane reconstructions of all emission lines and continuum data publicly available.

preprint2020arXiv

SHARP -- VI. Evidence for CO (1-0) molecular gas extended on kpc-scales in AGN star forming galaxies at high redshift

We present a study of the stellar host galaxy, CO (1$-$0) molecular gas distribution and AGN emission on 50 to 500 pc-scales of the gravitationally lensed dust-obscured AGN MG J0751+2716 and JVAS B1938+666 at redshifts 3.200 and 2.059, respectively. By correcting for the lensing distortion using a grid-based lens modelling technique, we spatially locate the different emitting regions in the source plane for the first time. Both AGN host galaxies have 300 to 500 pc-scale size and surface brightness consistent with a bulge/pseudo-bulge, and 2 kpc-scale AGN radio jets that are embedded in extended molecular gas reservoirs that are 5 to 20 kpc in size. The CO (1$-$0) velocity fields show structures possibly associated with discs (elongated velocity gradients) and interacting objects (off-axis velocity components). There is evidence for a decrement in the CO (1$-$0) surface brightness at the location of the host galaxy, which may indicate radiative feedback from the AGN, or offset star formation.We find CO-H$_2$ conversion factors of around $α_{\rm CO} = 1.5\pm0.5$ (K km s$^{-1}$ pc$^2$)$^{-1}$, molecular gas masses of $> 3\times10^{10}$ M$_{\odot}$, dynamical masses of $\sim 10^{11}$ M$_{\odot}$ and gas fractions of around 60 per cent. The intrinsic CO line luminosities are comparable to those of unobscured AGN and dusty star-forming galaxies at similar redshifts, but the infrared luminosities are lower, suggesting that the targets are less efficient at forming stars. Therefore, they may belong to the AGN feedback phase predicted by galaxy formation models, because they are not efficiently forming stars considering their large amount of molecular gas.

preprint2020arXiv

Systematic errors in strong gravitational lensing reconstructions, a numerical simulation perspective

We present the analysis of a sample of twenty-four SLACS-like galaxy-galaxy strong gravitational lens systems with a background source and deflectors from the Illustris-1 simulation. We study the degeneracy between the complex mass distribution of the lenses, substructures, the surface brightness distribution of the sources, and the time delays. Using a novel inference framework based on Approximate Bayesian Computation, we find that for all the considered lens systems, an elliptical and cored power-law mass density distribution provides a good fit to the data. However, the presence of cores in the simulated lenses affects most reconstructions in the form of a Source Position Transformation. The latter leads to a systematic underestimation of the source sizes by 50 per cent on average, and a fractional error in $H_{0}$ of around $25_{-19}^{+37}$ per cent. The analysis of a control sample of twenty-four lens systems, for which we have perfect knowledge about the shape of the lensing potential, leads to a fractional error on $H_{0}$ of $12_{-3}^{+6}$ per cent. We find no degeneracy between complexity in the lensing potential and the inferred amount of substructures. We recover an average total projected mass fraction in substructures of $f_{\rm sub}<1.7-2.0\times10^{-3}$ at the 68 per cent confidence level in agreement with zero and the fact that all substructures had been removed from the simulation. Our work highlights the need for higher-resolution simulations to quantify the lensing effect of more realistic galactic potentials better, and that additional observational constraint may be required to break existing degeneracies.

preprint2016arXiv

SHARP - II. Mass structure in strong lenses is not necessarily dark matter substructure: A flux ratio anomaly from an edge-on disc

Gravitational lens flux-ratio anomalies provide a powerful technique for measuring dark matter substructure in distant galaxies. However, before using these flux-ratio anomalies to test galaxy formation models, it is imperative to ascertain that the given anomalies are indeed due to the presence of dark matter substructure and not due to some other component of the lensing galaxy halo or to propagation effects. Here we present the case of CLASS~B1555+375, which has a strong radio-wavelength flux-ratio anomaly. Our high-resolution near-infrared Keck~II adaptive optics imaging and archival Hubble Space Telescope data reveal the lensing galaxy in this system to have a clear edge-on disc component that crosses directly over the pair of images that exhibit the flux-ratio anomaly. We find that simple models that include the disc can reproduce the cm-wavelength flux-ratio anomaly without requiring additional dark matter substructure. Although further studies are required, our results suggest the assumption that all flux-ratio anomalies are due to a population of dark matter sub-haloes may be incorrect, and analyses that do not account for the full complexity of the lens macro-model may overestimate the substructure mass fraction in massive lensing galaxies.

preprint2015arXiv

ALMA imaging of SDP.81 - I. A pixelated reconstruction of the far-infrared continuum emission

We present a sub-50 pc-scale analysis of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data. We model both the mass distribution of the gravitational lensing galaxy and the pixelated surface brightness distribution of the background source using a novel Bayesian technique that fits the data directly in visibility space. We find the 1 and 1.3 mm dust emission to be magnified by a factor of u_tot = 17.6+/-0.4, giving an intrinsic total star-formation rate of 315+/-60 M_sol/yr and a dust mass of 6.4+/-1.5*10^8 M_sol. The reconstructed dust emission is found to be non-uniform, but composed of multiple regions that are heated by both diffuse and strongly clumped star-formation. The highest surface brightness region is a ~1.9*0.7 kpc disk-like structure, whose small extent is consistent with a potential size-bias in gravitationally lensed starbursts. Although surrounded by extended star formation, with a density of 20-30+/-10 M_sol/yr/kpc^2, the disk contains three compact regions with densities that peak between 120-190+/-20 M_sol/yr/kpc^2. Such star-formation rate densities are below what is expected for Eddington-limited star-formation by a radiation pressure supported starburst. There is also a tentative variation in the spectral slope of the different star-forming regions, which is likely due to a change in the dust temperature and/or opacity across the source.

preprint2015arXiv

ALMA imaging of SDP.81 - II. A pixelated reconstruction of the CO emission lines

We present a sub-100 pc-scale analysis of the CO molecular gas emission and kinematics of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data and a visibility-plane lens reconstruction technique. We find clear evidence for an excitation dependent structure in the unlensed molecular gas distribution, with emission in CO (5-4) being significantly more diffuse and structured than in CO (8-7). The intrinsic line luminosity ratio is r_8-7/5-4 = 0.30 +/- 0.04, which is consistent with other low-excitation starbursts at z ~ 3. An analysis of the velocity fields shows evidence for a star-forming disk with multiple velocity components that is consistent with a merger/post-coalescence merger scenario, and a dynamical mass of M(< 1.56 kpc) = 1.6 +/- 0.6 x 10^10 M_sol . Source reconstructions from ALMA and the Hubble Space Telescope show that the stellar component is offset from the molecular gas and dust components. Together with Karl G. Jansky Very Large Array CO (1-0) data, they provide corroborative evidence for a complex ~2 kpc-scale starburst that is embedded within a larger ~15 kpc structure.

preprint2015arXiv

Strong gravitational lensing with the SKA

Strong gravitational lenses provide an important tool to measure masses in the distant Universe, thus testing models for galaxy formation and dark matter; to investigate structure at the Epoch of Reionization; and to measure the Hubble constant and possibly w as a function of redshift. However, the limiting factor in all of these studies has been the currently small samples of known gravitational lenses (~10^2). The era of the SKA will transform our understanding of the Universe with gravitational lensing, particularly at radio wavelengths where the number of known gravitational lenses will increase to ~10^5. Here we discuss the technical requirements, expected outcomes and main scientific goals of a survey for strong gravitational lensing with the SKA. We find that an all-sky (3pi sr) survey carried out with the SKA1-MID array at an angular resolution of 0.25-0.5 arcsec and to a depth of 3 microJy / beam is required for studies of galaxy formation and cosmology with gravitational lensing. In addition, the capability to carryout VLBI with the SKA1 is required for tests of dark matter and studies of supermassive black holes at high redshift to be made using gravitational lensing.

preprint2015arXiv

The Massive End of the Stellar Mass Function

We derive average flux corrections to the \texttt{Model} magnitudes of the Sloan Digital Sky Survey (SDSS) galaxies by stacking together mosaics of similar galaxies in bins of stellar mass and concentration. Extra flux is detected in the outer low surface brightness part of the galaxies, leading to corrections ranging from 0.05 to 0.32 mag for the highest stellar mass galaxies. We apply these corrections to the MPA-JHU (Max-Planck Institute for Astrophysics - John Hopkins University) stellar masses for a complete sample of half a million galaxies from the SDSS survey to derive a corrected galaxy stellar mass function at $z=0.1$ in the stellar mass range $9.5<\log(M_\ast/M_\odot)<12.0$. We find that the flux corrections and the use of the MPA-JHU stellar masses have a significant impact on the massive end of the stellar mass function, making the slope significantly shallower than that estimated by Li \& White (2009), but steeper than derived by Bernardi et al. (2013). This corresponds to a mean comoving stellar mass density of galaxies with stellar masses $\log(M_\ast/M_\odot) \ge 11.0$ that is a factor of 3.36 larger than the estimate by Li \& White (2009), but is 43\% smaller than reported by Bernardi et al. (2013).

preprint2015arXiv

The X-shooter Lens Survey - II. Sample presentation and spatially resolved kinematics

We present the X-shooter Lens Survey (XLENS) data. The main goal of XLENS is to disentangle the stellar and dark matter content of massive early-type galaxies (ETGs), through combined strong gravitational lensing, dynamics and spectroscopic stellar population studies. The sample consists of 11 lens galaxies covering the redshift range from $0.1$ to $0.45$ and having stellar velocity dispersions between $250$ and $380\,\mathrm{km}\,\mathrm{s}^{-1}$. All galaxies have multi-band, high-quality HST imaging. We have obtained long-slit spectra of the lens galaxies with X-shooter on the VLT. We are able to disentangle the dark and luminous mass components by combining lensing and extended kinematics data-sets, and we are also able to precisely constrain stellar mass-to-light ratios and infer the value of the low-mass cut-off of the IMF, by adding spectroscopic stellar population information. Our goal is to correlate these IMF parameters with ETG masses and investigate the relation between baryonic and non-baryonic matter during the mass assembly and structure formation processes. In this paper we provide an overview of the survey, highlighting its scientific motivations, main goals and techniques. We present the current sample, briefly describing the data reduction and analysis process, and we present the first results on spatially resolved kinematics.

preprint2014arXiv

Inference of the Cold Dark Matter substructure mass function at z=0.2 using strong gravitational lenses

We present the results of a search for galaxy substructures in a sample of 11 gravitational lens galaxies from the Sloan Lens ACS Survey. We find no significant detection of mass clumps, except for a luminous satellite in the system SDSS J0956+5110. We use these non-detections, in combination with a previous detection in the system SDSS J0946+1006, to derive constraints on the substructure mass function in massive early-type host galaxies with an average redshift z ~ 0.2 and an average velocity dispersion of 270 km/s. We perform a Bayesian inference on the substructure mass function, within a median region of about 32 kpc squared around the Einstein radius (~4.2 kpc). We infer a mean projected substructure mass fraction $f = 0.0076^{+0.0208}_{-0.0052}$ at the 68 percent confidence level and a substructure mass function slope $α$ < 2.93 at the 95 percent confidence level for a uniform prior probability density on alpha. For a Gaussian prior based on Cold Dark Matter (CDM) simulations, we infer $f = 0 .0064^{+0.0080}_{-0.0042}$ and a slope of $α$ = 1.90$^{+0.098}_{-0.098}$ at the 68 percent confidence level. Since only one substructure was detected in the full sample, we have little information on the mass function slope, which is therefore poorly constrained (i.e. the Bayes factor shows no positive preference for any of the two models).The inferred fraction is consistent with the expectations from CDM simulations and with inference from flux ratio anomalies at the 68 percent confidence level.

preprint2014arXiv

Parametrizing the Stellar Haloes of Galaxies

We study the stellar haloes of galaxies out to 70-100 kpc as a function of stellar mass and galaxy type by stacking aligned $r$ and $g$ band images from a sample of 45508 galaxies from SDSS DR9 in the redshift range $0.06\,\le\,z\,\le\,0.1$ and in the mass range $10^{10.0} M_{\odot} < M_{*} < 10^{11.4} M_{\odot}$r. We derive surface brightness profiles to a depth of almost $μ_r \sim 32 \,\mathrm{mag\,arcsec}^{-2}$. We find that the ellipticity of the stellar halo is a function of galaxy stellar mass and that the haloes of high concentration ($C > 2.6$) galaxies are more elliptical than those of low concentration ($C < 2.6$) galaxies. The $g$-$r$ colour profile of high concentration galaxies reveals that the $g$-$r$ colour of the stellar population in the stellar halo is bluer than in the main galaxy, and the colour of the stellar halo is redder for higher mass galaxies. We further demonstrate that the full two-dimensional surface intensity distribution of our galaxy stacks can only be fit through multi-component Sérsic models. Using the fraction of light in the outer component of the models as a proxy for the fraction of accreted stellar light, we show that this fraction is a function of stellar mass and galaxy type. For high concentration galaxies, the fraction of accreted stellar light rises from $30\%$ to $70\%$ for galaxies in the stellar mass range from $10^{10.0} M_{\odot}$ to $10^{11.4} M_{\odot}$. The fraction of accreted light is much smaller in low concentration systems, increasing from $2\%$ to $25\%$ over the same mass range. This work provides important constraints for the theoretical understanding of the formation of stellar haloes of galaxies.

preprint2014arXiv

The physical nature of the 8 o'clock arc based on near-IR IFU spectroscopy with SINFONI

We present an analysis of near-infrared integral field unit spectroscopy of the 8 o'clock arc, a gravitationally lensed Lyman break galaxy, taken with SINFONI. We explore the shape of the spatially-resolved Hβ profile and demonstrate that we can decompose it into three components that partially overlap (spatially) but are distinguishable when we include dynamical information. We use existing B and H imaging from the Hubble Space Telescope to construct a rigorous lens model using a Bayesian grid based lens modelling technique. We apply this lens model to the SINFONI data cube to construct the de-lensed Hβ line-flux velocity and velocity dispersion maps of the galaxy. We find that the 8 o'clock arc has a complex velocity field that is not simply explained by a single rotating disk. The Hβ profile of the galaxy shows a blue-shifted wing suggesting gas outflows of ~ 200 km s^{-1}. We confirm that the 8 o'clock arc lies on the stellar mass--oxygen abundance--SFR plane found locally, but it has nevertheless significantly different gas surface density (a factor of 2--4 higher) and electron density in the ionized gas (five times higher) from those in similar nearby galaxies, possibly indicating a higher density interstellar medium for this galaxy.

preprint2012arXiv

Gravitational detection of a low-mass dark satellite at cosmological distance

The mass-function of dwarf satellite galaxies that are observed around Local Group galaxies substantially differs from simulations based on cold dark matter: the simulations predict many more dwarf galaxies than are seen. The Local Group, however, may be anomalous in this regard. A massive dark satellite in an early-type lens galaxy at z = 0.222 was recently found using a new method based on gravitational lensing, suggesting that the mass fraction contained in substructure could be higher than is predicted from simulations. The lack of very low mass detections, however, prohibited any constraint on their mass function. Here we report the presence of a 1.9 +/- 0.1 x 10^8 M_sun dark satellite in the Einstein-ring system JVAS B1938+666 at z = 0.881, where M_sun denotes solar mass. This satellite galaxy has a mass similar to the Sagittarius galaxy, which is a satellite of the Milky Way. We determine the logarithmic slope of the mass function for substructure beyond the local Universe to be alpha = 1.1^+0.6_-0.4, with an average mass-fraction of f = 3.3^+3.6_-1.8 %, by combining data on both of these recently discovered galaxies. Our results are consistent with the predictions from cold dark matter simulations at the 95 per cent confidence level, and therefore agree with the view that galaxies formed hierarchically in a Universe composed of cold dark matter.

preprint2012arXiv

SHARP - I. A high-resolution multi-band view of the infra-red Einstein ring of JVAS B1938+666

We present new mass models for the gravitational lens system B1938+666, using multi-wavelength data acquired from Keck adaptive optics (AO) and Hubble Space Telescope (HST) observations. These models are the first results from the Strong-lensing at High Angular Resolution Program (SHARP), a project designed to study known quadruple-image and Einstein ring lenses using high-resolution imaging, in order to probe their mass distributions in unprecedented detail. Here, we specifically highlight differences between AO- and HST-derived lens models, finding that -- at least when the lens and source galaxies are both bright and red, and the system has a high degree of circular symmetry -- AO-derived models place significantly tighter constraints on model parameters. Using this improved precision, we infer important physical properties about the B1938+666 system, including the mass density slope of the lensing galaxy (gamma = 2.045), the projected dark matter mass fraction within the Einstein radius (M_dark/M_lens = 0.55), and the total magnification factor of the source galaxy (~ 13). Additionally, we measure an upper-limit constraint on luminous substructure (M_V > 16.2), based on the non-detection of bright satellite galaxies in all data sets. Finally, we utilize the improved image resolution of the AO data to reveal the presence of faint arcs outside of the primary Einstein ring. The positions and orientations of these arcs raise the intriguing possibility that B1938+666 has a second source galaxy, located at a more distant redshift. However, future work is needed to verify this hypothesis.

preprint2010arXiv

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.

preprint2010arXiv

The mass distribution of a moderate redshift galaxy group and brightest group galaxy from gravitational lensing and kinematics

The gravitational lens system CLASS B2108+213 has two radio-loud lensed images separated by 4.56 arcsec. The relatively large image separation implies that the lensing is caused by a group of galaxies. In this paper, new optical imaging and spectroscopic data for the lensing galaxies of B2108+213 and the surrounding field galaxies are presented. These data are used to investigate the mass and composition of the lensing structure. The redshift and stellar velocity dispersion of the main lensing galaxy (G1) are found to be z = 0.3648 +/- 0.0002 and sigma_v = 325 +/- 25 km/s, respectively. The optical spectrum of the lensed quasar shows no obvious emission or absorption features and is consistent with a BL Lac type radio source. However, the tentative detection of the G-band and Mg-b absorption lines, and a break in the spectrum of the host galaxy of the lensed quasar gives a likely source redshift of z = 0.67. Spectroscopy of the field around B2108+213 finds 51 galaxies at a similar redshift to G1, thus confirming that there is a much larger structure at z ~ 0.365 associated with this system. The width of the group velocity distribution is 694 +/- 93 km/s, but is non-Gaussian, implying that the structure is not yet viralized. The main lensing galaxy is also the brightest group member and has a surface brightness profile consistent with a typical cD galaxy. A lensing and dynamics analysis of the mass distribution, which also includes the newly found group members, finds that the logarithmic slope of the mass density profile is on average isothermal inside the Einstein radius, but steeper at the location of the Einstein radius. This apparent change in slope can be accounted for if an external convergence gradient, representing the underlying parent halo of the galaxy group, is included in the mass model.

preprint2009arXiv

Detection of a Dark Substructure through Gravitational Imaging

We report the detection of a dark substructure through direct gravitational imaging - undetected in the HST-ACS F814W image - in the gravitational lens galaxy of SLACS SDSSJ0946+1006 (the "Double Einstein Ring"). The detection is based on a Bayesian grid reconstruction of the two-dimensional surface density of the galaxy inside an annulus around its Einstein radius (few kpc). [...] We confirm this detection by modeling the system including a parametric mass model with a tidally truncated pseudo-Jaffe density profile; in that case the substructure mass is M_sub=(3.51+-0.15)x10^9 Msun, located at (-0.651+-0.038,1.040+-0.034)'', precisely where also the surface density map shows a strong convergence peak. [...] We set a lower limit of (M/L)_V}>=120 (Msun/L}_V,sun (3-sigma) inside a sphere of 0.3 kpc centred on the substructure (r_tidal=1.1kpc). The result is robust under substantial changes in the model and the data-set (e.g. PSF, pixel number and scale, source and potential regularization, rotations and galaxy subtraction). Despite being at the limits of detectability, it can therefore not be attributed to obvious systematic effects. Our detection implies a dark matter mass fraction at the radius of the inner Einstein ring of f_CDM=2.15^{+2.05}_{-1.25} percent (68 percent C.L) in the mass range 4x10^6 Msun to 4x10^9 Msun assuming alpha=1.9+-0.1 (with dN/dm ~ m^-alpha). Assuming a flat prior on alpha, between 1.0 and 3.0, increases this to f_CDM=2.56^{+3.26}_{-1.50} percent (68 percent C.L). The likelihood ratio is 0.51 between our best value (f_CDM=0.0215) and that from simulations (f_sim=0.003). Hence the inferred mass fraction, admittedly based on a single lens system, is large but still consistent with predictions. [...]

preprint2009arXiv

Statistics of mass substructure from strong gravitational lensing: quantifying the mass fraction and mass function

A Bayesian statistical formalism is developed to quantify the level at which the mass function slope (alpha) and the projected cumulative mass fraction (f) of (CDM) substructure in strong gravitational-lens galaxies, with arcs or Einstein rings, can be recovered as function of the lens-survey parameters and the detection threshold of the substructure mass. The method is applied to different sets of mock data to explore a range of observational limits: (i) the number of lens galaxies in the survey, (ii) the mass threshold, Mlow, for the detection of substructures and (iii) the uncertainty of the measured substructure masses. We explore two different priors on the mass function slope: a uniform prior and a Gaussian prior with alpha = 1.90+-0.1. With a substructure detection threshold Mlow=3x10^8 Msun, the number of lenses available now (n_l=30), a true dark-matter mass fraction in (CDM) substructure <=1.0% and a prior of alpha = 1.90+-0.1, we find that the upper limit of f can be constrained down to a level <=1.0% (95% CL). In the case of a Gaussian prior on alpha, it is always possible to set stringent constraints on both parameters. We also find that lowering the detection threshold has the largest impact on the ability to recover alpha, because of the (expected) steep mass-function slope. In the future, thanks to new surveys with telescopes, such as SKA, LSST and JDEM and follow-up telescopes with high-fidelity data, a significant increase in the number of known lenses will allow us to recover the satellite population in its completeness. For example, a sample of 200 lenses, equivalent in data-quality to the Sloan Lens ACS Survey and a detection threshold of 10^8 Msun, allows one to determine f=0.5+-0.1% (68% CL) and alpha=1.90+-0.2 (68% CL).