A Bayesian analysis of redshifted 21-cm HI signal and foregrounds: Simulations for LOFAR

Observations of the EoR with the 21-cm hyperfine emission of neutral hydrogen (HI) promise to open an entirely new window onto the formation of the first stars, galaxies and accreting black holes. In order to characterize the weak 21-cm signal, we need to develop imaging techniques which can reconstruct the extended emission very precisely. Here, we present an inversion technique for LOFAR baselines at NCP, based on a Bayesian formalism with optimal spatial regularization, which is used to reconstruct the diffuse foreground map directly from the simulated visibility data. We notice the spatial regularization de-noises the images to a large extent, allowing one to recover the 21-cm power-spectrum over a considerable $k_{perp}-k_{para}$ space in the range of $0.03,{rm Mpc^{-1}}<k_{perp}<0.19,{rm Mpc^{-1}}$ and $0.14,{rm Mpc^{-1}}<k_{para}<0.35,{rm Mpc^{-1}}$ without subtracting the noise power-spectrum. We find that, in combination with using the GMCA, a non-parametric foreground removal technique, we can mostly recover the spherically average power-spectrum within $2sigma$ statistical fluctuations for an input Gaussian random rms noise level of $60 , {rm mK}$ in the maps after 600 hrs of integration over a $10 , {rm MHz}$ bandwidth.

Measuring gravitational lens time delays using low-resolution radio monitoring observations

Obtaining lensing time delay measurements requires long-term monitoring campaigns with a high enough resolution (< 1 arcsec) to separate the multiple images. In the radio, a limited number of high-resolution interferometer arrays make these observations difficult to schedule. To overcome this problem, we propose a technique for measuring gravitational time delays which relies on monitoring the total flux density with low-resolution but high-sensitivity radio telescopes to follow the variation of the brighter image. This is then used to trigger high-resolution observations in optimal numbers which then reveal the variation in the fainter image. We present simulations to assess the efficiency of this method together with a pilot project observing radio lens systems with the Westerbork Synthesis Radio Telescope (WSRT) to trigger Very Large Array (VLA) observations. This new method is promising for measuring time delays because it uses relatively small amounts of time on high-resolution telescopes. This will be important because instruments that have high sensitivity but limited resolution, together with an optimum usage of followup high-resolution observations from appropriate radio telescopes may in the future be useful for gravitational lensing time delay measurements by means of this new method.

Ionospheric Power-Spectrum Tomography in Radio Interferometry

A tomographic method is described to quantify the three-dimensional power-spectrum of the ionospheric electron-density fluctuations based on radio-interferometric observations by a two-dimensional planar array. The method is valid to first-order Born approximation and might be applicable to correct observed visibilities for phase variations due to the imprint of the full three-dimensional ionosphere. It is shown that not the ionospheric electron density distribution is the primary structure to model in interferometry, but its autocorrelation function or equivalent its power-spectrum. An exact mathematical expression is derived that provides the three dimensional power-spectrum of the ionospheric electron-density fluctuations directly from a rescaled scattered intensity field and an incident intensity field convolved with a complex unit phasor that depends on the w-term and is defined on the full sky pupil plane. In the limit of a small field of view, the method reduces to the single phase screen approximation. Tomographic self-calibration can become important in high-dynamic range observations at low radio frequencies with wide-field antenna interferometers, because a three-dimensional ionosphere causes a spatially varying convolution of the sky, whereas a single phase screen results in a spatially invariant convolution. A thick ionosphere can therefore not be approximated by a single phase screen without introducing errors in the calibration process. By applying a Radon projection and the Fourier projection-slice theorem, it is shown that the phase-screen approach in three dimensions is identical to the tomographic method. Finally we suggest that residual speckle can cause a diffuse intensity halo around sources, due to uncorrectable ionospheric phase fluctuations in the short integrations, which could pose a fundamental limit on the dynamic range in long-integration images.

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.

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. [...]

The Structure & Dynamics of Massive Early-type Galaxies: On Homology, Isothermality and Isotropy inside one Effective Radius

Based on 58 SLACS strong-lens early-type galaxies with direct total-mass and stellar-velocity dispersion measurements, we find that inside one effective radius massive elliptical galaxies with M_eff >= 3x10^10 M_sun are well-approximated by a power-law ellipsoid with an average logaritmic density slope of <gamma_LD> = -dlog(rho_tot)/dlog(r)=2.085^{+0.025}_{-0.018} (random error on mean) for isotropic orbits with beta_r=0, +-0.1 (syst.) and sigma_gamma <= 0.20^{+0.04}_{-0.02} intrinsic scatter (all errors indicate the 68 percent CL). We find no correlation of gamma_LD with galaxy mass (M_eff), rescaled radius (i.e. R_einst/R_eff) or redshift, despite intrinsic differences in density-slope between galaxies. Based on scaling relations, the average logarithmic density slope can be derived in an alternative manner, fully independent from dynamics, yielding <gamma_SR>=1.959 +- 0.077. Agreement between the two values is reached for <beta_r> =0.45 +- 0.25, consistent with mild radial anisotropy. This agreement supports the robustness of our results, despite the increase in mass-to-light ratio with total galaxy mass: M_eff ~ L_{V,eff}^(1.363+-0.056). We conclude that massive early-type galaxies are structurally close-to homologous with close-to isothermal total density profiles (<=10 percent intrinsic scatter) and have at most some mild radial anisotropy. Our results provide new observational limits on galaxy formation and evolution scenarios, covering four Gyr look-back time.

Strong Gravitational Lensing as a Probe of Gravity, Dark-Matter and Super-Massive Black Holes

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

Bayesian Strong Gravitational-Lens Modeling on Adaptive Grids: Objective Detection of Mass Substructure in Galaxies

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 their 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 (...).