The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the partially overlapping BOSS and 2dFLenS surveys. Cross-correlations between galaxy positions and ellipticities are incorporated into the analysis, developing a hybrid model of non-linear scales that blends perturbative and non-perturbative approaches, and assessing signal contributions by astrophysical effects. All weak lensing signals are measured consistently via Fourier-space statistics that are insensitive to the survey mask and display low levels of mode mixing. The calibration of photometric redshift distributions and multiplicative gravitational shear bias is updated, and a more complete tally of residual calibration uncertainties is propagated into the likelihood. A dedicated suite of more than 20000 mocks is used to assess the performance of covariance models and to quantify the impact of survey geometry and spatial variations of survey depth on signals and their errors. The sampling distributions for the likelihood and the $\chi^2$ goodness-of-fit statistic are validated, with proposed changes to calculating the effective number of degrees of freedom. Standard weak lensing point estimates on $S_8=\sigma_8\,(\Omega_{\rm m}/0.3)^{1/2}$ derived from its marginal posterior are easily misinterpreted to be biased low, and an alternative estimator and associated credible interval are proposed. Known systematic effects pertaining to weak lensing modelling and inference are shown to bias $S_8$ by no more than 0.1 standard deviations, with the caveat that no conclusive validation data exist for models of intrinsic galaxy alignments. Compared to the previous KiDS analyses, $S_8$ constraints are expected to improve by 20% for weak lensing alone and by 29% for the joint analysis. [abridged]

We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our `gold-sample' of galaxies, with well calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of $\sim 6.2$ galaxies per square arcminute. We quantify the accuracy of the spatial, temporal and flux-dependent point-spread function (PSF) model, verifying that the model meets our requirements to induce less than a $0.1\sigma$ change in the inferred cosmic shear constraints on the clustering cosmological parameter $S_8 = \sigma_8\sqrt{\Omega_{\rm m}/0.3}$. Through a series of two-point null-tests we validate the shear estimates, finding no evidence for significant non-lensing B-mode distortions in the data. PSF residuals are detected in the highest-redshift bins, originating from object selection and/or weight bias. The amplitude is however shown to be sufficiently low and within our stringent requirements. With a shear-ratio null-test we verify the expected redshift scaling of the galaxy-galaxy lensing signal around luminous red galaxies. We conclude that the joint KiDS-1000 shear and photometric redshift calibration is sufficiently robust for combined-probe gravitational lensing and spectroscopic clustering analyses.

Weak lensing, which is the deflection of light by matter along the line of sight, has proven to be an efficient method for constraining models of structure formation and reveal the nature of dark energy. So far, most weak-lensing studies have focused on the shear field that can be measured directly from the ellipticity of background galaxies. However, within the context of forthcoming full-sky weak-lensing surveys such as Euclid, convergence maps (mass maps) offer an important advantage over shear fields in terms of cosmological exploitation. While it carry the same information, the lensing signal is more compressed in the convergence maps than in the shear field. This simplifies otherwise computationally expensive analyses, for instance, non-Gaussianity studies. However, the inversion of the non-local shear field requires accurate control of systematic effects caused by holes in the data field, field borders, shape noise, and the fact that the shear is not a direct observable (reduced shear). We present the two mass-inversion methods that are included in the official Euclid data-processing pipeline: the standard Kaiser & Squires method (KS), and a new mass-inversion method (KS+) that aims to reduce the information loss during the mass inversion. This new method is based on the KS method and includes corrections for mass-mapping systematic effects. The results of the KS+ method are compared to the original implementation of the KS method in its simplest form, using the Euclid Flagship mock galaxy catalogue. In particular, we estimate the quality of the reconstruction by comparing the two-point correlation functions and third- and fourth- order moments obtained from shear and convergence maps, and we analyse each systematic effect independently and simultaneously...

The nearby TRAPPIST-1 planetary system is an exciting target for characterizing the atmospheres of terrestrial planets. The planets e, f and g lie in the circumstellar habitable zone and could sustain liquid water on their surfaces. During the extended pre-main sequence phase of TRAPPIST-1, however, the planets may have experienced extreme water loss, leading to a desiccated mantle. The presence or absence of an ocean is challenging to determine with current and next generation telescopes. Therefore, we investigate whether indirect evidence of an ocean and/or a biosphere can be inferred from observations of the planetary atmosphere. We introduce a newly developed photochemical model for planetary atmospheres, coupled to a radiative-convective model and validate it against modern Earth, Venus and Mars. The coupled model is applied to the TRAPPIST-1 planets e and f, assuming different surface conditions and varying amounts of CO$_2$ in the atmosphere. As input for the model we use a constructed spectrum of TRAPPIST-1, based on near-simultaneous data from X-ray to optical wavelengths. We compute cloud-free transmission spectra of the planetary atmospheres and determine the detectability of molecular features using the Extremely Large Telescope (ELT) and the James Webb Space Telescope (JWST). We find that under certain conditions, the existence or non-existence of a biosphere and/or an ocean can be inferred by combining 30 transit observations with ELT and JWST within the K-band. A non-detection of CO could suggest the existence of an ocean, whereas significant CH$_4$ hints at the presence of a biosphere.

Context. The density split statistics in weak gravitational lensing analyses probes the correlation between regions of different (foreground) galaxy number densities and their weak lensing signal, measured by the shape distortion of background galaxies. Aims. In this paper, we reconsider density split statistics, by constructing a new angular filter function that is adapted to the expected relation between galaxy number density and shear pattern, in a way that the filter weighting the galaxy number density is matched to the filter that is used to quantify the shear signal. Methods. We use the results of numerical ray-tracing simulations, specifically through Millennium Simulation supplemented by a galaxy distribution based on a semi-analytic model, to construct a matched pair of adapted filter functions for the galaxy density and the tangential shear signal. We compare the performance of our new filter to the previously used top-hat filter, applying both to a different and independent set of numerical simulations (SLICS, cosmo-SLICS). Results. We show that the adapted filter yields a better correlation between the total matter and the galaxy distribution. Furthermore, the adapted filter provides a larger signal-to-noise ratio to constrain the bias between the total matter and the galaxy distribution, and we show that it is, in general, a more sensitive discriminator between different cosmologies, with the exception of cosmologies with very large $\sigma_8$ values. All analyses lead to the conclusion that our adapted filter should be favored in future density split statistic works.

Several semi-analytic models (SAMs) try to explain how galaxies form, evolve and interact inside the dark matter large-scale structure. These SAMs can be tested by comparing their predictions for galaxy-galaxy-galaxy-lensing (G3L), which is weak gravitational lensing around galaxy pairs, with observations. We evaluate the SAMs by Henriques et al. (2015; H15) and by Lagos et al. (2012; L12), implemented in the Millennium Run, by comparing their predictions for G3L to observations at smaller scales than previous studies and also for pairs of lens galaxies from different populations. We compare the G3L signal predicted by the SAMs to measurements in the overlap of the Galaxy And Mass Assembly survey (GAMA), the Kilo-Degree Survey (KiDS), and the VISTA Kilo-degree Infrared Galaxy survey (VIKING), splitting lens galaxies into two colour and five stellar-mass samples. Using an improved G3L estimator, we measure the three-point correlation of the matter distribution for mixed lens pairs with galaxies from different samples, and unmixed lens pairs with galaxies from the same sample. Predictions by the H15 SAM agree with the observations for all colour-selected and all but one stellar-mass-selected sample with 95% confidence. Deviations occur for lenses with stellar masses below $9.5h^{-2}\mathrm{M}_\odot$ at scales below $0.2h^{-1}\mathrm{Mpc}$. Predictions by the L12 SAM for stellar-mass selected samples and red galaxies are significantly higher than observed, while the predicted signal for blue galaxy pairs is too low. The L12 SAM predicts more pairs of small stellar-mass and red galaxies than the H15 SAM and the observations, as well as fewer pairs of blue galaxies. Likely explanations are different treatments of environmental effects by the SAMs and different models of the initial mass function. We conclude that G3L provides a stringent test for models of galaxy formation and evolution.

We revisit the X-ray properties of the main-sequence Hyades members and the relation between X-ray emission and stellar rotation. As input catalog for Hyades members, we combined three recent membership lists derived from Gaia DR2 data including the Hyades core and its tidal tails. We searched for X-ray detections from the ROSAT all-sky survey (RASS) and pointings from ROSAT, Chandra, and XMM-Newton of the Hyades members and adopted rotation periods derived from Kepler's K2 mission and other resources. We find an X-ray detection for 281 of 1066 bona fide main-sequence Hyades members and provide statistical upper limits for the undetected sources. F- and G-type stars have the highest detection fraction (72 %), while K- and M-type dwarfs have lower detection rates (22 %). The X-ray luminosities of the detected members range from ~$2\times10^{27}$ for late M-type dwarfs to ~$2\times10^{30}$ erg s$^{-1}$ for active binaries. The X-ray luminosity distribution functions formally differ for the members in the core and tidal tails, which is likely caused by a larger fraction of field stars in our Hyades tails sample. Compared to previous studies, our sample is slightly fainter in X-rays due to differences in the used Hyades membership list, furthermore, we extend the X-ray luminosity distribution to fainter luminosities. The X-ray activity of F- and G-type stars is well defined at $F_X/F_bol=10^{-5}$ and increases to later spectral types reaching the saturation limit for members later than spectral type M3. The X-ray flux varies by less than a factor of three between epochs for the 104 Hyades members with multiple epoch data. Rotation periods are found for 204 Hyades members, with about half of them being detected in X-rays. The activity-rotation-relation derived for the coeval Hyades members has properties very similar to those obtained by other authors investigating stars of different ages.

Young accreting stars drive outflows that collimate into jets, which can be seen hundreds of au from their driving sources. Accretion and outflow activity cease with system age, and it is believed that magneto-centrifugally launched disk winds are critical agents in regulating accretion through the protoplanetary disk. Protostellar jets are well studied in classical T~Tauri stars ($M_\star\lesssim2\,M_\odot$), while few nearby ($d\lesssim150\,$pc) intermediate-mass stars ($M_\star=2-10\,M_\odot$), known as Herbig Ae/Be stars, have detected jets. We report VLT/MUSE observations of the Herbig~Ae/Be star HD~100546 and the discovery of a protostellar jet. The jet is similar in appearance to jets driven by low-mass stars and compares well with the jet of HD~163296, the only other known optical jet from a nearby Herbig~Ae/Be star. We derive a (one-sided) mass-loss rate in the jet of $\log \dot{M}_{jet} \sim -9.5$ (in $M_\odot$\,yr$^{-1}$) and a ratio of outflow to accretion of roughly $3\times10^{-3}$, which is lower than that of CTTS jets. The discovery of the HD~100546 jet is particularly interesting because the protoplanetary disk around HD~100546 shows a large radial gap, spiral structure, and might host a protoplanetary system. A bar-like structure previously seen in H$\alpha$ with VLT/SPHERE shares the jet position angle, likely represents the base of the jet, and suggests a jet-launching region within about 2\,au. We conclude that the evolution of the disk at radii beyond a few au does not affect the ability of the system to launch jets.

The formation of low-mass stars in molecular clouds involves accretion disks and jets, which are of broad astrophysical interest. Accreting stars represent the closest examples of these phenomena. Star and planet formation are also intimately connected, setting the starting point for planetary systems like our own. The ultraviolet (UV) spectral range is particularly suited to study star formation, because virtually all relevant processes radiate at temperatures associated with UV emission processes or have strong observational signatures in the UV. In this review, we describe how UV observations provide unique diagnostics for the accretion process, the physical properties of the protoplanetary disk, and jets and outflows.

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Context. Future weak lensing surveys, such as the Euclid mission, will attempt to measure the shapes of billions of galaxies in order to derive cosmological information. These surveys will attain very low levels of statistical error, and systematic errors must be extremely well controlled. In particular, the point spread function (PSF) must be estimated using stars in the field, and recovered with high accuracy. Aims. The aims of this paper are twofold. Firstly, we took steps toward a nonparametric method to address the issue of recovering the PSF field, namely that of finding the correct PSF at the position of any galaxy in the field, applicable to Euclid. Our approach relies solely on the data, as opposed to parametric methods that make use of our knowledge of the instrument. Secondly, we studied the impact of imperfect PSF models on the shape measurement of galaxies themselves, and whether common assumptions about this impact hold true in an Euclid scenario. Methods. We extended the recently proposed resolved components analysis approach, which performs super-resolution on a field of under-sampled observations of a spatially varying, image-valued function. We added a spatial interpolation component to the method, making it a true 2-dimensional PSF model. We compared our approach to PSFEx, then quantified the impact of PSF recovery errors on galaxy shape measurements through image simulations. Results. Our approach yields an improvement over PSFEx in terms of the PSF model and on observed galaxy shape errors, though it is at present far from reaching the required Euclid accuracy. We also find that the usual formalism used for the propagation of PSF model errors to weak lensing quantities no longer holds in the case of an Euclid-like PSF. In particular, different shape measurement approaches can react differently to the same PSF modeling errors.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

Young stars accrete mass and angular momentum from their circumstellar disks. Some of them also drive outflows, which can be distinguished in optical forbidden emission lines (FELs). We analyze a sample of binary T Tauri stars observed with long-slit spectroscopy by the Hubble Space Telescope (HST), searching for spatially resolved outflows. We detect resolved [O I] emission in two cases out of twenty one. In DF Tau we resolve high and medium velocity outflows in a jet and counterjet out to 60 au. The outflows are accelerated within the inner 12 au and retain a constant speed thereafter. In UY Aur, we detect a blue- and a red-shifted outflow from UY Aur A, as well as a blue-shifted jet from UY Aur B. All of these features have been seen in [Fe II] with data taken ten years apart indicating that the underlying outflow pattern is stable on these timescales.

In 2015 we started the XMM-Newton monitoring of the young solar-like star Epsilon Eridani (440 Myr), one of the youngest solar-like stars with a known chromospheric CaII cycle. By analyzing the most recent Mount Wilson S-index CaII data of this star, we found that the chromospheric cycle lasts 2.92 +/- 0.02 yr, in agreement with past results. From the long-term X-ray lightcurve, we find clear and systematic X-ray variability of our target, consistent with the chromospheric CaII cycle. The average X-ray luminosity results to be 2 x 10^28 erg/s, with an amplitude that is only a factor 2 throughout the cycle. We apply a new method to describe the evolution of the coronal emission measure distribution of Epsilon Eridani in terms of solar magnetic structures: active regions, cores of active regions and flares covering the stellar surface at varying filling fractions. Combinations of these magnetic structures can describe the observed X-ray emission measure of Epsilon Eridani only if the solar flare emission measure distribution is restricted to events in the decay phase. The interpretation is that flares in the corona of Epsilon Eridani last longer than their solar counterparts. We ascribe this to the lower metallicity of Epsilon Eridani. Our analysis revealed also that the X-ray cycle of Epsilon Eridani is strongly dominated by cores of active regions. The coverage fraction of cores throughout the cycle changes by the same factor as the X-ray luminosity. The maxima of the cycle are characterized by a high percentage of covering fraction of the flares, consistent with the fact that flaring events are seen in the corresponding short-term X-ray lightcurves predominately at the cycle maxima. The high X-ray emission throughout the cycle of Epsilon Eridani is thus explained by the high percentage of magnetic structures on its surface.

The UV radiation field is a critical regulator of gas-phase chemistry in surface layers of disks around young stars. In an effort to understand the relationship between photocatalyzing UV radiation fields and gas emission observed at infrared and sub-mm wavelengths, we present an analysis of new and archival HST, Spitzer, ALMA, IRAM, and SMA data for five targets in the Lupus cloud complex and 14 systems in Taurus-Auriga. The HST spectra were used to measure LyA and FUV continuum fluxes reaching the disk surface, which are responsible for dissociating relevant molecular species (e.g. HCN, N2). Semi-forbidden C II] 2325 and UV-fluorescent H2 emission were also measured to constrain inner disk populations of C+ and vibrationally excited H2. We find a significant positive correlation between 14 micron HCN emission and fluxes from the FUV continuum and C II] 2325, consistent with model predictions requiring N2 photodissociation and carbon ionization to trigger the main CN/HCN formation pathways. We also report significant negative correlations between sub-mm CN emission and both C II] and FUV continuum fluxes, implying that CN is also more readily dissociated in disks with stronger FUV irradiation. No clear relationships are detected between either CN or HCN and LyA or UV-H2 emission. This is attributed to the spatial stratification of the various molecular species, which span several vertical layers and radii across the inner and outer disk. We expect that future observations with JWST will build on this work by enabling more sensitive IR surveys than were possible with Spitzer.

We present cosmological constraints from a joint cosmic shear analysis of the Kilo-Degree Survey (KV450) and the Dark Energy Survey (DES-Y1), conducted using Complete Orthogonal Sets of E/B-Integrals (COSEBIs). With COSEBIs we isolate any B-modes which have a non-cosmic shear origin and demonstrate the robustness of our cosmological E-mode analysis as no significant B-modes are detected. We highlight how COSEBIs are fairly insensitive to the amplitude of the non-linear matter power spectrum at high $k$-scales, mitigating the uncertain impact of baryon feedback in our analysis. COSEBIs, therefore, allow us to utilise additional small-scale information, improving the DES-Y1 joint constraints on $S_8=\sigma_8(\Omega_{\rm m}/0.3)^{0.5}$ and $\Omega_{\rm m}$ by $20\%$. Adopting a flat $\Lambda$CDM model we find $S_8=0.755^{+0.019}_{-0.021}$, which is in $3.2\sigma$ tension with the Planck Legacy analysis of the cosmic microwave background.

We present a semi-analytic model for the shear two-point correlation function of a cosmic shear survey with non-uniform depth. Ground-based surveys are subject to depth variations that primarily arise through varying atmospheric conditions. For a survey like the Kilo-Degree Survey (KiDS), we find that the measured depth variation increases the amplitude of the observed shear correlation function at the level of a few percent out to degree-scales, relative to the assumed uniform-depth case. The impact on the inferred cosmological parameters is shown to be insignificant for a KiDS-like survey. For next-generation cosmic shear experiments, however, we conclude that variable depth should be accounted for.

Galaxy-galaxy-galaxy lensing (G3L) is a powerful tool for constraining the three-point correlation between the galaxy and matter distribution and thereby models of galaxy evolution. We propose three improvements to current measurements of G3L: (i) a weighting of lens galaxies according to their redshift difference, (ii) adaptive binning of the three-point correlation function, and (iii) accounting for the effect of lens magnification by the cosmic large-scale structure. Improvement (i) is designed to improve the precision of the G3L measurement, whereas improvements (ii) and (iii) remove biases of the estimator. We further show how the G3L signal can be converted from angular into physical scales. The improvements were tested on simple mock data and simulated data based on the Millennium Run with an implemented semi-analytic galaxy model. Our improvements increase the signal-to-noise ratio by 35 % on average at angular scales between 0.1 arcmin and 10 arcmin and physical scales between $0.02$ and $2 \, h^{-1}\, \textrm{Mpc}$. They also remove the bias of the G3L estimator at angular scales below 1 arcmin, which was originally up to 40 %. The signal due to lens magnification is approximately 10 % of the total signal.

We present a weak lensing detection of filamentary structures in the cosmic web, combining data from the Kilo-Degree Survey, the Red Cluster Sequence Lensing Survey and the Canada-France-Hawaii Telescope Lensing Survey. The line connecting luminous red galaxies with a separation of $3 - 5\, h^{-1}\text{Mpc}$ is chosen as a proxy for the location of filaments. We measure the average weak lensing shear around $\sim$11,000 candidate filaments selected in this way from the Sloan Digital Sky Survey. After nulling the shear induced by the dark matter haloes around each galaxy, we report a $3.4\,\sigma$ detection of an anisotropic shear signal from the matter that connects them. Adopting a filament density profile, motivated from $N$-body simulations, the average density at the centre of these filamentary structures is found to be $15 \pm 4$ times the critical density.

We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey (VIKING). This is the first time that a full optical to near-infrared data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning $450~$deg$^2$, allows us to improve significantly the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and - most importantly - solidify our knowledge of the redshift distributions of the sources. Based on a flat $\Lambda$CDM model we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.737_{-0.036}^{+0.040}$ in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with $S_8$ differing by $2.3\sigma$. This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with Planck is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Magnification changes the observed number counts of galaxies on the sky. This biases the observed tangential shear profiles around galaxies, the so-called galaxy-galaxy lensing (GGL) signal, and the related excess mass profile. Correspondingly, inference of physical quantities, such as the mean mass profile of halos around galaxies, are affected by magnification effects. We use simulated shear and galaxy data of the Millennium Simulation to quantify the effect on shear and mass estimates from magnified lens and source number counts. The former are due to the large-scale matter distribution in the foreground of the lenses, the latter are caused by magnification of the source population by the matter associated with the lenses. The GGL signal is calculated from the simulations by an efficient fast-Fourier transform that can also be applied to real data. The numerical treatment is complemented by a leading-order analytical description of the magnification effects, which is shown to fit the numerical shear data well. We find the magnification effect is strongest for steep galaxy luminosity functions and high redshifts. For a lens redshift of $z_\mathrm{d}=0.83$, a limiting magnitude of $22\,\mathrm{mag}$ in the $r$-band and a source redshift of $z_\mathrm{s}=0.99$, we find that a magnification correction changes the shear profile up to $45\%$ and the mass is biased by up to $55 \%$. For medium-redshift galaxies the relative change in shear and mass is typically a few percent. As expected, the sign of the bias depends on the local slope of the lens luminosity function $\alpha_\mathrm{d}$, where the mass is biased low for $\alpha_\mathrm{d}<1$ and biased high for $\alpha_\mathrm{d}>1$. Whereas the magnification effect of sources is rarely than more $1\%$, the statistical power of future weak lensing surveys warrants correction for this effect.

Atmospheric escape is key to explaining the long-term evolution of planets in our Solar System and beyond, and in the interpretation of atmospheric measurements. Hydrodynamic escape is generally thought to be driven by the flux of extreme ultraviolet photons that the planet receives from its host star. Here, we show that the escape from planets orbiting hot stars proceeds through a different yet complementary process: drawing its energy from the intense near ultraviolet emission of the star that is deposited within an optically thin, high-altitude atmospheric layer of hydrogen excited into the lower state of the Balmer series. The ultra-hot exoplanet KELT-9b likely represents the first known instance of this Balmer-driven escape. In this regime of hydrodynamic escape, the near ultraviolet emission from the star is more important at determining the planet mass loss than the extreme ultraviolet emission, and uncertainties in the latter become less critical. Further, we predict that gas exoplanets around hot stars may experience catastrophic mass loss when they are less massive than 1--2 Jupiter masses and closer-in than KELT-9b, thereby challenging the paradigm that all large exoplanets are stable to atmospheric escape. We argue that extreme escape will affect the demographics of close-in exoplanets orbiting hot stars.

We present the curation and verification of a new combined optical and near infrared dataset for cosmology and astrophysics, derived from the combination of $ugri$-band imaging from the Kilo Degree Survey (KiDS) and $ZY\!J\!H\!K_{\rm s}$-band imaging from the VISTA Kilo degree Infrared Galaxy (VIKING) survey. This dataset is unrivaled in cosmological imaging surveys due to its combination of area ($458$ deg$^2$ before masking), depth ($r\le25$), and wavelength coverage ($ugriZY\!J\!H\!K_{\rm s}$). The combination of survey depth, area, and (most importantly) wavelength coverage allows significant reductions in systematic uncertainties (i.e. reductions of between 10 and 60\% in bias, outlier rate, and scatter) in photometric-to-spectroscopic redshift comparisons, compared to the optical-only case at photo-$z$ above $0.7$. The complementarity between our optical and NIR surveys means that over $80\%$ of our sources, across all photo-$z$, have significant detections (i.e. not upper limits) in our $8$ reddest bands. We derive photometry, photo-$z$, and stellar masses for all sources in the survey, and verify these data products against existing spectroscopic galaxy samples. We demonstrate the fidelity of our higher-level data products by constructing the survey stellar mass functions in 8 volume-complete redshift bins. We find that these photometrically derived mass functions provide excellent agreement with previous mass evolution studies derived using spectroscopic surveys. The primary data products presented in this paper are publicly available at http://kids.strw.leidenuniv.nl/.

CVSO 30 is a young, active, weak-line T Tauri star; it possibly hosts the only known planetary system with both a transiting hot-Jupiter and a cold-Jupiter candidate (CVSO 30 b and c). We analyzed archival ROSAT, Chandra, and XMM-Newton data to study the coronal emission in the system. According to our modeling, CVSO 30 shows a quiescent X-ray luminosity of about 8e29 erg/s. The X-ray absorbing column is consistent with interstellar absorption. XMM-Newton observed a flare, during which a transit of the candidate CVSO 30 b was expected, but no significant transit-induced variation in the X-ray flux is detectable. While the hot-Jupiter candidate CVSO 30 b has continuously been undergoing mass loss powered by the high-energy irradiation, we conclude that its evaporation lifetime is considerably longer than the estimated stellar age of 2.6 Myr.

In modern weak-lensing surveys, the common approach to correct for residual systematic biases in the shear is to calibrate shape measurement algorithms using simulations. These simulations must fully capture the complexity of the observations to avoid introducing any additional bias. In this paper we study the importance of faint galaxies below the observational detection limit of a survey. We simulate simplified Euclid VIS images including and excluding this faint population, and measure the shift in the multiplicative shear bias between the two sets of simulations. We measure the shear with three different algorithms: a moment-based approach, model fitting, and machine learning. We find that for all methods, a spatially uniform random distribution of faint galaxies introduces a shear multiplicative bias of the order of a few times $10^{-3}$. This value increases to the order of $10^{-2}$ when including the clustering of the faint galaxies, as measured in the Hubble Space Telescope Ultra-Deep Field. The magnification of the faint background galaxies due to the brighter galaxies along the line of sight is found to have a negligible impact on the multiplicative bias. We conclude that the undetected galaxies must be included in the calibration simulations with proper clustering properties down to magnitude 28 in order to reach a residual uncertainty on the multiplicative shear bias calibration of a few times $10^{-4}$, in line with the $2\times10^{-3}$ total accuracy budget required by the scientific objectives of the Euclid survey. We propose two complementary methods for including faint galaxy clustering in the calibration simulations.

The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope, specifically designed for measuring weak gravitational lensing by galaxies and large-scale structure. When completed it will consist of 1350 square degrees imaged in four filters (ugri). Here we present the fourth public data release which more than doubles the area of sky covered by data release 3. We also include aperture-matched ZYJHKs photometry from our partner VIKING survey on the VISTA telescope in the photometry catalogue. We illustrate the data quality and describe the catalogue content. Two dedicated pipelines are used for the production of the optical data. The Astro-WISE information system is used for the production of co-added images in the four survey bands, while a separate reduction of the r-band images using the theli pipeline is used to provide a source catalogue suitable for the core weak lensing science case. All data have been re-reduced for this data release using the latest versions of the pipelines. The VIKING photometry is obtained as forced photometry on the theli sources, using a re-reduction of the VIKING data that starts from the VISTA pawprints. Modifications to the pipelines with respect to earlier releases are described in detail. The photometry is calibrated to the Gaia DR2 G band using stellar locus regression. In this data release a total of 1006 square-degree survey tiles with stacked ugri images are made available, accompanied by weight maps, masks, and single-band source lists. We also provide a multi-band catalogue based on r-band detections, including homogenized photometry and photometric redshifts, for the whole dataset. Mean limiting magnitudes (5 sigma in a 2" aperture) are 24.23, 25.12, 25.02, 23.68 in ugri, respectively, and the mean r-band seeing is 0.70".

We analyse three public cosmic shear surveys; the Kilo-Degree Survey (KiDS-450), the Dark Energy Survey (DES-SV) and the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Adopting the COSEBIs statistic to cleanly and completely separate the lensing E-modes from the non-lensing B-modes, we detect B-modes in KiDS-450 and CFHTLenS at the level of about 2.7 $\sigma$. For DES- SV we detect B-modes at the level of 2.8 $\sigma$ in a non-tomographic analysis, increasing to a 5.5 $\sigma$ B-mode detection in a tomographic analysis. In order to understand the origin of these detected B-modes we measure the B-mode signature of a range of different simulated systematics including PSF leakage, random but correlated PSF modelling errors, camera-based additive shear bias and photometric redshift selection bias. We show that any correlation between photometric-noise and the relative orientation of the galaxy to the point-spread-function leads to an ellipticity selection bias in tomographic analyses. This work therefore introduces a new systematic for future lensing surveys to consider. We find that the B-modes in DES-SV appear similar to a superposition of the B-mode signatures from all of the systematics simulated. The KiDS-450 and CFHTLenS B-mode measurements show features that are consistent with a repeating additive shear bias.

Young stars accrete mass from a circumstellar disk, but at the same time disk and star eject outflows and jets. These outflows have an onion-like structure where the innermost and fastest layers are surrounded by increasingly lower velocity components. The outer layers are probably photo-evaporative and magnetocentrifugally launched disk winds, but the nature of the inner winds is still uncertain. Since the fastest components carry only a small fraction of the mass, they are best observed at high-energies (X-ray and UV) as the slower, more massive components do not reach plasma temperatures sufficient for relevant X-ray or UV emission. Outflows are the most likely way in which a star or its disk can shed angular momentum and allow accretion to proceed; thus we cannot understand the accretion and the rotation rate of young stars if we cannot solve the origin of the inner jet components. Stellar jets share characteristics with their counterparts in more massive astrophysical objects, such as stellar mass black holes and AGN, with the added benefit that young stars are found at much closer distances and thus scales not accessible in other types of objects can be resolved. To understand the origin and impact of the inner jets, sub-arcsecond imaging and spectroscopy in the UV and X-rays is required, together with theory and modelling to interpret existing and future observations.

All cool stars show magnetic activity, and X-ray emission is the hallmark of this activity. Gaining an understanding of activity aids us in answering fundamental questions about stellar astrophysics and in determining the impact of activity on the exoplanets that orbit these stars. Stellar activity is driven by magnetic fields, which are ultimately powered by convection and stellar rotation. However, the resulting dynamo properties heavily depend on the stellar interior structure and are far from being understood. X-ray radiation can evaporate exoplanet atmospheres and damage organic materials on the planetary surface, reducing the probability that life can form or be sustained, but also provides an important source of energy for prebiotic chemical reactions. Over the next two years, the TESS mission will deliver a catalog of the closest exoplanets, along with rotation periods and activity diagnostics for millions of stars, whether or not they have a planet. We propose to include all cool stars that are TESS targets and bright enough for Chandra observations, as determined by their detection and flux in the ROSAT all-sky survey (RASS), to the list of Chandra Cool Attitude Targets (CATs). For each target, the signal will be sufficient to fit the coronal plasma with at least two temperature components, and compare abundances of groups of elements with low, medium, or high first ionization potential. Similar to the known relation between X-ray luminosity $L_X$ and rotation period, we can correlate stellar properties with coronal temperatures and abundances to constrain models for stellar activity, coronal heating, and stellar dynamos. Detailed X-ray characterization for even a subset of planet-hosting systems would dramatically advance our knowledge of what impact these emissions have on orbiting planets.

We study projected underdensities in the cosmic galaxy density field known as 'troughs', and their overdense counterparts, which we call 'ridges'. We identify these regions using a bright sample of foreground galaxies from the photometric Kilo-Degree Survey (KiDS), specifically selected to mimic the spectroscopic Galaxy And Mass Assembly survey (GAMA). Using background galaxies from KiDS, we measure the weak gravitational lensing profiles of the troughs/ridges. We quantify the amplitude of their lensing strength $A$ as a function of galaxy density percentile rank $P$ and galaxy overdensity $\delta$, and find that the skewness in the galaxy density distribution is reflected in the total mass distribution measured by weak lensing. We interpret our results using the mock galaxy catalogue from the Marenostrum Institut de Ci\`encies de l'Espai (MICE) simulation, and find a good agreement with our observations. Using signal-to-noise weights derived from the Scinet LIghtCone Simulations (SLICS) mock catalogue we optimally stack the lensing signal of KiDS troughs with an angular radius $\theta_A$ = {5,10,15,20} arcmin, resulting in {16.8,14.9,10.13,7.55} $\sigma$ detections. Finally, we select troughs using a volume-limited sample of galaxies, split into two redshift bins between 0.1 < z < 0.3. For troughs/ridges with transverse comoving radius $R_A$ = 1.9 Mpc/h, we find no significant difference in the comoving Excess Surface Density as a function of $P$ and $\delta$ between the low- and high-redshift sample. Using the MICE and SLICS mocks we predict that trough and ridge evolution could be detected with gravitational lensing using deeper and wider lensing surveys, such as those from the Large Synoptic Survey Telescope and Euclid.

Close-in gas planets are subject to continuous photoevaporation that can erode their volatile envelopes. Today, ongoing mass loss has been confirmed in a few individual systems via transit observations in the ultraviolet spectral range. We demonstrate that the Ultraviolet/Optical Telescope (UVOT) onboard the Neil Gehrels Swift Observatory enables photometry to a relative accuracy of about 0.5% and present the first near-UV (200-270 nm, NUV) transit observations of WASP-121 b, a hot Jupiter with one of the highest predicted mass-loss rates. The data cover the orbital phases 0.85 to 1.15 with three visits. We measure a broad-band NUV transit depth of $2.10\pm0.29$%. While still consistent with the optical value of 1.55%, the NUV data indicate excess absorption of 0.55% at a 1.9$\sigma$ level. Such excess absorption is known from the WASP-12 system, and both of these hot Jupiters are expected to undergo mass loss at extremely high rates. With a CLOUDY simulation, we show that absorption lines of Fe II in a dense extended atmosphere can cause broad-band NUV absorption at the 0.5% level. Given the numerous lines of low-ionization metals, the NUV range is a promising tracer of photoevaporation in the hottest gas planets.

As part of the Mega MUSCLES Hubble Space Telescope (HST) Treasury program, we obtained time-series ultraviolet spectroscopy of the M2.5V star, GJ~674. During the FUV monitoring observations, the target exhibited several small flares and one large flare (E_FUV = 10^{30.75} ergs) that persisted over the entirety of a HST orbit and had an equivalent duration >30,000 sec, comparable to the highest relative amplitude event previously recorded in the FUV. The flare spectrum exhibited enhanced line emission from chromospheric, transition region, and coronal transitions and a blue FUV continuum with an unprecedented color temperature of T_c ~ 40,000+/-10,000 K. In this paper, we compare the flare FUV continuum emission with parameterizations of radiative hydrodynamic model atmospheres of M star flares. We find that the observed flare continuum can be reproduced using flare models but only with the ad hoc addition of hot, dense emitting component. This observation demonstrates that flares with hot FUV continuum temperatures and significant EUV/FUV energy deposition will continue to be of importance to exoplanet atmospheric chemistry and heating even as the host M dwarfs age beyond their most active evolutionary phases.

Using the same lens galaxies, the ratios of tangential shears for different source galaxy redshifts is equal to the ratios of their corresponding angular-diameter distances. This is the so-called shear-ratio test (SRT) and it is valid when effects induced by the intervening large-scale structure (LSS) can be neglected. The dominant LSS effect is magnification bias which, on the one hand, induces an additional shear, and on the other hand, causes a magnification of the lens population. Our objective is to quantify the magnification bias for the SRT and show an easy-to-apply mitigation strategy that does not rely on additional observations. We use ray-tracing data through the Millennium simulation to measure the influence of magnification on the SRT and test our mitigation strategy. Using the SRT as a null-test we find deviations from zero up to $10 \%$ for a flux-limited sample of lens galaxies, which is a strong function of lens redshift and the lens-source line-of-sight separation. Using our mitigation strategy we can improve the null-test by a factor of $\sim \!100$.

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.

Hot gas giant exoplanets can lose part of their atmosphere due to strong stellar irradiation, affecting their physical and chemical evolution. Studies of atmospheric escape from exoplanets have mostly relied on space-based observations of the hydrogen Lyman-{\alpha} line in the far ultraviolet which is strongly affected by interstellar absorption. Using ground-based high-resolution spectroscopy we detect excess absorption in the helium triplet at 1083 nm during the transit of the Saturn-mass exoplanet WASP-69b, at a signal-to-noise ratio of 18. We measure line blue shifts of several km/s and post transit absorption, which we interpret as the escape of part of the atmosphere trailing behind the planet in comet-like form. [Additional notes by authors: Furthermore, we provide upper limits for helium signals in the atmospheres of the exoplanets HD 209458b, KELT-9b, and GJ 436b. We investigate the host stars of all planets with detected helium signals and those of the three planets we derive upper limits for. In each case we calculate the X-ray and extreme ultraviolet flux received by these planets. We find that helium is detected in the atmospheres of planets (orbiting the more active stars and) receiving the larger amount of irradiation from their host stars.]

We present three transit observations of HD 189733 b obtained with the high-resolution spectrograph CARMENES at Calar Alto. A strong absorption signal is detected in the near-infrared He I triplet at 10830 \AA{} in all three transits. During mid-transit, the mean absorption level is $0.88\pm0.04$ % measured in a $\pm$10 km s$^{-1}$ range at a net blueshift of $-3.5\pm0.4$ km s$^{-1}$ (10829.84--10830.57 \AA{}). The absorption signal exhibits radial velocities of $+6.5\pm3.1$ km s$^{-1}$ and $-12.6\pm1.0$ km s$^{-1}$ during ingress and egress, respectively; measured in the planetary rest frame. We show that stellar activity related pseudo-signals interfere with the planetary atmospheric absorption signal. They could contribute as much as 80% of the observed signal and might also affect the radial velocity signature, but pseudo-signals are very unlikely to explain the entire signal. The observed line ratio between the two unresolved and the third line of the He I triplet is $2.8\pm0.2$, which strongly deviates from the value expected for an optically thin atmospheres. When interpreted in terms of absorption in the planetary atmosphere, this favors a compact helium atmosphere with an extent of only 0.2 planetary radii and a substantial column density on the order of $4\times 10^{12}$ cm$^{-2}$. The observed radial velocities can be understood either in terms of atmospheric circulation with equatorial superrotation or as a sign of an asymmetric atmospheric component of evaporating material. We detect no clear signature of ongoing evaporation, like pre- or post-transit absorption, which could indicate material beyond the planetary Roche lobe, or radial velocities in excess of the escape velocity. These findings do not contradict planetary evaporation, but only show that the detected helium absorption in HD 189733 b does not trace the atmospheric layers that show pronounced escape signatures.

The matter environment of galaxies is connected to the physics of galaxy formation and evolution. Utilising galaxy-galaxy-galaxy lensing as a direct probe, we map out the distribution of correlated surface mass-density around galaxy pairs for different lens separations in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We compare, for the first time, these so-called excess mass maps to predictions provided by a recent semi-analytic model, which is implanted within the dark-matter Millennium Simulation. We analyse galaxies with stellar masses between $10^9-10^{11}\,{\rm M}_\odot$ in two photometric redshift bins, for lens redshifts $z\lesssim0.6$, focusing on pairs inside groups and clusters. To allow us a better interpretation of the maps, we discuss the impact of chance pairs, i.e., galaxy pairs that appear close to each other in projection only. Our tests with synthetic data demonstrate that the patterns observed in the maps are essentially produced by correlated pairs that are close in redshift ($\Delta z\lesssim5\times10^{-3}$). We also verify the excellent accuracy of the map estimators. In an application to the galaxy samples in the CFHTLenS, we obtain a $3\sigma-6\sigma$ significant detection of the excess mass and an overall good agreement with the galaxy model predictions. There are, however, a few localised spots in the maps where the observational data disagrees with the model predictions on a $\approx3.5\sigma$ confidence level. Although we have no strong indications for systematic errors in the maps, this disagreement may be related to the residual B-mode pattern observed in the average of all maps. Alternatively, misaligned galaxy pairs inside dark matter halos or lensing by a misaligned distribution of the intra-cluster gas might also cause the unanticipated bulge in the distribution of the excess mass between lens pairs.

We present a weak gravitational lensing analysis of 815 square degree of $i$-band imaging from the Kilo-Degree Survey (KiDS-$i$-800). In contrast to the deep $r$-band observations, which take priority during excellent seeing conditions and form the primary KiDS dataset (KiDS-$r$-450), the complementary yet shallower KiDS-$i$-800 spans a wide range of observing conditions. The overlapping KiDS-$i$-800 and KiDS-$r$-450 imaging therefore provides a unique opportunity to assess the robustness of weak lensing measurements. In our analysis, we introduce two new `null' tests. The `nulled' two-point shear correlation function uses a matched catalogue to show that the calibrated KiDS-$i$-800 and KiDS-$r$-450 shear measurements agree at the level of $1 \pm 4$\%. We use five galaxy lens samples to determine a `nulled' galaxy-galaxy lensing signal from the full KiDS-$i$-800 and KiDS-$r$-450 surveys and find that the measurements agree to $7 \pm 5$\% when the KiDS-$i$-800 source redshift distribution is calibrated using either spectroscopic redshifts, or the 30-band photometric redshifts from the COSMOS survey.

(Abridged) The detection of forming planets in disks around young stars remains elusive, and state-of-the-art observational techniques provide somewhat ambiguous results. It has been reported that the pre-transitional T Tauri star LkCa 15 could host three planets; candidate planet b is in the process of formation, as inferred from its H$\alpha$ emission. However, a more recent work casts doubts on the planetary nature of the previous detections. We have observed LkCa 15 with ISIS/WHT. The spectrograph's slit was oriented towards the last reported position of LkCa 15 b (parallel direction) and 90degr from that (perpendicular). The photocenter and full width half maximum (FWHM) of the Gaussians fitting the spatial distribution at H$\alpha$ and the adjacent continuum were measured. A well-known binary (GU CMa) was used as a calibrator to test the spectro-astrometric performance of ISIS/WHT, recovering consistent photocenter and FWHM signals. However, the photocenter shift predicted for LkCa 15 b is not detected, but the FWHM in H$\alpha$ is broader than in the continuum for both slit positions. Our simulations show that the photocenter and FWHM observations cannot be explained simultaneously by an accreting planet. In turn, both spectro-astrometric observations are naturally reproduced from a roughly symmetric Halpha emitting region centered on the star and extent comparable to the orbit originally attributed to the planet at several au. The extended H$\alpha$ emission around LkCa 15 could be related to a variable disk wind, but additional multi-epoch data and detailed modeling are necessary to understand its physical nature. Spectro-astrometry in H$\alpha$ is able to test the presence of accreting planets and can be used as a complementary technique to survey planet formation in circumstellar disks.

We present a survey of far-ultraviolet (FUV; 1150 - 1450 Ang) emission line spectra from 71 planet-hosting and 33 non-planet-hosting F, G, K, and M dwarfs with the goals of characterizing their range of FUV activity levels, calibrating the FUV activity level to the 90 - 360 Ang extreme-ultraviolet (EUV) stellar flux, and investigating the potential for FUV emission lines to probe star-planet interactions (SPIs). We build this emission line sample from a combination of new and archival observations with the Hubble Space Telescope-COS and -STIS instruments, targeting the chromospheric and transition region emission lines of Si III, N V, C II, and Si IV. We find that the exoplanet host stars, on average, display factors of 5 - 10 lower UV activity levels compared with the non-planet hosting sample; this is explained by a combination of observational and astrophysical biases in the selection of stars for radial-velocity planet searches. We demonstrate that UV activity-rotation relation in the full F - M star sample is characterized by a power-law decline (with index $\alpha$ ~ -1.1), starting at rotation periods >~3.5 days. Using N V or Si IV spectra and a knowledge of the star's bolometric flux, we present a new analytic relationship to estimate the intrinsic stellar EUV irradiance in the 90 - 360 Ang band with an accuracy of roughly a factor of ~2. Finally, we study the correlation between SPI strength and UV activity in the context of a principal component analysis that controls for the sample biases. We find that SPIs are not a statistically significant contributor to the observed UV activity levels.

We present the first "clipped" cosmic shear measurement using data from the Kilo-Degree Survey (KiDS-450). "Clipping" transformations suppress the signal from the highest density, non-linear regions of cosmological fields. We demonstrate that these transformations improve constraints on $S_8=\sigma_8(\Omega_{\rm{m}}/0.3)^{0.5}$ when used in combination with conventional two-point statistics. For the KiDS-450 data, we find that the combined measurements improve the constraints on $S_8$ by 17%, compared to shear correlation functions alone. We determine the expectation value of the clipped shear correlation function using a suite of numerical simulations, and develop methodology to mitigate the impact of masking and shot noise. Future improvements in numerical simulations and mass reconstruction methodology will permit the precise calibration of clipped cosmic shear statistics such that clipping can become a standard tool in weak lensing analyses.

RW Aur is a binary system composed of two young, low-mass stars. The primary, RW Aur A, has undergone visual dimming events ($\Delta V =2-3$~mag) in 2011, 2014-16, and 2017-2018. Visual and IR observations indicate a gray absorber that moved into the line-of-sight. This dimming is also associated with changes in the outflow. In 2017, when the optical brightness was almost 2~mag below the long-term average we triggered a Chandra observation to measure the absorbing column density $N_\mathrm{H}$ and to constrain dust properties and the gas-to-dust ratio of the absorber. In 2017, the X-ray spectrum is more absorbed than it was in the optically bright state ($N_\mathrm{H} = (4\pm 1) \times 10^{23}\;\mathrm{cm}^{-2}$) and shows significantly more hot plasma than in X-ray observations taken before. Also, a new emission feature at $6.63\pm0.02$ keV (statistic) $\pm0.02$ keV (systematic) appeared indicating an Fe abundance an order of magnitude above Solar, in contrast with previous sub-Solar Fe abundance measurements. Comparing X-ray absorbing column density $N_\mathrm{H}$ and optical extinction $A_V$, we find that either the gas-to-dust ratio in the absorber is orders of magnitude higher than in the ISM or the absorber has undergone significant dust evolution. Given the high column density coupled with changes in the X-ray spectral shape, this absorber is probably located in the inner disk. We speculate that a break-up of planetesimals or a terrestrial planet could supply large grains causing gray absorption; some of these grains would be accreted and enrich the stellar corona with iron which could explain the inferred high abundance.

Classical T Tauri stars (CTTSs) accrete matter from the inner edge of their surrounding circumstellar disks. The impact of the accretion material on the stellar atmosphere results in a strong shock, which causes emission from the X-ray to the near-infrared (NIR) domain. Shock velocities of several 100 km s$^{-1}$ imply that the immediate post shock plasma emits mainly in X-rays. Indeed, two X-ray diagnostics, the so-called soft excess and the high densities observed in He-like triplets, differentiate CTTSs from their non-accreting siblings. However, accretion shock properties derived from X-ray diagnostics often contradict established ultraviolet (UV)-NIR accretion tracers and a physical model simultaneously explaining both, X-ray and UV-NIR accretion tracers, is not yet available. We present new XMM-Newton and Chandra grating observations of the CTTS T Tauri combined with UV and optical data. During all epochs, the soft excess is large and the densities derived from the O VII and Ne IX He-like triplets are compatible with coronal densities. This confirms that the soft X-ray emission cannot originate in accretion funnels that carry the bulk of the accretion rate despite T Tauri's large soft excess. Instead, we propose a model of radially density stratified accretion columns to explain the density diagnostics and the soft excess. In addition, accretion rate and X-ray luminosity are inversely correlated in T Tauri over several epochs. Such an anti-correlation has been observed in samples of stars. Hence the process causing it must be intrinsic to the accretion process, and we speculate that the stellar magnetic field configuration on the visibile hemisphere affects both the accretion rate and the coronal emission, eventually causing the observed anti-correlation.

Although timing variations in close binary systems have been studied for a long time, their underlying causes are still unclear. A possible explanation is the so-called Applegate mechanism, where a strong, variable magnetic field can periodically change the gravitational quadrupole moment of a stellar component, thus causing observable period changes. One of the systems exhibiting such strong orbital variations is the RS CVn binary HR 1099, whose activity cycle has been studied by various authors via photospheric and chromospheric activity indicators, resulting in contradicting periods. We aim at independently determining the magnetic activity cycle of HR 1099 using archival X-ray data to allow for a comparison to orbital period variations. Archival X-ray data from 80 different observations of HR 1099 acquired with 12 different X-ray facilities and covering almost four decades were used to determine X-ray fluxes in the energy range of 2-10 keV via spectral fitting and flux conversion. Via the Lomb-Scargle periodogram we analyze the resulting long-term X-ray light curve to search for periodicities. We do not detect any statistically significant periodicities within the X-ray data. An analysis of optical data of HR 1099 shows that the derivation of such periods is strongly dependent on the time coverage of available data, since the observed optical variations strongly deviate from a pure sine wave. We argue that this offers an explanation as to why other authors derive such a wide range of activity cycle periods based on optical data. We conclude that our analysis constitutes the longest stellar X-ray activity light curve acquired to date, yet the still rather sparse sampling of the X-ray data, along with stochastic flaring activity, does not allow for the independent determination of an X-ray activity cycle.

We measure the projected galaxy clustering and galaxy-galaxy lensing signals using the Galaxy And Mass Assembly (GAMA) survey and Kilo-Degree Survey (KiDS) to study galaxy bias. We use the concept of non-linear and stochastic galaxy biasing in the framework of halo occupation statistics to constrain the parameters of the halo occupation statistics and to unveil the origin of galaxy biasing. The bias function $\Gamma_{\text{gm}}(r_{\text{p}})$, where $r_{\text{p}}$ is the projected comoving separation, is evaluated using the analytical halo model from which the scale dependence of $\Gamma_{\text{gm}}(r_{\text{p}})$, and the origin of the non-linearity and stochasticity in halo occupation models can be inferred. Our observations unveil the physical reason for the non-linearity and stochasticity, further explored using hydrodynamical simulations, with the stochasticity mostly originating from the non-Poissonian behaviour of satellite galaxies in the dark matter haloes and their spatial distribution, which does not follow the spatial distribution of dark matter in the halo. The observed non-linearity is mostly due to the presence of the central galaxies, as was noted from previous theoretical work on the same topic. We also see that overall, more massive galaxies reveal a stronger scale dependence, and out to a larger radius. Our results show that a wealth of information about galaxy bias is hidden in halo occupation models. These models should therefore be used to determine the influence of galaxy bias in cosmological studies.

Planets that reside close-in to their host star are subject to intense high-energy irradiation. Extreme-ultraviolet (EUV) and X-ray radiation (together, XUV) is thought to drive mass loss from planets with volatile envelopes. We present $\textit{XMM-Newton}$ observations of six nearby stars hosting transiting planets in tight orbits (with orbital period, $P_\text{orb} < 10\,$d), wherein we characterise the XUV emission from the stars and subsequent irradiation levels at the planets. In order to reconstruct the unobservable EUV emission, we derive a new set of relations from Solar $\textit{TIMED/SEE}$ data that are applicable to the standard bands of the current generation of X-ray instruments. From our sample, WASP-80b and HD$\,$149026b experience the highest irradiation level, but HAT-P-11b is probably the best candidate for Ly$\,\alpha$ evaporation investigations because of the system's proximity to the Solar System. The four smallest planets have likely lost a greater percentage of their mass over their lives than their larger counterparts. We also detect the transit of WASP-80b in the near ultraviolet with the Optical Monitor on $\textit{XMM-Newton}$.

We present cosmological parameter constraints from a joint analysis of three cosmological probes: the tomographic cosmic shear signal in $\sim$450 deg$^2$ of data from the Kilo Degree Survey (KiDS), the galaxy-matter cross-correlation signal of galaxies from the Galaxies And Mass Assembly (GAMA) survey determined with KiDS weak lensing, and the angular correlation function of the same GAMA galaxies. We use fast power spectrum estimators that are based on simple integrals over the real-space correlation functions, and show that they are practically unbiased over relevant angular frequency ranges. We test our full pipeline on numerical simulations that are tailored to KiDS and retrieve the input cosmology. By fitting different combinations of power spectra, we demonstrate that the three probes are internally consistent. For all probes combined, we obtain $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}=0.800_{-0.027}^{+0.029}$, consistent with Planck and the fiducial KiDS-450 cosmic shear correlation function results. Marginalising over wide priors on the mean of the tomographic redshift distributions yields consistent results for $S_8$ with an increase of $28\%$ in the error. The combination of probes results in a $26\%$ reduction in uncertainties of $S_8$ over using the cosmic shear power spectra alone. The main gain from these additional probes comes through their constraining power on nuisance parameters, such as the galaxy intrinsic alignment amplitude or potential shifts in the redshift distributions, which are up to a factor of two better constrained compared to using cosmic shear alone, demonstrating the value of large-scale structure probe combination.

We compute the spherical-sky weak-lensing power spectrum of the shear and convergence. We discuss various approximations, such as flat-sky, and first- and second- order Limber equations for the projection. We find that the impact of adopting these approximations is negligible when constraining cosmological parameters from current weak lensing surveys. This is demonstrated using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We find that the reported tension with Planck Cosmic Microwave Background (CMB) temperature anisotropy results cannot be alleviated. For future large-scale surveys with unprecedented precision, we show that the spherical second-order Limber approximation will provide sufficient accuracy. In this case, the cosmic-shear power spectrum is shown to be in agreement with the full projection at the sub-percent level for l > 3, with the corresponding errors an order of magnitude below cosmic variance for all l. When computing the two-point shear correlation function, we show that the flat-sky fast Hankel transformation results in errors below two percent compared to the full spherical transformation. In the spirit of reproducible research, our numerical implementation of all approximations and the full projection are publicly available within the package nicaea at http://www.cosmostat.org/software/nicaea.

Disk warps around classical T Tauri stars (CTTS) can periodically obscure the central star for some viewing geometries. For these so-called AA Tau-like variables, the obscuring material is located in the inner disk and absorption spectroscopy allows one to characterize its dust and gas content. Since the observed emission from CTTS consists of several components (photospheric, accretion, jet, and disk emission), which can all vary with time, it is generally challenging to disentangling disk features from emission variability. Multi-epoch, flux-calibrated, broadband spectra provide us with the necessary information to cleanly separate absorption from emission variability. We applied this method to three epochs of VLT/X-Shooter spectra of the CTTS V354 Mon (CSI Mon-660) located in NGC 2264 and find that: (a) the accretion emission remains virtually unchanged between the three epochs; (b) the broadband flux evolution is best described by disk material obscuring part of the star, and (c) the Na and K gas absorption lines show only a minor increase in equivalent width during phases of high dust extinction. The limits on the absorbing gas column densities indicate a low gas-to-dust ratio in the inner disk, less than a tenth of the ISM value. We speculate that the evolutionary state of V354 Mon, rather old with a low accretion rate, is responsible for the dust excess through an evolution toward a dust dominated disk or through the fragmentation of larger bodies that drifted inward from larger radii in a still gas dominated disk.

Faint X-ray emission from hot plasma (T > 1 MK) has been detected extending outward a few arcseconds along the optically-delineated jets of some classical T Tauri stars including RY Tau. The mechanism and location where the jet is heated to X-ray temperatures is unknown. We present high spatial resolution HST far-ultraviolet long-slit observations of RY Tau with the slit aligned along the jet. The primary objective was to search for C IV emission from warm plasma at T(CIV) ~ 10$^{5}$ K within the inner jet (<1 arcsec) that cannot be fully-resolved by X-ray telescopes. Spatially-resolved C IV emission is detected in the blueshifted jet extending outward from the star to 1 arcsec and in the redshifted jet out to 0.5 arcsec. C IV line centroid shifts give a radial velocity in the blueshifted jet of -136 $\pm$ 10 km/s at an offset of 0.29 arcsec (39 au) and deceleration outward is detected. The deprojected jet speed is subject to uncertainties in the jet inclination but values >200 km/s are likely. The mass-loss rate in the blueshifted jet is at least 2.3 $\times$ 10$^{-9}$ M_sun yr$^{-1}$, consistent with optical determinations. We use the HST data along with optically-determined jet morphology to place meaningful constraints on candidate jet-heating models including a hot-launch model in which the jet is heated near the base to X-ray temperatures by an unspecified (but probably magnetic) process, and downstream heating from shocks or a putative jet magnetic field.

We present a study of molecular gas in the inner disk $\left(r < 20 \, \text{AU} \right)$ around RY Lupi, with spectra from HST-COS, HST-STIS, and VLT-CRIRES. We model the radial distribution of flux from hot gas in a surface layer between $r = 0.1-10$ AU, as traced by Ly$\alpha$-pumped H$_2$. The result shows H$_2$ emission originating in a ring centered at $\sim$3 AU that declines within $r < 0.1$ AU, which is consistent with the behavior of disks with dust cavities. An analysis of the H$_2$ line shapes shows that a two-component Gaussian profile $\left(\text{FWHM}_{broad, H_2} = 105 \pm 15 \, \text{km s}^{-1}; \, \text{FWHM}_{narrow, H_2} = 43 \pm 13 \, \text{km s}^{-1} \right)$ is statistically preferred to a single-component Gaussian. We interpret this as tentative evidence for gas emitting from radially separated disk regions $\left( \left \langle r_{broad, H_2} \right \rangle \sim 0.4 \pm 0.1 \, \text{AU}; \, \left \langle r_{narrow, H_2} \right \rangle \sim 3 \pm 2 \, \text{AU} \right)$. The 4.7 $\mu$m $^{12}$CO emission lines are also well fit by two-component profiles $\left( \left \langle r_{broad, CO} \right \rangle = 0.4 \pm 0.1 \, \text{AU}; \, \left \langle r_{narrow, CO} \right \rangle = 15 \pm 2 \, \text{AU} \right)$. We combine these results with 10 $\mu$m observations to form a picture of gapped structure within the mm-imaged dust cavity, providing the first such overview of the inner regions of a young disk. The HST SED of RY Lupi is available online for use in modeling efforts.

The central ambition of the modern time delay cosmography consists in determining the Hubble constant $H_0$ with a competitive precision. However, the tension with $H_0$ obtained from the Planck satellite for a spatially-flat $\Lambda$CDM cosmology suggests that systematic errors may have been underestimated. The most critical one probably comes from the degeneracy existing between lens models that was first formalized by the well-known mass-sheet transformation (MST). In this paper, we assess to what extent the source position transformation (SPT), a more general invariance transformation which contains the MST as a special case, may affect the time delays predicted by a model. To this aim we use pySPT, a new open-source python package fully dedicated to the SPT that we present in a companion paper. For axisymmetric lenses, we find that the time delay ratios between a model and its SPT-modified counterpart simply scale like the corresponding source position ratios, $\Delta \hat{t}/ \Delta t \approx \hat{\beta}/\beta$, regardless of the mass profile and the isotropic SPT. Similar behavior (almost) holds for non-axisymmetric lenses in the double image regime and for opposite image pairs in the quadruple image regime. In the latter regime, we also confirm that the time delay ratios are not conserved. In addition to the MST effects, the SPT-modified time delays deviate in general no more than a few percent for particular image pairs, suggesting that its impact on time-delay cosmography seems not be as crucial as initially suspected. We also reflected upon the relevance of the SPT validity criterion and present arguments suggesting that it should be reconsidered. Even though a new validity criterion would affect the time delays in a different way, we expect from numerical simulations that our conclusions will remain unchanged.

We perform a combined analysis of cosmic shear tomography, galaxy-galaxy lensing tomography, and redshift-space multipole power spectra (monopole and quadrupole) using 450 deg$^2$ of imaging data by the Kilo Degree Survey (KiDS) overlapping with two spectroscopic surveys: the 2-degree Field Lensing Survey (2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS). We restrict the galaxy-galaxy lensing and multipole power spectrum measurements to the overlapping regions with KiDS, and self-consistently compute the full covariance between the different observables using a large suite of $N$-body simulations. We methodically analyze different combinations of the observables, finding that galaxy-galaxy lensing measurements are particularly useful in improving the constraint on the intrinsic alignment amplitude (by 30%, positive at $3.5\sigma$ in the fiducial data analysis), while the multipole power spectra are useful in tightening the constraints along the lensing degeneracy direction (e.g. factor of two stronger matter density constraint in the fiducial analysis). The fully combined constraint on $S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.742 \pm 0.035$, which is an improvement by 20% compared to KiDS alone, corresponds to a $2.6\sigma$ discordance with Planck, and is not significantly affected by fitting to a more conservative set of scales. Given the tightening of the parameter space, we are unable to resolve the discordance with an extended cosmology that is simultaneously favored in a model selection sense, including the sum of neutrino masses, curvature, evolving dark energy, and modified gravity. The complementarity of our observables allows for constraints on modified gravity degrees of freedom that are not simultaneously bounded with either probe alone, and up to a factor of three improvement in the $S_8$ constraint in the extended cosmology compared to KiDS alone.

This paper is the first of a series of papers constraining cosmological parameters with weak lensing peak statistics using $\sim 450~\rm deg^2$ of imaging data from the Kilo Degree Survey (KiDS-450). We measure high signal-to-noise ratio (SNR: $\nu$) weak lensing convergence peaks in the range of $3<\nu<5$, and employ theoretical models to derive expected values. These models are validated using a suite of simulations. We take into account two major systematic effects, the boost factor and the effect of baryons on the mass-concentration relation of dark matter haloes. In addition, we investigate the impacts of other potential astrophysical systematics including the projection effects of large scale structures, intrinsic galaxy alignments, as well as residual measurement uncertainties in the shear and redshift calibration. Assuming a flat $\Lambda$CDM model, we find constraints for $S_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.5}=0.746^{+0.046}_{-0.107}$ according to the degeneracy direction of the cosmic shear analysis and $\Sigma_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.38}=0.696^{+0.048}_{-0.050}$ based on the derived degeneracy direction of our high-SNR peak statistics. The difference between the power index of $S_{\rm 8}$ and in $\Sigma_{\rm 8}$ indicates that combining the two probes has the potential to break the degeneracy in $\sigma_{\rm 8}$ and $\Omega_{\rm m}$. Our results are consistent with the cosmic shear tomographic correlation analysis of the same dataset and $\sim 2\sigma$ lower than the Planck 2016 results.

We study the statistics of peaks in a weak lensing reconstructed mass map of the first 450 square degrees of the Kilo Degree Survey. The map is computed with aperture masses directly applied to the shear field with an NFW-like compensated filter. We compare the peak statistics in the observations with that of simulations for various cosmologies to constrain the cosmological parameter $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$, which probes the ($\Omega_{\rm m}, \sigma_8$) plane perpendicularly to its main degeneracy. We estimate $S_8=0.750\pm0.059$, using peaks in the signal-to-noise range $0 \leq {\rm S/N} \leq 4$, and accounting for various systematics, such as multiplicative shear bias, mean redshift bias, baryon feedback, intrinsic alignment, and shear-position coupling. These constraints are $\sim25\%$ tighter than the constraints from the high significance peaks alone ($3 \leq {\rm S/N} \leq 4$) which typically trace single-massive halos. This demonstrates the gain of information from low-S/N peaks. However we find that including ${\rm S/N} < 0$ peaks does not add further information. Our results are in good agreement with the tomographic shear two-point correlation function measurement in KiDS-450. Combining shear peaks with non-tomographic measurements of the shear two-point correlation functions yields a $\sim20\%$ improvement in the uncertainty on $S_8$ compared to the shear two-point correlation functions alone, highlighting the great potential of peaks as a cosmological probe.

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.

HD 209458 is one of the benchmark objects in the study of hot Jupiter atmospheres and their evaporation through planetary winds. The expansion of the planetary atmosphere is thought to be driven by high-energy EUV and X-ray irradiation. We obtained new Chandra HRC-I data, which unequivocally show that HD 209458 is an X-ray source. Combining these data with archival XMM-Newton observations, we find that the corona of HD 209458 is characterized by a temperature of about 1 MK and an emission measure of 7e49 cm^-3, yielding an X-ray luminosity of 1.6e27 erg/s in the 0.124-2.48 keV band. HD 209458 is an inactive star with a coronal temperature comparable to that of the inactive Sun but a larger emission measure. At this level of activity, the planetary high-energy emission is sufficient to support mass-loss at a rate of a few times 1e10 g/s.

We present measurements of the weak gravitational lensing shear power spectrum based on $450$ sq. deg. of imaging data from the Kilo Degree Survey. We employ a quadratic estimator in two and three redshift bins and extract band powers of redshift auto-correlation and cross-correlation spectra in the multipole range $76 \leq \ell \leq 1310$. The cosmological interpretation of the measured shear power spectra is performed in a Bayesian framework assuming a $\Lambda$CDM model with spatially flat geometry, while accounting for small residual uncertainties in the shear calibration and redshift distributions as well as marginalising over intrinsic alignments, baryon feedback and an excess-noise power model. Moreover, massive neutrinos are included in the modelling. The cosmological main result is expressed in terms of the parameter combination $S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ yielding $S_8 = \ 0.651 \pm 0.058$ (3 z-bins), confirming the recently reported tension in this parameter with constraints from Planck at $3.2\sigma$ (3 z-bins). We cross-check the results of the 3 z-bin analysis with the weaker constraints from the 2 z-bin analysis and find them to be consistent. The high-level data products of this analysis, such as the band power measurements, covariance matrices, redshift distributions, and likelihood evaluation chains are available at http://kids.strw.leidenuniv.nl/

The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope. It aims to image 1500 square degrees in four filters (ugri). The core science driver is mapping the large-scale matter distribution in the Universe, using weak lensing shear and photometric redshift measurements. Further science cases include galaxy evolution, Milky Way structure, detection of high-redshift clusters, and finding rare sources such as strong lenses and quasars. Here we present the third public data release (DR3) and several associated data products, adding further area, homogenized photometric calibration, photometric redshifts and weak lensing shear measurements to the first two releases. A dedicated pipeline embedded in the Astro-WISE information system is used for the production of the main release. Modifications with respect to earlier releases are described in detail. Photometric redshifts have been derived using both Bayesian template fitting, and machine-learning techniques. For the weak lensing measurements, optimized procedures based on the THELI data reduction and lensfit shear measurement packages are used. In DR3 stacked ugri images, weight maps, masks, and source lists for 292 new survey tiles (~300 sq.deg) are made available. The multi-band catalogue, including homogenized photometry and photometric redshifts, covers the combined DR1, DR2 and DR3 footprint of 440 survey tiles (447 sq.deg). Limiting magnitudes are typically 24.3, 25.1, 24.9, 23.8 (5 sigma in a 2 arcsec aperture) in ugri, respectively, and the typical r-band PSF size is less than 0.7 arcsec. The photometric homogenization scheme ensures accurate colors and an absolute calibration stable to ~2% for gri and ~3% in u. Separately released are a weak lensing shear catalogue and photometric redshifts based on two different machine-learning techniques.

In weak gravitational lensing, weighted quadrupole moments of the brightness profile in galaxy images are a common way to estimate gravitational shear. We employ general adaptive moments (GLAM) to study causes of shear bias on a fundamental level and for a practical definition of an image ellipticity. The GLAM ellipticity has useful properties for any chosen weight profile: the weighted ellipticity is identical to that of isophotes of elliptical images, and in absence of noise and pixellation it is always an unbiased estimator of reduced shear. We show that moment-based techniques, adaptive or unweighted, are similar to a model-based approach in the sense that they can be seen as imperfect fit of an elliptical profile to the image. Due to residuals in the fit, moment-based estimates of ellipticities are prone to underfitting bias when inferred from observed images. The estimation is fundamentally limited mainly by pixellation which destroys information on the original, pre-seeing image. We give an optimized estimator for the pre-seeing GLAM ellipticity and quantify its bias for noise-free images. To deal with pixel noise, we consider a Bayesian approach where the posterior of the GLAM ellipticity can be inconsistent with the true ellipticity if we do not properly account for our ignorance about fit residuals. This underfitting bias is S/N-independent but changes with the pre-seeing brightness profile and the correlation or heterogeneity of pixel noise over the post-seeing image. Furthermore, when inferring a constant ellipticity or, more relevantly, constant shear from a source sample with a distribution of intrinsic properties (sizes, centroid positions, intrinsic shapes), an additional, now noise-dependent bias arises towards low S/N if incorrect priors for the intrinsic properties are used. We discuss the origin of this prior bias.

Characterizing the UV spectral energy distribution (SED) of an exoplanet host star is critically important for assessing its planet's potential habitability, particularly for M dwarfs as they are prime targets for current and near-term exoplanet characterization efforts and atmospheric models predict that their UV radiation can produce photochemistry on habitable zone planets different than on Earth. To derive ground-based proxies for UV emission for use when Hubble Space Telescope observations are unavailable, we have assembled a sample of fifteen early-to-mid M dwarfs observed by Hubble, and compared their non-simultaneous UV and optical spectra. We find that the equivalent width of the chromospheric Ca II K line at 3933 Angstroms, when corrected for spectral type, can be used to estimate the stellar surface flux in ultraviolet emission lines, including H I Lyman alpha. In addition, we address another potential driver of habitability: energetic particle fluxes associated with flares. We present a new technique for estimating soft X-ray and >10 MeV proton flux during far-UV emission line flares (Si IV and He II) by assuming solar-like energy partitions. We analyze several flares from the M4 dwarf GJ 876 observed with Hubble and Chandra as part of the MUSCLES Treasury Survey and find that habitable zone planets orbiting GJ 876 are impacted by large Carrington-like flares with peak soft X-ray fluxes >1e-3 W m-2 and possible proton fluxes ~100-1000 pfu, approximately four orders of magnitude more frequently than modern-day Earth.

We test extensions to the standard cosmological model with weak gravitational lensing tomography using 450 deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). In these extended cosmologies, which include massive neutrinos, nonzero curvature, evolving dark energy, modified gravity, and running of the scalar spectral index, we also examine the discordance between KiDS and cosmic microwave background measurements from Planck. The discordance between the two datasets is largely unaffected by a more conservative treatment of the lensing systematics and the removal of angular scales most sensitive to nonlinear physics. The only extended cosmology that simultaneously alleviates the discordance with Planck and is at least moderately favored by the data includes evolving dark energy with a time-dependent equation of state (in the form of the $w_0-w_a$ parameterization). In this model, the respective $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ constraints agree at the $1\sigma$ level, and there is `substantial concordance' between the KiDS and Planck datasets when accounting for the full parameter space. Moreover, the Planck constraint on the Hubble constant is wider than in LCDM and in agreement with the Riess et al. (2016) direct measurement of $H_0$. The dark energy model is moderately favored as compared to LCDM when combining the KiDS and Planck measurements, and remains moderately favored after including an informative prior on the Hubble constant. In both of these scenarios, marginalized constraints in the $w_0-w_a$ plane are discrepant with a cosmological constant at the $3\sigma$ level. Moreover, KiDS constrains the sum of neutrino masses to 4.0 eV (95% CL), finds no preference for time or scale dependent modifications to the metric potentials, and is consistent with flatness and no running of the spectral index. The analysis code is public at https://github.com/sjoudaki/kids450

We study the stellar-to-halo mass relation of central galaxies in the range 9.7<log_10(M_*/h^-2 M_sun)<11.7 and z<0.4, obtained from a combined analysis of the Kilo Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. We use ~100 deg^2 of KiDS data to study the lensing signal around galaxies for which spectroscopic redshifts and stellar masses were determined by GAMA. We show that lensing alone results in poor constraints on the stellar-to-halo mass relation due to a degeneracy between the satellite fraction and the halo mass, which is lifted when we simultaneously fit the stellar mass function. At M_sun>5x10^10 h^-2 M_sun, the stellar mass increases with halo mass as ~M_h^0.25. The ratio of dark matter to stellar mass has a minimum at a halo mass of 8x10^11 h^-1 M_sun with a value of M_h/M_*=56_-10^+16 [h]. We also use the GAMA group catalogue to select centrals and satellites in groups with five or more members, which trace regions in space where the local matter density is higher than average, and determine for the first time the stellar-to-halo mass relation in these denser environments. We find no significant differences compared to the relation from the full sample, which suggests that the stellar-to-halo mass relation does not vary strongly with local density. Furthermore, we find that the stellar-to-halo mass relation of central galaxies can also be obtained by modelling the lensing signal and stellar mass function of satellite galaxies only, which shows that the assumptions to model the satellite contribution in the halo model do not significantly bias the stellar-to-halo mass relation. Finally, we show that the combination of weak lensing with the stellar mass function can be used to test the purity of group catalogues.

We constrain the average halo ellipticity of ~2 600 galaxy groups from the Galaxy And Mass Assembly (GAMA) survey, using the weak gravitational lensing signal measured from the overlapping Kilo Degree Survey (KiDS). To do so, we quantify the azimuthal dependence of the stacked lensing signal around seven different proxies for the orientation of the dark matter distribution, as it is a priori unknown which one traces the orientation best. On small scales, the major axis of the brightest group/cluster member (BCG) provides the best proxy, leading to a clear detection of an anisotropic signal. In order to relate that to a halo ellipticity, we have to adopt a model density profile. We derive new expressions for the quadrupole moments of the shear field given an elliptical model surface mass density profile. Modeling the signal with an elliptical Navarro-Frenk-White (NFW) profile on scales < 250 kpc, which roughly corresponds to half the virial radius, and assuming that the BCG is perfectly aligned with the dark matter, we find an average halo ellipticity of e_h=0.38 +/- 0.12. This agrees well with results from cold-dark-matter-only simulations, which typically report values of e_h ~ 0.3. On larger scales, the lensing signal around the BCGs does not trace the dark matter distribution well, and the distribution of group satellites provides a better proxy for the halo's orientation instead, leading to a 3--4 sigma detection of a non-zero halo ellipticity at scales between 250 kpc and 750 kpc. Our results suggest that the distribution of stars enclosed within a certain radius forms a good proxy for the orientation of the dark matter within that radius, which has also been observed in hydrodynamical simulations.

We present new tests to identify stationary position-dependent additive shear biases in weak gravitational lensing data sets. These tests are important diagnostics for currently ongoing and planned cosmic shear surveys, as such biases induce coherent shear patterns that can mimic and potentially bias the cosmic shear signal. The central idea of these tests is to determine the average ellipticity of all galaxies with shape measurements in a grid in the pixel plane. The distribution of the absolute values of these averaged ellipticities can be compared to randomized catalogues; a difference points to systematics in the data. In addition, we introduce a method to quantify the spatial correlation of the additive bias, which suppresses the contribution from cosmic shear and therefore eases the identification of a position-dependent additive shear bias in the data. We apply these tests to the publicly available shear catalogues from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) and the Kilo Degree Survey (KiDS) and find evidence for a small but non-negligible residual additive bias at small scales. As this residual bias is smaller than the error on the shear correlation signal at those scales, it is highly unlikely that it causes a significant bias in the published cosmic shear results of CFHTLenS. In CFHTLenS, the amplitude of this systematic signal is consistent with zero in fields where the number of stars used to model the PSF is higher than average, suggesting that the position-dependent additive shear bias originates from undersampled PSF variations across the image.

Early-type galaxies provide unique tests for the predictions of the cold dark matter cosmology and the baryonic physics assumptions entering models for galaxy formation. In this work, we use the Illustris simulation to study correlations of three main properties of early-type galaxies, namely, the stellar orbital anisotropies, the central dark matter fractions and the central radial density slopes, as well as their redshift evolution since $z=1.0$. We find that lower-mass galaxies or galaxies at higher redshift tend to be bluer in rest-frame colour, have higher central gas fractions, and feature more tangentially anisotropic orbits and steeper central density slopes than their higher-mass or lower-redshift counterparts, respectively. The projected central dark matter fraction within the effective radius shows a very mild mass dependence but positively correlates with galaxy effective radii due to the aperture effect. The central density slopes obtained by combining strong lensing measurements with single aperture kinematics are found to differ from the true density slopes. We identify systematic biases in this measurement to be due to two common modelling assumptions, isotropic stellar orbital distributions and power-law density profiles. We also compare the properties of early-type galaxies in Illustris to those from existing galaxy and strong lensing surveys, we find in general broad agreement but also some tension, which poses a potential challenge to the stellar formation and feedback models adopted by the simulation.

A power-law density model, i.e., $\rho(r) \propto r^{-\gamma'}$ has been commonly employed in strong gravitational lensing studies, including the so-called time-delay technique used to infer the Hubble constant $H_0$. However, since the radial scale at which strong lensing features are formed corresponds to the transition from the dominance of baryonic matter to dark matter, there is no known reason why galaxies should follow a power law in density. The assumption of a power law artificially breaks the mass-sheet degeneracy, a well-known invariance transformation in gravitational lensing which affects the product of Hubble constant and time delay and can therefore cause a bias in the determination of $H_0$ from the time-delay technique. In this paper, we use the Illustris hydrodynamical simulations to estimate the amplitude of this bias, and to understand how it is related to observational properties of galaxies. Investigating a large sample of Illustris galaxies that have velocity dispersion $\sigma_{SIE}$>160 km/s at redshifts below $z=1$, we find that the bias on $H_0$ introduced by the power-law assumption can reach 20%-50%, with a scatter of $10\%-30\%$ (rms). However, we find that by selecting galaxies with an inferred power-law model slope close to isothermal, it is possible to reduce the bias on $H_0$ to <5%, and the scatter to <10%. This could potentially be used to form less biased statistical samples for $H_0$ measurements in the upcoming large survey era.

We investigate possible signatures of halo assembly bias for spectroscopically selected galaxy groups from the GAMA survey using weak lensing measurements from the spatially overlapping regions of the deeper, high-imaging-quality photometric KiDS survey. We use GAMA groups with an apparent richness larger than 4 to identify samples with comparable mean host halo masses but with a different radial distribution of satellite galaxies, which is a proxy for the formation time of the haloes. We measure the weak lensing signal for groups with a steeper than average and with a shallower than average satellite distribution and find no sign of halo assembly bias, with the bias ratio of $0.85^{+0.37}_{-0.25}$, which is consistent with the $\Lambda$CDM prediction. Our galaxy groups have typical masses of $10^{13} M_{\odot}/h$, naturally complementing previous studies of halo assembly bias on galaxy cluster scales.

We present the tomographic cross-correlation between galaxy lensing measured in the Kilo Degree Survey (KiDS-450) with overlapping lensing measurements of the cosmic microwave background (CMB), as detected by Planck 2015. We compare our joint probe measurement to the theoretical expectation for a flat $\Lambda$CDM cosmology, assuming the best-fitting cosmological parameters from the KiDS-450 cosmic shear and Planck CMB analyses. We find that our results are consistent within $1\sigma$ with the KiDS-450 cosmology, with an amplitude re-scaling parameter $A_{\rm KiDS} = 0.86 \pm 0.19$. Adopting a Planck cosmology, we find our results are consistent within $2\sigma$, with $A_{\it Planck} = 0.68 \pm 0.15$. We show that the agreement is improved in both cases when the contamination to the signal by intrinsic galaxy alignments is accounted for, increasing $A$ by $\sim 0.1$. This is the first tomographic analysis of the galaxy lensing -- CMB lensing cross-correlation signal, and is based on five photometric redshift bins. We use this measurement as an independent validation of the multiplicative shear calibration and of the calibrated source redshift distribution at high redshifts. We find that constraints on these two quantities are strongly correlated when obtained from this technique, which should therefore not be considered as a stand-alone competitive calibration tool.

We study intrinsic alignments (IA) of galaxy image shapes within the Illustris cosmic structure formation simulations. We investigate how IA correlations depend on observable galaxy properties such as stellar mass, apparent magnitude, redshift, and photometric type, and on the employed shape measurement method. The correlations considered include the matter density-intrinsic ellipticity (mI), galaxy density-intrinsic ellipticity (dI), gravitational shear-intrinsic ellipticity (GI), and intrinsic ellipticity-intrinsic ellipticity (II) correlations. We find stronger correlations for more massive and more luminous galaxies, as well as for earlier photometric types, in agreement with observations. Moreover, the correlations significantly depend on the choice of shape estimator, even if calibrated to serve as unbiased shear estimators. In particular, shape estimators that down-weight the outer parts of galaxy images produce much weaker IA signals on intermediate and large scales than methods employing flat radial weights. The expected contribution of intrinsic alignments to the observed ellipticity correlation in tomographic cosmic shear surveys may be below one percent or several percent of the full signal depending on the details of the shape measurement method. A comparison of our results to a tidal alignment model indicates that such a model is able to reproduce the IA correlations well on intermediate and large scales, provided the effect of varying galaxy density is correctly taken into account. We also find that the GI contributions to the observed ellipticity correlations could be inferred directly from measurements of galaxy density-intrinsic ellipticity correlations, except on small scales, where systematic differences between mI and dI correlations are large.

Testing predictions of semi-analytic models of galaxy evolution against observations help to understand the complex processes that shape galaxies. We compare predictions from the Garching and Durham models implemented on the Millennium Run with observations of galaxy-galaxy lensing (GGL) and galaxy-galaxy-galaxy lensing (G3L) for various galaxy samples with stellar masses in the range 0.5 < (M_* / 10^10 M_Sun) < 32 and photometric redshift range 0.2 < z < 0.6 in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We find that the predicted GGL and G3L signals are in qualitative agreement with CFHTLenS data. Quantitatively, the models succeed in reproducing the observed signals in the highest stellar mass bin (16 < ( M_* / 10^10 M_Sun) < 32) but show different degrees of tension for the other stellar mass samples. The Durham models are strongly excluded at the 95% confidence level by the observations as they largely over-predict the amplitudes of the GGL and G3L signals, probably because they predict too many satellite galaxies in massive halos.

We present The-wiZZ, an open source and user-friendly software for estimating the redshift distributions of photometric galaxies with unknown redshifts by spatially cross-correlating them against a reference sample with known redshifts. The main benefit of The-wiZZ is in separating the angular pair finding and correlation estimation from the computation of the output clustering redshifts allowing anyone to create a clustering redshift for their sample without the intervention of an "expert". It allows the end user of a given survey to select any sub-sample of photometric galaxies with unknown redshifts, match this sample's catalog indices into a value-added data file, and produce a clustering redshift estimation for this sample in a fraction of the time it would take to run all the angular correlations needed to produce a clustering redshift. We show results with this software using photometric data from the Kilo-Degree Survey (KiDS) and spectroscopic redshifts from the Galaxy and Mass Assembly (GAMA) survey and the Sloan Digital Sky Survey (SDSS). The results we present for KiDS are consistent with the redshift distributions used in a recent cosmic shear analysis from the survey. We also present results using a hybrid machine learning-clustering redshift analysis that enables the estimation of clustering redshifts for individual galaxies. The-wiZZ can be downloaded at http://github.com/morriscb/The-wiZZ/.

We present a re-analysis of the CFHTLenS weak gravitational lensing survey using Complete Orthogonal Sets of E/B-mode Integrals, known as COSEBIs. COSEBIs provide a complete set of functions to efficiently separate E-modes from B-modes and hence allow for robust and stringent tests for systematic errors in the data. This analysis reveals significant B-modes on large angular scales that were not previously seen using the standard E/B decomposition analyses. We find that the significance of the B-modes is enhanced when the data is split by galaxy type and analysed in tomographic redshift bins. Adding tomographic bins to the analysis increases the number of COSEBIs modes, which results in a less accurate estimation of the covariance matrix from a set of simulations. We therefore also present the first compressed COSEBIs analysis of survey data, where the COSEBIs modes are optimally combined based on their sensitivity to cosmological parameters. In this tomographic CCOSEBIs analysis we find the B-modes to be consistent with zero when the full range of angular scales are considered.

We present predictions for the galaxy-galaxy lensing profile from the EAGLE hydrodynamical cosmological simulation at redshift z=0.18, in the spatial range 0.02 < R/(Mpc/h) < 2, and for five logarithmically equi-spaced stellar mass bins in the range 10.3 < $\log_{10}$(Mstar/ $M_{\odot}$) < 11.8. We compare these excess surface density profiles to the observed signal from background galaxies imaged by the Kilo Degree Survey around spectroscopically confirmed foreground galaxies from the GAMA survey. Exploiting the GAMA galaxy group catalogue, the profiles of central and satellite galaxies are computed separately for groups with at least five members to minimise contamination. EAGLE predictions are in broad agreement with the observed profiles for both central and satellite galaxies, although the signal is underestimated at R$\approx$0.5-2 Mpc/h for the highest stellar mass bins. When central and satellite galaxies are considered simultaneously, agreement is found only when the selection function of lens galaxies is taken into account in detail. Specifically, in the case of GAMA galaxies, it is crucial to account for the variation of the fraction of satellite galaxies in bins of stellar mass induced by the flux-limited nature of the survey. We report the inferred stellar-to-halo mass relation and we find good agreement with recent published results. We note how the precision of the galaxy-galaxy lensing profiles in the simulation holds the potential to constrain fine-grained aspects of the galaxy-dark matter connection.

RW Aur is a young binary system showing strong signatures of a recent tidal encounter between the circumprimary disk and the secondary star. The primary star has recently undergone two major dimming events ($\Delta$mag $\approx$ 2 in V-band), whose origin is still under debate. To shed light on the mechanism leading to the dimming events, we study the extinction properties, accretion variability, and gas kinematics using absorption lines from the material obscuring star RW Aur A. We compare our moderate resolution X-Shooter spectra of the dim state of RW Aur A with other spectral observations. In particular, we analyse archival high resolution UVES spectra obtained during the bright state of the system, in order to track the evolution of the spectral properties across the second dimming event. The spectrum obtained during the dim state shows narrow absorption lines in the Na and K optical doublets, where the former is saturated. With a velocity of -60 km/s these lines indicate that during the dim state the disk wind is either enhanced, or significantly displaced into the line of sight. The photometric evolution across the dimming event shows a gray extinction, and is correlated with a significant reduction of the EW of all photospheric lines. Emission lines tracing accretion do not vary significantly across the dimming. We conclude that the dimming event is related to a major perturbation on the inner disk. We suggest that the inner disk is occulting (most of) the star, and thus its photosphere, but is not occulting the accretion regions within a few stellar radii. Since observations of the outer disk indicate that the disk is modestly inclined (45 - 60 deg), we propose that the inner disk might be warped by a yet unseen (sub-)stellar companion, which may also explain the 2.77 day periodic variability of the spectral lines.

We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). For a flat $\Lambda$CDM cosmology with a prior on $H_0$ that encompasses the most recent direct measurements, we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.745\pm0.039$. This result is in good agreement with other low redshift probes of large scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A $2.3$-$\sigma$ tension in $S_8$ and `substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved `self-calibrating' version of $lens$fit validated using an extensive suite of image simulations. Four-band $ugri$ photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov Chains are available at http://kids.strw.leidenuniv.nl.

We present the X-ray luminosity (L) versus dynamical mass (M) relation for 63 nearby clusters in the HIFLUGCS. The luminosity measurements are obtained based on ~1.3 Ms of clean XMM data and ROSAT pointed observations. The masses are estimated using optical spectroscopic redshifts of 13647 cluster galaxies in total. Given sufficient numbers of member galaxies in computing the dynamical masses, the L-M relations agree between the disturbed and undisturbed clusters. The cool-core clusters still dominate the scatter in the L-M relation even when a core corrected X-ray luminosity is used, which indicates that the scatter mainly reflects the structure formation history of the clusters. As shown by the clusters with a small number of redshifts, the dynamical masses can be underestimated leading to a biased scaling relation. To investigate the potential of spectroscopic surveys to follow up high-redshift galaxy clusters/groups observed in X-ray surveys for the identifications and mass calibrations, we carried out Monte-Carlo re-sampling of the cluster galaxy redshifts and calibrated the uncertainties of the redshift and dynamical mass estimates when only reduced numbers of galaxy redshifts per cluster are available. The re-sampling considers the SPIDERS and 4MOST configurations, designed for the follow-up of the eROSITA clusters, and was carried out for each cluster at the actual cluster redshift as well as at z=0.2, 0.4, 0.6, and 0.8. For following up very distant cluster/groups, we carried out the mass calibration based on the re-sampling with only 10zs/cluster, and redshift calibration based on the re-sampling with only 5zs/cluster and 10zs/cluster, respectively. Our results demonstrate the power of combining upcoming X-ray and optical spectroscopic surveys for mass calibration. The scatter in the dynamical mass estimates for the clusters with at least ten members is within 50%.

Strong gravitational lensing is regarded as the most precise technique to measure the mass in the inner region of galaxies or galaxy clusters. In particular, the mass within one Einstein radius can be determined with an accuracy of order of a few percent or better, depending on the image configuration. For other radii, however, degeneracies exist between galaxy density profiles, precluding an accurate determination of the enclosed mass. The source position transformation (SPT), which includes the well-known mass-sheet transformation (MST) as a special case, describes this degeneracy of the lensing observables in a more general way. In this paper we explore properties of an SPT, removing the MST to leading order, i.e., we consider degeneracies which have not been described before. The deflection field $\boldsymbol{\hat{\alpha}}(\boldsymbol{\theta})$ resulting from an SPT is not curl-free in general, and thus not a deflection that can be obtained from a lensing mass distribution. Starting from a variational principle, we construct lensing potentials that give rise to a deflection field $\boldsymbol{\tilde{\alpha}}$, which differs from $\boldsymbol{\hat{\alpha}}$ by less than an observationally motivated upper limit. The corresponding mass distributions from these 'valid' SPTs are studied: their radial profiles are modified relative to the original mass distribution in a significant and non-trivial way, and originally axi-symmetric mass distributions can obtain a finite ellipticity. These results indicate a significant effect of the SPT on quantitative analyses of lens systems. We show that the mass inside the Einstein radius of the original mass distribution is conserved by the SPT; hence, as is the case for the MST, the SPT does not affect the mass determination at the Einstein radius. [...]

Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91,195 galaxies, within 0.039 < z < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

The ultraviolet (UV) spectral energy distributions of low-mass (K- and M-type) stars play a critical role in the heating and chemistry of exoplanet atmospheres, but are not observationally well-constrained. Direct observations of the intrinsic flux of the Lyman alpha line (the dominant source of UV photons from low-mass stars) are challenging, as interstellar HI absorbs the entire line core for even the closest stars. To address the existing gap in empirical constraints on the UV flux of K and M dwarfs, the MUSCLES HST Treasury Survey has obtained UV observations of 11 nearby M and K dwarfs hosting exoplanets. This paper presents the Lyman alpha and extreme-UV spectral reconstructions for the MUSCLES targets. Most targets are optically inactive, but all exhibit significant UV activity. We use a Markov Chain Monte Carlo technique to correct the observed Lyman alpha profiles for interstellar absorption, and we employ empirical relations to compute the extreme-UV spectral energy distribution from the intrinsic Lyman alpha flux in ~100 {\AA} bins from 100-1170 {\AA}. The reconstructed Lyman alpha profiles have 300 km/s broad cores, while >1% of the total intrinsic Lyman alpha flux is measured in extended wings between 300 km/s to 1200 km/s. The Lyman alpha surface flux positively correlates with the MgII surface flux and negatively correlates with the stellar rotation period. Stars with larger Lyman alpha surface flux also tend to have larger surface flux in ions formed at higher temperatures, but these correlations remain statistically insignificant in our sample of 11 stars. We also present HI column density measurements for 10 new sightlines through the local interstellar medium.

We derive a new relation between cosological distances, valid in any (statistically) isotropic space-time and independent of cosmological parameters or even the validity of the field equation of General Relativity. In particular, this relation yields an equation between those distance ratios which are the geometrical factors determining the strength of the gravitational lensing effect of mass concentrations. Considering a combination of weak lensing shear ratios, based on lenses at two different redshifts, and sources at three different redshifts, we derive a relation between shear-ratio tests which must be identically satisfied. A redshift-dependent multiplicative bias in shear estimates will violate this relation, and thus can be probed by this generalized shear-ratio test. Combining the lensing effect for lenses at three different redshifts and three different source redshifts, a relation between shear ratios is derived which must be valid independent of a multiplicative bias. We propose these generalized shear-ratio tests as a diagnostic for the presence of systematics in upcoming weak lensing surveys.

Ground- and space-based planet searches employing radial velocity techniques and transit photometry have detected thousands of planet-hosting stars in the Milky Way. The chemistry of these atmospheres is controlled by the shape and absolute flux of the stellar spectral energy distribution, however, flux distributions of relatively inactive low-mass stars are poorly known at present. To better understand exoplanets orbiting low-mass stars, we have executed a panchromatic (X-ray to mid-IR) study of the spectral energy distributions of 11 nearby planet hosting stars, the {\it Measurements of the Ultraviolet Spectral Characteristics of Low-mass Exoplanetary Systems} (MUSCLES) Treasury Survey. The MUSCLES program consists of contemporaneous observations at X-ray, UV, and optical wavelengths. We show that energetic radiation (X-ray and ultraviolet) is present from magnetically active stellar atmospheres at all times for stars as late as M5. Emission line luminosities of \ion{C}{4} and \ion{Mg}{2} are strongly correlated with band-integrated luminosities. We find that while the slope of the spectral energy distribution, FUV/NUV, increases by approximately two orders of magnitude form early K to late M dwarfs ($\approx$~0.01~to~1), the absolute FUV and XUV flux levels at their corresponding habitable zone distances are constant to within factors of a few, spanning the range 10~--~70 erg cm$^{-2}$ s$^{-1}$ in the habitable zone. Despite the lack of strong stellar activity indicators in their optical spectra, several of the M dwarfs in our sample show spectacular flare emission in their UV light curves. Finally, we interpret enhanced $L(line)$/$L_{Bol}$ ratios for \ion{C}{4} and \ion{N}{5} as tentative observational evidence for the interaction of planets with large planetary mass-to-orbital distance ratios ($M_{plan}$/$a_{plan}$) with the transition regions of their host stars.

We present measurements of the galaxy bias $b$ and the galaxy-matter cross-correlation coefficient $r$ for the BOSS LOWZ luminous red galaxy sample. Using a new statistical weak lensing analysis of the Red Sequence Cluster Lensing Survey (RCSLenS) we find the bias properties of this sample to be higher than previously reported with $b=2.45^{+0.05}_{-0.05}$ and $r=1.64^{+0.17}_{-0.16}$ on scales between $3'$ and $20'$. We repeat the measurement for angular scales of $20'\leq \vartheta \leq70'$, which yields $b=2.39^{+0.07}_{-0.07}$ and $r=1.24^{+0.26}_{-0.25}$. This is the first application of a data compression analysis using a complete set of discrete estimators for galaxy-galaxy lensing and galaxy clustering. As cosmological data sets grow, our new method of data compression will become increasingly important in order to interpret joint weak lensing and galaxy clustering measurements and to estimate the data covariance. In future studies this formalism can be used as a tool to study the large-scale structure of the Universe to yield a precise determination of cosmological parameters.

Earth's atmosphere imprints a large number of telluric absorption and emission lines on astronomical spectra, especially in the near infrared, that need to be removed before analysing the affected wavelength regions. These lines are typically removed by comparison to A- or B-type stars used as telluric standards that themselves have strong hydrogen lines, which complicates the removal of telluric lines. We have developed a method to circumvent that problem. For our IDL software package tellrem we used a recent approach to model telluric absorption features with the line-by-line radiative transfer model (LBLRTM). The broad wavelength coverage of the X-Shooter at VLT allows us to expand their technique by determining the abundances of the most important telluric molecules H2O, O2, CO2, and CH4 from sufficiently isolated line groups. For individual observations we construct a telluric absorption model for most of the spectral range that is used to remove the telluric absorption from the object spectrum. We remove telluric absorption from both continuum regions and emission lines without systematic residuals for most of the processable spectral range; however, our method increases the statistical errors. The errors of the corrected spectrum typically increase by 10% for S/N~10 and by a factor of two for high-quality data (S/N~100), i.e. the method is accurate on the percent level. Modelling telluric absorption can be an alternative to the observation of standard stars for removing telluric contamination.

Gas planets in close proximity to their host stars experience photoevaporative mass loss. The energy-limited escape concept is generally used to derive estimates for the planetary mass-loss rates. Our photoionization hydrodynamics simulations of the thermospheres of hot gas planets show that the energy-limited escape concept is valid only for planets with a gravitational potential lower than $\log_\mathrm{10}\left( -\Phi_{\mathrm{G}}\right) < 13.11~$erg$\,$g$^{-1}$ because in these planets the radiative energy input is efficiently used to drive the planetary wind. Massive and compact planets with $\log_\mathrm{10}\left( -\Phi_{\mathrm{G}}\right) \gtrsim 13.6~$erg$\,$g$^{-1}$ exhibit more tightly bound atmospheres in which the complete radiative energy input is re-emitted through hydrogen Ly$\alpha$ and free-free emission. These planets therefore host hydrodynamically stable thermospheres. Between these two extremes the strength of the planetary winds rapidly declines as a result of a decreasing heating efficiency. Small planets undergo enhanced evaporation because they host expanded atmospheres that expose a larger surface to the stellar irradiation. We present scaling laws for the heating efficiency and the expansion radius that depend on the gravitational potential and irradiation level of the planet. The resulting revised energy-limited escape concept can be used to derive estimates for the mass-loss rates of super-Earth-sized planets as well as massive hot Jupiters with hydrogen-dominated atmospheres.

Absorption of high-energy radiation in planetary thermospheres is believed to lead to the formation of planetary winds. The resulting mass-loss rates can affect the evolution, particularly of small gas planets. We present 1D, spherically symmetric hydrodynamic simulations of the escaping atmospheres of 18 hot gas planets in the solar neighborhood. Our sample only includes strongly irradiated planets, whose expanded atmospheres may be detectable via transit spectroscopy. The simulations were performed with the PLUTO-CLOUDY interface, which couples a detailed photoionization and plasma simulation code with a general MHD code. We study the thermospheric escape and derive improved estimates for the planetary mass-loss rates. Our simulations reproduce the temperature-pressure profile measured via sodium D absorption in HD 189733 b, but show unexplained differences in the case of HD 209458 b. In contrast to general assumptions, we find that the gravitationally more tightly bound thermospheres of massive and compact planets, such as HAT-P-2 b are hydrodynamically stable. Compact planets dispose of the radiative energy input through hydrogen Ly$\alpha$ and free-free emission. Radiative cooling is also important in HD 189733 b, but it decreases toward smaller planets like GJ 436 b. The simulations show that the strong and cool winds of smaller planets mainly cause strong Ly$\alpha$ absorption but little emission. Compact and massive planets with hot, stable thermospheres cause small absorption signals but are strong Ly$\alpha$ emitters, possibly detectable with the current instrumentation. The absorption and emission signals provide a possible distinction between these two classes of thermospheres in hot gas planets. According to our results, WASP-80 and GJ 3470 are currently the most promising targets for observational follow-up aimed at detecting atmospheric Ly$\alpha$ absorption signals.

Young stars are surrounded by copious amounts of circumstellar material. Its composition, in particular its gas-to-dust ratio, is an important parameter. However, measuring this ratio is challenging, because gas mass estimates are often model dependent. X-ray absorption is sensitive to the gas along the line-of-sight while optical/near-IR extinction depends on the dust content. Therefore, the gas-to-dust ratio of an absorber is given by the ratio between X-ray and optical/near-IR extinction. We present three systems where we used X-ray and optical/near-IR data to constrain the gas-to-dust ratio of circumstellar material; from a dust-rich debris disk to gaseous protoplanetary disks.

The binary system RW Aur consists of two classical T~Tauri stars (CTTSs). The primary recently underwent its second observed major dimming event ($\Delta V\,\sim2\,$mag). We present new, resolved Chandra X-ray and UKIRT near-IR (NIR) data as well as unresolved optical photometry obtained in the dim state to study the gas and dust content of the absorber causing the dimming. The X-ray data show that the absorbing column density increased from $N_H<0.1\times10^{22}\,$cm$^{-2}$ during the bright state to $\approx2\times10^{22}\,$cm$^{-2}$ in the dim state. The brightness ratio between dim and bright state at optical to NIR wavelengths shows only a moderate wavelength dependence and the NIR color-color diagram suggests no substantial reddening. Taken together, this indicates gray absorption by large grains ($\gtrsim1\,\mu$m) with a dust mass column density of $\gtrsim2\times10^{-4}\,$g$\,$cm$^{-2}$. Comparison with $N_H$ shows that an absorber responsible for the optical/NIR dimming and the X-ray absorption is compatible with the ISM's gas-to-dust ratio, i.e., that grains grow in the disk surface layers without largely altering the gas-to-dust ratio. Lastly, we discuss a scenario in which a common mechanism can explain the long-lasting dimming in RW Aur and recently in AA Tau.

Context. In previous work, we developed a quasi-Gaussian approximation for the likelihood of correlation functions, which, in contrast to the usual Gaussian approach, incorporates fundamental mathematical constraints on correlation functions. The analytical computation of these constraints is only feasible in the case of correlation functions of one-dimensional random fields. Aims. In this work, we aim to obtain corresponding constraints in the case of higher-dimensional random fields and test them in a more realistic context. Methods. We develop numerical methods to compute the constraints on correlation functions which are also applicable for two- and three-dimensional fields. In order to test the accuracy of the numerically obtained constraints, we compare them to the analytical results for the one-dimensional case. Finally, we compute correlation functions from the halo catalog of the Millennium Simulation, check whether they obey the constraints, and examine the performance of the transformation used in the construction of the quasi-Gaussian likelihood. Results. We find that our numerical methods of computing the constraints are robust and that the correlation functions measured from the Millennium Simulation obey them. Despite the fact that the measured correlation functions lie well inside the allowed region of parameter space, i.e. far away from the boundaries of the allowed volume defined by the constraints, we find strong indications that the quasi-Gaussian likelihood yields a substantially more accurate description than the Gaussian one.

AA Tau is a well-studied, nearby classical T Tauri star, which is viewed almost edge-on. A warp in its inner disk periodically eclipses the central star, causing a clear modulation of its optical light curve. The system underwent a major dimming event beginning in 2011 caused by an extra absorber, which is most likely associated with additional disk material in the line of sight toward the central source. We present new XMM-Newton X-ray, Hubble Space Telescope FUV, and ground based optical and near-infrared data of the system obtained in 2013 during the long-lasting dim phase. The line width decrease of the fluorescent H$_2$ disk emission shows that the extra absorber is located at $r>1\,$au. Comparison of X-ray absorption ($N_H$) with dust extinction ($A_V$), as derived from measurements obtained one inner disk orbit (eight days) after the X-ray measurement, indicates that the gas-to-dust ratio as probed by the $N_H$ to $A_V$ ratio of the extra absorber is compatible with the ISM ratio. Combining both results suggests that the extra absorber, i.e., material at $r>1\,$au, has no significant gas excess in contrast to the elevated gas-to-dust ratio previously derived for material in the inner region ($\lesssim0.1\,$au).

The Kilo-Degree Survey (KiDS) is a multi-band imaging survey designed for cosmological studies from weak lensing and photometric redshifts. It uses the ESO VLT Survey Telescope with its wide-field camera OmegaCAM. KiDS images are taken in four filters similar to the SDSS ugri bands. The best-seeing time is reserved for deep r-band observations that reach a median 5-sigma limiting AB magnitude of 24.9 with a median seeing that is better than 0.7arcsec. Initial KiDS observations have concentrated on the GAMA regions near the celestial equator, where extensive, highly complete redshift catalogues are available. A total of 109 survey tiles, one square degree each, form the basis of the first set of lensing analyses, which focus on measurements of halo properties of GAMA galaxies. 9 galaxies per square arcminute enter the lensing analysis, for an effective inverse shear variance of 69 per square arcminute. Accounting for the shape measurement weight, the median redshift of the sources is 0.53. KiDS data processing follows two parallel tracks, one optimized for galaxy shape measurement (for weak lensing), and one for accurate matched-aperture photometry in four bands (for photometric redshifts). This technical paper describes how the lensing and photometric redshift catalogues have been produced (including an extensive description of the Gaussian Aperture and Photometry pipeline), summarizes the data quality, and presents extensive tests for systematic errors that might affect the lensing analyses. We also provide first demonstrations of the suitability of the data for cosmological measurements, and explain how the shear catalogues were blinded to prevent confirmation bias in the scientific analyses. The KiDS shear and photometric redshift catalogues, presented in this paper, are released to the community through http://kids.strw.leidenuniv.nl .

We use the first 100 sq. deg. of overlap between the Kilo-Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey to determine the galaxy halo mass of ~10,000 spectroscopically-confirmed satellite galaxies in massive ($M > 10^{13}h^{-1}{\rm M}_\odot$) galaxy groups. Separating the sample as a function of projected distance to the group centre, we jointly model the satellites and their host groups with Navarro-Frenk-White (NFW) density profiles, fully accounting for the data covariance. The probed satellite galaxies in these groups have total masses $\log M_{\rm sub} /(h^{-1}{\rm M}_\odot) \approx 11.7 - 12.2$ consistent across group-centric distance within the errorbars. Given their typical stellar masses, $\log M_{\rm \star,sat}/(h^{-2}{\rm M}_\odot) \sim 10.5$, such total masses imply stellar mass fractions of $M_{\rm \star,sat} /M_{\rm sub} \approx 0.04 h^{-1}$ . The average subhalo hosting these satellite galaxies has a mass $M_{\rm sub} \sim 0.015M_{\rm host}$ independent of host halo mass, in broad agreement with the expectations of structure formation in a $\Lambda$CDM universe.

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)

The Kilo-Degree Survey (KiDS) is an optical wide-field survey designed to map the matter distribution in the Universe using weak gravitational lensing. In this paper, we use these data to measure the density profiles and masses of a sample of $\sim \mathrm{1400}$ spectroscopically identified galaxy groups and clusters from the Galaxy And Mass Assembly (GAMA) survey. We detect a highly significant signal (signal-to-noise-ratio $\sim$ 120), allowing us to study the properties of dark matter haloes over one and a half order of magnitude in mass, from $M \sim 10^{13}-10^{14.5} h^{-1}\mathrm{M_{\odot}}$. We interpret the results for various subsamples of groups using a halo model framework which accounts for the mis-centring of the Brightest Cluster Galaxy (used as the tracer of the group centre) with respect to the centre of the group's dark matter halo. We find that the density profiles of the haloes are well described by an NFW profile with concentrations that agree with predictions from numerical simulations. In addition, we constrain scaling relations between the mass and a number of observable group properties. We find that the mass scales with the total r-band luminosity as a power-law with slope $1.16 \pm 0.13$ (1-sigma) and with the group velocity dispersion as a power-law with slope $1.89 \pm 0.27$ (1-sigma). Finally, we demonstrate the potential of weak lensing studies of groups to discriminate between models of baryonic feedback at group scales by comparing our results with the predictions from the Cosmo-OverWhelmingly Large Simulations (Cosmo-OWLS) project, ruling out models without AGN feedback.

Finding a sample of the most massive clusters with redshifts $z>0.6$ can provide an interesting consistency check of the $\Lambda$ cold dark matter ($\Lambda$CDM) model. Here we present results from our search for clusters with $0.6\lesssim z\lesssim1.0$ where the initial candidates were selected by cross-correlating the RASS faint and bright source catalogues with red galaxies from the Sloan Digital Sky Survey DR8. Our survey thus covers $\approx10,000\,\rm{deg^2}$, much larger than previous studies of this kind. Deeper follow-up observations in three bands using the William Herschel Telescope and the Large Binocular Telescope were performed to confirm the candidates, resulting in a sample of 44 clusters for which we present richnesses and red sequence redshifts, as well as spectroscopic redshifts for a subset. At least two of the clusters in our sample are comparable in richness to RCS2-$J$232727.7$-$020437, one of the richest systems discovered to date. We also obtained new observations with the Combined Array for Research in Millimeter Astronomy for a subsample of 21 clusters. For 11 of those we detect the Sunyaev-Zel'dovich effect signature. The Sunyaev-Zel'dovich signal allows us to estimate $M_{200}$ and check for tension with the cosmological standard model. We find no tension between our cluster masses and the $\Lambda$CDM model.

Exoplanet science is now in its full expansion, particularly after the CoRoT and Kepler space missions that led us to the discovery of thousands of extra-solar planets. The last decade has taught us that UV observations play a major role in advancing our understanding of planets and of their host stars, but the necessary UV observations can be carried out only by HST, and this is going to be the case for many years to come. It is therefore crucial to build a treasury data archive of UV exoplanet observations formed by a dozen "golden systems" for which observations will be available from the UV to the infrared. Only in this way we will be able to fully exploit JWST observations for exoplanet science, one of the key JWST science case.

We present an interface between the (magneto-) hydrodynamics code PLUTO and the plasma simulation and spectral synthesis code CLOUDY. By combining these codes, we constructed a new photoionization hydrodynamics solver: The PLUTO-CLOUDY Interface (TPCI), which is well suited to simulate photoevaporative flows under strong irradiation. The code includes the electromagnetic spectrum from X-rays to the radio range and solves the photoionization and chemical network of the 30 lightest elements. TPCI follows an iterative numerical scheme: First, the equilibrium state of the medium is solved for a given radiation field by CLOUDY, resulting in a net radiative heating or cooling. In the second step, the latter influences the (magneto-) hydrodynamic evolution calculated by PLUTO. Here, we validated the one-dimensional version of the code on the basis of four test problems: Photoevaporation of a cool hydrogen cloud, cooling of coronal plasma, formation of a Stroemgren sphere, and the evaporating atmosphere of a hot Jupiter. This combination of an equilibrium photoionization solver with a general MHD code provides an advanced simulation tool applicable to a variety of astrophysical problems.

Giant gas planets in close proximity to their host stars experience strong irradiation. In extreme cases photoevaporation causes a transonic, planetary wind and the persistent mass loss can possibly affect the planetary evolution. We have identified nine hot Jupiter systems in the vicinity of the Sun, in which expanded planetary atmospheres should be detectable through Lyman alpha transit spectroscopy according to predictions. We use X-ray observations with Chandra and XMM-Newton of seven of these targets to derive the high-energy irradiation level of the planetary atmospheres and the resulting mass loss rates. We further derive improved Lyman alpha luminosity estimates for the host stars including interstellar absorption. According to our estimates WASP-80 b, WASP-77 b, and WASP-43 b experience the strongest mass loss rates, exceeding the mass loss rate of HD 209458 b, where an expanded atmosphere has been confirmed. Furthermore, seven out of nine targets might be amenable to Lyman alpha transit spectroscopy. Finally, we check the possibility of angular momentum transfer from the hot Jupiters to the host stars in the three binary systems among our sample, but find only weak indications for increased stellar rotation periods of WASP-77 and HAT-P-20.

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

Lensing flux-ratio anomalies are most likely caused by gravitational lensing by small-scale dark matter structures. These anomalies offer the prospect of testing a fundamental prediction of the cold dark matter (CDM) cosmological model: the existence of numerous substructures that are too small to host visible galaxies. In two previous studies we found that the number of subhalos in the six high-resolution simulations of CDM galactic halos of the Aquarius project is not sufficient to account for the observed frequency of flux ratio anomalies seen in selected quasars from the CLASS survey. These studies were limited by the small number of halos used, their narrow range of masses (1-2E12 solar masses) and the small range of lens ellipticities considered. We address these shortcomings by investigating the lensing properties of a large sample of halos with a wide range of masses in two sets of high resolution simulations of cosmological volumes and comparing them to a currently best available sample of radio quasars. We find that, as expected, substructures do not change the flux-ratio probability distribution of image pairs and triples with large separations, but they have a significant effect on the distribution at small separations. For such systems, CDM substructures can account for a substantial fraction of the observed flux-ratio anomalies. For large close-pair separation systems, the discrepancies existing between the observed flux ratios and predictions from smooth halo models are attributed to simplifications inherent in these models which do not take account of fine details in the lens mass distributions.

Strong gravitational lensing of sources with different redshifts has been used to determine cosmological distance ratios, which in turn depend on the expansion history. Hence, such systems are viewed as potential tools for constraining cosmological parameters. Here we show that in lens systems with two distinct source redshifts, of which the nearest one contributes to the light deflection towards the more distant one, there exists an invariance transformation which leaves all strong lensing observables unchanged (except the product of time delay and Hubble constant), generalizing the well-known mass-sheet transformation in single plane lens systems. The transformation preserves the relative distribution of mass and light, so that a `mass-follows-light' assumption does not fix the MST. All time delays (from sources on both planes) scale with the same factor -- time-delay ratios are therefore invariant under the MST. Changing cosmological parameters, and thus distance ratios, is essentially equivalent to such a mass-sheet transformation. As an example, we discuss the double source plane system SDSSJ0946+1006, which has been recently studied by Collett and Auger, and show that variations of cosmological parameters within reasonable ranges lead to only a small mass-sheet transformation in both lens planes. Hence, the ability to extract cosmological information from such systems depends heavily on the ability to break the mass-sheet degeneracy.

Young stars accrete mass from circumstellar disks and in many cases, the accretion coincides with a phase of massive outflows, which can be highly collimated. Those jets emit predominantly in the optical and IR wavelength range. However, in several cases X-ray and UV observations reveal a weak but highly energetic component in those jets. X-rays are observed both from stationary regions close to the star and from knots in the jet several hundred AU from the star. In this article we show semi-analytically that a fast stellar wind which is recollimated by the pressure from a slower, more massive disk wind can have the right properties to power stationary X-ray emission. The size of the shocked regions is compatible with observational constraints. Our calculations support a wind-wind interaction scenario for the high energy emission near the base of YSO jets. For the specific case of DG Tau, a stellar wind with a mass loss rate of $5\cdot10^{-10}\;M_{\odot}\;\mathrm{yr}^{-1}$ and a wind speed of 800 km s$^{-1}$ reproduces the observed X-ray spectrum. We conclude that a stellar wind recollimation shock is a viable scenario to power stationary X-ray emission close to the jet launching point.

Future large scale cosmological surveys will provide huge data sets whose analysis requires efficient data compression. Calculating accurate covariances is extremely challenging with increasing number of statistics used. Here we introduce a formalism for achieving efficient data compression, based on a local expansion of statistical measures around a fiducial cosmological model. We specifically apply and test this approach for the case of cosmic shear statistics. We demonstrate the performance of our approach, using a Fisher analysis on cosmic shear tomography described in terms of E-/B-mode separating statistics (COSEBIs). We show that our data compression is highly effective in extracting essentially the full cosmological information from a strongly reduced number of observables. Specifically, the number of statistics needed decreases by at least one order of magnitude relative to the COSEBIs, which already compress the data substantially compared to the shear two-point correlation functions. The efficiency appears to be affected only slightly if a highly inaccurate covariance is used for defining the compressed statistics, showing the robustness of the method. We conclude that an efficient data compression is achievable and that the number of compressed statistics depends on the number of model parameters. In addition, we study how well band powers can be obtained from measuring shear correlation functions over a finite interval of separations. We show the strong limitations on the possibility to construct top-hat filters in Fourier space, for which the real-space analog has a finite support, yielding strong bounds on the accuracy of band power estimates. The error on an estimated band-power is larger for a narrower filter and a smaller angular range which for relevant cases can be as large as 10%.

We consider several aspects of the generalized multi-plane gravitational lens theory, in which light rays from a distant source are affected by several main deflectors, and in addition by the tidal gravitational field of the large-scale matter distribution in the Universe when propagating between the main deflectors. Specifically, we derive a simple expression for the time-delay function in this case, making use of the general formalism for treating light propagation in inhomogeneous spacetimes which leads to the characterization of distance matrices between main lens planes. Applying Fermat's principle, an alternative form of the corresponding lens equation is derived, which connects the impact vectors in three consecutive main lens planes, and we show that this form of the lens equation is equivalent to the more standard one. For this, some general relations for cosmological distance matrices are derived. The generalized multi-plane lens situation admits a generalized mass-sheet transformation, which corresponds to uniform isotropic scaling in each lens plane, a corresponding scaling of the deflection angle, and the addition of a tidal matrix (mass sheet plus external shear) to each main lens. We show that the time delay for sources in all lens planes scale with the same factor under this generalized mass-sheet transformation, thus precluding the use of time-delay ratios to break the mass-sheet transformation.

Classical T Tauri stars (CTTS) are young (< 10 Myr), cool stars that actively accrete matter from a disk. They show strong, broad and asymmetric, atomic FUV emission lines. Neither the width, nor the line profile is understood. Likely, different mechanisms influence the line profile; the best candidates are accretion, winds and stellar activity. We monitored the C IV 1548/1550 Ang doublet in the nearby, bright CTTS TW Hya with the Hubble Space Telescope Cosmic Origin Spectrograph (HST/COS) to correlate it with i) the cool wind, as seen in COS NUV Mg II line profiles, ii) the photometric period from joint ground-based monitoring, iii) the accretion rate as determined from the UV continuum, and iv) the Ha line profile from independent ground-based observations. The observations span 10 orbits distributed over a few weeks to cover the typical time scales of stellar rotation, accretion and winds. Here we describe a model with intrinsically asymmetric C IV lines.

HD 181068 is the only compact, triply eclipsing, hierarchical triple system containing a giant star known to date. With its central, highly-active G-type giant orbited by a close pair of main-sequence dwarfs, the system is ideal to study tidal interactions. We carried out a multiwavelength study to characterize the magnetic activity of the HD 181068 system. To this end, we obtained in- and out-of-eclipse X-ray snapshots with XMM-Newton and an optical spectrum, which we analyzed along with the Kepler light-curve. The primary giant shows strong quiescent X-ray emission at a level of 2e31 ergs, an S-index of 0.41 +/- 0.01, and marked white-light flares releasing up to 6e38 erg in the Kepler-band. During the second X-ray observation, we found a three-times elevated -- yet decaying -- level of X-ray emission, which might be due to an X-ray flare. The high level of magnetic activity is compatible with the previously reported absence of solar-like oscillations in the giant, whose atmosphere, however, undergoes tidally-induced oscillations imposed by the changing configuration of the dwarf-binary. We found that the driving force exciting these oscillations is comparable to the disturbances produced by a typical hot Jupiter, making the system a potential test bed to study the effects of tidal interactions also present in planetary systems.

Developing accurate analysis techniques to combine various probes of cosmology is essential to tighten constraints on cosmological parameters and to check for inconsistencies in our model of the Universe. In this paper we develop a joint analysis framework for six different second-order statistics calculated from three tracers of the dark matter density field, namely galaxy position, shear, and magnification. We extend a data compression scheme developed in the context of shear-shear statistics (the so-called COSEBIs) to the other five second-order statistics, thereby significantly reducing the number of data points in the joint data vector. We use CosmoLike, a newly developed software framework for joint likelihood analyses, to forecast parameter constraints for the Dark Energy Survey (DES). The simulated MCMCs cover a five dimensional cosmological parameter space comparing the information content of the individual probes to several combined probes (CP) data vectors. Given the significant correlations of these second-order statistics we model all cross terms in the covariance matrix; furthermore we go beyond the Gaussian covariance approximation and use the halo model to include higher order correlations of the density field. We find that adding magnification information (including cross probes with shear and clustering) noticeably increases the information content and that the correct modeling of the covariance (i.e., accounting for non-Gaussianity and cross terms) is essential for accurate likelihood contours from the CP data vector. We also identify several nulltests based on the degeneracy of magnification and shear statistics which can be used to quantify the contamination of data sets by astrophysical systematics and/or calibration issues.

Scaling properties of galaxy cluster observables with mass provide central insights into the processes shaping clusters. Calibrating proxies for cluster mass will be crucial to cluster cosmology with upcoming surveys like eROSITA and Euclid. The recent Planck results led to suggestions that X-ray masses might be biased low by $\sim\!40$ %, more than previously considered. We extend the direct calibration of the weak lensing -- X-ray mass scaling towards lower masses (as low as $1\!\times\!10^{14}\,\mathrm{M}_{\odot}$) in a sample representative of the $z\!\sim\!0.4$--$0.5$ population. We investigate the scaling of MMT/Megacam weak lensing (WL) masses for $8$ clusters at $0.39\!\leq\!z\!\leq\!0.80$ as part of the \emph{400d} WL programme with hydrostatic \textit{Chandra} X-ray masses as well as those based on the proxies, e.g. $Y_{\mathrm{X}}\!=\!T_{\mathrm{X}}M_{\mathrm{gas}}$. Overall, we find good agreement between WL and X-ray masses, with different mass bias estimators all consistent with zero. Subdividing the sample, we find the high-mass subsample to show no significant mass bias while for the low-mass subsample, there is a bias towards overestimated X-ray masses at the $\sim\!2\sigma$ level for some mass proxies. The overall scatter in the mass-mass scaling relations is surprisingly low. Neither observation can be traced back to the parameter settings in the WL analysis. We do not find evidence for a strong ($\sim\!40$ %) underestimate in the X-ray masses, as suggested to reconcile Planck cluster counts and cosmological constraints. For high-mass clusters, our measurements are consistent with studies in the literature. The mass dependent bias, significant at $\sim\!2\sigma$, may hint at a physically different cluster population (less relaxed clusters with more substructure and mergers); or it may be due to small number statistics.

Intrinsic galaxy alignments constitute the major astrophysical systematic of forthcoming weak gravitational lensing surveys but also yield unique insights into galaxy formation and evolution. We build analytic models for the distribution of galaxy shapes based on halo properties extracted from the Millennium Simulation, differentiating between early- and late-type galaxies as well as central galaxies and satellites. The resulting ellipticity correlations are investigated for their physical properties and compared to a suite of current observations. The best-faring model is then used to predict the intrinsic alignment contamination of planned weak lensing surveys. We find that late-type galaxy models generally have weak intrinsic ellipticity correlations, marginally increasing towards smaller galaxy separation and higher redshift. The signal for early-type models at fixed halo mass strongly increases by three orders of magnitude over two decades in galaxy separation, and by one order of magnitude from z=0 to z=2. The intrinsic alignment strength also depends strongly on halo mass, but not on galaxy luminosity at fixed mass, or galaxy number density in the environment. We identify models that are in good agreement with all observational data, except that all models over-predict alignments of faint early-type galaxies. The best model yields an intrinsic alignment contamination of a Euclid-like survey between 0.5-10% at z>0.6 and on angular scales larger than a few arcminutes. Cutting 20% of red foreground galaxies using observer-frame colours can suppress this contamination by up to a factor of two.

With 3rd-order statistics of gravitational shear it will be possible to extract valuable cosmological information from ongoing and future weak lensing surveys which is not contained in standard 2nd-order statistics, due to the non-Gaussianity of the shear field. Aperture mass statistics are an appropriate choice for 3rd-order statistics due to their simple form and their ability to separate E- and B-modes of the shear. However, it has been demonstrated that 2nd-order aperture mass statistics suffer from E-/B-mode mixing because it is impossible to reliably estimate the shapes of close pairs of galaxies. This finding has triggered developments of several new 2nd-order statistical measures for cosmic shear. Whether the same developments are needed for 3rd-order shear statistics is largely determined by how severe this E-/B-mixing is for 3rd-order statistics. We test 3rd-order aperture mass statistics against E-/B-mode mixing, and find that the level of contamination is well-described by a function of $\theta/\theta_{\rm min}$, with $\theta_{\rm min}$ being the cut-off scale. At angular scales of $\theta > 10 \;\theta_{\rm min}$, the decrease in the E-mode signal due to E-/B-mode mixing is smaller than 1 percent, and the leakage into B-modes is even less. For typical small-scale cut-offs this E-/B-mixing is negligible on scales larger than a few arcminutes. Therefore, 3rd-order aperture mass statistics can safely be used to separate E- and B-modes and infer cosmological information, for ground-based surveys as well as forthcoming space-based surveys such as Euclid.

Interstellar reddening corrections are necessary to reconstruct the intrinsic spectral energy distributions (SEDs) of accreting protostellar systems. The stellar SED determines the heating and chemical processes that can occur in circumstellar disks. Measurement of neutral hydrogen absorption against broad Lyman-$\alpha$ emission profiles in young stars can be used to obtain the total H I column density (N(H I)) along the line of sight. We measure N(H I) with new and archival ultraviolet observations from the Hubble Space Telescope ($HST$) of 31 classical T Tauri and Herbig Ae/Be stars. The H I column densities range from log$_{10}$(N(H I)) $\approx 19.6 - 21.1$, with corresponding visual extinctions of A$_{V}$ $= 0.02 - 0.72$ mag, assuming an R$_{V}$ of 3.1. We find that the majority of the H I absorption along the line of sight likely comes from interstellar rather than circumstellar material. Extinctions derived from new $HST$ blue-optical spectral analyses, previous IR and optical measurements, and new X-ray column densities on average overestimate the interstellar extinction toward young stars compared to the N(H I) values by $\sim 0.6$ mag. We discuss possible explanations for this discrepancy in the context of a protoplanetary disk geometry.

The light travel time differences in strong gravitational lensing systems allows an independent determination of the Hubble constant. This method has been successfully applied to several lens systems. The formally most precise measurements are, however, in tension with the recent determination of $H_0$ from the Planck satellite for a spatially flat six-parameters $\Lambda CDM$ cosmology. We reconsider the uncertainties of the method, concerning the mass profile of the lens galaxies, and show that the formal precision relies on the assumption that the mass profile is a perfect power law. Simple analytical arguments and numerical experiments reveal that mass-sheet like transformations yield significant freedom in choosing the mass profile, even when exquisite Einstein rings are observed. Furthermore, the characterization of the environment of the lens does not break that degeneracy which is not physically linked to extrinsic convergence. We present an illustrative example where the multiple imaging properties of a composite (baryons + dark matter) lens can be extremely well reproduced by a power-law model having the same velocity dispersion, but with predictions for the Hubble constant that deviate by $\sim 20%$. Hence we conclude that the impact of degeneracies between parametrized models have been underestimated in current $H_0$ measurements from lensing, and need to be carefully reconsidered.

One of the best-studied jets from all young stellar objects is the jet of DG Tau, which we imaged in the FUV with HST for the first time. These high spatial resolution images were obtained with long-pass filters and allow us to construct images tracing mainly molecular hydrogen and C IV emission. We find that the H2 emission appears as a limb-brightened cone with additional emission close to the jet axis. The length of the rims is about 0.3 arcsec or 42 AU (proj.) before their brightness strongly drops, and the opening angle is about 90 deg. Comparing our FUV data with near-IR data we find that the fluorescent H2 emission likely traces the outer, cooler part of the disk wind while an origin of the H2 emission in the surface layers (atmosphere) of the (flared) disk is unlikely. Furthermore, the spatial shape of the H2 emission shows little variation over six years which suggests that the outer part of the disk wind is rather stable and probably not associated with the formation of individual knots. The C IV image shows that the emission is concentrated towards the jet axis. We find no indications for additional C IV emission at larger distances, which strengthens the association with the X-ray emission observed to originate within the DG Tau jet.

Context. Whenever correlation functions are used for inference about cosmological parameters in the context of a Bayesian analysis, the likelihood function of correlation functions needs to be known. Usually, it is approximated as a multivariate Gaussian, though this is not necessarily a good approximation. Aims. We show how to calculate a better approximation for the probability distribution of correlation functions, which we call "quasi-Gaussian". Methods. Using the exact univariate PDF as well as constraints on correlation functions previously derived, we transform the correlation functions to an unconstrained variable for which the Gaussian approximation is well justified. From this Gaussian in the transformed space, we obtain the quasi-Gaussian PDF. The two approximations for the probability distributions are compared to the "true" distribution as obtained from simulations. Additionally, we test how the new approximation performs when used as likelihood in a toy-model Bayesian analysis. Results. The quasi-Gaussian PDF agrees very well with the PDF obtained from simulations; in particular, it provides a significantly better description than a straightforward copula approach. In a simple toy-model likelihood analysis, it yields noticeably different results than the Gaussian likelihood, indicating its possible impact on cosmological parameter estimation.

The main obstacle for gravitational lensing to determine accurate masses of deflectors, or to determine precise estimates for the Hubble constant, is the degeneracy of lensing observables with respect to the mass-sheet transformation (MST). The MST is a global modification of the mass distribution which leaves all image positions, shapes and flux ratios invariant, but which changes the time delay. Here we show that another global transformation of lensing mass distributions exists which almost leaves image positions and flux ratios invariant, and of which the MST is a special case. Whereas for axi-symmetric lenses this source position transformation exactly reproduces all strong lensing observables, it does so only approximately for more general lens situations. We provide crude estimates for the accuracy with which the transformed mass distribution can reproduce the same image positions as the original lens model, and present an illustrative example of its performance. This new invariance transformation most likely is the reason why the same strong lensing information can be accounted for with rather different mass models.

The statistical properties of the ellipticities of galaxy images depend on how galaxies form and evolve, and therefore constrain models of galaxy morphology, which are key to the removal of the intrinsic alignment contamination of cosmological weak lensing surveys, as well as to the calibration of weak lensing shape measurements. We construct such models based on the halo properties of the Millennium Simulation and confront them with a sample of 90,000 galaxies from the COSMOS Survey, covering three decades in luminosity and redshifts out to z=2. The ellipticity measurements are corrected for effects of point spread function smearing, spurious image distortions, and measurement noise. Dividing galaxies into early, late, and irregular types, we find that early-type galaxies have up to a factor of two lower intrinsic ellipticity dispersion than late-type galaxies. None of the samples shows evidence for redshift evolution, while the ellipticity dispersion for late-type galaxies scales strongly with absolute magnitude at the bright end. The simulation-based models reproduce the main characteristics of the intrinsic ellipticity distributions although which model fares best depends on the selection criteria of the galaxy sample. We observe fewer close-to-circular late-type galaxy images in COSMOS than expected for a sample of randomly oriented circular thick disks and discuss possible explanations for this deficit.

Young A and B stars, the so-called Herbig Ae/Be stars (HAeBe), are surrounded by an active accretion disk and drive outflows. We study the jet HH 409, which is launched from the HAeBe star HD 163296, using new and archival observations from Chandra and HST/STIS. In X-rays we can show that the central source is not significantly extended. The approaching jet, but not the counter-jet, is detected in Ly alpha. In addition, there is red-shifted Ly alpha emission extended in the same direction as the jet, that is also absent in the counter-jet. We can rule out an accretion or disk-wind origin for this feature. In the optical we find the knots B and B2 in the counter-jet. Knot B has been observed previously, so we can derive its proper motion of 0.37+-0.01 arcsec/yr. Its electron density is 3000/cm^3, thus the cooling time scale is a few months only, so the knot needs to be reheated continuously. The shock speed derived from models of H alpha and forbidden emission lines (FELs) decreased from 50 km/s in 1999 to 30 km/s in 2011 because the shock front loses energy as it travels along the jet. Knot B2 is observed at a similar position in 2011 as knot B was in 1999, but shows a lower ionization fraction and higher mass loss rate, proving variations in the jet launching conditions.

The active young protostar DG Tau has an extended jet that has been well studied at radio, optical, and X-ray wavelengths. We report sensitive new VLA full-polarization observations of the core and jet between 5 GHz and 8 GHz. Our high angular resolution observation at 8 GHz clearly shows an unpolarized inner jet with a size 42 AU (0.35") extending along a position angle similar to the optical-X ray outer jet. Using our nearly coeval 2012 VLA observations, we find a spectral-index=+0.46+/-0.05, which combined with the lack of polarization, is consistent with bremsstrahlung (free-free) emission, with no evidence for a non-thermal coronal component. By identifying the end of the radio jet as the optical depth unity surface, and calculating the resulting emission measure, we find our radio results are in agreement with previous optical line studies of electron density and consequent mass-loss rate. We also detect a weak radio knot at 5 GHz located 7" from the base of the jet, coincident with the inner radio knot detected by Rodriguez et al. (2012) in 2009 but at lower surface brightness. We interpret this as due to expansion of post-shock ionized gas in the three years between observations.

We study the cosmological information of weak lensing (WL) peaks, focusing on two other statistics besides their abundance: the stacked tangential-shear profiles and the peak-peak correlation function. We use a large ensemble of simulated WL maps with survey specifications relevant to future missions like Euclid and LSST, to explore the three peak probes. We find that the correlation function of peaks with high signal-to-noise (S/N) measured from fields of size 144 sq. deg. has a maximum of ~0.3 at an angular scale ~10 arcmin. For peaks with smaller S/N, the amplitude of the correlation function decreases, and its maximum occurs on smaller angular scales. We compare the peak observables measured with and without shape noise and find that for S/N~3 only ~5% of the peaks are due to large-scale structures, the rest being generated by shape noise. The covariance matrix of the probes is examined: the correlation function is only weakly covariant on scales < 30 arcmin, and slightly more on larger scales; the shear profiles are very correlated for theta > 2 arcmin, with a correlation coefficient as high as 0.7. Using the Fisher-matrix formalism, we compute the cosmological constraints for {Om_m, sig_8, w, n_s} considering each probe separately, as well as in combination. We find that the correlation function of peaks and shear profiles yield marginalized errors which are larger by a factor of 2-4 for {Om_m, sig_8} than the errors yielded by the peak abundance alone, while the errors for {w, n_s} are similar. By combining the three probes, the marginalized constraints are tightened by a factor of ~2 compared to the peak abundance alone, the least contributor to the error reduction being the correlation function. This work therefore recommends that future WL surveys use shear peaks beyond their abundance in order to constrain the cosmological model.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

The gravitational lensing effect provides various ways to study the mass environment of galaxies. We investigate how galaxy-galaxy(-galaxy) lensing can be used to test models of galaxy formation and evolution. We consider two semi-analytic galaxy formation models based on the Millennium Run N-body simulation: the Durham model by Bower et al. (2006) and the Garching model by Guo et al. (2011). We generate mock lensing observations for the two models, and then employ Fast Fourier Transform methods to compute second- and third-order aperture statistics in the simulated fields for various galaxy samples. We find that both models predict qualitatively similar aperture signals, but there are large quantitative differences. The Durham model predicts larger amplitudes in general. In both models, red galaxies exhibit stronger aperture signals than blue galaxies. Using these aperture measurements and assuming a linear deterministic bias model, we measure relative bias ratios of red and blue galaxy samples. We find that a linear deterministic bias is insufficient to describe the relative clustering of model galaxies below ten arcmin angular scales. Dividing galaxies into luminosity bins, the aperture signals decrease with decreasing luminosity for brighter galaxies, but increase again for fainter galaxies. This increase is likely an artifact due to too many faint satellite galaxies in massive group and cluster halos predicted by the models. Our study shows that galaxy-galaxy(-galaxy) lensing is a sensitive probe of galaxy evolution.

We present the first direct measurement of the galaxy-matter bispectrum as a function of galaxy luminosity, stellar mass, and SED type. Our analysis uses a galaxy-galaxy-galaxy lensing technique (G3L), on angular scales between 9 arcsec to 50 arcmin, to quantify (i) the excess surface mass density around galaxy pairs (excess mass hereafter) and (ii) the excess shear-shear correlations around single galaxies, both of which yield a measure of two types of galaxy-matter bispectra. We apply our method to the state-of-the-art Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), spanning 154 square degrees. This survey allows us to detect a significant change of the bispectra with lens properties. Measurements for lens populations with distinct redshift distributions become comparable by a newly devised normalisation technique. That will also aid future comparisons to other surveys or simulations. A significant dependence of the normalised G3L statistics on luminosity within -23<M_r<-18 and stellar mass within 5x10^9 M_sol<M_star<2x10^11 M_sol is found (h=0.73). Both bispectra exhibit a stronger signal for more luminous lenses or those with higher stellar mass (up to a factor 2-3). This is accompanied by a steeper equilateral bispectrum for more luminous or higher stellar mass lenses for the excess mass. Importantly, we find the excess mass to be very sensitive to galaxy type as recently predicted with semi-analytic galaxy models: luminous (M_r<-21) late-type galaxies show no detectable signal, while all excess mass detected for luminous galaxies seems to be associated with early-type galaxies. We also present the first observational constraints on third-order stochastic galaxy biasing parameters.

Third-order galaxy-galaxy lensing (G3L) is a next generation galaxy-galaxy lensing technique that either measures the excess shear about lens pairs or the excess shear-shear correlations about lenses. It is clear that these statistics assess the three-point correlations between galaxy positions and projected matter density. For future applications of these novel statistics, we aim at a more intuitive understanding of G3L to isolate the main features that possibly can be measured. We construct a toy model ("isolated lens model"; ILM) for the distribution of galaxies and associated matter to determine the measured quantities of the two G3L correlation functions and traditional galaxy-galaxy lensing (GGL) in a simplified context. The ILM presumes single lens galaxies to be embedded inside arbitrary matter haloes that, however, are statistically independent ("isolated") from any other halo or lens position. In the ILM, the average mass-to-galaxy number ratio of clusters of any size cannot change. GGL and galaxy clustering alone cannot distinguish an ILM from any more complex scenario. The lens-lens-shear correlator in combination with second-order statistics enables us to detect deviations from a ILM, though. This can be quantified by a difference signal defined in the paper. We demonstrate with the ILM that this correlator picks up the excess matter distribution about galaxy pairs inside clusters. The lens-shear-shear correlator is sensitive to variations among matter haloes. In principle, it could be devised to constrain the ellipticities of haloes, without the need for luminous tracers, or maybe even random halo substructure. [Abridged]

Protostellar jets are tightly connected to the accretion process and regulate the angular momentum balance of accreting star-disk systems. The DG Tau jet is one of the best-studied protostellar jets and contains plasma with temperatures ranging over three orders of magnitude within the innermost 50 AU of the jet. We present new Hubble Space Telescope (HST) far ultraviolet (FUV) long-slit spectra spatially resolving the C IV emission (T~1e5 K) from the jet for the first time, and quasi-simultaneous HST observations of optical forbidden emission lines ([O I], [N II], [S II] and [O III]) and fluorescent H2 lines. The C IV emission peaks at 42 AU from the stellar position and has a FWHM of 52 AU along the jet. Its deprojected velocity of around 200 km/s decreases monotonically away from the driving source. In addition, we compare our HST data with the X-ray emission from the DG Tau jet. We investigate the requirements to explain the data by an initially hot jet compared to local heating. Both scenarios indicate a mass loss by the T~1e5 K jet of ~1e-9 Msun/year, i.e., between the values for the lower temperature jet (T~1e4 K) and the hotter X-ray emitting part (T>1e6 K). However, a simple initially hot wind requires a large launching region (~1 AU), and we therefore favor local heating.

Evolution in the mass function of galaxy clusters sensitively traces both the expansion history of the Universe and cosmological structure formation. Robust cluster mass determinations are a key ingredient for a reliable measurement of this evolution, especially at high redshift. Weak gravitational lensing is a promising tool for, on average, unbiased mass estimates. This weak lensing project aims at measuring reliable weak lensing masses for a complete X-ray selected sample of 36 high redshift (0.35<z<0.9) clusters. The goal of this paper is to demonstrate the robustness of the methodology against commonly encountered problems, including pure instrumental effects, the presence of bright (8--9 mag) stars close to the cluster centre, ground based measurements of high-z (z~0.8) clusters, and the presence of massive unrelated structures along the line-sight. We select a subsample of seven clusters observed with MMT/Megacam. Instrumental effects are checked in detail by cross-comparison with an archival CFHT/MegaCam observation. We derive mass estimates for seven clusters by modelling the tangential shear with an NFW profile, in two cases with multiple components to account for projected structures in the line-of-sight. We firmly detect lensing signals from all seven clusters at more than $3.5\sigma$ and determine their masses, ranging from $10^{14} M_{\odot}$ to $10^{15} M_{\odot}$, despite the presence of nearby bright stars. We retrieve the lensing signal of more than one cluster in the CL 1701+6414 field, while apparently observing CL 1701+6414 through a massive foreground filament. We also find a multi-peaked shear signal in CL 1641+4001. Shear structures measured in the MMT and CFHT images of CL 1701+6414 are highly correlated.

We studied one of the most X-ray luminous cluster of galaxies in the REFLEX survey, RXC J1504.1-0248 (hereafter R1504; z=0.2153), using XMM-Newton X-ray imaging spectroscopy, VLT/VIMOS optical spectroscopy, and WFI optical imaging. The mass distributions were determined using both the so-called hydrostatic method with X-ray imaging spectroscopy and the dynamical method with optical spectroscopy, respectively, which yield M^{H.E.}_{500}=(5.81+/-0.49)*1.e14Msun and M^{caustic}_{500}=(4.17+/-0.42)*1e14Msun. According to recent calibrations, the richness derived mass estimates closely agree with the hydrostatic and dynamical mass estimates. The line-of-sight velocities of spectroscopic members reveal a group of galaxies with high-velocities (>1000 km/s) at a projected distance of about r^{H.E.}_{500}=(1.18+/-0.03) Mpc south-east of the cluster centroid, which is also indicated in the X-ray two-dimensional (2-D) temperature, density, entropy, and pressure maps. The dynamical mass estimate is 80% of the hydrostatic mass estimate at r^{H.E.}_{500}. It can be partially explained by the ~20% scatter in the 2-D pressure map that can be propagated into the hydrostatic mass estimate. The uncertainty in the dynamical mass estimate caused by the substructure of the high velocity group is ~14%. The dynamical mass estimate using blue members is 1.23 times that using red members. The global properties of R1504 obey the observed scaling relations of nearby clusters, although its stellar-mass fraction is rather low.

Young stellar objects often drive powerful bipolar outflows which evolve on time scales of a few years. An increasing number of these outflows has been detected in X-rays implying the existence of million degree plasma almost co-spatial with the lower temperature gas observed in the optical and near-IR. The details of the heating and cooling processes of the X-ray emitting part of these so-called Herbig-Haro objects are still ambiguous, e.g., whether the cooling is dominated by expansion, radiation or thermal conduction. We present a second epoch Chandra observation of the first X-ray detected Herbig-Haro object (HH 2) and derive the proper-motion of the X-ray emitting plasma and its cooling history. We argue that the most likely explanation for the constancy of the X-ray luminosity, the alignment with the optical emission and the proper-motion is that the cooling is dominated by radiative losses leading to cooling times exceeding a decade. We explain that a strong shock caused by fast material ramming into slower gas in front of it about ten years ago can explain the X-ray emission while being compatible with the available multi-wavelength data of HH 2.

We explore the possibility of detecting primordial non-Gaussianity of the local type using weak lensing peak counts. We measure the peak abundance in sets of simulated weak lensing maps corresponding to three models f_NL={0, +100, -100}. Using survey specifications similar to those of Euclid and without assuming any knowledge of the lens and source redshifts, we find the peak functions of the non-Gaussian models with f_NL=+-100 to differ by up to 15% from the Gaussian peak function at the high-mass end. For the assumed survey parameters, the probability of fitting an f_NL=0 peak function to the f_NL=+-100 peak functions is less than 0.1%. Assuming the other cosmological parameters known, f_NL can be measured with an error \Delta f_NL ~ 13. It is therefore possible that future weak lensing surveys like Euclid and LSST may detect primordial non-Gaussianity from the abundance of peak counts, and provide complementary information to that obtained from the cosmic microwave background.

We present a new method to extract cosmological constraints from weak lensing (WL) peak counts, which we denote as `the hierarchical algorithm'. The idea of this method is to combine information from WL maps sequentially smoothed with a series of filters of different size, from the largest down to the smallest, thus increasing the cosmological sensitivity of the resulting peak function. We compare the cosmological constraints resulting from the peak abundance measured in this way and the abundance obtained by using a filter of fixed size, which is the standard practice in WL peak studies. For this purpose, we employ a large set of WL maps generated by ray-tracing through N-body simulations, and the Fisher matrix formalism. We find that if low-S/N peaks are included in the analysis (S/N ~ 3), the hierarchical method yields constraints significantly better than the single-sized filtering. For a large future survey such as Euclid or LSST, combined with information from a CMB experiment like Planck, the results for the hierarchical (single-sized) method are: \Delta n=0.0039 (0.004); \Delta \Omega m=0.002 (0.0045); \Delta \sigma 8=0.003 (0.006); \Delta w=0.019 (0.0525). This forecast is conservative, as we assume no knowledge of the redshifts of the lenses, and consider a single broad bin for the redshifts of the sources. If only peaks with S/N >= 6 are considered, then there is little difference between the results of the two methods. We also examine the statistical properties of the hierarchical peak function: Its covariance matrix has off-diagonal terms for bins with S/N <= 6 and aperture mass of M < 3 x 1e+14 Ms/h, the higher bins being largely uncorrelated and therefore well described by a Poisson distribution.

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.

Decomposing the shear signal into E and B-modes properly, i.e. without leakage of B-modes into the E-mode signal and vice versa, has been a long-standing problem in weak gravitational lensing. At the two-point level this problem was resolved by developing the so-called ring statistics, and later the COSEBIs; however, extending these concepts to the three-point level is far from trivial. Currently used methods to decompose three-point shear correlation functions (3PCFs) into E- and B-modes require knowledge of the 3PCF down to arbitrary small scales. This implies that the 3PCF needs to be modeled on scales smaller than the minimum separation of 2 galaxies and subsequently will be biased towards the model, or, in the absence of a model, the statistics is affected by E/B-mode leakage (or mixing). In this paper we derive a new third-order E/B-mode statistic that performs the decomposition using the 3PCF only on a finite interval, and thereby is free of any E/B-mode leakage while at the same time relying solely on information from the data. In addition, we relate this third-order ring statistics to the convergence field, thereby enabling a fast and convenient calculation of this statistic from numerical simulations. We note that our new statistics should be applicable to corresponding E/B-mode separation problems in the CMB polarization field.

Context. Gravitational lensing is one of the leading tools in understanding the dark side of the Universe. The need for accurate, efficient and effective methods which are able to extract this information along with other cosmological parameters from cosmic shear data is ever growing. COSEBIs, Complete Orthogonal Sets of E-/B-Integrals, is a recently developed statistical measure that encompasses the complete E-/B-mode separable information contained in the shear correlation functions measured on a finite angular range. Aims. The aim of the present work is to test the properties of this newly developed statistics for a higher-dimensional parameter space and to generalize and test it for shear tomography. Methods. We use Fisher analysis to study the effectiveness of COSEBIs. We show our results in terms of figure-of-merit quantities, based on Fisher matrices. Results. We find that a relatively small number of COSEBIs modes is always enough to saturate to the maximum information level. This number is always smaller for 'logarithmic COSEBIs' than for 'linear COSEBIs', and also depends on the number of redshift bins, the number and choice of cosmological parameters, as well as the survey characteristics. Conclusions. COSEBIs provide a very compact way of analyzing cosmic shear data, i.e., all the E-/B-mode separable second-order statistical information in the data is reduced to a small number of COSEBIs modes. Furthermore, with this method the arbitrariness in data binning is no longer an issue since the COSEBIs modes are discrete. Finally, the small number of modes also implies that covariances, and their inverse, are much more conveniently obtainable, e.g., from numerical simulations, than for the shear correlation functions themselves.

It is believed that the global baryon content of clusters of galaxies is representative of the matter distribution of the universe, and can, therefore, be used to reliably determine the matter density parameter Omega_m. This assumption is challenged by the growing evidence from optical and X-ray observations that the total baryon mass fraction increases towards rich clusters. In this context, we investigate the dependence of stellar, and total baryon mass fractions as a function of mass. To do so, we used a subsample of nineteen clusters extracted from the X-ray flux limited sample HIFLUGCS that have available DR-7 Sloan Digital Sky Survey (SDSS) data. From the optical analysis we derived the stellar masses. Using XMM-Newton we derived the gas masses. Then, adopting a scaling relation we estimate the total masses. Adding the gas and the stellar mass fractions we obtain the total baryonic content that we find to increase with cluster mass, reaching 7-year Wilkinson Microwave Anisotropy Probe (WMAP-7) prediction for clusters with M_500 = 1.6 x 10^{15} M_sun. We observe a decrease of the stellar mass fraction (from 4.5% to ~1.0%) with increasing total mass where our findings for the stellar mass fraction agree with previous studies. This result suggests a difference in the number of stars formed per unit of halo mass, though with a large scatter for low-mass systems. That is, the efficiency of star formation varies on cluster scale that lower mass systems are likely to have higher star formation efficiencies. It follows immediately that the dependence of the stellar mass fraction on total mass results in an increase of the mass-to-light ratio from lower to higher mass systems. We also discuss the consequences of these results in the context of determining the cosmic matter density parameter Omega_m.

[Abridged] We seek approximations to the cosmic shear covariance that are as easy to use as the common approximations based on normal statistics, but yield more accurate covariance matrices and parameter errors. We derive expressions for the cosmic shear covariance under the assumption that the underlying convergence field follows log-normal statistics. We also derive a simplified version of this log-normal approximation. We use numerical simulations of weak lensing to study how well the normal, log-normal, and simplified log-normal approximations as well as empirical corrections to the normal approximation proposed in the literature reproduce shear covariances for cosmic shear surveys. We find that the normal approximation substantially underestimates the cosmic shear covariances and the inferred parameter confidence regions, in particular for surveys with small fields of view and large galaxy densities, but also for very wide surveys. In contrast, the log-normal approximation yields more realistic covariances and confidence regions, but also requires evaluating slightly more complicated expressions. However, the simplified log-normal approximation, although as simple as the normal approximation, yields confidence regions that are almost as accurate as those obtained from the log-normal approximation. The empirical corrections to the normal approximation do not yield more accurate covariances and confidence regions than the (simplified) log-normal approximation. Moreover, they fail to produce positive-semidefinite data covariance matrices in certain cases, rendering them unusable for parameter estimation. The simplified log-normal approximation should be used in favour of the normal approximation for parameter estimation and parameter error forecasts. Any approximation to the cosmic shear covariance should ensure a positive-semidefinite data covariance matrix.

We have investigated the baryon-mass content in a subsample of 19 clusters of galaxies extracted from the X-ray flux-limited sample HIFLUGCS according to their positions in the sky. For these clusters, we measured total masses and characteristic radii on the basis of a rich optical spectroscopic data set, the physical properties of the intracluster medium (ICM) using XMM-Newton and ROSAT X-ray data, and total (galaxy) stellar masses utilizing the SDSS DR7 multi-band imaging. The observed (hot) gas-mass fractions are almost constant in this mass range. We confirm that the stellar mass fraction decreases as the total mass increases and shows (20+/-4)% scatter; in addition, we show that it decreases as the central entropy increases. The latter behavior supports a twofold interpretation where heating from merging quenches the star-formation activity of galaxies in massive systems, and feedback from supernovae and/or radio galaxies drives a significant amount of gas to the regions beyond r_{500} or, alternatively, a substantially large amount of intracluster light (ICL) is associated with galaxies in nonrelaxed systems. Furthermore, less massive clusters are confirmed to host less gas per unit total mass; however, they exhibit higher mass fractions in metals, so that their ICM is more metal-rich. This again supports the interpretation that in the potential wells of low-mass systems the star-formation efficiency of galaxies was high or, alternatively, some gas is missing from the hot phase of the ICM. The former hypothesis is preferred as the main driver of the mass-dependent metal enrichment since the total mass-to-optical luminosity ratio increases as the total mass increases.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

Context: Two-point correlation functions are used throughout cosmology as a measure for the statistics of random fields. When used in Bayesian parameter estimation, their likelihood function is usually replaced by a Gaussian approximation. However, this has been shown to be insufficient. Aims: For the case of Gaussian random fields, we search for an exact probability distribution of correlation functions, which could improve the accuracy of future data analyses. Methods: We use a fully analytic approach, first expanding the random field in its Fourier modes, and then calculating the characteristic function. Finally, we derive the probability distribution function using integration by residues. We use a numerical implementation of the full analytic formula to discuss the behaviour of this function. Results: We derive the univariate and bivariate probability distribution function of the correlation functions of a Gaussian random field, and outline how higher joint distributions could be calculated. We give the results in the form of mode expansions, but in one special case we also find a closed-form expression. We calculate the moments of the distribution and, in the univariate case, we discuss the Edgeworth expansion approximation. We also comment on the difficulties in a fast and exact numerical implementation of our results, and on possible future applications.

To investigate and specify the statistical properties of cosmological fields with particular attention to possible non-Gaussian features, accurate formulae for the bispectrum and the bispectrum covariance are required. The bispectrum is the lowest-order statistic providing an estimate for non-Gaussianities of a distribution, and the bispectrum covariance depicts the errors of the bispectrum measurement and their correlation on different scales. Currently, there do exist fitting formulae for the bispectrum and an analytical expression for the bispectrum covariance, but the former is not very accurate and the latter contains several intricate terms and only one of them can be readily evaluated from the power spectrum of the studied field. Neglecting all higher-order terms results in the Gaussian approximation of the bispectrum covariance. We study the range of validity of this Gaussian approximation for two-dimensional non-Gaussian random fields. For this purpose, we simulate Gaussian and non-Gaussian random fields, the latter represented by log-normal fields and obtained directly from the former by a simple transformation. From the simulated fields, we calculate the power spectra, the bispectra, and the covariance from the sample variance of the bispectra, for different degrees of non-Gaussianity \alpha, which is equivalent to the skewness on a given angular scale \theta g. We find that the Gaussian approximation provides a good approximation for \alpha<0.6 and a reasonably accurate approximation for \alpha< 1, both on scales >8\theta g. Using results from cosmic shear simulations, we estimate that the cosmic shear convergence fields are described by \alpha<0.7 at \theta g~4". We therefore conclude that the Gaussian approximation for the bispectrum covariance is likely to be applicable in ongoing and future cosmic shear studies.

With the growing interest in and ability of using weak lensing studies to probe the non-Gaussian properties of the matter density field, there is an increasing need for the study of suitable statistical measures, e.g. shear three-point statistics. In this paper we establish the relations between the three-point configuration space shear and convergence statistics, which are an important missing link between different weak lensing three-point statistics and provide an alternative way of relating observation and theory. The method we use also allows us to derive the relations between other two- and three-point correlation functions. We show the consistency of the relations obtained with already established results and demonstrate how they can be evaluated numerically. As a direct application, we use these relations to formulate the condition for E/B-mode decomposition of lensing three-point statistics, which is the basis for constructing new three-point statistics which allow for exact E/B-mode separation. Our work applies also to other two-dimensional polarization fields such as that of the Cosmic Microwave Background.

Protostars are actively accreting matter and they drive spectacular, dynamic outflows, which evolve on timescales of years. X-ray emission from these jets has been detected only in a few cases and little is known about its time evolution. We present a new Chandra observation of L1551 IRS 5's jet in the context of all available X-ray data of this object. Specifically, we perform a spatially resolved spectral analysis of the X-ray emission and find that (a) the total X-ray luminosity is constant over almost one decade, (b) the majority of the X-rays appear to be always located close to the driving source, (c) there is a clear trend in the photon energy as a function of the distance to the driving source indicating that the plasma is cooler at larger distances and (d) the X-ray emission is located in a small volume which is unresolved perpendicular to the jet axis by Chandra. A comparison of our X-ray data of the L1551 IRS 5 jet both with models as well as X-ray observations of other protostellar jets shows that a base/standing shock is a likely and plausible explanation for the apparent constancy of the observed X-ray emission. Internal shocks are also consistent with the observed morphology if the supply of jet material by the ejection of new blobs is sufficiently constant. We conclude that the study of the X-ray emission of protostellar jet sources allows us to diagnose the innermost regions close to the acceleration region of the outflows.

Context: Understanding the evolution of the dark matter halos of galaxies after they become part of a cluster is essential for understanding the evolution of these satellite galaxies. Aims: We investigate the potential of galaxy-galaxy lensing to map the halo density profiles of galaxies in clusters. Methods: We propose a method that separates the weak-lensing signal of the dark-matter halos of galaxies in clusters from the weak-lensing signal of the cluster's main halo. Using toy cluster models as well as ray-tracing through N-body simulations of structure formation along with semi-analytic galaxy formation models, we test the method and assess its performance. Results: We show that with the proposed method, one can recover the density profiles of the cluster galaxy halos in the range 30 - 300 kpc. Using the method, we find that weak-lensing signal of cluster member galaxies in the Millennium Simulation is well described by an Navarro-Frenk-White (NFW) profile. In contrast, non-singular isothermal mass distribution (like PIEMD) model provide a poor fit. Furthermore, we do not find evidence for a sharp truncation of the galaxy halos in the range probed by our method. Instead, there is an observed overall decrease of the halo mass profile of cluster member galaxies with increasing time spent in the cluster. This trend, as well as the presence or absence of a truncation radius, should be detectable in future weak-lensing surveys like the Dark Energy Survey (DES) or the Large Synoptic Survey Telescope (LSST) survey. Such surveys should also allow one to infer the mass-luminosity relation of cluster galaxies with our method over two decades in mass. Conclusions: It is possible to recover in a non-parametric way the mass profile of satellite galaxies and their dark matter halos in future surveys, using our proposed weak lensing method.

We identify and study a previously unknown systematic effect on cosmic shear measurements, caused by the selection of galaxies used for shape measurement, in particular the rejection of close (blended) galaxy pairs. We use ray-tracing simulations based on the Millennium Simulation and a semi-analytical model of galaxy formation to create realistic galaxy catalogues. From these, we quantify the bias in the shear correlation functions by comparing measurements made from galaxy catalogues with and without removal of close pairs. A likelihood analysis is used to quantify the resulting shift in estimates of cosmological parameters. The filtering of objects with close neighbours (a) changes the redshift distribution of the galaxies used for correlation function measurements, and (b) correlates the number density of sources in the background with the density field in the foreground. This leads to a scale-dependent bias of the correlation function of several percent, translating into biases of cosmological parameters of similar amplitude. This makes this new systematic effect potentially harmful for upcoming and planned cosmic shear surveys. As a remedy, we propose and test a weighting scheme that can significantly reduce the bias.

We report on new X-ray observations of the classical T Tauri star DG Tau. DG Tau drives a collimated bi-polar jet known to be a source of X-ray emission perhaps driven by internal shocks. The rather modest extinction permits study of the jet system to distances very close to the star itself. Our initial results presented here show that the spatially resolved X-ray jet has been moving and fading during the past six years. In contrast, a stationary, very soft source much closer (~ 0.15-0.2") to the star but apparently also related to the jet has brightened during the same period. We report accurate temperatures and absorption column densities toward this source, which is probably associated with the jet base or the jet collimation region.

Three-point weak lensing statistics provide cosmic information complementary to that of two-point statistics. However, both statistics suffer from intrinsic-shear alignment, which is one of their limiting systematics. The nulling technique is a model-independent method developed to eliminate intrinsic-shear alignment at the two-point level. In this paper we demonstrate that the nulling technique can also be naturally generalized to the three-point level, controlling the corresponding GGI systematics. We show that under the assumption of exact redshift information the intrinsic-shear alignment contamination can be completely eliminated. To show how well the nulling technique performs on data with limited redshift information, we apply the nulling technique to three-point weak lensing statistics from a fictitious survey analogous to a typical future deep imaging survey, in which the three-point intrinsic-shear alignment systematics is generated from a power-law toy model. Using 10 redshift bins, the nulling technique leads to a factor of 10 suppression of the GGI/GGG ratio, and reduces the bias on cosmological parameters to less than the original statistical error. More detailed redshift information allowing for finer redshift bins leads to better reduction of bias. The information loss during the nulling procedure doubles the statistical error on cosmological parameters. A comparison of the nulling technique with an unconditioned compression of the data suggests that part of the information loss can be retained by considering higher-order nulling weights during the nulling procedure. A combined analysis of two- and three-point statistics confirms that the information contained in them is of comparable size and is complementary to each other, both before and after nulling.

Cosmic shear is considered one of the most powerful methods for studying the properties of Dark Energy in the Universe. As a standard method, the two-point correlation functions $xi_\pm(theta)$ of the cosmic shear field are used as statistical measures for the shear field. In order to separate the observed shear into E- and B-modes, the latter being most likely produced by remaining systematics in the data set and/or intrinsic alignment effects, several statistics have been defined before. Here we aim at a complete E-/B-mode decomposition of the cosmic shear information contained in the $xi_\pm$ on a finite angular interval. We construct two sets of such E-/B-mode measures, namely Complete Orthogonal Sets of E-/B-mode Integrals (COSEBIs), characterized by weight functions between the $xi_\pm$ and the COSEBIs which are polynomials in $theta$ or polynomials in $ln(theta)$, respectively. Considering the likelihood in cosmological parameter space, constructed from the COSEBIs, we study their information contents. We show that the information grows with the number of COSEBI modes taken into account, and that an asymptotic limit is reached which defines the maximum available information in the E-mode component of the $xi_\pm$. We show that this limit is reached the earlier (i.e., for a smaller number of modes considered) the narrower the angular range is over which $xi_\pm$ are measured, and it is reached much earlier for logarithmic weight functions. For example, for $xi_\pm$ on the interval $1'\le \theta\le 400'$, the asymptotic limit for the parameter pair $(Omega_m, sigma_8)$ is reached for $\sim 25$ modes in the linear case, but already for 5 modes in the logarithmic case. The COSEBIs form a natural discrete set of quantities, which we suggest as method of choice in future cosmic shear likelihood analyses.

Aims. The predictions of the ellipticity of the dark matter halos from models of structure formation are notoriously difficult to test with observations. A direct measurement would give important constraints on the formation of galaxies, and its effect on the dark matter distribution in their halos. Here we show that galaxy-galaxy flexion provides a direct and potentially powerful method for determining the ellipticity of (an ensemble of) elliptical lenses. Methods. We decompose the spin-1 flexion into a radial and a tangential component. Using the ratio of tangential-to- radial flexion, which is independent of the radial mass profile, the mass ellipticity can be estimated. Results. An estimator for the ellipticity of the mass distribution is derived and tested with simulations. We show that the estimator is slightly biased. We quantify this bias, and provide a method to reduce it. Furthermore, a parametric fitting of the flexion ratio and orientation provides another estimate for the dark halo ellipticity, which is more accurate for individual lenses Overall, galaxy-galaxy flexion appears as a powerful tool for constraining the ellipticity of mass distributions.

The intrinsic alignment of galaxies constitutes the major astrophysical source of systematic errors in surveys of weak gravitational lensing by the large-scale structure. We discuss the principles, summarise the implementation, and highlight the performance of two model-independent methods that control intrinsic alignment signals in weak lensing data: the nulling technique which eliminates intrinsic alignments to ensure unbiased constraints on cosmology, and the boosting technique which extracts intrinsic alignments and hence allows one to further study this contribution. Making only use of the characteristic dependence on redshift of the signals, both approaches are robust, but reduce the statistical power due to the similar redshift scaling of intrinsic alignment and lensing signals.

Intrinsic alignments constitute the major astrophysical systematic for cosmological weak lensing surveys. We present a purely geometrical method with which one can study gravitational shear-intrinsic ellipticity correlations directly in weak lensing data. Linear combinations of second-order cosmic shear measures are constructed such that the intrinsic alignment signal is boosted while suppressing the contribution by gravitational lensing. We then assess the performance of a specific parametrisation of the weights entering these linear combinations for three representative survey models. Moreover a relation between this boosting technique and the intrinsic alignment removal via nulling is derived. For future all-sky weak lensing surveys with photometric redshift information the boosting technique yields statistical errors on model parameters of intrinsic alignments whose order of magnitude is compatible with current constraints determined from indirect measurements. Parameter biases due to a residual cosmic shear signal are negligible in case of quasi-spectroscopic redshifts and remain sub-dominant for typical values of the photometric redshift scatter. We find good agreement between the performance of the intrinsic alignment removal based on the boosting technique and standard nulling methods, possibly indicating a fundamental limit in the separation of lensing and intrinsic alignment signals.

Galaxy clusters are considered as excellent probes for cosmology. For that purpose, their mass needs to be measured and their structural properties needs to be understood. We propose a method for galaxy cluster mass reconstruction which combines information from strong lensing, weak lensing shear and flexion. We extend the weak lensing analysis to the inner parts of the cluster and, in particular, improve the resolution of substructure. We use simulations to show that the method recovers the mass density profiles of the cluster. We find that the weak lensing flexion is sensitive to substructure. After combining the flexion data into the joint weak and strong lensing analysis, we can resolve the cluster properties with substructures.

AU Mic is a young, nearby X-ray active M-dwarf with an edge-on debris disk. Debris disk are the successors of the gaseous disks usually surrounding pre-main sequence stars which form after the first few Myrs of their host stars' lifetime, when - presumably - also the planet formation takes place. Since X-ray transmission spectroscopy is sensitive to the chemical composition of the absorber, features in the stellar spectrum of AU Mic caused by its debris disk can in principle be detected. The upper limits we derive from our high resolution Chandra LETGS X-ray spectroscopy are on the same order as those from UV absorption measurements, consistent with the idea that AU Mic's debris disk possesses an inner hole with only a very low density of sub-micron sized grains or gas.

We present a tomographic cosmological weak lensing analysis of the HST COSMOS Survey. Applying our lensing-optimized data reduction, principal component interpolation for the ACS PSF, and improved modelling of charge-transfer inefficiency, we measure a lensing signal which is consistent with pure gravitational modes and no significant shape systematics. We carefully estimate the statistical uncertainty from simulated COSMOS-like fields obtained from ray-tracing through the Millennium Simulation. We test our pipeline on simulated space-based data, recalibrate non-linear power spectrum corrections using the ray-tracing, employ photometric redshifts to reduce potential contamination by intrinsic galaxy alignments, and marginalize over systematic uncertainties. We find that the lensing signal scales with redshift as expected from General Relativity for a concordance LCDM cosmology, including the full cross-correlations between different redshift bins. For a flat LCDM cosmology, we measure sigma_8(Omega_m/0.3)^0.51=0.75+-0.08 from lensing, in perfect agreement with WMAP-5, yielding joint constraints Omega_m=0.266+0.025-0.023, sigma_8=0.802+0.028-0.029 (all 68% conf.). Dropping the assumption of flatness and using HST Key Project and BBN priors only, we find a negative deceleration parameter q_0 at 94.3% conf. from the tomographic lensing analysis, providing independent evidence for the accelerated expansion of the Universe. For a flat wCDM cosmology and prior w in [-2,0], we obtain w<-0.41 (90% conf.). Our dark energy constraints are still relatively weak solely due to the limited area of COSMOS. However, they provide an important demonstration for the usefulness of tomographic weak lensing measurements from space. (abridged)

The mass function of galaxy clusters at high redshifts is a particularly useful probe to learn about the history of structure formation and constrain cosmological parameters. We aim at deriving reliable masses for a high-redshift, high-luminosity sample of clusters of galaxies selected from the 400d survey of X-ray selected clusters. Here, we will focus on a particular object, CL0030+2618 at z=0.50 Using deep imaging in three passbands with the MEGACAM instrument at MMT, we show that MEGACAM is well-suited for measuring gravitational shear. We detect the weak lensing signal of CL0030+2618 at 5.8 sigma significance, using the aperture mass technique. Furthermore, we find significant tangential alignment of galaxies out to ~10 arcmin or >2r_200 distance from the cluster centre. The weak lensing centre of CL0030+2618 agrees with several X-ray measurements and the position of the brightest cluster galaxy. Finally, we infer a weak lensing virial mass of M_200=7.5 10^{14} M_sun for CL0030+2618. Despite complications by a tentative foreground galaxy group in the line of sight, the X-ray and weak lensing estimates for CL0030+2618 are in remarkable agreement. This study paves the way for the largest weak lensing survey of high-redshift galaxy clusters to date.

To probe cosmological fields beyond the Gaussian level, three-point statistics can be used, all of which are related to the bispectrum. Hence, measurements of CMB anisotropies, galaxy clustering, and weak gravitational lensing alike have to rely upon an accurate theoretical background concerning the bispectrum and its noise properties. If only small portions of the sky are considered, it is often desirable to perform the analysis in the flat-sky limit. We aim at a formal, detailed derivation of the bispectrum covariance in the flat-sky approximation, focusing on a pure two-dimensional Fourier-plane approach. We define an unbiased estimator of the bispectrum, which takes the average over the overlap of annuli in Fourier space, and compute its full covariance. The outcome of our formalism is compared to the flat-sky spherical harmonic approximation in terms of the covariance, the behavior under parity transformations, and the information content. We introduce a geometrical interpretation of the averaging process in the estimator, thus providing an intuitive understanding. Contrary to foregoing work, we find a difference by a factor of two between the covariances of the Fourier-plane and the spherical harmonic approach. We argue that this discrepancy can be explained by the differing behavior with respect to parity. However, in an exemplary analysis it is demonstrated that the Fisher information of both formalisms agrees to high accuracy. Via the geometrical interpretation we are able to link the normalization in the bispectrum estimator to the area enclosed by the triangle configuration at consideration as well as to the Wigner symbol, which leads to convenient approximation formulae for the covariances of both approaches.

Cosmic shear is regarded one of the most powerful probes to reveal the properties of dark matter and dark energy. To fully utilize its potential, one has to be able to control systematic effects down to below the level of the statistical parameter errors. Particularly worrisome in this respect is intrinsic alignment, causing considerable parameter biases via correlations between the intrinsic ellipticities of galaxies and the gravitational shear, which mimic lensing. In an earlier work we have proposed a nulling technique that downweights this systematic, only making use of its well-known redshift dependence. We assess the practicability of nulling, given realistic conditions on photometric redshift information. For several simplified intrinsic alignment models and a wide range of photometric redshift characteristics we calculate an average bias before and after nulling. Modifications of the technique are introduced to optimize the bias removal and minimize the information loss by nulling. We demonstrate that one of the presented versions is close to optimal in terms of bias removal, given high quality of photometric redshifts. For excellent photometric redshift information, i.e. at least 10 bins with a small dispersion, a negligible fraction of catastrophic outliers, and precise knowledge about the redshift distributions, one version of nulling is capable of reducing the shear-intrinsic ellipticity contamination by at least a factor of 100. Alternatively, we describe a robust nulling variant which suppresses the systematic signal by about 10 for a very broad range of photometric redshift configurations. Irrespective of the photometric redshift quality, a loss of statistical power is inherent to nulling, which amounts to a decrease of the order 50% in terms of our figure of merit.

The gravitational potential of clusters of galaxies acts as a cosmic telescope allowing us to find and study galaxies at fainter limits than otherwise possible and thus probe closer to the epoch of formation of the first galaxies. We use the Bullet Cluster 1E0657-56 (z = 0.296) as a case study, because its high mass and merging configuration makes it one of the most efficient cosmic telescopes we know. We develop a new algorithm to reconstruct the gravitational potential of the Bullet Cluster, based on a non-uniform adaptive grid, combining strong and weak gravitational lensing data derived from deep HST/ACS F606W-F775W-F850LP and ground-based imaging. We exploit this improved mass map to study z~5-6 Lyman Break Galaxies (LBGs), which we detect as dropouts. One of the LBGs is multiply imaged, providing a geometric confirmation of its high redshift, and is used to further improve our mass model. We quantify the uncertainties in the magnification map reconstruction in the intrinsic source luminosity, and in the volume surveyed, and show that they are negligible compared to sample variance when determining the luminosity function of high-redshift galaxies. With shallower and comparable magnitude limits to HUDF and GOODS, the Bullet cluster observations, after correcting for magnification, probe deeper into the luminosity function of the high redshift galaxies than GOODS and only slightly shallower than HUDF. We conclude that accurately focused cosmic telescopes are the most efficient way to sample the bright end of the luminosity function of high redshift galaxies and - in case they are multiply imaged - confirm their redshifts.

The central part of the massive star-forming region Cepheus A contains several radio sources which indicate multiple outflow phenomena, yet the driving sources of the individual outflows have not been identified. We present a high-resolution Chandra observation of this region that shows the presence of bright X-ray sources, consistent with active pre-main sequence stars, while the strong absorption hampers the detection of less luminous objects. A new source has been discovered located on the line connecting H_2 emission regions at the eastern and western parts of Cepheus A. This source could be the driving source of HH 168. We present a scenario relating the observed X-ray and radio emission.

X-ray emission from massive stellar outflows has been detected in several cases. We present a Chandra observation of HH 168 and show that the soft X-ray emission from a plasma of 0.55 keV within HH 168 is diffuse. The X-ray emission is observed on two different scales: Three individual, yet extended, regions are embedded within a complex of low X-ray surface brightness. Compared to the bow shock the emission is displaced against the outflow direction. We show that there is no significant contribution from young stellar objects (YSOs) and discuss several shock scenarios that can produce the observed signatures. We establish that the X-ray emission of HH 168 is excited by internal shocks in contrast to simple models, which expect the bow shock to be the most X-ray luminous.

(abridged) We study the validity of the approximation of a Gaussian cosmic shear likelihood. We estimate the true likelihood for a fiducial cosmological model from a large set of ray-tracing simulations and investigate the impact of non-Gaussianity on cosmological parameter estimation. We investigate how odd the recently reported very low value of $\sigma_8$ really is as derived from the \textit{Chandra} Deep Field South (CDFS) using cosmic shear by taking the non-Gaussianity of the likelihood into account as well as the possibility of biases coming from the way the CDFS was selected. We find that the cosmic shear likelihood is significantly non-Gaussian. This leads to both a shift of the maximum of the posterior distribution and a significantly smaller credible region compared to the Gaussian case. We re-analyse the CDFS cosmic shear data using the non-Gaussian likelihood. Assuming that the CDFS is a random pointing, we find $\sigma_8=0.68_{-0.16}^{+0.09}$ for fixed $\Omega_{\rm m}=0.25$. In a WMAP5-like cosmology, a value equal to or lower than this would be expected in $\approx 5%$ of the times. Taking biases into account arising from the way the CDFS was selected, which we model as being dependent on the number of haloes in the CDFS, we obtain $\sigma_8 = 0.71^{+0.10}_{-0.15}$. Combining the CDFS data with the parameter constraints from WMAP5 yields $\Omega_{\rm m} = 0.26^{+0.03}_{-0.02}$ and $\sigma_8 = 0.79^{+0.04}_{-0.03}$ for a flat universe.

Commonly used methods to decompose E- and B-modes in cosmic shear, namely the aperture mass dispersion and the E/B-mode shear correlation function, suffer from incomplete knowledge of the two-point correlation function (2PCF) on very small and/or very large scales. The ring statistics, the most recently developed cosmic shear measure, improves on this issue and is able to decompose E- and B-modes using a 2PCF measured on a finite interval. First, we improve on the ring statistics' filter function with respect to the signal-to-noise ratio. Second, we examine the ability of the ring statistics to constrain cosmology and compare the results to cosmological constraints obtained with the aperture mass dispersion. Third, we use the ring statistics to measure a cosmic shear signal from CFHTLS (Canada-France-Hawaii Telescope Legacy Survey) data. We consider a scale-dependent filter function for the ring statistics which improves its signal-to-noise ratio. In addition, we show that there exist filter functions which decompose E- and B-modes using a finite range of 2PCFs (EB-statistics) and have higher S/N ratio than the ring statistics. However, we find that data points of the latter are significantly less correlated than data points of the aperture mass dispersion and the EB-statistics. As a consequence the ring statistics is an ideal tool to identify remaining systematics accurately as a function of angular scale. We use the 2PCF of the latest CFHTLS analysis and therefrom calculate the ring statistics and its error bars.

Weak gravitational lensing is one of the most promising tools to investigate the equation-of-state of dark energy. In order to obtain reliable parameter estimations for current and future experiments, a good theoretical understanding of dark matter clustering is essential. Of particular interest is the statistical precision to which weak lensing observables, such as cosmic shear correlation functions, can be determined. We construct a fitting formula for the non-Gaussian part of the covariance of the lensing power spectrum. The Gaussian contribution to the covariance, which is proportional to the lensing power spectrum squared, and optionally shape noise can be included easily by adding their contributions. Starting from a canonical estimator for the dimensionless lensing power spectrum, we model first the covariance in the halo model approach including all four halo terms for one fiducial cosmology and then fit two polynomials to the expression found. On large scales, we use a first-order polynomial in the wave-numbers and dimensionless power spectra that goes asymptotically towards $1.1 C_{pt}$ for $\ell \to 0$, i.e., the result for the non-Gaussian part of the covariance using tree-level perturbation theory. On the other hand, for small scales we employ a second-order polynomial in the dimensionless power spectra for the fit. We obtain a fitting formula for the non-Gaussian contribution of the convergence power spectrum covariance that is accurate to 10% for the off-diagonal elements, and to 5% for the diagonal elements, in the range $50 \lesssim \ell \lesssim 5000$ and can be used for single source redshifts $z_{s} \in [0.5,2.0]$ in WMAP5-like cosmologies.

We present an equation of state and the composition of low-density supernova matter composed of light nuclei with mass number A \le 13. We work within the quasiparticle gas model, which accounts for bound states with decay time scales larger than the relevant time scale of supernova and protoneutron star evolution. The mean-field contribution is included in terms of Skyrme density functional. Deuterons, tritons, and 3H(e) nuclei appear in matter in concentrations that are substantially higher than those of heavier nuclei. We calculate the critical temperature of deuteron condensation in such matter, and demonstrate that the appearance of clusters substantially lowers the critical temperature.

We show that correlation functions have to satisfy contraint relations, owing to the non-negativity of the power spectrum of the underlying random process. Specifically, for any statistically homogeneous and (for more than one spatial dimension) isotropic random field with correlation function $\xi(x)$, we derive inequalities for the correlation coefficients $r_n\equiv \xi(n x)/\xi(0)$ (for integer $n$) of the form $r_{n{\rm l}}\le r_n\le r_{n{\rm u}}$, where the lower and upper bounds on $r_n$ depend on the $r_j$, with $j<n$. Explicit expressions for the bounds are obtained for arbitrary $n$. These constraint equations very significantly limit the set of possible correlation functions. For one particular example of a fiducial cosmic shear survey, we show that the Gaussian likelihood ellipsoid has a significant spill-over into the forbidden region of correlation functions, rendering the resulting best-fitting model parameters and their error region questionable, and indicating the need for a better description of the likelihood function. We conduct some simple numerical experiments which explicitly demonstrate the failure of a Gaussian description for the likelihood of $\xi$. Instead, the shape of the likelihood function of the correlation coefficients appears to follow approximately that of the shape of the bounds on the $r_n$, even if the Gaussian ellipsoid lies well within the allowed region. For more than one spatial dimension of the random field, the explicit expressions of the bounds on the $r_n$ are not optimal. We outline a geometrical method how tighter bounds may be obtained in principle. We illustrate this method for a few simple cases; a more general treatment awaits future work.

This paper studies the relative spatial distribution of red-sequence and blue-cloud galaxies, and their relation to the dark matter distribution in the COMBO-17 survey as function of scale down to z~1. We measure the 2nd-order auto- and cross-correlation functions of galaxy clustering and express the relative biasing by using aperture statistics. Also estimated is the relation between the galaxies and the dark matter distribution exploiting galaxy-galaxy lensing (GGL). All observables are further interpreted in terms of a halo model. To fully explain the galaxy clustering cross-correlation function with a halo model, we need to introduce a new parameter,R, that describes the statistical relation between numbers of red and blue galaxies within the same halo. We find that red and blue galaxies are clearly differently clustered, a significant evolution of the relative clustering with redshift was not found. There is evidence for a scale-dependence of relative biasing. The relative clustering, the GGL and, with some tension, the galaxy numbers can be explained consistently within a halo model. For the cross-correlation function one requires a HOD variance that becomes Poisson even for relatively small occupancy numbers. For our sample, this rules out with high confidence a "Poisson satellite" scenario as found in semi-analytical models. Red galaxies have to be concentrated towards the halo centre, either by a central red galaxy or by a concentration parameter above that for dark matter.The value of R depends on the presence or absence of central galaxies: If no central galaxies or only red central galaxies are allowed, R is consistent with zero, whereas a positive correlation $R=+0.5\pm0.2$ is needed if both blue and red galaxies can have central galaxies.[ABRIDGED]

(abridged) We study the accuracy of various approximations to cosmic shear and weak galaxy-galaxy lensing and investigate effects of Born corrections and lens-lens coupling. We use ray-tracing through the Millennium Simulation to calculate various cosmic-shear and galaxy-galaxy-lensing statistics. We compare the results from ray-tracing to semi-analytic predictions. We find: (i) The linear approximation provides an excellent fit to cosmic-shear power spectra as long as the actual matter power spectrum is used as input. Common fitting formulae, however, strongly underestimate the cosmic-shear power spectra. Halo models provide a better fit to cosmic shear-power spectra, but there are still noticeable deviations. (ii) Cosmic-shear B-modes induced by Born corrections and lens-lens coupling are at least three orders of magnitude smaller than cosmic-shear E-modes. Semi-analytic extensions to the linear approximation predict the right order of magnitude for the B-mode. Compared to the ray-tracing results, however, the semi-analytic predictions may differ by a factor two on small scales and also show a different scale dependence. (iii) The linear approximation may under- or overestimate the galaxy-galaxy-lensing shear signal by several percent due to the neglect of magnification bias, which may lead to a correlation between the shear and the observed number density of lenses. We conclude: (i) Current semi-analytic models need to be improved in order to match the degree of statistical accuracy expected for future weak-lensing surveys. (ii) Shear B-modes induced by corrections to the linear approximation are not important for future cosmic-shear surveys. (iii) Magnification bias can be important for galaxy-galaxy-lensing surveys.

We present the CFHTLS-Archive-Research Survey (CARS). It is a virtual multi-colour survey based on public archive images from the CFHT-Legacy-Survey. Our main scientific interests in CARS are optical searches for galaxy clusters from low to high redshift and their subsequent study with photometric and weak-gravitational lensing techniques. As a first step of the project we present multi-colour catalogues from 37 sq. degrees of the CFHTLS-Wide component. Our aims are to create astrometrically and photometrically well calibrated co-added images. Second goal are five-band (u*, g', r', i', z') multi-band catalogues with an emphasis on reliable estimates for object colours. These are subsequently used for photometric redshift estimates. The article explains in detail data processing, multi-colour catalogue creation and photometric redshift estimation. Furthermore we apply a novel technique, based on studies of the angular galaxy cross-correlation function, to quantify the reliability of photo-z's. The accuracy of our high-confidence photo-z sample (10-15 galaxies per sq. arcmin) is estimated to $\sigma_{\Delta_z/(1+z)}\approx 0.04-0.05$ up to i'<24 with typically only 1-3% outliers. Interested users can obtain access to our data by request to the authors.

In cosmic shear likelihood analyses the covariance is most commonly assumed to be constant in parameter space. Therefore, when calculating the covariance matrix (analytically or from simulations), its underlying cosmology should not influence the likelihood contours. We examine whether the aforementioned assumption holds and quantify how strong cosmic shear covariances vary within a reasonable parameter range. Furthermore, we examine the impact on likelihood contours when assuming different cosmologies in the covariance. We find that covariances vary significantly within the considered parameter range (Omega_m=[0.2;0.4], sigma_8=[0.6;1.0]) and that this has a non-negligible impact on the size of likelihood contours. This impact increases with increasing survey size, increasing number density of source galaxies, decreasing ellipticity noise, and when using non-Gaussian covariances. To improve on the assumption of a constant covariance we present two methods. The adaptive covariance is the most accurate method, but it is computationally expensive. To reduce the computational costs we give a scaling relation for covariances. As a second method we outline the concept of an iterative likelihood analysis. Here, we additionally account for non-Gaussianity using a ray-tracing covariance derived from the Millennium simulation.

To render cosmic shear an astronomical tool of high precision, it is essential to eliminate systematic effects upon its signal, one of the most significant ones being correlations between the gravitational shear and the intrinsic ellipticity of source galaxies. Because of the crudeness of current models of intrinsic alignment, we have developed a model-independent, purely geometrical method to remove shear-ellipticity correlations. We eliminate the contributions to a tomographic cosmic shear signal that may be subject to contamination by shear-ellipticity correlations, making use of the characteristic dependence of these correlations on redshift. By introducing an appropriately chosen weight function to the lensing efficiency that nulls signals stemming from certain distances, new second-order measures of cosmic shear, free of intrinsic alignment, can be constructed. The loss of information induced by this nulling technique and the subsequent degradation of constraints on cosmological parameters is quantified in a likelihood analysis. Concerning the construction of optimal weight functions, good agreement is achieved between all approaches considered. For a survey divided into 20 redshift bins we find that the area of credible regions increases by 20% up to about 50% after the application of nulling, depending on the cosmological parameters considered. We demonstrate that due to an optimization of the weight functions nearly all information is contained in a small subset of the new second-order measures. The use of a significantly smaller number of redshift bins than 20 for the nulling considerably degrades parameter constraints under conservative assumptions, emphasizing the need for detailed redshift information.

The eclipsing brown-dwarf binary system 2MASS J05352184-0546085 is a case sui generis. For the first time, it allows a detailed analysis of the individual properties of young brown dwarfs, in particular, masses, and radii, and the temperature ratio of the system components can be determined accurately. The system shows a "temperature reversal" with the more massive component being the cooler one, and both components are found to be active. We analyze X-ray images obtained by Chandra and XMM-Newton containing 2MASS J05352184-0546085 in their respective field of view. The Chandra observatory data show a clear X-ray source at the position of 2MASS J05352184-0546085, whereas the XMM-Newton data suffer from contamination from other nearby sources, but are consistent with the Chandra detection. No indications of flaring activity are found in either of the observations (together about 70 ks), and we thus attribute the observed flux to quiescent emission. With an X-ray luminosity of 3*10^{28} erg/s we find an L_X/L_{bol}-ratio close to the saturation limit of 10^{-3} and an L_{X}/L_{H\alpha}-ratio consistent with values obtained from low-mass stars. The X-ray detection of 2MASS J05352184-0546085 reported here provides additional support for the interpretation of the temperature reversal in terms of magnetically suppressed convection, and suggests that the activity phenomena of young brown dwarfs resemble those of their more massive counterparts.

The classical T Tauri star DG Tau shows all typical signatures of X-ray activity and, in particular, harbors a resolved X-ray jet. We demonstrate that its soft and hard X-ray components are separated spatially by approximately 0.2 arcsec by deriving the spatial offset between both components from the event centroids of the soft and hard photons utilizing the intrinsic energy-resolution of the Chandra ACIS-S detector. We also demonstrate that this offset is physical and cannot be attributed to an instrumental origin or to low counting statistics. Furthermore, the location of the derived soft X-ray emission peak coincides with emission peaks observed for optical emission lines, suggesting that both, soft X-rays and optical emission, have the same physical origin.

We investigate how strong gravitational lensing in the concordance LCDM cosmology is affected by the stellar mass in galaxies. We extend our previous studies, based on ray-tracing through the Millennium Simulation, by including the stellar components predicted by galaxy formation models. We find that the inclusion of these components greatly enhances the probability for strong lensing compared to a `dark matter only' universe. The identification of the `lenses' associated with strong-lensing events reveals that the stellar mass of galaxies (i) significantly enhances the strong-lensing cross-sections of group and cluster halos, and (ii) gives rise to strong lensing in smaller halos, which would not produce noticeable effects in the absence of the stars. Even if we consider only image splittings >10 arcsec, the luminous matter can enhance the strong-lensing optical depths by up to a factor of 2.

The galaxy cluster RX J1347-1145 is one of the most X-ray luminous and most massive clusters known. Its extreme mass makes it a prime target for studying issues addressing cluster formation and cosmology. In this paper we present new high-resolution HST/ACS and Chandra X-ray data. The high resolution and sensitivity of ACS enabled us to detect and quantify several new multiply imaged sources, we now use a total of eight for the strong lensing analysis. Combining this information with shape measurements of weak lensing sources in the central regions of the cluster, we derive a high-resolution, absolutely-calibrated mass map. This map provides the best available quantification of the total mass of the central part of the cluster to date. We compare the reconstructed mass with that inferred from the new Chandra X-ray data, and conclude that both mass estimates agree extremely well in the observed region, namely within 400 / h_70 kpc of the cluster center. In addition we study the major baryonic components (gas and stars) and hence derive the dark matter distribution in the center of the cluster. We find that the dark matter and baryons are both centered on the BCG within the uncertainties (alignment is better than <10 kpc). We measure the corresponding 1-D profiles and find that dark matter distribution is consistent with both NFW and cored profiles, indicating that a more extended radial analysis is needed to pinpoint the concentration parameter, and hence the inner slope of the dark matter profile.

In recent years cosmic shear, the weak gravitational lensing effect by the large-scale structure of the Universe, has proven to be one of the observational pillars on which the cosmological concordance model is founded. Several cosmic shear statistics have been developed in order to analyze data from surveys. For the covariances of the prevalent second-order measures we present simple and handy formulae, valid under the assumptions of Gaussian density fluctuations and a simple survey geometry. We also formulate these results in the context of shear tomography, i.e. the inclusion of redshift information, and generalize them to arbitrary data field geometries. We define estimators for the E- and B-mode projected power spectra and show them to be unbiased in the case of Gaussianity and a simple survey geometry. From the covariance of these estimators we demonstrate how to derive covariances of arbitrary combinations of second-order cosmic shear measures. We then recalculate the power spectrum covariance for general survey geometries and examine the bias thereby introduced on the estimators for exemplary configurations. Our results for the covariances are considerably simpler than and analytically shown to be equivalent to the real-space approach presented in the first paper of this series. We find good agreement with other numerical evaluations and confirm the general properties of the covariance matrices. The studies of the specific survey configurations suggest that our simplified covariances may be employed for realistic survey geometries to good approximation.

We investigate the impact of the observed correlation between a galaxies shape and its surrounding density field on the measurement of third order weak lensing shear statistics. Using numerical simulations, we estimate the systematic error contribution to a measurement of the third order moment of the aperture mass statistic (GGG) from three-point intrinsic ellipticity correlations (III), and the three-point coupling between the weak lensing shear experienced by distant galaxies and the shape of foreground galaxies (GGI and GII). We find that third-order weak lensing statistics are typically more strongly contaminated by these physical systematics compared to second-order shear measurements, contaminating the measured three-point signal for moderately deep surveys with a median redshift z_m ~ 0.7 by ~ 15%. It has been shown that accurate photometric redshifts will be crucial to correct for this effect, once a model and the redshift dependence of the effect can be accurately constrained. To this end we provide redshift-dependent fitting functions to our results and propose a new tool for the observational study of intrinsic galaxy alignments. For a shallow survey with z_m ~ 0.4 we find III to be an order of magnitude larger than the expected cosmological GGG shear signal. Compared to the two-point intrinsic ellipticity correlation which is similar in amplitude to the two-point shear signal at these survey depths, third order statistics therefore offer a promising new way to constrain models of intrinsic galaxy alignments. Early shallow data from the next generation of very wide weak lensing surveys will be optimal for this type of study.

High resolution HST/ACS images of the galaxy cluster RX J1347-1145 have enabled us to identify several new multiple image candidates in the cluster, including a 5 image system with a central image. The multiple images allow us to construct an accurate 2-dimensional mass map of the central part of the cluster. The modelling of the cluster mass includes the most prominent cluster galaxies modelled as truncated isothermal spheres and a smooth halo component that is described with 2 parametric profiles. The mass reconstruction is done using a Markov chain Monte Carlo method that provides us with a total projected mass density as well as estimates for the parameters of interest and their respective errors. The mass profile is in reasonable agreement with previous mass estimates based on the X-ray emission from the hot intra-cluster gas, however the X-ray mass estimates are systematically lower than what we obtain with gravitational lensing.

The weak gravitational lensing effect, small coherent distortions of galaxy images by means of a gravitational tidal field, can be used to study the relation between the matter and galaxy distribution. In this context, weak lensing has so far only been used for considering a second-order correlation function that relates the matter density and galaxy number density as a function of separation. We implement two new, third-order correlation functions that have recently been suggested in the literature, and apply them to the Red-Sequence Cluster Survey. We demonstrate that it is possible, even with already existing data, to make significant measurements of third-order lensing correlations. We develop an optimised computer code for the correlation functions. To test its reliability a set of tests are performed. The correlation functions are transformed to aperture statistics, which allow easy tests for remaining systematics in the data. In order to further verify the robustness of our measurement, the signal is shown to vanish when randomising the source ellipticities. Finally, the lensing signal is compared to crude predictions based on the halo-model. On angular scales between roughly 1 arcmin and 11 arcmin a significant third-order correlation between two lens positions and one source ellipticity is found. We discuss this correlation function as a novel tool to study the average matter environment of pairs of galaxies. Correlating two source ellipticities and one lens position yields a less significant but nevertheless detectable signal on a scale of 4 arcmin. Both signals lie roughly within the range expected by theory which supports their cosmological origin.[ABRIDGED]

We investigate strong gravitational lensing in the concordance $\Lambda$CDM cosmology by carrying out ray-tracing along past light cones through the Millennium Simulation, the largest simulation of cosmic structure formation ever carried out. We extend previous ray-tracing methods in order to take full advantage of the large volume and the excellent spatial and mass resolution of the simulation. As a function of source redshift we evaluate the probability that an image will be highly magnified, will be highly elongated or will be one of a set of multiple images. We show that such strong lensing events can almost always be traced to a single dominant lensing object and we study the mass and redshift distribution of these primary lenses. We fit analytic models to the simulated dark halos in order to study how our optical depth measurements are affected by the limited resolution of the simulation and of the lensing planes that we construct from it. We conclude that such effects lead us to underestimate total strong-lensing cross sections by about 15 percent. This is smaller than the effects expected from our neglect of the baryonic components of galaxies. Finally we investigate whether strong lensing is enhanced by material in front of or behind the primary lens. Although strong lensing lines-of-sight are indeed biased towards higher than average mean densities, this additional matter typically contributes only a few percent of the total surface density.

In weak gravitational lensing, the image distortion caused by shear measures the projected tidal gravitational field of the deflecting mass distribution. To lowest order, the shear is proportional to the mean image ellipticity. If the image sizes are not small compared to the scale over which the shear varies, higher-order distortions occur, called flexion. For ordinary weak lensing, the observable quantity is not the shear, but the reduced shear, owing to the mass-sheet degeneracy. Likewise, the flexion itself is unobservable. Rather, higher-order image distortions measure the reduced flexion, i.e., derivatives of the reduced shear. We derive the corresponding lens equation in terms of the reduced flexion and calculate the resulting relation between brightness moments of source and image. Assuming an isotropic distribution of source orientations, estimates for the reduced shear and flexion are obtained; these are then tested with simulations. In particular, the presence of flexion affects the determination of the reduced shear. The results of these simulations yield the amount of bias of the estimators, as a function of the shear and flexion. We point out and quantify a fundamental limitation of the flexion formalism, in terms of the product of reduced flexion and source size. If this product increases above the derived threshold, multiple images of the source are formed locally, and the formalism breaks down. Finally, we show how a general (reduced) flexion field can be decomposed into its four components: two of them are due to a shear field, carrying an E- and B-mode in general. The other two components do not correspond to a shear field; they can also be split up into corresponding E- and B-modes.

We introduce an optimized data vector of cosmic shear measures (N). This data vector has high information content, is not sensitive against B-mode contamination and only shows small correlation between data points of different angular scales. We show that a data vector of the two-point correlation function (2PCF) in general contains more information on cosmological parameters compared to a data vector of the aperture mass dispersion. Reason for this is the fact that <M_ap^2> lacks the information of the convergence power spectrum (P_\kappa) on large angular scales, which is contained in the 2PCF data vector. Therefore we create a combined data vector N, which retains the advantages of <M_ap^2> and in addition is also sensitive to the large-scale information of P_\kappa. We compare the information content of the three data vectors by performing a detailed likelihood analysis and use ray-tracing simulations to derive the covariance matrices. In the last part of the paper we contaminate all data vectors with B-modes on small angular scales and examine their robustness against this contamination.The combined data vector strongly improves constraints on cosmological parameters compared to <M_ap^2>. Although, in case of a pure E-mode signal the information content of the 2PCF is higher, in the more realistic case where B-modes are present the 2PCF data vector is strongly contaminated and yields biased cosmological parameter estimates. N shows to be robust against this contamination. Furthermore the individual data points of N show a much smaller correlation compared to the 2PCF leading to an almost diagonal covariance matrix.

The mass function of galaxy clusters is an important cosmological probe. Differences in the selection method could potentially lead to biases when determining the mass function. From the optical and X-ray data of the XMM-Newton Follow-Up Survey, we obtained a sample of galaxy cluster candidates using weak gravitational lensing, the optical Postman matched filter method, and a search for extended X-ray sources. We developed our weak-lensing search criteria by testing the performance of the aperture mass statistic on realistic ray-tracing simulations matching our survey parameters and by comparing two filter functions. We find that the dominant noise source for our survey is shape noise at almost all significance levels and that spurious cluster detections due to projections of large-scale structures are negligible, except possibly for highly significantly detected peaks. Our full cluster catalog has 155 cluster candidates, 116 found with the Postman matched filter, 59 extended X-ray sources, and 31 shear selected potential clusters. Most of these cluster candidates were not previously known. The present catalog will be a solid foundation for studying possible selection effects in either method.

This is the first paper of a series describing our measurement of weak lensing by large-scale structure using archival observations from the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST). In this work we present results from a pilot study testing the capabilities of the ACS for cosmic shear measurements with early parallel observations and presenting a re-analysis of HST/ACS data from the GEMS survey and the GOODS observations of the Chandra Deep Field South (CDFS). We describe our new correction scheme for the time-dependent ACS PSF based on observations of stellar fields. This is currently the only technique which takes the full time variation of the PSF between individual ACS exposures into account. We estimate that our PSF correction scheme reduces the systematic contribution to the shear correlation functions due to PSF distortions to < 2*10^{-6} for galaxy fields containing at least 10 stars. We perform a number of diagnostic tests indicating that the remaining level of systematics is consistent with zero for the GEMS and GOODS data confirming the success of our PSF correction scheme. For the parallel data we detect a low level of remaining systematics which we interpret to be caused by a lack of sufficient dithering of the data. Combining the shear estimate of the GEMS and GOODS observations using 96 galaxies arcmin^{-2} with the photometric redshift catalogue of the GOODS-MUSIC sample, we determine a local single field estimate for the mass power spectrum normalisation sigma_{8,CDFS}=0.52^{+0.11}_{-0.15} (stat) +/- 0.07 (sys) (68% confidence assuming Gaussian cosmic variance) at fixed Omega_m=0.3 for a LambdaCDM cosmology. We interpret this exceptionally low estimate to be due to a local under-density of the foreground structures in the CDFS.

We introduce our publicly available Wide-Field-Imaging reduction pipeline THELI. The procedures applied for the efficient pre-reduction and astrometric calibration are presented. A special emphasis is put on the methods applied to the photometric calibration. As a test case the reduction of optical data from the ESO Deep Public Survey including the WFI-GOODS data is described. The end-products of this project are now available via the ESO archive Advanced Data Products section.

AIMS. The maximum-likelihood method is the standard approach to obtain model fits to observational data and the corresponding confidence regions. We investigate possible sources of bias in the log-likelihood function and its subsequent analysis, focusing on estimators of the inverse covariance matrix. Furthermore, we study under which circumstances the estimated covariance matrix is invertible. METHODS. We perform Monte-Carlo simulations to investigate the behaviour of estimators for the inverse covariance matrix, depending on the number of independent data sets and the number of variables of the data vectors. RESULTS. We find that the inverse of the maximum-likelihood estimator of the covariance is biased, the amount of bias depending on the ratio of the number of bins (data vector variables), P, to the number of data sets, N. This bias inevitably leads to an -- in extreme cases catastrophic -- underestimation of the size of confidence regions. We report on a method to remove this bias for the idealised case of Gaussian noise and statistically independent data vectors. Moreover, we demonstrate that marginalisation over parameters introduces a bias into the marginalised log-likelihood function. Measures of the sizes of confidence regions suffer from the same problem. Furthermore, we give an analytic proof for the fact that the estimated covariance matrix is singular if P>N.

Aims. The clustering properties of a large sample of U-dropouts are investigated and compared to very precise results for B-dropouts from other studies to identify a possible evolution from z=4 to z=3. Methods. A population of ~8800 candidates for star-forming galaxies at z=3 is selected via the well-known Lyman-break technique from a large optical multicolour survey (the ESO Deep Public Survey). The selection efficiency, contamination rate, and redshift distribution of this population are investigated by means of extensive simulations. Photometric redshifts are estimated for every Lyman-break galaxy (LBG) candidate from its UBVRI photometry yielding an empirical redshift distribution. The measured angular correlation function is deprojected and the resulting spatial correlation lengths and slopes of the correlation function of different subsamples are compared to previous studies. Results. By fitting a simple power law to the correlation function we do not see an evolution in the correlation length and the slope from other studies at z=4 to our study at z=3. In particular, the dependence of the slope on UV-luminosity similar to that recently detected for a sample of B-dropouts is confirmed also for our U-dropouts. For the first time number statistics for U-dropouts are sufficient to clearly detect a departure from a pure power law on small scales down to ~2" reported by other groups for B-dropouts.

ESO and ESA agreed to establish a number of Working Groups to explore possible synergies between these two major European astronomical institutions. This Working Group's mandate was to concentrate on fundamental questions in cosmology, and the scope for tackling these in Europe over the next ~15 years. One major resulting recommendation concerns the provision of new generations of imaging survey, where the image quality and near-IR sensitivity that can be attained only in space are naturally matched by ground-based imaging and spectroscopy to yield massive datasets with well-understood photometric redshifts (photo-z's). Such information is essential for a range of new cosmological tests using gravitational lensing, large-scale structure, clusters of galaxies, and supernovae. Great scope in future cosmology also exists for ELT studies of the intergalactic medium and space-based studies of the CMB and gravitational waves; here the synergy is less direct, but these areas will remain of the highest mutual interest to the agencies. All these recommended facilities will produce vast datasets of general applicability, which will have a tremendous impact on broad areas of astronomy.

Aims. One of the main probes for systematic errors in the cosmic shear signal are the division of the shear field into E- and B-mode shear, where gravitational lensing only produces the former. As shown in a recent note, all currently used E-/B-mode separation methods for the shear correlation functions xi_pm require them to be measured to arbitrarily small and/or large separations which is of course not feasible in practice. Methods. We derive second-order shear statistics which provide a clean separation into E- and B-modes from measurements of xi_pm(theta) over a finite interval only. We call these new statistics the circle and ring statistics, respectively; the latter is obtained by an integral over the former. The mathematical properties of these new shear statistics are obtained, as well as specific expressions for applying them to observed data. Results. It is shown that an E-/B-mode separation can be performed on measurements of xi_pm over a finite interval in angular separation, using the ring statistics. We furthermore generalize this result to derive the most general class of second-order shear statistics which provide a separation of E- and B-mode shear on a given angular interval theta_min <= theta <= theta_max. Our results will be of practical use particularly for future cosmic shear surveys where highly precise measurements of the shear will become available and where control of systematics will be mandatory.

[ABRIDGED] The weak gravitational lensing effect is used to infer matter density fluctuations within the field-of-view of the Garching-Bonn Deep Survey (GaBoDS). This information is employed for a statistical comparison of the galaxy distribution to the total matter distribution. The result of this comparison is expressed by means of the linear bias factor, b, the ratio of density fluctuations, and the correlation factor $r$ between density fluctuations. The total galaxy sample is divided into three sub-samples using R-band magnitudes and the weak lensing analysis is applied separately for each sub-sample. Together with the photometric redshifts from the related COMBO-17 survey we estimate the typical mean redshifts of these samples with $\bar{z}=0.35, 0.47, 0.61$, respectively. For all three samples, a slight galaxy anti-bias, b~0.8+-0.1, on scales of a few Mpc/h is found; the bias factor shows evidence for a slight scale-dependence. The correlation between galaxy and (dark) matter distribution is high, r~0.6+-0.2, indicating a non-linear or/and stochastic biasing relation between matter and galaxies. Between the three samples no significant evolution with redshift is found.

Aims. We quantify the mixing of the measured cosmic-shear E- and B-modes caused by the lack of shear-correlation measurements on small and large scales, arising from a lack of close projected galaxy pairs and the finite field size, respectively. Methods. We calculate the aperture-mass statistics <M_{ap, perp}^2> and the E-/B-mode shear-correlation functions xi_{E, B +/-} where small- and large-scale cutoffs are taken into account. We assess the deviation of the obtained E-mode to the true E-mode and the introduction of a spurious B-mode. Results. The measured aperture-mass dispersion is underestimated by more than 10% on scales smaller than 12 times the lower cutoff. For a precise measurement of the E- and B-modes at the percent level using a combination of xi_{E, B +} and xi_{E, B -}, a field as large as 7 (2.4) degrees is necessary for ground-based (space-based) observations.

The aim of the present work is the construction of a mass-selected galaxy cluster sample based on weak gravitational lensing methods. This sample will be subject to spectroscopic follow-up observations. We apply the mass aperture statistics and a derivative of it to 19 square degrees of high quality, single colour wide field imaging data obtained with the WFI@MPG/ESO 2.2m telescope. For the statistics a family of filter functions is used that approximates the expected tangential radial shear profile and thus allows for the efficient detection of mass concentrations. We identify 158 possible mass concentrations. This is the first time that such a large and blindly selected sample is published. 72 of the detections are associated with concentrations of bright galaxies. For about 22 of those we found spectra in the literature, indicating or proving that the galaxies seen are indeed spatially concentrated. 15 of those were previously known to be clusters or have meanwhile been secured as such. We currently follow-up a larger number of them spectroscopically to obtain deeper insight into their physical properties. The remaining 55% of the possible mass concentrations found are not associated with any optical light, or could not be classified unambiguously. We show that those "dark" detections are to a significant degree due to noise, and appear preferentially in shallow data.

Aims. We present a cosmic shear analysis and data validation of 15 square degree high-quality R-band data of the Garching-Bonn Deep Survey obtained with the Wide Field Imager of the MPG/ESO 2.2m telescope. Methods. We measure the two-point shear correlation functions to calculate the aperture mass dispersion. Both statistics are used to perform the data quality control. Combining the cosmic shear signal with a photometric redshift distribution of a galaxy sub-sample obtained from two square degree of UBVRI-band observations of the Deep Public Survey we determine constraints for the matter density Omega_m, the mass power spectrum normalisation sigma_8 and the dark energy density Omega_Lambda in the magnitude interval R in [21.5,24.5]. In this magnitude interval the effective number density of source galaxies is n=12.5/sq. arcmin, and their mean redshift is z_m=0.78. To estimate the posterior likelihood we employ the Monte Carlo Markov Chain method. Results. Using the aperture mass dispersion we obtain for the mass power spectrum normalisation sigma_8=0.80 +- 0.10 (1 sigma statistical error) at a fixed matter density Omega_m=0.30 assuming a flat universe with negligible baryon content and marginalising over the Hubble parameter and the uncertainties in the fitted redshift distribution.

We outline a novel satellite mission concept, DEMON, aimed at advancing our comprehension of both dark matter and dark energy, taking full advantage of two complementary methods: weak lensing and the statistics of galaxy clusters. We intend to carry out a 5000 sqdeg combined IR, optical and X-ray survey with galaxies up to a redshift of z~2 in order to determine the shear correlation function. We will also find ~100000 galaxy clusters, making it the largest survey of this type to date. The DEMON spacecraft will comprise one IR/optical and eight X-ray telescopes, coupled to multiple cameras operating at different frequency bands. To a great extent, the technology employed has already been partially tested on ongoing missions, therefore ensuring improved reliability.

Aims. In this paper the optical data of the ESO Deep-Public-Survey observed with the Wide Field Imager and reduced with the THELI pipeline are described. Methods. Here we present 63 fully reduced and stacked images. The astrometric and photometric calibrations are discussed and the properties of the images are compared to images released by the ESO Imaging Survey team covering a subset of our data. Results. These images are publicly released to the community. Our main scientific goals with this survey are to study the high-redshift universe by optically pre-selecting high-redshift objects from imaging data and to use VLT instruments for follow-up spectroscopy as well as weak lensing applications.

We present weak lensing mass reconstructions for the 20 high-redshift clusters i n the ESO Distant Cluster Survey. The weak lensing analysis was performed on deep, 3-color optical images taken with VLT/FORS2, using a composite galaxy catalog with separate shape estimators measured in each passband. We find that the EDisCS sample is composed primarily of clusters that are less massive than t hose in current X-ray selected samples at similar redshifts, but that all of the fields are likely to contain massive clusters rather than superpositions of low mass groups. We find that 7 of the 20 fields have additional massive structures which are not associated with the clusters and which can affect the weak lensing mass determination. We compare the mass measurements of the remaining 13 clusters with luminosity measurements from cluster galaxies selected using photometric redshifts and find evidence of a dependence of the cluster mass-to-light ratio with redshift. Finally we determine the noise level in the shear measurements for the fields as a function of exposure time and seeing and demonstrate that future ground-based surveys which plan to perform deep optical imaging for use in weak lensing measurements must achieve point-spread functions smaller than a median of 0.6" FWHM.

As a short introduction to the astronomy session, the response of the community to the Call for Themes issued by ESA and the specific themes selected by the Astronomy Working Group are briefly presented in connection with the four grand themes finally selected for the ESA Science Programme. They are placed in the context of the main discoveries of the past decade and the astronomy projects currently in their development or definition phase. Finally, possible strategies for their implementation are summarised.

Application of the aperture mass (Map-) statistics provides a weak lensing method for the detection of cluster-sized dark matter halos. We present a new aperture filter function and maximise the effectiveness of the Map-statistics to detect cluster-sized halos using analytical models. We then use weak lensing mock catalogues generated from ray-tracing through N-body simulations, to analyse the effect of image treatment on the expected number density of halos. Using the Map-statistics, the aperture radius is typically several arcminutes, hence the aperture often lies partly outside a data field, consequently the signal-to-noise ratio of a halo detection decreases. We study these border effects analytically and by using mock catalogues. We find that the expected number density of halos decreases by a factor of two if the size of a field is comparable to the diameter of the aperture used. We finally report on the results of a weak lensing cluster search applying the Map-statistics to 50 randomly selected fields which were observed with FORS1 at the VLT. Altogether the 50 VLT fields cover an area of 0.64 square degrees. The I-band images were taken under excellent seeing conditions (average seeing 0.6 arcsec.) which results in a high number density of galaxies used for the weak lensing analysis (26/sq.arcmin). In five of the VLT fields, we detect a significant Map-signal which coincides with an overdensity of the light distribution. These detections are thus excellent candidates for shear-selected clusters.

This review forms the Weak Lensing part of the Saas-Fee Advanced Course on Gravitational Lensing. It describes the basicsm applications and results of weak lensing. Contents: (1) Introduction (2) The principles of weak gravitational lensing (3) Observational issues and challenges (4) Clusters of galaxies: Introduction, and strong lensing (5) Mass reconstructions from weak lensing (6) Cosmic shear -- lensing by the LSS (7) Large-scale structure lensing: results (8) The mass of, and associated with galaxies (9) Additional issues in cosmic shear (10) Concluding remarks.

We present the ESO Distant Cluster Survey (EDisCS) a survey of 20 fields containing distant galaxy clusters with redshifts ranging from 0.4 to almost 1.0. Candidate clusters were chosen from among the brightest objects identified in the Las Campanas Distant Cluster Survey, half with estimated redshift z_est~0.5 and half with z_est~0.8. They were confirmed by identifying red sequences in moderately deep two colour data from VLT/FORS2. For confirmed candidates we have assembled deep three-band optical photometry using VLT/FORS2, deep near-infrared photometry in one or two bands using NTT/SOFI, deep optical spectroscopy using VLT/FORS2, wide field imaging in two or three bands using the ESO Wide Field Imager, and HST/ACS mosaic images for 10 of the most distant clusters. This first paper presents our sample and the VLT photometry we have obtained. We present images, colour-magnitude diagrams and richness estimates for our clusters, as well as giving redshifts and positions for the brightest cluster members. Subsequent papers will present our infrared photometry, spectroscopy, HST and wide-field imaging, as well as a wealth of further analysis and science results. Our reduced data become publicly available as these papers are accepted.

We present predictions for cosmological parameter constraints from combined measurements of second- and third-order statistics of cosmic shear. We define the generalized third-order aperture mass statistics <M_ap^3> and show that it contains much more information about the bispectrum of the projected matter density than the skewness of the aperture mass. From theoretical models as well as from LCDM ray-tracing simulations, we calculate <M_ap^2> and <M_ap^3> and their dependence on cosmological parameters. The covariances including shot noise and cosmic variance of M_ap^2, M_ap^3 and their cross-correlation are calculated using ray-tracing simulations. We perform an extensive Fisher matrix analysis, and for various combinations of cosmological parameters, we predict 1-sigma-errors corresponding to measurements from a deep 29 square degree cosmic shear survey. Although the parameter degeneracies can not be lifted completely, the (linear) combination of second- and third-order aperture mass statistics reduces the errors significantly. The strong degeneracy between Omega_m and sigma_8, present for all second-order cosmic shear measures, is diminished substantially, whereas less improvement is found for the near-degenerate pair consisting of the shape parameter Gamma and the spectral index n_s. Uncertainties in the source galaxy redshift z_0 increase the errors of all other parameters.

We present first results of our search for high-redshift galaxies in deep CCD mosaic images. As a pilot study for a larger survey, very deep images of the Chandra Deep Field South (CDFS), taken withWFI@MPG/ESO2.2m, are used to select large samples of 1070 U-band and 565 B-band dropouts with the Lyman-break method. The data of these Lyman-break galaxies are made public as an electronic table. These objects are good candidates for galaxies at z~3 and z~4 which is supported by their photometric redshifts. The distributions of apparent magnitudes and the clustering properties of the two populations are analysed, and they show good agreement to earlier studies. We see no evolution in the comoving clustering scale length from z~3 to z~4. The techniques presented here will be applied to a much larger sample of U-dropouts from the whole survey in near future.

We present a weak lensing analysis of the double cluster system Abell 222 and Abell 223. The lensing reconstruction shows evidence for a possible dark matter filament connecting both clusters. The case for a filamentary connection between A 222/223 is supported by an analysis of the galaxy density and X-ray emission between the clusters. Using the results of N-body simulations, we try to develop a criterion that separates this system into cluster and filament regions. The aim is to find a technique that allows the quantification of the significance of (weak lensing) filament candidates in close pairs of clusters. While this mostly fails, the aperture quadrupole statistics (Schneider & Bartelmann 1997) shows some promise in this area. The cluster masses determined from weak lensing in this system are considerably lower than those previously determined from spectroscopic and X-ray observations (Dietrich et al. 2002; Proust et al. 2000; Wang et al. 1999). Additionally, we report the serendipitous weak lensing detection of a previously unknown cluster in the field of this double cluster system.

We have shown that the cluster-mass reconstruction method which combines strong and weak gravitational lensing data, developed in the first paper in the series, successfully reconstructs the mass distribution of a simulated cluster. In this paper we apply the method to the ground-based high-quality multi-colour data of RX J1347.5-1145, the most X-ray luminous cluster to date. A new analysis of the cluster core on very deep, multi-colour data analysis of VLT/FORS data reveals many more arc candidates than previously known for this cluster. The combined strong and weak lensing reconstruction confirms that the cluster is indeed very massive. If the redshift and identification of the multiple-image system as well as the redshift estimates of the source galaxies used for weak lensing are correct, we determine the enclosed cluster mass in a cylinder to M(<360 h^-1 kpc)= (1.2 +/- 0.3) 10^15 Msun. In addition the reconstructed mass distribution follows the distribution found with independent methods (X-ray measurements, SZ). With higher resolution (e.g. HST imaging data) more reliable multiple imaging information can be obtained and the reconstruction can be improved to accuracies greater than what is currently possible with weak and strong lensing techniques.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, the mass-sheet degeneracy transformation has limited its success. We present a novel method for a cluster mass reconstruction which combines weak and strong lensing information on common scales and can, as a consequence, break the mass-sheet degeneracy. We extend the weak lensing formalism to the inner parts of the cluster and combine it with the constraints from multiple image systems. We demonstrate the feasibility of the method with simulations, finding an excellent agreement between the input and reconstructed mass also on scales within and beyond the Einstein radius. Using a single multiple image system and photometric redshift information of the background sources used for weak and strong lensing analysis, we find that we are effectively able to break the mass-sheet degeneracy, therefore removing one of the main limitations on cluster mass estimates. We conclude that with high resolution (e.g. HST) imaging data the method can more accurately reconstruct cluster masses and their profiles than currently existing lensing techniques.

We study the relation of the three-point cosmic shear statistics to the third-order statistical properties of the underlying convergence, expressed in terms of its bispectrum. Explicit relations for the natural components of the shear three-point correlation function in terms of the bispectrum are derived. The behavior of the correlation function under parity transformation is obtained and found to agree with previous results. We find that in contrast to the two-point shear correlation function, the three-point function at a given angular scale \theta is not affected by power in the bispectrum on much larger scales. These relations are then inverted to obtain the bispectrum in terms of the three-point shear correlator; two different expressions, corresponding to different natural components of the shear correlator, are obtained and can be used to separate E and B-mode shear contributions. These relations allow us to explicitly show that correlations containing an odd power of B-mode shear vanish for parity-symmetric fields. Generalizing a recent result by Jarvis et al., we derive expressions for the third-order aperture measures, employing multiple angular scales, in terms of the (natural components of the) three-point shear correlator and show that they contain essentially all the information about the underlying bispectrum. We discuss the many useful features these (generalized) aperture measures have that makes them convenient for future analyses of the skewness of the cosmic shear field (and any other polar field, such as the polarization of the Cosmic Microwave Background). (Abridged)

Galaxy-galaxy lensing (GGL) measures the 2-point cross-correlation between galaxies and mass in the Universe. In this work we seek to generalise this effect by considering the third-order correlations between galaxies and mass: galaxy-galaxy-galaxy lensing. Third-order correlations in the cosmic shear field have recently been reported in the VIRMOS-DESCART and CTIO surveys. Such data should also be ideal for measuring galaxy-galaxy-galaxy lensing. Indeed, the effects of these higher-order correlations may have already been detected in recent studies of galaxy-galaxy lensing. In this paper we lay out the basic relations for third-order cross correlations and their projections and introduce a new set of scale dependent third-order bias parameters. We define three new observables: two galaxy-shear-shear correlation functions, $G_{\pm}$, and a galaxy-galaxy-shear correlation, ${\cal G}$. We relate these to the various projected cross-bispectra and give practical estimators for their measurement. We note that the observational signature of these correlators is simply the excess shear-shear correlation measured about foreground galaxies (for $G_{\pm}$) and the average tangential shear around foreground galaxy pairs (for $\cal G$). Finally we derive expressions for the third order aperture mass statistics in terms of both the cross-bispectra and the real-space correlation functions. (Abridged)

Weak lensing surveys provide the possibility of identifying dark matter halos based on their total matter content rather than just the luminous matter content. On the basis of two sets of observations carried out with the CFHT, Erben et al. (2000) presented the first candidate dark clump, i.e. a dark matter concentration identified by its significant weak lensing signal without a corresponding galaxy overdensity or X-ray emission. We present a set of HST mosaic observations which confirms the presence of an alignment signal at the dark clump position. The signal strength, however, is weaker than in the ground-based data. It is therefore still unclear whether the signal is caused by a lensing mass or is just a chance alignment. We also present Chandra observations of the dark clump, which fail to reveal any significant extended emission. A comparison of the ellipticity measurements from the space-based HST data and the ground-based CFHT data shows a remarkable agreement on average, demonstrating that weak lensing studies from high-quality ground-based observations yield reliable results.

We present our image processing system for the reduction of optical imaging data from multi-chip cameras. In the framework of the Garching Bonn Deep Survey (GaBoDS; Schirmer et al. 2003) consisting of about 20 square degrees of high-quality data from WFI@MPG/ESO 2.2m, our group developed an imaging pipeline for the homogeneous and efficient processing of this large data set. Having weak gravitational lensing as the main science driver, our algorithms are optimised to produce deep co-added mosaics from individual exposures obtained from empty field observations. However, the modular design of our pipeline allows an easy adaption to different scientific applications. Our system has already been ported to a large variety of optical instruments and its products have been used in various scientific contexts. In this paper we give a thorough description of the algorithms used and a careful evaluation of the accuracies reached. This concerns the removal of the instrumental signature, the astrometric alignment, photometric calibration and the characterisation of final co-added mosaics. In addition we give a more general overview on the image reduction process and comment on observing strategies where they have significant influence on the data quality.

We investigate how galaxy-galaxy lensing measurements depend on the knowledge of redshifts for lens and source galaxies. Galaxy-galaxy lensing allows one to study dark matter halos of galaxies statistically using weak gravitational lensing. Redshift information is required to reliably distinguish foreground lens galaxies from background source galaxies and to convert the measured shear into constraints on the lens model. Without spectroscopy or multi-colour information, redshifts can be drawn from independently estimated probability distributions. The COMBO-17 survey provides redshifts for both lens and source galaxies. It thus offers the unique possibility to do this investigation with observational data. We find that it is of great importance to know the redshifts of individual lens galaxies in order to constrain the properties of their dark matter halos. Whether the redshifts are derived from $UBVRI$ or the larger number of filters available in COMBO-17 is not very important. In contrast, knowledge of individual source redshifts improves the measurements only very little over the use of statistical source redshift distributions as long as the source redshift distribution is known accurately.

We present a measurement of mass estimates for dark matter halos around galaxies from the COMBO-17 survey using weak gravitational lensing. COMBO-17 is particularly useful for this kind of investigation because it covers observations in 17 optical filters from which accurate photometric redshifts and spectral classification for objects with $R<24$ are derived. This allows us to select lens and source galaxies from their redshifts and to thus avoid any uncertainties from estimates of the source redshift distribution. We study galaxy lenses at redshifts $z_\mathrm{d}=0.2-0.7$ by fitting the normalization of either singular isothermal spheres (SIS) or Navarro-Frenk-White (NFW) profiles to the whole lens sample; we then consider halos around blue and red subsamples separately. We also constrain the scaling of halo mass with light. For the NFW model, we find virial masses $M_\mathrm{vir}^*=3.9^{+3.3}_{-2.4}\times 10^{11}h^{-1}M_{\sun}$ (1-$\sigma$) for blue and $M_\mathrm{vir}^*=7.1^{+7.1}_{-3.8}\times 10^{11}h^{-1}M_{\sun}$ for red galaxies of $L_\star=10^{10}h^{-2}L_{\sun}$, respectively. The derived mass-to-light scaling relations suggest that the mass-to-light ratio might decrease with increasing luminosity for blue galaxies but increase with increasing luminosity for red galaxies. However, these differences between blue and red galaxies are only marginally significant and both subsamples are consistent with having the same mass-to-light ratio at all luminosities. Finally, we compare our results to those obtained from the Red-Sequence Cluster Survey (RCS) and the Sloan Digital Sky Survey (SDSS). Taking differences in the actual modelling into account, we find very good agreement with these surveys.

We measured the bias and correlation factor of galaxies with respect to the dark matter using the aperture statistics including the aperture mass from weak gravitational lensing. The analysis was performed for three galaxy samples selected by R-band magnitudes; the median redshifts of the samples were $\ave{z}=0.34,0.49$ and 0.65, respectively. The brightest sample has the strongest peak in redshift and can therefore be most accurately represented by a single redshift. The data used is the GaBoDS, and the COMBO-17 survey for an accurate estimate of the redshift distribution of the galaxies. Assuming the currently favoured $\Lambda\rm CDM$ model as cosmology, we obtained values for the linear stochastic galaxy-dark matter bias on angular scales $1^\prime\le\theta\le20^\prime$. Averaging the measurements for the bias over the range $4^\prime\le \theta_{\rm ap}\le 18^\prime$, weighting with the cosmic variance error, yields $b=0.89(5), 0.79(5), 0.89(5)$, respectively ($1\sigma$). As for the correlation factor $r$ we see, as Hoekstra et al. (2002), a slight increase to $r=1$ in the last angular bin from an almost constant value on smaller scales; the weighted average here over the same range as before is $r=0.8(1), 0.8(1), 0.5(1)$, respectively (abridged).

We study the dark matter halos of galaxies with galaxy-galaxy lensing using the COMBO-17 survey. This survey offers an unprecedented data set for studying lens galaxies at z=0.2-0.7 including redshift information and spectral classification from 17 optical filters for objects brighter than R=24. So far, redshifts and classification for the lens galaxies have mainly been available for local surveys like the Sloan Digital Sky Survey (SDSS). Further, redshifts for the source galaxies have typically not been available at all but had to be estimated from redshift probability distribution which -- for faint surveys -- even had to be extrapolated. To study the dark matter halos we parametrize the lens galaxies as singular isothermal spheres (SIS) or by Navarro-Frenk-White (NFW) profiles. In both cases we find a dependence of the velocity dispersion or virial radius, respectively, on lens luminosity and colour. For the SIS model, we are able to reproduce the Tully-Fisher/Faber-Jackson relation on a scale of 150h^-1kpc. For the NFW profile we also calculate virial masses, mass-to-light ratios and rotation velocities. Finally, we investigate differences between the three survey fields used here.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, weak lensing mass reconstructions are plagued by the so-called mass-sheet degeneracy--the surface mass density \kappa of the cluster can be determined only up to a degeneracy transformation \kappa \to \kappa' = \lambda \kappa + (1 -\lambda), where \lambda is an arbitrary constant. This transformation fundamentally limits the accuracy of cluster mass determinations if no further assumptions are made. We discuss here a possibility to break the mass-sheet degeneracy in weak lensing mass maps using distortion and redshift information of background galaxies. Compared to other techniques proposed in the past, it does not rely on any assumptions on cluster potential and does not make use of weakly constrained information (such as the source number counts, used in the magnification effect).Our simulations show that we are effectively able to break the mass-sheet degeneracy for supercritical lenses and that for undercritical lenses the mass-sheet degeneracy is very difficult to be broken, even under idealised conditions.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, weak lensing mass reconstructions are plagued by the so-called mass-sheet degeneracy--the surface mass density \kappa of the cluster can be determined only up to a degeneracy transformation \kappa \to \kappa' = \lambda \kappa + (1 -\lambda), where \lambda is an arbitrary constant. This transformation fundamentally limits the accuracy of cluster mass determinations if no further assumptions are made. We describe here a method to break the mass-sheet degeneracy in weak lensing mass maps using distortion and redshift information of background galaxies and illustrate this by two simple toy models. Compared to other techniques proposed in the past, it does not rely on any assumptions on cluster potential; it can be easily applied to non-parametric mass-reconstructions and no assumptions on boundary conditions have to be made. In addition it does not make use of weakly constrained information (such as the source number counts, used in the magnification effect). Our simulations show that we are effectively able to break the mass-sheet degeneracy for supercritical lenses, but that for undercritical lenses the mass-sheet degeneracy is very difficult to be broken, even under idealised conditions.

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.

We present a weak lensing analysis and comparison to optical and X-ray maps of the close pair of massive clusters A222/223. Indications for a filamentary connection between the clusters are found and discussed.

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 present simulations of a cosmic shear survey and show how the survey geometry influences the accuracy of determination of cosmological parameters. We numerically calculate the full covariance matrices Cov of two-point statistics of cosmic shear, based on the expressions derived in the first paper of this series. The individual terms are compared for two survey geometries with large and small cosmic variance. We use analyses based on maximum likelihood of Cov and the Fisher information matrix in order to derive expected constraints on cosmological parameters. As an illustrative example, we simulate various survey geometries consisting of 300 individual fields of 13'x13' size, placed (semi-)randomly into patches which are assumed to be widely separated on the sky and therefore uncorrelated. Using the aperture mass statistics, the optimum survey consists of 10 patches with 30 images in each patch. If \Omega_m, \sigma_8 and \Gamma are supposed to be extracted from the data, the minimum variance bounds on these three parameters are 0.17, 0.25 and 0.04 respectively. These variances raise slightly when the initial power spectrum index n_s is also to be determined from the data. The cosmological constant is only poorly constrained.

We report the first confirmation of colour-selected galaxy cluster candidates by means of weak gravitational lensing. Significant lensing signals were identified in the course of the shear-selection programme of dark matter haloes in the Garching-Bonn Deep Survey, which currently covers 20 square degrees of deep, high-quality imaging data on the southern sky. The detection was made in a field that was previously covered by the ESO Imaging Survey (EIS) in 1997. A highly significant shear-selected mass-concentration perfectly coincides with the richest EIS cluster candidate at z~0.2, thus confirming its cluster nature. Several other shear patterns in the field can also be identified with cluster candidates, one of which could possibly be part of a filament at z~0.45.

Following the detection of a cosmic shear signal at the 30'' scale using archival parallel data from the STIS CCD camera onboard HST in Haemmerle et al. (2002), we analyzed a larger data set obtained from an HST GO pure parallel program. Although this data set is considerably larger than the one analyzed previously, we do not obtain a significant detection of the cosmic shear signal. The potential causes of this null result are the multiple systematics that plague the STIS CCD data, and in particular the degradation of the CCD charge transfer efficiency after 4 years in space.

We consider the parity transformation and the consequences of parity invariance on the $n$-point correlation function of the cosmic shear, or any other polar (or `spin-2') field. The decomposition of the shear field into E- and B-modes then yields the result that any correlation function which contains an odd number of B-mode shear components vanishes for parity-invariant random shear fields. This result is at variance with some statements in the literature, and the reasons for this is explained. We argue that the parity invariance provides a very specific diagnostic of systematic effects in shear data.

We perform a weak lensing analysis on a previously reported dark-lens candidate on STIS Parallel data (Miralles et al. 2002). New VLT-data indicate that the reported signal originates from a small number of galaxies tangentially aligned towards the center of the STIS field but no signature for an extended mass distribution is found. We argue that we should be able to detect a massive cluster ($M\geq 3.2\times 10^{14}M_{\odot}$) through its lensing signal up to a redshift of $z\approx 0.6$ with our data. Also the double image hypothesis of two galaxies with very similar morphologies in the STIS data is ruled out with colour information.

Gravitational lensing has become one of the most interesting tools to study the mass distribution in the Universe. Since gravitational light deflection is independent of the nature and state of the matter, it is ideally suited to investigate the distribution of all (and thus also of dark) matter in the Universe. Lensing results have now become available over a wide range of scales, from the search for MACHOs in the Galactic halo, to the mass distribution in galaxies and clusters of galaxies, and the statistical properties of the large-scale matter distribution in the Universe. Here, after introducing the concepts of strong and weak lensing, several applications are outlined, from strong lensing by galaxies, to strong and weak lensing by clusters and the lensing properties of the large-scale structure.

In a series of papers (Lombardi & Schneider 2001, 2002) we studied in detail the statistical properties of an interpolation technique widely used in astronomy. In particular, we considered the average interpolated map and its covariance under the hypotheses that the map is obtained by smoothing unbiased measurements of an unknown field, and that the measurements are uniformly distributed on the sky. In this paper we generalize the results obtained to the case of observations carried out only on a finite field and distributed on the field with a non-uniform density. These generalizations, which are required in many astronomically relevant cases, still allow an exact, analytical solution of the problem. We also consider a number of properties of the interpolated map, and provide asymptotic expressions for the average map and the two-point correlation function which are valid at high densities.

The Garching-Bonn Deep Survey (GaBoDS) is a virtual 12 square degree cosmic shear and cluster lensing survey, conducted with the WFI@2.2m MPG/ESO telescope at La Silla. It consists of shallow, medium and deep random fields taken in R-band in subarcsecond seeing conditions at high galactic latitude. A substantial amount of the data was taken from the ESO archive, by means of a dedicated ASTROVIRTEL program. In the present work we describe the main characteristics and scientific goals of GaBoDS. Our strategy for mining the ESO data archive is introduced, and we comment on the Wide Field Imager data reduction as well. In the second half of the paper we report on clusters of galaxies found in the background of NGC 300, a random archival field. We use weak gravitational lensing and the red cluster sequence method for the selection of these objects. Two of the clusters found were previously known and already confirmed by spectroscopy. Based on the available data we show that there is significant evidence for substructure in one of the clusters, and an increasing fraction of blue galaxies towards larger cluster radii. Two other mass peaks detected by our weak lensing technique coincide with red clumps of galaxies. We estimate their redshifts and masses.

We use galaxy-galaxy lensing to study the statistical properties of the dark matter halos of galaxies. Our data is taken from the COMBO-17 survey which has imaged 1 square degree in 17 optical filters giving accurate photometric redshifts and spectral classification down to R=24. This allows us to study lens galaxies at $z=0.2-0.7$. Shapes are measured from deep $R$-band images taken at the best seeing conditions ($0\farcs 8$ PSF). We model the lens galaxies as singular isothermal spheres (SIS). Investigating the dependence of the velocity dispersion on the luminosities of the lens galaxies we are able to reproduce the Tully-Fisher/Faber-Jackson relation. Further we find a larger lensing signal for early-type galaxies than for the late-types. However, when testing the dependence of the velocity dispersion on radius we see clear deviations from the SIS model. The velocity dispersion is first rising and then declining again. This is exactly what is expected for lens galaxies that can be modeled by Navarro-Frenk-White (NFW) profiles.

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.

Previous work on three-point statistics of cosmic shear has mainly concentrated on the convergence, or on aperture measures of the shear. However, as has become clear recently for the two-point statistics of cosmic shear, the basic quantity that should be used is the correlation function: first, it is much easier to measure from observational data, since it is immune against complicated geometries of data fields; second, all other (linear) 2-pt statistics can be expressed as integrals over the correlation function. The situation is the same for the 3-pt statistics. However, in contrast to the 2-pt correlation function, the invariants of the shear 3-pt correlation function have not been employed yet. Here we consider the transformation properties of the shear 3-pt correlation function under rotations. We show that there are four complex linear combinations of components of the 3-pt correlation function, which we shall call `natural components', since they are multiplied just by a phase factor for arbitrary rotations, but do not mix. In particular, their moduli are invariant under rotations and thus (non-linear) invariants of the 3-pt correlation function. In terms of these natural components, the invariance of the statistical properties of the shear field under parity transformations are easily obtained. Our results do not apply only to cosmic shear, but also to other quantities with the same mathematical properties -- that of a polar. For example, practically every relation derived here applies also to the polarization of the cosmic microwave background radiation.

The present status of weak lensing analyses of clusters of galaxies and of cosmic shear surveys are presented and discussed. We focus on the impact of very large telescopes on present-day and future surveys and compare their potential with HST or wide field 4 meter telescopes.

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.

We present mass reconstructions from weak lensing for the galaxy clusters A1835 and A2204 over 34'x34' fields using data from the ESO/MPG Wide Field Imager. Using a background galaxy population of 22<R<25.5 we detect the gravitational shear of A1835 at 8.8 sigma significance, and obtain best-fit mass profiles of sigma_v=1233^{+66}_{-70} km/s for a singular isothermal sphere model and r_{200}=1550 h^{-1} kpc, c=2.96 for a `universal' CDM profile. Using a color-selected background galaxy population of 22<R<25.8 we detect the gravitational shear of A2204 at 7.2 sigma significance, and obtain best-fit mass profiles of sigma_v=1035^{+65}_{-71} km/s for a SIS model and r_{200}=1310 h^{-1} km/s, c=6.3 for a `universal' CDM profile. The gravitational shear at distances greater than 10' is significantly detected for both clusters. The best fit weak lensing cluster masses agree well with both X-ray and dynamical mass measurements, although the central concentration of A1835 is much lower in the weak lensing mass profile than that measured by recent Chandra results. We suggest that this lower concentration is most likely a combination of contamination of the 'background' galaxy population with cluster dwarf galaxies and the effect of a prolate or tri-axial cluster core with the major axis lying near the plane of the sky. We also detect a number of additional structures at moderate significance, some of which appear to be sub-haloes associated with the clusters.

Weak gravitational lensing by the large scale structure can be used to probe the dark matter distribution in the Universe directly and thus to probe cosmological models. The recent detection of cosmic shear by several groups has demonstrated the feasibility of this new mode of observational cosmology. In the currently most extensive analysis of cosmic shear, it was found that the shear field contains unexpected modes, so-called B-modes, which are thought to be unaccountable for by lensing. B-modes can in principle be generated by an intrinsic alignment of galaxies from which the shear is measured, or may signify some remaining systematics in the data reduction and analysis. In this paper we show that B-modes in fact {\it are produced} by lensing itself. The effect comes about through the clustering of source galaxies, which in particular implies an angular separation-dependent clustering in redshift. After presenting the theory of the decomposition of a general shear field into E- and B-modes, we calculate their respective power spectra and correlation functions for a clustered source distribution. Numerical and analytical estimates of the relative strength of these two modes show that the resulting B-mode is very small on angular scales larger than a few arcminutes, but its relative contribution rises quickly towards smaller angular scales, with comparable power in both modes at a few arcseconds. The relevance of this effect with regard to the current cosmic shear surveys is discussed.

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.

Interpolation techniques play a central role in Astronomy, where one often needs to smooth irregularly sampled data into a smooth map. In a previous article (Lombardi & Schneider 2001), we have considered a widely used smoothing technique and we have evaluated the expectation value of the smoothed map under a number of natural hypotheses. Here we proceed further on this analysis and consider the variance of the smoothed map, represented by a two-point correlation function. We show that two main sources of noise contribute to the total error budget and we show several interesting properties of these two noise terms. The expressions obtained are also specialized to the limiting cases of low and high densities of measurements. A number of examples are used to show in practice some of the results obtained.

We derive in this paper expressions for the covariance matrix of the cosmic shear two-point correlation functions which are readily applied to any survey geometry. Furthermore, we consider the more special case of a simple survey geometry which allows us to obtain approximations for the covariance matrix in terms of integrals which are readily evaluated numerically. These results are then used to study the covariance of the aperture mass dispersion which has been employed earlier in quantitative cosmic shear analyses. We show that the aperture mass dispersion, measured at two different angular scales, quickly decorrelates with the ratio of the scales. Inverting the relation between the shear two-point correlation functions and the power spectrum of the underlying projected matter distribution, we construct estimators for the power spectrum and for the band powers, and show that they yields accurate approximations; in particular, the correlation between band powers at different wave numbers is quite weak. The covariance matrix of the shear correlation function is then used to investigate the expected accuracy of cosmological parameter estimates from cosmic shear surveys. Depending on the use of prior information, e.g. from CMB measurements, cosmic shear can yield very accurate determinations of several cosmological parameters, in particular the normalization $\sigma_8$ of the power spectrum of the matter distribution, the matter density parameter $\Omega_{\rm m}$, and the shape parameter $\Gamma$.

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 report on the marginal detection of cosmic shear on sub-arcmin scales with archive data from the STIS camera on board HST. For the measurement 121 galaxy fields with a field of view of 51'' x 51'' are used to obtain an rms cosmic shear of ~ 4% with 1.5sigma significance. This value is consistent with groundbased results obtained by other groups on larger scales, and with theoretical predictions for a standard LambdaCDM cosmology. To show the suitability of STIS for weak shear measurements we carefully investigated the stability of the PSF. We demonstrate that small temporal changes do not affect the cosmic shear measurement by more than ~ 10%. We also discuss the influence of various weighting and selection schemes for the galaxy ellipticities.

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.

We report the serendipituous discovery of a conspicuous alignment of galaxies in a field obtained through the STIS Parallel Shear Survey. This project collects randomly distributed 50'' x 50'' fields to investigate the cosmic shear effect on this scale. Analyzing the parallel observations having the Seyfert galaxy NGC625 as primary target, we recognized over the whole field of view a strong apparent tangential alignment of galaxy ellipticities towards the image center. The field shows several arclet-like features typical for images of massive galaxy clusters, but no obvious over-density of bright foreground galaxies. We also find a multiple image candidate. On the basis of the possible strong and weak lensing effect within the data, we discuss whether this could be compatible with a massive halo with no clear optical counterpart.

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.

The parameter-free reconstruction of the surface-mass density of clusters of galaxies is one of the principal applications of weak gravitational lensing. From the observable ellipticities of images of background galaxies, the tidal gravitational field (shear) of the mass distribution is estimated, and the corresponding surface mass density is constructed. The noise of the resulting mass map is investigated here, generalizing previous work which included mainly the noise due to the intrinsic galaxy ellipticities. Whereas this dominates the noise budget if the lens is very weak, other sources of noise become important, or even dominant, for the medium-strong lensing regime close to the center of clusters. In particular, shot noise due to a Poisson distribution of galaxy images, and increased shot noise owing to the correlation of galaxies in angular position and redshift, can yield significantly larger levels of noise than that from the intrinsic ellipticities only. We estimate the contributions from these various effects for two widely used smoothing operations, showing that one of them effectively removes the Poisson and the correlation noises related to angular positions of galaxies. Noise sources due to the spread in redshift of galaxies are still present in the optimized estimator and are shown to be relevant in many cases. We show how (even approximate) redshift information can be profitably used to reduce the noise in the mass map. The dependence of the various noise terms on the relevant parameters (lens redshift, strength, smoothing length, redshift distribution of background galaxies) are explicitly calculated and simple estimates are provided.

Since the Universe is inhomogeneous on scales well below the Hubble radius, light bundles from distant galaxies are deflected and distorted by the tidal gravitational field of the large-scale matter distribution as they propagate through the Universe. Two-point statistical measures of the observed ellipticities, like the dispersion within a finite aperture or the ellipticity cross-correlation, can be related to the power spectrum of the large-scale structure. The measurement of cosmic shear (especially on small angular scales) can thus be used to constrain cosmological parameters and to test cosmological structure formation in the non-linear regime, without any assumptions about the relation between luminous and dark matter. In this paper we will present preliminary cosmic shear measurements on sub-arcminute scales, obtained from archival STIS parallel data. The high angular resolution of HST, together with the sensitivity and PSF-stability of STIS, allows us to measure cosmic shear along many independent lines-of-sight. Ongoing STIS parallel observations, currently being carried out in the frame of a big GO program (8562+9248), will greatly increase the number of available useful fields and will enable us to measure cosmic shear with higher accuracy on sub-arcminute scales.

The measurement of cosmic shear requires deep imaging with high image quality on many lines of sight to sample the statistics of large-scale structure. The expected distortion of galaxy images by cosmic shear on the STIS angular scale is a few percent, therefore the PSF anisotropy has to be understood and controlled to an accuracy better than 1%. In this poster we present the analysis of the PSF of STIS and a preliminary cosmic shear measurement using archival data from the STIS pure parallel program to show that the STIS camera on-board HST is well suited for our project. The data reduction and catalog production are described in an accompanying paper (astro-ph/0102330).

In June 1997, parallel observations using the Space Telescope Imaging Spectrograph (STIS) on the HST started to be taken in substantial numbers along many different lines-of-sight. We are using the imaging data to investigate the distortion of background galaxies by the gravitational field of the large scale matter distribution, also known as Cosmic Shear. This poster presents the data and the catalog production that leads to the cosmic shear result presented in the poster "First Cosmic Shear results from STIS parallel program archive data" (Haemmerle et al., this conference). The data is publicly available also at http://www.stecf.org/projects/shear .

We present a mass profile for A1689 from 0.13 h^{-1} Mpc < r < 2 h^{-1} Mpc from a weak lensing analysis of a 34'x34' R-band image from the ESO/MPG Wide Field Imager. We detect the gravitational shearing of a 23<R<25.5 background galaxy population even at the edge of the image with a 4-sigma significance, and find a two-dimensional mass reconstruction has a 13.4-sigma significance mass peak centered on the brightest cluster galaxy. This peak is well fit by both a 1028+-35 km/s singular isothermal sphere and a r200 = 1.28 Mpc, c = 6 universal CDM profile, although the universal CDM profile provides a better fit with 95.5% confidence. These mass measurements are lower than most of those derived by other means and we discuss possible reasons for weak lensing providing an underestimate of the true mass of the cluster. We find that the correction factors needed to reconcile the weak lensing mass models with the strong lensing Einstein radius would result is a much larger fraction of faint stars and foreground and cluster dwarf galaxies in the 23<R<25.5 object catalog than is seen in other fields.

We present the weak lensing mass reconstructions of three clusters of galaxies from data taken on the ESO Wide Field Imager. We detect a lensing signal in all three clusters over a range of 1' < r < 17' (~0.1 h^-1 Mpc < r < 2 h^-1 Mpc). We measure the best fitting SIS and NFW profiles for each cluster and show that in two of the three clusters the NFW profile is preferred to the SIS with marginal significance (~95%).

Following the second HST servicing mission in 1997 when the STIS instrument was installed and the capability for parallel observations was enhanced, a substantial archive of non-proprietary parallel data has been accumulating. In this paper, we discuss the use of unfiltered STIS imaging data for a project that requires deep observations along as many independent lines-of-sight as possible. We have developed a technique to determine which datasets in the archive can safely be co-added together and have developed an iterative co-addition technique which enabled us to produce 498 high-quality, deep images. The principal motivation for this work is to measure the Cosmic Shear on small angular scales and a value derived from these data will be presented in a subsequent paper. A valuable by-product of this work is a set of high quality combined fields which can be used for other projects. The data are publicly available at http://www.stecf.org/projects/shear/

We study an estimator for smoothing irregularly sampled data into a smooth map. The estimator has been widely used in astronomy, owing to its low level of noise; it involves a weight function -- or smoothing kernel -- w(\theta). We show that this estimator is not unbiased, in the sense that the expectation value of the smoothed map is not the underlying process convolved with $w$, but a convolution with a modified kernel w_eff(\theta). We show how to calculate w_eff for a given kernel w and investigate its properties. In particular, it is found that (1) w_eff is normalized, (2) has a shape `similar' to the original kernel w, (3) converges to w in the limit of high number density of data points, and (4) reduces to a top-hat filter in the limit of very small number density of data points. Hence, although the estimator is biased, the bias is well understood analytically, and since w_eff has all the desired properties of a smoothing kernel, the estimator is in fact very useful. We present explicit examples for several filter functions which are commonly used, and provide a series expression valid in the limit of large density of data points.

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.

With the recent detection of cosmic shear, the most challenging effect of weak gravitational lensing has been observed. The main difficulties for this detection were the need for a large amount of high quality data and the control of systematics during the gravitational shear measurement process, in particular those coming from the Point Spread Function anisotropy. In this paper we perform detailed simulations with the state-of-the-art algorithm developed by Kaiser, Squires and Broadhurst (KSB) to measure gravitational shear. We show that for realistic PSF profiles the KSB algorithm can recover any shear amplitude in the range $0.012 < |\gammavec |<0.32$ with a relative, systematic error of $10-15%$. We give quantitative limits on the PSF correction method as a function of shear strength, object size, signal-to-noise and PSF anisotropy amplitude, and we provide an automatic procedure to get a reliable object catalog for shear measurements out of the raw images.

We report a measurement of cosmic shear correlations using an effective area of 6.5 sq. deg. of the VIRMOS deep imaging survey in progress at the Canada-France-Hawaii Telescope. We measured various shear correlation functions, the aperture mass statistic and the top-hat smoothed variance of the shear with a detection significance exceeding 12 sigma for each of them. We present results on angular scales from 3 arc-seconds to half a degree. The consistency of different statistical measures is demonstrated and confirms the lensing origin of the signal through tests that rely on the scalar nature of the gravitational potential. For Cold Dark Matter models we find $\sigma_8 \Omega_0^{0.6}=0.43^{+0.04}_{-0.05}$ at the 95% confidence level. The measurement over almost three decades of scale allows to discuss the effect of the shape of the power spectrum on the cosmological parameter estimation. The degeneracy on sigma_8-Omega_0 can be broken if priors on the shape of the linear power spectrum (that can be parameterized by Gamma) are assumed. For instance, with Gamma=0.21 and at the 95% confidence level, we obtain 0.6<sigma_8<1.1 and 0.2<Omega_0<0.5 for open models, and sigma_8>0.65 and Omega_0<0.4 for flat (Lambda-CDM) models. From the tangential/radial modes decomposition we can set an upper limit on the intrinsic shape alignment, which was recently suggested as a possible contribution to the lensing signal. Within the error bars, there is no detection of intrinsic shape alignment for scales larger than 1'.

Gravitational weak shear produced by large-scale structures of the universe induces a correlated ellipticity distribution of distant galaxies. The amplitude and evolution with angular scale of the signal depend on cosmological models and can be inverted in order to constrain the power spectrum and the cosmological parameters. We present our recent analysis of 50 uncorrelated VLT fields and the very first constrains on ($\Omega_m,\sigma_8$) and the nature of primordial fluctuations based on the join analysis of present-day cosmic shear surveys.

In this contribution the ongoing effort to build a statistical sample of clusters of galaxies over a wide range of redshifts to study the evolution of clusters and member galaxies is reviewed. The starting point for this project has been the list of candidate clusters identified from the I-band EIS-WIDE survey data. Since the completion of this survey, new optical/infrared observations have become available and have been used to confirm some of these candidates using the photometric data alone or in combination with the results of follow-up spectroscopic observations. Our preliminary results show that the yield of real physical associations from the original catalog is conservatively >60% and that a large sample of clusters in the southern hemisphere, extending to high-redshifts, is within reach.

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.

The recent detections of cosmic shear signal announced by several groups have demonstrated the feasibility of this challenging program and convinced astro- nomers of its potential for cosmology. Cosmic shear analysis demands to handle Gigabytes of data in order to probe several square degrees in subarcsecond deep imaging mode. The success of these surveys is sensitive to the designs of the observation strategy, the organization of the data reduction pipelines and the links of the data base with surveys like X-ray or spectroscopic follow up. We describe the cosmic shear surveys we have carried out at the VLT and at CFHT and the way we handle this huge data set in a more general context including the VIRMOS and the XMM-LSS surveys, and the future CMB surveys.

We observed with the camera FORS1 on the VLT (UT1, ANTU) 50 randomly selected fields and analyzed the cosmic shear inside circular apertures with diameter ranging from 0.5 to 5.0 arc-minutes. The images were obtained in optimal conditions by using the Service Observing proposed by ESO on the VLT, which enabled us to get a well-defined and homogeneous set of data. The 50 fields cover a 0.65 square-degrees area spread over more than 1000 square-degrees which provides a sample ideal for minimizing the cosmic variance. Using the same techniques as in Van Waerbeke et al. 2000, we measured the cosmic shear signal and investigated the systematics of the VLT sample. We find a significant excess of correlations between galaxy ellipticities on those angular scales. The amplitude and the shape of the correlation as function of angular scale are remarkably similar to those reported so far. Using our combined VLT and CFHT data and adding the results published by other teams we put the first joint constraints on Omega_0 and sigma_8 using cosmic shear surveys. From a deduced average of the redshift of the sources the combined data are consistent with sigma_8\simeq (0.59 \pm 0.03) Omega_0^{-0.47} (for a CDM power spectrum and $\Gamma=0.21$), in excellent agreement with the cosmological constraints obtained from the local cluster abundance.

The first detections of cosmic shear signal reported recently by 4 independent groups cover angular scales between one and 10 arcmin. On those scales, the cosmic shear is a signature of non-linear perturbations, like groups and clusters of galaxies. We present the results obtained by our team on CFHT and on the VLT and discuss its impact for the analysis of cluster abundances and cosmology

We present results of our cosmic shear survey in progress at the Canada- France-Hawaii Telescope (CFHT). We analysed 1.7 sq. degrees of data, out of which we measured a significant correlation of galaxy shape-orientation over several arcmin scale. We present measurements of the variance of the shear $<\gamma^2>$ and of the correlation functions $<e_r(\theta)e_r(0)>$, $<e_t(\theta)e_t(0)>$. The signal is consistent with gravitational lensing by large scale structures predictions. The level of residual systematics after PSF correction is small compared to our measured signal. We outline several evolutions of our work by using additional data covering 9 sq. degrees.

We report on the detection of cosmic shear on angular scales of 1.3-6.5 arcmin using 45 independent empty fields observed with the Very Large Telescope (VLT). This result confirms previous measurements obtained with the CFH Telescope at the same angular scales. We present the data analysis and the first preliminary results.

We report the detection of a significant (5.5 sigma) excess of correlations between galaxy ellipticities at scales ranging from 0.5 to 3.5 arc-minutes. This detection of a gravitational lensing signal by large-scale structure was made using a composite high quality imaging survey of 6300 arcmin^2 obtained at the Canada France Hawaii Telescope (CFHT) with the UH8K and CFH12K panoramic CCD cameras. The amplitude of the excess correlation is 2.2\pm 0.2 % at 1 arcmin scale, in agreement with theoretical predictions of the lensing effect induced by large-scale structure.We provide a quantitative analysis of systematics which could contribute to the signal and show that the net effect is small and can be corrected for. We show that the measured ellipticity correlations behave as expected for a gravitational shear signal. The relatively small size of our survey precludes tight constraints on cosmological models. However the data are in favor of cluster normalized cosmological models, and marginally reject Cold Dark Matter models with (Omega=0.3, sigma_8<0.6) or (Omega=1, sigma_8=1). The detection of cosmic shear demonstrates the technical feasibility of using weak lensing surveys to measure dark matter clustering and the potential for cosmological parameter measurements, in particular with upcoming wide field CCD cameras.

The aperture mass has been shown in a series of recent publications to be a useful quantitative tool for weak lensing studies, ranging from cosmic shear to the detection of a mass-selected sample of dark matter haloes. Quantitative analytical predictions for the aperture mass have been based on a number of simplifying assumptions. In this paper, we test the reliability of these assumptions and the quality of the analytic approximations, using ray-tracing simulations through a cosmological density field generated by very large N-body simulations. We find that those analytic predictions which take into account the non-linear evolution of the matter distribution, such as the dispersion of the aperture mass and the halo abundance, are surprisingly accurately reproduced with our numerical results, whereas the predictions for the skewness, based on quasi-linear theory, are rather imprecise. In particular, we verify numerically that the probability distribution of the aperture mass decreases exponentially for values much larger than the rms. Given the good overall agreement, comparisons between the observed distribution of the aperture mass and the theoretical values provide a powerful tool for testing cosmological models.

We review theory and applications of weak gravitational lensing. After summarising Friedmann-Lemaitre cosmological models, we present the formalism of gravitational lensing and light propagation in arbitrary space-times. We discuss how weak-lensing effects can be measured. The formalism is then applied to reconstructions of galaxy-cluster mass distributions, gravitational lensing by large-scale matter distributions, QSO-galaxy correlations induced by weak lensing, lensing of galaxies by galaxies, and weak lensing of the cosmic microwave background. Contents: Introduction - Cosmological Background - Gravitational Light Deflection - Principles of Weak Gravitational Lensing - Weak Lensing by Galaxy Clusters - Weak Cosmological Lensing - QSO Magnification Bias and Large-Scale Structure - Galaxy-Galaxy Lensing - The Impact of Weak Gravitational Lensing on the Microwave Background Radiation - Summary and Outlook

I shall outline the basic principles and some observational aspects of weak gravitational lensing, and discuss several applications of this powerful tool in observational cosmology. It will be explained why the applications have been restricted to optical observations up to now, and why SKA is going to change this. I conclude with a few general remarks on a comparison between SKA and the NGST, both being facilities which will provide a tremendous step forward in radio and near-IR astronomy, respectively, into completely unknown territory.

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.

A weak-lensing analysis of wide-field $V$- and $I$-band images centered on the cluster Abell 1942 has uncovered a mass concentration $\sim 7$ arcminutes South of the cluster center. A statistical analysis shows that the detections are highly significant. No strong concentration of bright galaxies is seen at the position of the mass concentration, though a slight galaxy number overdensity and a weak extended X-ray source are present about 1' away from its center. From the spatial dependence of the tangential alignment around the center of the mass concentration, we inferred a lower bound on the mass inside a sphere of radius $0.5 h^{-1}$\ts Mpc of $1\times 10^{14}h^{-1}M_\odot$, much higher than crude mass estimates based on X-ray data. No firm conclusion can be inferred about the nature of the clump. If it were a high-redshift cluster, the weak X-ray flux would indicate that it had an untypically low X-ray luminosity for its mass; if the X-ray emission were physically unrelated to the mass concentration, this conclusion would be even stronger. The search for massive halos by weak lensing enables us for the first time to select halos based on their mass properties only and to detect new types of objects, e.g., dark halos. The mass concentration in the field of A1942 may be the first example of such a halo.

Due to gravitational instability, an initially Gaussian density field develops non-Gaussian features as the Universe evolves. The most prominent non-Gaussian features are massive haloes, visible as clusters of galaxies. The distortion of high-redshift galaxy images due to the tidal gravitational field of the large-scale matter distribution, called cosmic shear, can be used to investigate the statistical properties of the LSS. In particular, non-Gaussian properties of the LSS will lead to a non-Gaussian distribution of cosmic-shear statistics. The aperture mass ($M_{\rm ap}$) statistics, recently introduced as a measure for cosmic shear, is particularly well suited for measuring these non-Gaussian properties. In this paper we calculate the highly non-Gaussian tail of the aperture mass probability distribution, assuming Press-Schechter theory for the halo abundance and the `universal' density profile of haloes as obtained from numerical simulations. We find that for values of $M_{\rm ap}$ much larger than its dispersion, this probability distribution is closely approximated by an exponential, rather than a Gaussian. We determine the amplitude and shape of this exponential for various cosmological models and aperture sizes, and show that wide-field imaging surveys can be used to distinguish between some of the currently most popular cosmogonies. Our study here is complementary to earlier cosmic-shear investigations which focussed more on two-point statistical properties.

We present radio observations at three frequencies and contemporaneous optical monitoring of the peculiar BL Lac object AO 0235+164. During a three-week campaign with the VLA we observed intraday variability in this source and found a distinct peak which can be identified throughout the radio frequencies and tentatively connected to the R-band variations. This event is characterized by unusual properties: its strength increases, and its duration decreases with wavelength, and it peaks earlier at 20cm than at 3.6 and 6cm. We discuss several generic models (a `standard' shock-in-jet model, a precessing beam, free-free-absorption in a foreground screen, interstellar scattering, and gravitational microlensing), and explore whether they can account for our observations. Most attempts at explaining the data on 0235+164 require an extremely small source size, which can be reconciled with the $10^{12}$ K inverse Compton limit only when the Doppler factor of the bulk flow is of order 100. However, none of the models is completely satisfactory, and we suggest that the observed variability is due to a superposition of intrinsic and propagation effects.

We demonstrate that the aperture mass as a measure for cosmic shear closely approximates (to better than ~5%) the scaled and shifted power spectrum of the projected mass density. This cosmological weak-lensing information can thus be used to directly infer the projected matter power spectrum with high accuracy. As an application, we show that aperture-mass observations can be used to constrain the cosmic density parameter and the power-spectrum amplitude. We show that, for a particular example, it should be possible to constrain Omega_0 to within ~ +/- 27%, and sigma_8 to within ~ +/- 8% using weak-shear data on a square-shaped field of 8 degrees side length.

The developments summarized with the name ``weak gravitational lensing'' have led to exciting possibilities to study the (statistical properties of the) dark matter distribution in the Universe. Concentrating on those aspects which require deep wide-field imaging surveys, we describe the basic principles of weak lensing and discuss its future applications in view of NGST (a) to determine the statistical properties of the dark matter halos of individual galaxies, (b) to determine the mass and the mass profile of very low-mass clusters and groups at medium redshift and/or of more massive clusters at very high redshift, and (c) to measure the power spectrum of the matter distribution in the Universe in the non-linear regime, thereby also obtaining a mass-selected sample of halos and providing a means to investigate the scale- and redshift dependence of the bias `factor'.

The distortion of the images of faint high-redshift galaxies can be used to probe the intervening mass distribution. This weak gravitational lensing effect has been used recently to study the (projected) mass distribution of several clusters at intermediate and high redshifts. In addition, the weak lensing effect can be employed to detect (dark) matter concentrations in the Universe, based on their mass properties alone. Thus it is feasible to obtain a mass-selected sample of `clusters', and thereby probe the full range of their mass-to-light ratios. We study the expected number density of such haloes which can be detected in ongoing and future deep wide-field imaging surveys, using the number density of haloes as predicted by the Press-Schechter theory, and modeling their mass profile by the `universal' density profile found by Navarro, Frenk & White. We find that in all cosmological models considered, the number density of haloes with a signal-to-noise ratio larger than 5 exceeds 10 per square degree. With the planned MEGACAM imaging survey of $\sim 25 deg^2$, it will be easily possible to distinguish between the most commonly discussed cosmological parameter sets.

In a previous paper we investigated means for constraining the mass distribution of cluster galaxies by weak lensing. We concluded that a comprehensive method should treat the lensing effects of individual cluster galaxies and those resulting from a general cluster component simultaneously. To this end we now develop a non-parametric maximum likelihood cluster reconstruction algorithm that can implicitly take into account the presence of cluster galaxies. The method includes an entropy-like regularization prescription and directly uses the ellipticities of individual source galaxy images as observables rather than relying on an averaged ellipticity field. The mass distribution of cluster galaxies is described by parametrized models. For each set of galaxy parameters the cluster reconstruction algorithm allows to determine the best representation of the global underlying cluster component that is consistent with the presence of the cluster galaxies and the observed image ellipticities of background galaxies. Tests with simulations yielded convincing and robust results. We applied the method to a WFPC2 image of the cluster Cl0939+4713 and obtained a detection of the lensing effects of luminous elliptical cluster galaxies. We consider this application as a successful test of our technique. However, the small size of the image we analysed does not yet allow to draw strong conclusions regarding the mass distribution of cluster galaxies.

We present a new method for reconstructing two-dimensional mass maps of galaxy clusters from the image distortion of background galaxies. In contrast to most previous approaches, which directly convert locally averaged image ellipticities to mass maps (direct methods), our entropy-regularized maximum-likelihood method is an inverse approach. Albeit somewhat more expensive computationally, our method allows high spatial resolution in those parts of the cluster where the lensing signal is strong enough. Furthermore, it allows to straightforwardly incorporate additional constraints, such as magnification information or strong-lensing features. Using synthetic data, we compare our new approach to direct methods and find indeed a substantial improvement especially in the reconstruction of mass peaks. The main differences to previously published inverse methods are discussed.

A new method for the reconstruction of the projected mass distribution of clusters of galaxies from the image distortion of background galaxies is discussed. This method is essentially equivalent to the one we developed previously, i.e., the noise-filtering method, but has several practical advantages: (1) It is much easier to implement; (2) it can be easily applied to wide-field images, since the constraints on the number of gridpoints are much weaker than for the previous method, and (3) it can be easily generalized to more complicated field geometries, such as that of the Wide Field Planetary Camera 2 (WFPC2) onboard HST. We have tested the performance of our new inversion method (for which a FORTRAN-77 implementation is available from the authors) using simulated data, demonstrating that it fares very favourably.

The correlation between cosmic shear as measured by the image distortion of high-redshift galaxies and the number counts of foreground galaxies is calculated. For a given power spectrum of the cosmic density fluctuations, this correlation is proportional to the bias factor, which can thus directly be measured. In addition, this correlation provides a first-order measure of cosmic shear and is therefore easier to observe than quadratic measures hitherto proposed. Analytic approximations show that the expected signal-to-noise ratio of the correlation is large, so that a significant detection is possible with a moderate amount of data; in particular, it is predicted that the ongoing ESO Imaging Survey (EIS) will be able to detect this correlation on scales of $\sim 10'$ at a 3-$\sigma$ level, and at with higher significance on smaller angular scales.

We consider here a new statistical measure for cosmic shear, the aperture mass Map, which is defined as a spatially filtered projected density field and which can be measured directly from the image distortions of high-redshift galaxies. By selecting an appropriate spatial filter function, the dispersion of the aperture mass is a convolution of the power spectrum of the projected density field with a narrow kernel, so that <Map**2(theta)> provides a well localized estimate of the power spectrum at wavenumbers s~5/theta. We calculate <Map**2(theta)> for various cosmological models, using the fully non-linear power spectrum of the cosmic density fluctuations. The non-linear evolution yields a significant increase of <Map**2> relative to the linear growth on scales below ~0.5 degree. The third-order moment of Map can be used to define a skewness, which is a measure of the non-Gaussianity of the density field. We present the first calculation of the skewness of the shear in the frame of quasi-linear theory of structure growth. We show that it yields a sensitive measure of the cosmological model; in particular, it is independent of the normalization of the power spectrum. Several practical estimates for <Map**2> are constructed and their dispersions calculated. We show that one sq.deg. of a high-quality image is sufficient to detect the cosmic shear with the Map-statistics on scales below ~10 arcmin, and to estimate its amplitude with an accuracy of ~30% on scales below ~5 arcmin. (abbreviated)

We discuss whether one should expect that multiply imaged QSOs can be understood with `simple' lens models which contain a handful of parameters. Whereas for many lens systems such simple mass models yield a remarkably good description of the observed properties, there are some systems which are notoriously difficult to understand quantitatively. We argue that at least in one case (B 1422+231) these difficulties are not due to a `wrong' parametrization of the lens model, but that the discrepancy between observed and model-predicted flux ratios are due to substructure in the lens. Similar to microlensing for optical fluxes, such substructure can distort also the radio flux ratios predicted by `simple' mass models, in particular for highly magnified images, without appreciably changing image positions. Substructure also does not change the time delay significantly and therefore has little effect on the determination of the Hubble constant using time delays. We quantify these statements with several simple scenarios for substructure, and propose a strategy to model lens systems in which substructure is suspected.

Analysing the weak lensing distortions of the images of faint background galaxies provides a means to constrain the average mass distribution of cluster galaxies and potentially to test the extent of their dark matter haloes as a function of the density of their environment. The observable image distortions are a consequence of the interplay between the effects of a global cluster mass distribution and the perturbations due to individual cluster galaxies. Starting from a reconstruction of the cluster mass distribution with conventional techniques, we apply a maximum likelihood method to infer the average properties of an ensemble of cluster galaxies. From simulations this approach is found to be reliable as long as the galaxies including their dark matter haloes only contribute a small fraction to the total mass of the system. If their haloes are extended, the galaxies contain a substantial mass fraction. In this case our method is still applicable in the outer regions of clusters, where the surface mass density is low, but yields biased estimates of the parameters describing the mass profiles of the cluster galaxies in the central part of the cluster. In that case it will be necessary to resort to more sophisticated strategies by modelling cluster galaxies and an underlying global mass distribution simultaneously. We conclude that galaxy-galaxy lensing in clusters provides a unique means to probe the presence and extent of dark haloes of cluster galaxies.

The gravitational lens equation resulting from a single (non-linear) mass concentration (the main lens) plus inhomogeneities of the large-scale structure is shown to be strictly equivalent to the single-plane gravitational lens equation without the cosmological perturbations. The deflection potential (and, by applying the Poisson equation, also the mass distribution) of the equivalent single-plane lens is derived. If the main lens is described by elliptical isopotential curves plus a shear term, the equivalent single-plane lens will be of the same form. Due to the equivalence shown, the determination of the Hubble constant from time delay measurements is affected by the same mass-sheet invariance transformation as for the single-plane lens. If the lens strength is fixed (e.g., by measuring the velocity dispersion of stars in the main lens), the determination of $H_0$ is affected by inhomogeneous matter between us and the lens. The orientation of the mass distribution relative to the image positions is the same for the cosmological lens situation and the single-plane case. In particular this implies that cosmic shear cannot account for a misalignment of the observed galaxy orientation relative to the best-fitting lens model.

We present evidence for a coherent shear signal in a field containing the z=1.2 radio-source PKS1508-05. Since there were no intervening mass concentrations known before targeting this field, we interpret this signal as due to weak lensing by large-scale structure. This result is the outcome of a re-analysis of the observations of Fort et al. (1996) in the fields of three high redshift QSO/radio-sources. Several tests of the robustness of the signal were performed: subdivision of the field and tests using simulated data with randomized ellipticity orientation; using the pixel-to-pixel autocorrelation function method to measure the shear; using the correlation function of the galaxy ellipticities with bootstrap resampling to estimate errors and to find the typical angular separation that dominates the signal. We find that the field of PKS1508-05 contains a robust shear signal that is coherent over the 2' by 2' field. The signal is detected at a high level of significance (99.6%) using the tests described above. The amplitude of the shear on a scale of 1' is about 3%, which is consistent with theoretical expectations from weak lensing by large-scale structure.

The projected mass of a gravitational lens inside (circular) apertures can be derived from the measured shear inside an annulus which is caused by the tidal field of the deflecting mass distribution. Here we show that also the multipoles of the two-dimensional mass distribution can be derived from the shear in annuli. We derive several expressions for these mass multipole moments in terms of the shear, which allow large flexibility in the choice of a radial weight function. In contrast to determining multipole moments from weak-lensing mass reconstructions, this approach allows to quantify the signal-to-noise ratio of the multipole moments directly from the observed galaxy ellipticities, and thus to estimate the significance of the multipole detection. Radial weight functions can therefore be chosen such as to optimize the significance of the detection given an assumed radial mass profile. Application of our formulae to numerically simulated clusters demonstrates that the quadrupole moment of realistic cluster models can be detected with high signal-to-noise ratio S/N; in about 85 per cent of the simulated cluster fields S/N >~ 3. We also show that the shear inside a circular annulus determines multipole moments inside and outside the annulus. This is relevant for clusters whose central region is too bright to allow the observation of the shear of background galaxies, or which extend beyond the CCD. We also generalize the aperture mass equation to the case of `radial' weight functions which are constant on arbitrarily-shaped curves which are not necessarily self-similar.

A new method for measuring the shear induced by gravitational light deflection is proposed. It is based on analyzing the anisotropy induced in the auto-correlation function (ACF) of the extragalactic background light which is produced by very faint distant galaxies. The ACF can be measured `locally', and its anisotropy is caused by the tidal gravitational field of the deflecting mass distribution in the foreground of these faint background galaxies. Since the method does not require individual galaxy detection, it can be used to measure the shear of extremely faint galaxies which are not detectable individually, but are present in the noise. The shear estimated from the ACF of the noise provides an independent measurement which can be compared to the shear obtained from the distortion of individual galaxy images. Combining these two independent estimates clearly increases the sensitivity of shear measurements. In addition, our method may allow to determine the local magnification caused by the deflector if the auto-correlation function is caused by a large number density of faint galaxies; in this case, the intrinsic ACF may provide a `standard source' with respect to which shear and magnification can be obtained. Applications to real and synthetic data are shown and the feasibility of our new method is demonstrated. In particular, we present the shear maps obtained with our method for the double QSO 2345+007 and the cluster Cl0024+16 and compare them to published shear maps.

In this paper we explore a quantitative and efficient method to constrain the halo properties of distant galaxy populations through ``galaxy--galaxy" lensing and show that the mean masses and sizes of halos can be estimated accurately, without excessive data requirements. Specifically, we propose a maximum-likelihood analysis which takes full account of the actual image ellipticities, positions and apparent magnitudes. We apply it to simulated observations, using the same model for the lensing galaxy population as in BBS, where the galaxy halos are described by isothermal spheres with velocity dispersion $\sigma$, truncated at a radius $s$. Both parameters are assumed to scale with the luminosity of the galaxy. The best fitting values are then determined with the maximum-likelihood analysis. We explore two different observing strategies, (a) taking deep images (e.g., with HST) on small fields, and (b) using shallower images on larger fields. We find that $\sigma_*$ can be determined to $\lesssim$10\% accuracy if a sample of about 5000 galaxies with measured ellipticities are available, down to $R\lesssim 23$. Alternatively, the same accuracy in the determination of $\sigma_*$ can be achieved from about ten, moderately deep WFPC2 fields, on which galaxy shapes can be measured to about $R\sim 25$ and for which ground-based images are available on which the WFPC2 fields are centered. Firm lower limits can be set on the radial extent of the halo, but the maximal halo extent is poorly constrained. We show that the likelihood approach can also be used to constrain other parameters of the galaxy population, such as the Tully-Fischer index, or the mean redshift of the galaxies as a function of apparent magnitude. Finally we show how multi-color information, constraining the redshift of individual galaxies, can dramatically

We propose a new technique to measure the time delay of radio-loud gravitational lens systems, which does not rely on the excessive use of interferometric observations. Instead, the method is based on single-dish flux density monitoring of the (unresolved) lens system's total lightcurve, combined with additional interferometric measurements of the flux density ratio at a few epochs during that monitoring period. The basic idea of the method is to reconstruct the individual image lightcurves from the observed total lightcurve by assuming a range of potential values for the time delay and the magnification ratio of the images. It is then possible to single out the correct reconstruction, and therefore determine the time delay, by checking the consistency of the reconstructed individual lightcurves with the additional interferometric observations. We performed extensive numerical simulations of synthetic lightcurves to investigate the dependence of the performance of this method on various parameters which are involved in the problem. Probably the most promising candidates for applying the method (and also for determining the Hubble constant) are lens systems consisting of multiply imaged compact sources and an Einstein ring, such as B0218+357 from which some of the parameters used for our simulations were adopted.

We develop a new statistical method to reanalyse angular correlations between background QSOs and foreground galaxies that are supposed to be a consequence of dark matter inhomogeneities acting as weak gravitational lenses. The method is based on a weighted average over the galaxy positions and is optimized to distinguish between a random distribution of galaxies around QSOs and a distribution which follows an assumed QSO-galaxy two-point correlation function, by choosing an appropriate weight function. With simulations we demonstrate that this weighted average is slightly more significant than Spearman's rank-order test which was used in previous investigations. In particular, the advantages of the weighted average show up if the two-point correlation function is weak. We then reanalyze the correlation between high-redshift 1-Jansky QSOs and IRAS galaxies, taken from the IRAS Faint Source Catalog; these samples were analyzed previously using Spearman's rank-order test. In agreement with the previous work, we find moderate to strong correlations between these two samples; considering the angular two-point correlation function of these samples, we find a typical scale of order $5'$ from which most of the correlation signal derives. However, the statistical significance of the correlation changes with the redshift slices of the QSO sample one considers. Comparing with simple theoretical estimates of the expected correlation, we find that the signal we derive is considerably stronger than expected. On the other hand, recent direct verifications of the overdensity of matter in the line-of-sight to high-redshift radio QSOs obtained from the shear field around these sources, indicates that the observed association can be attributed to a gravitational lens effect.

The double QSO2345+007 comprises two optical components separated by 7.1\ arcseconds and is the most prominent `dark matter' gravitational lens candidate. Despite being known for more than a decade, optical spectroscopy and imaging have been unable to determine whether this double QSO is a binary QSO or a gravitational lens system. In this note we report a deep VLA observation of this system, yielding a map with a noise level of 8.5\ts $\mu$Jy per beam. We have a $4\sigma$ detection of a radio source within one arcsecond of the optical position of the brighter A-component of the QSO, but no significant detection of any radio counterpart of the B component. Given that the flux ratio in the optical waveband is $\sim$3--4, the gravitational lens hypothesis for this system would predict a radio flux of image B of $\lesssim 10$\ts $\mu$Jy. The non-detection of the B component is thus consistent with, but does not prove, the lens interpretation.

The distortion of images of faint background galaxies by (weak) gravitational lensing can be used to measure the mass distribution of the deflector. The image distortions can be used to define a weighted mean of the mass inside a circular aperture, as was first suggested by Kaiser. The aperture mass can be used to {\it detect} dark matter concentrations. Keeping in mind that wide-field cameras will become increasingly available, this method can be used to search for mass concentrations on wide-field images. To do this, the aperture mass measure is generalized to account for different weighting functions. For each such weighting function, a signal-to-noise ratio can be calculated. For an assumed mass profile of the density concentrations, the weighting function can be chosen such as to maximize the resulting signal-to-noise ratio. Numerical simulations which adopt parameters characteristic of 4-m class telescopes are then used to show that dark halos with a velocity dispersion in excess of $\sim 600$\ts km/s can be reliably detected as significant peaks in the signal-to-noise map. The effects of seeing and an anisotropic PSF are then investigated and shown to be less important than might be feared. It is thus suggested that the method of aperture mass measures developed here can be used to obtain a mass-selected sample of dark halos, in contrast to flux-selected samples. Shear fields around high-redshift bright QSOs as detected by Fort et al. provide a first successful application of this strategy. The simplicity of the method allows its routine application to wide-field images of sufficient depth and image quality.

In a series of previous papers we have considered the reconstruction of the surface mass density of a cluster of galaxies from images of lensed faint background galaxies. We showed that the reconstructed surface mass density is not uniquely determined, but that there exists a global invariance transformation that leaves the shape of the images of the lensed galaxies unchanged. Because of this, only lower limits on the total mass of a cluster can be derived if no further informations besides image ellipticities are used. Throughout these papers we used the simplifying assumption that all sources are at the same redshift. In this paper we account for a redshift distribution of the faint galaxies, and in particular, some of these galaxies can lie in front of the cluster or can be cluster members. We show how the mass distribution of a cluster of galaxies can be obtained from images of these faint galaxies, if the redshift distribution of these galaxies is known. We demonstrate that for the reconstruction of non-critical clusters we need less information on the redshift distribution of the galaxies, i.e., we only need to know two or three moments of the distribution. We show that the mean mass density across the data field is still a free variable, i.e., there remains a global invariance transformation of the resulting mass density field. For non-critical clusters we can derive the transformation explicitly and it is similar to that derived previously for the case of a single redshift of the sources. We discuss several theoretical ideas to break the mass degeneracy; of those considered, we find that only the magnification effect on the number density of galaxy images can be used successfully in practice.

The image distortions of high-redshift galaxies caused by gravitational light deflection of foreground clusters of galaxies can be used to reconstruct the two-dimensional surface mass density of these clusters. We apply an unbiased parameter-free reconstruction technique to the cluster CL0939+4713 (Abell 851), observed with the WFPC2 on board of the HST. We demonstrate that a single deep WFPC2 observation can be used for cluster mass reconstruction despite its small field of view and the irregular shape of the data field (especially for distant clusters). For CL0939, we find a strong correlation between the reconstructed mass distribution and the bright cluster galaxies indicating that mass follows light on average. The detected anti-correlation between the faint galaxies and the reconstructed mass is most likely an effect of the magnification (anti) bias, which was detected previously in the cluster A1689. Because of the high redshift of CL0939 (z_d=0.41), the redshift distribution of the lensed, faint galaxies has to be accounted for in the reconstruction technique. We derive an approximate global transformation for the surface mass density which leaves the mean image ellipticities invariant, resulting in an uncertainty in the normalization of the mass. From the non-negativity of the surface mass density, we derive lower limits on the mass inside the observed field. We can break the invariance transformation for the mass using the magnification effect on the observed number density of the background galaxies. Assuming a mean redshift of 0.8 and a fraction of 15 % (20 %) of cluster galaxies in the observed galaxy sample with 23<R<25.5 we obtain for the mass inside the field approximately 5 (7) x 10**{14}/ h_{50} solar masses which corresponds to M/L ~ 100 (140) h_{50}.

It is pointed out that there are now known four very close pairs of QSOs with separations < 5 arcsec and very different redshifts. Several estimates of the probability that they are accidental configurations range between 10^{-7} and 3.5 x 10^{-3}. We conclude either that this is further evidence that QSOs have significant non-cosmological redshift components, or that thave significant non-cosmological redshift components, or that the pairs must be explained by gravitational lensing.

We present a novel method to recontruct the mass distribution of galaxy clusters from their gravitational lens effect on background galaxies. The method is based on a least-chisquare fit of the two-dimensional gravitational cluster potential. The method combines information from shear and magnification by the cluster lens and is designed to easily incorporate possible additional information. We describe the technique and demonstrate its feasibility with simulated data. Both the cluster morphology and the total cluster mass are well reproduced.

The last decade has seen an enormous increase of activity in the field of gravitational lensing, mainly driven by improvements of observational capabilities. I will review the basics of gravitational lens theory, just enough to understand the rest of this contribution, and will then concentrate on several of the main applications in cosmology. Cluster lensing, and weak lensing, will constitute the main part of this review.

Gravitational light deflection can distort the images of distant sources by its tidal effects. The population of faint blue galaxies is at sufficiently high redshift so that their images are distorted near foreground clusters, with giant luminous arcs being the most spectacular evidence for this effect. Much weaker distortions, however, can observationally be detected by a statistical analysis of the numerous faint galaxy images, as first demonstrated by Tyson, Valdes \& Wenk. This distortion effect can be used as a {\it quantitativetool} for the reconstruction of the surface mass density of galaxy clusters with appropriate redshifts, as was demonstrated by Kaiser \&Squires. They have derived an explicit equation for this surface mass density in terms of its tidal field. The reconstruction formula by Kaiser \& Squires must be modified because of two effects: in its original form it applies only to weak lenses, and hence must be generalized to account for stronger lensing effects. Second, due to the nature of the inversion formula, it produces boundary artefacts (or biases) if applied to real data which are confined to a finite field on the sky. We discuss several possibilities to obtain inversion formulae which are exact for ideal data on a finite data field (CCD). We demonstrate that there exists an infinite number of such finite-field inversion formulae, which differ in their sensitivity to observational effects (such as noise, intrinsic ellipticities of the sources, etc.). We show that, using two simple conditions, one can uniquely specify a finite-field formula which in a well-defined sense minimizes the sensitivity to observational effects. We then use synthetic data to compare the quality of our new reconstruction method with that of previous

The reconstruction of the density profile in clusters of galaxies from the distortion of the images of faint background galaxies is reconsidered. The inversion formula of Kaiser \& Squires is known to provide a quantitative way to perform this reconstruction; however, the practical application of this formula faces two problems of principle (besides problems related to the analysis of the observational data): (1) the shear distribution of a lens cannot be inferred from the distortion of images, but only a combination of shear and surface mass density can be observed. (2) The inversion formula is exact only if one assumes observational data on the whole lens plane, whereas in reality, the size of the data field is limited by the size of the CCD. We have considered a possible solution to the first problem in a previous paper. Here we consider the second problem. It is shown that the application of the inversion formula to a finite data field induces systemmatic boundary effects. An alternative inversion formula is derived, based on some recently published results by Kaiser. We demonstrate, using synthetic data, that this new inversion formula which does not require an extrapolation of the data beyond the observed region, yields results which are comparable with those from the Kaiser \& Squires inversion in their `noise levels', but lack the systemmatic boundary effects.

In a recent paper by Rodrigues-Williams \& Hogan (RH94), a correlation between high-redshift, optically-selected QSOs and Zwicky clusters was reported at a very high significance level. Due to the fairly bright flux threshold of the cluster sample, these correlations cannot be interpreted as being due to an environmental effect of the clusters on the quasar activity. The most likely interpretation employed in RH94 was the effect of gravitational lensing by the foreground clusters, though the required magnification to explain the observed correlations has to be considerably higher than obtained from simple mass models for the clusters. We have repeated the analysis of RH94 using a different QSO sample, namely radio quasars and radio galaxies from the 1-Jy sample. In accordance with RH, we detect a statistically significant correlation between Zwicky clusters and 1-Jy sources with intermediate redshift ($z\approx 1$), but fail to detect significant effects for higher-redshift sources. In addition, we detect a highly significant underdensity of low-redshift radio sources around Zwicky clusters, for which an environmental interpretation seems to be most reasonable. Our result for the overdensity of $z\sim 1$ sources is in good agreement with previous results and can possibly be interpreted as a lensing effect, though we have not tried to quantitatively investigate a lensing model. Other interpretations, such as patchy dust obscuration, can not explain the observed effect and its tendency with redshift.

The weak distortions of high-redshift galaxies caused by gravitational light deflection near clusters of galaxies can be used to reconstruct the projected (two-dimensional) surface mass density of intermediate redshift clusters. This technique, pioneered by Tyson, Valdes & Wenk, and Kaiser & Squires, is reconsidered in the present paper, where we generalize the inversion equation found by Kaiser & Squires (KS) in several respect. Adopting a different smoothing procedure for the discreetly sampled data (individual galaxy images), we effectively reduce the shot noise in the KS procedure. In particular we show that the best density reconstructions are obtained if the smoothing scale is adopted to the `strength of the signal', which yields a better resolution near the center of the cluster where the distortions are strongest. Furthermore, we point out the importance of boundary effects and demonstrate their disastrous impact on rectangular data fields (CCDs) with large side ratio. Most important, however, is the generalization of the KS method to critical clusters, i.e., to such clusters which are capable of producing multiple images and giant luminous arcs. The corresponding modifications of the inversion procedure are severe; in particular, the resulting inversion equation is much more difficult to solve. As we pointed out in a previous paper (Schneider & Seitz), there exists a local degeneracy if the cluster is critical. We have developed an iteration procedure to solve the inversion equation, which we demonstrate to yield a very accurate

The distortion of images of faint, high-redshift galaxies by light deflection at foreground clusters of galaxies can be used to determine the (projected) mass distribution of the clusters. In the case of strong distortions, which lead to arcs in clusters, the position of the arc and/or its radius of curvature yields an estimate for the total mass inside a circle traced out by the arc. Weak distortions, which can be observed to much larger angular separations, can be used to determine the mass profile. In the case of weak distortions, an approximation which identifies the observed distortion with the shear produced by the lens can be made; this (linear) approximation breaks down, however, if one wants to probe the center of the cluster, i.e., approach the region within which giant arcs can be formed. The methods developed hitherto, the most advanced of which is due to Kaiser \& Squires, rely on the linear approximation and thus cannot yield reliable results for regions close to the center of the cluster. The purpose of this paper is to provide a theoretical basis to extend the Kaiser \& Squires method into the nonlinear regime, thereby making it more powerful, since the nonlinear distortions provide strong constraints on the mass profile, which in combination with the weak distortions should yield more reliable cluster inversions. We discuss the statistical properties of the observable image ellipticities and provide several methods to determine the local distortion by the lens from observed galaxy images, some of which rely

We reformulate the transport equation which determines the size, shape and orientation of infinitesimal light beams in arbitrary spacetimes. The behaviour of such light beams near vertices and conjugate points is investigated, with special attention to the singular behaviour of the optical scalars. We then specialize the general transport equation to the case of an approximate metric of an inhomogeneous universe, which is a Friedmann metric `on average' with superposed isolated weak matter inhomogeneities. In a series of well-defined approximations, the equations of gravitational lens theory are derived. Finally, we derive a relative optical focusing equation which describes the focusing of light beams relative to the case that the beam is unaffected by matter inhomogeneities in the universe, from which it follows immediately that no beam can be focused less than one which is unaffected by matter clumps, before it propagates through its first conjugate point.

We announce the detection of correlations on angular scales of $\ga10'$ between optically bright, high-redshift, radio-loud QSOs with diffuse X-ray emission seen by ROSAT in the {\it All-Sky Survey}. These correlations reach significance levels of up to $99.8\%$. A comparison of the results with a sample of control fields, bootstrapping analyses, and Kolmogorov-Smirnov tests provide unambiguous support for the statistical significance of the correlations found. We argue that the detected enhanced diffuse X-ray emission is in the foreground of the QSOs, and that it is probably due to galaxy clusters which magnify the QSOs by their gravitational lensing effect, thereby giving rise to a magnification bias in the background source sample. A comparison of the results presented below with correlations previously found between the same QSO sample and either Lick or IRAS galaxies provides further evidence for this interpretation, and identifies positions in the sky where weak gravitational lensing may be detected by searching for coherent distortions of background galaxy images.

We demonstrate that high-redshift 1-Jansky radio sources are correlated with foreground IRAS galaxies on angular scales of several 10 arc minutes on a statistical significance level of up to 99.8\%. Based on a theoretical analysis published earlier, we propose to interpret these correlations in terms of gravitational lensing by large- scale structures.

We construct a simple elliptical gravitational lens model for the quadruple lens B1422+231 and show that the details of the configuration cannot be easily understood in terms of this model; in particular, the flux ratios of the images are hard to reproduce. This qualitatively verifies the results from a different lens model constructed for the same object by Hogg and Blandford.