The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in $H_0$. We quantify any potential effect of the MST with flexible mass models that are maximally degenerate with H0. Our calculation is based on a new hierarchical approach in which the MST is only constrained by stellar kinematics. The approach is validated on hydrodynamically simulated lenses. We apply the method to the TDCOSMO sample of 7 lenses (6 from H0LiCOW) and measure $H_0=74.5^{+5.6}_{-6.1}$ km s$^{-1}{\rm Mpc}^{-1}$. In order to further constrain the deflector mass profiles, we then add imaging and spectroscopy for 33 strong gravitational lenses from the SLACS sample. For 9 of the SLAC lenses we use resolved kinematics to constrain the stellar anisotropy. From the joint analysis of the TDCOSMO+SLACS sample, we measure $H_0=67.4^{+4.1}_{-3.2}$ km s$^{-1}{\rm Mpc}^{-1}$, assuming that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. While our new analysis does not statistically invalidate the mass profile assumptions by H0LiCOW, and thus their $H_0$ measurement relying on those, it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on $H_0$ derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data, chiefly spatially resolved kinematics of lens galaxies.

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

We report new high-quality galaxy scale strong lens candidates found in the Kilo Degree Survey data release 4 using Machine Learning. We have developed a new Convolutional Neural Network (CNN) classifier to search for gravitational arcs, following the prescription by \cite{2019MNRAS.484.3879P} and using only $r-$band images. We have applied the CNN to two "predictive samples": a Luminous red galaxy (LRG) and a "bright galaxy" (BG) sample ($r<21$). We have found 286 new high probability candidates, 133 from the LRG sample and 153 from the BG sample. We have then ranked these candidates based on a value that combines the CNN likelihood to be a lens and the human score resulting from visual inspection (P-value) and we present here the highest 82 ranked candidates with P-values $\ge 0.5$. All these high-quality candidates have obvious arc or point-like features around the central red defector. Moreover, we define the best 26 objects, all with scores P-values $\ge 0.7$ as a "golden sample" of candidates. This sample is expected to contain very few false positives and thus it is suitable for follow-up observations. The new lens candidates come partially from the the more extended footprint adopted here with respect to the previous analyses, partially from a larger predictive sample (also including the BG sample). These results show that machine learning tools are very promising to find strong lenses in large surveys and more candidates that can be found by enlarging the predictive samples beyond the standard assumption of LRGs. In the future, we plan to apply our CNN to the data from next-generation surveys such as the Large Synoptic Survey Telescope, Euclid, and the Chinese Space Station Optical Survey.

Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, which includes mass-dependent star formation rates and ionising escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high-$z$ galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionisation such as ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at $z\approx9.1$ are excluded at $\gtrsim2\sigma$ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below $\sim2.6$ K across an IGM neutral fraction range of $0.15 \lesssim \bar{x}_{H{\scriptscriptstyle I}} \lesssim 0.6$. Note, these limits are only obtained from 141 hrs of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21CMMC in combining 21-cm observations with other high-$z$ probes to constrain the astrophysics of galaxies.

Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $\Lambda$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $\sigma$. We find that modifications of $\Lambda$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

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 $\alpha_{\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.

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

We use nine different galaxy formation scenarios in ten cosmological simulation boxes from the EAGLE suite of {\Lambda}CDM hydrodynamical simulations to assess the impact of feedback mechanisms in galaxy formation and compare these to observed strong gravitational lenses. To compare observations with simulations, we create strong lenses with $M_\star$ > $10^{11}$ $M_\odot$ with the appropriate resolution and noise level, and model them with an elliptical power-law mass model to constrain their total mass density slope. We also obtain the mass-size relation of the simulated lens-galaxy sample. We find significant variation in the total mass density slope at the Einstein radius and in the projected stellar mass-size relation, mainly due to different implementations of stellar and AGN feedback. We find that for lens selected galaxies, models with either too weak or too strong stellar and/or AGN feedback fail to explain the distribution of observed mass-density slopes, with the counter-intuitive trend that increasing the feedback steepens the mass density slope around the Einstein radius ($\approx$ 3-10 kpc). Models in which stellar feedback becomes inefficient at high gas densities, or weaker AGN feedback with a higher duty cycle, produce strong lenses with total mass density slopes close to isothermal (i.e. -d log({\rho})/d log(r) $\approx$ 2.0) and slope distributions statistically agreeing with observed strong lens galaxies in SLACS and BELLS. Agreement is only slightly worse with the more heterogeneous SL2S lens galaxy sample. Observations of strong-lens selected galaxies thus appear to favor models with relatively weak feedback in massive galaxies.

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at $1.42\,$GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen (Fialkov & Barkana 2019). We use the upper limit on the 21-cm signal from the Epoch of Reionization ($z=9.1$) based on $141\,$hours of observations with LOFAR (Mertens et al. 2020) to evaluate the contribution of the high redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming halos, we find (at 68% C.L.) that the cosmological radio background can be at most 0.7% of the CMB at $1.42\,$GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of $0.1-0.7$% over the CMB (at $1.42\,$GHz) is still allowed and could explain the EDGES Low-Band detection. We also constrain the thermal and ionization state of the gas at $z = 9.1$, and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES High-Band data, LOFAR data disfavour scenarios with inefficient X-ray sources and cases where the Universe was ionized by massive halos only.

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.

Ultra-compact massive galaxies UCMGs, i.e. galaxies with stellar masses $M_{*} > 8 \times 10^{10} M_{\odot}$ and effective radii $R_{e} < 1.5$ kpc, are very rare systems, in particular at low and intermediate redshifts. Their origin as well as their number density across cosmic time are still under scrutiny, especially because of the paucity of spectroscopically confirmed samples. We have started a systematic census of UCMG candidates within the ESO Kilo Degree Survey, together with a large spectroscopic follow-up campaign to build the largest possible sample of confirmed UCMGs. This is the third paper of the series and the second based on the spectroscopic follow-up program. Here, we present photometrical and structural parameters of 33 new candidates at redshifts $0.15 \lesssim z \lesssim 0.5$ and confirm 19 of them as UCMGs, based on their nominal spectroscopically inferred $M_{*}$ and $R_{e}$. This corresponds to a success rate of $\sim 58\%$, nicely consistent with our previous findings. The addition of these 19 newly confirmed objects, allows us to fully assess the systematics on the system selection, and finally reduce the number density uncertainties. Moreover, putting together the results from our current and past observational campaigns and some literature data, we build the largest sample of UCMGs ever collected, comprising 92 spectroscopically confirmed objects at $0.1 \lesssim z \lesssim 0.5$. This number raises to 116, allowing for a $3\sigma$ tolerance on the $M_{*}$ and $R_{e}$ thresholds for the UCMG definition. For all these galaxies we have estimated the velocity dispersion values at the effective radii which have been used to derive a preliminary mass-velocity dispersion correlation.

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$\sigma$ upper limit of $\Delta^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power.

We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space.

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.

We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $\Lambda$CDM cosmology with a uniform prior on $\Omega_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII).

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $\Lambda$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1\sigma$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3\sigma$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $\Lambda$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $\Lambda$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

This chapter addresses limitations to current 21-cm signal detection instruments, being related to the instrument, environment, signal processing or science, and what lies beyond the current horizon for 21-cm science, especially in the 2030s and beyond. We address how to overcome current challenges and drive the field forward, not only approaching a detection of the 21-cm signal but to a full characterisation of its parameterspace, in particular, probing an increasingly larger volume of k-modes (spatially and in redshifts). We also will shortly touch upon the kinds of questions that could drive such future endeavours.

The local expansion rate of the Universe is parametrized by the Hubble constant, $H_0$, the ratio between recession velocity and distance. Different techniques lead to inconsistent estimates of $H_0$. Observations of Type Ia supernovae (SNe) can be used to measure $H_0$, but this requires an external calibrator to convert relative distances to absolute ones. We use the angular diameter distance to strong gravitational lenses as a suitable calibrator, which is only weakly sensitive to cosmological assumptions. We determine the angular diameter distances to two gravitational lenses, $810^{+160}_{-130}$ and $1230^{+180}_{-150}$~Mpc, at redshifts of $z=0.295$ and $0.6304$. Using these absolute distances to calibrate 740 previously-measured relative distances to SNe, we measure the Hubble constant to be $H_0=82.4^{+8.4}_{-8.3} ~{\rm km\,s^{-1}\,Mpc^{-1}}$.

Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI2033-4723 ($z_{\rm src} = 1.662$, $z_{\rm lens}$ = 0.6575), one of the core targets of the H0LICOW project for which we present cosmological inferences in a companion paper (Rusu et al. 2019). We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens ($\sigma_{\rm {los}}= 250^{+15}_{-21}$ km/s) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to $i < 23$ mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalog, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify 2 galaxy groups and 3 galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al. (2019).

We analyse the mass density distribution in the centres of galaxies across five orders of magnitude in mass range. Using high-quality spiral galaxy rotation curves and infrared photometry from SPARC, we conduct a systematic study of their central dark matter fraction ($f_{\rm DM}$) and their mass density slope ($\alpha$), within their effective radius. We show that lower-mass spiral galaxies are more dark matter dominated and have more shallow mass density slopes when compared with more massive galaxies, which have density profiles closer to isothermal. Low-mass ($M_{\star} \lesssim 10^{10}\, \rm M_{\odot}$) gas-rich spirals span a wide range of \fdm\ values, but systematically lower than in gas-poor systems of similar mass. With increasing galaxy mass, the values of \fdm\ decrease and the density profiles steepen. In the most massive late-type gas-poor galaxies, a possible flattening of these trends is observed. When comparing these results to massive ($M_{\star} \gtrsim 10^{10}M_{\odot}$) elliptical galaxies from SPIDER and to dwarf ellipticals from SMACKED, these trends result to be inverted. Hence, the values of both $f_{\rm DM}$ and $\alpha$, as a function of $M_{\star}$, exhibit a U-shape trend. At a fixed stellar mass, the mass density profiles in dwarf ellipticals are steeper than in spirals. These trends can be understood by stellar feedback from a more prolonged star formation period in spirals, causing a transformation of the initial steep density cusp to a more shallow profile via differential feedback efficiency by supernovae, and by galaxy mergers or AGN feedback in higher-mass galaxies.

In this paper, we exploit the gravitational lensing effect to detect proper motion in the highly magnified gravitationally lensed source MG B2016+112. We find positional shifts up to 6 mas in the lensed images by comparing two Very Long Baseline Interferometric (VLBI) radio observations at 1.7 GHz that are separated by 14.359 years, and provide an astrometric accuracy of the order of tens of $\mu$as. From lens modelling, we exclude a shift in the lensing galaxy as the cause of the positional change of the lensed images, and we assign it to the background source. The source consists of four sub-components separated by $\sim 175$ pc, with proper motion of the order of tens $\mu$as yr$^{-1}$ for the two components at highest magnification ($\mu\sim350$) and of the order of a few mas yr$^{-1}$ for the two components at lower magnification ($\mu\sim2$). We propose single AGN and dual AGN scenarios to explain the source plane. Although, the latter interpretation is supported by the archival multi-wavelength properties of the object. In this case, MG B2016+112 would represent the highest redshift dual radio-loud AGN system discovered thus far, and would support the merger interpretation for such systems. Also, given the low probability ($\sim10^{-5}$) of detecting a dual AGN system that is also gravitationally lensed, if confirmed, this would suggest that such dual AGN systems must be more abundant in the early Universe than currently thought.

We analyse the accuracy of radio interferometric gridding of visibilities with the aim to quantify the Epoch of Reionization (EoR) 21-cm power spectrum bias caused by gridding, ultimately to determine the suitability of different imaging algorithms and gridding settings for 21-cm power spectrum analysis. We simulate realistic LOFAR data, and construct power spectra with convolutional gridding and w-stacking, w-projection, image domain gridding and without w-correction. These are compared against directly Fourier transformed data. The influence of oversampling, kernel size, w-quantization, kernel windowing function and image padding are quantified. The gridding excess power is measured with a foreground subtraction strategy, for which foregrounds have been subtracted using Gaussian progress regression, as well as with a foreground avoidance strategy. Constructing a power spectrum that has a bias significantly lower compared to the expected EoR signals is possible with the tested methods, but requires a kernel oversampling factor > 4000 and, when using w-correction, > 500 w-quantization levels. These values are higher than typical values used for imaging, but are computationally feasible. The kernel size and padding factor parameters are less crucial. Among the tested methods, image domain gridding shows the highest accuracy with the lowest imaging time. LOFAR 21-cm power spectrum results are not affected by gridding. Image domain gridding is overall the most suitable algorithm for 21-cm EoR experiments, including for future SKA EoR analyses. Nevertheless, convolutional gridding with tuned parameters results in sufficient accuracy. This holds also for w-stacking for wide-field imaging. The w-projection algorithm is less suitable because of the kernel oversampling requirements, and a faceting approach is unsuitable due to the resulting spatial discontinuities.

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

Due to ionosphere absorption and the interference by natural and artificial radio emissions, ground observation of the sky at the decameter or longer is very difficult. This unexplored part of electromagnetic spectrum has the potential of great discoveries, notably in the study of cosmic dark ages and dawn, but also in heliophysics and space weather, planets, cosmic ray and neutrinos, pulsar and interstellar medium, extragalactic radio sources, and even SETI. At a forum organized by the International Space Science Institute-Beijing (ISSI-BJ), we discussed the prospect of opening up this window for astronomical observations by using a constellation of small or micro-satellites. We discussed the past experiments and the current ones such as the low frequency payload on Chang'e-4 mission lander, relay satellite and the Longjiang satellite, and also the future DSL mission, which is a linear array on lunar orbit which can make synthesized map of the whole sky as well as measure the global spectrum. We also discuss the synergy with other experiments, including ground global experiments such as EDGES, SARAS, SCI-HI and High-z, PRIZM/Albatros, ground imaging facillities such as LOFAR and MWA, and space experiments such as SUNRISE, DARE/DAPPER and PRATUSH. We also discussed some technical aspects of the DSL concept.

Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz frequency range) using 14 hours of data obtained with the LOFAR-Low Band Antenna (LBA) array. We also demonstrate the application of a multiple pointing calibration technique to calibrate the LOFAR-LBA dual-pointing observations centred on the North Celestial Pole and the radio galaxy 3C220.3. We observe an unexplained excess of $\sim 30-50\%$ in Stokes $I$ noise compared to Stokes $V$ for the two observed fields, which decorrelates on $\gtrsim 12$ seconds and might have a physical origin. We show that enforcing smoothness of gain errors along frequency direction during calibration reduces the additional variance in Stokes $I$ compared Stokes $V$ introduced by the calibration on sub-band level. After subtraction of smooth foregrounds, we achieve a $2\sigma$ upper limit on the 21-cm power spectrum of $\Delta_{21}^2 < (14561\,\text{mK})^2$ at $k\sim 0.038\,h\,\text{cMpc}^{-1}$ and $\Delta_{21}^2 < (14886\,\text{mK})^2$ at $k\sim 0.038 \,h\,\text{cMpc}^{-1}$ for the 3C220 and NCP fields respectively and both upper limits are consistent with each other. The upper limits for the two fields are still dominated by systematics on most $k$ modes.

We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper.

Interferometry of the cosmic 21-cm signal is set to revolutionise our understanding of the Epoch of Reionisation (EoR) and the Cosmic Dawn (CD). The culmination of ongoing efforts will be the upcoming Square Kilometre Array (SKA), which will provide tomography of the 21-cm signal from the first billion years of our Universe. Using a galaxy formation model informed by high-$z$ luminosity functions, here we forecast the accuracy with which the first phase of SKA-low (SKA1-low) can constrain the properties of the unseen galaxies driving the astrophysics of the EoR and CD. We consider three observing strategies: (i) deep (1000h on a single field); (ii) medium-deep (100hr on 10 independent fields); and (iii) shallow (10hr on 100 independent fields). Using the 21-cm power spectrum as a summary statistic, and conservatively only using the 21-cm signal above the foreground wedge, we predict that all three observing strategies should recover astrophysical parameters to a fractional precision of $\sim 0.1$ -- 10 per cent. The reionisation history is recovered to an uncertainty of $\Delta z \lesssim 0.1$ (1$\sigma$) for the bulk of its duration. The medium-deep strategy, balancing thermal noise against cosmic variance, results in the tightest constraints, slightly outperforming the deep strategy. The shallow observational strategy performs the worst, with up to a $\sim 10$ -- 60 per cent increase in the recovered uncertainty. We note, however, that non-Gaussian summary statistics, tomography, as well as unbiased foreground removal would likely favour the deep strategy.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

WD J2356$-$209 is a cool metal-polluted white dwarf whose visible spectrum is dominated by a strong and broad sodium feature. Although discovered nearly two decades ago, no detailed and realistic analysis of this star had yet been realized. In the absence of atmosphere models taking into account the nonideal high-density effects arising at the photosphere of WD J2356$-$209, the origin of its unique spectrum had remained nebulous. We use the cool white dwarf atmosphere code presented in the first paper of this series to finally reveal the secrets of this peculiar object and details about the planetesimal that polluted its atmosphere. Thanks to the improved input physics of our models, we find a solution that is in excellent agreement with the photometric observations and the visible spectrum. Our solution reveals that the photosphere of WD J2356$-$209 has a number density ratio of $\log\,{\rm Na/Ca}= 1.0 \pm 0.2$, which is the highest ever found in a white dwarf. Since we do not know how long ago the accretion episode stopped (if it has), we cannot precisely determine the composition nor the mass of the accreted planetesimal. Nevertheless, all scenarios considered indicate that its composition is incompatible with that of chondrite-like material and that its mass was at least $10^{21}\,{\rm g}$.

Calibration is a key step in the signal processing pipeline of any radio astronomical instrument. The required sky, ionospheric and instrumental models for this step can suffer from various kinds of incompleteness. In this paper we analyze several important calibration methods, ignoring for now the ionosphere. The aim is to use established statistical and signal processing tools to provide a generic method to assess calibratability of an instrument. We show how currently popular calibration techniques differ in their assumptions and also discuss their theoretical commonalities. We also study the effect of only using a sub-set of baselines on the calibration and provide theoretical methods to analyze excess noise and biases that it might introduce. In order to simplify the physical interpretation of the results, we introduce a new signal processing model which is capable of modeling instrumental direction dependent effects and spectral smoothness of the individual receiver gain within a beam-formed station. The statistical properties of this model are then studied by deriving the Cram\'er--Rao bound (CRB). We finally define a mathematical framework for calibratability of an instrument based on the model used which is generic and can be used to study different instruments. These theoretical results are then verified using numerical simulations.

Quadruply lensed quasars are extremely rare objects, but incredibly powerful cosmological tools. Only few dozen are known in the whole sky. Here we present the spectroscopic confirmation of two new quadruplets WG0214-2105 and WG2100-4452 discovered by Agnello & Spiniello (2018) within the Dark Energy Survey (DES) public footprints. We have conducted spectroscopic follow-up of these systems with the Southern African Large Telescope as part of a program that aims at confirming the largest possible number of optically selected strong gravitational lensing systems in the Equatorial and Southern Hemisphere. For both systems, we present the spectra for the sources and deflectors that allowed us to estimate the source redshifts and unambiguously confirm their lensing nature. For the brighter deflector (WG2100-4452), we measure the stellar velocity dispersion from the spectrum. We also obtain photometry for both lenses, directly from DES multi-band images, isolating the lens galaxies from the quasar images. One of the quadruplets, WG0214-2105, was also observed by Pan-STARRS, allowing us to estimate the apparent brightness of each quasar image at two different epochs, and thus to find evidence for flux variability. This result could suggest a microlensing event for the faintest components, although intrinsic variability cannot be excluded with only two epochs. Finally, we present simple lens models for both quadruplets, obtaining Einstein radii, SIE velocity dispersions, ellipticities, and position angles of the lens systems, as well as time delay predictions assuming a concordance cosmological model.

We present a new sample of galaxy-scale strong gravitational-lens candidates, selected from 904 square degrees of Data Release 4 of the Kilo-Degree Survey (KiDS), i.e., the "Lenses in the Kilo-Degree Survey" (LinKS) sample. We apply two Convolutional Neural Networks (ConvNets) to $\sim88\,000$ colour-magnitude selected luminous red galaxies yielding a list of 3500 strong-lens candidates. This list is further down-selected via human inspection. The resulting LinKS sample is composed of 1983 rank-ordered targets classified as "potential lens candidates" by at least one inspector. Of these, a high-grade subsample of 89 targets is identified with potential strong lenses by all inspectors. Additionally, we present a collection of another 200 strong lens candidates discovered serendipitously from various previous ConvNet runs. A straightforward application of our procedure to future Euclid or LSST data can select a sample of $\sim3000$ lens candidates with less than 10 per cent expected false positives and requiring minimal human intervention.

We investigate the implications of interference detection for experiments that are pursuing a detection of the redshifted 21-cm signals from the Epoch of Reionization. Interference detection causes samples to be sporadically flagged and rejected. As a necessity to reduce the data volume, flagged samples are typically (implicitly) interpolated during time or frequency averaging or uv-gridding. This so-far unexplored systematic biases the 21-cm power spectrum, and it is important to understand this bias for current 21-cm experiments as well as the upcoming SKA Epoch of Reionization experiment. We analyse simulated data using power spectrum analysis and Gaussian process regression. We find that the combination of flagging and averaging causes tiny spectral fluctuations, resulting in `flagging excess power'. This excess power does not substantially average down over time and, without extra mitigation techniques, can exceed the power of realistic models of the 21-cm reionization signals in LOFAR observations. We mitigate the bias by i) implementing a novel way to average data using a Gaussian-weighted interpolation scheme; ii) using unitary instead of inverse-variance weighting of visibilities; and iii) using low-resolution forward modelling of the data. After these modifications, which have been integrated in the LOFAR EoR processing pipeline, the excess power reduces by approximately three orders of magnitude, and is no longer preventing a detection of the 21-cm signals.

We use a sample of 17 strong gravitational lens systems from the BELLS GALLERY survey to quantify the amount of low-mass dark matter haloes within the lensing galaxies and along their lines of sight, and to constrain the properties of dark matter. Based on a detection criterion of 10$\sigma$, we report no significant detection in any of the lenses. Using the sensitivity function at the 10-$\sigma$ level, we have calculated the predicted number of detectable cold dark matter (CDM) line-of-sight haloes to be $\mu_{l} = 1.17\pm1.08$, in agreement with our null detection. Assuming a detection sensitivity that improved to the level implied by a 5-$\sigma$ threshold, the expected number of detectable line-of-sight haloes rises to $\mu_l = 9.0\pm3.0$. Whilst the current data find zero detections at this sensitivity level (which has a probability of P$^{{\rm5}\sigma}_{{\rm CDM}}(n_{\rm det}=0)$=0.0001 and would be in strong tension with the CDM framework), we find that such a low detection threshold leads to many spurious detections and non-detections and therefore the current lack of detections is unreliable and requires data with improved sensitivity. Combining this sample with a subsample of 11 SLACS lenses, we constrain the half-mode mass to be $\log$(M$_{\rm hm}) < 12.26$ at the 2-$\sigma$ level. The latter is consistent with resonantly produced sterile neutrino masses m$_{\rm s} < 0.8$ keV at any value of the lepton asymmetry at the 2-$\sigma$ level.

Convolutional Neural Networks (ConvNets) are one of the most promising methods for identifying strong gravitational lens candidates in survey data. We present two ConvNet lens-finders which we have trained with a dataset composed of real galaxies from the Kilo Degree Survey (KiDS) and simulated lensed sources. One ConvNet is trained with single \textit{r}-band galaxy images, hence basing the classification mostly on the morphology. While the other ConvNet is trained on \textit{g-r-i} composite images, relying mostly on colours and morphology. We have tested the ConvNet lens-finders on a sample of 21789 Luminous Red Galaxies (LRGs) selected from KiDS and we have analyzed and compared the results with our previous ConvNet lens-finder on the same sample. The new lens-finders achieve a higher accuracy and completeness in identifying gravitational lens candidates, especially the single-band ConvNet. Our analysis indicates that this is mainly due to improved simulations of the lensed sources. In particular, the single-band ConvNet can select a sample of lens candidates with $\sim40\%$ purity, retrieving 3 out of 4 of the confirmed gravitational lenses in the LRG sample. With this particular setup and limited human intervention, it will be possible to retrieve, in future surveys such as Euclid, a sample of lenses exceeding in size the total number of currently known gravitational lenses.

We present results from an ongoing investigation using the Kilo Degree Survey (KiDS) on the VLT Survey Telescope (VST) to provide a census of ultra-compact massive galaxies (UCMGs), defined as galaxies with stellar masses $M_{\rm \star} > 8 \times 10^{10} \rm M_{\odot}$ and effective radii $R_{\rm e} < 1.5\,\rm kpc$. UCMGs, which are expected to have undergone very few merger events, provide a unique view on the accretion history of the most massive galaxies in the Universe. Over an effective sky area of nearly 330 square degrees, we select UCMG candidates from KiDS multi-colour images, which provide high quality structural parameters, photometric redshifts and stellar masses. Our sample of $\sim 1000$ photometrically selected UCMGs at $z < 0.5$ represents the largest sample of UCMG candidates assembled to date over the largest sky area. In this paper we present the first effort to obtain their redshifts using different facilities, starting with first results for 28 candidates with redshifts $z < 0.5$, obtained at NTT and TNG telescopes. We confirmed, as bona fide UCMGs, 19 out of the 28 candidates with new redshifts. A further 46 UCMG candidates are confirmed with literature spectroscopic redshifts (35 at $z < 0.5$), bringing the final cumulative sample of spectroscopically-confirmed lower-z UCMGs to 54 galaxies, which is the largest sample at redshifts below $0.5$. We use these spectroscopic redshifts to quantify systematic errors in our photometric selection, and use these to correct our UCMG number counts. We finally compare the results to independent datasets and simulations.

The 21-cm signal of neutral hydrogen - emitted during the Epoch of Reionization - promises to be an important source of information for the study of the infant universe. However, its detection is impossible without sufficient mitigation of other strong signals in the data, which requires an accurate knowledge of the instrument. Using the result of instrument calibration, a large part of the contaminating signals are removed and the resulting residual data is further analyzed in order to detect the 21-cm signal. Direction dependent calibration (DDC) can strongly affect the 21-cm signal, however, its effect has not been precisely quantified. In the analysis presented here we show how to exactly calculate what part of the 21-cm signal is removed as a result of the DDC. We also show how a-priori information about the frequency behavior of the instrument can be used to reduce signal suppression. The theoretical results are tested using a realistic simulation based on the LOFAR setup. Our results show that low-order smooth gain functions (e.g. polynomials) over a bandwidth of ~10\,MHz - over which the signal is expected to be stationary - is sufficient to allow for calibration with limited, quantifiable, signal suppression in its power spectrum. We also show mathematically and in simulations that more incomplete sky models lead to larger 21-cm signal suppression, even if the gain models are enforced to be fully smooth. This result has immediate consequences for current and future radio telescopes with non-identical station beams, where DDC might be necessary (e.g. SKA-low).

Time-delay strong lensing provides a unique way to directly measure the Hubble constant ($H_{0}$). The precision of the $H_{0}$ measurement depends on the uncertainties in the time-delay measurements, the mass distribution of the main deflector(s), and the mass distribution along the line of sight. Tie and Kochanek (2018) have proposed a new microlensing effect on time delays based on differential magnification of the coherent accretion disc variability of the lensed quasar. If real, this effect could significantly broaden the uncertainty on the time delay measurements by up to $30\%$ for lens systems such as PG1115+080, which have relatively short time delays and monitoring over several different epochs. In this paper we develop a new technique that uses the time-delay ratios and simulated microlensing maps within a Bayesian framework in order to limit the allowed combinations of microlensing delays and thus to lessen the uncertainties due to the proposed effect. We show that, under the assumption of Tie and Kochanek (2018), the uncertainty on the time-delay distance ($D_{\Delta t}$, which is proportional to 1/$H_{0}$) of short time-delay ($\sim18$ days) lens, PG1115+080, increases from $\sim7\%$ to $\sim10\%$ by simultaneously fitting the three time-delay measurements from the three different datasets across twenty years, while in the case of long time-delay ($\sim90$ days) lens, the microlensing effect on time delays is negligible as the uncertainty on $D_{\Delta t}$ of RXJ1131-1231 only increases from $\sim2.5\%$ to $\sim2.6\%$.

Radio detection of Ultra High Energetic Cosmic Rays and Neutrinos (UHECRv) which hit the Moon has been investigated in recent years. In preparation for near-future lunar science missions, we discuss technical requirements for radio experiments onboard lunar orbiters or on a lunar lander. We also develop an analysis of UHECRv aperture by including UHECv events occurring in the sub-layers of lunar regolith. It is verified that even using a single antenna onboard lunar orbiters or a few meters above the Moon's surface, dozens of lunar UHECRv events are detectable for one-year of observation at energy levels of $10^{18}$ eV to $10^{23}$ eV. Furthermore, it is shown that an antenna 3 meters above the Moon's surface could detect lower energy lunar UHECR events at the level of $10^{15}$ eV to $10^{18}$ eV which might not be detectable from lunar orbiters or ground-based observations.

We present milliarcsecond (mas) angular resolution observations of the gravitationally lensed radio source MG J0751+2716 (at z=3.2) obtained with global Very Long Baseline Interferometry (VLBI) at 1.65 GHz. The background object is highly resolved in the tangential and radial directions, showing evidence of both compact and extended structure across several gravitational arcs that are 200 to 600~mas in size. By identifying compact sub-components in the multiple images, we constrain the mass distribution of the foreground z=0.35 gravitational lens using analytic models for the main deflector [power-law elliptical mass model; $\rho(r) \propto r^{-\gamma}$, where $\gamma=2$ corresponds to isothermal] and for the members of the galaxy group. Moreover, our mass models with and without the group find an inner mass-density slope steeper than isothermal for the main lensing galaxy, with $\gamma_1 = 2.08 \pm 0.02$ and $\gamma_2 = 2.16 \pm 0.02$ at the 4.2$\sigma$ level and 6.8$\sigma$ level, respectively, at the Einstein radius ($b_1 = 0.4025 \pm 0.0008$ and $b_2 = 0.307 \pm 0.002$ arcsec, respectively). We find randomly distributed image position residuals of about 3 mas, which are much larger that the measurement errors ($40$ $\mu$as on average). This suggests that at the mas level, the assumption of a smooth mass distribution fails, requiring additional structure in the model. However, given the environment of the lensing galaxy, it is not clear whether this extra mass is in the form of sub-haloes within the lens or along the line of sight, or from a more complex halo for the galaxy group.

New methods have been recently developed to search for strong gravitational lenses, in particular lensed quasars, in wide-field imaging surveys. Here, we compare the performance of three different, morphology- and photometry- based methods to find lens candidates over the Kilo-Degree Survey (KiDS) DR3 footprint (440 deg$^2$). The three methods are: i) a multiplet detection in KiDS-DR3 and/or Gaia-DR1, ii) direct modeling of KiDS cutouts and iii) positional offsets between different surveys (KiDS-vs-Gaia, Gaia-vs-2MASS), with purpose-built astrometric recalibrations. The first benchmark for the methods has been set by the recovery of known lenses. We are able to recover seven out of ten known lenses and pairs of quasars observed in the KiDS DR3 footprint, or eight out of ten with improved selection criteria and looser colour pre-selection. This success rate reflects the combination of all methods together, which, taken individually, performed significantly worse (four lenses each). One movelty of our analysis is that the comparison of the performances of the different methods has revealed the pros and cons of the approaches and, most of all, the complementarities. We finally provide a list of high-grade candidates found by one or more methods, awaiting spectroscopic follow-up for confirmation. Of these, KiDS 1042+0023 is to our knowledge the first confirmed lensed quasar from KiDS, exhibiting two quasar spectra at the same source redshift at either sides of a red galaxy, with uniform flux-ratio $f\approx1.25$ over the wavelength range $0.45\mu\mathrm{m}<\lambda<0.75\mu\mathrm{m}.$

In this paper we introduce the SEAGLE (i.e. Simulating EAGLE LEnses) program, that approaches the study of galaxy formation through strong gravitational lensing, using a suite of high-resolution hydrodynamic simulations, Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. We introduce the simulation and analysis pipeline and present the first set of results from our analysis of early-type galaxies. We identify and extract an ensemble of simulated lens galaxies and use the GLAMER ray-tracing lensing code to create mock lenses similar to those observed in the SLACS and SL2S surveys, using a range of source parameters and galaxy orientations, including observational effects such as the Point-Spread-Function (PSF), pixelization and noise levels, representative of single-orbit observations with the Hubble Space Telescope (HST) using the ACS-F814W filter. We subsequently model these mock lenses using the code LENSED, treating them in the same way as observed lenses. We also estimate the mass model parameters directly from the projected surface mass density of the simulated galaxy, using an identical mass model family. We perform a three-way comparison of all the measured quantities with real lenses. We find the average total density slope of EAGLE lenses, $t=2.26\; (0.25\; \rm{rms})$ to be higher than SL2S, $t=2.16$ or SLACS, $t=2.08$. We find a very strong correlation between the external shear ($\gamma$) and the complex ellipticity ($\epsilon$), with $\gamma \sim \epsilon/4$. This correlation indicates a degeneracy in the lens mass modeling. We also see a dispersion between lens modeling and direct fitting results, indicating systematical biases.

Detecting and characterizing the Epoch of Reionization and Cosmic Dawn via the redshifted 21-cm hyperfine line of neutral hydrogen will revolutionize the study of the formation of the first stars, galaxies, black holes and intergalactic gas in the infant Universe. The wealth of information encoded in this signal is, however, buried under foregrounds that are many orders of magnitude brighter. These must be removed accurately and precisely in order to reveal the feeble 21-cm signal. This requires not only the modeling of the Galactic and extra-galactic emission, but also of the often stochastic residuals due to imperfect calibration of the data caused by ionospheric and instrumental distortions. To stochastically model these effects, we introduce a new method based on `Gaussian Process Regression' (GPR) which is able to statistically separate the 21-cm signal from most of the foregrounds and other contaminants. Using simulated LOFAR-EoR data that include strong instrumental mode-mixing, we show that this method is capable of recovering the 21-cm signal power spectrum across the entire range $k = 0.07 - 0.3 \ \rm{h\, cMpc^{-1}}$. The GPR method is most optimal, having minimal and controllable impact on the 21-cm signal, when the foregrounds are correlated on frequency scales $\gtrsim 3$\,MHz and the rms of the signal has $\sigma_{\mathrm{21cm}} \gtrsim 0.1\,\sigma_{\mathrm{noise}}$. This signal separation improves the 21-cm power-spectrum sensitivity by a factor $\gtrsim 3$ compared to foreground avoidance strategies and enables the sensitivity of current and future 21-cm instruments such as the {\sl Square Kilometre Array} to be fully exploited.

Contamination due to foregrounds (Galactic and Extra-galactic), calibration errors and ionospheric effects pose major challenges in detection of the cosmic 21 cm signal in various Epoch of Reionization (EoR) experiments. We present the results of a pilot study of a field centered on 3C196 using LOFAR Low Band (56-70 MHz) observations, where we quantify various wide field and calibration effects such as gain errors, polarized foregrounds, and ionospheric effects. We observe a `pitchfork' structure in the 2D power spectrum of the polarized intensity in delay-baseline space, which leaks into the modes beyond the instrumental horizon (EoR/CD window). We show that this structure largely arises due to strong instrumental polarization leakage ($\sim30\%$) towards {Cas\,A} ($\sim21$ kJy at 81 MHz, brightest source in northern sky), which is far away from primary field of view. We measure an extremely small ionospheric diffractive scale ($r_{\text{diff}} \approx 430$ m at 60 MHz) towards {Cas\,A} resembling pure Kolmogorov turbulence compared to $r_{\text{diff}} \sim3 - 20$ km towards zenith at 150 MHz for typical ionospheric conditions. This is one of the smallest diffractive scales ever measured at these frequencies. Our work provides insights in understanding the nature of aforementioned effects and mitigating them in future Cosmic Dawn observations (e.g. with SKA-low and HERA) in the same frequency window.

Having an accurate calibration method is crucial for any scientific research done by a radio telescope. The next generation radio telescopes such as the Square Kilometre Array (SKA) will have a large number of receivers which will produce exabytes of data per day. In this paper we propose new direction-dependent and independent calibration algorithms that, while requiring much less storage during calibration, converge very fast. The calibration problem can be formulated as a non-linear least square optimization problem. We show that combining a block-LDU decomposition with Gauss-Newton iterations produces systems of equations with convergent matrices. This allows significant reduction in complexity per iteration and very fast converging algorithms. We also discuss extensions to direction-dependent calibration. The proposed algorithms are evaluated using simulations.

For subspace estimation with an unknown colored noise, Factor Analysis (FA) is a good candidate for replacing the popular eigenvalue decomposition (EVD). Finding the unknowns in factor analysis can be done by solving a non-linear least square problem. For this type of optimization problems, the Gauss-Newton (GN) algorithm is a powerful and simple method. The most expensive part of the GN algorithm is finding the direction of descent by solving a system of equations at each iteration. In this paper we show that for FA, the matrices involved in solving these systems of equations can be diagonalized in a closed form fashion and the solution can be found in a computationally efficient way. We show how the unknown parameters can be updated without actually constructing these matrices. The convergence performance of the algorithm is studied via numerical simulations.

Constraining the sub-galactic matter-power spectrum on 1-10 kpc scales would make it possible to distinguish between the concordance $\Lambda$CDM model and various alternative dark-matter models due to the significantly different levels of predicted mass structure. Here, we demonstrate a novel approach to observationally constrain the population of overall law-mass density fluctuations in the inner regions of massive elliptical lens galaxies, based on the power spectrum of the associated surface-brightness perturbations observable in highly magnified galaxy-scale Einstein rings and gravitational arcs. The application of our method to the SLACS lens system SDSS J0252+0039 results in the following limits (at the 99 per cent confidence level) on the dimensionless convergence-power spectrum (and the associated standard deviation in aperture mass): $\Delta^{2}_{\delta\kappa}<1$ ($\sigma_{AM}< 0.8 \times 10^8 M_\odot$) on 0.5-kpc scale, $\Delta^{2}_{\delta\kappa}<0.1$ ($\sigma_{AM}< 1 \times 10^8 M_\odot$) on 1-kpc scale and $\Delta^{2}_{\delta\kappa}<0.01$ ($\sigma_{AM}< 3 \times 10^8 M_\odot$) on 3-kpc scale. The estimated effect of CDM sub-haloes lies considerably below these first observational upper-limit constraints on the level of inhomogeneities in the projected total mass distribution of galactic haloes. Future analysis for a larger sample of galaxy-galaxy strong lens systems will narrow down these constraints and rule out all cosmological models predicting a significantly larger level of clumpiness on these critical sub-galactic scales.

Direct detection of the Epoch of Reionization via the redshifted 21-cm line will have unprecedented implications on the study of structure formation in the early Universe. To fulfill this promise current and future 21-cm experiments will need to detect the weak 21-cm signal over foregrounds several order of magnitude greater. This requires accurate modeling of the galactic and extragalactic emission and of its contaminants due to instrument chromaticity, ionosphere and imperfect calibration. To solve for this complex modeling, we propose a new method based on Gaussian Process Regression (GPR) which is able to cleanly separate the cosmological signal from most of the foregrounds contaminants. We also propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere. Using this method, chromatic effects causing the so-called "wedge" are effectively eliminated (i.e. deconvolved) in the cylindrical ($k_{\perp}, k_{\parallel}$) power spectrum.

We present the full sample of 118 galaxy-scale strong-lens candidates in the Sloan Lens ACS (SLACS) Survey for the Masses (S4TM) Survey, which are spectroscopically selected from the final data release of the Sloan Digital Sky Survey. Follow-up Hubble Space Telescope (HST) imaging observations confirm that 40 candidates are definite strong lenses with multiple lensed images. The foreground lens galaxies are found to be early-type galaxies (ETGs) at redshifts 0.06 to 0.44, and background sources are emission-line galaxies at redshifts 0.22 to 1.29. As an extension of the SLACS Survey, the S4TM Survey is the first attempt to preferentially search for strong-lens systems with relatively lower lens masses than those in the pre-existing strong-lens samples. By fitting HST data with a singular isothermal ellipsoid model, we find total projected mass within the Einstein radius of the S4TM strong-lens sample ranges from $3 \times10^{10} M_{\odot}$ to $2 \times10^{11} M_{\odot}$. In [Shu15], we have derived the total stellar mass of the S4TM lenses to be $5 \times10^{10} M_{\odot}$ to $1 \times10^{12} M_{\odot}$. Both total enclosed mass and stellar mass of the S4TM lenses are on average almost a factor of 2 smaller than those of the SLACS lenses, which also represent typical mass scales of the current strong-lens samples. The extended mass coverage provided by the S4TM sample can enable a direct test, with the aid of strong lensing, for transitions in scaling relations, kinematic properties, mass structure, and dark-matter content trends of ETGs at intermediate-mass scales as noted in previous studies.

We present ALMA 2-mm continuum and CO (2-1) spectral line imaging of the gravitationally lensed z=0.654 star-forming/quasar composite RX J1131-1231 at 240-400 mas angular resolution. The continuum emission is found to be compact and coincident with the optical emission, whereas the molecular gas forms a complete Einstein ring, which shows strong differential magnification. The de-lensed source structure is determined on 400-pc resolution using a visibility-fitting lens modelling technique. The reconstructed molecular gas velocity-field is consistent with a rotating disk with a maximum rotational velocity of 280 km/s. From dynamical model fitting we find an enclosed mass M(r<5 kpc)=(1.46+/-0.31)*10^11 M_sol. The molecular gas distribution is highly structured, with clumps that are co-incident with higher gas velocity dispersion regions 40-50 km/s and with the intensity peaks in the optical emission, which are associated with sites of on-going turbulent star-formation. The peak in the CO (2-1) distribution is not co-incident with the AGN, where there is a paucity of molecular gas emission, possibly due to radiative feedback from the central engine. The intrinsic molecular gas luminosity is L'_CO=(1.2+/-0.3)*10^10 K km/s pc^2 and the inferred gas mass is M(H2)=(8.3+/-3.0)*10^10 M_sol, which given its dynamical mass is consistent with a CO-H2 conversion factor of alpha = 5.5+/-2.0 M_solar(K km/s pc^2)^-1. This suggests that the star-formation efficiency is dependent on the host galaxy morphology as opposed to the nature of the AGN. The far-infrared continuum spectral energy distribution shows evidence for heated dust, equivalent to an obscured star-formation rate of SFR=69^(+41)_(-25)*(7.3/u_IR)M_sol/yr, which demonstrates the composite star-forming/AGN nature of this system. RX J1131-1231

New and upcoming radio interferometers will produce unprecedented amounts of data that demand extremely powerful computers for processing. This is a limiting factor due to the large computational power and energy costs involved. Such limitations restrict several key data processing steps in radio interferometry. One such step is calibration where systematic errors in the data are determined and corrected. Accurate calibration is an essential component in reaching many scientific goals in radio astronomy and the use of consensus optimization that exploits the continuity of systematic errors across frequency significantly improves calibration accuracy. In order to reach full consensus, data at all frequencies need to be calibrated simultaneously. In the SKA regime, this can become intractable if the available compute agents do not have the resources to process data from all frequency channels simultaneously. In this paper, we propose a multiplexing scheme that is based on the alternating direction method of multipliers (ADMM) with cyclic updates. With this scheme, it is possible to simultaneously calibrate the full dataset using far fewer compute agents than the number of frequencies at which data are available. We give simulation results to show the feasibility of the proposed multiplexing scheme in simultaneously calibrating a full dataset when a limited number of compute agents are available.

The initial mass function (IMF) is a key ingredient in many studies of galaxy formation and evolution. Although the IMF is often assumed to be universal, there is continuing evidence that it is not universal. Spectroscopic studies that derive the IMF of the unresolved stellar populations of a galaxy often assume that this spectrum can be described by a single stellar population (SSP). To alleviate these limitations, in this paper we have developed a unique hierarchical Bayesian framework for modelling composite stellar populations (CSPs). Within this framework we use a parameterized IMF prior to regulate a direct inference of the IMF. We use this new framework to determine the number of SSPs that is required to fit a set of realistic CSP mock spectra. The CSP mock spectra that we use are based on semi-analytic models and have an IMF that varies as a function of stellar velocity dispersion of the galaxy. Our results suggest that using a single SSP biases the determination of the IMF slope to a higher value than the true slope, although the trend with stellar velocity dispersion is overall recovered. If we include more SSPs in the fit, the Bayesian evidence increases significantly and the inferred IMF slopes of our mock spectra converge, within the errors, to their true values. Most of the bias is already removed by using two SSPs instead of one. We show that we can reconstruct the variable IMF of our mock spectra for signal-to-noise ratios exceeding $\sim$75.

Direct detection of the Epoch of Reionization (EoR) via the red-shifted 21-cm line will have unprecedented implications on the study of structure formation in the infant Universe. To fulfill this promise, current and future 21-cm experiments need to detect this weak EoR signal in the presence of foregrounds that are several orders of magnitude larger. This requires extreme noise control and improved wide-field high dynamic-range imaging techniques. We propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere, or equivalently in the $uvw$-domain. The method uses the one-to-one relation between spherical waves and spherical harmonics (SpH). It consistently handles signals from the entire sky, and does not require a $w$-term correction. The spherical-harmonics coefficients represent the sky-brightness distribution and the visibilities in the $uvw$-domain, and provide a direct estimate of the spatial power spectrum. Using these spectrally-smooth SpH coefficients, bright foregrounds can be removed from the signal, including their side-lobe noise, which is one of the limiting factors in high dynamics range wide-field imaging. Chromatic effects causing the so-called "wedge" are effectively eliminated (i.e. deconvolved) in the cylindrical ($k_{\perp}, k_{\parallel}$) power spectrum, compared to a power spectrum computed directly from the images of the foreground visibilities where the wedge is clearly present. We illustrate our method using simulated LOFAR observations, finding an excellent reconstruction of the input EoR signal with minimal bias.

We introduce a novel technique, called "granulometry", to characterize and recover the mean size and the size distribution of HII regions from 21-cm tomography. The technique is easy to implement, but places the previously not very well defined concept of morphology on a firm mathematical foundation. The size distribution of the cold spots in 21-cm tomography can be used as a direct tracer of the underlying probability distribution of HII region sizes. We explore the capability of the method using large-scale reionization simulations and mock observational data cubes while considering capabilities of SKA1-low and a future extension to SKA2. We show that the technique allows the recovery of the HII region size distribution with a moderate signal-to-noise ratio from wide-field imaging ($\rm SNR\lesssim3$), for which the statistical uncertainty is sample variance dominated. We address the observational requirements on the angular resolution, the field-of-view, and the thermal noise limit for a successful measurement. To achieve a full scientific return from 21-cm tomography and to exploit a synergy with 21-cm power spectra, we suggest an observing strategy using wide-field imaging (several tens of square degrees) by an interferometric mosaicking/multi-beam observation with additional intermediate baselines (~2-4 km).

In the last decade the detection of individual massive dark matter sub-halos has been possible using potential correction formalism in strong gravitational lens imaging. Here we propose a statistical formalism to relate strong gravitational lens surface brightness anomalies to the lens potential fluctuations arising from dark matter distribution in the lens galaxy. We consider these fluctuations as a Gaussian random field in addition to the unperturbed smooth lens model. This is very similar to weak lensing formalism and we show that in this way we can measure the power spectrum of these perturbations to the potential. We test the method by applying it to simulated mock lenses of different geometries and by performing an MCMC analysis of the theoretical power spectra. This method can measure density fluctuations in early type galaxies on scales of 1-10 kpc at typical rms-levels of a percent, using a single lens system observed with the Hubble Space Telescope with typical signal-to-noise ratios obtained in a single orbit.

We study the dark matter (DM) assembly in the central regions of massive early-type galaxies up to $z\sim 0.65$. We use a sample of $\sim 3800$ massive ($\log M_{\rm \star}/M_{\rm \odot} > 11.2$) galaxies with photometry and structural parameters from 156 sq. deg. of the Kilo Degree Survey, and spectroscopic redshifts and velocity dispersions from SDSS. We obtain central total-to-stellar mass ratios, $M_{\rm dyn}/M_{\rm \star}$, and DM fractions, by determining dynamical masses, $M_{\rm dyn}$, from Jeans modelling of SDSS aperture velocity dispersions and stellar masses, $M_{\rm \star}$, from KiDS galaxy colours. We first show how the central DM fraction correlates with structural parameters, mass and density proxies, and demonstrate that most of the local correlations are still observed up to $z \sim 0.65$; at fixed $M_{\rm \star}$, local galaxies have larger DM fraction, on average, than their counterparts at larger redshift. We also interpret these trends with a non universal Initial Mass Function (IMF), finding a strong evolution with redshift, which contrast independent observations and is at odds with the effect of galaxy mergers. For a fixed IMF, the galaxy assembly can be explained, realistically, by mass and size accretion, which can be physically achieved by a series of minor mergers. We reproduce both the $R_{\rm e}-M_{\rm \star}$ and $M_{\rm dyn}/M_{\rm \star}-M_{\rm \star}$ evolution with stellar and dark mass changing at a different rate. This result suggests that the main progenitor galaxy is merging with less massive systems, characterized by a smaller $M_{\rm dyn}/M_{\rm \star}$, consistently with results from halo abundance matching.

The volume of data that will be produced by new-generation surveys requires automatic classification methods to select and analyze sources. Indeed, this is the case for the search for strong gravitational lenses, where the population of the detectable lensed sources is only a very small fraction of the full source population. We apply for the first time a morphological classification method based on a Convolutional Neural Network (CNN) for recognizing strong gravitational lenses in $255$ square degrees of the Kilo Degree Survey (KiDS), one of the current-generation optical wide surveys. The CNN is currently optimized to recognize lenses with Einstein radii $\gtrsim 1.4$ arcsec, about twice the $r$-band seeing in KiDS. In a sample of $21789$ colour-magnitude selected Luminous Red Galaxies (LRG), of which three are known lenses, the CNN retrieves 761 strong-lens candidates and correctly classifies two out of three of the known lenses. The misclassified lens has an Einstein radius below the range on which the algorithm is trained. We down-select the most reliable 56 candidates by a joint visual inspection. This final sample is presented and discussed. A conservative estimate based on our results shows that with our proposed method it should be possible to find $\sim100$ massive LRG-galaxy lenses at $z\lsim 0.4$ in KiDS when completed. In the most optimistic scenario this number can grow considerably (to maximally $\sim$2400 lenses), when widening the colour-magnitude selection and training the CNN to recognize smaller image-separation lens systems.

Verlinde derived gravity as an emergent force from the information flow, through two-dimensional surfaces and recently, by a priori postulating the entanglement of information in 3D space, he derived the effect of the gravitational potential from dark matter (DM) as the entropy displacement of dark energy by baryonic matter. In Emergent Gravity (EG) this apparent DM depends only on the baryonic mass distribution and the present-day value of the Hubble parameter. In this paper we test the EG proposition, formalized by Verlinde for a spherical and isolated mass distribution, using the central velocity dispersion, $\sigma$ and the light distribution in a sample of 4260 massive and local early-type galaxies (ETGs) from the SPIDER sample. Our results remain unaltered if we consider the sample of 807 roundest field galaxies. We derive the predictions by EG for the stellar mass-to-light ratio (M/L) and the Initial Mass Function (IMF), and compare them with the same inferences derived from a) DM-based models, b) MOND and c) stellar population models. We demonstrate that, consistently with a classical Newtonian framework with a DM halo component, or alternative theories of gravity as MOND, the central dynamics can be fitted if the IMF is assumed non-universal. The results can be interpreted with a IMF lighter than a standard Chabrier at low-$\sigma$, and bottom-heavier IMFs at larger $\sigma$. We find lower, but still acceptable, stellar M/L in EG theory, if compared with the DM-based NFW model and with MOND. The results from EG are comparable to what is found if the DM haloes are adiabatically contracted and with expectations from spectral gravity-sensitive features. If the strain caused by the entropy displacement would be not maximal, as adopted in the current formulation, then the dynamics of ETGs could be reproduced with larger M/L. (abridged)

Leakage of polarized Galactic diffuse emission into total intensity can potentially mimic the 21-cm signal coming from the epoch of reionization (EoR), as both of them might have fluctuating spectral structure. Although we are sensitive to the EoR signal only in small fields of view, chromatic sidelobes from further away can contaminate the inner region. Here, we explore the effects of leakage into the 'EoR window' of the cylindrically averaged power spectra (PS) within wide fields of view using both observation and simulation of the 3C196 and NCP fields, two observing fields of the LOFAR-EoR project. We present the polarization PS of two one-night observations of the two fields and find that the NCP field has higher fluctuations along frequency, and consequently exhibits more power at high-$k_\parallel$ that could potentially leak to Stokes $I$. Subsequently, we simulate LOFAR observations of Galactic diffuse polarized emission based on a model to assess what fraction of polarized power leaks into Stokes $I$ because of the primary beam. We find that the rms fractional leakage over the instrumental $k$-space is $0.35\%$ in the 3C196 field and $0.27\%$ in the NCP field, and it does not change significantly within the diameters of $15^\circ$, $9^\circ$ and $4^\circ$. Based on the observed PS and simulated fractional leakage, we show that a similar level of leakage into Stokes $I$ is expected in the 3C196 and NCP fields, and the leakage can be considered to be a bias in the PS.

Flux ratio anomalies in quasar lenses can be attributed to dark matter substructure surrounding the lensing galaxy and, thus, used to constrain the substructure mass fraction. Previous applications of this approach infer a substructure abundance that potentially in tension with the predictions of a $\Lambda$CDM cosmology. However, the assumption that all flux ratio anomalies are due to substructure is a strong one, and alternative explanations have not been fully investigated. Here, we use new high-resolution near-IR Keck~II adaptive optics imaging for the lens system CLASS B0712+472 to perform pixel-based lens modelling for this system and, in combination with new VLBA radio observations, show that the inclusion of the disc in the lens model can explain the flux ratio anomalies without the need for dark matter substructures. The projected disc mass comprises 16% of the total lensing mass within the Einstein radius and the total disc mass is $1.79 \times 10^{10} M_{sun}$. The case of B0712+472 adds to the evidence that not all flux ratio anomalies are due to dark subhaloes, and highlights the importance of taking the effects of baryonic structures more fully into account in order to obtain an accurate measure of the substructure mass fraction.

This paper is the first in a series of papers describing the impact of antenna instrumental artefacts on the 21-cm cosmology experiments to be carried out by the low frequency instrument (SKA1-LOW) of the Square Kilometre Array telescope (SKA), i.e., the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The smoothness of the passband response of the current log-periodic antenna being developed for the SKA1-LOW is analyzed using numerical electromagnetic simulations. The amplitude variations over the frequency range are characterized using low-order polynomials defined locally, in order to study the impact of the passband smoothness in the instrument calibration and CD/EoR Science. A solution is offered to correct a fast ripple found at 60~MHz during a test campaign at the SKA site at the Murchison Radio-astronomy Observatory, Western Australia in September 2015 with a minor impact on the telescope's performance and design. A comparison with the Hydrogen Epoch of Reionization Array antenna is also shown demonstrating the potential use of the SKA1-LOW antenna for the Delay Spectrum technique to detect the EoR.

We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $\Delta^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-$\sigma$) excess variance and a best 2-$\sigma$ upper limit of $\Delta^2_{\rm 21} < (79.6 {\rm mK})^2$ at $k=0.053$$h$cMpc$^{-1}$ in the range $z=$9.6-10.6. The excess variance decreases when optimizing the smoothness of the direction- and frequency-dependent gain calibration, and with increasing the completeness of the sky model. It is likely caused by (i) residual side-lobe noise on calibration baselines, (ii) leverage due to non-linear effects, (iii) noise and ionosphere-induced gain errors, or a combination thereof. Further analyses of the excess variance will be discussed in forthcoming publications.

The radio light-curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the cm-wave and mm-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the cm-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing by compact plasma structures in the Galactic interstellar medium, but the precise nature of these plasma structures remains unknown. Since the strength of a plasma lens evolves with wavelength $\lambda$ as $\lambda^2$, the presence of correlated variations at over a wide wavelength range casts doubt on the canonical ESE interpretation for J1415+1320. In this paper, we critically examine the evidence for plasma lensing in J1415+1320. We compute limits on the lensing strength, and the associated free-free opacity of the putative plasma lenses. We compare the observed and model ESE light curves, and also derive a lower limit on the lens distance based on the effects of parallax due to the Earth's orbit around the Sun. We conclude that plasma lensing is not a viable interpretation for J1415+1320's light curves and that symmetric U-shaped features in the radio light curves of extragalactic sources do not present {\em prima facie} evidence for ESEs. The methodology presented here is generic enough to be applicable to any plasma lensing candidate.

We present a new measurement of the Hubble Constant H0 and other cosmological parameters based on the joint analysis of three multiply-imaged quasar systems with measured gravitational time delays. First, we measure the time delay of HE0435-1223 from 13-year light curves obtained as part of the COSMOGRAIL project. Companion papers detail the modeling of the main deflectors and line of sight effects, and how these data are combined to determine the time-delay distance of HE 0435-1223. Crucially, the measurements are carried out blindly with respect to cosmological parameters in order to avoid confirmation bias. We then combine the time-delay distance of HE0435-1223 with previous measurements from systems B1608+656 and RXJ1131-1231 to create a Time Delay Strong Lensing probe (TDSL). In flat $\Lambda$CDM with free matter and energy density, we find $H_0$ = 71.9 +2.4 -3.0 km/s/Mpc and $\Omega_{\Lambda}$ = 0.62 +0.24 -0.35 . This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H0. We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The TDSL and Planck joint constraints are $H_0$ = 69.2 +1.4 -2.2 km/s/Mpc, $\Omega_{\Lambda}$ = 0.70 +0.01 -0.01 and $\Omega_k$ = 0.003 +0.004 -0.006 in open $\Lambda$CDM and $H_0$ = 79.0 +4.4 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.02 -0.03 and $w$ = -1.38 +0.14 -0.16 in flat $w$CDM. Combined with Planck and Baryon Acoustic Oscillation data, when relaxing the constraints on the numbers of relativistic species we find $N_{eff}$ = 3.34 +0.21 -0.21 and when relaxing the total mass of neutrinos we find 0.182 eV. In an open $w$CDM in combination with Planck and CMB lensing we find $H_0$ = 77.9 +5.0 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.03 -0.03, $\Omega_k$ = -0.003 +0.004 -0.004 and $w$ = -1.37 +0.18 -0.23.

Galaxies located in the environment or on the line of sight towards gravitational lenses can significantly affect lensing observables, and can lead to systematic errors on the measurement of $H_0$ from the time-delay technique. We present the results of a systematic spectroscopic identification of the galaxies in the field of view of the lensed quasar HE0435-1223, using the W. M. Keck, Gemini and ESO-Very Large telescopes. Our new catalog triples the number of known galaxy redshifts in the vicinity of the lens, expanding to 102 the number of measured redshifts for galaxies separated by less than 3 arcmin from the lens. We complement our catalog with literature data to gather redshifts up to 15 arcmin from the lens, and search for galaxy groups or clusters projected towards HE0435-1223. We confirm that the lens is a member of a small group that includes at least 12 galaxies, and find 8 other group candidates near the line of sight of the lens. The flexion shift, namely the shift of lensed images produced by high order perturbation of the lens potential, is calculated for each galaxy/group and used to identify which objects produce the largest perturbation of the lens potential. This analysis demonstrates that i) at most three of the five brightest galaxies projected within 12 arcsec of the lens need to be explicitly used in the lens models, and ii) the groups can be treated in the lens model as an external tidal field (shear) contribution.

Strong gravitational lens systems with time delays between the multiple images allow measurements of time-delay distances, which are primarily sensitive to the Hubble constant that is key to probing dark energy, neutrino physics, and the spatial curvature of the Universe, as well as discovering new physics. We present H0LiCOW ($H_0$ Lenses in COSMOGRAIL's Wellspring), a program that aims to measure $H_0$ with $<3.5\%$ uncertainty from five lens systems (B1608+656, RXJ1131-1231, HE0435-1223, WFI2033-4723 and HE1104-1805). We have been acquiring (1) time delays through COSMOGRAIL and Very Large Array monitoring, (2) high-resolution Hubble Space Telescope imaging for the lens mass modeling, (3) wide-field imaging and spectroscopy to characterize the lens environment, and (4) moderate-resolution spectroscopy to obtain the stellar velocity dispersion of the lenses for mass modeling. In cosmological models with one-parameter extension to flat $\Lambda$CDM, we expect to measure $H_0$ to $<3.5\%$ in most models, spatial curvature $\Omega_{\rm k}$ to 0.004, $w$ to 0.14, and the effective number of neutrino species to 0.2 (1$\sigma$ uncertainties) when combined with current CMB experiments. These are, respectively, a factor of $\sim15$, $\sim2$, and $\sim1.5$ tighter than CMB alone. Our data set will further enable us to study the stellar initial mass function of the lens galaxies, and the co-evolution of supermassive black holes and their host galaxies. This program will provide a foundation for extracting cosmological distances from the hundreds of time-delay lenses that are expected to be discovered in current and future surveys.

We present an exploration of the mass structure of a sample of 12 strongly lensed massive, compact early-type galaxies at redshifts $z\sim0.6$ to provide further possible evidence for their inside-out growth. We obtain new ESI/Keck spectroscopy and infer the kinematics of both lens and source galaxies, and combine these with existing photometry to construct (a) the fundamental plane (FP) of the source galaxies and (b) physical models for their dark and luminous mass structure. We find their FP to be tilted towards the virial plane relative to the local FP, and attribute this to their unusual compactness, which causes their kinematics to be totally dominated by the stellar mass as opposed to their dark matter; that their FP is nevertheless still inconsistent with the virial plane implies that both the stellar and dark structure of early-type galaxies is non-homologous. We also find the intrinsic scatter of their FP to be comparable to the local value, indicating that variations in the stellar mass structure outweight variations in the dark halo in the central regions of early-type galaxies. Finally, we show that inference on the dark halo structure -- and, in turn, the underlying physics -- is sensitive to assumptions about the stellar initial mass function (IMF), but that physically-motivated assumptions about the IMF imply haloes with sub-NFW inner density slopes, and may present further evidence for the inside-out growth of compact early-type galaxies via minor mergers and accretion.

Strong gravitational lenses with measured time delays between the multiple images allow a direct measurement of the time-delay distance to the lens, and thus a measure of cosmological parameters, particularly the Hubble constant, $H_{0}$. We present a blind lens model analysis of the quadruply-imaged quasar lens HE 0435-1223 using deep Hubble Space Telescope imaging, updated time-delay measurements from the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL), a measurement of the velocity dispersion of the lens galaxy based on Keck data, and a characterization of the mass distribution along the line of sight. HE 0435-1223 is the third lens analyzed as a part of the $H_{0}$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) project. We account for various sources of systematic uncertainty, including the detailed treatment of nearby perturbers, the parameterization of the galaxy light and mass profile, and the regions used for lens modeling. We constrain the effective time-delay distance to be $D_{\Delta t} = 2612_{-191}^{+208}~\mathrm{Mpc}$, a precision of 7.6%. From HE 0435-1223 alone, we infer a Hubble constant of $H_{0} = 73.1_{-6.0}^{+5.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$ assuming a flat $\Lambda$CDM cosmology. The cosmographic inference based on the three lenses analyzed by H0LiCOW to date is presented in a companion paper (H0LiCOW Paper V).

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

We present a new sample of strong gravitational lens systems where both the foreground lenses and background sources are early-type galaxies. Using imaging from HST/ACS and Keck/NIRC2, we model the surface brightness distributions and show that the sources form a distinct population of massive, compact galaxies at redshifts $0.4 \lesssim z \lesssim 0.7$, lying systematically below the size-mass relation of the global elliptical galaxy population at those redshifts. These may therefore represent relics of high-redshift red nuggets or their partly-evolved descendants. We exploit the magnifying effect of lensing to investigate the structural properties, stellar masses and stellar populations of these objects with a view to understanding their evolution. We model these objects parametrically and find that they generally require two S\'ersic components to properly describe their light profiles, with one more spheroidal component alongside a more envelope-like component, which is slightly more extended though still compact. This is consistent with the hypothesis of the inside-out growth of these objects via minor mergers. We also find that the sources can be characterised by red-to-blue colour gradients as a function of radius which are stronger at low redshift -- indicative of ongoing accretion -- but that their environments generally appear consistent with that of the general elliptical galaxy population, contrary to recent suggestions that these objects are predominantly associated with clusters.

Recent studies based on the integrated light of distant galaxies suggest that the initial mass function (IMF) might not be universal. Variations of the IMF with galaxy type and/or formation time may have important consequences for our understanding of galaxy evolution. We have developed a new stellar population synthesis (SPS) code specifically designed to reconstruct the IMF. We implement a novel approach combining regularization with hierarchical Bayesian inference. Within this approach we use a parametrized IMF prior to regulate a direct inference of the IMF. This direct inference gives more freedom to the IMF and allows the model to deviate from parametrized models when demanded by the data. We use Markov Chain Monte Carlo sampling techniques to reconstruct the best parameters for the IMF prior, the age, and the metallicity of a single stellar population. We present our code and apply our model to a number of mock single stellar populations with different ages, metallicities, and IMFs. When systematic uncertainties are not significant, we are able to reconstruct the input parameters that were used to create the mock populations. Our results show that if systematic uncertainties do play a role, this may introduce a bias on the results. Therefore, it is important to objectively compare different ingredients of SPS models. Through its Bayesian framework, our model is well-suited for this.

Accurate and precise measurements of the Hubble constant are critical for testing our current standard cosmological model and revealing possibly new physics. With Hubble Space Telescope (HST) imaging, each strong gravitational lens system with measured time delays can allow one to determine the Hubble constant with an uncertainty of $\sim 7\%$. Since HST will not last forever, we explore adaptive-optics (AO) imaging as an alternative that can provide higher angular resolution than HST imaging but has a less stable point spread function (PSF) due to atmospheric distortion. To make AO imaging useful for time-delay-lens cosmography, we develop a method to extract the unknown PSF directly from the imaging of strongly lensed quasars. In a blind test with two mock data sets created with different PSFs, we are able to recover the important cosmological parameters (time-delay distance, external shear, lens mass profile slope, and total Einstein radius). Our analysis of the Keck AO image of the strong lens system RXJ1131-1231 shows that the important parameters for cosmography agree with those based on HST imaging and modeling within 1-$\sigma$ uncertainties. Most importantly, the constraint on the model time-delay distance by using AO imaging with $0.045"$resolution is tighter by $\sim 50\%$ than the constraint of time-delay distance by using HST imaging with $0.09"$when a power-law mass distribution for the lens system is adopted. Our PSF reconstruction technique is generic and applicable to data sets that have multiple nearby point sources, enabling scientific studies that require high-precision models of the PSF.

Leakage of diffuse polarized emission into Stokes I caused by the polarized primary beam of the instrument might mimic the spectral structure of the 21-cm signal coming from the epoch of reionization (EoR) making their separation difficult. Therefore, understanding polarimetric performance of the antenna is crucial for a successful detection of the EoR signal. Here, we have calculated the accuracy of the nominal model beam of LOFAR in predicting the leakage from Stokes I to Q, U by comparing them with the corresponding leakage of compact sources actually observed in the 3C295 field. We have found that the model beam has errors of less than or equal to 10% on the predicted levels of leakage of ~1% within the field of view, i. e. if the leakage is taken out perfectly using this model the leakage will reduce to $10^{-3}$ of the Stokes I flux. If similar levels of accuracy can be obtained in removing leakage from Stokes Q, U to I, we can say, based on the results of our previous paper, that the removal of this leakage using this beam model would ensure that the leakage is well below the expected EoR signal in almost the whole instrumental k-space of the cylindrical power spectrum. We have also shown here that direction dependent calibration can remove instrumentally polarized compact sources, given an unpolarized sky model, very close to the local noise level.

Based on spectroscopy and multiband wide-field observations of the gravitationally lensed quasar HE 0435-1223, we determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for this system. We measure the under/overdensity of the line of sight towards the lens system and compare it to the average line of sight throughout the universe, determined by using the CFHTLenS as a control field. Aiming to constrain $\kappa_\mathrm{ext}$ as tightly as possible, we determine under/overdensities using various combinations of relevant informative weighing schemes for the galaxy counts, such as projected distance to the lens, redshift, and stellar mass. We then convert the measured under/overdensities into a $\kappa_\mathrm{ext}$ distribution, using ray-tracing through the Millennium Simulation. We explore several limiting magnitudes and apertures, and account for systematic and statistical uncertainties relevant to the quality of the observational data, which we further test through simulations. Our most robust estimate of $\kappa_\mathrm{ext}$ has a median value $\kappa^\mathrm{med}_\mathrm{ext} = 0.004$ and a standard deviation of $\sigma_\kappa = 0.025$. The measured $\sigma_\kappa$ corresponds to $2.5\%$ uncertainty on the time delay distance, and hence the Hubble constant $H_0$ inference from this system. The median $\kappa^\mathrm{med}_\mathrm{ext}$ value is robust to $\sim0.005$ (i.e. $\sim0.5\%$ on $H_0$) regardless of the adopted aperture radius, limiting magnitude and weighting scheme, as long as the latter incorporates galaxy number counts, the projected distance to the main lens, and a prior on the external shear obtained from mass modeling. The availability of a well-constrained $\kappa_\mathrm{ext}$ makes \hequad\ a valuable system for measuring cosmological parameters using strong gravitational lens time delays.

LOFAR is the LOw Frequency Radio interferometer ARray located at mid-latitude ($52^{\circ} 53'N$). Here, we present results on ionospheric structures derived from 29 LOFAR nighttime observations during the winters of 2012/2013 and 2013/2014. We show that LOFAR is able to determine differential ionospheric TEC values with an accuracy better than 1 mTECU over distances ranging between 1 and 100 km. For all observations the power law behavior of the phase structure function is confirmed over a long range of baseline lengths, between $1$ and $80$ km, with a slope that is in general larger than the $5/3$ expected for pure Kolmogorov turbulence. The measured average slope is $1.89$ with a one standard deviation spread of $0.1$. The diffractive scale, i.e. the length scale where the phase variance is $1\, \mathrm{rad^2}$, is shown to be an easily obtained single number that represents the ionospheric quality of a radio interferometric observation. A small diffractive scale is equivalent to high phase variability over the field of view as well as a short time coherence of the signal, which limits calibration and imaging quality. For the studied observations the diffractive scales at $150$ MHz vary between $3.5$ and $30\,$ km. A diffractive scale above $5$ km, pertinent to about $90 \%$ of the observations, is considered sufficient for the high dynamic range imaging needed for the LOFAR Epoch of Reionization project. For most nights the ionospheric irregularities were anisotropic, with the structures being aligned with the Earth magnetic field in about $60\%$ of the observations.

The redshifted 21 cm line of neutral hydrogen is a promising probe of the Epoch of Reionization (EoR). However, its detection requires a thorough understanding and control of the systematic errors. We study two systematic biases observed in the LOFAR EoR residual data after calibration and subtraction of bright discrete foreground sources. The first effect is a suppression in the diffuse foregrounds, which could potentially mean a suppression of the 21 cm signal. The second effect is an excess of noise beyond the thermal noise. The excess noise shows fluctuations on small frequency scales, and hence it can not be easily removed by foreground removal or avoidance methods. Our analysis suggests that sidelobes of residual sources due to the chromatic point spread function and ionospheric scintillation can not be the dominant causes of the excess noise. Rather, both the suppression of diffuse foregrounds and the excess noise can occur due to calibration with an incomplete sky model containing predominantly bright discrete sources. We show that calibrating only on bright sources can cause suppression of other signals and introduce an excess noise in the data. The levels of the suppression and excess noise depend on the relative flux of sources which are not included in the model with respect to the flux of modeled sources. We discuss possible solutions such as using only long baselines to calibrate the interferometric gain solutions as well as simultaneous multi-frequency calibration along with their benefits and shortcomings.

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.

Compact radio sources sometimes exhibit intervals of large, rapid changes in their flux-density, due to lensing by interstellar plasma crossing the line-of-sight. A novel survey program has made it possible to discover these "Extreme Scattering Events" (ESEs) in real time, resulting in a high-quality dynamic spectrum of an ESE observed in PKS 1939-315. Here we present a method for determining the column-density profile of a plasma lens, given only the dynamic radio spectrum of the lensed source, under the assumption that the lens is either axisymmetric or totally anisotropic. Our technique relies on the known, strong frequency dependence of the plasma refractive index in order to determine how points in the dynamic spectrum map to positions on the lens. We apply our method to high-frequency (4.2-10.8 GHz) data from the Australia Telescope Compact Array of the PKS 1939-315 ESE. The derived electron column-density profiles are very similar for the two geometries we consider, and both yield a good visual match to the data. However, the fit residuals are substantially above the noise level, and deficiencies are evident when we compare the predictions of our model to lower-frequency (1.6-3.1 GHz) data on the same ESE, thus motivating future development of more sophisticated inversion techniques.

The observed 21-cm signal from the epoch of reionization will be distorted along the line-of-sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus potentially be used to extract information about the sources of reionization. In this paper, we study a collection of simulated reionization scenarios assuming different models for the sources of reionization. We show that the 21-cm anisotropy is best measured by the quadrupole moment of the power spectrum. We find that, unless the properties of the reionization sources are extreme in some way, the quadrupole moment evolves very predictably as a function of global neutral fraction. This predictability implies that redshift-space distortions are not a very sensitive tool for distinguishing between reionization sources. However, the quadrupole moment can be used as a model-independent probe for constraining the reionization history. We show that such measurements can be done to some extent by first-generation instruments such as LOFAR, while the SKA should be able to measure the reionization history using the quadrupole moment of the power spectrum to great accuracy.

Visibility scintillation resulting from wave propagation through the turbulent ionosphere can be an important sources of noise at low radio frequencies ($\nu\lesssim 200$ MHz). Many low frequency experiments are underway to detect the power spectrum of brightness temperature fluctuations of the neutral-hydrogen $21$-cm signal from the Epoch of Reionization (EOR: $12\gtrsim z\gtrsim 7$, $100\lesssim \nu \lesssim 175$ MHz). In this paper, we derive scintillation noise power-spectra in such experiments while taking into account the effects of typical data processing operations such as self-calibration and Fourier synthesis. We find that for minimally redundant arrays such as LOFAR and MWA, scintillation noise is of the same order of magnitude as thermal noise, has a spectral coherence dictated by stretching of the snapshot $uv$-coverage with frequency, and thus is confined to the well known wedge-like structure in the cylindrical ($2$-dimensional) power spectrum space. Compact, fully redundant ($d_{\rm core}\lesssim r_{\rm F} \approx 300$ m at $150$ MHz) arrays such as HERA and SKA-LOW (core) will be scintillation noise dominated at all baselines, but the spatial and frequency coherence of this noise will allow it to be removed along with spectrally smooth foregrounds.

In this paper, we consider random phase fluctuations imposed during wave propagation through a turbulent plasma (e.g. ionosphere) as a source of additional noise in interferometric visibilities. We derive expressions for visibility variance for the wide field of view case (FOV$\sim10$ deg) by computing the statistics of Fresnel diffraction from a stochastic plasma, and provide an intuitive understanding. For typical ionospheric conditions (diffractive scale $\sim 5-20$ km at $150$ MHz), we show that the resulting ionospheric `scintillation noise' can be a dominant source of uncertainty at low frequencies ($\nu \lesssim 200$ MHz). Consequently, low frequency widefield radio interferometers must take this source of uncertainty into account in their sensitivity analysis. We also discuss the spatial, temporal, and spectral coherence properties of scintillation noise that determine its magnitude in deep integrations, and influence prospects for its mitigation via calibration or filtering.

This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. We have used the high band antennas (HBA) of LOFAR to image this region and Rotation Measure (RM) synthesis to unravel the distribution of polarized structures in Faraday depth. The brightness temperature of the detected Galactic emission is $5-15~{\rm K}$ in polarized intensity and covers the range from -3 to +8 ${\rm rad~m^{-2}}$ in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of $+0.5~{\rm rad~m^{-2}}$ running from north to south, right through the centre of the field and parallel to the Galactic plane. There is also an interesting system of linear depolarization canals conspicuous in an image showing the peaks of Faraday spectra. We used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. For the first time, we see some common morphology in the RM cubes made at 150 and 350~{\rm MHz}. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies and previous LOFAR observations. Based on our results, we determined physical parameters of the ISM and proposed a simple model that may explain the observed distribution of the intervening magneto-ionic medium. The mean line-of-sight magnetic field component, $B_\parallel$, is determined to be $0.3\pm0.1~{\rm \mu G}$ and its spatial variation across the 3C196 field is $0.1~{\rm \mu G}$. The filamentary structure is probably an ionized filament in the ISM, located somewhere within the Local Bubble. This filamentary structure shows an excess in thermal electron density ($n_e B_\parallel>6.2~{\rm cm^{-3}\mu G}$) compared to its surroundings.

The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru's Hyper Suprime Cam (HSC), we investigate the observability of such cross-power spectrum with these two instruments, which are both planning to observe the ELAIS-N1 field at z=6.6. In this paper we use N-body + radiative transfer (both for continuum and Lyalpha photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm-galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power spectrum expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum show clear anti-correlation on scales larger than ~ 60 h$^{-1}$ Mpc (corresponding to k ~ 0.1 h Mpc$^{-1}$), with levels of significance p=0.04 at z=6.6 and p=0.048 at z=7.3. On smaller scales, instead, the signal is completely contaminated.

With the first phase of the Square Kilometre Array (SKA1) entering into its final pre-construction phase, we investigate how best to maximise its scientific return. Specifically, we focus on the statistical measurement of the 21 cm power spectrum (PS) from the epoch of reionization (EoR) using the low frequency array, SKA1-low. To facilitate this investigation we use the recently developed MCMC based EoR analysis tool 21CMMC (Greig & Mesinger). In light of the recent 50 per cent cost reduction, we consider several different SKA core baseline designs, changing: (i) the number of antenna stations; (ii) the number of dipoles per station; and also (iii) the distribution of baseline lengths. We find that a design with a reduced number of dipoles per core station (increased field of view and total number of core stations), together with shortened baselines, maximises the recovered EoR signal. With this optimal baseline design, we investigate three observing strategies, analysing the trade-off between lowering the instrumental thermal noise against increasing the field of view. SKA1-low intends to perform a three tiered observing approach, including a deep 100 deg$^{2}$ at 1000 h, a medium-deep 1000 deg$^{2}$ at 100 h and a shallow 10,000 deg$^{2}$ at 10 h survey. We find that the three observing strategies result in comparable ($\lesssim$ per cent) constraints on our EoR astrophysical parameters. This is contrary to naive predictions based purely on the total signal-to-noise, thus highlighting the need to use EoR parameter constraints as a figure of merit, in order to maximise scientific returns with next generation interferometers.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structures along the line of sight, whereas the Planck dust polarization maps probe the orientation of the sky projected magnetic field component. Hence, no apparent correlation between the two is expected. Here we report a surprising, yet clear, correlation between the filamentary structures, detected with LOFAR, and the magnetic field orientation, probed by the Planck satellite. This finding points to a common, yet unclear, physical origin of the two measurements in this specific area in the sky. A number of follow-up multi- frequency studies are proposed to shed light on this unexpected finding.

The Kilo-Degree Survey (KiDS) is an optical wide-field imaging survey carried out with the VLT Survey Telescope and the OmegaCAM camera. KiDS will image 1500 square degrees in four filters (ugri), and together with its near-infrared counterpart VIKING will produce deep photometry in nine bands. Designed for weak lensing shape and photometric redshift measurements, the core science driver of the survey is mapping the large-scale matter distribution in the Universe back to a redshift of ~0.5. Secondary science cases are manifold, covering topics such as galaxy evolution, Milky Way structure, and the detection of high-redshift clusters and quasars. KiDS is an ESO Public Survey and dedicated to serving the astronomical community with high-quality data products derived from the survey data, as well as with calibration data. Public data releases will be made on a yearly basis, the first two of which are presented here. For a total of 148 survey tiles (~160 sq.deg.) astrometrically and photometrically calibrated, coadded ugri images have been released, accompanied by weight maps, masks, source lists, and a multi-band source catalog. A dedicated pipeline and data management system based on the Astro-WISE software system, combined with newly developed masking and source classification software, is used for the data production of the data products described here. The achieved data quality and early science projects based on the data products in the first two data releases are reviewed in order to validate the survey data. Early scientific results include the detection of nine high-z QSOs, fifteen candidate strong gravitational lenses, high-quality photometric redshifts and galaxy structural parameters for hundreds of thousands of galaxies. (Abridged)

We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.

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.

We have started a systematic search of strong lens candidates in the ESO public survey KiDS based on the visual inspection of massive galaxies in the redshift range $0.1<z<0.5$. As a pilot program we have inspected 100 sq. deg., which overlap with SDSS and where there are known lenses to use as a control sample. Taking advantage of the superb image quality of VST/OmegaCAM, the colour information and accurate model subtracted images, we have found 18 new lens candidates, for which spectroscopic confirmation will be needed to confirm their lensing nature and study the mass profile of the lensing galaxies.

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 $2\sigma$ 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.

Experiments designed to measure the redshifted 21~cm line from the Epoch of Reionization (EoR) are challenged by strong astrophysical foreground contamination, ionospheric distortions, complex instrumental response and other different types of noise (e.g. radio frequency interference). The astrophysical foregrounds are dominated by diffuse synchrotron emission from our Galaxy. Here we present a simulation of the Galactic emission used as a foreground module for the LOFAR- EoR key science project end-to-end simulations. The simulation produces total and polarized intensity over $10^\circ \times 10^\circ$ maps of the Galactic synchrotron and free-free emission, including all observed characteristics of the emission: spatial fluctuations of amplitude and spectral index of the synchrotron emission, together with Faraday rotation effects. The importance of these simulations arise from the fact that the Galactic polarized emission could behave in a manner similar to the EoR signal along the frequency direction. As a consequence, an improper instrumental calibration will give rise to leakages of the polarized to the total signal and mask the desired EoR signal. In this paper we address this for the first time through realistic simulations.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

We conduct a detailed lensing, dynamics and stellar population analysis of nine massive lens early-type galaxies (ETGs) from the X-Shooter Lens Survey (XLENS). Combining gravitational lensing constraints from HST imaging with spatially-resolved kinematics and line-indices constraints from VLT X-Shooter (XSH) spectra, we infer the low-mass slope and the low cut-off mass of the stellar Initial Mass Function (IMF): $x_{250}=2.37^{+0.12}_{-0.12}$ and $M_{{\rm low}, 250}= 0.131^{+0.023}_{-0.026}\, M_{\odot}$, respectively, for a reference point with $\sigma \equiv 250\, {{\rm kms}}^{-1}$ and R$_{{\rm eff}} \equiv 10$ kpc. All the XLENS systems are consistent with an IMF slope steeper than Milky Way-like. We find no significant correlations between IMF slope and any other quantity, except for an anti-correlation between total dynamical mass density and low-mass IMF slope at the 87% CL [$dx/d\log(\rho)$ = $ -0.19^{+0.15}_{-0.15}$]. This anti-correlation is consistent with the low redshift lenses found by Smith et al. (2015) that have high velocity dispersions and high stellar mass densities but surprisingly shallow IMF slopes.

Detection of the 21-cm signal coming from the epoch of reionization (EoR) is challenging especially because, even after removing the foregrounds, the residual Stokes $I$ maps contain leakage from polarized emission that can mimic the signal. Here, we discuss the instrumental polarization of LOFAR and present realistic simulations of the leakages between Stokes parameters. From the LOFAR observations of polarized emission in the 3C196 field, we have quantified the level of polarization leakage caused by the nominal model beam of LOFAR, and compared it with the EoR signal using power spectrum analysis. We found that at 134--166 MHz, within the central 4$^\circ$ of the field the $(Q,U)\rightarrow I$ leakage power is lower than the EoR signal at $k<0.3$ Mpc$^{-1}$. The leakage was found to be localized around a Faraday depth of 0, and the rms of the leakage as a fraction of the rms of the polarized emission was shown to vary between 0.2-0.3\%, both of which could be utilized in the removal of leakage. Moreover, we could define an `EoR window' in terms of the polarization leakage in the cylindrical power spectrum above the PSF-induced wedge and below $k_\parallel\sim 0.5$ Mpc$^{-1}$, and the window extended up to $k_\parallel\sim 1$ Mpc$^{-1}$ at all $k_\perp$ when 70\% of the leakage had been removed. These LOFAR results show that even a modest polarimetric calibration over a field of view of $\lesssim 4^\circ$ in the future arrays like SKA will ensure that the polarization leakage remains well below the expected EoR signal at the scales of 0.02-1 Mpc$^{-1}$.

We perform a direct comparison of two state-of-the art single stellar population (SSP) models that have been used to demonstrate the non-universality of the low-mass end of the Initial Mass Function (IMF) slope. The two public versions of the SSP models are restricted to either solar abundance patterns or solar metallicity, too restrictive if one aims to disentangle elemental enhancements, metallicity changes and IMF variations in massive early-type galaxies (ETGs) with star formation histories different from the solar neighborhood. We define response functions (to metallicity and \alpha-abundance) to extend the parameter space of each set of models. We compare these extended models with a sample of Sloan Digital Sky Survey (SDSS) ETGs spectra with varying velocity dispersions. We measure equivalent widths of optical IMF-sensitive stellar features to examine the effect of the underlying model assumptions and ingredients, such as stellar libraries or isochrones, on the inference of the IMF slope down to ~0.1 solar masses. We demonstrate that the steepening of the low-mass end of the Initial Mass Function (IMF) based on a non-degenerate set of spectroscopic optical indicators is robust against the choice of the stellar population model. Although the models agree in a relative sense (i.e. both imply more bottom-heavy IMFs for more massive systems), we find non-negligible differences on the absolute values of the IMF slope inferred at each velocity dispersion by using the two different models. In particular, we find large inconsistency in the quantitative predictions of IMF slope variations and abundance patterns when sodium lines are used. We investigate the possible reasons for these inconsistencies.

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.

We present observational results from a new Hubble Space Telescope (HST) Snapshot program to extend the methods of the Sloan Lens ACS (SLACS) Survey to lower lens-galaxy masses. We discover 40 new galaxy-scale strong lenses, which we supplement with 58 previously discovered SLACS lenses. In addition, we determine the posterior PDFs of the Einstein radius for 33 galaxies (18 new and 15 from legacy SLACS data) based on single lensed images. We find a less-than-unity slope of $0.64\pm0.06$ for the $\log_{10} {\sigma}_*$-$\log_{10} {\sigma}_{\rm SIE}$ relation, which corresponds to a 6-$\sigma$ evidence that the total mass-density profile of early-type galaxies varies systematically in the sense of being shallower at higher lens-galaxy velocity dispersions. The trend is only significant when single-image systems are considered, highlighting the importance of including both "lenses" and "non-lenses" for an unbiased treatment of the lens population when extending to lower mass ranges. By scaling simple stellar population models to the HST I-band data, we identify a strong trend of increasing dark-matter fraction at higher velocity dispersions, which can be alternatively interpreted as a trend in the stellar initial mass function (IMF) normalization. Consistent with previous findings and the suggestion of a non-universal IMF, we find that a Salpeter IMF is ruled out for galaxies with velocity dispersion less than $180$ km/s. Considered together, our mass-profile and dark-matter-fraction trends with increasing galaxy mass could both be explained by an increasing relative contribution on kiloparsec scales from a dark-matter halo with a spatial profile more extended than that of the stellar component.

We provide an overview of 21cm tomography of the Cosmic Dawn and Epoch of Reionization as possible with SKA-Low. We show why tomography is essential for studying CD/EoR and present the scales which can be imaged at different frequencies for the different phases of SKA- Low. Next we discuss the different ways in which tomographic data can be analyzed. We end with an overview of science questions which can only be answered by tomography, ranging from the characterization of individual objects to understanding the global processes shaping the Universe during the CD/EoR

The early Universe is a precious probe of the birth of primordial objects, first star formation events and consequent production of photons and heavy elements. Higher-order corrections to the cosmological linear perturbation theory predicts the formation of coherent supersonic gaseous streaming motions at decoupling time. These bulk flows impact the gas cooling process and determine a cascade effect on the whole baryon evolution. By analytical estimates and N-body hydrodynamical chemistry numerical simulations including atomic and molecular evolution, gas cooling, star formation, feedback effects and metal spreading for individual species from different stellar populations according to the proper yields and lifetimes, we discuss the role of these primordial bulk flows at the end of the dark ages and their detectable impacts during the first Gyr in view of the upcoming SKA mission. Early bulk flows can inhibit molecular gas cooling capabilities, suppressing star formation, metal spreading and the abundance of small primordial galaxies in the infant Universe. This can determine a delay in the re-ionization process and in the heating of neutral hydrogen making the observable HI signal during cosmic evolution patchier and noisier. The planned SKA mission will represent a major advance over existing instruments, since it will be able to probe the effects on HI 21cm at z ~ 6-20 and on molecular line emissions from first collapsing sites at z ~ 20-40. Therefore, it will be optimal to address the effects of primordial streaming motions on early baryon evolution and to give constraints on structure formation in the first Gyr.

We observed six strongly lensed, radio-loud quasars (MG 0414+0534, CLASS B0712+472, JVAS B1030+074, CLASS B1127+385, CLASS B1152+199, and JVAS B1938+666) in order to identify systems suitable for measuring cosmological parameters using time delays between their multiple images. These systems are in standard two and four image configurations, with B1938 having a faint secondary pair of images. Two separate monitoring campaigns were carried out using the VLA and upgraded JVLA. Lightcurves were extracted for each individual lensed image and analyzed for signs of intrinsic variability. While it was not possible to measure time delays from these data, $\chi^2$-based and structure function tests found evidence for variability in a majority of the lightcurves. B0712 and B1030 had particularly strong variations, exhibiting linear flux trends. These results show that most of these systems should be targeted with followup monitoring campaigns, especially B0712 and B1030.

The combination of strong gravitational lensing and stellar kinematics provides a powerful and robust method to investigate the mass and dynamical structure of early-type galaxies. We demonstrate this approach by analysing two massive ellipticals from the XLENS Survey for which both high-resolution HST imaging and X-Shooter spectroscopic observations are available. We adopt a flexible axisymmetric two-component mass model for the lens galaxies, consisting of a generalised NFW dark halo and a realistic self-gravitating stellar mass distribution. For both systems, we put constraints on the dark halo inner structure and flattening, and we find that they are dominated by the luminous component within one effective radius. By comparing the tight inferences on the stellar mass from the combined lensing and dynamics analysis with the values obtained from stellar population studies, we conclude that both galaxies are characterised by a Salpeter-like stellar initial mass function.

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<\nu<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the global history of the EoR. Signal variance is the integral of the signal's power spectrum, and it is expected to be measured with a high significance. We demonstrate this through results from a simulation and parameter estimation pipeline developed for the Low Frequency Array (LOFAR)-EoR experiment. We show that LOFAR should be able to detect the EoR in 600 hours of integration using the variance statistic. Additionally, the redshift ($z_r$) and duration ($\Delta z$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $\Delta z = 0.43$ to test the pipeline. We are able to detect the simulated signal with a significance of 4 standard deviations and extract the EoR parameters as $z_r = 7.72^{+0.37}_{-0.18}$ and $\Delta z = 0.53^{+0.12}_{-0.23}$ in 600 hours, assuming that systematic errors can be adequately controlled. We further show that the significance of detection and constraints on EoR parameters can be improved by measuring the cross-variance of the signal by cross-correlating consecutive redshift bins.

Herschel Space Observatory photometry and extensive multiwavelength followup have revealed that the powerful radio galaxy 3C 220.3 at z=0.685 acts as a gravitational lens for a background submillimeter galaxy (SMG) at z=2.221. At an observed wavelength of 1mm, the SMG is lensed into three distinct images. In the observed near infrared, these images are connected by an arc of 1.8" radius forming an Einstein half-ring centered near the radio galaxy. In visible light, only the arc is apparent. 3C 220.3 is the only known instance of strong galaxy-scale lensing by a powerful radio galaxy not located in a galaxy cluster and therefore it offers the potential to probe the dark matter content of the radio galaxy host. Lens modeling rejects a single lens, but two lenses centered on the radio galaxy host A and a companion B, separated by 1.5", provide a fit consistent with all data and reveal faint candidates for the predicted fourth and fifth images. The model does not require an extended common dark matter halo, consistent with the absence of extended bright X-ray emission on our Chandra image. The projected dark matter fractions within the Einstein radii of A (1.02") and B (0.61") are about 0.4 +/- 0.3 and 0.55 +/- 0.3. The mass to i-band light ratios of A and B, M/L ~ 8 +/- 4 Msun/Lsun, appear comparable to those of radio-quiet lensing galaxies at the same redshift in the CASTLES, LSD, and SL2S samples. The lensed SMG is extremely bright with observed f(250um) = 440mJy owing to a magnification factor mu~10. The SMG spectrum shows luminous, narrow CIV 154.9nm emission, revealing that the SMG houses a hidden quasar in addition to a violent starburst. Multicolor image reconstruction of the SMG indicates a bipolar morphology of the emitted ultraviolet (UV) light suggestive of cones through which UV light escapes a dust-enshrouded nucleus.

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 $\alpha$ < 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 $\alpha$ = 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.

Under the assumption of a flat Lambda-CDM cosmology, recent data from the Planck satellite point toward a Hubble constant that is in tension with that measured by gravitational lens time delays and by the local distance ladder. Prosaically, this difference could arise from unknown systematic uncertainties in some of the measurements. More interestingly -- if systematics were ruled out -- resolving the tension would require a departure from the flat Lambda-CDM cosmology, introducing for example a modest amount of spatial curvature, or a non-trivial dark energy equation of state. To begin to address these issues, we present here an analysis of the gravitational lens RXJ1131-1231 that is improved in one particular regard: we examine the issue of systematic error introduced by an assumed lens model density profile. We use more flexible gravitational lens models with baryonic and dark matter components, and find that the exquisite Hubble Space Telescope image with thousands of intensity pixels in the Einstein ring and the stellar velocity dispersion of the lens contain sufficient information to constrain these more flexible models. The total uncertainty on the time-delay distance is 6.6% for a single system. We proceed to combine our improved time-delay distance measurements with the WMAP9 and Planck posteriors. In an open Lambda-CDM model, the data for RXJ1131-1231 in combination with Planck favor a flat universe with Omega_k = 0.00+0.01/-0.02 (68% CI). In a flat wCDM model, the combination of RXJ1131-1231 and Planck yields w = -1.52+0.19/-0.20 (68% CI).

Galaxy masses play a fundamental role in our understanding of structure formation models. This review addresses the variety and reliability of mass estimators that pertain to stars, gas, and dark matter. The different sections on masses from stellar populations, dynamical masses of gas-rich and gas-poor galaxies, with some attention paid to our Milky Way, and masses from weak and strong lensing methods, all provide review material on galaxy masses in a self-consistent manner.

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.

We present an algorithm using Principal Component Analysis (PCA) to subtract galaxies from imaging data, and also two algorithms to find strong, galaxy-scale gravitational lenses in the resulting residual image. The combined method is optimized to find full or partial Einstein rings. Starting from a pre-selection of potential massive galaxies, we first perform a PCA to build a set of basis vectors. The galaxy images are reconstructed using the PCA basis and subtracted from the data. We then filter the residual image with two different methods. The first uses a curvelet (curved wavelets) filter of the residual images to enhance any curved/ring feature. The resulting image is transformed in polar coordinates, centered on the lens galaxy center. In these coordinates, a ring is turned into a line, allowing us to detect very faint rings by taking advantage of the integrated signal-to-noise in the ring (a line in polar coordinates). The second way of analysing the PCA-subtracted images identifies structures in the residual images and assesses whether they are lensed images according to their orientation, multiplicity and elongation. We apply the two methods to a sample of simulated Einstein rings, as they would be observed with the ESA Euclid satellite in the VIS band. The polar coordinates transform allows us to reach a completeness of 90% and a purity of 86%, as soon as the signal-to-noise integrated in the ring is higher than 30, and almost independent of the size of the Einstein ring. Finally, we show with real data that our PCA-based galaxy subtraction scheme performs better than traditional subtraction based on model fitting to the data. Our algorithm can be developed and improved further using machine learning and dictionary learning methods, which would extend the capabilities of the method to more complex and diverse galaxy shapes.

We investigate the optical spectral region of spectra of 1000 stars searching for IMF-sensitive features to constrain the low-mass end of the initial mass function (IMF) slope in elliptical galaxies. We use the MILES stellar library in the wavelength range [3500-7500] {\AA} to select indices that are sensitive to cool dwarf stars and that only weakly depend on age and metallicity. In this wavelength range, the response of a change in the effective temperature of the cool red giant population is similar to the response of a change in the number of dwarf stars in the galaxy. We investigate the degeneracy between IMF variation and Teff,RGB and show that it is possible to break it with the new IMF indicators defined here. In particular, we define a CaH1 index around {\lambda}6380 {\AA} that arises purely from cool dwarfs, does not strongly depend on age and is anti-correlated with [a/Fe]. This index allows the determination of the the IMF slope when combined with TiO lines and age- and metallicity-dependent features. The use of several indicators is crucial to break degeneracies between IMF, age, abundance pattern and Teff of the RGB population. We measure line-index strengths of the optical IMF indicators in the SSP models and compare these with index strengths of the same features in stacked SSDSS early-type galaxy spectra with varying velocity dispersions. Using 10 different indicators, we find a clear trend of a steepening IMF with increasing velocity dispersion described by the equation x=(2.3+-0.1)*log{\sigma}200+(2.13+-0.15), where x is the IMF slope and {\sigma}200 is the central stellar velocity dispersion measured in units of 200km/s. We compute the 'IMF mismatch' parameter as the ratio of stellar mass-to-light ratio predicted from the x-{\sigma}200 relation to that inferred from SSP models assuming a Salpeter IMF and find good agreement with independent published results.

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

We conduct a detailed investigation of the properties of the stellar initial mass function (IMF) in two massive early-type lens galaxies with velocity dispersions of sigma ~245 km/s and sigma ~325 km/s, for which both HST imaging and X-Shooter spectra are available. We compare the inferences obtained from two fully independent methods: (i) a combined gravitational lensing and stellar dynamics (L&D) analysis of the data sets employing self-consistent axisymmetric models, and (ii) a spectroscopic simple stellar population (SSP) analysis of optical line-strength indices, assuming single power-law IMFs. The results from the two approaches are found to be in agreement within the 1-sigma uncertainties. Both galaxies are consistent with having a Salpeter IMF (power-law slope of x = 2.35), which is strongly favoured over a Chabrier IMF (x = 1.8), with probabilities inferred from the joint analysis of 89% and 99%, respectively. Bottom-heavy IMFs significantly steeper than Salpeter (x >= 3.0) are ruled out with decisive evidence (Bayes factor B > 1000) for both galaxies, as they exceed the total mass derived from the L&D constraints. Our analysis allows, for the first time, the inference of the low-mass cut-off of the IMF (M_low). Combining the joint L&D and SSP analyses of both galaxies, we infer an IMF slope of x = 2.22 +/- 0.14, consistent with Salpeter IMF, and a low-mass limit M_low = 0.13 +/- 0.03 M_sun, just above the hydrogen burning limit.

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios.

The redshifted 21 cm brightness distribution from neutral hydrogen is a promising probe into the cosmic dark ages, cosmic dawn, and re-ionization. LOFAR's Low Band Antennas (LBA) may be used in the frequency range 45 MHz to 85 MHz (30>z>16) to measure the sky averaged redshifted 21 cm brightness temperature as a function of frequency, or equivalently, cosmic redshift. These low frequencies are affected by strong Galactic foreground emission that is observed through frequency dependent ionospheric and antenna beam distortions which lead to chromatic mixing of spatial structure into spectral structure. Using simple models, we show that (i) the additional antenna temperature due to ionospheric refraction and absorption are at a \sim 1% level--- 2 to 3 orders of magnitude higher than the expected 21 cm signal, and have an approximate \nu^{-2} dependence, (ii) ionospheric refraction leads to a knee-like modulation on the sky spectrum at \nu\approx 4\times plasma frequency. Using more realistic simulations, we show that in the measured sky spectrum, more than 50% of the 21 cm signal variance can be lost to confusion from foregrounds and chromatic effects. We conclude that foregrounds and chromatic mixing may not be subtracted as generic functions of frequency as previously thought, but must rather be carefully modeled using additional priors and interferometric measurements.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

Using a combination of N-body simulations, semi-analytic models and radiative transfer calculations, we have estimated the theoretical cross power spectrum between galaxies and the 21cm emission from neutral hydrogen during the epoch of reionization. In accordance with previous studies, we find that the 21cm emission is initially correlated with halos on large scales (> 30 Mpc), anti-correlated on intermediate (~ 5 Mpc), and uncorrelated on small (< 3 Mpc) scales. This picture quickly changes as reionization proceeds and the two fields become anti-correlated on large scales. The normalization of the cross power spectrum can be used to set constraints on the average neutral fraction in the intergalactic medium and its shape can be a tool to study the topology of reionization. When we apply a drop-out technique to select galaxies and add to the 21cm signal the noise expected from the LOFAR telescope, we find that while the normalization of the cross power spectrum remains a useful tool for probing reionization, its shape becomes too noisy to be informative. On the other hand, for a Lyalpha Emitter (LAE) survey both the normalization and the shape of the cross power spectrum are suitable probes of reionization. A closer look at a specific planned LAE observing program using Subaru Hyper-Suprime Cam reveals concerns about the strength of the 21cm signal at the planned redshifts. If the ionized fraction at z ~ 7 is lower that the one estimated here, then using the cross power spectrum may be a useful exercise given that at higher redshifts and neutral fractions it is able to distinguish between two toy models with different topologies.

High precision cosmological distance measurements towards individual objects such as time delay gravitational lenses or type Ia supernovae are affected by weak lensing perturbations by galaxies and groups along the line of sight. In time delay gravitational lenses, "external convergence," kappa, can dominate the uncertainty in the inferred distances and hence cosmological parameters. In this paper we attempt to reconstruct kappa, due to line of sight structure, using a simple halo model. We use mock catalogues from the Millennium Simulation, and calibrate and compare our reconstructed P(kappa) to ray-traced kappa "truth" values; taking into account realistic observational uncertainties. We find that the reconstruction of kappa provides an improvement in precision of ~50% over galaxy number counts. We find that the lowest-kappa lines of sight have the best constrained P(kappa). In anticipation of large future samples of lenses, we find that selecting the third of the systems with the highest precision kappa estimates gives a sample of unbiased time delay distance measurements with just ~1% uncertainty due to line of sight external convergence effects. Photometric data are sufficient to pre-select the best-constrained lines of sight, and can be done before investment in light-curve monitoring. Conversely, we show that selecting lines of sight with high external shear could, with the reconstruction model presented, induce biases of up to 1% in time delay distance. We find that a major potential source of systematic error is uncertainty in the high mass end of the stellar mass-halo mass relation; this could introduce ~2% biases on the time-delay distance if completely ignored. We suggest areas for the improvement of this general analysis framework (including more sophisticated treatment of high mass structures) that should allow yet more accurate cosmological inferences to be made.

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line of sight convergence k_ext. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with k_ext obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of k_ext from sigma_k $0.06 to ~0.04 for very overdense lines of sight like that of the system B1608+656. For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to an uncertainty of ~<3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain k_ext almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from 0.065 using only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

Strong gravitational lenses with measured time delays between the multiple images and models of the lens mass distribution allow a one-step determination of the time-delay distance, and thus a measure of cosmological parameters. We present a blind analysis of the gravitational lens RXJ1131-1231 incorporating (1) the newly measured time delays from COSMOGRAIL, (2) archival HST imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323+/-20km/s based on Keck spectroscopy, and (4) a characterization of the line-of-sight structures via observations of the lens' environment and ray tracing through the Millennium Simulation. Our blind analysis is designed to prevent experimenter bias. The joint analysis of the data sets allows a time-delay distance measurement to 6% precision that takes into account all known systematic uncertainties. In combination with the WMAP7 data set in flat wCDM cosmology, our unblinded cosmological constraints for RXJ1131-1231 are: H_0=80.0+5.8/-5.7km/s/Mpc, OmegaDE=0.79+/-0.03 and w=-1.25+0.17/-0.21. We find the results to be statistically consistent with those from the analysis of the gravitational lens B1608+656. The joint constraints from the two lenses and WMAP7 are H_0=75.2+4.4/-4.2km/s/Mpc, OmegaDE=0.76+0.02/-0.03 and w=-1.14+0.17/-0.20 in flat wCDM, and H_0=73.1+2.4/-3.6km/s/Mpc, OmegaL=0.75+0.01/-0.02 and OmegaK=0.003+0.005/-0.006 in open LCDM. Time-delay lenses constrain especially tightly the Hubble constant (5.7% and 4.0% respectively in wCDM and open LCDM) and curvature of the universe. They complement well other cosmological probes, and provide an independent check of unknown systematics. Our measurement of the Hubble constant is completely independent of those based on the local distance ladder method, providing an important consistency check of the standard cosmological model and of general relativity.

One of the key science projects of the Low-Frequency Array (LOFAR) is the detection of the cosmological signal coming from the Epoch of Reionization (EoR). Here we present the LOFAR EoR Diagnostic Database (LEDDB) that is used in the storage, management, processing and analysis of the LOFAR EoR observations. It stores referencing information of the observations and diagnostic parameters extracted from their calibration. This stored data is used to ease the pipeline processing, monitor the performance of the telescope and visualize the diagnostic parameters which facilitates the analysis of the several contamination effects on the signals. It is implemented with PostgreSQL and accessed through the psycopg2 python module. We have developed a very flexible query engine, which is used by a web user interface to access the database, and a very extensive set of tools for the visualization of the diagnostic parameters through all their multiple dimensions.

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

The accurate and precise removal of 21-cm foregrounds from Epoch of Reionization redshifted 21-cm emission data is essential if we are to gain insight into an unexplored cosmological era. We apply a non-parametric technique, Generalized Morphological Component Analysis or GMCA, to simulated LOFAR-EoR data and show that it has the ability to clean the foregrounds with high accuracy. We recover the 21-cm 1D, 2D and 3D power spectra with high accuracy across an impressive range of frequencies and scales. We show that GMCA preserves the 21-cm phase information, especially when the smallest spatial scale data is discarded. While it has been shown that LOFAR-EoR image recovery is theoretically possible using image smoothing, we add that wavelet decomposition is an efficient way of recovering 21-cm signal maps to the same or greater order of accuracy with more flexibility. By comparing the GMCA output residual maps (equal to the noise, 21-cm signal and any foreground fitting errors) with the 21-cm maps at one frequency and discarding the smaller wavelet scale information, we find a correlation coefficient of 0.689, compared to 0.588 for the equivalently smoothed image. Considering only the central 50% of the maps, these coefficients improve to 0.905 and 0.605 respectively and we conclude that wavelet decomposition is a significantly more powerful method to denoise reconstructed 21-cm maps than smoothing.

Recent work has suggested that the stellar initial mass function (IMF) is not universal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or morphological type. In this paper, we investigate variations of the IMF within individual galaxies. For this purpose, we use strong lensing and gas kinematics to measure independently the normalisation of the IMF of the bulge and disk components of a sample of 5 massive spiral galaxies with substantial bulge components taken from the SWELLS survey. We find that the stellar mass of the bulges are tightly constrained by the lensing and kinematic data. A comparison with masses based on stellar population synthesis models fitted to optical and near infrared photometry favors a Salpeter-like normalisation of the IMF. Conversely, the disk masses are less well constrained due to degeneracies with the dark matter halo, but are consistent with Milky Way type IMFs in agreement with previous studies. The disks are submaximal at 2.2 disk scale lengths, but due to the contribution of the bulges, the galaxies are baryon dominated at 2.2 disk scale lengths. Globally, our inferred IMF normalisation is consistent with that found for early-type galaxies of comparable stellar mass (> 10^11 M_sun). Our results suggest a non-universal IMF within the different components of spiral galaxies, adding to the well-known differences in stellar populations between disks and bulges.

We discuss the feasibility of the detection of the 21cm forest in the diffuse IGM with the radio telescope LOFAR. The optical depth to the 21cm line has been derived using simulations of reionization which include detailed radiative transfer of ionizing photons. We find that the spectra from reionization models with similar total comoving hydrogen ionizing emissivity but different frequency distribution look remarkably similar. Thus, unless the reionization histories are very different from each other (e.g. a predominance of UV vs. x-ray heating) we do not expect to distinguish them by means of observations of the 21cm forest. Because the presence of a strong x-ray background would make the detection of 21cm line absorption impossible, the lack of absorption could be used as a probe of the presence/intensity of the x-ray background and the thermal history of the universe. Along a random line of sight LOFAR could detect a global suppression of the spectrum from z>12, when the IGM is still mostly neutral and cold, in contrast with the more well-defined, albeit broad, absorption features visible at lower redshift. Sharp, strong absorption features associated with rare, high density pockets of gas could be detected also at z~7 along preferential lines of sight.

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (\sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it hard to disentangle the two on such large scales. The poor sensitivity on small scales on the one hand, and the foregrounds effect on large scales on the other hand, make direct imaging of the Epoch of Reionization of the Universe very difficult, and detection of the signal therefore is expected to be statistical.Despite these hurdles, in this paper we argue that for many reionization scenarios low resolution images could be obtained from the expected data. This is because at the later stages of the process one still finds very large pockets of neutral regions in the IGM, reflecting the clustering of the large-scale structure, which stays strong up to scales of \sim 120 comoving Mpc/h (\sim 1 degree). The coherence of the emission on those scales allows us to reach sufficient S/N (\sim 3) so as to obtain reionization 21 cm images. Such images will be extremely valuable for answering many cosmological questions but above all they will be a very powerful tool to test our control of the systematics in the data. The existence of this typical scale (\sim 120 comoving Mpc/h) also argues for designing future EoR experiments, e.g., with SKA, with a field of view of at least 4 degree.

We measure equivalent widths (EW) - focussing on two unique features (NaI and TiO2) of low-mass stars (<0.3M\odot) - for luminous red galaxy spectra from the the Sloan Digital Sky Survey (SDSS) and X-Shooter Lens Survey (XLENS) in order to study the low-mass end of the initial mass function (IMF). We compare these EWs to those derived from simple stellar population models computed with different IMFs, ages, [{\alpha}/Fe], and elemental abundances. We find that models are able to simultaneously reproduce the observed NaD {\lambda}5895 and Na I {\lambda}8190 features for lower-mass (\sim {\sigma}\ast) early-type galaxies (ETGs) but deviate increasingly for more massive ETGs, due do strongly mismatching NaD EWs. The TiO2 {\lambda}6230 and the Na I {\lambda}8190 features together appear to be a powerful IMF diagnostic, with age and metallicity effects orthogonal to the effect of IMF. We find that both features correlate strongly with galaxy velocity dispersion. The XLENS ETG (SDSSJ0912+0029) and an SDSS ETG (SDSSJ0041-0914) appear to require both an extreme dwarf-rich IMF and a high sodium enhancement ([Na/Fe] = +0.4). In addition, lensing constraints on the total mass of the XLENS system within its Einstein radius limit a bottom-heavy IMF with a power-law slope to x \leq 3.0 at the 90% C.L. We conclude that NaI and TiO features, in comparison with state-of-the-art SSP models, suggest a mildly steepening IMF from Salpeter (dn/dm \propto m-x with x = 2.35) to x \approx 3.0 for ETGs in the range {\sigma} = 200 - 335 km s-1.

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.

Dispersion in the interstellar medium is a well known phenomenon that follows a simple relationship, which has been used to predict the time delay of dispersed radio pulses since the late 1960s. We performed wide-band simultaneous observations of four pulsars with LOFAR (at 40-190 MHz), the 76-m Lovell Telescope (at 1400 MHz) and the Effelsberg 100-m Telescope (at 8000 MHz) to test the accuracy of the dispersion law over a broad frequency range. In this paper we present the results of these observations which show that the dispersion law is accurate to better than 1 part in 100000 across our observing band. We use this fact to constrain some of the properties of the ISM along the line-of-sight and use the lack of any aberration or retardation effects to determine upper limits on emission heights in the pulsar magnetosphere. We also discuss the effect of pulse profile evolution on our observations, and the implications that it could have for precision pulsar timing projects such as the detection of gravitational waves with pulsar timing arrays.

We introduce a new implementation of the FastICA algorithm on simulated LOFAR EoR data with the aim of accurately removing the foregrounds and extracting the 21-cm reionization signal. We find that the method successfully removes the foregrounds with an average fitting error of 0.5 per cent and that the 2D and 3D power spectra are recovered across the frequency range. We find that for scales above several PSF scales the 21-cm variance is successfully recovered though there is evidence of noise leakage into the reconstructed foreground components. We find that this blind independent component analysis technique provides encouraging results without the danger of prior foreground assumptions.

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

Stars and dark matter account for most of the mass of early-type galaxies, but uncertainties in the stellar population and the dark matter profile make it challenging to distinguish between the two components. Nevertheless, precise observations of stellar and dark matter are extremely valuable for testing the many models of structure formation and evolution. We present a measurement of the stellar mass and inner slope of the dark matter halo of a massive early-type galaxy at $z=0.222$. The galaxy is the foreground deflector of the double Einstein ring gravitational lens system SDSSJ0946+1006, also known as the "Jackpot". By combining the tools of lensing and dynamics, we first constrain the mean slope of the total mass density profile ($\rho_{\rm{tot}}\propto r^{-\gamma'}$) within the radius of the outer ring to be $\gamma' = 1.98\pm0.02\pm0.01$. Then we obtain a bulge-halo decomposition, assuming a power-law form for the dark matter halo. Our analysis yields $\gamma_{\rm{DM}} = 1.7\pm0.2$ for the inner slope of the dark matter profile, in agreement with theoretical findings on the distribution of dark matter in ellipticals, and a stellar mass from lensing and dynamics $M_*^{\rm{LD}} = 5.5_{-1.3}^{+0.4}\times10^{11}M_\Sun$. By comparing this measurement with stellar masses inferred from stellar population synthesis fitting we find that a Salpeter IMF provides a good description of the stellar population of the lens while the probability of the IMF being heavier than Chabrier is 95%. Our data suggest that growth by accretion of small systems from a compact red nugget is a plausible formation scenario for this object.

Measuring the relative mass contributions of luminous and dark matter in spiral galaxies is important for understanding their formation and evolution. The combination of a galaxy rotation curve and strong lensing is a powerful way to break the disk-halo degeneracy that is inherent in each of the methods individually. We present an analysis of the 10-image radio spiral lens B1933+503 at z_l=0.755, incorporating (1) new global VLBI observations, (2) new adaptive-optics assisted K-band imaging, (3) new spectroscopic observations for the lens galaxy rotation curve and the source redshift. We construct a three-dimensionally axisymmetric mass distribution with 3 components: an exponential profile for the disk, a point mass for the bulge, and an NFW profile for the halo. The mass model is simultaneously fitted to the kinematics and the lensing data. The NFW halo needs to be oblate with a flattening of a/c=0.33^{+0.07}_{-0.05} to be consistent with the radio data. This suggests that baryons are effective at making the halos oblate near the center. The lensing and kinematics analysis probe the inner ~10 kpc of the galaxy, and we obtain a lower limit on the halo scale radius of 16 kpc (95% CI). The dark matter mass fraction inside a sphere with a radius of 2.2 disk scale lengths is f_{DM,2.2}=0.43^{+0.10}_{-0.09}. The contribution of the disk to the total circular velocity at 2.2 disk scale lengths is 0.76^{+0.05}_{-0.06}, suggesting that the disk is marginally submaximal. The stellar mass of the disk from our modeling is log_{10}(M_{*}/M_{sun}) = 11.06^{+0.09}_{-0.11} assuming that the cold gas contributes ~20% to the total disk mass. In comparison to the stellar masses estimated from stellar population synthesis models, the stellar initial mass function of Chabrier is preferred to that of Salpeter by a probability factor of 7.2.

We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged]

We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral to lenticular. For each system, we compare the total mass inside the critical curve inferred from gravitational lens modelling to the stellar mass inferred from stellar population synthesis (SPS) models, computing the stellar mass fraction f* = M(SPS)/M(lens). We find that, for the lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) leads to estimates of f* that exceed 1. This is unphysical, and provides strong evidence against the Salpeter IMF being valid for these systems. Taking the lower mass end of the SWELLS sample sigma(SIE) < 230 km/s, we find that the IMF is lighter (in terms of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with the Chabrier IMF and IMFs between the two. This result is consistent with previous studies of spiral galaxies based on independent techniques. In combination with recent studies of massive early-type galaxies that have favoured a heavier Salpeter-like IMF, this result strengthens the evidence against a universal stellar IMF.

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.

This paper presents the full VLT/VIMOS-IFU data set and related data products from an ESO Large Programme with the observational goal of obtaining two-dimensional kinematic data of early-type lens galaxies, out to one effective radius. The sample consists of 17 early-type galaxies (ETG) selected from the SLACS gravitational-lens survey. The galaxies cover the redshift range from 0.08 to 0.35 and have stellar velocity dispersions between 200 and 350 km/s. This programme is complemented by a similar observational programme on Keck, using long-slit spectroscopy. In combination with multi-band imaging data, the kinematic data provide stringent constraints on the inner mass profiles of ETGs beyond the local universe. Our Large Programme thus extends studies of nearby early-type galaxies (e.g. SAURON/ATLAS3D) by an order of magnitude in distance and toward higher masses. We provide an overview of our observational strategy, the data products (luminosity-weighted spectra and Hubble Space Telescope images) and derived products (i.e. two-dimensional fields of velocity dispersions and streaming motions) that have been used in a number of published and forthcoming lensing, kinematic and stellar-population studies.These studies also pave the way for future studies of early-type galaxies at z=1 with the upcoming extremely large telescopes.

When an ultra-high energy neutrino or cosmic ray strikes the Lunar surface a radio-frequency pulse is emitted. We plan to use the LOFAR radio telescope to detect these pulses. In this work we propose an efficient trigger implementation for LOFAR optimized for the observation of short radio pulses.

We present the first results from the X-shooter Lens Survey (XLENS): an analysis of the massive early-type galaxy SDSS J1148+1930 at redshift z=0.444. We combine its extended kinematic profile -derived from spectra obtained with X-shooter on the ESO VLT- with strong gravitational lensing and multi-color information derived from SDSS images. Our main results are (i) The luminosity-weighted stellar velocity dispersion is \sigma *(<Reff) = 351\pm10 km/s, extracted from a rectangular aperture of 1.8"x1.6", more accurate and considerably lower than a previously published value of \sim450km/s. (ii) A single-component (stellar plus dark) mass model of the lens galaxy yields a logarithmic total-density slope of {\gamma}' = 1.72+0.05-0.06. (iii) The projected stellar mass fraction, derived solely from lensing, is f*(<Reff) = 0.19+0.04-0.09 inside the Einstein radius for a Hernquist profile and no anisotropy. The dark-matter fraction inside the effective radius f_DM(<Reff) = 0.60+0.15\pm0.1. (iv) Based on the SDSS colors, we find f*(<Reff)=0.17 \pm 0.06 for a Salpeter IMF and f*(<Reff)=0.07 \pm 0.02 for a Chabrier IMF. The lensing and dynamics constraints on the stellar mass fraction agree well with those independently derived from the SDSS colors for a Salpeter IMF, preferred over Chabrier IMF at variance with standard results for lower mass galaxies. Dwarf-rich IMFs in the lower mass range of 0.1-0.7 solar mass, with {\alpha} >= 3 (with dN/dM \propto M^(-{\alpha})) -such as that recently suggested for massive ETGs- are excluded at the > 90 % C.L. and in some cases violate the total lensing-derived mass limit. We conclude that this very massive ETG is dark-matter dominated inside one effective radius, consistent with the trend recently found from massive SLACS galaxies, with a IMF normalization consistent with Salpeter.

The relative contribution of baryons and dark matter to the inner regions of spiral galaxies provides critical clues to their formation and evolution, but it is generally difficult to determine. For spiral galaxies that are strong gravitational lenses, however, the combination of lensing and kinematic observations can be used to break the disk-halo degeneracy. In turn, such data constrain fundamental parameters such as i) the mass density profile slope and axis ratio of the dark matter halo, and by comparison with dark matter-only numerical simulations the modifications imposed by baryons; ii) the mass in stars and therefore the overall star formation efficiency, and the amount of feedback; iii) by comparison with stellar population synthesis models, the normalization of the stellar initial mass function. In this first paper of a series, we present a sample of 16 secure, 1 probable, and 6 possible strong lensing spiral galaxies, for which multi-band high-resolution images and rotation curves were obtained using the Hubble Space Telescope and Keck-II Telescope as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS). The sample includes 8 newly discovered secure systems. [abridged] We find that the SWELLS sample of secure lenses spans a broad range of morphologies (from lenticular to late-type spiral), spectral types (quantified by Halpha emission), and bulge to total stellar mass ratio (0.22-0.85), while being limited to M_*>10^{10.5} M_sun. The SWELLS sample is thus well-suited for exploring the relationship between dark and luminous matter in a broad range of galaxies. We find that the deflector galaxies obey the same size-mass relation as that of a comparison sample of elongated non-lens galaxies selected from the SDSS survey. We conclude that the SWELLS sample is consistent with being representative of the overall population of high-mass high-inclination disky galaxies.

We exploit the strong lensing effect to explore the properties of intrinsically faint and compact galaxies at intermediate redshift, at the highest possible resolution at optical wavelengths. Our sample consists of 46 strongly-lensed emission line galaxies discovered by the Sloan Lens ACS (SLACS) Survey. The galaxies have been imaged at high resolution with HST in three bands (V_HST, I_814 and H_160), allowing us to infer their size, luminosity, and stellar mass using stellar population synthesis models. Lens modeling is performed using a new fast and robust code, klens, which we test extensively on real and synthetic non-lensed galaxies, and also on simulated galaxies multiply-imaged by SLACS- like galaxy-scale lenses. Our tests show that our measurements of galaxy size, flux, and Sersic index are robust and accurate, even for objects intrinsically smaller than the HST point spread function. The median magnification is 8.8, with a long tail that extends to magnifications above 40. Modeling the SLACS sources reveals a population of galaxies with colors and Sersic indices (median n ~ 1) consistent with the objects detected in the field with HST in the GEMS survey, but that are (typically) ~ 2 magnitudes fainter and ~ 5 times smaller in apparent size. The closest analog are ultracompact emission line galaxies identified by HST grism surveys. The lowest mass galaxies in our sample are comparable to the brightest Milky Way satellites in stellar mass (10^7 solar masses) and have well-determined half light radii of 0."05 (~0.3 kpc).

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

We combine in a self-consistent way the constraints from both gravitational lensing and stellar kinematics to perform a detailed investigation of the internal mass distribution, amount of dark matter, and dynamical structure of the 16 early-type lens galaxies from the SLACS Survey, at z = 0.08 - 0.33, for which both HST/ACS and NICMOS high-resolution imaging and VLT VIMOS IFU spectroscopy are available. Based on this data set, we analyze the inner regions of the galaxies, i.e. typically within one (3D) effective radius r_e, under the assumption of axial symmetry and by constructing dynamical models supported by two-integral stellar DFs. For all systems, the total mass density distribution is found to be well approximated by a simple power-law: this profile is on average slightly super-isothermal, with a logarithmic slope <gamma'> = 2.074^{+0.043}_{-0.041} (68% CL) and an intrinsic scatter 0.144^{+0.055}_{-0.014}, and is fairly round, with an average axial ratio <q> = 0.77+/-0.04. The lower limit for the dark matter fraction (fDM) inside r_e ranges, in individual systems, from nearly zero to almost a half, with a median value of 12%. By including stellar masses derived from SPS models with a Salpeter IMF, we obtain an average fDM = 31%. The fDM rises to 61% if, instead, a Chabrier IMF is assumed. For both IMFs, the dark matter fraction increases with the total mass of the galaxy (3-sigma correlation). Based on the intrinsic angular momentum parameter calculated from our models, we find that the galaxies can be divided into two dynamically distinct groups, which are shown to correspond to the usual classes of the slow and fast rotators. Overall, the SLACS systems are structurally and dynamically very similar to their nearby counterparts, indicating that the inner regions of early-type galaxies have undergone little, if any, evolution since redshift z ~ 0.35. (Abridged)

The degeneracy among the disk, bulge and halo contributions to galaxy rotation curves prevents an understanding of the distribution of baryons and dark matter in disk galaxies. In an attempt to break this degeneracy, we present an analysis of the spiral galaxy strong gravitational lens SDSS J2141-0001, discovered as part of the SLACS survey. We present new Hubble Space Telescope multicolor imaging, gas and stellar kinematics data derived from long-slit spectroscopy, and K-band LGS adaptive optics imaging, both from the Keck telescopes. We model the galaxy as a sum of concentric axisymmetric bulge, disk and halo components and infer the contribution of each component, using information from gravitational lensing and gas kinematics. This analysis yields a best-fitting total (disk plus bulge) stellar mass of log_{10}(Mstar/Msun) = 10.99(+0.11,-0.25). The photometric data combined with stellar population synthesis models yield log_{10}(Mstar/Msun) = 10.97\pm0.07, and 11.21\pm0.07 for the Chabrier and Salpeter IMFs, respectively. Accounting for the expected gas fraction of \simeq 20% reduces the lensing plus kinematics stellar mass by 0.10\pm0.05 dex, resulting in a Bayes factor of 11.9 in favor of a Chabrier IMF. The dark matter halo is roughly spherical, with minor to major axis ratio q_{halo}=0.91(+0.15,-0.13). The dark matter halo has a maximum circular velocity of V_{max}=276(+17,-18) km/s, and a central density parameter of log_{10}\Delta_{V/2}=5.9(+0.9,-0.5). This is higher than predicted for uncontracted dark matter haloes in LCDM cosmologies, log_{10}\Delta_{V/2}=5.2, suggesting that either the halo has contracted in response to galaxy formation, or that the halo has a higher than average concentration. At 2.2 disk scale lengths the dark matter fraction is f_{DM}=0.55(+0.20,-0.15), suggesting that SDSS J2141-0001 is sub-maximal.

We compare galaxy number counts in HST/ACS fields containing moderate-redshift (0.2<z<1.0) strong gravitational lenses with those in two control samples: (1) the first square degree of the COSMOS survey, comprising 259 ACS fields and (2) 20 "pure parallel" fields randomly located on the sky. Through a Bayesian analysis we determine the expectation values (mu_0) and confidence levels of the underlying number counts for a range of apertures and magnitude bins. Our analysis has produced the following results: (i) We infer that our control samples are not consistent, with the number counts in the COSMOS sample being significantly higher than in the pure parallel sample for 21 <= F814W <= 23. This result matches those found in previous analyses of COSMOS data using different techniques. (ii) We find that small-size apertures, centered on strong lenses, are overdense compared with randomly placed apertures in the control samples, even compared to the COSMOS sample. Correcting for the local clustering of elliptical galaxies, based on the average two-point correlation function reduces this overdensity to the 1-2 sigma level. Thus, the overdensity of galaxies seen along a typical line of sight to a lens can be explained mostly by the natural clustering of galaxies, rather than being due to lenses lying along otherwise biased lines of sight. However, a larger sample of lenses is needed to see if the remaining bias persists when the lens-field uncertainties are smaller. (iii) There is considerable scatter in the lines of sight to _individual_ lens systems, but quantities that are linearly dependent on the external convergence (e.g., H_0) should become unbiased if the extra galaxies that cause the bias can be accounted for in the lens models either through direct modeling or via an informed prior on the external convergence. The number counts can used to set such an informed prior.

The aim of the new generation of radio synthesis arrays such as LOFAR and SKA is to achieve much higher sensitivity, resolution and frequency coverage than what is available now, especially at low frequencies. To accomplish this goal, the accuracy of the calibration techniques used is of considerable importance. Moreover, since these telescopes produce huge amounts of data, speed of convergence of calibration is a major bottleneck. The errors in calibration are due to system noise (sky and instrumental) as well as the estimation errors introduced by the calibration technique itself, which we call solver noise. We define solver noise as the distance between the optimal solution (the true value of the unknowns, uncorrupted by the system noise) and the solution obtained by calibration. We present the Space Alternating Generalized Expectation Maximization (SAGE) calibration technique, which is a modification of the Expectation Maximization algorithm, and compare its performance with the traditional Least Squares calibration based on the level of solver noise introduced by each technique. For this purpose, we develop statistical methods that use the calibrated solutions to estimate the level of solver noise. The SAGE calibration algorithm yields very promising results both in terms of accuracy and speed of convergence. The comparison approaches we adopt introduce a new framework for assessing the performance of different calibration schemes.

Tseliakhovich & Hirata recently discovered that higher-order corrections to the cosmological linear-perturbation theory lead to supersonic coherent baryonic flows just after recombination (i.e. z~1020), with rms velocities of ~30 km/s relative to the underlying dark-matter distribution, on comoving scales of ~<3 Mpc/h. To study the impact of these coherent flows we performed high-resolution N-body plus SPH simulations in boxes of 5.0 and 0.7 Mpc/h, for bulk-flow velocities of 0 (as reference), 30 and 60 km/s. The simulations follow the evolution of cosmic structures by taking into account detailed, primordial, non-equilibrium gas chemistry (i.e. H, He, H2, HD, HeH, etc.), cooling, star formation, and feedback effects from stellar evolution. We find that these bulk flows suppress star formation in low-mass haloes (i.e. Mvir<~10^8Msun until z~13, lower the abundance of the first objects by ~1%-20%, and, as consequence, delay cosmic star formation history by ~2e7yr. The gas fractions in individual objects can change up to a factor of two at very early times. Coherent bulk flow, therefore, has implications for (i) the star-formation in the lowest-mass haloes (e.g. dSphs), (ii) the start of reionization by suppressing it in some patches of the Universe, and (iii) the heating (i.e. spin temperature) of neutral hydrogen. We speculate that the patchy nature of reionization and heating on several Mpc scales could lead to enhanced differences in the HI spin-temperature, giving rise to stronger variations in the HI brightness temperatures during the late dark ages.

The Hubble constant value is currently known to 10% accuracy unless assumptions are made for the cosmology (Sandage et al. 2006). Gravitational lens systems provide another probe of the Hubble constant using time delay measurements. However, current investigations of ~20 time delay lenses, albeit of varying levels of sophistication, have resulted in different values of the Hubble constant ranging from 50-80 km/s/Mpc. In order to reduce uncertainties, more time delay measurements are essential together with better determined mass models (Oguri 2007, Saha et al. 2006). We propose a more efficient technique for measuring time delays which does not require regular monitoring with a high-resolution interferometer array. The method uses double image and long-axis quadruple lens systems in which the brighter component varies first and dominates the total flux density. Monitoring the total flux density with low-resolution but high sensitivity radio telescopes provides the variation of the brighter image and is used to trigger high-resolution observations which can then be used to see the variation in the fainter image. We present simulations of this method together with a pilot project using the WSRT (Westerbork Radio Synthesis Telescope) to trigger VLA (Very Large Array) 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 many SKA pathfinder telescopes, such as MeerKAT (Karoo Array Telescope) and ASKAP (Australian Square Kilometre Array Pathfinder), have high sensitivity but limited resolution.

In the coming years a new insight into galaxy formation and the thermal history of the Universe is expected to come from the detection of the highly redshifted cosmological 21 cm line. The cosmological 21 cm line signal is buried under Galactic and extragalactic foregrounds which are likely to be a few orders of magnitude brighter. Strategies and techniques for effective subtraction of these foreground sources require a detailed knowledge of their structure in both intensity and polarization on the relevant angular scales of 1-30 arcmin. We present results from observations conducted with the Westerbork telescope in the 140-160 MHz range with 2 arcmin resolution in two fields located at intermediate Galactic latitude, centred around the bright quasar 3C196 and the North Celestial Pole. They were observed with the purpose of characterizing the foreground properties in sky areas where actual observations of the cosmological 21 cm line could be carried out. The polarization data were analysed through the rotation measure synthesis technique. We have computed total intensity and polarization angular power spectra. Total intensity maps were carefully calibrated, reaching a high dynamic range, 150000:1 in the case of the 3C196 field. [abridged]

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We use stellar masses, photometry, lensing, and velocity dispersions to investigate empirical correlations for the final sample of 73 early-type lens galaxies (ETGs) from the SLACS survey. The traditional correlations (Fundamental Plane [FP] and its projections) are consistent with those found for non-lens galaxies, supporting the thesis that SLACS lens galaxies are representative of massive ETGs. The addition of strong lensing estimates of the total mass allows us to gain further insights into their internal structure: i) the mean slope of the total mass density profile is <gamma'> = 2.078+/-0.027 with an intrinsic scatter of 0.16+/-0.02; ii) gamma' correlates with effective radius and central mass density, in the sense that denser galaxies have steeper profiles; iii) the dark matter fraction within reff/2 is a monotonically increasing function of galaxy mass and size; iv) the dimensional mass M_dim is proportional to the total mass, and both increase more rapidly than stellar mass M*; v) the Mass Plane (MP), obtained by replacing surface brightness with surface mass density in the FP, is found to be tighter and closer to the virial relation than the FP and the M*P, indicating that the scatter of those relations is dominated by stellar population effects; vi) we construct the Fundamental Hyper-Plane by adding stellar masses to the MP and find the M* coefficient to be consistent with zero and no residual intrinsic scatter. Our results demonstrate that the dynamical structure of ETGs is not scale invariant and that it is fully specified by the total mass, r_eff, and sigma. Although the basic trends can be explained qualitatively in terms of varying star formation efficiency as a function of halo mass and as the result of dry and wet mergers, reproducing quantitatively the observed correlations and their tightness may be a significant challenge for galaxy formation models.

We use stellar dynamics, strong lensing, stellar population synthesis models, and weak lensing shear measurements to constrain the dark matter (DM) profile and stellar mass in a sample of 53 massive early-type galaxies. We explore three DM halo models (unperturbed Navarro Frenk & White [NFW] halos and the adiabatic contraction models of Blumenthal and Gnedin) and impose a model for the relationship between the stellar and virial mass (i.e., a relationship for the star-formation efficiency as a function of halo mass). We show that, given our model assumptions, the data clearly prefer a Salpeter-like initial mass function (IMF) over a lighter IMF (e.g., Chabrier or Kroupa), irrespective of the choice of DM halo. In addition, we find that the data prefer at most a moderate amount of adiabatic contraction (Blumenthal adiabatic contraction is strongly disfavored) and are only consistent with no adiabatic contraction (i.e., a NFW halo) if a mass-dependent IMF is assumed, in the sense of a more massive normalization of the IMF for more massive halos.

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.

We study the internal dynamical structure of the early-type lens galaxy SDSS J0728+3835 at z = 0.206. The analysis is based on two-dimensional kinematic maps extending out to 1.7 effective radii obtained from Keck spectroscopy, on lensing geometry and on stellar mass estimates obtained from multiband Hubble Space Telescope imaging. The data are modelled under the assumptions of axial symmetry supported by a two-integral distribution function (DF), by applying the combined gravitational lensing and stellar dynamics code CAULDRON, and yielding high-quality constraints for an early-type galaxy at cosmological redshifts. Modelling the total density profile as a power-law of the form rho_tot ~ 1/r^{gamma}, we find that it is nearly isothermal (logarithmic slope gamma = 2.08^{+0.04}_{-0.02}), and quite flattened (axial ratio q = 0.60^{+0.08}_{-0.03}). The galaxy is mildly anisotropic (delta = 0.08 +/- 0.02) and shows a fair amount of rotational support, in particular towards the outer regions. We determine a dark matter fraction lower limit of 28 per cent within the effective radius. The stellar contribution to the total mass distribution is close to maximal for a Chabrier initial mass function (IMF), whereas for a Salpeter IMF the stellar mass exceeds the total mass within the galaxy inner regions. We find that the combination of a NFW dark matter halo with the maximally rescaled luminous profile provides a remarkably good fit to the total mass distribution over a broad radial range. Our results confirm and expand the findings of the SLACS survey for early-type galaxies of comparable velocity dispersion (sigma_SDSS = 214 +/- 11 km/s). The internal structure of J0728 is consistent with that of local early-type galaxies of comparable velocity dispersion as measured by the SAURON project, suggesting lack of evolution in the past two billion years.

One of the aims of the Low Frequency Array (LOFAR) Epoch of Reionization (EoR) project is to measure the power spectrum of variations in the intensity of redshifted 21-cm radiation from the EoR. The sensitivity with which this power spectrum can be estimated depends on the level of thermal noise and sample variance, and also on the systematic errors arising from the extraction process, in particular from the subtraction of foreground contamination. We model the extraction process using realistic simulations of the cosmological signal, the foregrounds and noise, and so estimate the sensitivity of the LOFAR EoR experiment to the redshifted 21-cm power spectrum. Detection of emission from the EoR should be possible within 360 hours of observation with a single station beam. Integrating for longer, and synthesizing multiple station beams within the primary (tile) beam, then enables us to extract progressively more accurate estimates of the power at a greater range of scales and redshifts. We discuss different observational strategies which compromise between depth of observation, sky coverage and frequency coverage. A plan in which lower frequencies receive a larger fraction of the time appears to be promising. We also study the nature of the bias which foreground fitting errors induce on the inferred power spectrum, and discuss how to reduce and correct for this bias. The angular and line-of-sight power spectra have different merits in this respect, and we suggest considering them separately in the analysis of LOFAR data.

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.

Strong gravitational lens systems with measured time delays between the multiple images provide a method for measuring the "time-delay distance" to the lens, and thus the Hubble constant. We present a Bayesian analysis of the strong gravitational lens system B1608+656, incorporating (i) new, deep Hubble Space Telescope (HST) observations, (ii) a new velocity dispersion measurement of 260+/-15 km/s for the primary lens galaxy, and (iii) an updated study of the lens' environment. When modeling the stellar dynamics of the primary lens galaxy, the lensing effect, and the environment of the lens, we explicitly include the total mass distribution profile logarithmic slope gamma' and the external convergence kappa_ext; we marginalize over these parameters, assigning well-motivated priors for them, and so turn the major systematic errors into statistical ones. The HST images provide one such prior, constraining the lens mass density profile logarithmic slope to be gamma'=2.08+/-0.03; a combination of numerical simulations and photometric observations of the B1608+656 field provides an estimate of the prior for kappa_ext: 0.10 +0.08/-0.05. This latter distribution dominates the final uncertainty on H_0. Compared with previous work on this system, the new data provide an increase in precision of more than a factor of two. In combination with the WMAP 5-year data set, we find that the B1608+656 data set constrains the curvature parameter to be -0.031 < Omega_k < 0.009 (95% CL), a level of precision comparable to that afforded by the current Type Ia SNe sample. Asserting a flat spatial geometry, we find that, in combination with WMAP, H_0 = 69.7 +4.9/-5.0 km/s/Mpc and w=-0.94 +0.17/-0.19 (68% CL), suggesting that the observations of B1608+656 constrain w as tightly as do the current Baryon Acoustic Oscillation data. (abridged)

We construct a mass model for the spiral lens galaxy 2237+0305, at redshift z_l=0.04, based on gravitational-lensing constraints, HI rotation, and new stellar-kinematic information, based on data taken with the ESI spectrograph on the 10m Keck-II Telescope. High resolution rotation curves and velocity dispersion profiles along two perpendicular directions, close to the major and minor axes of the lens galaxy, were obtained by fitting the Mgb-Fe absorption line region. The stellar rotation curve rises slowly and flattens at r~1.5" (~1.1 kpc). The velocity dispersion profile is approximately flat. A combination of photometric, kinematic and lensing information is used to construct a mass model for the four major mass components of the system -- the dark matter halo, disc, bulge, and bar. The best-fitting solution has a dark matter halo with a logarithmic inner density slope of gamma=0.9+/-0.3 for rho_DM propto r^-gamma, a bulge with M/L_B=6.6+/-0.3 Upsilon_odot, and a disc with M/L_B =1.2+/-0.3 Upsilon_odot, in agreement with measurements of late-type spirals. The bulge dominates support in the inner regions where the multiple images are located and is therefore tightly constrained by the observations. The disc is sub-maximal and contributes 45+/-11 per cent of the rotational support of the galaxy at 2.2r_d. The halo mass is (2.0+/-0.6) x 10^12 M_odot, and the stellar to virial mass ratio is 7.0+/-2.3 per cent, consistent with typical galaxies of the same mass.

We determine an absolute calibration of the initial mass function (IMF) of early-type galaxies, by studying a sample of 56 gravitational lenses identified by the SLACS Survey. Under the assumption of standard Navarro, Frenk & White dark matter halos, a combination of lensing, dynamical, and stellar population synthesis models is used to disentangle the stellar and dark matter contribution for each lens. We define an "IMF mismatch" parameter \alpha=M*(L+D)/M*(SPS) as the ratio of stellar mass inferred by a joint lensing and dynamical models (M*(L+D)) to the current stellar mass inferred from stellar populations synthesis models (M*(SPS)). We find that a Salpeter IMF provides stellar masses in agreement with those inferred by lensing and dynamical models (<\log \alpha>=0.00+-0.03+-0.02), while a Chabrier IMF underestimates them (<\log \alpha>=0.25+-0.03+-0.02). A tentative trend is found, in the sense that \alpha appears to increase with galaxy velocity dispersion. Taken at face value, this result would imply a non universal IMF, perhaps dependent on metallicity, age, or abundance ratios of the stellar populations. Alternatively, the observed trend may imply non-universal dark matter halos with inner density slope increasing with velocity dispersion. While the degeneracy between the two interpretations cannot be broken without additional information, the data imply that massive early-type galaxies cannot have both a universal IMF and universal dark matter halos.

We present the current photometric dataset for the Sloan Lens ACS (SLACS) Survey, including HST photometry from ACS, WFPC2, and NICMOS. These data have enabled the confirmation of an additional 15 grade `A' (certain) lens systems, bringing the number of SLACS grade `A' lenses to 85; including 13 grade `B' (likely) systems, SLACS has identified nearly 100 lenses and lens candidates. Approximately 80% of the grade `A' systems have elliptical morphologies while ~10% show spiral structure; the remaining lenses have lenticular morphologies. Spectroscopic redshifts for the lens and source are available for every system, making SLACS the largest homogeneous dataset of galaxy-scale lenses to date. We have developed a novel Bayesian stellar population analysis code to determine robust stellar masses with accurate error estimates. We apply this code to deep, high-resolution HST imaging and determine stellar masses with typical statistical errors of 0.1 dex; we find that these stellar masses are unbiased compared to estimates obtained using SDSS photometry, provided that informative priors are used. The stellar masses range from 10^10.5 to 10^11.8 M$_\odot$ and the typical stellar mass fraction within the Einstein radius is 0.4, assuming a Chabrier IMF. The ensemble properties of the SLACS lens galaxies, e.g. stellar masses and projected ellipticities, appear to be indistinguishable from other SDSS galaxies with similar stellar velocity dispersions. This further supports that SLACS lenses are representative of the overall population of massive early-type galaxies with M* >~ 10^11 M$_\odot$, and are therefore an ideal dataset to investigate the kpc-scale distribution of luminous and dark matter in galaxies out to z ~ 0.5.

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

SMM J04542-0301 is an extended (~1 arcmin) sub-mm source located near the core of the cluster MS0451.6-0305. It has been suggested that part of its emission arises from the interaction between a LBG and two EROs at z~2.9 that are multiply-imaged. However, the dramatic resolution difference between the sub-mm map and the optical/NIR images make it difficult to confirm this hypothesis. In this paper, we present a deep (~ 10 microJy/beam), high resolution (~2 arcsec) 1.4 GHz radio map of the cluster core, in which we have identified 6 sources located within SMM J04542-0301. The strong lensing effect in the radio data has been quantified by constructing a new lens model of the cluster. The brightest and most extended of these sources (RJ) is located in the middle of the sub-mm emission, and has no obvious counterpart in the optical/NIR. Three other detections (E1, E2 and E3) seem to be associated with the images of one of the EROs. The last two detections (CR1 and CR2), for which no optical/NIR counterpart have been found, seem to constitute two relatively compact emitting regions embedded in a ~5 arcsec extended radio source located at the position of the sub-mm peak. The presence of this extended component can only be explained if it is being produced by a lensed region of dust obscured star formation in the center of the merger. A comparison between the radio and sub-mm data at the same resolution suggests that E1, E2, E3, CR1 and CR2 are associated with the sub-mm emission. The radio observations provide strong observational evidence in favor of the merger hypothesis. However, the question if RJ is also contributing to the observed sub-mm emission remains open. These results illustrate the promising prospects for radio interferometry and strong gravitational lensing to study the internal structure of SMGs.

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

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.

The Universe's Epoch of Reionization (EoR) can be studied using a number of observational probes that provide complementary or corroborating information. Each of these probes suffers from its own systematic and statistical uncertainties. It is therefore useful to consider the mutual information that these data sets contain. In this paper we present a cross-correlation study between the kinetic Sunyaev-Zel'dovich effect (kSZ) -- produced by the scattering of CMB photons off free electrons produced during the reionization process -- and the cosmological 21cm signal -- which reflects the neutral hydrogen content of the Universe, as a function of redshift. The study is carried out using a simulated reionization history in 100Mpc/h scale N-body simulations with radiative transfer. In essence we find that the two probes anti-correlate. The significance of the anti-correlation signal depends on the extent of the reionization process, wherein extended histories result in a much stronger signal compared to instantaneous cases. Unfortunately however, once the primary CMB fluctuations are included into our simulation they serve as a source of large correlated noise that renders the cross-correlation signal insignificant, regardless of the reionization scenario.

We present the first results from a series of observations conducted with the Westerbork telescope in the 140--160 MHz range with a 2 arcmin resolution aimed at characterizing the properties of the foregrounds for epoch of reionization experiments. For the first time we have detected fluctuations in the Galactic diffuse emission on scales greater than 13 arcmin at 150 MHz, in the low Galactic latitude area known as Fan region. Those fluctuations have an $rms$ of 14 K. The total intensity power spectrum shows a power--law behaviour down to $\ell \sim 900$ with slope $\beta^I_\ell = -2.2 \pm 0.3$. The detection of diffuse emission at smaller angular scales is limited by residual point sources. We measured an $rms$ confusion noise of $\sim$3 mJy beam$^{-1}$. Diffuse polarized emission was also detected for the first time at this frequency. The polarized signal shows complex structure both spatially and along the line of sight. The polarization power spectrum shows a power--law behaviour down to $\ell \sim 2700$ with slope $\beta^P_\ell = -1.65 \pm 0.15$. The $rms$ of polarization fluctuations is 7.2 K on 4 arcmin scales. By extrapolating the measured spectrum of total intensity emission, we find a contamination on the cosmological signal of $\delta T= \sqrt{\ell (\ell+1) C^I_\ell / 2\pi} \sim 5.7$ K on 5 arcmin scales and a corresponding $rms$ value of $\sim$18.3 K at the same angular scale. The level of the polarization power spectrum is $\delta T \sim 3.3$ K on 5 arcmin scales. Given its exceptionally bright polarized signal, the Fan region is likely to represent an upper limit on the sky brightness at moderate and high Galactic latitude.

An obstacle to the detection of redshifted 21cm emission from the epoch of reionization (EoR) is the presence of foregrounds which exceed the cosmological signal in intensity by orders of magnitude. We argue that in principle it would be better to fit the foregrounds non-parametrically - allowing the data to determine their shape - rather than selecting some functional form in advance and then fitting its parameters. Non-parametric fits often suffer from other problems, however. We discuss these before suggesting a non-parametric method, Wp smoothing, which seems to avoid some of them. After outlining the principles of Wp smoothing we describe an algorithm used to implement it. We then apply Wp smoothing to a synthetic data cube for the LOFAR EoR experiment. The performance of Wp smoothing, measured by the extent to which it is able to recover the variance of the cosmological signal and to which it avoids leakage of power from the foregrounds, is compared to that of a parametric fit, and to another non-parametric method (smoothing splines). We find that Wp smoothing is superior to smoothing splines for our application, and is competitive with parametric methods even though in the latter case we may choose the functional form of the fit with advance knowledge of the simulated foregrounds. Finally, we discuss how the quality of the fit is affected by the frequency resolution and range, by the characteristics of the cosmological signal and by edge effects.

We present the first detailed analysis of the mass and dynamical structure of a sample of six early-type lens galaxies, selected from the Sloan Lens ACS Survey, in the redshift range 0.08 < z < 0.33. Both Hubble Space Telescope (HST)/ACS high-resolution imaging and VLT VIMOS integral field spectroscopy are available for these systems. The galaxies are modelled - under the assumptions of axial symmetry and two-integral stellar distribution function - by making use of the CAULDRON code, which self-consistently combines gravitational lensing and stellar dynamics, and is fully embedded within the framework of Bayesian statistics. The principal results of this study are: (i) all galaxies in the sample are well described by a simple axisymmetric power-law profile for the total density, with a logarithmic slope gamma very close to isothermal (<gamma> = 1.98 +/- 0.05 and an intrinsic spread close to 5 per cent) showing no evidence of evolution over the probed range of redshift; (ii) the axial ratio of the total density distribution is rounder than 0.65 and in all cases, except for a fast rotator, does not deviate significantly from the flattening of the intrinsic stellar distribution; (iii) the dark matter fraction within the effective radius has a lower limit of about 15 to 30 per cent; (iv) the sample galaxies are only mildly anisotropic, with delta <= 0.16; (v) the physical distinction among slow and fast rotators, quantified by the v/sigma ratio and the intrinsic angular momentum, is already present at z > 0.1. Altogether, early-type galaxies at z = 0.08 - 0.33 are found to be markedly smooth and almost isothermal systems, structurally and dynamically very similar to their nearby counterparts. (Abridged)

We apply the joint lensing and dynamics code for the analysis of early-type galaxies, CAULDRON, to a rotating N-body stellar system with dark matter halo which significantly violates the two major assumptions of the method, i.e. axial symmetry supported by a two-integral distribution function. The goal is to study how CAULDRON performs in an extreme case, and to determine which galaxy properties can still be robustly recovered. Three data sets, corresponding to orthogonal lines of sight, are generated from the N-body system and analysed with the identical procedure followed in the study of real lens galaxies, adopting an axisymmetric power-law total density distribution. We find that several global properties of the N-body system are recovered with remarkable accuracy, despite the fact that the adopted power-law model is too simple to account for the lack of symmetry of the true density distribution. In particular, the logarithmic slope of the total density distribution is robustly recovered to within less than 10 per cent (with the exception of the ill-constrained very inner regions), the inferred angle-averaged radial profile of the total mass closely follows the true distribution, and the dark matter fraction of the system (inside the effective radius) is correctly determined within ~ 10 per cent of the total mass. Unless the line of sight direction is almost parallel to the total angular momentum vector of the system, reliably recovered quantities also include the angular momentum, the V/sigma ratio, and the anisotropy parameter delta. We conclude that the CAULDRON code can be safely and effectively applied to real early-type lens galaxies, providing reliable information also for systems that depart significantly from the method's assumptions.

MG 2016+112 is a quadruply imaged lens system with two complete images A and B and a pair of merging partial images in region C as seen in the radio. The merging images are found to violate the expected mirror symmetry. This indicates an astrometric anomaly which could only be of gravitational origin and could arise due to substructure in the environment or line-of-sight of the lens galaxy. We present new high resolution multi-frequency VLBI observations at 1.7, 5 and 8.4 GHz. Three new components are detected in the new VLBI imaging of both the lensed images A and B. The expected opposite parity of the lensed images A and B was confirmed due to the detection of non-collinear components. Furthermore, the observed properties of the newly detected components are inconsistent with the predictions of previous mass models. We present new scenarios for the background quasar which are consistent with the new observations. We also investigate the role of the satellite galaxy situated at the same redshift as the main lensing galaxy. Our new mass models demonstrate quantitatively that the satellite galaxy is the primary cause of the astrometric anomaly found in region C. The detected satellite is consistent with the abundance of subhaloes expected in the halo from cold dark matter (CDM) simulations. However, the fraction of the total halo mass in the satellite as computed from lens modeling is found to be higher than that predicted by CDM simulations.

We outline the case for a comprehensive wide and deep survey ultimately targeted at obtaining 21-cm HI line emission spectroscopic observations of more than a billion galaxies to redshift z=1.5 and greater over half the sky. This survey provides a database of galaxy redshifts, HI gas masses, and galaxy rotation curves that would enable a wide range of science, including fundamental cosmology and studies of Dark Energy. This science requires the next generation of radio arrays, which are being designed under the umbrella of the Square Kilometer Array (SKA) project. We present a science roadmap, extending to 2020 and beyond, that would enable this ambitious survey. We also place this survey in the context of other multi-wavelength surveys.

Detecting redshifted 21cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, and we suggest that using the correct statistics may tease out signatures of reionization. We generate mock datacubes simulating the output of the Low Frequency Array (LOFAR) EoR experiment. These cubes combine realistic models for Galactic and extragalactic foregrounds and the noise with three different simulations of the CS. We fit out the foregrounds, which are smooth in the frequency direction, to produce residual images in each frequency band. We denoise these images and study the skewness of the one-point distribution in the images as a function of frequency. We find that, under sufficiently optimistic assumptions, we can recover the main features of the redshift evolution of the skewness in the 21cm signal. We argue that some of these features - such as a dip at the onset of reionization, followed by a rise towards its later stages - may be generic, and give us a promising route to a statistical detection of reionization.

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.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

A number of experiments are set to measure the 21-cm signal of neutral hydrogen from the Epoch of Reionization (EoR). The common denominator of these experiments are the large data sets produced, contaminated by various instrumental effects, ionospheric distortions, RFI and strong Galactic and extragalactic foregrounds. In this paper, the first in a series, we present the Data Model that will be the basis of the signal analysis for the LOFAR (Low Frequency Array) EoR Key Science Project (LOFAR EoR KSP). Using this data model we simulate realistic visibility data sets over a wide frequency band, taking properly into account all currently known instrumental corruptions (e.g. direction-dependent gains, complex gains, polarization effects, noise, etc). We then apply primary calibration errors to the data in a statistical sense, assuming that the calibration errors are random Gaussian variates at a level consistent with our current knowledge based on observations with the LOFAR Core Station 1. Our aim is to demonstrate how the systematics of an interferometric measurement affect the quality of the calibrated data, how errors correlate and propagate, and in the long run how this can lead to new calibration strategies. We present results of these simulations and the inversion process and extraction procedure. We also discuss some general properties of the coherency matrix and Jones formalism that might prove useful in solving the calibration problem of aperture synthesis arrays. We conclude that even in the presence of realistic noise and instrumental errors, the statistical signature of the EoR signal can be detected by LOFAR with relatively small errors. A detailed study of the statistical properties of our data model and more complex instrumental models will be considered in the future.

Strong gravitational lensing is a powerful technique for probing galaxy mass distributions and for measuring cosmological parameters. We present a pixelated approach to modeling simultaneously the lens potential and source intensity of strong gravitational lens systems with extended source-intensity distributions. For systems with sources of sufficient extent such that the separate lensed images are connected by intensity measurements, the accuracy in the reconstructed potential is solely limited by the quality of the data. We apply this potential reconstruction technique to deep HST observations of B1608+656, a four-image gravitational lens system formed by a pair of interacting lens galaxies. We present a comprehensive Bayesian analysis of the system that takes into account the extended source-intensity distribution, dust extinction, and the interacting lens galaxies. Our approach allows us to compare various models of the components of the lens system, which include the point-spread function (PSF), dust, lens galaxy light, source-intensity distribution, and lens potential. Using optimal combinations of the PSF, dust, and lens galaxy light models, we successfully reconstruct both the lens potential and the extended source-intensity distribution of B1608+656. The resulting reconstruction can be used as the basis of a measurement of the Hubble constant. We use our reconstruction of the gravitational potential to study the relative distribution of mass and light in the lensing galaxies. We find that the mass-to-light ratio for the primary lens galaxy is (2.0+/-0.2)h M_{\sun} L_{B,\sun}^{-1} within the Einstein radius 3.9 h^{-1} kpc, in agreement with what is found for noninteracting lens galaxies at the same scales. (Abridged)

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

The combination of two-dimensional kinematics and gravitational lens modelling permits detailed reconstruction of the phase-space structure of early-type galaxies and sets constraints on the dark-matter distribution in their inner regions. We describe a project which combines integral-field spectroscopy from an ESO Large Programme using VIMOS on the VLT with deep HST ACS and NICMOS images to study a sample of 17 early-type lens galaxies at redshifts between 0.1 and 0.3, drawn from the Sloan Lens ACS survey (SLACS).

Radio Interferometry is an essential method for astronomical observations. Self-calibration techniques have increased the quality of the radio astronomical observations (and hence the science) by orders of magnitude. Recently, there is a drive towards sensor arrays built using inexpensive hardware and distributed over a wide area acting as radio interferometers. Calibration of such arrays poses new problems in terms of computational cost as well as in performance of existing calibration algorithms. We consider the application of the Space Alternating Generalized Expectation Maximization (SAGE) \cite{Fess94} algorithm for calibration of radio interferometric arrays. Application to real data shows that this is an improvement over existing calibration algorithms that are based on direct, deterministic non linear optimization. As presented in this paper, we can improve the computational cost as well as the quality of the calibration using this algorithm.

We study the relation between the internal structure of early-type galaxies and their environment using 70 strong gravitational lenses from the Sloan ACS Lens Survey. The Sloan database is used to determine two measures of overdensity of galaxies around each lens: the projected number density of galaxies inside the tenth nearest neighbor (\Sigma_{10}) and within a cone of radius one h^{-1} Mpc (D_1). Our main results are: 1) The average overdensity is somewhat larger than unity, consistent with lenses preferring overdense environments as expected for massive early-type galaxies (12/70 lenses are in known groups/clusters). 2) The distribution of overdensities is indistinguishable from that of "twin" non-lens galaxies selected from SDSS to have the same redshift and stellar velocity dispersion \sigma_*. Thus, within our errors, lens galaxies are an unbiased population, and the SLACS results can be generalized to the overall population of early-type galaxies. 3) Typical contributions from external mass distribution are no more than a few per cent, reaching 10-20% (~0.05-0.10 external convergence) only in the most extreme overdensities. 4) No significant correlation between overdensity and slope of the mass density profile of the lens is found. 5) Satellite galaxies (those with a more luminous companion) have marginally steeper mass density profiles than central galaxies (those without). This result suggests that tidal stripping may affect the mass structure of early-type galaxies down to kpc scales probed by strong lensing, when they fall into larger structures [ABRIDGED].

We present an efficient method to generate large simulations of the Epoch of Reionization (EoR) without the need for a full 3-dimensional radiative transfer code. Large dark-matter-only simulations are post-processed to produce maps of the redshifted 21cm emission from neutral hydrogen. Dark matter haloes are embedded with sources of radiation whose properties are either based on semi-analytical prescriptions or derived from hydrodynamical simulations. These sources could either be stars or power-law sources with varying spectral indices. Assuming spherical symmetry, ionized bubbles are created around these sources, whose radial ionized fraction and temperature profiles are derived from a catalogue of 1-D radiative transfer experiments. In case of overlap of these spheres, photons are conserved by redistributing them around the connected ionized regions corresponding to the spheres. The efficiency with which these maps are created allows us to span the large parameter space typically encountered in reionization simulations. We compare our results with other, more accurate, 3-D radiative transfer simulations and find excellent agreement for the redshifts and the spatial scales of interest to upcoming 21cm experiments. We generate a contiguous observational cube spanning redshift 6 to 12 and use these simulations to study the differences in the reionization histories between stars and quasars. Finally, the signal is convolved with the LOFAR beam response and its effects are analyzed and quantified. Statistics performed on this mock data set shed light on possible observational strategies for LOFAR.

Future high redshift 21-cm experiments will suffer from a high degree of contamination, due both to astrophysical foregrounds and to non-astrophysical and instrumental effects. In order to reliably extract the cosmological signal from the observed data, it is essential to understand very well all data components and their influence on the extracted signal. Here we present simulated astrophysical foregrounds datacubes and discuss their possible statistical effects on the data. The foreground maps are produced assuming 5 deg x 5 deg windows that match those expected to be observed by the LOFAR Epoch-of-Reionization (EoR) key science project. We show that with the expected LOFAR-EoR sky and receiver noise levels, which amount to ~52 mK at 150 MHz after 300 hours of total observing time, a simple polynomial fit allows a statistical reconstruction of the signal. We also show that the polynomial fitting will work for maps with realistic yet idealised instrument response, i.e., a response that includes only a uniform uv coverage as a function of frequency and ignores many other uncertainties. Polarized galactic synchrotron maps that include internal polarization and a number of Faraday screens along the line of sight are also simulated. The importance of these stems from the fact that the LOFAR instrument, in common with all current interferometric EoR experiments has an instrumentally polarized response.

We use a sample of 53 massive early-type strong gravitational lens galaxies with well-measured redshifts (ranging from z=0.06 to 0.36) and stellar velocity dispersions (between 175 and 400 km/s) from the Sloan Lens ACS (SLACS) Survey to derive numerous empirical scaling relations. The ratio between central stellar velocity dispersion and isothermal lens-model velocity dispersion is nearly unity within errors. The SLACS lenses define a fundamental plane (FP) that is consistent with the FP of the general population of early-type galaxies. We measure the relationship between strong-lensing mass M_lens within one-half effective radius (R_e/2) and the dimensional mass variable M_dim = G^-1 sigma_e2^2 R_e/2 to be log_10 [M_lens/10^11 M_Sun] = (1.03 +/- 0.04) log_10 [M_dim/10^11 M_Sun] + (0.54 +/- 0.02) (where sigma_e2 is the projected stellar velocity dispersion within R_e/2). The near-unity slope indicates that the mass-dynamical structure of massive elliptical galaxies is independent of mass, and that the "tilt" of the SLACS FP is due entirely to variation in total (luminous plus dark) mass-to-light ratio with mass. Our results imply that dynamical masses serve as a good proxies for true masses in massive elliptical galaxies. Regarding the SLACS lenses as a homologous population, we find that the average enclosed 2D mass profile goes as log_10 [M(<R)/M_dim] = (1.10 +/- 0.09) log_10 [R/R_e] + (0.85 +/- 0.03), consistent with an isothermal (flat rotation curve) model when de-projected into 3D. This measurement is inconsistent with the slope of the average projected aperture luminosity profile at a confidence level greater than 99.9%, implying a minimum dark-matter fraction of f_DM = 0.38 +/- 0.07 within one effective radius. (abridged)

We present the definitive data for the full sample of 131 strong gravitational lens candidates observed with the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope by the Sloan Lens ACS (SLACS) Survey. All targets were selected for higher-redshift emission lines and lower-redshift continuum in a single Sloan Digital Sky Survey (SDSS) spectrum. The foreground galaxies are primarily of early-type morphology, with redshifts from approximately 0.05 to 0.5 and velocity dispersions from 160 km/s to 400 km/s; the faint background emission-line galaxies have redshifts ranging from about 0.2 to 1.2. We confirm 70 systems showing clear evidence of multiple imaging of the background galaxy by the foreground galaxy, as well as an additional 19 systems with probable multiple imaging. For 63 clear lensing systems, we present singular isothermal ellipsoid and light-traces-mass gravitational lens models fitted to the ACS imaging data. These strong-lensing mass measurements are supplemented by magnitudes and effective radii measured from ACS surface-brightness photometry and redshifts and velocity dispersions measured from SDSS spectroscopy. These data constitute a unique resource for the quantitative study of the inter-relations between mass, light, and kinematics in massive early-type galaxies. We show that the SLACS lens sample is statistically consistent with being drawn at random from a parent sample of SDSS galaxies with comparable spectroscopic parameters and effective radii, suggesting that the results of SLACS analyses can be generalized to the massive early-type population.

The dynamic spectrum of a radio pulsar is an in-line digital hologram of the ionised interstellar medium. It has previously been demonstrated that such holograms permit image reconstruction, in the sense that one can determine an approximation to the complex electric field values as a function of Doppler-shift and delay, but to date the quality of the reconstructions has been poor. Here we report a substantial improvement in the method which we have achieved by simultaneous optimisation of the thousands of coefficients that describe the electric field. For our test spectrum of PSR B0834+06 we find that the model provides an accurate representation of the data over the full 63 dB dynamic range of the observations: residual differences between model and data are noise-like. The advent of interstellar holography enables detailed quantitative investigation of the interstellar radio-wave propagation paths for a given pulsar at each epoch of observation; we illustrate this using our test data which show the scattering material to be structured and highly anisotropic. The temporal response of the medium exhibits a scattering tail out to beyond 100 microsec and a pulse arrival time measurement at this frequency and this epoch of observation would be affected by a mean delay of 15 microsec due to multipath propagation in the interstellar medium.

We present the first results of a combined VLT VIMOS integral-field unit and Hubble Space Telescope (HST)/ACS study of the early-type lens galaxy SDSS J2321-097 at z=0.0819, extending kinematic studies to a look-back time of 1 Gyr. This system, discovered in the Sloan Lens ACS Survey (SLACS), has been observed as part of a VLT Large Programme with the goal of obtaining two-dimensional stellar kinematics of 17 early-type galaxies to z~0.35 and Keck spectroscopy of an additional dozen lens systems. Bayesian modelling of both the surface brightness distribution of the lensed source and the two-dimensional measurements of velocity and velocity dispersion has allowed us, under the only assumptions of axisymmetry and a two-integral stellar distribution function (DF) for the lens galaxy, to dissect this galaxy in three dimensions and break the classical mass--anisotropy, mass-sheet and inclination--oblateness degeneracies. Our main results are that the galaxy (i) has a total density profile well described by a single power-law rho propto r^{-gamma'} with gamma'=2.06^{+0.03}_{-0.06}; (ii) is a very slow rotator (specific stellar angular momentum parameter lambda_R = 0.075); (iii) shows only mild anisotropy (delta ~ 0.15); and (iv) has a dark matter contribution of ~30 per cent inside the effective radius. Our first results from this large combined imaging and spectroscopic effort with the VLT, Keck and HST show that the structure of massive early-type galaxies beyond the local Universe can now be studied in great detail using the combination of stellar kinematics and gravitational lensing. Extending these studies to look-back times where evolutionary effects become measurable holds great promise for the understanding of formation and evolution of early-type galaxies.

We report the discovery of two concentric Einstein rings around the gravitational lens SDSSJ0946+1006, as part of the Sloan Lens ACS Survey. The main lens is at redshift zl=0.222, while the inner ring (1) is at zs1=0.609 and Einstein radius $Re_1=1.43\pm0.01"$. The wider image separation ($Re_2=2.07\pm 0.02"$) of the outer ring (2) implies that it is at higher redshift. Its detection in the F814W filter implies zs2<6.9. The configuration can be well described by a total density profile $\rho_{tot} ~ r^-g'$ with $g'=2.00\pm0.03$ and velocity dispersion $\sigma_{SIE}=287\pm5\kms$. [...] We consider whether this configuration can be used to constrain cosmological parameters exploiting angular distance ratios entering the lens equations. Constraints for SDSSJ0946+1006, are uninteresting due to the sub-optimal lens and source redshifts. We then consider the perturbing effect of the mass associated with Ring 1 building a double lens plane compound lens model. This introduces minor changes to the mass of the main lens and allows to estimate the mass of Ring 1 $(\sigma_{SIE,s1}=94\pm30\kms)$. We examine the prospects of doing cosmography with a sample of 50 double lenses, expected from future space based surveys such as DUNE or JDEM. Taking full account of the model uncertainties, such a sample could be used to measure $\Omega_m$ and $w$ with 10% accuracy, for a flat cosmology.

The gravitational lens system CLASS B2108+213 has two lensed images separated by 4.56 arcsec. Such a wide image separation suggests that the lens is either a massive galaxy, or is composed of a group of galaxies. To investigate the structure of the lensing potential we have carried out new high resolution imaging of the two lensed images at 1.7 GHz with the VLBA and at 5 GHz with global VLBI. Compact and extended emission is detected from the two lensed images, which provides additional constraints to the lensing mass model. We find that the data are consistent with either a single lensing galaxy, or a two galaxy lens model that takes account of a nearby companion to the main lensing galaxy within the Einstein radius of the system. However, for an ensemble of global power-law mass models, those with density profiles steeper than isothermal are a better fit. The best-fitting profile for a single spherical mass model has a slope of $\gamma=$~2.45$_{-0.18}^{+0.19}$. The system also has a third radio component which is coincident with the main lensing galaxy. This component is detected at milli-arcsecond scales for the first time by the 1.7 GHz VLBA and 5 GHz global VLBI imaging. However, the third radio component is found not to be consistent with a core lensed image because the radio spectrum differs from the two lensed images, and its flux-density is too high when compared to what is expected from simple mass models with a variable power-law density profile and/or a reasonable core radius. Furthermore, 1.4 GHz imaging of the system with the MERLIN finds extended lobe emission on either side of the main lensing galaxy. Therefore, the radio emission from the third radio component is almost certainly from an AGN within the main lensing galaxy, which is classified as an FR I type radio source.

We combine high-resolution images in four optical/infra-red bands, obtained with the laser guide star adaptive optics system on the Keck Telescope and with the Hubble Space Telescope, to study the gravitational lens system SDSSJ0737+3216 (lens redshift 0.3223, source redshift 0.5812). We show that (under favorable observing conditions) ground-based images are comparable to those obtained with HST in terms of precision in the determination of the parameters of both the lens mass distribution and the background source. We also quantify the systematic errors associated with both the incomplete knowledge of the PSF, and the uncertain process of lens galaxy light removal, and find that similar accuracy can be achieved with Keck LGSAO as with HST. We then exploit this well-calibrated combination of optical and gravitational telescopes to perform a multi-wavelength study of the source galaxy at 0.01" effective resolution. We find the Sersic index to be indicative of a disk-like object, but the measured half-light radius (0.59+-0.007+-0.1 kpc) and stellar mass (2.0+-1.0+-0.8e9Msun) place it more than three sigma away from the local disk size-mass relation. The SDSSJ0737+3216 source has the characteristics of the most compact faint blue galaxies studied, and has comparable size and mass to dwarf early-type galaxies in the local universe. With the aid of gravitational telescopes to measure individual objects' brightness profiles to 10% accuracy, the study of the high-redshift size-mass relation may be extended by an order of magnitude or more beyond existing surveys at the low-mass end, thus providing a new observational test of galaxy formation models.

Gravitational lensing and stellar dynamics are two independent methods, based solely on gravity, to study the mass distributions of galaxies. Both methods suffer from degeneracies, however, that are difficult to break. In this paper, we present a new framework that self-consistently unifies gravitational lensing and stellar dynamics. This approach breaks some of classical degeneracies that have limited their individual usage, in particular in the study of high-redshift galaxies. The methodology is based on the premise that, for any given galaxy potential, the mapping of both the unknown lensed source brightness distribution and the stellar phase-space distribution function on to the photometric and kinematic observables, can be cast as a single set of coupled linear equations. This set of linear equations is solved, maximizing the likelihood penalty function. The evidence penalty function, as derived from Bayesian statistics, subsequently allows the best potential-model parameters to be found and potential-model families, or other model assumptions (e.g. PSF), to be quantitatively ranked. We have implemented a fast algorithm that solves for the maximum-likelihood pixelized lensed source brightness distribution and the two-integral stellar phase-space distribution function f(E, L_z), assuming axisymmetric potentials. To make the method practical, we have devised a new Monte-Carlo approach to Schwarzschild's orbital superposition method, based on the superposition of two-integral (E and L_z) toroidal components, to find the maximum-likelihood two-integral distribution function in a matter of seconds in any axisymmetric potential. The non-linear parameters of the potential are subsequently found through a hybrid MCMC and Simplex optimization of the evidence. (Abridged)

We combine strong-lensing masses with SDSS stellar velocity dispersions and HST-ACS effective (half-light) radii for 36 lens galaxies from the Sloan Lens ACS (SLACS) Survey to study the mass dependence of mass-dynamical structure in early-type galaxies. We find that over a 180--390 km/s range in velocity dispersion, structure is independent of lensing mass to within 5%. This result suggests a systematic variation in the total (i.e., luminous plus dark matter) mass-to-light ratio as the origin of the tilt of the fundamental plane (FP) scaling relationship between galaxy size, velocity dispersion, and surface brightness. We construct the FP of the lens sample, which we find to be consistent with the FP of the parent SDSS early-type galaxy population, and present the first observational correlation between mass-to-light ratio and residuals about the FP. Finally, we re-formulate the FP in terms of surface mass density rather than surface brightness. By removing the complexities of stellar-population effects, this mass-plane formulation will facilitate comparison to numerical simulations and possible use as a cosmological distance indicator.

We present a weak gravitational lensing analysis of 22 early-type strong lens galaxies, based on deep HST images obtained as part of the Sloan Lens ACS Survey. Using the most advanced techniques to control systematic uncertainties related to the variable PSF and charge transfer efficiency of the ACS, we detect weak lensing signal out to 300 kpc/h. We analyze blank control fields from the COSMOS survey in the same manner, inferring that the residual systematic uncertainty in the tangential shear is <0.3%. A joint strong and weak lensing analysis shows that the average total mass density profile is consistent with isothermal over two decades in radius (3-300 kpc/h, approximately 1-100 Reff). This finding extends by over an order of magnitude in radius previous results, based on strong lensing and/or stellar dynamics, that luminous and dark component ``conspire'' to form an isothermal mass distribution. In order to disentangle the contributions of luminous and dark matter, we fit a two-component mass model (R^1/4 + NFW) to the weak and strong lensing constraints. It provides a good fit to the data with only two free parameters; i) the average stellar mass-to-light ratio M_*/L_V=4.48 +- 0.46 hMo/Lo, in agreement with that expected for an old stellar population; ii) the average virial mass-to-light ratio M_vir/L_V = 246+101-87 hMo/Lo. [abridged]

The existence of flux-ratio anomalies between fold and cusp images in galaxy-scale strong-lens systems has led to an interpretation based on the presence of a high mass-fraction of cold-dark-matter (CDM) substructures around galaxies, as predicted by numerical N-body simulations. The flux-ratio anomaly is particularly evident in the radio-loud quadruple gravitational lens system CLASS B2045+265. In this paper, new high-resolution radio, optical, and infrared imaging of B2045+265 is presented which sheds more light on this anomaly and its possible causes. First, deep Very Long Baseline Array observations show very compact images, possibly with a hint of a jet, but with no evidence for differential scattering or scatter broadening. Second, optical and infrared observations with the Hubble Space Telescope and through Adaptive-Optics imaging with the W. M. Keck Telescope, show a previously undiscovered object -- interpreted as a (tidally disrupted) dwarf satellite based on its colours and slight extension -- between the main lens galaxy and the three anomalous flux-ratio images. Third, colour variations in the early-type lens galaxy indicate recent star-formation, possibly the result of secondary infall of gas-rich satellites. A population of such galaxies around the lens system could explain the previously discovered strong [O II] emission. However, spiral structure and/or normal star formation in the lens galaxy cannot be excluded. In light of these new data, we propose a lens model for the system, including the observed dwarf satellite, which reproduces all positional and flux-ratio constraints, without the need for additional CDM substructure. [abridged]

LOFAR, the Low Frequency Array, is a next-generation radio telescope that is being built in Northern Europe and expected to be fully operational at the end of this decade. It will operate at frequencies from 15 to 240 MHz (corresponding to wavelengths of 20 to 1.2 m). Its superb sensitivity, high angular resolution, large field of view and flexible spectroscopic capabilities will represent a dramatic improvement over previous facilities at these wavelengths. As such, LOFAR will carry out a broad range of fundamental astrophysical studies. The design of LOFAR has been driven by four fundamental astrophysical applications: (i) The Epoch of Reionisation, (ii) Extragalactic Surveys and their exploitation to study the formation and evolution of clusters, galaxies and black holes, (iii) Transient Sources and their association with high energy objects such as gamma ray bursts, and (iv) Cosmic Ray showers and their exploitation to study the origin of ultra-high energy cosmic rays. In this conference the foreseen LOFAR work on the epoch of reionisation has been covered by de Bruyn and on cosmic ray showers by Falcke. During this contribution we will first present the LOFAR project with an emphasis on the challenges faced when carrying out sensitive imaging at low radio frequencies. Subsequently, we will discuss LOFAR's capabilities to survey the low-frequency radio sky. Main aims for the planned surveys are studies of z>6 radio galaxies, diffuse emission associated with distant clusters and starbursting galaxies at z>2.

Previous authors have reported the detection of intrinsically faint sub-mm emission lensed by the cluster MS0451.6-0305. They suggest that this emission arises from a merging system composed of a Ly-break galaxy and a pair of extremely red objects which are multiply-imaged in the optical/NIR observations. Since the submm emission presents an unusually large angular extent (~1 arcmin), the possible radio emission asociatted with that system can help to identify optical/NIR counterparts due to the higher spatial resolution and astrometric accuracy of the radio observations. Archive VLA data (BnA configuration at 1.4 GHz) was reduced and analysed. A simple lens model was constructed to aid the interpretation of the radio and pre-existing sub-mm and optical/NIR data. We present a 1.4 GHz map of the central region of MS0451.6-0305 and report the detection of gravitationally lensed radio emission, coincident with the previously discovered sub-mm lensed emission. The overall morphology and scale of the radio and sub-mm emission are strikingly similar, extending ~1 arcmin across the sky. This observation strongly suggests that the radio and sub-mm emission arise from the same sources. Preliminary estimates of the total S_850 microns/S_1.4 GHz flux density ratio appear to be consistent with that expected from distant star forming galaxies. The radio emission is resolved into 7 distinct components, and the overall structure can be explained, using a simple lens model, with three multiply-imaged radio sources at z~2.9. One of these sources is predicted to lie in the middle of the previously mentioned system in the source plane, suggesting that it is related to the intense star formation generated during the merging process.

LOFAR (Low Frequency Array) is an innovative radio telescope optimized for the frequency range 30-240 MHz. The telescope is realized as a phased aperture array without any moving parts. Digital beam forming allows the telescope to point to any part of the sky within a second. Transient buffering makes retrospective imaging of explosive short-term events possible. The scientific focus of LOFAR will initially be on four key science projects (KSPs): 1) detection of the formation of the very first stars and galaxies in the universe during the so-called epoch of reionization by measuring the power spectrum of the neutral hydrogen 21-cm line (Shaver et al. 1999) on the ~5' scale; 2) low-frequency surveys of the sky with of order $10^8$ expected new sources; 3) all-sky monitoring and detection of transient radio sources such as gamma-ray bursts, x-ray binaries, and exo-planets (Farrell et al. 2004); and 4) radio detection of ultra-high energy cosmic rays and neutrinos (Falcke & Gorham 2003) allowing for the first time access to particles beyond 10^21 eV (Scholten et al. 2006). Apart from the KSPs open access for smaller projects is also planned. Here we give a brief description of the telescope.

We derive Fundamental Plane parameters of 15 early-type lens galaxies identified by the Sloan Lens ACS (SLACS) Survey. The size of the sample allows us to investigate for the first time the distribution of lens galaxies in the FP space. After correcting for evolution, we find that lens galaxies occupy a subset of the local FP. The edge-on projection (approximately M vs M/L) is indistinguishable from that of normal early-type galaxies. However -- within the fundamental plane -- the lens galaxies appear to concentrate at the edge of the region populated by normal early-type galaxies. We show that this is a result of our selection procedure (approximately velocity dispersion sigma>240km/s). We conclude that SLACS lenses are a fair sample of high velocity dispersion early-type galaxies. By comparing the central stellar velocity dispersion that of the best fit lens model, we find <f_SIE>=<sigma/sigma_SIE>= =1.01+-0.02 with 0.065 rms scatter. We conclude that within the Einstein radii the SLACS lenses are very well approximated by isothermal ellipsoids, requiring a fine tuning of the stellar and dark matter distribution (bulge-halo ``conspiracy''). Interpreting the offset from the local FP in terms of evolution of the stellar mass-to-light ratio, we find for the SLACS lenses d log M/L_B/dz=-0.69+-0.08 (rms 0.11) consistent with the rate found for field early-type galaxies and with a scenario where most of the stars were formed at high redshift (>2) with secondary episodes of star formation providing less than ~10% of the stellar mass below z=1. We discuss star formation history and structural homogeneity in the context of formation mechanisms such as collisionless (``dry'') mergers. [Abridged]

(Abridged) We present a joint gravitational lensing and stellar dynamical analysis of fifteen massive field early-type galaxies, selected from the Sloan Lens (SLACS) Survey. The following numerical results are found: (i) A joint-likelihood gives an average logarithmic density slope for the total mass density of 2.01 (+0.02/-0.03) (68 perecnt C.L). inside the Einstein radius. (ii) The average position-angle difference between the light distribution and the total mass distribution is found to be 0+-3 degrees, setting an upper limit of <= 0.035 on the average external shear. (iii) The average projected dark-matter mass fraction is inferred to be 0.25+-0.06 inside R_E, using the stellar mass-to-light ratios derived from the Fundamental Plane as priors. (iv) Combined with results from the LSD Survey, we find no significant evolution of the total density slope inside one effective radius: a linear fit gives d\gamma'/dz = 0.23+-0.16 (1-sigma) for the range z=0.08-1.01. The small scatter and absence of significant evolution in the inner density slopes suggest a collisional scenario where gas and dark matter strongly couple during galaxy formation, leading to a total mass distribution that rapidly converge to dynamical isothermality.

We report the serendipitous discovery of two strong gravitational lens candidates (ACS J160919+6532 and ACS J160910+6532) in deep images obtained with the Advanced Camera for Surveys on the Hubble Space Telescope, each less than 40 arcsec from the previously known gravitational lens system CLASS B1608+656. The redshifts of both lens galaxies have been measured with Keck and Gemini: one is a member of a small galaxy group at z~0.63, which also includes the lensing galaxy in the B1608+656 system, and the second is a member of a foreground group at z~0.43. By measuring the effective radii and surface brightnesses of the two lens galaxies, we infer their velocity dispersions based on the passively evolving Fundamental Plane (FP) relation. Elliptical isothermal lens mass models are able to explain their image configurations within the lens hypothesis, with a velocity dispersion compatible with that estimated from the FP for a reasonable source-redshift range. Based on the large number of massive early-type galaxies in the field and the number-density of faint blue galaxies, the presence of two additional lens systems around CLASS B1608+656 is not unlikely in hindsight. Gravitational lens galaxies are predominantly early-type galaxies, which are clustered, and the lensed quasar host galaxies are also clustered. Therefore, obtaining deep high-resolution images of the fields around known strong lens systems is an excellent method of enhancing the probability of finding additional strong gravitational lens systems.

The Sloan Lens ACS (SLACS) Survey is an efficient Hubble Space Telescope Snapshot imaging survey for new galaxy-scale strong gravitational lenses. The targeted lens candidates are selected spectroscopically from within the Sloan Digital Sky Survey (SDSS) database of galaxy spectra for having multiple nebular emission lines at a redshift significantly higher than that of the SDSS target galaxy. In this paper, we present a catalog of 19 newly discovered gravitational lenses, along with 9 other observed candidate systems that are either possible lenses, non-lenses, or non-detections. The survey efficiency is thus >=68%. We also present Gemini and Magellan IFU data for 9 of the SLACS targets, which further support the lensing interpretation. A new method for the effective subtraction of foreground galaxy images to reveal faint background features is presented. We show that the SLACS lens galaxies have colors and ellipticities typical of the spectroscopic parent sample from which they are drawn (SDSS luminous red galaxies and quiescent main-sample galaxies), but are somewhat brighter and more centrally concentrated. Several explanations for the latter bias are suggested. The SLACS survey provides the first statistically significant and homogeneously selected sample of bright early-type lens galaxies, furnishing a powerful probe of the structure of early-type galaxies within the half-light radius. The high confirmation rate of lenses in the SLACS survey suggests consideration of spectroscopic lens discovery as an explicit science goal of future spectroscopic galaxy surveys (abridged).

Strong gravitational lensing and stellar dynamics provide two complementary and orthogonal constraints on the density profiles of galaxies. Based on spherically symmetric, scale-free, mass models, it is shown that the combination of both techniques is powerful in breaking the mass-sheet and mass-anisotropy degeneracies. Second, observational results are presented from the Lenses Structure & Dynamics (LSD) Survey and the Sloan Lens ACS (SLACS) Survey collaborations to illustrate this new methodology in constraining the dark and stellar density profiles, and mass structure, of early-type galaxies to redshifts of unity.

We propose the novel method of ``gravitational imaging'' to detect and quantify luminous and dark-matter substructure in gravitational-lens galaxies. The method utilizes highly-magnified Einstein rings and arcs as sensitive probes of small perturbations in the lens potential (due to the presence of mass substructure), reconstructing the gravitational lens potential non-parametrically. Numerical simulations show that the implemented algorithm can reconstruct the smooth mass distribution of a typical lens galaxy - exhibiting reasonable signal-to-noise Einstein rings - as well as compact substructure with masses as low as M_sub~10^-3 M_lens, if present. ``Gravitational imaging'' of pure dark-matter substructure around massive galaxies can provide a new window on the standard cold-dark-matter paradigm, using very different physics than ground-based direct-detection experiments, and probe the hierarchical structure-formation model which predicts this substructure to exist in great abundance.

We report the discovery of a new gravitational lens system from the CLASS survey, CLASS B0631+519. VLA, MERLIN and VLBA observations show a doubly-imaged radio core, a doubly-imaged lobe and a second lobe that is probably quadruply-imaged. The maximum image separation is 1.16 arcseconds. The VLBA resolves the most magnified image of the flat-spectrum radio core into a number of sub-components spread across approximately 20 milli-arcseconds. Optical and near-infrared imaging with the ACS and NICMOS cameras on the HST show that there are two galaxies along the line of sight to the lensed source, as was previously discovered by optical spectroscopy. The nearer galaxy at z=0.0896 is a small blue irregular, while the more distant galaxy at z=0.6196 is an elliptical type and appears to contribute most of the lensing effect. The host galaxy of the lensed source is visible in NICMOS imaging as a set of arcs that form an almost complete Einstein ring. Mass modelling using non-parametric techniques can reproduce the ring and indicates that the irregular galaxy has a (localised) effect on the flux density distribution in the Einstein ring at the 5-10% level.

Gravitational lensing and stellar dynamics provide two complementary, nearly orthogonal, constraints on the mass distribution of early-type lens galaxies. This allows the luminous and dark-matter distribution in higher-redshift (z>0.1) galaxies to be studied beyond the limitations of each individual method. Two surveys have been initiated to compile a large sample of early-type galaxies suitable to lensing and dynamical studies: (1) The Lenses Structure & Dynamics (LSD) Survey and (2) the Sloan Lens ACS (SLACS) Survey. Using spherically symmetric mass models, I illustrated how lensing and dynamical constraints can be used to measure the ``effective'' density slope (gamma') of galaxies inside their Einstein radii and estimate the typical error on this determination. The main results from the LSD survey thus far are: (i) Massive (typically >L_*) early-type galaxies at z=0.5-1 contain a significant fraction f_CDM=0.4-0.7 of dark matter inside their Einstein radii. [The null-hypothesis, f_CDM=0, is excluded at the >99% in all analyzed systems.] (ii) The inner CDM density slope is gamma_CDM=1.3^{+0.2}_{-0.3} (68% CL). (iii) The total density slope gamma' = 1.9 +- 0.1 (with 0.3 rms scatter in the sample). The intrinsic scatter of 15% in gamma' is consistent with local dynamical studies and can lead to a 30% rms scatter in inferred values of H_0 from lens time-delays, when purely isothermal mass models are assumed. Hence, the common practice to assume that lens galaxies are perfectly isothermal should be abandoned, especially in cases where this assumption is critical.

We present observations of CLASS B2108+213, the widest separation gravitational lens system discovered by the Cosmic Lens All-Sky Survey. Radio imaging using the VLA at 8.46 GHz and MERLIN at 5 GHz shows two compact components separated by 4.56 arcsec with a faint third component in between which we believe is emission from a lensing galaxy. 5-GHz VLBA observations reveal milliarcsecond-scale structure in the two lensed images that is consistent with gravitational lensing. Optical emission from the two lensed images and two lensing galaxies within the Einstein radius is detected in Hubble Space Telescope imaging. Furthermore, an optical gravitational arc, associated with the strongest lensed component, has been detected. Surrounding the system are a number of faint galaxies which may help explain the wide image separation. A plausible mass distribution model for CLASS B2108+213 is also presented.

(abridged) We present Gemini integral-field unit (IFU) spectroscopy and Hubble Space Telescope (HST) F435W- and F814W-band images of a new four-image gravitational lens, SDSSJ140228.22+632133.3, obtained as part of an HST Snapshot program designed to expand the sample of known gravitational lenses amenable to detailed photometric, lensing, and dynamical studies. The lens is a smooth elliptical galaxy at a redshift of z_l = 0.2046 +/- 0.0001 with a Sloan r-band magnitude of 17.00 +/- 0.05 and a stellar velocity dispersion of 267 +/- 17 km s^-1, obtained from its SDSS spectrum. Multiple emission lines place the quadruply-imaged source at a redshift of z_s = 0.4814 +/- 0.0001. The best-fitting singular isothermal ellipsoid lens model gives an Einstein radius b = 1.35" +/- 0.05" (or [4.9 +/- 0.2] h_65^-1 kpc), corresponding to a total mass of (30.9 +/- 2.3) x 10^10 h_65^-1 M_sun within the critical curve. In combination with HST photometry this gives a rest-frame B-band mass-to-light ratio of (8.1 +/- 0.7) h_65 times solar within the same region. The lens model predicts a luminosity-weighted stellar dispersion within the 3"-diameter SDSS aperture of sigma_* ~= 270 km s^-1, in good agreement with the observed value. Using the model to de-lens the four lensed images yields a source with a smooth, monotonically-decreasing brightness distribution.

The advent of new observational facilities in the last two decades has allowed the rapid discovery and high-resolution optical imaging of many strong lens systems from galaxy to cluster scales, as well as their spectroscopic follow-up. Radio telescopes have played the dominant role in the systematic detection of dozens of new arcsec-scale lens systems. For the future, we expect nothing less! The next major ground- and space-based facilities, especially the Square Kilometer Array can discover tens of thousands of new lens systems in large sky surveys. For optical imaging and spectroscopic follow-up a strong synergy with planned optical facilities is needed. Here, we discuss the field where strong gravitational lensing is expected to play the dominant role and where SKA can have a major impact: The study of the internal mass structure and evolution of galaxies and clusters to z~1. In addition, studies of more exotic phenomena are contemplated. For example, milli- and microlensing can provide a way to measure the mass-functions of stars and CDM substructure at cosmological distances. All-sky radio monitoring will also rapidly develop the field of time-domain lensing.

Compact groups of galaxies recently have been discovered in association with several strong gravitational lens systems. These groups provide additional convergence to the lensing potential and thus affect the value of H_0 derived from the systems. Lens system time delays are now being measured with uncertainties of only a few percent or better. Additionally, vast improvements are being made in incorporating observational constraints such as Einstein ring structures and stellar velocity dispersions into the lens models. These advances are reducing the uncertainties on H_0 to levels at which the the effects of associated galaxy groups may contribute significantly to the overall error budget. We describe a dedicated multiwavelength program, using Keck, HST, and Chandra, to find such groups and measure their properties. We present, as a case study, results obtained from observations of the CLASS lens system B1608+656 and discuss the implications for the value of H_0 derived from this system.

We report the discovery of a transient source in the central regions of galaxy cluster Abell 267. The object, which we call "PALS1", was found in a survey aimed at identifying highly-magnified Lyman-break galaxies in the fields of intervening rich clusters. At discovery, the source had U>24.7 (2-sigma; AB), g=21.96, and very blue g-r and r-i colors; i.e., PALS1 was a U-band drop-out, characteristic of star-forming galaxies and quasars at z~3. However, three months later the source had faded by more than three magnitudes. Further observations showed a continued decline in luminosity, to R>26.4 seven months after discovery. Though the apparent brightness is suggestive of a supernova at roughly the cluster redshift, we show that the photometry and light curve argue against any known type of supernova at any redshift. The spectral energy distribution and location near the center of a galaxy cluster are consistent with the hypothesis that PALS1 is a gravitationally-lensed transient at z~3.3. If this interpretation is correct, the source is magnified by a factor of 4-7 and two counterimages are predicted. Our lens model predicts time delays between the three images of 1-10 years and that we have witnessed the final occurrence of the transient. The intense luminosity (M(AB) ~ -23.5 after correcting for lensing) and blue UV continuum (implying T>50,000 K) argue the source may have been a flare resulting from the tidal disruption of a star by a 10^6-10^8 solar-mass black hole. Regardless of its physical nature, PALS1 highlights the importance of monitoring regions of high magnification in galaxy clusters for distant, time-varying phenomena.

We present a study of the lens properties of quadruply imaged systems, lensed by numerically simulated galaxies. We investigate a simulated elliptical and disc galaxy drawn from high resolution simulations of galaxy formation in a concordance $\Lambda$CDM universe. The simulations include the effects of gas dynamics, star formation and feedback processes. Flux-ratio anomalies observed in strong gravitational lensing potentially provide an indicator for the presence of mass substructure in lens galaxies as predicted from CDM simulations. We particularly concentrate on the prediction that, for an ideal cusp caustic, the sum of the signed magnifications of the three highly magnified images should vanish when the source approaches the cusp. Strong violation of this cusp relation indicates the presence of substructure, regardless of the global, smooth mass model of the lens galaxy. We draw the following conclusions: (1) the level of substructure present in simulations produces violations of the cusp relation comparable to those observed, (2) higher-order catastrophes (e.g. swallowtails) can also cause changes of the order of 0.6 in the cusp relation as predicted by a smooth model, (3) the flux anomaly distribution depends on the image parity and flux and both the brightest minimum and saddle-point images are more affected by substructure than the fainter images. In addition, the brightest saddle point is demagnified w.r.t. the brightest minimum. Our results are fully numerical and properly include all mass scales, without making semi-analytic assumptions. They are ultimately limited by the mass resolution of single particles in the simulation determined by current computational limits, however show that our results are not affected by shot-noise due to the finite number of particles.

(Shortened) The combination of lensing and stellar dynamics breaks the mass-anisotropy degeneracy and provides stringent constraints on the mass distribution in early-type (E/S0) galaxies out to z~1. We present the combined results from the five field E/S0 lens galaxies at z=0.5-1.0 analyzed as part of the LSD Survey. We find: (i) Constant M/L models are ruled out at >99% CL for all five E/S0s. The projected dark-matter mass fractions inside the Einstein (effective) radius is f_DM=0.37-0.72 (0.15-0.65) for isotropic models. (ii) The average power-law slope of the total mass distribution is <\gamma'>=1.75+-0.10 for isotropic models with 0.20 rms scatter. The ratio between the observed central stellar velocity dispersion and that from the best-fit SIE lens model is <f_SIE>=<\sigma/\sigma_SIE>=0.87+-0.04 with 0.08 rms. Considering that \gamma'>2 and f_SIE>1 have been reported for other systems, we conclude that there is a significant intrinsic scatter in the density slopes of E/S0s (rms \~15%). Hence, the isothermal approximation is not sufficiently accurate for applications that depend critically on the slope of the mass density profile (i.e. measuring H_0). (iii) The inner power-law slope of the dark-matter halo is constrained to be <\gamma>=1.3(+0.2/-0.4) (68% CL) for the isotropic model or an upper limit of \gamma<0.6, if the galaxies are radially anisotropic (r_i=R_e). This is consistent with numerical simulations only for an isotropic velocity ellipsoid and if baryonic collapse and star-formation do not steepen dark-matter density profiles. (iv) The average stellar M/L evolves as d\log(M_*/L_B)/dz =-0.72+-0.10, obtained via the FP. Based on lensing and dynamics we find d\log(M_*/L_B)/dz=-0.75+-0.17, indicating that the M/L ratio evolution for our sample of field E/S0s is faster than those in clusters.

We present the optical spectra of four newly discovered gravitational lenses from the Cosmic Lens All-Sky Survey (CLASS). These observations were carried out using the Low Resolution Imaging Spectrograph on the W. M. Keck-I Telescope as part of a program to study galaxy-scale gravitational lenses. From our spectra we found the redshift of the background source in CLASS B0128+437 (z_s=3.1240+-0.0042) and the lensing galaxy redshifts in CLASS B0445+123 (z_l=0.5583+-0.0003) and CLASS B0850+054 (z_l=0.5883+-0.0006). Intriguingly, we also discovered that CLASS B0631+519 may have two lensing galaxies (z_l,1=0.0896+-0.0001, z_l,2=0.6196+-0.0004). We also found a single unidentified emission line from the lensing galaxy in CLASS B0128+437 and the lensed source in CLASS B0850+054. We find the lensing galaxies in CLASS B0445+123 and CLASS B0631+519 (l,2) to be early-type galaxies with Einstein Radii of 2.8-3.0 h^{-1} kpc. The deflector in CLASS B0850+054 is a late-type galaxy with an Einstein Radius of 1.6 h^{-1} kpc.

We present a kinematic model of the neutral-hydrogen in the spiral galaxy of the lens system PKS1830-211, based on a Multi Element Radio-Linked Interferometer (MERLIN) 1.4-GHz radio map and the integrated and redshifted 21-cm hydrogen absorption-line profile as measured with the Westerbork Synthesis Radio Telescope (WSRT). Degeneracies in the models do not allow a unique determination of the kinematic center and forthcoming deeper Hubble Space Telescope observations with the Advanced Camera for Surveys (ACS) are required to break this degeneracy. Even so, we measure the inclination of the hydrogen disk: i=(17-32) deg., indicating a close-to face-on spiral galaxy. The optical depth increases with radius over the extent of the Einstein ring, suggesting HI depletion towards the lens center. The latter could be due to star formation or conversion of HI in to molecular hydrogen because of a higher metalicity/dust content in the galaxy center. The neutral hydrogen optical depth gives N_HI= 2x10^21 cm^-2 at r=5.0 h_{70}^-1 kpc in the disk (T_s=100K), comparable to local spiral galaxies. Our study shows that planned new radio telescopes (i.e. ALMA, LOFAR and SKA) -- which will discover large numbers of similar lens systems -- are powerful new tools to probe the internal structure, kinematics and evolution of spiral galaxies in detail to z>=1, as well as their neutral-hydrogen (and molecular) content, thereby complementing studies of HI emission from spirals at z<=1 and damped Ly-alpha systems to z>>1.

We describe the first results from two observational projects aimed at measuring the amount and spatial distribution of dark matter in distant early-type galaxies (E/S0s) and clusters of galaxies. At the galaxy scale, the Lenses Structure and Dynamics (LSD) Survey is gathering kinematic data for distant (up to $z\sim1$) E/S0s that are gravitational lenses. A joint lensing and dynamical analysis constrains the fraction of dark matter within the Einstein radius, the mass-to-light ratio of the stellar component, and the total slope of the mass density profile. These properties and their evolution with redshift are briefly discussed in terms of the formation and evolution of E/S0 galaxies and measurement of the Hubble Constant from gravitational time delay systems. At the cluster scale -- after careful removal of the stellar component with a joint lensing and dynamical analysis -- systems with giant radial arcs can be used to measure precisely the inner slope of the dark matter halo. An HST search for radial arcs and the analysis of a first sample are briefly discussed in terms of the universal dark matter halos predicted by CDM simulations.

In this paper we present new and improved data on 38 cool white dwarfs identified by Oppenheimer et al. 2001 (OHDHS) as candidate dark halo objects. Using the high-res spectra obtained with LRIS, we measure radial velocities for 13 WDs that show an H alpha line. We show that the knowledge of RVs decreases the UV-plane velocities by only 6%. The radial velocity sample has a W-velocity dispersion of sig_W = 59 km/s--in between the values associated with the thick disk and the stellar halo. We also see indications for the presence of two populations by analyzing the velocities in the UV plane. In addition, we present CCD photometry for half of the sample, and with it recalibrate the photographic photometry of the remaining WDs. Using the new photometry in standard bands, and by applying the appropriate color-magnitude relations for H and He atmospheres, we obtain new distance estimates. New distances of the WDs that were not originally selected as halo candidates yield 13 new candidates. On average, new distances produce velocities in the UV plane that are larger by 10%, with already fast objects gaining more. Using the new data, while applying the same UV-velocity cut (94 km/s) as in OHDHS, we find a density of cool WDs of 1.7e-4 pc^-3, confirming the value of OHDHS. In addition, we derive the density as a function of the UV-velocity cutoff. The density (corrected for losses due to higher UV cuts) starts to flatten out at 150 km/s (0.4e-4 pc^-3), and is minimized (thus minimizing a possible non-halo contamination) at 190 km/s (0.3e-4 pc^-3). These densities are in a rough agreement with the estimates for the stellar halo WDs, corresponding to a factor of 1.9 and 1.4 higher values.

We present a refined gravitational lens model of the four-image lens system B1608+656 based on new and improved observational constraints: (i) the three independent time-delays and flux-ratios from VLA observations, (ii) the radio-image positions from VLBA observations, (iii) the shape of the deconvolved Einstein Ring from optical and infrared HST images, (iv) the extinction-corrected lens-galaxy centroids and structural parameters, and (v) a stellar velocity dispersion, sigma_ap=247+-35 km/s, of the primary lens galaxy (G1), obtained from an echelle spectrum taken with the Keck--II telescope. The lens mass model consists of two elliptical mass distributions with power-law density profiles and an external shear, totaling 22 free parameters, including the density slopes which are the key parameters to determine the value of H_0 from lens time delays. This has required the development of a new lens code that is highly optimized for speed. The minimum-chi^2 model reproduces all observations very well, including the stellar velocity dispersion and the shape of the Einstein Ring. A combined gravitational-lens and stellar dynamical analysis leads to a value of the Hubble Constant of H_0=75(+7/-6) km/s/Mpc (68 percent CL; Omega_m=0.3, Omega_Lambda=0.7. The non-linear error analysis includes correlations between all free parameters, in particular the density slopes of G1 and G2, yielding an accurate determination of the random error on H_0. The lens galaxy G1 is ~5 times more massive than the secondary lens galaxy (G2), and has a mass density slope of gamma_G1=2.03(+0.14/-0.14) +- 0.03 (68 percent CL) for rho~r^-gamma', very close to isothermal (gamma'=2). (Abridged)

The Lenses Structure & Dynamics (LSD) Survey aims at studying the internal structure of luminous and dark matter - as well as their evolution - of field early-type (E/SO) galaxies to z~1. In particular, E/S0 lens galaxies are studied by combining gravitational lensing, photometric and kinematic data obtained with ground-based (VLA/Keck/VLT) and space-based telescopes (HST). Here, we report on preliminary results from the LSD Survey, in particular on (i) the constraints set on the luminous and dark-matter distributions in the inner several R_eff of E/S0 galaxies, (ii) the evolution of their stellar component and (iii) the constraints set on the value of H_0 from time-delay systems by combining lensing and kinematic data to break degeneracies in gravitational-lens models.

We present flux-ratio curves of the fold and cusp (i.e. close multiple) images of six JVAS/CLASS gravitational lens systems. The data were obtained over a period of 8.5 months in 2001 with the Multi-Element Radio-Linked Interferometer Network (MERLIN) at 5-GHz with 50 mas resolution, as part of a MERLIN Key-Project. Even though the time delays between the fold and cusp images are small (<~1 day) compared to the time-scale of intrinsic source variability, all six lens systems show evidence that suggests the presence of extrinsic variability. In particular, the cusp images of B2045+265 -- regarded as the strongest case of the violation of the cusp relation (i.e. the sum of the magnifications of the three cusp images add to zero) -- show extrinsic variations in their flux-ratios up to ~40 percent peak-to-peak on time scales of several months. Its low Galactic latitude of b=-10 degree and a line-of-sight toward the Cygnus superbubble region suggest that Galactic scintillation is the most likely cause. The cusp images of B1422+231 at b=+69 degree do not show strong extrinsic variability. Galactic scintillation can therefore cause significant scatter in the cusp and fold relations of some radio lens systems (up to 10 percent rms), even though these relations remain violated when averaged over a <~1 year time baseline.

A deep spectrum taken with the Echelle Spectrograph and Imager (ESI) at the Keck II Telescope as part of the Lenses Structure and Dynamics (LSD) Survey reveals the redshifts of the extremely red source of the radio Einstein Ring in the gravitational lens system MG1549+305 ($z_{\rm s}=1.170\pm 0.001$) and an intermediate redshift lensed spiral galaxy ($z_{\rm G2}=0.604\pm 0.001$). The source redshift allows us to determine the mass of the SB0 lens galaxy enclosed by the Einstein Radius ($R_{\rm E}=1\farcs15\pm0\farcs05$) $M_{\rm E}$$\equiv$$M(<R_{\rm E}) = 8.4\pm0.7\times 10^{10} h_{65}^{-1}$ M$_\odot$. This corresponds to a Singular Isothermal Ellipsoid (SIE) velocity dispersion $\sigma_{\rm SIE}=214\pm5$ \kms, in good agreement with the measured stellar velocity dispersion $\sigma=227\pm18$ \kms (Leh\'ar et al. 1996). The mass-to-light ratio within the Einstein Radius ($\sim$1.4 effective radii) is $10\pm1 h_{65}$ \mlu. This is only marginally larger than typical stellar mass-to-light ratios of local early-type galaxies, indicating that dark matter is not likely to be dominant inside the Einstein Radius.

We report the final results of the search for gravitationally lensed flat-spectrum radio sources found in the combination of CLASS (Cosmic Lens All Sky Survey) and JVAS (Jodrell-Bank VLA Astrometric Survey). VLA observations of 16,503 sources have been made, resulting in the largest sample of arcsec-scale lens systems available. Contained within the 16,503 sources is a complete sample of 11,685 sources having two-point spectral indices between 1.4 and 5 GHz flatter than -0.5 and 5 GHz flux densities $\geq$30 mJy. A subset of 8,958 sources form a well-defined statistical sample suitable for analysis of the lens statistics. We describe the systematic process by which 149 candidate lensed sources were picked from the statistical sample on the basis of possessing multiple compact components in the 0.2 arcsec-resolution VLA maps. Candidates were followed up with 0.05 arcsec resolution MERLIN and 0.003 arcsec VLBA observations at 5 GHz and rejected as lens systems if they failed well-defined surface brightness and/or morphological tests. Maps for all the candidates can be found on the World Wide Web at http://www.jb.man.ac.uk/research/gravlens/index.html We summarize the properties of each of the 22 gravitational lens systems in JVAS/CLASS. Twelve are double-image systems, nine are four-image systems and one is a six-image system. Thirteen constitute a statistically well-defined sample giving a point-source lensing rate of 1:690$\pm$190. The interpretation of the results in terms of the properties of the lensing galaxy population and cosmological parameters will be published elsewhere. (Abridged)

The Cosmic Lens All-Sky Survey (CLASS) is an international collaborative program which has obtained high-resolution radio images of over 10000 flat-spectrum radio sources in order to create the largest and best studied statistical sample of radio-loud gravitationally lensed systems. With this survey, combined with detailed studies of the lenses found therein, constraints can be placed on the expansion rate, matter density, and dark energy (e.g. cosmological constant, quintessence) content of the Universe that are complementary to and independent of those obtained through other methods. CLASS is aimed at identifying lenses where multiple images are formed from compact flat-spectrum radio sources, which should be easily identifiable in the radio maps. Because CLASS is radio-based, dust obscuration in lensing galaxies is not a factor, and the relative insensitivity of the instrument to environmental conditions leads to nearly uniform sensitivity and resolution over the entire survey. In four observing seasons from 1994-1999, CLASS has observed 13783 radio sources with the VLA at 8.4 GHz at 0.2 arcsecond resolution. When combined with the JVAS survey, the CLASS sample contains over 16,000 images. A complete sample of 11685 flat-spectrum sources was observed, selected from GB6 catalogue at 4.85 GHz and the NVSS at 1.4 GHz. So far, CLASS has found 16 new gravitational lens systems, and the JVAS/CLASS survey contains a total of 22 lenses. (Abridged)

We report the discovery of a new gravitational lens system from the CLASS survey. Radio observations with the VLA, the WSRT and MERLIN show that the radio source B0850+054 is comprised of two compact components with identical spectra, a separation of 0.7 arcsec and a flux density ratio of 6:1. VLBA observations at 5 GHz reveal structures that are consistent with the gravitational lens hypothesis. The brighter of the two images is resolved into a linear string of at least six sub-components whilst the weaker image is radially stretched towards the lens galaxy. UKIRT K-band imaging detects an 18.7 mag extended object, but the resolution of the observations is not sufficient to resolve the lensed images and the lens galaxy. Mass modelling has not been possible with the present data and the acquisition of high-resolution optical data is a priority for this system.

A new two-image gravitational lens system has been discovered as a result of the Cosmic Lens All-Sky Survey (CLASS). Radio observations with the VLA, MERLIN and the VLBA at increasingly higher resolutions all show two components with a flux density ratio of ~7:1 and a separation of 1.34". Both components are compact and have the same spectral index. Followup observations made with the VLA at 8.4 GHz show evidence of a feature to the south-east of the brighter component and a corresponding extension of the weaker component to the north-west. Optical observations with the WHT show ~1.7" extended emission aligned in approximately the same direction as the separation between the radio components with an R-band magnitude of 21.8 +/- 0.4.

We have measured the kinematic profile of the early-type (E/S0) lens galaxy in the system 0047-281 (z=0.485) with the {\sl Echelle Spectrograph and Imager} (ESI) on the W.M. Keck--II Telescope, as part of the {\sl Lenses Structure and Dynamics (LSD) Survey}. The central velocity dispersion is \sigma=229\pm 15 \kms, and the dispersion profile is nearly flat to beyond one effective radius (R_e). Surface photometry of the lens galaxy is measured from {\sl Hubble Space Telescope} images. From the offset from the local Fundamental Plane (FP), we measure an evolution of the effective mass-to-light ratio of \Delta \log M/L_B=-0.37\pm0.06 between z=0 and z=0.485, consistent with the observed evolution of field E/S0 galaxies. (We assume h_{65}=1,\Omega_{\rm m}=0.3 and \Omega_\Lambda=0.7 throughout.) Gravitational lens models provide a mass of M_{\rmE}=(4.06\pm0.20)\times 10^{11} h_{65}^{-1}M_\odot inside the Einstein radius of R_{\rm E}=(8.70\pm0.07) h_{65}^{-1} kpc. This allows us to break the degeneracy between velocity anisotropy and density profile, typical of dynamical models for E/S0 galaxies. We find that constant M/L model, even with strongly tangential anisotropy of the stellar velocity ellipsoid, are excluded at >99.9%CL. The total mass distribution inside R_{\rm E} can be described by a single power-law density profile, \rho_t\propto r^{-\gamma'}, with an effective slope \gamma'=1.90^{+0.05}_{-0.23} (68%CL; \pm0.1 systematic error). Two-component models yield an upper limit (68% CL) of \gamma\le 1.55(1.12) on the power-law slope of the dark-matter density profile and a projected dark-matter mass fraction of 0.41(0.54)^{+0.15}_{-0.05}(^{+0.09}_{-0.06}) (68% CL) inside R_{\rm E}, for Osipkov--Merritt models with anisotropy radius r_i=\infty(R_e).

We combine lensing, stellar kinematic and mass-to-light ratio constraints to build a two-component (luminous plus dark) mass model of the early-type lens galaxy in PG1115+080. We find a total mass density profile steeper than r^-2, effectively rho r^-gamma' with gamma'=2.35+-0.1+-0.05 (random + systematic). The stellar mass fraction is f*=0.67(+0.20/-0.25)+-0.03 inside the Einstein Radius (R_E 1.2 effective radii). The dynamical mass model breaks the degeneracies in the mass profile of the lens galaxy and allows us to obtain a value of the Hubble Constant that is no longer dominated by systematic errors: H_0=59(+12/-7)+-3 km/s/Mpc (68 percent C.L.; Omega_m=0.3, Omega_L=0.7). The offset of PG1115+080 from the Fundamental Plane might indicate deviations from homology of the mass profile of some early-type galaxies.

We derive constraints on cosmological parameters and the properties of the lensing galaxies from gravitational lens statistics based on the final Cosmic Lens All Sky Survey (CLASS) data. For a flat universe with a classical cosmological constant, we find that the present matter fraction of the critical density is $\Omega_{\rm m}=0.31^{+0.27}_{-0.14}$ (68%) $^{+0.12}_{-0.10}$ (systematic). For a flat universe with a constant equation of state for dark energy $w = p_x({pressure})/\rho_x({energy density})$, we find $w < -0.55^{+0.18}_{-0.11}$ (68%).

The MACHO and OGLE collaborations have argued that the three longest duration bulge microlensing events are likely caused by nearby black holes, given the small velocities measured with microlensing parallax and nondetection of the lenses. However, these events may be due to lensing by more numerous lower mass stars at greater distances. We find a posteriori probabilities of 76%, 16%, and 4% that the three longest events are black holes, assuming a Salpeter initial mass function (IMF) and a 40 M cutoff for neutron star progenitors; the numbers depend strongly on the assumed mass function but favor a black hole for the longest event for most standard IMFs. The longest events (>600 days) have an a priori 26% probability of being black holes for a standard mass function. We propose a new technique for measuring the lens mass function using the mass distribution of long events measured with the Advanced Camera for Surveys on the Hubble Space Telescope, the Very Large Telescope Interferometer, the Space Interferometry Mission, or the Global Astrometric Interferometer for Astrophysics.

The gravitational lens CLASS B1608+656 is the only four-image lens system for which all three independent time delays have been measured. This makes the system an excellent candidate for a high-quality determination of H_0 at cosmological distances. However, the original measurements of the time delays had large (12-20%) uncertainties, due to the low level of variability of the background source during the monitoring campaign. In this paper, we present results from two additional VLA monitoring campaigns. In contrast to the ~5% variations seen during the first season of monitoring, the source flux density changed by 25-30% in each of the subsequent two seasons. We analyzed the combined data set from all three seasons of monitoring to improve significantly the precision of the time delay measurements; the delays are consistent with those found in the original measurements, but the uncertainties have decreased by factors of two to three. We combined the delays with revised isothermal mass models to derive a measurement of H_0. Depending on the positions of the galaxy centroids, which vary by up to 0.1 arcsec in HST images obtained with different filters, we obtain H_0 = 61-65 km/s/Mpc, for (Omega_M,Omega_L) = (0.3,0.7). The value of H_0 decreases by 6% if (Omega_M,Omega_L) = (1.0,0.0). The formal uncertainties on H_0 due to the time delay measurements are +/- 1 (+/- 2) km/s/Mpc for the 1-sigma (2-sigma) confidence limits. Thus, the systematic uncertainties due to the lens model, which are on the order of +/- 15 km/s/Mpc, now dominate the error budget for this system. In order to improve the measurement of H_0 with this lens, new models that incorporate the constraints provided by stellar dynamics and the optical/infrared Einstein ring seen in HST images must be developed.

The Lenses Structure and Dynamics (LSD) Survey aims at measuring the luminous and dark matter distribution of high redshift early-type galaxies by combining gravitational lens analysis with newly determined velocity dispersion profiles. Here we describe the first results from the LSD survey, a measurement of the internal structure of the lens galaxy in MG2016+112 (z=1.004). The relevance of this measurement to the cosmological model is briefly discussed, in particular in the context of E/S0 galaxy formation and evolution, the existence and shape of universal dark matter profiles, and the measurement of the Hubble constant from gravitational time delay.

[Abridged] We combine our measurements of the velocity dispersion and the surface brightness profile of the lens galaxy D in the system MG2016+112 (z=1.004) with constraints from gravitational lensing to study its internal mass distribution. We find that: (i) dark matter accounts for >50% of the total mass within the Einstein radius (99% CL), excluding at the 8-sigma level that mass follows light inside the Einstein radius with a constant mass-to-light ratio (M/L). (ii) the total mass distribution inside the Einstein radius is well-described by a density profile ~r^-gamma' with an effective slope gamma'=2.0+-0.1+-0.1, including random and systematic uncertainties. (iii) The offset of galaxy D from the local Fundamental Plane independently constrains the stellar M/L, and matches the range derived from our models, leading to a more stringent lower limit of >60% on the fraction of dark matter within the Einstein radius (99%CL). Under the assumption of adiabatic contraction, the inner slope of the dark matter halo before the baryons collapsed is gamma_i<1.4 (68 CL), marginally consistent with the highest-resolution cold dark matter simulations that indicate gamma_i~1.5. This might indicate that either adiabatic contraction is a poor description of E/S0 formation or that additional processes play a role as well. Indeed, the apparently isothermal density distribution inside the Einstein radius, is not a natural outcome of adiabatic contraction models, where it appears to be a mere coincidence. By contrast, we argue that isothermality might be the result of a stronger coupling between luminous and dark-matter, possibly the result of (incomplete) violent relaxation processes. Hence, we conclude that galaxy D appears already relaxed 8 Gyr ago.

We present a direct measurement of the stellar velocity dispersion of the early-type lens galaxy D in the system MG2016+112 (z=1.004), determined from a spectrum obtained with the Echelle Spectrograph and Imager (ESI) on the W.M. Keck-II Telescope, as part of the Lenses Structure and Dynamics (LSD) Survey. We find a velocity dispersion of sigma_ap=304+-27 km/s inside an effective circular aperture with a radius of 0.65'', corresponding to a central velocity dispersion of sigma=328+-32 km/s. From a Hubble Space Telescope F160W-band image, we measure the effective radius and effective surface brightness in order to determine the offset of the lens galaxy with respect to the local Fundamental Plane. The offset corresponds to an evolution of the rest-frame effective mass-to-light ratio of Delta log(M/L_B)=-0.62+-0.08 from z=0 to z=1.004. By interpreting colors and offset of the FP with two independent stellar population synthesis models, we obtain a single-burst equivalent age of 2.8+-0.8 Gyr (i.e. z_f>1.9) and a supersolar metallicity of log[Z/Z_sun]=0.25+-0.25. The lens galaxy is therefore a massive elliptical dominated by an old and metal rich stellar population at z>1. The excellent agreement of the stellar velocity dispersion with that predicted from recent lens models confirms that the angular separation of the multiple images of the background QSO is predominantly due to the lens galaxy, and not to a massive ``dark cluster'', in agreement with recent weak lensing and X-ray observations. However, the significant overdensity of galaxies in the field might indicate that this system is a proto-cluster, in formation around galaxy D, responsible for the ~10% external shear inferred from the strong lens models.

A single-screen model of the gravitational lens system 2016+112 is proposed, that explains recent Hubble Space Telescope} (HST) infrared (NICMOS-F160W) observations and new high-resolution European VLBI Network (EVN) 5-GHz radio observations, presented in this paper. In particular, we find that a massive `dark' structure at the lens position, previously suggested by X-ray, optical and spectroscopic observations of the field around 2016+112, is not necessarily required to accommodate the strong lensing constraints. A massive structure to the north-west of the lens system, suggested from a weak-lensing analysis of the field, is included in the model. The lensed source is an X-ray bright active galaxy at z=3.273 with a central bright optical continuum core and strong narrow emission lines, suggestive of a type-II quasar. The EVN 5-GHz radio maps show a radio-jet structure with at least two compact subcomponents. We propose that the diamond caustic crosses the counter-jet of the radio source, so that part of the counter-jet, host galaxy and narrow-line emission regions are quadruply imaged. The remainder of the radio source, including the core, is doubly imaged. Our lens model predicts a very high magnification (mu~300) at the bightness peaks of the inner two radio components of complex C. If the jet exhibits relativistic velocities on micro-arsecond scales, it might result in apparent hyperluminal motion. However, the lack of strong radio variability and the peaked radio spectrum imply that these motions need not be present in the source. Our model furthermore implies that the optical spectrum of C' can only show features of the AGN and its host galaxy.

We analyse the kinematics of the entire spectroscopic sample of 99 recently discovered high proper-motion white dwarfs by Oppenheimer et al. using a maximum-likelihood analysis, and discuss the claim that the high-velocity white dwarfs are members of a halo population with a local density at least ten times greater than traditionally assumed. We argue that the observations, as reported, are consistent with the presence of an almost undetected thin disc plus a thick disc, with densities as conventionally assumed. In addition, there is a kinematically distinct, flattened, halo population at the more than 99% confidence level. Surprisingly, the thick disc and halo populations are indistinguishable in terms of luminosity, color and apparent age (1-10 Gyr). Adopting a bimodal, Schwarzschild model for the local velocity ellipsoid, with the ratios sigma_U:sigma_V:sigma_W=1:2/3:1/2, we infer radial velocity dispersions of sigma_U=62(+8/-10) km/s and 150(+80/-40) km/s (90% C.L.) for the local thick disc and halo populations, respectively. The thick disc result agrees with the empirical relation between asymmetric drift and radial velocity dispersion, inferred from local stellar populations. The local thick-disc plus halo density of white dwarfs is n^{td+h}=(1.9+-0.5)x10^-3 pc^-3 (90% C.L.), of which n^{h}=1.1(+2.1/-0.7)x10^-4 pc^-3 (90% C.L.) belongs to the halo, a density about five times higher than previously thought. (Abridged)

We examine the claim by Oppenheimer et al. (2001) that the local halo density of white dwarfs is an order of magnitude higher than previously thought. As it stands, the observational data support the presence of a kinematically distinct population of halo white dwarfs at the >99% confidence level. A maximum-likelihood analysis gives a radial velocity dispersion of sigma^h_U=150(+80/-40) km/s and an asymmetric drift of v_a^h=176(+102/-80) km/s, for a Schwarzschild velocity distribution function with sigma_U:sigma_V:sigma_W=1:2/3:1/2. Halo white dwarfs have a local number density of 1.1(+2.1/-0.7)x10^-4 pc-3, which amounts to 0.8(+1.6/-0.5) per cent of the nominal local dark-matter halo density and is 5.0(+9.5/-3.2) times higher and thus only marginally in agreement with previous estimates (all errors indicate the 90% C.L.). We discuss several direct consequences of this white-dwarf population (e.g. microlensing) and postulate a potential mechanism to eject young white dwarfs from the disc to the halo, through the orbital instabilities in triple or multiple stellar systems.

Recent developments in the spectropolarimetric study of compact objects, specifically black holes (stellar and massive) and neutron stars are reviewed. The lectures are organized around five topics: disks, jets, outflows, neutron stars and black holes. They emphasize physical mechanisms and are intended to bridge the gap between the fundamentals of polarimetry and the phenomenology of observed cosmic sources of polarized radiation, as covered by the other lecturers. There has been considerable recent progress in spectropolarimetry from radio through optical frequencies and this is producing some unique diagnostics of the physical conditions around compact objects. It is argued that there is a great need to develop a correspondingly sensitive polarimetric capability at ultraviolet through gamma-ray energies.

We report the discovery of a new two-image gravitational lens system from the Cosmic Lens All-Sky Survey, CLASS B2319+051. Radio imaging with the Very Large Array (VLA) and Multi-Element Radio-Linked Interferometer Network (MERLIN) shows two compact components with a flux density ratio of 5:1, separated by 1.36 arcsec. Observations with the Very Long Baseline Array (VLBA) resolve each of the radio components into a pair of parity-reversed subcomponents. Hubble Space Telescope (HST) observations with the Near-Infrared Camera and Multi-Object Spectrometer (NICMOS) show a bright elliptical galaxy (G1) coincident with the radio position, and a second irregular galaxy (G2) 3.4 arcsec to the northwest. Previous spectroscopic studies have indicated that these galaxies are at different redshifts: z(G1) = 0.624, z(G2) = 0.588. Infrared counterparts to the lensed radio components are not detected in the NICMOS image, and the source redshift has not yet been determined. Preliminary mass modeling based on the VLBA subcomponent data indicates that the lensing potential includes a strong external shear contribution. A VLA monitoring program is currently being undertaken to measure the differential time delay.

The declining lightcurve of the optical afterglow of gamma-ray burst GRB000301C showed rapid variability with one particularly bright feature at about t-t0=3.8 days. This event was interpreted as gravitational microlensing by Garnavich, Loeb & Stanek (2000) and subsequently used to derive constraints on the structure of the GRB optical afterglow. In this paper, we use these structural parameters to calculate the probability of such a microlensing event in a realistic scenario, where all compact objects in the universe are associated with observable galaxies. For GRB000301C at a redshift of z=2.04, the `a posteriori' probability for a microlensing event with an amplitude of dm>=0.95 mag (as observed) is 0.7% (2.7%) for the most plausible scenario of a flat Lambda-dominated FRW universe with Omega_m=0.3 and a fraction f_*=0.2 (1.0) of dark-matter in the form of compact objects. If we lower the magnification threshold to dm>=0.10 mag, the probabilities for microlensing events of GRB afterglows increase to 17% (57%). We emphasise that this low probability for a microlensing signature of almost a magnitude does not exclude that the observed event in the afterglow lightcurve of GRB000301C was caused by microlensing, especially in light of the fact that a galaxy was found within 2 arcsec from the GRB. In that case, however, a more robust upper limit on the `a posteriori' probability of ~5% is found. It does show, however, that it will not be easy to create a large sample of strong GRB afterglow microlensing events for statistical studies of their physical conditions on micro-arcsec scales.

One of the main objectives of the Cosmic Lens All-Sky Survey (CLASS) collaboration has been to find gravitational lens (GL) systems at radio wavelengths that are suitable for the determination of time delays between image pairs. The survey is now near completion and at least 18 GL systems have been found. Here, I will discuss our efforts to measure time delays from several of these systems with the ultimate aim of constraining the Hubble Constant (H0). Thus far three CLASS GL systems (i.e. B0218+357, B1600+434 and B1608+656) have yielded measurements of time delays, from which values of H0=60-70 km/s/Mpc have been estimated. Although most GL systems give similar values of H0, statistical and systematic uncertainties are still considerable. To reduce these uncertainties, I will shortly mention two monitoring programs that we are undertaking to (re)measure time delays in 15 CLASS GL systems and address several important issues for the near future.

Strongly correlated non-intrinsic variability between 5 and 8.5 GHz has been observed in one of the lensed images of the gravitational lens B1600+434. These non-intrinsic (i.e. `external') variations are interpreted as `radio-microlensing' of relativistic micro-arsec-scale jet components in the source at a redshift of z=1.59 by massive compact objects in the halo of the edge-on disk lens galaxy at z=0.41. We shortly summarize these observations and discuss several new observational and theoretical programs to investigate this new phenomenon in more detail.

High resolution MERLIN observations of a newly-discovered four-image gravitational lens system, B0128+437, are presented. The system was found after a careful re-analysis of the entire CLASS dataset. The MERLIN observations resolve four components in a characteristic quadruple-image configuration; the maximum image separation is 542 mas and the total flux density is 48 mJy at 5 GHz. A best-fit lens model with a singular isothermal ellipsoid results in large errors in the image positions. A significantly improved fit is obtained after the addition of a shear component, suggesting that the lensing system is more complex and may consist of multiple deflectors. The integrated radio spectrum of the background source indicates that it is a GigaHertz-Peaked Spectrum (GPS) source. It may therefore be possible to resolve structure within the radio images with deep VLBI observations and thus better constrain the lensing mass distribution.

We present the discovery of CLASS B0739+366, a new gravitational lens system from the Cosmic Lens All-Sky Survey. Radio imaging of the source with the Very Large Array (VLA) shows two compact components separated by $0\farcs54$, with a flux density ratio of $\sim$ 6:1. High-resolution follow-up observations using the Very Long Baseline Array (VLBA) at 1.7 GHz detect weak, parity-reversed jet emission from each of the radio components. Hubble Space Telescope NICMOS F160W observations detect infrared counterparts to the lensed images, as well as an extended object between them which we identify as the lensing galaxy. Redshifts for the galaxy and lensed source have not yet been obtained. For typical lens and source redshifts of $z=0.5$ and $z=1.5$, respectively, preliminary mass modeling predicts a time delay of $\sim7h^{-1}$ days in a flat $\Omega_{M}=1.0$ universe. The small predicted time delay and weak radio components will make CLASS B0739+366 a challenging target for Hubble constant determination.

We present radio, optical, near-infrared and spectroscopic observations of the source B0827+525. We consider this source as the best candidate from the Cosmic Lens All-Sky Survey (CLASS) for a `dark lens' system or binary radio-loud quasar. The system consists of two radio components with somewhat different spectral indices, separated by 2.815 arcsec. VLBA observations show that each component has substructure on a scale of a few mas. A deep K-band exposure with the W.M.Keck-II Telescope reveals emission near both radio components. The K-band emission of the weaker radio component appears extended, whereas the emission from the brighter radio component is consistent with a point source. Hubble Space Telescope F160W-band observations with the NICMOS instrument confirms this. A redshift of 2.064 is found for the brighter component, using the LRIS instrument on the W.M.Keck-II Telescope. The probability that B0827+525 consists of two unrelated compact flat-spectrum radio sources is ~3%, although the presence of similar substructure in both component might reduce this. We discuss two scenarios to explain this system: (i) CLASS B0827+525 is a `dark lens' system or (ii) B0827+525 is a binary radio-loud quasar. B0827+525 has met all criteria that thus far have in 100% of the cases confirmed a source as an indisputable gravitational lens system. Despite this, no lens galaxy has been detected with m_F160W<=23 mag. Hence, we might have found the first binary radio-loud quasar. At this moment, however, we feel that the `dark lens' hypothesis cannot yet be fully excluded.

Recent observations indicate that many if not all galaxies host massive central black holes. In this paper we explore the influence of black holes on the lensing properties. We model the lens as an isothermal ellipsoid with a finite core radius plus a central black hole. We show that the presence of the black hole substantially changes the critical curves and caustics. If the black hole mass is above a critical value, then it will completely suppress the central images for all source positions. Realistic central black holes likely have masses below this critical value. Even in such sub-critical cases, the black hole can suppress the central image when the source is inside a zone of influence, which depends on the core radius and black hole mass. In the sub-critical cases, an additional image may be created by the black hole in some regions, which for some radio lenses may be detectable with high-resolution and large dynamic-range VLBI maps. The presence of central black holes should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.

The complex ten-component gravitational lens system B1933+503 has been monitored with the VLA during the period February to June 1998 with a view to measuring the time delay between the four compact components and hence to determine the Hubble parameter. Here we present the results of an `A' configuration 8.4-GHz monitoring campaign which consists of 37 epochs with an average spacing of 2.8 days. The data have yielded light curves for the four flat-spectrum radio components (components 1, 3, 4 and 6). We observe only small flux density changes in the four flat-spectrum components which we do not believe are predominantly intrinsic to the source. Therefore the variations do not allow us to determine the independent time delays in this system. However, the data do allow us to accurately determine the flux density ratios between the four flat-spectrum components. These will prove important as modelling constraints and could prove crucial in future monitoring observations should these data show only a monotonic increase or decrease in the flux densities of the flat-spectrum components.

First, we review the current status of the detection of strong `external' variability in the CLASS gravitational B1600+434, focusing on the 1998 VLA 8.5-GHz and 1998/9 WSRT multi-frequency observations. We show that this data can best be explained in terms of radio-microlensing. We then proceed to show some preliminary results from our new multi-frequency VLA monitoring program, in particular the detection of a strong feature (~30%) in the light curve of the lensed image which passes predominantly through the dark-matter halo of the lens galaxy. We tentatively interpret this event, which lasted for several weeks, as a radio-microlensing caustic crossing, i.e. the superluminal motion of a micro-arcsec-scale jet-component in the lensed source over a single caustic in the magnification pattern, that has been created by massive compact objects along the line-of-sight to the lensed image.

We present the first unambiguous case of external variability of a radio gravitational lens, CLASS B1600+434. The VLA 8.5-GHz difference light curve of the lensed images, taking the proper time-delay into account, shows the presence of external variability with 14.6-sigma confidence. We investigate two plausible causes of this external variability: scattering by the ionized component of the Galactic interstellar medium and microlensing by massive compact objects in the bulge/disk and halo of the lens galaxy. Based on the tight relation between the modulation-index and variability time scale and the quantitative difference between the light curves of both lensed images, we conclude that the observed short-term variability characteristics of the lensed images are incompatible with scintillation in our Galaxy. This conclusion is strongly supported by multi-frequency WSRT observations at 1.4 and 5 GHz, which are in strong disagreement with predictions based on the scintillation hypothesis. ... On the other hand, a single superluminal jet-component, having an apparent velocity 9<=(v_app/c)<=26, a radius of 2-5 micro-arcsec and containing 5-11% of the observed 8.5-GHz source flux density, can reproduce the observed modulation-indices and variability time scale at 8.5 GHz, when it is microlensed by compact objects in the lens galaxy. It also reproduces the frequency-dependence of the modulation-indices, determined from the independent WSRT 1.4 and 5-GHz observations. ... The only conclusion fully consistent with the data gathered thus far is that we have indeed detected radio microlensing. The far reaching consequence of this statement is that a significant fraction of the mass in the dark-matter halo at ~6 kpc (h=0.65) above the lens-galaxy disk in B1600+434 consists of massive compact objects. [abridged]