We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray-tracing through light-cones, and use the aperture mass statistic for the analysis. The peaks detected on the maps reflect the cumulative signal from massive bound objects and general large-scale structure. We present the first study of weak lensing peaks in simulations that include both baryonic physics and massive neutrinos (summed neutrino mass $M_{\nu} =$ 0.06, 0.12, 0.24, and 0.48 eV assuming normal hierarchy), so that the uncertainty due to physics beyond the gravity of dark matter can be factored into constraints on cosmological models. Assuming a fiducial model of baryonic physics, we also investigate the correlation between peaks and massive haloes, over a range of summed neutrino mass values. As higher neutrino mass tends to suppress the formation of massive structures in the Universe, the halo mass function and lensing peak counts are therefore modified as a function of $M_{\nu}$. Over most of the S/N range, the impact of fiducial baryonic physics is greater (less) than neutrinos for 0.06 and 0.12 (0.24 and 0.48) eV models. Both baryonic physics and massive neutrinos should be accounted for when deriving cosmological parameters from weak lensing observations.

Weak gravitational lensing depends on the integrated mass along the line of sight. Baryons contribute to the mass distribution of galaxy clusters and the resulting mass estimates from lensing analysis. We use the cosmo-OWLS suite of hydrodynamic simulations to investigate the impact of baryonic processes on the bias and scatter of weak lensing mass estimates of clusters. These estimates are obtained by fitting NFW profiles to mock data using MCMC techniques. In particular, we examine the difference in estimates between dark matter-only runs and those including various prescriptions for baryonic physics. We find no significant difference in the mass bias when baryonic physics is included, though the overall mass estimates are suppressed when feedback from AGN is included. For lowest-mass systems for which a reliable mass can be obtained ($M_{200} \approx 2 \times 10^{14}$ $M_{\odot}$), we find a bias of $\approx -10$ per cent. The magnitude of the bias tends to decrease for higher mass clusters, consistent with no bias for the most massive clusters which have masses comparable to those found in the CLASH and HFF samples. For the lowest mass clusters, the mass bias is particularly sensitive to the fit radii and the limits placed on the concentration prior, rendering reliable mass estimates difficult. The scatter in mass estimates between the dark matter-only and the various baryonic runs is less than between different projections of individual clusters, highlighting the importance of triaxiality.

For more than two decades, the Navarro, Frenk, and White (NFW) model has stood the test of time; it has been used to describe the distribution of mass in galaxy clusters out to their outskirts. Stacked weak lensing measurements of clusters are now revealing the distribution of mass out to and beyond their virial radii, where the NFW model is no longer applicable. In this study we assess how well the parameterised Diemer & Kravstov (DK) density profile describes the characteristic mass distribution of galaxy clusters extracted from cosmological simulations. This is determined from stacked synthetic lensing measurements of the 50 most massive clusters extracted from the Cosmo-OWLS simulations, using the Dark Matter Only run and also the run that most closely matches observations. The characteristics of the data reflect the Weighing the Giants survey and data from the future Large Synoptic Survey Telescope (LSST). In comparison with the NFW model, the DK model favored by the stacked data, in particular for the future LSST data, where the number density of background galaxies is higher. The DK profile depends on the accretion history of clusters which is specified in the current study. Eventually however subsamples of galaxy clusters with qualities indicative of disparate accretion histories could be studied.

Abell 2146 ($z=0.232$) is a massive galaxy cluster currently undergoing a spectacular merger in the plane of the sky with a bullet-like morphology. It was the first system in which both the bow and upstream shock fronts were detected at X-ray wavelengths (Mach$\sim2$), yet deep GMRT 325 MHz observations failed to detect extended radio emission associated with the cluster as is typically seen in such systems. We present new, multi-configuration $1-2$ GHz Karl G. Jansky Very Large Array (VLA) observations of Abell 2146 totalling 16 hours of observations. These data reveal for the first time the presence of an extended ($\approx850$ kpc), faint radio structure associated with Abell 2146. The structure appears to harbour multiple components, one associated with the upstream shock which we classify as a radio relic and one associated with the subcluster core which is consisted as being a radio halo bounded by the bow shock. The newly detected structures have some of the lowest radio powers detected thus far in any cluster ($P_{\rm 1.4 GHz, halo}=2.4\pm0.2\times10^{23}$ W Hz$^{-1}$ and $P_{\rm 1.4 GHz, relic}=2.2\pm0.2\times10^{23}$ W Hz$^{-1}$). The flux measurement of the halo, as well as its morphology, also suggest that the halo was recently created ($\approx0.3$ Gyr after core passage), consistent with the dynamical state of the cluster. These observations demonstrate the capacity of the upgraded VLA to detect extremely faint and extended radio structures. Based on these observations, we predict that many more radio relics and radio halos in merging clusters should be detected by future radio facilities such as the Square Kilometre Array (SKA).

The recent discovery of gravitational waves (GWs) from stellar-mass binary black holes (BBHs) provided direct evidence of the existence of these systems. BBH lenses would have gravitational microlensing signatures that are distinct from single-lens signals. We apply Bayesian statistics to examine the distinguishability of BBH microlensing events from single-lens events under ideal observing conditions, using the photometric capabilities of the Korean Microlensing Telescope Network. Given one year of observations, a source star at the Galactic Centre, a GW150914-like BBH lens (total mass 65M$_\odot$, mass ratio 0.8) at half that distance, and an impact parameter of 0.4 Einstein radii, we find that binarity is detectable for BBHs with separations down to 0.0250 Einstein radii, which is nearly 3.5 times greater than the maximum separation for which such BBHs would merge within the age of the Universe. Microlensing searches are thus sensitive to more widely separated BBHs than GW searches, perhaps allowing the discovery of BBH populations produced in different channels of binary formation.

Abell 2146 ($z$ = 0.232) consists of two galaxy clusters undergoing a major merger. The system was discovered in previous work, where two large shock fronts were detected using the $\textit{Chandra X-ray Observatory}$, consistent with a merger close to the plane of the sky, caught soon after first core passage. A weak gravitational lensing analysis of the total gravitating mass in the system, using the distorted shapes of distant galaxies seen with ACS-WFC on $\textit{Hubble Space Telescope}$, is presented. The highest peak in the reconstruction of the projected mass is centred on the Brightest Cluster Galaxy (BCG) in Abell 2146-A. The mass associated with Abell 2146-B is more extended. Bootstrapped noise mass reconstructions show the mass peak in Abell 2146-A to be consistently centred on the BCG. Previous work showed that BCG-A appears to lag behind an X-ray cool core; although the peak of the mass reconstruction is centred on the BCG, it is also consistent with the X-ray peak given the resolution of the weak lensing mass map. The best-fit mass model with two components centred on the BCGs yields $M_{200}$ = 1.1$^{+0.3}_{-0.4}$$\times$10$^{15}$M$_{\odot}$ and 3$^{+1}_{-2}$$\times$10$^{14}$M$_{\odot}$ for Abell 2146-A and Abell 2146-B respectively, assuming a mass concentration parameter of $c=3.5$ for each cluster. From the weak lensing analysis, Abell 2146-A is the primary halo component, and the origin of the apparent discrepancy with the X-ray analysis where Abell 2146-B is the primary halo is being assessed using simulations of the merger.

Abell 2146 consists of two galaxy clusters that have recently collided close to the plane of the sky, and it is unique in showing two large shocks on $\textit{Chandra X-ray Observatory}$ images. With an early stage merger, shortly after first core passage, one would expect the cluster galaxies and the dark matter to be leading the X-ray emitting plasma. In this regard, the cluster Abell 2146-A is very unusual in that the X-ray cool core appears to lead, rather than lag, the Brightest Cluster Galaxy (BCG) in their trajectories. Here we present a strong lensing analysis of multiple image systems identified on $\textit{Hubble Space Telescope}$ images. In particular, we focus on the distribution of mass in Abell 2146-A in order to determine the centroid of the dark matter halo. We use object colours and morphologies to identify multiple image systems; very conservatively, four of these systems are used as constraints on a lens mass model. We find that the centroid of the dark matter halo, constrained using the strongly lensed features, is coincident with the BCG, with an offset of $\approx$ 2 kpc between the centres of the dark matter halo and the BCG. Thus from the strong lensing model, the X-ray cool core also leads the centroid of the dark matter in Abell 2146-A, with an offset of $\approx$ 30 kpc.

We present a dynamical analysis of the merging galaxy cluster system Abell 2146 using spectroscopy obtained with the Gemini Multi-Object Spectrograph on the Gemini North telescope. As revealed by the Chandra X-ray Observatory, the system is undergoing a major merger and has a gas structure indicative of a recent first core passage. The system presents two large shock fronts, making it unique amongst these rare systems. The hot gas structure indicates that the merger axis must be close to the plane of the sky and that the two merging clusters are relatively close in mass, from the observation of two shock fronts. Using 63 spectroscopically determined cluster members, we apply various statistical tests to establish the presence of two distinct massive structures. With the caveat that the system has recently undergone a major merger, the virial mass estimate is M_vir = 8.5 +4.3 -4.7 x 10 ^14 M_sol for the whole system, consistent with the mass determination in a previous study using the Sunyaev-Zeldovich signal. The newly calculated redshift for the system is z = 0.2323. A two-body dynamical model gives an angle of 13-19 degrees between the merger axis and the plane of the sky, and a timescale after first core passage of 0.24-0.28 Gyr.

We determine the concentration-mass relation of 19 X-ray selected galaxy clusters from the CLASH survey in theories of gravity that directly modify the lensing potential. We model the clusters as NFW haloes and fit their lensing signal, in the Cubic Galileon and Nonlocal gravity models, to the lensing convergence profiles of the clusters. We discuss a number of important issues that need to be taken into account, associated with the use of nonparametric and parametric lensing methods, as well as assumptions about the background cosmology. Our results show that the concentration and mass estimates in the modified gravity models are, within the errorbars, the same as in $\Lambda$CDM. This result demonstrates that, for the Nonlocal model, the modifications to gravity are too weak at the cluster redshifts, and for the Galileon model, the screening mechanism is very efficient inside the cluster radius. However, at distances $\sim \left[2-20\right] {\rm Mpc}/h$ from the cluster center, we find that the surrounding force profiles are enhanced by $\sim20-40\%$ in the Cubic Galileon model. This has an impact on dynamical mass estimates, which means that tests of gravity based on comparisons between lensing and dynamical masses can also be applied to the Cubic Galileon model.

We present X-shooter observations of CASSOWARY 20 (CSWA 20), a star-forming (SFR ~6 Msol/yr) galaxy at z=1.433, magnified by a factor of 11.5 by the gravitational lensing produced by a massive foreground galaxy at z=0.741. We analysed the integrated physical properties of the HII regions of CSWA 20 using temperature- and density-sensitive emission lines. We find the abundance of oxygen to be ~1/7 of solar, while carbon is ~50 times less abundant than in the Sun. The unusually low C/O ratio may be an indication of a particularly rapid timescale of chemical enrichment. The wide wavelength coverage of X-shooter gives us access to five different methods for determining the metallicity of CSWA 20, three based on emission lines from HII regions and two on absorption features formed in the atmospheres of massive stars. All five estimates are in agreement, within the factor of ~2 uncertainty of each method. The interstellar medium of CSWA 20 only partially covers the star-forming region as viewed from our direction; in particular, absorption lines from neutrals and first ions are exceptionally weak. We find evidence for large-scale outflows of the interstellar medium (ISM) with speeds of up 750 km/s, similar to the values measured in other high-z galaxies sustaining much higher rates of star formation.

(Abridged) We quantify the bias and scatter in galaxy cluster masses and concentrations derived from an idealised mock weak gravitational lensing (WL) survey, and their effect on the cluster mass-concentration relation. For this, we simulate WL distortions on a population of background galaxies due to a large (~3000) sample of galaxy cluster haloes extracted from the Millennium Simulation at z~0.2. This study takes into account the influence of shape noise, cluster substructure and asphericity as well as correlated large-scale structure, but not uncorrelated large-scale structure along the line of sight and observational effects. We find a small, but non-negligble, negative median bias in both mass and concentration at a level of ~5%, the exact value depending both on cluster mass and radial survey range. Both the mass and concentration derived from WL show considerable scatter about their true values. This scatter has, even for the highest mass clusters of M200 > 10^14.8 M_sun, a level of ~30% and ~20% for concentration and mass respectively and increases strongly with decreasing cluster mass. For a typical survey analysing 30 galaxies per arcmin^2 over a radial range from 30" to 15' from the cluster centre, the derived M200-c relation has a slope and normalisation too low compared to the underlying true (3D) relation by ~40% and ~15% respectively. The scatter and bias in mass are shown to reflect a departure at large radii of the true WL shear/matter distribution of the simulated clusters from the NFW profile adopted in modelling the mock observations. Orientation of the triaxial cluster haloes dominates the concentration scatter (except at low masses, where galaxy shape noise becomes dominant), while the bias in c is mostly due to substructure within the virial radius.

We present a new 400 ks Chandra X-ray observation of the merging galaxy cluster Abell 2146. This deep observation reveals detailed structure associated with the major merger event including the Mach M=2.3+/-0.2 bow shock ahead of the dense, ram pressure stripped subcluster core and the first known example of an upstream shock in the ICM (M=1.6+/-0.1). By measuring the electron temperature profile behind each shock front, we determine the timescale for the electron population to thermally equilibrate with the shock-heated ions. We find that the temperature profile behind the bow shock is consistent with the timescale for Coulomb collisional equilibration and the postshock temperature is lower than expected for instant shock-heating of the electrons. Although like the Bullet cluster the electron temperatures behind the upstream shock front are hotter than expected, favouring the instant heating model, the uncertainty on the temperature values is greater here and there is significant substructure complicating the interpretation. We also measured the width of each shock front and the contact discontinuity on the leading edge of the subcluster core to investigate the suppression of transport processes in the ICM. The upstream shock is ~440 kpc in length but appears remarkably narrow over this distance with a best-fit width of only 6^{+5}_{-3} kpc compared with the mean free path of 23+/-5 kpc. The leading edge of the subcluster core is also narrow with an upper limit on the width of only 2 kpc separating the cool, multiphase gas at 0.5-2 keV from the shock-heated surrounding ICM at ~6 keV. The strong suppression of diffusion and conduction across this edge suggests a magnetic draping layer may have formed around the subcluster core.[abridged]

We propose convergent plasma lenses, possibly from current sheets, as a generic solution to strong interstellar scattering. These lenses resolve the overpressure problem by geometric alignment as noted by Goldreich and Shridhar (2006). They further quantitatively explain properties of extreme scattering events, and pulsar parabolic arcs. This model makes quantitative predictions testable by VLBI on scattering events. It differs conceptually from previous models by occurring through rare, localized underdense sheets. Such sheets are thermally and kinematically stable, and could be consequences of reconnection. The apparent diffractive effects are a result of coherent interference of refractive images. We propose that these lenses can be used for precision distance determination to pulsars, enabling accurate gravity source localization.

Using WHT OASIS integral field unit observations, we report the discovery of a thin plume of ionised gas extending from the brightest cluster galaxy in Abell 2146 to the sub-cluster X-ray cool core which is offset from the BCG by ~37 kpc. The plume is greater than 15 kpc long and less than 3 kpc wide. This plume is unique in that the cluster it is situated in is currently undergoing a major galaxy cluster merger. The brightest cluster galaxy is unusually located behind the X-ray shock front and in the wake of the ram pressure stripped X-ray cool core and evidence for recent disruption to the BCG is observed. We examine the gas and stellar morphology, the gas kinematics of the BCG and their relation to the X-ray gas. We propose that a causal link between the ionised gas plume and the offset X-ray cool core provides the simplest explanation for the formation of the plume. An interaction or merger between the BCG and another cluster galaxy is probably the cause of the offset.

A generic property of the cuspy simulated virialized halos in cold dark matter cosmogenies is that their concentration is inversely correlated with their mass. This behavior has also been confirmed in observations, although differences in the exact form and dispersion of this so-called mass-concentration relationship have been reported. Some observational studies of massive halos suggest that they are statistically over-concentrated with respect to the expectations of Lambda-CDM. Here we investigate the impact that various published mass-concentration relationships, both from simulations and derived from observations, would have on other cosmological observables, in particular considering upcoming surveys. We find that an integral measure of lensing shear, such as counts of peaks from halos, is very sensitive to the relationship between mass and concentration at fixed sigma-8, and the disparity between some reported fits is much larger than the impact of uncertainty in sigma-8 itself. We also briefly assess the impact of baryonic physics on cluster scale observables, using state-of-the-art simulations, concluding that it is unlikely to give rise to the high concentrations reported for some clusters.

In a recent publication, the flexion aperture mass statistic was found to provide a robust and effective method by which substructure in galaxy clusters might be mapped. Moreover, we suggested that the masses and mass profile of structures might be constrained using this method. In this paper, we apply the flexion aperture mass technique to HST ACS images of Abell 1689. We demonstrate that the flexion aperture mass statistic is sensitive to small-scale structures in the central region of the cluster. While the central potential is not constrained by our method, due largely to missing data in the central 0.5$^\prime$ of the cluster, we are able to place constraints on the masses and mass profiles of prominent substructures. We identify 4 separate mass peaks, and use the peak aperture mass signal and zero signal radius in each case to constrain the masses and mass profiles of these substructures. The three most massive peaks exhibit complex small-scale structure, and the masses indicated by the flexion aperture mass statistic suggest that these three peaks represent the dominant substructure component of the cluster ($\sim 7\times 10^{14}h^{-1}M_\odot$). Their complex structure indicates that the cluster -- far from being relaxed -- may have recently undergone a merger. The smaller, subsidiary peak is located coincident with a group of galaxies within the cluster, with mass $\sim 1\times10^{14}h^{-1}M_\odot$. These results are in excellent agreement with previous substructure studies of this cluster.

We propose a new ray-tracing algorithm to measure the weak lensing shear and convergence fields directly from N-body simulations. We calculate the deflection of the light rays lensed by the 3-D mass density field or gravitational potential along the line of sight on a grid-by-grid basis, rather than using the projected 2-D lens planes. Our algorithm uses simple analytic formulae instead of numerical integrations in the computation of the projected density field along the line of sight, and so is computationally efficient, accurate and straightforward to implement. This will prove valuable in the interpretation of data from the next generation of surveys that will image many thousands of square degrees of sky.

We show that simple models of scalar-field dark energy leave a generic enhancement in the weak-lensing power spectrum when compared to the LCDM prediction. In particular, we calculate the linear-scale enhancement in the convergence (or cosmic-shear) power spectrum for two best-fit models of scalar-field dark energy, namely, the Ratra-Peebles and SUGRA-type quintessence. Our calculations are based on linear perturbation theory, using gauge-invariant variables with carefully defined adiabatic initial conditions. We find that geometric effects enhance the lensing power spectrum on a broad range of scales, whilst the clustering of dark energy gives rise to additional power on large scales. The dark-energy power spectrum for these models are also explicitly obtained. On degree scales, the total enhancement may be as large as 30-40% for sources at redshift ~1. We argue that there are realistic prospects for detecting such an enhancement using the next generation of large telescopes.

Gravitational lensing induces significant errors in the measured distances to high-redshift standard candles and standard sirens such as type-Ia supernovae, gamma-ray bursts, and merging supermassive black hole binaries. There will therefore be a significant benefit from correcting for the lensing error by using independent and accurate estimates of the lensing magnification. We investigate how accurately the magnification can be inferred from convergence maps reconstructed from galaxy shear and flexion data. We employ ray-tracing through the Millennium Simulation to simulate lensing observations in large fields, and perform a weak-lensing reconstruction on these fields. We identify optimal ways to filter the reconstructed convergence maps and to convert them to magnification maps. We find that a shear survey with 100 galaxies/arcmin^2 can help to reduce the lensing-induced distance errors for standard candles/sirens at redshifts z=1.5 (z=5) on average by 20% (10%), whereas a futuristic survey with shear and flexion estimates from 500 galaxies/arcmin^2 yields much larger reductions of 50% (35%). For redshifts z>=3, a further improvement by 5% can be achieved, if the individual redshifts of the galaxies are used in the reconstruction. Moreover, the reconstruction allows one to identify regions for which the convergence is low, and in which an error reduction by up to 75% can be achieved.

Local non-Gaussianity causes correlations between large scale perturbation modes and the small scale power. The large-scale CMB signal has contributions from the integrated Sachs Wolfe (ISW) effect, which does not correlate with the small scale power. If this ISW contribution can be removed, the sensitivity to local non-Gaussianity is improved. Gravitational lensing and galaxy counts can be used to trace the ISW contribution; in particular we show that the CMB lensing potential is highly correlated with the ISW signal. We construct a nearly-optimal estimator for the local non-Gaussianity parameter $\fnl$ and investigate to what extent we can use this to decrease the variance on ${\fnl}$. We show that the variance can be decreased by up to $20\%$ at Planck sensitivity using galaxy counts. CMB lensing is a good bias-independent ISW tracer for future more sensitive observations, though the fractional decrease in variance is small if good polarization data is also available.

The sightline to the brighter member of the gravitationally lensed quasar pair UM 673A,B intersects a damped Lyman-alpha system (DLA) at z = 1.62650 which, because of its low redshift, has not been recognised before. Our high quality echelle spectra of the pair, obtained with HIRES on the Keck I telescope, show a drop in neutral hydrogen column density N(H I) by a factor of at least 400 between UM 673A and B, indicating that the DLA's extent in this direction is much less than the 2.7 kpc separation between the two sightlines at z = 1.62650. By reassessing this new case together with published data on other quasar pairs, we conclude that the typical size (radius) of DLAs at these redshifts is R ~ (5 +/- 3) kpc, smaller than previously realised. Highly ionized gas associated with the DLA is more extended, as we find only small differences in the C IV absorption profiles between the two sightlines. Coincident with UM 673B, we detect a weak and narrow Ly-alpha emission line which we attribute to star formation activity at a rate SFR >~ 0.2 M_solar/yr. The DLA in UM 673A is metal-poor, with an overall metallicity Z_DLA ~ 1/30 Z_solar, and has a very low internal velocity dispersion. It exhibits some apparent peculiarities in its detailed chemical composition, with the elements Ti, Ni, and Zn being deficient relative to Fe by factors of 2-3. The [Zn/Fe] ratio is lower than those measured in any other DLA or Galactic halo star, presumably reflecting somewhat unusual previous enrichment by stellar nucleosynthesis. We discuss the implications of these results for the nature of the galaxy hosting the DLA.

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 investigate the impact of AGN feedback on the gravitational lensing properties of a sample of galaxy clusters with masses in the range 10^14 - 10^15 Msol, using state-of-the-art simulations. Adopting a ray-tracing algorithm, we compute the cross-section of giant arcs from clusters simulated with dark-matter only physics (DM), dark matter plus gas with cooling and star formation (CSF), and dark matter plus gas with cooling, star formation and AGN feedback (CSFBH). Once AGN feedback is included, baryonic physics boosts the strong lensing cross-section by much less than previously estimated using clusters simulated with only cooling and star formation. For a cluster with a virial mass of 7.4 X 10^14 Msol, inclusion of baryonic physics without feedback can boost the cross-section by as much as a factor of 3, in agreement with previous studies, whereas once AGN feedback is included this maximal figure falls to a factor of 2 at most. Typically, clusters simulated with DM and CSFBH physics have similar cross-sections for the production of giant arcs. We also investigate how baryonic physics affects the weak lensing properties of the simulated clusters by fitting NFW profiles to synthetic weak lensing data sets using a Markov Chain Monte Carlo approach, and by performing non-parametric mass reconstructions. Without the inclusion of AGN feedback, measured concentration parameters can be much larger than those obtained with AGN feedback, which are similar to the dark-matter only case.

In a recent publication, an aperture mass statistic for gravitational flexion was derived and shown to be effective, at least with simulated data, in detecting massive structures and substructures within clusters of galaxies. Further, it was suggested that the radius at which the flexion aperture mass signal falls to zero might allow for estimation of the mass or density profile of the structures detected. In this paper, we more fully explore this possibility, considering the behaviour both of the peak signal and the zero-signal contours for two mass models--the singular isothermal sphere and Navarro-Frenk-White profiles--under varying aperture size, filter shape and mass concentration parameter. We demonstrate the effectiveness of the flexion aperture mass statistic in discriminating between mass profiles and concentration parameters, and in providing an accurate estimate of the mass of the lens, to within a factor of 1.5 or better. In addition, we compare the aperture mass method to a direct nonparametric reconstruction of the convergence from flexion measurements. We demonstrate that the aperture mass technique is much better able to constrain the shape of the central density profile, obtains much finer angular resolution in reconstructions, and does not suffer from ambiguity in the normalisation of the signal, in contrast to the direct method.

The problem of detecting dark matter filaments in the cosmic web is considered. Weak lensing is an ideal probe of dark matter, and therefore forms the basis of particularly promising detection methods. We consider and develop a number of weak lensing techniques that could be used to detect filaments in individual or stacked cluster fields, and apply them to synthetic lensing data sets in the fields of clusters from the Millennium Simulation. These techniques are multipole moments of the shear and convergence, mass reconstruction, and parameterized fits to filament mass profiles using a Markov Chain Monte Carlo approach. In particular, two new filament detection techniques are explored (multipole shear filters and Markov Chain Monte Carlo mass profile fits), and we outline the quality of data required to be able to identify and quantify filament profiles. We also consider the effects of large scale structure on filament detection. We conclude that using these techniques, there will be realistic prospects of detecting filaments in data from future space-based missions. The methods presented in this paper will be of great use in the identification of dark matter filaments in future surveys.

Upcoming ground and space based observatories such as the DES, the LSST, the JDEM concepts and the SKA, promise to dramatically increase the size of strong gravitational lens samples. A significant fraction of the systems are expected to be time delay lenses. Many of the existing lensing degeneracies become less of an issue with large samples since the distributions of a number of parameters are predictable, and can be incorporated into an analysis, thus helping to lessen the degeneracy. Assuming a mean galaxy density profile that does not evolve with redshift, a Lambda-CDM cosmology, and Gaussian distributions for bulk parameters describing the lens and source populations, we generate synthetic lens catalogues and examine the relationship between constraints on the Omega_m - Omega_Lambda plane and H_0 with increasing lens sample size. We find that, with sample sizes of ~400 time delay lenses, useful constraints can be obtained for Omega_m and Omega_Lambda with approximately similar levels of precision as from the best of other methods. In addition, sample sizes of ~100 time delay systems yield estimates of H_0 with errors of only a couple of percent, exceeding the level of precision from current best estimates such as the HST Key Project. We note that insufficient prior knowledge of the lens samples employed in the analysis, via under or overestimates in the mean values of the sample distributions, results in broadening of constraints. This highlights the need for sound prior knowledge of the sample before useful cosmological constraints can be obtained from large time delay samples (abridged).

Gravitational flexion has recently been introduced as a technique by which one can map out and study substructure in clusters of galaxies. Previous analyses involving flexion have measured the individual galaxy-galaxy flexion signal, or used either parametric techniques or a KSB-type inversion to reconstruct the mass distribution in Abell 1689. In this paper, we present an aperture mass statistic for flexion, and apply it to the lensed images of background galaxies obtained by ray-tracing simulations through a simple analytic mass distribution and through a galaxy cluster from the Millennium simulation. We show that this method is effective at detecting and accurately tracing structure within clusters of galaxies on sub-arcminute scales with high signal-to-noise even using a moderate background source number density and image resolution. In addition, the method provides much more information about both the overall shape and the small-scale structure of a cluster of galaxies than can be achieved through a weak lensing mass reconstruction using gravitational shear data. Lastly, we discuss how the zero-points of the aperture mass might be used to infer the masses of structures identified using this method.

Accurate measurement of the cluster mass function is a crucial element in efforts to constrain structure formation models, the normalisation of the matter power spectrum and the cosmological matter density, and the nature and evolution of dark energy. Large weak lensing surveys of ~20,000 galaxy clusters and groups will be key tools in the observational pursuit of that goal. These weak lensing studies often proceed by stacking the lensing signals of many clusters and groups binned by mass-correlated observables such as richness and luminosity; typically such analyses ignore the triaxial structure of dark matter halos on the assumption that the averaging of many shear signals within each mass bin makes its effects (as large as factors of two in mass model parameter estimates from individual clusters) negligible. We test this assumption and find that triaxiality can bias 3D virial mass estimates compared to those for a spherical population by a few percent if suboptimal mass estimators are used. This bias affects not only direct lensing constraints on the mass function but can also affect the scatter and normalization of the mass-observable relations derived from lensing that are so crucial to constraining the cluster mass function with large samples. However, we demonstrate that a careful choice of mass estimator can remove the bias very effectively if the lensing signals from a sufficient number of triaxial halos are averaged together, and further quantify that sufficient number for adequate shape averaging. We thus show that by choosing observable bins to contain an adequate number of halos and by utilizing a carefully chosen 3D mass estimator stacked weak-lensing analyses can give unbiased constraints on the triaxial mass function.

Measuring the 3D distribution of mass on galaxy cluster scales is a crucial test of the LCDM model, providing constraints on the nature of dark matter. Recent work investigating mass distributions of individual galaxy clusters (e.g. Abell 1689) using weak and strong gravitational lensing has revealed potential inconsistencies between the predictions of structure formation models relating halo mass to concentration and those relationships as measured in massive clusters. However, such analyses employ simple spherical halo models while a growing body of work indicates that triaxial 3D halo structure is both common and important in parameter estimates. We recently introduced a Markov Chain Monte Carlo (MCMC) method to fit fully triaxial models to weak lensing data that gives parameter and error estimates that fully incorporate the true shape uncertainty present in nature. In this paper we apply that method to weak lensing data obtained with the ESO/MPG Wide-Field Imager for galaxy clusters A1689, A1835, and A2204, under a range of Bayesian priors derived from theory and from independent X-ray and strong lensing observations. For Abell 1689, using a simple strong lensing prior we find marginalized mean parameter values M_200 = (0.83 +- 0.16)x10^15 M_solar/h and C=12.2 +- 6.7, which are marginally consistent with the mass-concentration relation predicted in LCDM. The large error contours that accompany our triaxial parameter estimates more accurately represent the true extent of our limited knowledge of the structure of galaxy cluster lenses, and make clear the importance of combining many constraints from other theoretical, lensing (strong, flexion), or other observational (X-ray, SZ, dynamical) data to confidently measure cluster mass profiles. (Abridged)

Measuring the 3D distribution of mass on galaxy cluster scales is a crucial test of the LCDM model, providing constraints on the behaviour of dark matter. Recent work investigating mass distributions of individual galaxy clusters (e.g. Abell 1689) using weak and strong gravitational lensing has revealed potential inconsistencies between the predictions of structure formation models relating halo mass to concentration and those relationships as measured in massive clusters. However, such analyses employ simple spherical halo models while a growing body of work indicates that triaxial 3D halo structure is both common and important in parameter estimates. The very strong assumptions about the symmetry of the lensing halo implied with circular or perturbative elliptical NFW models are not physically motivated and lead to incorrect parameter estimates with significantly underestimated error bars. We here introduce a Markov Chain Monte Carlo (MCMC) method to fit fully triaxial models to weak lensing data that gives parameter and error estimates that fully incorporate the true uncertainty present in nature. Applying the MCMC triaxial fitting method to a population of NFW triaxial lenses drawn from the shape distribution of structure formation simulations, we find that including triaxiality cannot explain a population of massive, highly concentrated clusters within the framework of LCDM, but easily explains rare cases of apparently massive, highly concentrated, very efficient lensing clusters. Our MCMC triaxial NFW fitting method is easily expandable to include constraints from additional data types, and its application returns model parameters and errors that more accurately capture the true (and limited) extent of our knowledge of the structure of galaxy cluster lenses. (abridged)

We present discovery images, together with follow-up imaging and spectroscopy, of two large separation gravitational lenses found by our survey for wide arcs (the CASSOWARY). The survey exploits the multicolor photometry of the Sloan Digital Sky Survey to find multiple blue components around red galaxies. CASSOWARY~2 (or "the Cheshire Cat") is composed of two massive early-type galaxies at z = 0.426 and 0.432 respectively lensing two background sources, the first a star-forming galaxy at z = 0.97 and the second a high redshift galaxy (z> 1.4). There are at least three images of the former source and probably four or more of the latter, arranged in two giant arcs. The mass enclosed within the larger arc of radius 11 arcsecs is about 33 x 10^{12} solar masses. CASSOWARY~3 comprises an arc of three bright images of a z = 0.725 source, lensed by a foreground elliptical at z = 0.274. The radius of the arc is about 4 arcsecs and the enclosed mass is 2.5 x 10^{12} solar masses. Together with earlier discoveries like the Cosmic Horseshoe and the 8 O'Clock Arc, these new systems, with separations intermediate between the arcsecond separation lenses of typical strong galaxy lensing and the larger separation cluster lenses, probe the very high end of the galaxy mass function.

The Hubble constant can be constrained using the time delays between multiple images of gravitationally lensed sources. In some notable cases, typical lensing analyses assuming isothermal galaxy density profiles produce low values for the Hubble constant, inconsistent with the result of the HST Key Project (72 +- 8 km/s/Mpc). Possible systematics in the values of the Hubble constant derived from galaxy lensing systems can result from a number of factors, e.g. neglect of environmental effects, assumption of isothermality, or contamination by line-of-sight structures. One additional potentially important factor is the triaxial structure of the lensing galaxy halo; most lens models account for halo shape simply by perturbing the projected spherical lensing potential, an approximation that is often necessary but that is inadequate at the levels of triaxiality predicted in the CDM paradigm. To quantify the potential error introduced by this assumption in estimates of the Hubble parameter, we strongly lens a distant galaxy through a sample of triaxial softened isothermal halos and use an MCMC method to constrain the lensing halo profile and the Hubble parameter from the resulting multiple image systems. We explore the major degeneracies between the Hubble parameter and several parameters of the lensing model, finding that without a way to accurately break these degeneracies accurate estimates of the Hubble parameter are not possible. Crucially, we find that triaxiality does not significantly bias estimates of the Hubble constant, and offer an analytic explanation for this behaviour in the case of isothermal profiles. Neglected triaxial halo shape cannot contribute to the low Hubble constant values derived in a number of galaxy lens systems.

Cosmic shear constrains cosmology by exploiting the apparent alignments of pairs of galaxies due to gravitational lensing by intervening mass clumps. However galaxies may become (intrinsically) aligned with each other, and with nearby mass clumps, during their formation. This effect needs to be disentangled from the cosmic shear signal to place constraints on cosmology. We use the linear intrinsic alignment model as a base and compare it to an alternative model and data. If intrinsic alignments are ignored then the dark energy equation of state is biased by ~50 per cent. We examine how the number of tomographic redshift bins affects uncertainties on cosmological parameters and find that when intrinsic alignments are included two or more times as many bins are required to obtain 80 per cent of the available information. We investigate how the degradation in the dark energy figure of merit depends on the photometric redshift scatter. Previous studies have shown that lensing does not place stringent requirements on the photometric redshift uncertainty, so long as the uncertainty is well known. However, if intrinsic alignments are included the requirements become a factor of three tighter. These results are quite insensitive to the fraction of catastrophic outliers, assuming that this fraction is well known. We show the effect of uncertainties in photometric redshift bias and scatter. Finally we quantify how priors on the intrinsic alignment model would improve dark energy constraints.

We report the discovery of an almost complete Einstein ring of diameter 10" in Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5). Spectroscopic data from the 6m telescope of the Special Astrophysical Observatory reveals that the deflecting galaxy has a line-of-sight velocity dispersion in excess of 400 km/s and a redshift of 0.444, whilst the source is a star-forming galaxy with a redshift of 2.379. From its color and luminosity, we conclude that the lens is an exceptionally massive Luminous Red Galaxy (LRG) with a mass within the Einstein radius of 5 x 10^12 solar masses. This remarkable system provides a laboratory for probing the dark matter distribution in LRGs at distances out to 3 effective radii, and studying the properties of high redshift star-forming galaxies.

We show that the existence of cosmic strings can be strongly constrained by the next generation of gravitational lensing surveys at radio frequencies. We focus on cosmic string loops, which simulations suggest would be far more numerous than long (horizon-sized) strings. Using simple models of the loop population and minimal assumptions about the lensing cross section per loop, we estimate the optical depth to lensing and show that extant radio surveys such as CLASS have already ruled out a portion of the cosmic string model parameter space. Future radio interferometers, such as LOFAR and especially SKA, may constrain $G\mu/c^2 < 10^{-9}$ in some regions of parameter space, outperforming current constraints from pulsar timing and the CMB by up to two orders of magnitude. This method relies on direct detections of cosmic strings, and so is less sensitive to the theoretical uncertainties in string network evolution that weaken other constraints.

Measuring the distribution of mass on galaxy cluster scales is a crucial test of the Lambda CDM model, providing constraints on the nature of dark matter. Recent work investigating mass distributions of individual galaxy clusters using gravitational lensing has illuminated potential inconsistencies between the predictions of structure formation models relating halo mass to concentration and those relationships as measured in massive clusters. However, such analyses typically employ only simple spherical halo models with canonical NFW slopes, while the halos formed in simulations show a range of more complex features. Here we investigate the impact of such expected deviations from the canonical NFW halo profile on mass and parameter estimation using weak gravitational lensing on massive cluster scales. The most important deviation is halo triaxiality because it is impossible even with fiducial weak lensing data to fully resolve the three-dimensional structure of the halo due to lensing's sensitivity only to projected mass. Significant elongation of the halo along the line of sight can cause the mass and concentration to be overestimated by as much as 50% and by a factor of 2, respectively, while foreshortening has the opposite effect. Additionally, triaxial halos in certain orientations are much better lenses than their spherical counterparts of the same mass, indicating that clusters chosen for study because of evident lensing are likely to be drawn from the high-triaxiality end of the halo shape distribution. While the effects of triaxiality alone may not be enough to fully explain the very high concentrations reported for some clusters, they go a long way in easing the tensions between observations and the predictions of the cold dark matter paradigm. (Abridged)

We continue to see a range of values for the Hubble constant obtained from gravitationally lensed multiple image time delays when assuming an isothermal lens despite a robust value from the HST key project (72 +- 8 km s^-1 Mpc^-1. One explanation is that there is a variation in Hubble constant values due to a fundamental heterogeneity in lens galaxies present in groups. Our goal is to see if a variety of group interactions between the most massive group members can result in significant changes in the galaxy density profiles over the scale probed by strong lensing (<15 kpc). While stripping of the outer parts of the halo can be expected, the impact on inner regions where the luminous component is important is less clear in the context of lensing, though still crucial, as a steepened density profile within this inner region allows these lens systems to be consistent with current HST/WMAP estimates on H_0. We employ the particle-mesh code SUPERBOX to carry out the group interaction simulations. We simulate interactions between group members, comparing the density profile for the satellite before and after interaction for the mass range of 10^11 to 10^13 M_sun. Our investigations show a significant steepening of the density profile in the region of < 5-20 kpc, i.e. that which dominates strong lensing in lens galaxies. Additionally, the steepening in the inner region is transient in nature, with consecutive interactions returning the profile to an isothermal state within a timeframe of ~ 0.5 - 2.0 Gyr. This factor may help explain why lens galaxies that produce lower values of H_0 (i.e. those with possibly steeper profiles) are far fewer in number than those which agree with both the HST key project value for H_0 and isothermality, since one would have to observe the lens galaxy during this transient steepened phase.

The distribution of mass on galaxy cluster scales is an important test of structure formation scenarios, providing constraints on the nature of dark matter itself. Several techniques have been used to probe the mass distributions of clusters, sometimes yielding results which are discrepant, or at odds with clusters formed in simulations - for example giving NFW concentration parameters much higher than expected in the standard CDM model. In addition, the velocity fields of some well studied galaxy clusters reveal the presence of several structures close to the line-of-sight, often not dynamically bound to the cluster itself. We investigate what impact such neighbouring but unbound massive structures would have on the determination of cluster profiles using weak gravitational lensing. Depending on its concentration and mass ratio to the primary halo, one secondary halo close to the line-of-sight can cause the estimated NFW concentration parameter to be significantly higher than that of the primary halo, and also cause the estimated mass to be biased high. Although it is difficult to envisage how this mechanism alone could yield concentrations as high as reported for some clusters, multiple haloes close to the line-of-sight, such as in the case of Abell 1689, can substantially increase the concentration parameter estimate. Together with the fact that clusters are triaxial, and that including baryonic physics also leads to an increase in the concentration of a dark matter halo, the tension between observations and the standard CDM model is eased. If the alignment with the secondary structure is imprecise, then the estimated concentration parameter can also be even lower than that of the primary halo, reinforcing the importance of identifying structures in cluster fields.

Gravitational lensing time delays depend upon the Hubble constant and the density distribution of the lensing galaxies. This allows one to either model the lens and estimate the Hubble constant, or to use a prior on the Hubble constant from other studies and investigate what the preferred density distribution is. Some studies have required compact dark matter halos (constant M/L ratio) in order to reconcile gravitational lenses with the HST/WMAP value of the Hubble constant (72 +/- 8 km/s /Mpc and 72 +/- 5 km/s /Mpc, respectively). This is in direct contradiction with X-ray, stellar dynamical, and weak lensing studies, which all point towards extended halos and isothermal density profiles. In this work, we examine an up-to-date sample of 13 lensing galaxies resulting in a data set consisting of 21 time delays. We select systems in which there is a single primary lensing galaxy (e.g. excluding systems undergoing mergers). Analysis is performed using analytic models based upon a powerlaw density profile (rho \propto r^-n) of which the isothermal profile is a special case (n = 2). This yields a value of n = 2.11+/-0.12 (3sigma) for the mean profile when modeling with a prior on the Hubble constant, which is only consistent with isothermality within 3 sigma. Note that this is a formal error from our calculations, and does not include the impact of sample selection or simplifications in the lens modeling. We conclude that time delays are a useful probe of density profiles, in particular as a function of the environment in which the lens resides, when combined with a prior on the Hubble constant.

Gravitational lensing can be used to directly constrain the projected density profile of galaxy clusters. We discuss possible future constraints from lensing of the CMB temperature and polarization, and compare to results from galaxy weak lensing. We model the moving lens and kinetic SZ signals that confuse the temperature CMB lensing when cluster velocities and angular momenta are unknown, and show how they degrade parameter constraints. The CMB polarization cluster lensing signal is ~1 micro-Kelvin for massive clusters and challenging to detect; however it should be significantly cleaner than the temperature signal and may provide the most robust constraints at low noise levels. Galaxy lensing is likely to be much better for constraining cluster masses at low redshift, but for clusters at redshift z >~ 1 future CMB lensing observations may be able to do better.

Cosmic shear leads to a correlation of the observed ellipticities of galaxies, an effect which is used to place constraints on cosmological parameters, and to explore the evolution of dark matter and dark energy in the universe. However, a possible systematic contaminant of the lensing signal is intrinsic galaxy alignment, with a correlation length of a few Mpc. Hirata & Seljak (2004) have recently demonstrated that for some models of intrinsic distortions, there may also be a cross-correlation between the intrinsic and lensing signals, which may dominate the intrinsic signal, and suppress the lensing power spectrum by several tens of percent. Unlike the pure intrinsic signal, this new term cannot be accounted for by neglecting or down-weighting pairs of galaxies which are physically close. Extending the correlation function tomography method of King & Schneider (2003) we illustrate how the impact of both intrinsic and cross-correlations can be significantly reduced, in the context of surveys with photometric redshift information. For a ground-based cosmic shear survey of ~100 sq. degrees with photometric redshifts, even in the presence of systematic contaminants at the level considered here, cosmological models degenerate in the Omega_m-sigma_8 plane can be distinguished well in excess of the 3-sigma level.

Weak gravitational lensing can be used to directly measure the mass along a line-of-sight without any dependence on the dynamical state of the mass, and thus can be used to measure the masses of clusters even if they are not relaxed. One common technique used to measure cluster masses is fitting azimuthally-averaged gravitational shear profiles with a spherical mass model. In this paper we quantify how asphericity and projected substructure in clusters can affect the virial mass and concentration measured with this technique by simulating weak lensing observations on 30 independent lines-of-sights through each of four high-resolution N-body cluster simulations. We find that the variations in the measured virial mass and concentration are of a size similar to the error expected in ideal weak lensing observations and are correlated, but that the virial mass and concentration of the mean shear profile agree well with that measured in three dimensional models of the clusters. The dominant effect causing the variations is the proximity of the line-of-sight to the major axis of the 3-D cluster mass distribution, with projected substructure only causing minor perturbations in the measured concentration. Finally we find that the best-fit "universal" CDM models used to fit the shear profiles over-predict the surface density of the clusters due to the cluster mass density falling off faster than the r^{-3} model assumption.

Cosmological weak lensing by the large scale structure of the Universe, cosmic shear, is coming of age as a powerful probe of the parameters describing the cosmological model and matter power spectrum. It complements CMB studies, by breaking degeneracies and providing a cross-check. An important measure of the cosmic shear signal are the shear correlation functions; these can be directly calculated from data, and compared with theoretical expectations for different cosmological models and matter power spectra. We present a Monte Carlo method to quickly simulate mock cosmic shear surveys. One application of this method is in the determination of the full covariance matrix for the correlation functions; this includes redshift binning and is applicable to arbitrary survey geometries. Terms arising from shot noise and cosmic variance (dominant on small and large scales respectively) are accounted for naturally. As an illustration of the use of such covariance matrices, we consider to what degree confidence regions on parameters are tightened when redshift binning is employed. The parameters considered are those commonly discussed in cosmic shear analyses - the matter density parameter, dark energy density parameter (classical cosmological constant), power spectrum normalisation and shape parameter. We incorporate our covariance matrices into a likelihood treatment, and also use the Fisher formalism to explore a larger region of parameter space (abridged).

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)

During the past few years, secure detections of cosmic shear have been obtained, manifest in the correlation of the observed ellipticities of galaxies. Constraints have already been placed on cosmological parameters, such as the normalisation of the matter power spectrum sigma_8. One possible systematic contaminant of the lensing correlation signal arises from intrinsic galaxy alignment, which is still poorly constrained. Unlike lensing, intrinsic correlations only pertain to galaxies with small physical separations, the correlation length being a few Mpc. We present a new method that harnesses this property, and isolates the lensing and intrinsic components of the galaxy ellipticity correlation function using measurements between different redshift slices. The observed signal is approximated by a set of template functions, making no strong assumptions about the amplitude or correlation length of any intrinsic alignment. We also show that the near-degeneracy between the matter density parameter Omega_m and sigma_8 can be lifted using correlation function tomography, even in the presence of an intrinsic alignment signal.

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

Cosmological weak lensing gives rise to correlations in the ellipticities of faint galaxies. This cosmic shear signal depends upon the matter power spectrum, thus providing a means to constrain cosmological parameters. It has recently been proposed that intrinsic alignments arising at the epoch of galaxy formation can also contribute significantly to the observed correlations, the amplitude increasing with decreasing survey depth. Here we consider the two-point shear correlation function, and demonstrate that photometric redshift information can be used to suppress the intrinsic signal; at the same time Poisson noise is increased, due to a decrease in the effective number of galaxy pairs. The choice to apply such a redshift-depending weighting will depend on the characteristics of the survey in question. In surveys with a mean z of about 1, although the lensing signal dominates, the measurement error bars may soon become smaller than the intrinsic alignment signal; hence, in order not to be dominated by systematics, redshift information in cosmic shear statistics will become a necessity. We discuss various aspects of this.

Gravitationally lensed systems with multiply imaged quasars are an excellent tool for studying the properties of distant galaxies. In particular, they provide the most accurate mass measures for the lensing galaxy. The system B1422+231 is a well studied example of a quadruply imaged quasar, with high-quality VLBI data available. Very accurate data on image positions, fluxes and deconvolved image sizes provide good constraints for lensing models. We discuss here the failure of smooth models in fitting the data. Since it is intuitively clear that the mass of a lens galaxy is not a smooth entity, we have investigated how deviation from a smooth model can influence lensing phenomena, especially the image flux ratios. To explore expectations about the level of substructure in galaxies and its influence on strong lensing, N-body simulations of a model galaxy are employed. By using the mass distribution of this model galaxy as a lens, synthetic data sets of different four image system configurations are generated. Their analysis can possibly provide evidence for the presence and strong influence of substructure in the primary lens galaxy.

In this work we investigate the gravitationally lensed system B1422+231. High--quality VLBI image positions, fluxes and shapes as well as an optical HST lens galaxy position are used. First, two simple and smooth models for the lens galaxy are applied to fit observed image positions and fluxes; no even remotely acceptable model was found. Such models also do not accurately reproduce the image shapes. In order to fit the data successfully, mass substructure has to be added to the lens, and its level is estimated. To explore expectations about the level of substructure in galaxies and its influence on strong lensing, N-body simulation results of a model galaxy are employed. By using the mass distribution of this model galaxy as a lens, synthetic data sets of different four image system configurations are generated and simple lens models are again applied to fit them. The difficulties in fitting these lens systems turn out to be similar to the case of some real gravitationally lensed systems, thus possibly providing evidence for the presence and strong influence of substructure in the primary lens galaxy.

We present the first detection of weak gravitational shear at infrared wavelengths, using observations of the lensing cluster Abell 1689, taken with the SofI camera on the ESO-NTT telescope. The imprint of cluster lenses on the shapes of the background galaxy population has previously been harnessed at optical wavelengths, and this gravitational shear signal enables cluster mass distributions to be probed, independent of whether the matter is luminous or dark. At near-infrared wavelengths, the spectrophotometric properties of galaxies facilitate a clean selection of background objects for use in the lensing analysis. A finite-field mass reconstruction and application of the aperture mass (Map) statistic are presented. The probability that the peak of the Map detection S/N~5, arises from a chance alignment of background sources is only ~4.5*10^-7. The velocity dispersion of the best-fit singular isothermal sphere model for the cluster is sigma_1D=1030^{+70}_{-80} km/s, and we find a K-band mass-to-light ratio of ~40 M_solar/L_solar inside a 0.44 Mpc radius.

Here we apply a recently developed maximum likelihood method for determining best-fit parameterised lens models to observations of the rich lensing cluster Abell 1689. The observations that we use were taken with the ESO/MPG Wide Field Imager. The wide field-of-view enables us to use the weakly lensed images of faint background objects on an unsurpassed range of scales, 0.12Mpc/h<R<1.8Mpc/h from the cluster centre, to determine best-fit models for the 1-parameter singular isothermal sphere (SIS), 2-parameter general power-law and NFW models, and 3-parameter singular isothermal ellipsoid (SIE). The best-fit SIS has an Einstein radius of 0.37arcmin (0.043Mpc/h) i.e. a velocity dispersion of 1028(+35,-42)km/s in an Omega=1.0, Lambda=0.0 cosmology. For the best-fit NFW profile, the virial radius r_200=1.14Mpc/h and the concentration parameter c=4.7, giving a virial mass M_200=5.7x10**14 M_solar/h. At q=0.88, the slope of the best-fit power-law model is slightly flatter than isothermal (q=1.0), indicating that the galaxies most important in the fitting procedure lie inside the scale radius r_s. By fitting an SIE, the deviation of the projected mass distribution from circular symmetry is evident, with a best-fit axial ratio f=0.74.

Integral measures of weak gravitational lensing by dark-matter haloes, like the aperture mass, are sensitive to different physical halo properties dependent on the halo mass density profile. For isothermal profiles, the relation between aperture mass and virial mass is steeper than for haloes with the universal NFW profile. Consequently, the halo mass range probed by the aperture mass is much wider for NFW than for isothermal haloes. We use recent modifications to the Press-Schechter mass function in CDM models normalised to the local abundance of rich clusters, to predict the properties of the halo sample expected to be accessible with the aperture mass technique. While ~10 haloes should be detected per square degree if the haloes have NFW profiles, their number density is lower by approximately an order of magnitude if they have isothermal profiles. These results depend only very mildly on the cosmological background model. We conclude that counts of haloes with a significant weak-lensing signal are a powerful discriminator between different dark-matter profiles.

In Schneider, King & Erben (2000) we developed likelihood techniques to compare the constraints on cluster mass profiles that can be obtained using the shear and magnification information. This work considered circularly symmetric power-law models for clusters at fairly low redshifts where the redshift distribution of source galaxies could be neglected. Here this treatment is extended to encompass NFW profiles which are a good description of clusters from cosmological N-body simulations, and NFW clusters at higher redshifts where the influence of various scenarios for the knowledge of the redshift distribution are examined. Since in reality the overwhelming majority of clusters have ellipsoidal rather than spherical profiles, the singular isothermal ellipsoid (SIE) is investigated. We also briefly consider the impact of substructure on such a likelihood analysis. In general, we find that the shear information provides a better constraint on the NFW profile under consideration, so this becomes the focus of what follows. The ability to differentiate between the NFW and power-law profiles strongly depends on the size of the data field, and on the number density of galaxies for which an ellipticity can be measured. For higher redshift NFW profiles, there is very little reduction (~1.5%) in the dispersion of parameter estimates when spectroscopic redshifts, as opposed to photometric redshift estimates, are available for the galaxies used in the lensing analysis.

We present observations of microlensing event MACHO-98-BLG-35 which reached a peak magnification factor of almost 80. These observations by the Microlensing Planet Search (MPS) and the MOA Collaborations place strong constraints on the possible planetary system of the lens star and show intriguing evidence for a low mass planet with a mass fraction $4\times 10^{-5} \leq \epsilon \leq 2\times 10^{-4}$. A giant planet with $\epsilon = 10^{-3}$ is excluded from 95% of the region between 0.4 and 2.5 $R_E$ from the lens star, where $R_E$ is the Einstein ring radius of the lens. This exclusion region is more extensive than the generic "lensing zone" which is $0.6 - 1.6 R_E$. For smaller mass planets, we can exclude 57% of the "lensing zone" for $\epsilon = 10^{-4}$ and 14% of the lensing zone for $\epsilon = 10^{-5}$. The mass fraction $\epsilon = 10^{-5}$ corresponds to an Earth mass planet for a lensing star of mass $\sim 0.3 \msun$. A number of similar events will provide statistically significant constraints on the prevalence of Earth mass planets. In order to put our limits in more familiar terms, we have compared our results to those expected for a Solar System clone averaging over possible lens system distances and orientations. We find that such a system is ruled out at the 90% confidence level. A copy of the Solar System with Jupiter replaced by a second Saturn mass planet can be ruled out at 70% confidence. Our low mass planetary signal (few Earth masses to Neptune mass) is significant at the $4.5\sigma$ confidence level. If this planetary interpretation is correct, the MACHO-98-BLG-35 lens system constitutes the first detection of a low mass planet orbiting an ordinary star without gas giant planets.

We present a 9 million star color-magnitude diagram (9M CMD) of the LMC bar. The 9M CMD reveals a complex superposition of different age and metallicity stellar populations, with important stellar evolutionary phases occurring over 3 orders of magnitude in number density. First, we count the non-variable supergiants, the associated Cepheids, and measure the effective temperatures defining the instability strip. Lifetime predictions of stellar evolution theory are tested, with implications for the origin of low-luminosity Cepheids. The highly-evolved AGB stars have a bimodal distribution in brightness, which we interpret as discrete old populations (>1 Gyr). The faint AGB may be metal-poor and very old. We identify the clusters NGC 411 and M3 as templates for the admixture of old stellar populations. However, there are indications that the old and metal-poor field population has a red HB morphology: the RR Lyraes lie on the red edge of the instability strip, the AGB-bump is very red, and the ratio of AGB-bump stars to RR Lyraes is quite large. If the HB second parameter is age, the old and metal-poor field population likely formed after the oldest clusters. Lifetime predictions of stellar evolution theory lead us to associate a significant fraction of the red HB clump giants with the same old and metal-poor population producing the RR Lyraes and the AGB-bump. In this case, compared to the age-dependent luminosity predictions of stellar evolution theory, the red HB clump is too bright relative to the RR Lyraes and AGB-bump. Last, the surface density profile of RR Lyraes is fit by an exponential, favoring a disk-like rather than spheroidal distribution. We conclude that the age of the LMC disk is probably similar to the age of the Galactic disk. (ABRIDGED)

A massive foreground cluster lens changes the shapes (shear effect) and number density (magnification effect) of the faint background galaxy population. In this paper we investigate how the shear, magnification and combined information can be used to constrain cluster mass profiles in the weak lensing regime. We develop maximum likelihood techniques to make quantitative predictions for each of the methods. Our analytic results are checked against Monte Carlo simulations. In general, we find that the shear method is superior to the magnification method. However, the magnification information complements the shear information if the former has accurate external calibration. For the magnification method, we discuss the effects of random and systematic uncertainties in the background galaxy counts.

We present the observations of the binary lensing event MACHO-98-SMC-1 conducted at the Mt.~Stromlo 74" telescope by the Microlensing Planet Search (MPS) collaboration. The MPS data constrain the first caustic crossing to have occurred after 1998 June 5.55 UT and thus directly rule out one of the two fits presented by the PLANET collaboration (model II). This substantially reduces the uncertainty in the the relative proper motion estimations of the lens object. We perform joint binary microlensing fits of the MPS data together with the publicly available data from the EROS, MACHO/GMAN and OGLE collaborations. We also study the binary lens fit parameters previously published by the PLANET and MACHO/GMAN collaborations by using them as initial values for $\chi^2$ minimization. Fits based on the PLANET model I appear to be in conflict with the GMAN-CTIO data. From our best fit, we find that the lens system has a proper motion of $\mu = 1.3\pm 0.2 \kmsk$ with respect to the source, which implies that the lens system is most likely to be located in the Small Magellanic Cloud strengthening the conclusion of previous reports.

We present photometric observations and analysis of the second microlensing event detected towards the Small Magellanic Cloud (SMC), MACHO Alert 98-SMC-1. This event was detected early enough to allow intensive observation of the lightcurve. These observations revealed 98-SMC-1 to be the first caustic crossing, binary microlensing event towards the Magellanic Clouds to be discovered in progress. Frequent coverage of the evolving lightcurve allowed an accurate prediction for the date of the source crossing out of the lens caustic structure. The caustic crossing temporal width, along with the angular size of the source star, measures the proper motion of the lens with respect to the source, and thus allows an estimate of the location of the lens. Lenses located in the Galactic halo would have a velocity projected to the SMC of v^hat ~1500 km/s, while an SMC lens would typically have v^hat ~60 km/s. We have performed a joint fit to the MACHO/GMAN data presented here, including recent EROS data of this event. These joint data are sufficient to constrain the time for the lens to move an angle equal to the source angular radius; 0.116 +/- 0.010 days. We estimate a radius for the lensed source of 1.4 +/- 0.1 R_sun. This yields a projected velocity of v^hat = 84 +/- 9 km/s. Only 0.15% of halo lenses would be expected to have a v^hat value at least as small as this, while 31% of SMC lenses would be expected to have v^hat as large as this. This implies that the lensing system is more likely to reside in the SMC than in the Galactic halo.

We report the discovery of a new double image gravitational lens system B1030+074 which was found during the Jodrell Bank - VLA Astrometric Survey (JVAS). We have collected extensive radio data on the system using the VLA, MERLIN, the EVN and the VLBA and optical observations using WFPC2 on the HST. The lensed images are separated by 1.56 arcseconds and their flux density ratio at centimetric wavelengths is approximately 14:1 although the ratio is slightly frequency dependent and the images appear to be time variable. The HST pictures show both the lensed images and the lensing galaxy close to the weaker image. The lensing galaxy has substructure which could be a spiral arm or an interacting galaxy.

We report the discovery, using NICMOS on the Hubble Space Telescope, of an arcsecond-diameter Einstein ring in the gravitational lens system B1938+666. The lensing galaxy is also detected, and is most likely an early-type. Modelling of the ring is presented and compared with the radio structure from MERLIN maps. We show that the Einstein ring is consistent with the gravitational lensing of an extended infrared component, centred between the two radio components.

The extreme bolometric luminosity of IRAS F10214+4724, and in particular the huge mass in molecular gas inferred from the CO line fluxes have led to suggestions that this is a giant galaxy in the process of formation. An arc-like structure and the closeness of a second object suggest, however, that gravitational lensing may be responsible for its anomalously high luminosity and mass. In this paper, we use an optical spectrum taken in conditions of 0.7-arcsec seeing to provide further evidence that F10214+4724 is a gravitationally lensed system. We measure tentative redshifts of 0.896 and 0.899 for galaxies projected $\approx 1$ and $\approx 3$ arcsec from IRAS F10214+4724 respectively. Identifying the former as the lensing galaxy we obtain a mass:light ratio consistent with those derived for other lenses, and find that its luminosity is consistent with the velocity dispersion deduced from the radius of the Einstein ring. If lensed, our models suggest magnification by a factor $\sim 10$, and hence an intrinsic bolometric luminosity for F10214+4724 similar to local ULIRGs.