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-plan
e 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.
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 gala
xy 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.
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 p
ossibly 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.
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 Ein
stein 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.
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 gr
id 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. [...]
We introduce a new adaptive and fully Bayesian grid-based method to model strong gravitational lenses with extended images. The primary goal of this method is to quantify the level of luminous and dark-mass substructure in massive galaxies, through t
heir effect on highly-magnified arcs and Einstein rings. The method is adaptive on the source plane, where a Delaunay tessellation is defined according to the lens mapping of a regular grid onto the source plane. The Bayesian penalty function allows us to recover the best non-linear potential-model parameters and/or a grid-based potential correction and to objectively quantify the level of regularization for both the source and the potential. In addition, we implement a Nested-Sampling technique to quantify the errors on all non-linear mass model parameters -- ... -- and allow an objective ranking of different potential models in terms of the marginalized evidence. In particular, we are interested in comparing very smooth lens mass models with ones that contain mass-substructures. The algorithm has been tested on a range of simulated data sets, created from a model of a realistic lens system. One of the lens systems is characterized by a smooth potential with a power-law density profile, twelve include a NFW dark-matter substructure of different masses and at different positions and one contains two NFW dark substructures with the same mass but with different positions. Reconstruction of the source and of the lens potential for all of these systems shows the method is able, in a realistic scenario, to identify perturbations with masses >=10^7 solar mass when located on the Einstein ring. For positions both inside and outside of the ring, masses of at least 10^9 solar mass are required (...).
