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.

Measuring cosmic shear in wide-field imaging surveys requires accurate knowledge of the redshift distribution of all sources. The clustering-redshift technique exploits the angular cross-correlation of a target galaxy sample with unknown redshifts and a reference sample with known redshifts, and is an attractive alternative to colour-based methods of redshift calibration. We test the performance of such clustering redshift measurements using mock catalogues that resemble the Kilo-Degree Survey (KiDS). These mocks are created from the MICE simulation and closely mimic the properties of the KiDS source sample and the overlapping spectroscopic reference samples. We quantify the performance of the clustering redshifts by comparing the cross-correlation results with the true redshift distributions in each of the five KiDS photometric redshift bins. Such a comparison to an informative model is necessary due to the incompleteness of the reference samples at high redshifts. Clustering mean redshifts are unbiased at $|\Delta z|<0.006$ under these conditions. The redshift evolution of the galaxy bias can be reliably mitigated at this level of precision using auto-correlation measurements and self-consistency relations, and will not become a dominant source of systematic error until the arrival of Stage-IV cosmic shear surveys. Using redshift distributions from a direct colour-based estimate instead of the true redshift distributions as a model for comparison with the clustering redshifts increases the biases in the mean to up to $|\Delta z|\sim0.04$. This indicates that the interpretation of clustering redshifts in real-world applications will require more sophisticated (parameterised) models of the redshift distribution in the future. If such better models are available, the clustering-redshift technique promises to be a highly complementary alternative to other methods of redshift calibration.

We simultaneously present constraints on the stellar-to-halo mass relation for central and satellite galaxies through a weak lensing analysis of spectroscopically classified galaxies. Using overlapping data from the 4th data release of the Kilo-Degree Survey (KiDS), and the Galaxy And Mass Assembly survey (GAMA), we find that satellite galaxies are hosted by halo masses that are $0.53 \pm 0.39$ dex (68% confidence, $3\sigma$ detection) smaller than those of central galaxies of the same stellar mass (for a stellar mass of $\log(M_{\star}/M_{\odot}) = 10.6$). This is consistent with galaxy formation models, whereby infalling satellite galaxies are preferentially stripped of their dark matter. We find consistent results with similar uncertainties when comparing constraints from a standard azimuthally averaged galaxy-galaxy lensing analysis and a two-dimensional likelihood analysis of the full shear field. As the latter approach is somewhat biased due to the lens incompleteness, and does not provide any improvement in precision when applied to actual data, we conclude that stacked tangential shear measurements are best-suited for studies of the galaxy-halo connection.

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

The physics of gravity on cosmological scales affects both the rate of assembly of large-scale structure, and the gravitational lensing of background light through this cosmic web. By comparing the amplitude of these different observational signatures, we can construct tests that can distinguish General Relativity from its potential modifications. We use the latest weak gravitational lensing dataset from the Kilo-Degree Survey, KiDS-1000, in conjunction with overlapping galaxy spectroscopic redshift surveys BOSS and 2dFLenS, to perform the most precise existing amplitude-ratio test. We measure the associated E_G statistic with 15-20% errors, in five dz = 0.1 tomographic redshift bins in the range 0.2 < z < 0.7, on projected scales up to 100 Mpc/h. The scale-independence and redshift-dependence of these measurements are consistent with the theoretical expectation of General Relativity in a Universe with matter density Omega_m = 0.27 +/- 0.04. We demonstrate that our results are robust against different analysis choices, including schemes for correcting the effects of source photometric redshift errors, and compare the performance of angular and projected galaxy-galaxy lensing statistics.

We present updated cosmological constraints for the KiDS+VIKING-450 cosmic shear dataset (KV450), estimated using redshift distributions and photometric samples defined using self organising maps (SOMs). Our fiducial analysis finds marginal posterior constraints of $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.716^{+0.043}_{-0.038}$; smaller than, but fully consistent with, previous work using this dataset ($|\Delta S_8| = 0.023$). We analyse additional samples and redshift distributions constructed in three ways: excluding certain spectroscopic surveys during redshift calibration, excluding lower-confidence spectroscopic redshifts in redshift calibration, and considering only photometric sources which are jointly calibrated by at least three spectroscopic surveys. In all cases, the method utilised here proves robust: we find a maximal deviation from our fiducial analysis of $|\Delta S_8| \leq 0.009$ for all samples defined and analysed using our SOM. Our largest shift in $S_8$ is found when calibrating redshift distributions without the DEEP2 spectroscopic subset, where we find $S_8=0.707_{-0.042}^{+0.046}$. This difference with respect to the fiducial is both significantly smaller than, and in the opposite direction to, the equivalent shift from previous work. No sample analysed in this work results in a meaningful positive shift in $S_8$ with respect to our fiducial constraints. These results suggest that our improved cosmological parameter estimates are insensitive to pathological misrepresentation of photometric sources by the spectroscopy used for direct redshift calibration, and therefore that this systematic effect cannot be responsible for the observed difference between $S_8$ estimates made with KV450 and Planck CMB probes.

We present a combined tomographic weak gravitational lensing analysis of the Kilo Degree Survey (KV450) and the Dark Energy Survey (DES-Y1). We homogenize the analysis of these two public cosmic shear datasets by adopting consistent priors and modeling of nonlinear scales, and determine new redshift distributions for DES-Y1 based on deep public spectroscopic surveys. Adopting these revised redshifts results in a $0.8\sigma$ reduction in the DES-inferred value for $S_8$, which decreases to a $0.5\sigma$ reduction when including a systematic redshift calibration error model from mock DES data based on the MICE2 simulation. The combined KV450 + DES-Y1 constraint on $S_8 = 0.762^{+0.025}_{-0.024}$ is in tension with the Planck 2018 constraint from the cosmic microwave background at the level of $2.5\sigma$. This result highlights the importance of developing methods to provide accurate redshift calibration for current and future weak lensing surveys.

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

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.

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.

Dust lanes bisect the plane of a typical edge-on spiral galaxy as a dark optical absorption feature. Their appearance is linked to the gravitational stability of spiral disks; the fraction of edge-on galaxies that displays a dust lane is a direct indicator of the typical vertical balance between gravity and turbulence; a balance struck between the energy input from star-formation and the gravitational pull into the plane of the disk. Based on morphological classifications by the Galaxy~Zoo project on the Kilo-Degree Survey (KiDS) imaging data in the Galaxy and Mass Assembly (GAMA) fields, we explore the relation of dust lanes to the galaxy characteristics, most of which were determined using the {\sc magphys} spectral energy distribution fitting tool: stellar mass, total and specific star-formation rates, and several parameters describing the cold dust component. We find that the fraction of dust lanes does depend on the stellar mass of the galaxy; they start to appear at $M^* \sim 10^9 M_\odot$. A dust lane also implies strongly a dust mass of at least $10^5 M_\odot$, but otherwise does not correlate with cold dust mass parameters of the {\sc magphys} spectral energy distribution analysis, nor is there a link with star-formation rate, specific or total. Dust lane identification does not depend on disk ellipticity (disk thickness) or Sersic profile but correlates with bulge morphology; a round bulge favors dust lane votes. The central component along the line of sight that produces the dust lane is not associated with either one of the components fit by {\sc magphys}, the cold diffuse component or the localized, heated component in HII regions, but a mix of these two.

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

Intrinsic galaxy alignments are a source of bias for weak lensing measurements as well as a tool for understanding galaxy formation and evolution. In this work, we measure the alignment of shapes of satellite galaxies, in galaxy groups, with respect to the brightest group galaxy (BGG), as well as alignments of the BGG shape with the satellite positions, using the highly complete Galaxy And Mass Assembly (GAMA) spectroscopic survey and deep imaging from the Kilo Degree Survey. We control systematic errors with dedicated image simulations and measure accurate shapes using the DEIMOS shape measurement method. We find a significant satellite radial alignment signal, which vanishes at large separations from the BGG. We do not identify any strong trends of the signal with galaxy absolute magnitude or group mass. The alignment signal is dominated by red satellites. We also find that the outer regions of galaxies are aligned more strongly than their inner regions, by varying the radial weight employed during the shape measurement process. This behaviour is evident for both red and blue satellites. BGGs are also found to be aligned with satellite positions, with this alignment being stronger when considering the innermost satellites, using red BGGs and the shape of the outer region of the BGG. Lastly, we measure the global intrinsic alignment signal in the GAMA sample for two different radial weight functions and find no significant difference.

We revisit the performance and biases of the two-dimensional approach to galaxy-galaxy lensing. This method exploits the information for the actual positions and ellipticities of source galaxies, rather than using only the ensemble properties of statistically equivalent samples. We compare the performance of this method with the traditionally used one-dimensional tangential shear signal on a set of mock data that resemble the current state-of-the-art weak lensing surveys. We find that under idealised circumstances the confidence regions of joint constraints for the amplitude and scale parameters of the NFW model in the two-dimensional analysis can be more than three times tighter than the one-dimensional results. Moreover, this improvement depends on the lens number density and it is larger for higher densities. We compare the method against the results from the hydrodynamical EAGLE simulation in order to test for possible biases that might arise due to lens galaxies being missed, and find that the method is able to return unbiased estimates of halo masses when compared to the true properties of the EAGLE galaxies. Because of its advantage in high galaxy density areas, the method is especially suitable for studying the properties of satellite galaxies in clusters of galaxies.

The KiDS Strongly lensed QUAsar Detection project (KiDS-SQuaD) aims at finding as many previously undiscovered gravitational lensed quasars as possible in the Kilo Degree Survey. This is the second paper of this series where we present a new, automatic object classification method based on machine learning technique. The main goal of this paper is to build a catalogue of bright extragalactic objects (galaxies and quasars), from the KiDS Data Release 4, with a minimum stellar contamination, preserving the completeness as much as possible, to then apply morphological methods to select reliable gravitationally lensed (GL) quasar candidates. After testing some of the most used machine learning algorithms, decision trees based classifiers, we decided to use CatBoost, that was specifically trained with the aim of creating a sample of extragalactic sources as clean as possible from stars. We discuss the input data, define the training sample for the classifier, give quantitative estimates of its performances, and finally describe the validation results with Gaia DR2, AllWISE, and GAMA catalogues. We have built and make available to the scientific community the KiDS Bright EXtraGalactic Objects catalogue (KiDS-BEXGO), specifically created to find gravitational lenses. This is made of $\approx6$ millions of sources classified as quasars ($\approx 200\,000$) and galaxies ($\approx 5.7$M), up to $r<22^m$. From this catalog we selected 'Multiplets': close pairs of quasars or galaxies surrounded by at least one quasar, presenting the 12 most reliable gravitationally lensed quasar candidates, to demonstrate the potential of the catalogue, which will be further explored in a forthcoming paper. We compared our search to the previous one, presented in the first paper from this series, showing that employing a machine learning method decreases the stars-contaminators within the GL candidates.

This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and is intended to clarify what LSST observations are needed in order to maximally enhance the combined science output of LSST and Euclid. To facilitate LSST planning we provide a range of possible LSST surveys with clear metrics based on the improvement in the Dark Energy figure of merit (FOM). To provide a quantifiable metric we present five survey options using only between 0.3 and 3.8% of the LSST 10 year survey. We also provide information so that the LSST DDF cadence can possibly be matched to those of \emph{Euclid} in common deep fields, SXDS, COSMOS, CDFS, and a proposed new LSST deep field (near the Akari Deep Field South). Co-coordination of observations from the Large Synoptic Survey Telescope (LSST) and Euclid will lead to a significant number of synergies. The combination of optical multi-band imaging from LSST with high resolution optical and near-infrared photometry and spectroscopy from \emph{Euclid} will not only improve constraints on Dark Energy, but provide a wealth of science on the Milky Way, local group, local large scale structure, and even on first galaxies during the epoch of reionization. A detailed paper has been published on the Dark Energy science case (Rhodes et al.) by a joint LSST/Euclid working group as well as a white paper describing LSST/Euclid/WFIRST synergies (Jain et al.), and we will briefly describe other science cases here. A companion white paper argues the general science case for an extension of the LSST footprint to the north at airmass < 1.8, and we support the white papers for southern extensions of the LSST survey.

WAVES is designed to study the growth of structure, mass and energy on scales of ~1 kpc to ~10 Mpc over a 7 Gyr timeline. On the largest length scales (1-10 Mpc) WAVES will measure the structures defined by groups, filaments and voids, and their emergence over recent times. Comparisons with bespoke numerical simulations will be used to confirm, refine or refute the Cold Dark Matter paradigm. At intermediate length scales (10 kpc-1 Mpc) WAVES will probe the size and mass distribution of galaxy groups, as well as the galaxy merger rates, in order to directly measure the assembly of dark matter halos and stellar mass. On the smallest length scales (1-10 kpc) WAVES will provide accurate distance and environmental measurements to complement high-resolution space-based imaging to study the mass and size evolution of galaxy bulges, discs and bars. In total, WAVES will provide a panchromatic legacy dataset of ~1.6 million galaxies, firmly linking the very low ($z < 0.1$) and intermediate ($z \sim 0.8$) redshift Universe.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

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.

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

Exploiting the full statistical power of future cosmic shear surveys will necessitate improvements to the accuracy with which the gravitational lensing signal is measured. We present a framework for calibrating shear with image simulations that demonstrates the importance of including realistic correlations between galaxy morphology, size and more importantly, photometric redshifts. This realism is essential so that selection and shape measurement biases can be calibrated accurately for a tomographic cosmic shear analysis. We emulate Kilo-Degree Survey (KiDS) observations of the COSMOS field using morphological information from {\it Hubble} Space Telescope imaging, faithfully reproducing the measured galaxy properties from KiDS observations of the same field. We calibrate our shear measurements from lensfit, and find through a range of sensitivity tests that lensfit is robust and unbiased within the allowed 2 per cent tolerance of our study. Our results show that the calibration has to be performed by selecting the tomographic samples in the simulations, consistent with the actual cosmic shear analysis, because the joint distributions of galaxy properties are found to vary with redshift. Ignoring this redshift variation could result in misestimating the shear bias by an amount that exceeds the allowed tolerance. To improve the calibration for future cosmic shear analyses, it will be essential to also correctly account for the measurement of photometric redshifts, which requires simulating multi-band observations.

We directly constrain the non-linear alignment (NLA) model of intrinsic galaxy alignments, analysing the most representative and complete flux-limited sample of spectroscopic galaxies available for cosmic shear surveys. We measure the projected galaxy position-intrinsic shear correlations and the projected galaxy clustering signal using high-resolution imaging from the Kilo Degree Survey (KiDS) overlapping with the GAMA spectroscopic survey, and data from the Sloan Digital Sky Survey. Separating samples by colour, we make no significant detection of blue galaxy alignments, constraining the blue galaxy NLA amplitude $A_{\textrm{IA}}^{\textrm{B}}=0.21^{+0.37}_{-0.36}$ to be consistent with zero. We make robust detections ($\sim9\sigma$) for red galaxies, with $A_{\textrm{IA}}^{\textrm{R}}=3.18^{+0.47}_{-0.46}$, corresponding to a net radial alignment with the galaxy density field, and we find no evidence for any scaling of alignments with galaxy luminosity. We provide informative priors for current and future weak lensing surveys, an improvement over de facto wide priors that allow for unrealistic levels of intrinsic alignment contamination. For a colour-split cosmic shear analysis of the final KiDS survey area, we forecast that our priors will improve the constraining power on $S_{8}$ and the dark energy equation of state $w_{0}$, by up to $62\%$ and $51\%$, respectively. Our results indicate, however, that the modelling of red/blue-split galaxy alignments may be insufficient to describe samples with variable central/satellite galaxy fractions.

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.

The outer regions of galaxies are more susceptible to the tidal interactions that lead to intrinsic alignments of galaxies. The resulting alignment signal may therefore depend on the passband if the colours of galaxies vary spatially. To quantify this, we measured the shapes of galaxies with spectroscopic redshifts from the GAMA survey using deep gri imaging data from the KiloDegree Survey. The performance of the moment-based shape measurement algorithm DEIMOS was assessed using dedicated image simulations, which showed that the ellipticities could be determined with an accuracy better than 1% in all bands. Additional tests for potential systematic errors did not reveal any issues. We measure a significant difference of the alignment signal between the g,r and i-band observations. This difference exceeds the amplitude of the linear alignment model on scales below 2 Mpc/h. Separating the sample into central/satellite and red/blue galaxies, we find that that the difference is dominated by red satellite galaxies.

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

We use the overlap between multiband photometry of the Kilo-Degree Survey (KiDS) and spectroscopic data based on the Sloan Digital Sky Survey (SDSS) and Galaxy And Mass Assembly (GAMA) to infer the colour-magnitude relation of red-sequence galaxies. We then use this inferred relation to select luminous red galaxies (LRGs) in the redshift range of $0.1<z<0.7$ over the entire KiDS Data Release 3 footprint. We construct two samples of galaxies with different constant comoving densities and different luminosity thresholds. The selected red galaxies have photometric redshifts with typical photo-z errors of $\sigma_z \sim 0.014 (1+z)$ that are nearly uniform with respect to observational systematics. This makes them an ideal set of galaxies for lensing and clustering studies. As an example, we use the KiDS-450 cosmic shear catalogue to measure the mean tangential shear signal around the selected LRGs. We detect a significant weak lensing signal for lenses out to $z \sim 0.7$.

We present a new measurement of $E_{\rm G}$, which combines measurements of weak gravitational lensing, real-space galaxy clustering and redshift space distortions. This statistic was proposed as a consistency test of General Relativity (GR) that is insensitive to linear, deterministic galaxy bias and the matter clustering amplitude. We combine deep imaging data from KiDS with overlapping spectroscopy from 2dFLenS, BOSS DR12 and GAMA and find $E_{\rm G}(\overline{z}=0.267)=0.43 \pm 0.13$ (GAMA), $E_{\rm G}(\overline{z}=0.305)=0.27 \pm 0.08$ (LOWZ+2dFLOZ) and $E_{\rm G}(\overline{z}=0.554)=0.26 \pm 0.07$ (CMASS+2dFHIZ). We demonstrate that the existing tension in the value of the matter density parameter hinders the robustness of this statistic as solely a test of GR. We find that our $E_{\rm G}$ measurements, as well as existing ones in the literature, favour a lower matter density cosmology than the Cosmic Microwave Background. For a flat $\Lambda$CDM Universe and assuming GR, we find $\Omega_{\rm m}(z=0)=0.25\pm0.03$. With this paper we publicly release the 2dFLenS dataset at: \url{http://2dflens.swin.edu.au}.

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.

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

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.

In the hierarchical two-phase formation scenario, the extended halos of early type galaxies (ETGs) are expected to have different physical properties from those of the galaxies' central regions. This work aims at characterizing the kinematic properties of ETG halos using planetary nebulae (PNe) as tracers, which allow us to overcome the limitations of absorption line spectroscopy of continuum at low surface brightness. We present two-dimensional velocity and velocity dispersion fields for 33 ETGs, including both fast (FRs) and slow rotators (SRs), making this the largest kinematic survey to-date of extragalactic PNe. The velocity fields extend out to a median 5.6 effective radii (Re), with a range [3Re-13Re] for the PN.S ETGs. We complemented the PN kinematics with absorption line data from the literature. We find that ETGs typically show a kinematic transition between inner regions and halo. Estimated transition radii in units of Re anticorrelate with stellar mass. SRs have increased but still modest rotational support at large radii, while most of the FRs show a decrease in rotation, due to the fading of the stellar disk in the outer, more slowly rotating spheroid. 30% of the FRs are dominated by rotation also at large radii. Most ETGs have flat or slightly falling halo velocity dispersion profiles, but 15% of the sample have steeply falling profiles. 40% of the FRs show kinematic twists, misalignments, or rotation along two axes, indicating that they turn from oblate at the center to triaxial in the halo, consistently with the distribution of isophote twists from extended photometry. ETGs have more diverse kinematic properties in their halos than in the central regions. The observed kinematic transition to the halo and its dependence on stellar mass is consistent with LambdaCDM simulations and supports a two-phase formation scenario. [abridged]

We have obtained structural parameters of about 340,000 galaxies from the Kilo Degree Survey (KiDS) in 153 square degrees of data release 1, 2 and 3. We have performed a seeing convolved 2D single S\'ersic fit to the galaxy images in the 4 photometric bands (u, g, r, i) observed by KiDS, by selecting high signal-to-noise ratio (S/N > 50) systems in every bands. We have classified galaxies as spheroids and disc-dominated by combining their spectral energy distribution properties and their S\'ersic index. Using photometric redshifts derived from a machine learning technique, we have determined the evolution of the effective radius, \Re\ and stellar mass, \mst, versus redshift, for both mass complete samples of spheroids and disc-dominated galaxies up to z ~ 0.6. Our results show a significant evolution of the structural quantities at intermediate redshift for the massive spheroids ($\mbox{Log}\ M_*/M_\odot>11$, Chabrier IMF), while almost no evolution has found for less massive ones ($\mbox{Log}\ M_*/M_\odot < 11$). On the other hand, disc dominated systems show a milder evolution in the less massive systems ($\mbox{Log}\ M_*/M_\odot < 11$) and possibly no evolution of the more massive systems. These trends are generally consistent with predictions from hydrodynamical simulations and independent datasets out to redshift z ~ 0.6, although in some cases the scatter of the data is large to drive final conclusions. These results, based on 1/10 of the expected KiDS area, reinforce precedent finding based on smaller statistical samples and show the route toward more accurate results, expected with the the next survey releases.

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

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.

We present a machine-learning photometric redshift analysis of the Kilo-Degree Survey Data Release 3, using two neural-network based techniques: ANNz2 and MLPQNA. Despite limited coverage of spectroscopic training sets, these ML codes provide photo-zs of quality comparable to, if not better than, those from the BPZ code, at least up to zphot<0.9 and r<23.5. At the bright end of r<20, where very complete spectroscopic data overlapping with KiDS are available, the performance of the ML photo-zs clearly surpasses that of BPZ, currently the primary photo-z method for KiDS. Using the Galaxy And Mass Assembly (GAMA) spectroscopic survey as calibration, we furthermore study how photo-zs improve for bright sources when photometric parameters additional to magnitudes are included in the photo-z derivation, as well as when VIKING and WISE infrared bands are added. While the fiducial four-band ugri setup gives a photo-z bias $\delta z=-2e-4$ and scatter $\sigma_z<0.022$ at mean z = 0.23, combining magnitudes, colours, and galaxy sizes reduces the scatter by ~7% and the bias by an order of magnitude. Once the ugri and IR magnitudes are joined into 12-band photometry spanning up to 12 $\mu$, the scatter decreases by more than 10% over the fiducial case. Finally, using the 12 bands together with optical colours and linear sizes gives $\delta z<4e-5$ and $\sigma_z<0.019$. This paper also serves as a reference for two public photo-z catalogues accompanying KiDS DR3, both obtained using the ANNz2 code. The first one, of general purpose, includes all the 39 million KiDS sources with four-band ugri measurements in DR3. The second dataset, optimized for low-redshift studies such as galaxy-galaxy lensing, is limited to r<20, and provides photo-zs of much better quality than in the full-depth case thanks to incorporating optical magnitudes, colours, and sizes in the GAMA-calibrated photo-z derivation.

We aim to investigate the galaxy environment in GAMA Galaxy Groups Catalogue (G3C) using a volume-limited galaxy sample from the Kilo Degree Survey Data Release 3. The k-Nearest Neighbour technique is adapted to take into account the probability density functions (PDFs) of photometric redshifts in our calculations. This algorithm was tested on simulated KiDS tiles, showing its capability of recovering the relation between galaxy colour, luminosity and local environment. The characterization of the galaxy environment in G3C groups shows systematically steeper density contrasts for more massive groups. The red galaxy fraction gradients in these groups is evident for most of group mass bins. The density contrast of red galaxies is systematically higher at group centers when compared to blue galaxy ones. In addition, distinct group center definitions are used to show that our results are insensitive to center definitions. These results confirm the galaxy evolution scenario which environmental mechanisms are responsible for a slow quenching process as galaxies fall into groups and clusters, resulting in a smooth observed colour gradients in galaxy systems.

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.

The decomposition of the rotation curve of galaxies into contribution from the disc and dark halo remains uncertain and depends on the adopted mass to light ratio (M/L) of the disc. Given the vertical velocity dispersion of stars and disc scale height, the disc surface mass density and hence the M/L can be estimated. We address a conceptual problem with previous measurements of the scale height and dispersion. When using this method, the dispersion and scale height must refer to the same population of stars. The scale height is obtained from near-IR studies of edge-on galaxies and is weighted towards older kinematically hotter stars, whereas the dispersion obtained from integrated light in the optical bands includes stars of all ages. We aim to extract the dispersion for the hotter stars, so that it can then be used with the correct scale height to obtain the disc surface mass density. We use a sample of planetary nebulae (PNe) as dynamical tracers in the face-on galaxy NGC 628. We extract two different dispersions from its velocity histogram -- representing the older and younger PNe. We also present complementary stellar absorption spectra in the inner regions of this galaxy and use a direct pixel fitting technique to extract the two components. Our analysis concludes that previous studies, which do not take account of the young disc, underestimate the disc surface mass density by a factor of ~ 2. This is sufficient to make a maximal disc for NGC 628 appear like a submaximal disc.

We study the scaling relations between the baryonic content and total mass of groups of galaxies, as these systems provide a unique way to examine the role of non-gravitational processes in structure formation. Using Planck and ROSAT data, we conduct detailed comparisons of the stacked thermal Sunyaev-Zel'dovich (tSZ) effect and X-ray scaling relations of galaxy groups found in the the Galaxy And Mass Assembly (GAMA) survey and the BAHAMAS hydrodynamical simulation. We use weak gravitational lensing data from the Kilo Degree Survey (KiDS) to determine the average halo mass of the studied systems. We analyse the simulation in the same way, using realistic weak lensing, X-ray, and tSZ synthetic observations. Furthermore, to keep selection biases under control, we employ exactly the same galaxy selection and group identification procedures to the observations and simulation. Applying this careful comparison, we find that the simulations are in agreement with the observations, particularly with regards to the scaling relations of the lensing and tSZ results. This finding demonstrates that hydrodynamical simulation have reached the level of realism that is required to interpret observational survey data and study the baryon physics within dark matter haloes, where analytical modelling is challenging. Finally, using simulated data, we demonstrate that our observational processing of the X-ray and tSZ signals is free of significant biases. We find that our optical group selection procedure has, however, some room for improvement.

We describe data release 3 (DR3) of the Galaxy And Mass Assembly (GAMA) survey. The GAMA survey is a spectroscopic redshift and multi-wavelength photometric survey in three equatorial regions each of 60.0 deg^2 (G09, G12, G15), and two southern regions of 55.7 deg^2 (G02) and 50.6 deg^2 (G23). DR3 consists of: the first release of data covering the G02 region and of data on H-ATLAS sources in the equatorial regions; and updates to data on sources released in DR2. DR3 includes 154809 sources with secure redshifts across four regions. A subset of the G02 region is 95.5% redshift complete to r<19.8 over an area of 19.5 deg^2, with 20086 galaxy redshifts, that overlaps substantially with the XXL survey (X-ray) and VIPERS (redshift survey). In the equatorial regions, the main survey has even higher completeness (98.5%), and spectra for about 75% of H-ATLAS filler targets were also obtained. This filler sample extends spectroscopic redshifts, for probable optical counterparts to H-ATLAS sub-mm sources, to 0.8 mag deeper (r<20.6) than the GAMA main survey. There are 25814 galaxy redshifts for H-ATLAS sources from the GAMA main or filler surveys. GAMA DR3 is available at the survey website (www.gama-survey.org/dr3/).

Weak gravitational lensing analyses are fundamentally limited by the intrinsic, non-Gaussian distribution of galaxy shapes. We explore alternative statistics for samples of ellipticity measurements that are unbiased, efficient, and robust. We take the non-linear mapping of gravitational shear and the effect of noise into account. We then discuss how the distribution of individual galaxy shapes in the observed field of view can be modeled by fitting Fourier modes to the shear pattern directly. We simulated samples of galaxy ellipticities, using both theoretical distributions and real data for ellipticities and noise. We determined the possible bias $\Delta e$, the efficiency $\eta$ and the robustness of the least absolute deviations, the biweight, and the convex hull peeling estimators, compared to the canonical weighted mean. Using these statistics for regression, we have shown the applicability of direct Fourier mode fitting. These estimators can be unbiased in the absence of noise, and decrease noise bias by more than $\sim 30\%$. The convex hull peeling estimator distribution is centered around the underlying shear, and its bias least affected by noise. The least absolute deviations estimator to be the most efficient estimator in almost all cases, except in the Gaussian case, where it's still competitive ($0.83<\eta <5.1$) and therefore robust. These results hold when fitting Fourier modes, where amplitudes of variation in ellipticity are determined to the order of $10^{-3}$. The peak of the ellipticity distribution is a direct tracer of the underlying shear and unaffected by noise, and we have shown that estimators that are sensitive to a central cusp perform more efficiently, potentially reducing uncertainties by more than $50\%$ and significantly decreasing noise bias.

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.

The search for minor bodies in the solar system promises insights into its formation history. Wide imaging surveys offer the opportunity to serendipitously discover and identify these traces of planetary formation and evolution. We aim to present a method to acquire position, photometry, and proper motion measurements of solar system objects in surveys using dithered image sequences. The application of this method on the Kilo-Degree Survey is demonstrated. Optical images of 346 square degree fields of the sky are searched in up to four filters using the AstrOmatic software suite to reduce the pixel to catalog data. The solar system objects within the acquired sources are selected based on a set of criteria depending on their number of observation, motion, and size. The Virtual Observatory SkyBoT tool is used to identify known objects. We observed 20,221 SSO candidates, with an estimated false-positive content of less than 0.05%. Of these SSO candidates, 53.4% are identified by SkyBoT. KiDS can detect previously unknown SSOs because of its depth and coverage at high ecliptic latitude, including parts of the Southern Hemisphere. Thus we expect the large fraction of the 46.6% of unidentified objects to be truly new SSOs. Our method is applicable to a variety of dithered surveys such as DES, LSST, and Euclid. It offers a quick and easy-to-implement search for solar system objects. SkyBoT can then be used to estimate the completeness of the recovered sample.

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.

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

In recent years, multiple studies have reported substantial populations of large, low-surface-brightness galaxies in local galaxy clusters. Various theories that aim to explain the presence of such ultra-diffuse galaxies (UDGs) have since been proposed. A key question that will help to differentiate between models is whether UDGs have counterparts in lower-mass host haloes, and what their abundance as a function of halo mass is. In this study we extend our previous study of UDGs in galaxy clusters to galaxy groups. We measure the abundance of UDGs in 325 spectroscopically-selected groups from the Galaxy And Mass Assembly (GAMA) survey. We make use of the overlapping imaging from the ESO Kilo-Degree Survey (KiDS), from which we can identify galaxies with mean surface brightnesses within their effective radii down to ~25.5 mag arcsec$^{-2}$ in the r-band. We are able to measure a significant overdensity of UDGs (with sizes r_eff > 1.5 kpc) in galaxy groups down to M200=10^12 Msun, a regime where approximately only 1 in 10 groups contains a UDG that we can detect. We combine measurements of the abundance of UDGs in haloes that cover three orders of magnitude in halo mass, finding that their numbers scale quite steeply with halo mass; N_UDG (R<R200) $\propto$ M200^(1.11+/-0.07). To better interpret this, we also measure the mass-richness relation for brighter galaxies down to $M^*_r$+2.5 in the same GAMA groups, and find a much shallower relation of N_Bright (R<R200) $\propto$ M200^(0.78+/-0.05). This shows that UDGs are relatively more abundant, compared to bright galaxies, in massive clusters than in groups. We discuss implications, but whether this difference is related to a higher destruction rate of UDGs in groups, or whether massive haloes have a positive effect on their formation, is not yet clear.

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.

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 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 measure the cross-correlation between Fermi-LAT gamma-ray photons and over 1000 deg$^2$ of weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Red Cluster Sequence Lensing Survey (RCSLenS), and the Kilo Degree Survey (KiDS). We present the first measurement of tomographic weak lensing cross-correlations and the first application of spectral binning to cross-correlations between gamma rays and weak lensing. The measurements are performed using an angular power spectrum estimator while the covariance is estimated using an analytical prescription. We verify the accuracy of our covariance estimate by comparing it to two internal covariance estimators. Based on the non-detection of a cross-correlation signal, we derive constraints on weakly interacting massive particle (WIMP) dark matter. We compute exclusion limits on the dark matter annihilation cross-section $\langle\sigma_\rm{ann} v \rangle$, decay rate $\Gamma_\rm{dec}$, and particle mass $m_\rm{DM}$. We find that in the absence of a cross-correlation signal, tomography does not significantly improve the constraining power of the analysis. Assuming a strong contribution to the gamma-ray flux due to small-scale clustering of dark matter and accounting for known astrophysical sources of gamma rays, we exclude the thermal relic cross-section for particle masses of $m_\rm{DM}\lesssim 20$ GeV.

We describe the selection of galaxies targeted in eight low redshift clusters (APMCC0917, A168, A4038, EDCC442, A3880, A2399, A119 and A85; $0.029 < z < 0.058$) as part of the Sydney-AAO Multi-Object integral field Spectrograph Galaxy Survey (SAMI-GS). We have conducted a redshift survey of these clusters using the AAOmega multi-object spectrograph on the 3.9m Anglo-Australian Telescope. The redshift survey is used to determine cluster membership and to characterise the dynamical properties of the clusters. In combination with existing data, the survey resulted in 21,257 reliable redshift measurements and 2899 confirmed cluster member galaxies. Our redshift catalogue has a high spectroscopic completeness ($\sim 94\%$) for $r_{\rm petro} \leq 19.4$ and clustercentric distances $R< 2\rm{R}_{200}$. We use the confirmed cluster member positions and redshifts to determine cluster velocity dispersion, $\rm{R}_{200}$, virial and caustic masses, as well as cluster structure. The clusters have virial masses $14.25 \leq {\rm log }({\rm M}_{200}/\rm{M}_{\odot}) \leq 15.19$. The cluster sample exhibits a range of dynamical states, from relatively relaxed-appearing systems, to clusters with strong indications of merger-related substructure. Aperture- and PSF-matched photometry are derived from SDSS and VST/ATLAS imaging and used to estimate stellar masses. These estimates, in combination with the redshifts, are used to define the input target catalogue for the cluster portion of the SAMI-GS. The primary SAMI-GS cluster targets have $R< \rm{R}_{200}$, velocities $|v_{\rm pec}| < 3.5\sigma_{200}$ and stellar masses $9.5 \leq {\rm log(M}^*_{approx}/\rm{M}_{\odot}) \leq 12$. Finally, we give an update on the SAMI-GS progress for the cluster regions.

Current optical imaging surveys for cosmology are covering large areas of sky. To exploit the statistical power of these surveys for weak lensing measurements requires shape measurement methods with subpercent systematic errors. We introduce a new weak lensing shear measurement algorithm, Shear Nulling after PSF Gaussianisation (SNAPG), designed to avoid the noise biases that affect most other methods. SNAPG operates on images that have been convolved with a kernel that renders the Point Spread Function (PSF) a circular Gaussian, and uses weighted second moments of the sources. The response of such second moments to a shear of the pre-seeing galaxy image can be predicted analytically, allowing us to construct a shear nulling scheme that finds the shear parameters for which the observed galaxies are consistent with an unsheared, isotropically oriented population of sources. The inverse of this nulling shear is then an estimate of the gravitational lensing shear. We identify the uncertainty of the estimated centre of each galaxy as the source of noise bias, and incorporate an approximate estimate of the centroid covariance into the scheme. We test the method on extensive suites of simulated galaxies of increasing complexity, and find that it is capable of shear measurements with multiplicative bias below 0.5 percent.

Verlinde (2016) proposed that the observed excess gravity in galaxies and clusters is the consequence of Emergent Gravity (EG). In this theory the standard gravitational laws are modified on galactic and larger scales due to the displacement of dark energy by baryonic matter. EG gives an estimate of the excess gravity (described as an apparent dark matter density) in terms of the baryonic mass distribution and the Hubble parameter. In this work we present the first test of EG using weak gravitational lensing, within the regime of validity of the current model. Although there is no direct description of lensing and cosmology in EG yet, we can make a reasonable estimate of the expected lensing signal of low redshift galaxies by assuming a background LambdaCDM cosmology. We measure the (apparent) average surface mass density profiles of 33,613 isolated central galaxies, and compare them to those predicted by EG based on the galaxies' baryonic masses. To this end we employ the ~180 square degrees overlap of the Kilo-Degree Survey (KiDS) with the spectroscopic Galaxy And Mass Assembly (GAMA) survey. We find that the prediction from EG, despite requiring no free parameters, is in good agreement with the observed galaxy-galaxy lensing profiles in four different stellar mass bins. Although this performance is remarkable, this study is only a first step. Further advancements on both the theoretical framework and observational tests of EG are needed before it can be considered a fully developed and solidly tested theory.

After its first implementation in 2003 the Astro-WISE technology has been rolled out in several European countries and is used for the production of the KiDS survey data. In the multi-disciplinary Target initiative this technology, nicknamed WISE technology, has been further applied to a large number of projects. Here, we highlight the data handling of other astronomical applications, such as VLT-MUSE and LOFAR, together with some non-astronomical applications such as the medical projects Lifelines and GLIMPS, the MONK handwritten text recognition system, and business applications, by amongst others, the Target Holding. We describe some of the most important lessons learned and describe the application of the data-centric WISE type of approach to the Science Ground Segment of the Euclid satellite.

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.

We present a new training set for estimating empirical photometric redshifts of galaxies, which was created as part of the 2dFLenS project. This training set is located in a 700 sq deg area of the KiDS South field and is randomly selected and nearly complete at r<19.5. We investigate the photometric redshift performance obtained with ugriz photometry from VST-ATLAS and W1/W2 from WISE, based on several empirical and template methods. The best redshift errors are obtained with kernel-density estimation, as are the lowest biases, which are consistent with zero within statistical noise. The 68th percentiles of the redshift scatter for magnitude-limited samples at r<(15.5, 17.5, 19.5) are (0.014, 0.017, 0.028). In this magnitude range, there are no known ambiguities in the colour-redshift map, consistent with a small rate of redshift outliers. In the fainter regime, the KDE method produces p(z) estimates per galaxy that represent unbiased and accurate redshift frequency expectations. The p(z) sum over any subsample is consistent with the true redshift frequency plus Poisson noise. Further improvements in redshift precision at r<20 would mostly be expected from filter sets with narrower passbands to increase the sensitivity of colours to small changes in redshift.

We develop a statistical estimator to infer the redshift probability distribution of a photometric sample of galaxies from its angular cross-correlation in redshift bins with an overlapping spectroscopic sample. This estimator is a minimum variance weighted quadratic function of the data: a quadratic estimator. This extends and modifies the methodology presented by McQuinn & White (2013). The derived source redshift distribution is degenerate with the source galaxy bias, which must be constrained via additional assumptions. We apply this estimator to constrain source galaxy redshift distributions in the Kilo-Degree imaging survey through cross-correlation with the spectroscopic 2-degree Field Lensing Survey, presenting results first as a binned step-wise distribution in the range z < 0.8, and then building a continuous distribution using a Gaussian process model. We demonstrate the robustness of our methodology using mock catalogues constructed from N-body simulations, and comparisons with other techniques for inferring the redshift distribution.

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 Red-sequence Cluster Lensing Survey (RCSLenS), an application of the methods developed for the Canada France Hawaii Telescope Lensing Survey (CFHTLenS) to the ~785deg$^2$, multi-band imaging data of the Red-sequence Cluster Survey 2 (RCS2). This project represents the largest public, sub-arcsecond seeing, multi-band survey to date that is suited for weak gravitational lensing measurements. With a careful assessment of systematic errors in shape measurements and photometric redshifts we extend the use of this data set to allow cross-correlation analyses between weak lensing observables and other data sets. We describe the imaging data, the data reduction, masking, multi-colour photometry, photometric redshifts, shape measurements, tests for systematic errors, and a blinding scheme to allow for more objective measurements. In total we analyse 761 pointings with r-band coverage, which constitutes our lensing sample. Residual large-scale B-mode systematics prevent the use of this shear catalogue for cosmic shear science. The effective number density of lensing sources over an unmasked area of 571.7deg$^2$ and down to a magnitude limit of r~24.5 is 8.1 galaxies per arcmin$^2$ (weighted: 5.5 arcmin$^{-2}$) distributed over 14 patches on the sky. Photometric redshifts based on 4-band griz data are available for 513 pointings covering an unmasked area of 383.5 deg$^2$ We present weak lensing mass reconstructions of some example clusters as well as the full survey representing the largest areas that have been mapped in this way. All our data products are publicly available through CADC at http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/community/rcslens/query.html in a format very similar to the CFHTLenS data release.

We present the 2-degree Field Lensing Survey (2dFLenS), a new galaxy redshift survey performed at the Anglo-Australian Telescope. 2dFLenS is the first wide-area spectroscopic survey specifically targeting the area mapped by deep-imaging gravitational lensing fields, in this case the Kilo-Degree Survey. 2dFLenS obtained 70,079 redshifts in the range z < 0.9 over an area of 731 sq deg, and is designed to extend the datasets available for testing gravitational physics and promote the development of relevant algorithms for joint imaging and spectroscopic analysis. The redshift sample consists first of 40,531 Luminous Red Galaxies (LRGs), which enable analyses of galaxy-galaxy lensing, redshift-space distortion, and the overlapping source redshift distribution by cross-correlation. An additional 28,269 redshifts form a magnitude-limited (r < 19.5) nearly-complete sub-sample, allowing direct source classification and photometric-redshift calibration. In this paper, we describe the motivation, target selection, spectroscopic observations, and clustering analysis of 2dFLenS. We use power spectrum multipole measurements to fit the redshift-space distortion parameter of the LRG sample in two redshift ranges 0.15 < z < 0.43 and 0.43 < z < 0.7 as beta = 0.49 +/- 0.15 and beta = 0.26 +/- 0.09, respectively. These values are consistent with those obtained from LRGs in the Baryon Oscillation Spectroscopic Survey. 2dFLenS data products will be released via our website http://2dflens.swin.edu.au.

With the aim of paving the road for future accurate astrometry with MICADO at the European-ELT, we performed an astrometric study using two different but complementary approaches to investigate two critical components that contribute to the total astrometric accuracy. First, we tested the predicted improvement in the astrometric measurements with the use of an atmospheric dispersion corrector (ADC) by simulating realistic images of a crowded Galactic globular cluster. We found that the positional measurement accuracy should be improved by up to ~2 mas with the ADC, making this component fundamental for high-precision astrometry. Second, we analysed observations of a globular cluster taken with the only currently available Multi-Conjugate Adaptive Optics assisted camera, GeMS/GSAOI at Gemini South. Making use of previously measured proper motions of stars in the field of view, we were able to model the distortions affecting the stellar positions. We found that they can be as large as ~200 mas, and that our best model corrects them to an accuracy of ~1 mas. We conclude that future astrometric studies with MICADO requires both an ADC and an accurate modelling of distortions to the field of view, either through an a-priori calibration or an a-posteriori correction.

MICADO will equip the E-ELT with a first light capability for diffraction limited imaging at near-infrared wavelengths. The instrument's observing modes focus on various flavours of imaging, including astrometric, high contrast, and time resolved. There is also a single object spectroscopic mode optimised for wavelength coverage at moderately high resolution. This contribution provides an overview of the key functionality of the instrument, outlining the scientific rationale for its observing modes. The interface between MICADO and the adaptive optics system MAORY that feeds it is summarised. The design of the instrument is discussed, focussing on the optics and mechanisms inside the cryostat, together with a brief overview of the other key sub-systems.

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 abundance of compact, massive, early-type galaxies (ETGs) provides important constraints to galaxy formation scenarios. Thanks to the area covered, depth, excellent spatial resolution and seeing, the ESO Public optical Kilo Degree Survey (KiDS), carried out with the VLT Survey Telescope (VST), offers a unique opportunity to conduct a complete census of the most compact galaxies in the Universe. This paper presents a first census of such systems from the first 156 square degrees of KiDS. Our analysis relies on g-, r-, and i-band effective radii ($R_{\rm e}$), derived by fitting galaxy images with PSF-convolved S\'ersic models, high-quality photometric redshifts, $z_{\rm phot}$, estimated from machine learning techniques, and stellar masses, $M_{\rm \star}$, calculated from KiDS aperture photometry. After massiveness ($M_{\rm \star} > 8 \times 10^{10}\, \rm M_{\odot}$) and compactness ($R_{\rm e} < 1.5 \, \rm kpc$ in g-, r- and i-bands) criteria are applied, a visual inspection of the candidates plus near-infrared photometry from VIKING-DR1 are used to refine our sample. The final catalog, to be spectroscopically confirmed, consists of 92 systems in the redshift range $z \sim 0.2-0.7$. This sample, which we expect to increase by a factor of ten over the total survey area, represents the first attempt to select massive super-compact ETGs (MSCGs) in KiDS. We investigate the impact of redshift systematics in the selection, finding that this seems to be a major source of contamination in our sample. A preliminary analysis shows that MSCGs exhibit negative internal colour gradients, consistent with a passive evolution of these systems. We find that the number density of MSCGs is only mildly consistent with predictions from simulations at $z>0.2$, while no such system is found at $z < 0.2$.

We measure the weak gravitational lensing shear power spectra and their cross-power in two photometric redshift bins from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). The measurements are performed directly in multipole space in terms of adjustable band powers. For the extraction of the band powers from the data we have implemented and extended a quadratic estimator, a maximum likelihood method that allows us to readily take into account irregular survey geometries, masks, and varying sampling densities. We find the 68 per cent credible intervals in the $\sigma_8$-$\Omega_{\rm m}$-plane to be marginally consistent with results from $Planck$ for a simple five parameter $\Lambda$CDM model. For the projected parameter $S_8 \equiv \sigma_8(\Omega_{\rm m}/0.3)^{0.5}$ we obtain a best-fitting value of $S_8 = 0.768_{-0.039}^{+0.045}$. This constraint is consistent with results from other CFHTLenS studies as well as the Dark Energy Survey. Our most conservative model, including modifications to the power spectrum due to baryon feedback and marginalization over photometric redshift errors, yields an upper limit on the total mass of three degenerate massive neutrinos of $\Sigma m_\nu < 4.53 \, {\rm eV}$ at 95 per cent credibility, while a Bayesian model comparison does not favour any model extension beyond a simple five parameter $\Lambda$CDM model. Combining the shear likelihood with $Planck$ breaks the $\sigma_8$-$\Omega_{\rm m}$-degeneracy and yields $\sigma_8=0.818 \pm 0.013$ and $\Omega_{\rm m} = 0.300 \pm 0.011$ which is fully consistent with results from $Planck$ alone.

We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models.

Dark energy may be the first sign of new fundamental physics in the Universe, taking either a physical form or revealing a correction to Einsteinian gravity. Weak gravitational lensing and galaxy peculiar velocities provide complementary probes of General Relativity, and in combination allow us to test modified theories of gravity in a unique way. We perform such an analysis by combining measurements of cosmic shear tomography from the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) with the growth of structure from the WiggleZ Dark Energy Survey and the Six-degree-Field Galaxy Survey (6dFGS), producing the strongest existing joint constraints on the metric potentials that describe general theories of gravity. For scale-independent modifications to the metric potentials which evolve linearly with the effective dark energy density, we find present-day cosmological deviations in the Newtonian potential and curvature potential from the prediction of General Relativity to be (Delta Psi)/Psi = 0.05 \pm 0.25 and (Delta Phi)/Phi = -0.05 \pm 0.3 respectively (68 per cent CL).

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 present the results of our first year of quasar search in the on-going ESO public Kilo Degree Survey (KiDS) and VISTA Kilo-Degree Infrared Galaxy (VIKING) surveys. These surveys are among the deeper wide-field surveys that can be used to uncovered large numbers of z~6 quasars. This allows us to probe a more common population of z~6 quasars that is fainter than the well-studied quasars from the main Sloan Digital Sky Survey. From this first set of combined survey catalogues covering ~250 deg^2 we selected point sources down to Z_AB=22 that had a very red i-Z (i-Z>2.2) colour. After follow-up imaging and spectroscopy, we discovered four new quasars in the redshift range 5.8<z<6.0. The absolute magnitudes at a rest-frame wavelength of 1450 A are between -26.6 < M_1450 < -24.4, confirming that we can find quasars fainter than M^*, which at z=6 has been estimated to be between M^*=-25.1 and M^*=-27.6. The discovery of 4 quasars in 250 deg^2 of survey data is consistent with predictions based on the z~6 quasar luminosity function. We discuss various ways to push the candidate selection to fainter magnitudes and we expect to find about 30 new quasars down to an absolute magnitude of M_1450=-24. Studying this homogeneously selected faint quasar population will be important to gain insight into the onset of the co-evolution of the black holes and their stellar hosts.

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.

We present a weak gravitational lensing analysis of supergroup SG1120$-$1202, consisting of four distinct X-ray-luminous groups, that will merge to form a cluster comparable in mass to Coma at $z=0$. These groups lie within a projected separation of 1 to 4 Mpc and within $\Delta v=550$ km s$^{-1}$ and form a unique protocluster to study the matter distribution in a coalescing system. Using high-resolution {\em HST}/ACS imaging, combined with an extensive spectroscopic and imaging data set, we study the weak gravitational distortion of background galaxy images by the matter distribution in the supergroup. We compare the reconstructed projected density field with the distribution of galaxies and hot X-ray emitting gas in the system and derive halo parameters for the individual density peaks. We show that the projected mass distribution closely follows the locations of the X-ray peaks and associated brightest group galaxies. One of the groups that lies at slightly lower redshift ($z\approx 0.35$) than the other three groups ($z\approx 0.37$) is X-ray luminous, but is barely detected in the gravitational lensing signal. The other three groups show a significant detection (up to $5 \sigma$ in mass), with velocity dispersions between $355^{+55}_{-70}$ and $530^{+45}_{-55}$ km s$^{-1}$ and masses between $0.8^{+0.4}_{-0.3} \times 10^{14}$ and $1.6^{+0.5}_{-0.4}\times 10^{14} h^{-1} M_{\odot}$, consistent with independent measurements. These groups are associated with peaks in the galaxy and gas density in a relatively straightforward manner. Since the groups show no visible signs of interaction, this supports the picture that we are catching the groups before they merge into a cluster.

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.

We present weak lensing and X-ray analysis of 12 low mass clusters from the CFHTLenS and XMM-CFHTLS surveys. We combine these systems with high-mass systems from CCCP and low-mass systems from COSMOS to obtain a sample of 70 systems, spanning over two orders of magnitude in mass. We measure core-excised Lx-Tx, M-Lx and M-Tx scaling relations and include corrections for observational biases. By providing fully bias corrected relations, we give the current limitations for Lx and Tx as cluster mass proxies. We demonstrate that Tx benefits from a significantly lower intrinsic scatter at fixed mass than Lx. By studying the residuals of the bias corrected relations, we show for the first time using weak lensing masses that galaxy groups seem more luminous and warmer for their mass than clusters. This implies a steepening of the M-Lx and M-Tx relations at low masses. We verify the inferred steepening using a different high mass sample from the literature and show that variance between samples is the dominant effect leading to discrepant scaling relations. We divide our sample into subsamples of merging and relaxed systems, and find that mergers may have enhanced scatter in lensing measurements, most likely due to stronger triaxiality and more substructure. For the Lx-Tx relation, which is unaffected by lensing measurements, we find the opposite trend in scatter. We also explore the effects of X-ray cross-calibration and find that Chandra calibration leads to flatter Lx-Tx and M-Tx relations than XMM-Newton.

We present new constraints on the relationship between galaxies and their host dark matter halos, measured from the location of the peak of the stellar-to-halo mass ratio (SHMR), up to the most massive galaxy clusters at redshift $z\sim0.8$ and over a volume of nearly 0.1~Gpc$^3$. We use a unique combination of deep observations in the CFHTLenS/VIPERS field from the near-UV to the near-IR, supplemented by $\sim60\,000$ secure spectroscopic redshifts, analysing galaxy clustering, galaxy-galaxy lensing and the stellar mass function. We interpret our measurements within the halo occupation distribution (HOD) framework, separating the contributions from central and satellite galaxies. We find that the SHMR for the central galaxies peaks at $M_{\rm h, peak} = 1.9^{+0.2}_{-0.1}\times10^{12} M_{\odot}$ with an amplitude of $0.025$, which decreases to $\sim0.001$ for massive halos ($M_{\rm h} > 10^{14} M_{\odot}$). Compared to central galaxies only, the total SHMR (including satellites) is boosted by a factor 10 in the high-mass regime (cluster-size halos), a result consistent with cluster analyses from the literature based on fully independent methods. After properly accounting for differences in modelling, we have compared our results with a large number of results from the literature up to $z=1$: we find good general agreement, independently of the method used, within the typical stellar-mass systematic errors at low to intermediate mass (${M}_{\star} < 10^{11} M_{\odot}$) and the statistical errors above. We have also compared our SHMR results to semi-analytic simulations and found that the SHMR is tilted compared to our measurements in such a way that they over- (under-) predict star formation efficiency in central (satellite) galaxies.

We study the density profile and shape of the Galactic halo using deep multicolour images from the MENeaCS and CCCP projects, over 33 fields selected to avoid overlap with the Galactic plane. Using multicolour selection and PSF homogenization techniques we obtain catalogues of F stars (near-main sequence turnoff stars) out to Galactocentric distances up to 60kpc. Grouping nearby lines of sight, we construct the stellar density profiles through the halo in eight different directions by means of photometric parallaxes. Smooth halo models are then fitted to these profiles. We find clear evidence for a steepening of the density profile power law index around R=20 kpc, from -2.50 +- 0.04 to -4.85 +- 0.04, and for a flattening of the halo towards the poles with best-fit axis ratio 0.63 +- 0.02. Furthermore, we cannot rule out a mild triaxiality (w>=0.8). We recover the signatures of well-known substructure and streams that intersect our lines of sight. These results are consistent with those derived from wider but shallower surveys, and augur well for upcoming, wide-field surveys of comparable depth to our pencil beam surveys.

This paper presents the first application of 3D cosmic shear to a wide-field weak lensing survey. 3D cosmic shear is a technique that analyses weak lensing in three dimensions using a spherical harmonic approach, and does not bin data in the redshift direction. This is applied to CFHTLenS, a 154 square degree imaging survey with a median redshift of 0.7 and an effective number density of 11 galaxies per square arcminute usable for weak lensing. To account for survey masks we apply a 3D pseudo-Cl approach on weak lensing data, and to avoid uncertainties in the highly non-linear regime, we separately analyse radial wave numbers k<=1.5h/Mpc and k<=5.0h/Mpc, and angular wavenumbers l~400-5000. We show how one can recover 2D and tomographic power spectra from the full 3D cosmic shear power spectra and present a measurement of the 2D cosmic shear power spectrum, and measurements of a set of 2-bin and 6-bin cosmic shear tomographic power spectra; in doing so we find that using the 3D power in the calculation of such 2D and tomographic power spectra from data naturally accounts for a minimum scale in the matter power spectrum. We use 3D cosmic shear to constrain cosmologies with parameters OmegaM, OmegaB, sigma8, h, ns, w0, wa. For a non-evolving dark energy equation of state, and assuming a flat cosmology, lensing combined with WMAP7 results in h=0.78+/-0.12, OmegaM=0.252+/-0.079, sigma8=0.88+/-0.23 and w=-1.16+/-0.38 using only scales k<=1.5h/Mpc. We also present results of lensing combined with first year Planck results, where we find no tension with the results from this analysis, but we also find no significant improvement over the Planck results alone. We find evidence of a suppression of power compared to LCDM on small scales 1.5 < k < 5.0 h/Mpc in the lensing data, which is consistent with predictions of the effect of baryonic feedback on the matter power spectrum.

We present the cluster mass-richness scaling relation calibrated by a weak lensing analysis of >18000 galaxy cluster candidates in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). Detected using the 3D-Matched-Filter cluster-finder of Milkeraitis et al., these cluster candidates span a wide range of masses, from the small group scale up to $\sim10^{15} M_{\odot}$, and redshifts 0.2 $\lesssim z\lesssim$ 0.9. The total significance of the stacked shear measurement amounts to 54$\sigma$. We compare cluster masses determined using weak lensing shear and magnification, finding the measurements in individual richness bins to yield 1$\sigma$ compatibility, but with magnification estimates biased low. This first direct mass comparison yields important insights for improving the systematics handling of future lensing magnification work. In addition, we confirm analyses that suggest cluster miscentring has an important effect on the observed 3D-MF halo profiles, and we quantify this by fitting for projected cluster centroid offsets, which are typically $\sim$ 0.4 arcmin. We bin the cluster candidates as a function of redshift, finding similar cluster masses and richness across the full range up to $z \sim$ 0.9. We measure the 3D-MF mass-richness scaling relation $M_{200} = M_0 (N_{200} / 20)^\beta$. We find a normalization $M_0 \sim (2.7^{+0.5}_{-0.4}) \times 10^{13} M_{\odot}$, and a logarithmic slope of $\beta \sim 1.4 \pm 0.1$, both of which are in 1$\sigma$ agreement with results from the magnification analysis. We find no evidence for a redshift-dependence of the normalization. The CFHTLenS 3D-MF cluster catalogue is now available at cfhtlens.org.

Galaxy-galaxy weak lensing is a direct probe of the mean matter distribution around galaxies. The depth and sky coverage of the CFHT Legacy Survey yield statistically significant galaxy halo mass measurements over a much wider range of stellar masses ($10^{8.75}$ to $10^{11.3} M_{\odot}$) and redshifts ($0.2 < z < 0.8$) than previous weak lensing studies. At redshift $z \sim 0.5$, the stellar-to-halo mass ratio (SHMR) reaches a maximum of $4.0\pm0.2$ percent as a function of halo mass at $\sim 10^{12.25} M_{\odot}$. We find, for the first time from weak lensing alone, evidence for significant evolution in the SHMR: the peak ratio falls as a function of cosmic time from $4.5 \pm 0.3$ percent at $z \sim 0.7$ to $3.4 \pm 0.2$ percent at $z \sim 0.3$, and shifts to lower stellar mass haloes. These evolutionary trends are dominated by red galaxies, and are consistent with a model in which the stellar mass above which star formation is quenched "downsizes" with cosmic time. In contrast, the SHMR of blue, star-forming galaxies is well-fit by a power law that does not evolve with time. This suggests that blue galaxies form stars at a rate that is balanced with their dark matter accretion in such a way that they evolve along the SHMR locus. The redshift dependence of the SHMR can be used to constrain the evolution of the galaxy population over cosmic time.

Higher-order, non-Gaussian aspects of the large-scale structure carry valuable information on structure formation and cosmology, which is complementary to second-order statistics. In this work we measure second- and third-order weak-lensing aperture-mass moments from CFHTLenS and combine those with CMB anisotropy probes. The third moment is measured with a significance of $2\sigma$. The combined constraint on $\Sigma_8 = \sigma_8 (\Omega_{\rm m}/0.27)^\alpha$ is improved by 10%, in comparison to the second-order only, and the allowed ranges for $\Omega_{\rm m}$ and $\sigma_8$ are substantially reduced. Including general triangles of the lensing bispectrum yields tighter constraints compared to probing mainly equilateral triangles. Second- and third-order CFHTLenS lensing measurements improve Planck CMB constraints on $\Omega_{\rm m}$ and $\sigma_8$ by 26% for flat $\Lambda$CDM. For a model with free curvature, the joint CFHTLenS-Planck result is $\Omega_{\rm m} = 0.28 \pm 0.02$ (68% confidence), which is an improvement of 43% compared to Planck alone. We test how our results are potentially subject to three astrophysical sources of contamination: source-lens clustering, the intrinsic alignment of galaxy shapes, and baryonic effects. We explore future limitations of the cosmological use of third-order weak lensing, such as the nonlinear model and the Gaussianity of the likelihood function.

Aims. This is the second in a series of papers that attempt to unveil the kinematic structure of the Galactic bulge through studying radial velocities and proper motions. We report here ~15000 new proper motions for three low foreground-extinction off-axis fields of the Galactic bulge. Methods. Proper motions were derived from a combination of Hubble Space Telescope Wide Field Planetary Camera 2 (WFPC2) and Advanced Camera for Surveys (ACS) images taken 8 and 9 years apart, and ACS observations taken 9 and 10 years apart, and they reach accuracies better than 0.9 mas/yr for more than ~10000 objects with magnitudes F814W < 24. Results. The proper motion distributions in these fields are similar to those of Galactic minor axis bulge fields. We observe the rotation of main sequence stars below the turn-off within the Galactic bulge, as in the minor axis fields. Conclusions. Our stellar proper motions measurements show a significant bulge rotation for fields as far from the galactic plane as b=-8.

Studying quasars at the highest redshifts can constrain models of galaxy and black hole formation, and it also probes the intergalactic medium in the early universe. Optical surveys have to date discovered more than 60 quasars up to z~6.4, a limit set by the use of the z-band and CCD detectors. Only one z>6.4 quasar has been discovered, namely the z=7.08 quasar ULAS J1120+0641, using near-infrared imaging. Here we report the discovery of three new z>6.4 quasars in 332 square degrees of the Visible and Infrared Survey Telescope for Astronomy Kilo-degree Infrared Galaxy (VIKING) survey, thus extending the number from 1 to 4. The newly discovered quasars have redshifts of z=6.60, 6.75, and 6.89. The absolute magnitudes are between -26.0 and -25.5, 0.6-1.1 mag fainter than ULAS J1120+0641. Near-infrared spectroscopy revealed the MgII emission line in all three objects. The quasars are powered by black holes with masses of ~(1-2)x10^9 M_sun. In our probed redshift range of 6.44<z<7.44 we can set a lower limit on the space density of supermassive black holes of \rho(M_BH>10^9 M_sun) > 1.1x10^(-9) Mpc^(-3). The discovery of three quasars in our survey area is consistent with the z=6 quasar luminosity function when extrapolated to z~7. We do not find evidence for a steeper decline in the space density of quasars with increasing redshift from z=6 to z=7.

We use data from two CFHT-MegaCam photometric pencil-beam surveys in the g' and the r' bands to measure distances to the Sagittarius, the Palomar 5 and the Orphan stream. We show that, using a cross-correlation algorithm to detect the turnoff point of the main sequence, it is possible to overcome the main limitation of a two-bands pencil-beam survey, namely the lack of adjacent control-fields that can be used to subtract the foreground and background stars to enhance the signal on the colour-magnitude diagrams (CMDs). We describe the cross-correlation algorithm and its implementation. We combine the resulting main sequence turnoff points with theoretical isochrones to derive photometric distances to the streams. Our results (31 detections on the Sagittarius stream and one each for the Palomar 5 and the Orphan streams) confirm the findings by previous studies, expand the distance trend for the Sagittarius faint southern branch and, for the first time, trace the Sagittarius faint branch of the northern-leading arm out to 56 kpc. In addition, they show evidence for new substructure: we argue that these detections trace the continuation of the Sagittarius northern-leading arm into the southern hemisphere, and find a nearby branch of the Sagittarius trailing wrap in the northern hemisphere.

We present a study of the relation between dark matter halo mass and the baryonic content of host galaxies, quantified via luminosity and stellar mass. Our investigation uses 154 deg2 of Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) lensing and photometric data, obtained from the CFHT Legacy Survey. We employ a galaxy-galaxy lensing halo model which allows us to constrain the halo mass and the satellite fraction. Our analysis is limited to lenses at redshifts between 0.2 and 0.4. We express the relationship between halo mass and baryonic observable as a power law. For the luminosity-halo mass relation we find a slope of 1.32+/-0.06 and a normalisation of 1.19+0.06-0.07x10^13 h70^-1 Msun for red galaxies, while for blue galaxies the best-fit slope is 1.09+0.20-0.13 and the normalisation is 0.18+0.04-0.05x10^13 h70^-1 Msun. Similarly, we find a best-fit slope of 1.36+0.06-0.07 and a normalisation of 1.43+0.11-0.08x10^13 h70^-1 Msun for the stellar mass-halo mass relation of red galaxies, while for blue galaxies the corresponding values are 0.98+0.08-0.07 and 0.84+0.20-0.16x10^13 h70^-1 Msun. For red lenses, the fraction which are satellites tends to decrease with luminosity and stellar mass, with the sample being nearly all satellites for a stellar mass of 2x10^9 h70^-2 Msun. The satellite fractions are generally close to zero for blue lenses, irrespective of luminosity or stellar mass. This, together with the shallower relation between halo mass and baryonic tracer, is a direct confirmation from galaxy-galaxy lensing that blue galaxies reside in less clustered environments than red galaxies. We also find that the halo model, while matching the lensing signal around red lenses well, is prone to over-predicting the large-scale signal for faint and less massive blue lenses. This could be a further indication that these galaxies tend to be more isolated than assumed. [abridged]

We present a quantitative analysis of the largest contiguous maps of projected mass density obtained from gravitational lensing shear. We use data from the 154 deg2 covered by the Canada-France-Hawaii Telescope Lensing Survey. Our study is the first attempt to quantitatively characterize the scientific value of lensing maps, which could serve in the future as a complementary approach to the study of the dark universe with gravitational lensing. We show that mass maps contain unique cosmological information beyond that of traditional two-points statistical analysis techniques. Using a series of numerical simulations, we first show how, reproducing the CFHTLenS observing conditions, gravitational lensing inversion provides a reliable estimate of the projected matter distribution of large scale structure. We validate our analysis by quantifying the robustness of the maps with various statistical estimators. We then apply the same process to the CFHTLenS data. We find that the 2-points correlation function of the projected mass is consistent with the cosmological analysis performed on the shear correlation function discussed in the CFHTLenS companion papers. The maps also lead to a significant measurement of the third order moment of the projected mass, which is in agreement with analytic predictions, and to a marginal detection of the fourth order moment. Tests for residual systematics are found to be consistent with zero for the statistical estimators we used. A new approach for the comparison of the reconstructed mass map to that predicted from the galaxy distribution reveals the existence of giant voids in the dark matter maps as large as 3 degrees on the sky. Our analysis shows that lensing mass maps can be used for new techniques such as peak statistics and the morphological analysis of the projected dark matter distribution.

We present the stellar mass functions (SMFs) of star-forming and quiescent galaxies from observations of 10 rich clusters in the Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS) in the redshift range 0.86<z<1.34. We compare our results with field measurements at similar redshifts using data from a Ks-band selected catalogue of the COSMOS/UltraVISTA field. We construct a Ks-band selected multi-colour catalogue for the clusters in 11 photometric bands covering u-8um, and estimate photometric redshifts and stellar masses using SED fitting techniques. To correct for interlopers in our cluster sample, we use the deep spectroscopic component of GCLASS, which contains spectra for 1282 identified cluster and field galaxies taken with Gemini/GMOS. Both the photometric and spectroscopic samples are sufficiently deep that we can probe the SMF down to masses of 10^10 Msun. We distinguish between star-forming and quiescent galaxies using the rest-frame U-V versus V-J diagram, and find that the best-fitting Schechter parameters alpha and M* are similar within the uncertainties for these galaxy types within the different environments. However, there is a significant difference in the shape and normalisation of the total SMF between the clusters and the field sample. This difference in the total SMF is primarily a reflection of the increased fraction of quiescent galaxies in high-density environments. We apply a simple quenching model that includes components of mass- and environment-driven quenching, and find that in this picture 45% of the star-forming galaxies, which normally would be forming stars in the field, are quenched by the cluster. If galaxies in clusters and the field quench their star formation via different mechanisms, these processes have to conspire in such a way that the shapes of the quiescent and star-forming SMF remain similar in these different environments.

The stellar kinematics of the spheroids and discs of S0 galaxies contain clues to their formation histories. Unfortunately, it is difficult to disentangle the two components and to recover their stellar kinematics in the faint outer parts of the galaxies using conventional absorption line spectroscopy. This paper therefore presents the stellar kinematics of six S0 galaxies derived from observations of planetary nebulae (PNe), obtained using the Planetary Nebula Spectrograph. To separate the kinematics of the two components, we use a maximum-likelihood method that combines the discrete kinematic data with a photometric component decomposition. The results of this analysis reveal that: the discs of S0 galaxies are rotationally supported; however, the amount of random motion in these discs is systematically higher than in comparable spiral galaxies; and the S0s lie around one magnitude below the Tully--Fisher relation for spiral galaxies, while their spheroids lie nearly one magnitude above the Faber--Jackson relation for ellipticals. All of these findings are consistent with a scenario in which spirals are converted into S0s through a process of mild harassment or "pestering," with their discs somewhat heated and their spheroid somewhat enhanced by the conversion process. In such a scenario, one might expect the properties of S0s to depend on environment. We do not see such an effect in this fairly small sample, although any differences would be diluted by the fact that the current location does not necessarily reflect the environment in which the transformation occurred. Similar observations of larger samples probing a broader range of environments, coupled with more detailed modelling of the transformation process to match the wide range of parameters that we have shown can now be measured, should take us from these first steps to the definitive answer as to how S0 galaxies form.

The origins of S0 galaxies remain obscure, with various mechanisms proposed for their formation, likely depending on environment. These mechanisms would imprint different signatures in the galaxies' stellar kinematics out to large radii, offering a method for distinguishing between them. We aim to study a sample of six S0 galaxies from a range of environments, and use planetary nebulae (PNe) as tracers of their stellar populations out to very large radii, to determine their kinematics in order to understand their origins. Using a special-purpose instrument, the Planetary Nebula Spectrograph, we observe and extract PNe catalogues for these six systems*. We show that the PNe have the same spatial distribution as the starlight, that the numbers of them are consistent with what would be expected in a comparable old stellar population in elliptical galaxies, and that their kinematics join smoothly onto those derived at smaller radii from conventional spectroscopy. The high-quality kinematic observations presented here form an excellent set for studying the detailed kinematics of S0 galaxies, in order to unravel their formation histories. We find that PNe are good tracers of stellar kinematics in these systems. We show that the recovered kinematics are largely dominated by rotational motion, although with significant random velocities in most cases.

We present a finely-binned tomographic weak lensing analysis of the Canada-France-Hawaii Telescope Lensing Survey, CFHTLenS, mitigating contamination to the signal from the presence of intrinsic galaxy alignments via the simultaneous fit of a cosmological model and an intrinsic alignment model. CFHTLenS spans 154 square degrees in five optical bands, with accurate shear and photometric redshifts for a galaxy sample with a median redshift of zm =0.70. We estimate the 21 sets of cosmic shear correlation functions associated with six redshift bins, each spanning the angular range of 1.5<theta<35 arcmin. We combine this CFHTLenS data with auxiliary cosmological probes: the cosmic microwave background with data from WMAP7, baryon acoustic oscillations with data from BOSS, and a prior on the Hubble constant from the HST distance ladder. This leads to constraints on the normalisation of the matter power spectrum sigma_8 = 0.799 +/- 0.015 and the matter density parameter Omega_m = 0.271 +/- 0.010 for a flat Lambda CDM cosmology. For a flat wCDM cosmology we constrain the dark energy equation of state parameter w = -1.02 +/- 0.09. We also provide constraints for curved Lambda CDM and wCDM cosmologies. We find the intrinsic alignment contamination to be galaxy-type dependent with a significant intrinsic alignment signal found for early-type galaxies, in contrast to the late-type galaxy sample for which the intrinsic alignment signal is found to be consistent with zero.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

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.

We present cosmological constraints from 2D weak gravitational lensing by the large-scale structure in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) which spans 154 square degrees in five optical bands. Using accurate photometric redshifts and measured shapes for 4.2 million galaxies between redshifts of 0.2 and 1.3, we compute the 2D cosmic shear correlation function over angular scales ranging between 0.8 and 350 arcmin. Using non-linear models of the dark-matter power spectrum, we constrain cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. The best constraints from lensing alone are obtained for the small-scale density-fluctuations amplitude sigma_8 scaled with the total matter density Omega_m. For a flat LambdaCDM model we obtain sigma_8(Omega_m/0.27)^0.6 = 0.79+-0.03. We combine the CFHTLenS data with WMAP7, BOSS and an HST distance-ladder prior on the Hubble constant to get joint constraints. For a flat LambdaCDM model, we find Omega_m = 0.283+-0.010 and sigma_8 = 0.813+-0.014. In the case of a curved wCDM universe, we obtain Omega_m = 0.27+-0.03, sigma_8 = 0.83+-0.04, w_0 = -1.10+-0.15 and Omega_K = 0.006+0.006-0.004. We calculate the Bayesian evidence to compare flat and curved LambdaCDM and dark-energy CDM models. From the combination of all four probes, we find models with curvature to be at moderately disfavoured with respect to the flat case. A simple dark-energy model is indistinguishable from LambdaCDM. Our results therefore do not necessitate any deviations from the standard cosmological model.

The OmegaCAM wide-field optical imager is the sole instrument on the VLT Survey Telescope at ESO's Paranal Observatory. The instrument, as well as the telescope, have been designed for surveys with very good, natural seeing-limited image quality over a 1 square degree field. OmegaCAM was commissioned in 2011 and has been observing three ESO Public Surveys in parallel since October 15, 2011. We use the Astro-WISE information system to monitor the calibration of the observatory and to produce the Kilo Degree Survey (KiDS). Here we describe the photometric monitoring procedures in Astro-WISE and give a first impression of OmegaCAM's photometric behavior as a function of time. The long-term monitoring of the observatory goes hand in hand with the KiDS survey production in Astro-WISE. KiDS is observed under partially non-photometric conditions. Based on the first year of OmegaCAM operations it is expected that a $\sim 1%-2%$ photometric homogeneity will be achieved for KiDS.

The Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) comprises deep multi-colour (u*g'r'i'z') photometry spanning 154 square degrees, with accurate photometric redshifts and shape measurements. We demonstrate that the redshift probability distribution function summed over galaxies provides an accurate representation of the galaxy redshift distribution accounting for random and catastrophic errors for galaxies with best fitting photometric redshifts z_p < 1.3. We present cosmological constraints using tomographic weak gravitational lensing by large-scale structure. We use two broad redshift bins 0.5 < z_p <= 0.85 and 0.85 < z_p <= 1.3 free of intrinsic alignment contamination, and measure the shear correlation function on angular scales in the range ~1-40 arcmin. We show that the problematic redshift scaling of the shear signal, found in previous CFHTLS data analyses, does not afflict the CFHTLenS data. For a flat Lambda-CDM model and a fixed matter density Omega_m=0.27, we find the normalisation of the matter power spectrum sigma_8=0.771 \pm 0.041. When combined with cosmic microwave background data (WMAP7), baryon acoustic oscillation data (BOSS), and a prior on the Hubble constant from the HST distance ladder, we find that CFHTLenS improves the precision of the fully marginalised parameter estimates by an average factor of 1.5-2. Combining our results with the above cosmological probes, we find Omega_m=0.2762 \pm 0.0074 and sigma_8=0.802 \pm 0.013.

The spatial variation of the colour of a galaxy may introduce a bias in the measurement of its shape if the PSF profile depends on wavelength. We study how this bias depends on the properties of the PSF and the galaxies themselves. The bias depends on the scales used to estimate the shape, which may be used to optimise methods to reduce the bias. Here we develop a general approach to quantify the bias. Although applicable to any weak lensing survey, we focus on the implications for the ESA Euclid mission. Based on our study of synthetic galaxies we find that the bias is a few times 10^-3 for a typical galaxy observed by Euclid. Consequently, it cannot be neglected and needs to be accounted for. We demonstrate how one can do so using spatially resolved observations of galaxies in two filters. We show that HST observations in the F606W and F814W filters allow us to model and reduce the bias by an order of magnitude, sufficient to meet Euclid's scientific requirements. The precision of the correction is ultimately determined by the number of galaxies for which spatially-resolved observations in at least two filters are available. We use results from the Millennium Simulation to demonstrate that archival HST data will be sufficient for the tomographic cosmic shear analysis with the Euclid dataset.

A likelihood-based method for measuring weak gravitational lensing shear in deep galaxy surveys is described and applied to the Canada-France-Hawaii Telescope (CFHT) Lensing Survey (CFHTLenS). CFHTLenS comprises 154 sq deg of multicolour optical data from the CFHT Legacy Survey, with lensing measurements being made in the i' band to a depth i'(AB)<24.7, for galaxies with signal-to-noise ratio greater than about 10. The method is based on the lensfit algorithm described in earlier papers, but here we describe a full analysis pipeline that takes into account the properties of real surveys. The method creates pixel-based models of the varying point spread function (PSF) in individual image exposures. It fits PSF-convolved two-component (disk plus bulge) models, to measure the ellipticity of each galaxy, with bayesian marginalisation over model nuisance parameters of galaxy position, size, brightness and bulge fraction. The method allows optimal joint measurement of multiple, dithered image exposures, taking into account imaging distortion and the alignment of the multiple measurements. We discuss the effects of noise bias on the likelihood distribution of galaxy ellipticity. Two sets of image simulations that mirror the observed properties of CFHTLenS have been created, to establish the method's accuracy and to derive an empirical correction for the effects of noise bias.

We present data products from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). CFHTLenS is based on the Wide component of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). It encompasses 154 deg^2 of deep, optical, high-quality, sub-arcsecond imaging data in the five optical filters u^*g'r'i'z'. The article presents our data processing of the complete CFHTLenS data set. We were able to obtain a data set with very good image quality and high-quality astrometric and photometric calibration. Our external astrometric accuracy is between 60-70 mas with respect to SDSS data and the internal alignment in all filters is around 30 mas. Our average photometric calibration shows a dispersion on the order of 0.01 to 0.03 mag for g'r'i'z' and about 0.04 mag for u^* with respect to SDSS sources down to i <= 21. In the spirit of the CFHTLS all our data products are released to the astronomical community via the Canadian Astronomy Data Centre. We give a description and how-to manuals of the public products which include image pixel data, source catalogues with photometric redshift estimates and all relevant quantities to perform weak lensing studies.

We present the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) that accurately determines a weak gravitational lensing signal from the full 154 square degrees of deep multi-colour data obtained by the CFHT Legacy Survey. Weak gravitational lensing by large-scale structure is widely recognised as one of the most powerful but technically challenging probes of cosmology. We outline the CFHTLenS analysis pipeline, describing how and why every step of the chain from the raw pixel data to the lensing shear and photometric redshift measurement has been revised and improved compared to previous analyses of a subset of the same data. We present a novel method to identify data which contributes a non-negligible contamination to our sample and quantify the required level of calibration for the survey. Through a series of cosmology-insensitive tests we demonstrate the robustness of the resulting cosmic shear signal, presenting a science-ready shear and photometric redshift catalogue for future exploitation.

Context: The mass-metallicity relationship (MMR) of star-forming galaxies is well-established, however there is still some disagreement with respect to its exact shape and its possible dependence on other observables. Aims: We measure the MMR in the Galaxy And Mass Assembly (GAMA) survey. We compare our measured MMR to that measured in the Sloan Digital Sky Survey (SDSS) and study the dependence of the MMR on various selection criteria to identify potential causes for disparities seen in the literature. Methods: We use strong emission line ratio diagnostics to derive oxygen abundances. We then apply a range of selection criteria for the minimum signal-to-noise in various emission lines, as well as the apparent and absolute magnitude to study variations in the inferred MMR. Results: The shape and position of the MMR can differ significantly depending on the metallicity calibration and selection used. After selecting a robust metallicity calibration amongst those tested, we find that the mass-metallicity relation for redshifts 0.061< z<0.35 in GAMA is in reasonable agreement with that found in the SDSS despite the difference in the luminosity range probed. Conclusions: In view of the significant variations of the MMR brought about by reasonable changes in the sample selection criteria and method, we recommend that care be taken when comparing the MMR from different surveys and studies directly. We also conclude that there could be a modest level of evolution over 0.06<z<0.35 within the GAMA sample.

We use the GAMA I dataset combined with GALEX, SDSS and UKIDSS imaging to construct the low-redshift (z<0.1) galaxy luminosity functions in FUV, NUV, ugriz, and YJHK bands from within a single well constrained volume of 3.4 x 10^5 (Mpc/h)^{3}. The derived luminosity distributions are normalised to the SDSS DR7 main survey to reduce the estimated cosmic variance to the 5 per cent level. The data are used to construct the cosmic spectral energy distribution (CSED) from 0.1 to 2.1 \mum free from any wavelength dependent cosmic variance for both the elliptical and non-elliptical populations. The two populations exhibit dramatically different CSEDs as expected for a predominantly old and young population respectively. Using the Driver et al. (2008) prescription for the azimuthally averaged photon escape fraction, the non-ellipticals are corrected for the impact of dust attenuation and the combined CSED constructed. The final results show that the Universe is currently generating (1.8 +/- 0.3) x 10^{35} h W Mpc^{-3} of which (1.2 +/- 0.1) x 10^{35} h W Mpc^{-3} is directly released into the inter-galactic medium and (0.6 +/- 0.1) x 10^{35} h W Mpc^{-3} is reprocessed and reradiated by dust in the far-IR. Using the GAMA data and our dust model we predict the mid and far-IR emission which agrees remarkably well with available data. We therefore provide a robust description of the pre- and post dust attenuated energy output of the nearby Universe from 0.1micron to 0.6mm. The largest uncertainty in this measurement lies in the mid and far-IR bands stemming from the dust attenuation correction and its currently poorly constrained dependence on environment, stellar mass, and morphology.

Strategic Plan for Astronomy in the Netherlands 2011 - 2020, written by the Netherlands Committee for Astronomy (NCA), on behalf of the excellence research school in astronomy NOVA, (combining the university astronomy institutes of the universities of Amsterdam, Groningen, Leiden and Nijmegen), the NWO division of Physical Sciences, the Netherlands Institute for Radio Astronomy ASTRON and the Netherlands Institute for Space Research SRON. The Strategic plan outlines the scientific priorities for Dutch astronomy in the next decade; the instrumentation effort required to address these priorities, and the connection between astronomical instrumentation and technology development and fundamental technological R&D; the financial contours needed to realise the priorities; and the role of Dutch astronomy in education and outreach. The Strategic Plan also includes a retrospective on the achievements since the last Strategic Plan (2000) and a forward look beyond 2020.

The Kilo Degree Survey (KiDS) is a 1500 square degree optical imaging survey with the recently commissioned OmegaCAM wide-field imager on the VLT Survey Telescope (VST). A suite of data products will be delivered to the European Southern Observatory (ESO) and the community by the KiDS survey team. Spread over Europe, the KiDS team uses Astro-WISE to collaborate efficiently and pool hardware resources. In Astro-WISE the team shares, calibrates and archives all survey data. The data-centric architectural design realizes a dynamic 'live archive' in which new KiDS survey products of improved quality can be shared with the team and eventually the full astronomical community in a flexible and controllable manner.

We measure the two-point angular correlation function of a sample of 4,289,223 galaxies with r < 19.4 mag from the Sloan Digital Sky Survey as a function of photometric redshift, absolute magnitude and colour down to M_r - 5log h = -14 mag. Photometric redshifts are estimated from ugriz model magnitudes and two Petrosian radii using the artificial neural network package ANNz, taking advantage of the Galaxy and Mass Assembly (GAMA) spectroscopic sample as our training set. The photometric redshifts are then used to determine absolute magnitudes and colours. For all our samples, we estimate the underlying redshift and absolute magnitude distributions using Monte-Carlo resampling. These redshift distributions are used in Limber's equation to obtain spatial correlation function parameters from power law fits to the angular correlation function. We confirm an increase in clustering strength for sub-L* red galaxies compared with ~L* red galaxies at small scales in all redshift bins, whereas for the blue population the correlation length is almost independent of luminosity for ~L* galaxies and fainter. A linear relation between relative bias and log luminosity is found to hold down to luminosities L~0.03L*. We find that the redshift dependence of the bias of the L* population can be described by the passive evolution model of Tegmark & Peebles (1998). A visual inspection of a random sample of our r < 19.4 sample of SDSS galaxies reveals that about 10 per cent are spurious, with a higher contamination rate towards very faint absolute magnitudes due to over-deblended nearby galaxies. We correct for this contamination in our clustering analysis.

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

The detailed structure of the Galactic bulge still remain uncertain. The strong difficulties of obtaining observations of stars in the Galactic bulge have hindered the acquisition of a kinematic representation for the inner kpc of the Milky Way. The observation of the 3-d kinematics in several low foreground extinction windows can solve this problem. We have developed a new technique, which combines precise stellar HST positions and proper motions with integral field spectroscopy, in order to obtain reliable 3-d stellar kinematics in crowded fields of the Galactic center. In addition, we present results using the new techniques for six fields in our project. A significant vertex deviation has been found in some of the fields in agreement with previous determinations. This result confirms the presence of a stellar bar in the Galactic bulge.

Here we present the results of various approaches to measure accurate colours and photometric redshifts (photo-z's) from wide-field imaging data. We use data from the Canada-France-Hawaii-Telescope Legacy Survey (CFHTLS) which have been re- processed by the CFHT Lensing Survey (CFHTLenS) team in order to carry out a number of weak gravitational lensing studies. An emphasis is put on the correction of systematic effects in the photo-z's arising from the different Point Spread Functions (PSF) in the five optical bands. Different ways of correcting these effects are discussed and the resulting photo-z accuracies are quantified by comparing the photo-z's to large spectroscopic redshift (spec-z) data sets. Careful homogenisation of the PSF between bands leads to increased overall accuracy of photo-z's. The gain is particularly pronounced at fainter magnitudes where galaxies are smaller and flux measurements are affected more by PSF-effects. We also study possible re- calibrations of the photometric zeropoints (ZPs) with the help of galaxies with known spec-z's. We find that if PSF-effects are properly taken into account, a re-calibration of the ZPs becomes much less important suggesting that previous such re-calibrations described in the literature could in fact be mostly corrections for PSF-effects rather than corrections for real inaccuracies in the ZPs. The implications of this finding for future surveys like KiDS, DES, LSST, or Euclid are mixed. On the one hand, ZP re-calibrations with spec-z's might not be as accurate as previously thought. On the other hand, careful PSF homogenisation might provide a way out and yield accurate, homogeneous photometry without the need for full spectroscopic coverage. This is the first paper in a series describing the technical aspects of CFHTLenS. (abridged)

The Kilo Degree Survey (KiDS) is a 1500 square degree optical imaging survey with the recently commissioned OmegaCAM wide-field imager on the VLT Survey Telescope (VST). A suite of data products will be delivered to ESO and the community by the KiDS survey team. Spread over Europe, the KiDS team uses Astro-WISE to collaborate efficiently and pool hardware resources. In Astro-WISE the team shares, calibrates and archives all survey data. The data-centric architectural design realizes a dynamic 'live archive' in which new KiDS survey products of improved quality can be shared with the team and eventually the full astronomical community in a flexible and controllable manner

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image analysis challenge aimed at the development of algorithms to analyze astronomical images. Specifically, the challenge is to measure varying image distortions in the presence of a variable convolution kernel, pixelization and noise. This is the second in a series of challenges set to the astronomy, computer science and statistics communities, providing a structured environment in which methods can be improved and tested in preparation for planned astronomical surveys. GREAT10 extends upon previous work by introducing variable fields into the challenge. The "Galaxy Challenge" involves the precise measurement of galaxy shape distortions, quantified locally by two parameters called shear, in the presence of a known convolution kernel. Crucially, the convolution kernel and the simulated gravitational lensing shape distortion both now vary as a function of position within the images, as is the case for real data. In addition, we introduce the "Star Challenge" that concerns the reconstruction of a variable convolution kernel, similar to that in a typical astronomical observation. This document details the GREAT10 Challenge for potential participants. Continually updated information is also available from http://www.greatchallenges.info.

Galaxy and Mass Assembly (GAMA) is a project to study galaxy formation and evolution, combining imaging data from ultraviolet to radio with spectroscopic data from the AAOmega spectrograph on the Anglo-Australian Telescope. Using data from phase 1 of GAMA, taken over three observing seasons, and correcting for various minor sources of incompleteness, we calculate galaxy luminosity functions (LFs) and their evolution in the ugriz passbands. At low redshift, z < 0.1, we find that blue galaxies, defined according to a magnitude-dependent but non-evolving colour cut, are reasonably well fit over a range of more than ten magnitudes by simple Schechter functions in all bands. Red galaxies, and the combined blue-plus-red sample, require double power-law Schechter functions to fit a dip in their LF faintward of the characteristic magnitude M* before a steepening faint end. This upturn is at least partly due to dust-reddened disk galaxies. We measure evolution of the galaxy LF over the redshift range 0.002 < z < 0.5 using both a parametric fit and by measuring binned LFs in redshift slices. The characteristic luminosity L* is found to increase with redshift in all bands, with red galaxies showing stronger luminosity evolution than blue galaxies. The comoving number density of blue galaxies increases with redshift, while that of red galaxies decreases, consistent with prevailing movement from blue cloud to red sequence. As well as being more numerous at higher redshift, blue galaxies also dominate the overall luminosity density beyond redshifts z = 0.2. At lower redshifts, the luminosity density is dominated by red galaxies in the riz bands, by blue galaxies in u and g.

We use 14-year baseline images obtained with the Wide Field Planetary Camera 2 on board the Hubble Space telescope to derive a proper motion for one of the Milky Way's most distant dwarf spheroidal companions, Leo II, relative to an extragalactic background reference frame. Astrometric measurements are performed in the effective point spread function (ePSF) formalism using our own developed code. An astrometric reference grid is defined using 3,224 stars that are members of Leo II that are brighter than magnitude 25 in the F814W band. We identify 17 compact extra-galactic sources, for which we measure a systemic proper motion relative to this stellar reference grid. We derive a proper motion [\mu_{\alpha},\mu_{\delta}]=[+104+/-113,-33+/-151] microarcseconds/yr for Leo II in the heliocentric reference frame. Though marginally detected, the proper motion yields constraints on the orbit of Leo II. Given a distance of 230 Kpc and a heliocentric radial velocity +79 km/s, and after subtraction of the solar motion, our measurement indicates a total orbital motion 266.1+/-128.7 km/s in the Galactocentric reference frame, with a radial component +21.5+/-4.3 km/s and tangential component 265.2+/-129.4 km/s. The small radial component indicates that Leo II either has a low-eccentricity orbit, or is currently close to perigalacticon or apogalacticon distance. We see evidence for systematic errors in the astrometry of the extragalactic sources which, while close to being point sources, are slightly resolved in the HST images. We argue that more extensive observations at later epochs will be necessary to better constrain the proper motion of Leo II. We provide a detailed catalog of the stellar and extragalactic sources identified in the HST data which should provide a solid early-epoch reference for future astrometric measurements.

This paper describes the first catalogue of photometrically-derived stellar mass estimates for intermediate-redshift (z < 0.65) galaxies in the Galaxy And Mass Assembly (GAMA) spectroscopic redshift survey. These masses, as well as the full set of ancillary stellar population parameters, will be made public as part of GAMA data release 2. Although the GAMA database does include NIR photometry, we show that the quality of our stellar population synthesis fits is significantly poorer when these NIR data are included. Further, for a large fraction of galaxies, the stellar population parameters inferred from the optical-plus-NIR photometry are formally inconsistent with those inferred from the optical data alone. This may indicate problems in our stellar population library, or NIR data issues, or both; these issues will be addressed for future versions of the catalogue. For now, we have chosen to base our stellar mass estimates on optical photometry only. In light of our decision to ignore the available NIR data, we examine how well stellar mass can be constrained based on optical data alone. We use generic properties of stellar population synthesis models to demonstrate that restframe colour alone is in principle a very good estimator of stellar mass-to-light ratio, M*/Li. Further, we use the observed relation between restframe (g-i) and M*/Li for real GAMA galaxies to argue that, modulo uncertainties in the stellar evolution models themselves, (g-i) colour can in practice be used to estimate M*/Li to an accuracy of < ~0.1 dex. This 'empirically calibrated' (g-i)-M*/Li relation offers a simple and transparent means for estimating galaxies' stellar masses based on minimal data, and so provides a solid basis for other surveys to compare their results to z < ~0.4 measurements from GAMA.

We investigate the properties of satellite galaxies that surround isolated hosts within the redshift range 0.01 < z < 0.15, using data taken as part of the Galaxy And Mass Assembly survey. Making use of isolation and satellite criteria that take into account stellar mass estimates, we find 3514 isolated galaxies of which 1426 host a total of 2998 satellites. Separating the red and blue populations of satellites and hosts, using colour-mass diagrams, we investigate the radial distribution of satellite galaxies and determine how the red fraction of satellites varies as a function of satellite mass, host mass and the projected distance from their host. Comparing the red fraction of satellites to a control sample of small neighbours at greater projected radii, we show that the increase in red fraction is primarily a function of host mass. The satellite red fraction is about 0.2 higher than the control sample for hosts with 11.0 < log M_* < 11.5, while the red fractions show no difference for hosts with 10.0 < log M_* < 10.5. For the satellites of more massive hosts the red fraction also increases as a function of decreasing projected distance. Our results suggest that the likely main mechanism for the quenching of star formation in satellites hosted by isolated galaxies is strangulation.

Using the complete GAMA-I survey covering ~142 sq. deg. to r=19.4, of which ~47 sq. deg. is to r=19.8, we create the GAMA-I galaxy group catalogue (G3Cv1), generated using a friends-of-friends (FoF) based grouping algorithm. Our algorithm has been tested extensively on one family of mock GAMA lightcones, constructed from Lambda-CDM N-body simulations populated with semi-analytic galaxies. Recovered group properties are robust to the effects of interlopers and are median unbiased in the most important respects. G3Cv1 contains 14,388 galaxy groups (with multiplicity >= 2$), including 44,186 galaxies out of a possible 110,192 galaxies, implying ~40% of all galaxies are assigned to a group. The similarities of the mock group catalogues and G3Cv1 are multiple: global characteristics are in general well recovered. However, we do find a noticeable deficit in the number of high multiplicity groups in GAMA compared to the mocks. Additionally, despite exceptionally good local spatial completeness, G3Cv1 contains significantly fewer compact groups with 5 or more members, this effect becoming most evident for high multiplicity systems. These two differences are most likely due to limitations in the physics included of the current GAMA lightcone mock. Further studies using a variety of galaxy formation models are required to confirm their exact origin.

We introduce a novel method to measure the masses of galaxy clusters at high redshift selected from optical and IR Spitzer data via the red-sequence technique. Lyman-break galaxies are used as a well understood, high-redshift background sample allowing mass measurements of lenses at unprecedented high redshifts using weak lensing magnification. By stacking a significant number of clusters at different redshifts with average masses of ~1-3x10^14M_sun, as estimated from their richness, we can calibrate the normalisation of the mass-richness relation. With the current data set (area: 6 deg^2) we detect a magnification signal at the >3-sigma level. There is good agreement between the masses estimated from the richness of the clusters and the average masses estimated from magnification, albeit with large uncertainties. We perform tests that suggest the absence of strong systematic effects and support the robustness of the measurement. This method - when applied to larger data sets in the future - will yield an accurate calibration of the mass-observable relations at z>~1 which will represent an invaluable input for cosmological studies using the galaxy cluster mass function and astrophysical studies of cluster formation. Furthermore this method will probably be the least expensive way to measure masses of large numbers of z>1 clusters detected in future IR-imaging surveys.

The stellar initial mass function (IMF) describes the distribution in stellar masses produced from a burst of star formation. For more than fifty years, the implicit assumption underpinning most areas of research involving the IMF has been that it is universal, regardless of time and environment. We measure the high-mass IMF slope for a sample of low-to-moderate redshift galaxies from the Galaxy And Mass Assembly survey. The large range in luminosities and galaxy masses of the sample permits the exploration of underlying IMF dependencies. A strong IMF-star formation rate dependency is discovered, which shows that highly star forming galaxies form proportionally more massive stars (they have IMFs with flatter power-law slopes) than galaxies with low star formation rates. This has a significant impact on a wide variety of galaxy evolution studies, all of which rely on assumptions about the slope of the IMF. Our result is supported by, and provides an explanation for, the results of numerous recent explorations suggesting a variation of or evolution in the IMF.

Context: In 2.2m MPG-ESO/WFI data we discovered several mass peaks through weak lensing, forming a possible supercluster at redshift 0.45. Through multi-colour wide-field imaging with CFHT/Megaprime and INT/WFC we identify early-type galaxies and trace the supercluster network with them. Through EMMI/NTT multi-object spectroscopy we verify the initial shear-selected cluster candidates. Using weak lensing we obtain mass estimates for the supercluster centre and the filaments. Results: We identified the centre of the SCL2243-0935 supercluster, MACS J2243-0935, which was found independently by Ebeling et al. (2010). 13 more clusters or overdensities are embedded in a filamentary network, half of them are already spectroscopically confirmed. Three (5-15) Mpc filaments are detected, and we estimate the global size of SCL2243 to 45x15x50 Mpc, making it one of the largest superclusters known at intermediate redshifts. Weak lensing yields r_200=(2.06+/-0.13) Mpc and M_200=(1.54+/-0.29)x10^15 M_sun for MACS J2243 with M/L=428+/-82, very similar to results from size-richness cluster scaling relations. Integrating the weak lensing surface mass density over the supercluster network (defined by increased i-band luminosity or g-i colours), we find (1.53+/-1.01)x10^15 M_sun and M/L=305+/-201 for the three main filaments, consistant with theoretical predictions. The filaments' projected surface mass density is 0.007-0.012, corresponding to 10-100 times the critical density. The greatly varying density of the cosmic web is also reflected in the mean colour of galaxies. Conclusions: SCL2243 is significantly larger and much more richly structured than other known superclusters such as A901/902 or MS0302 studied with weak lensing before. It is a text-book supercluster with little contamination along the line of sight, making it a perfect sandbox for testing new techniques probing the cosmic web.

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

Boxy and peanut-shaped bulges are seen in about half of edge-on disc galaxies. Comparisons of the photometry and major-axis gas and stellar kinematics of these bulges to simulations of bar formation and evolution indicate that they are bars viewed in projection. If the properties of boxy bulges can be entirely explained by assuming they are bars, then this may imply that their hosts are pure disc galaxies with no classical bulge. A handful of these bulges, including that of the Milky Way, have been observed to rotate cylindrically, i.e. with a mean stellar velocity independent of height above the disc. In order to assess whether such behaviour is ubiquitous in boxy bulges, and whether a pure disc interpretation is consistent with their stellar populations, we have analysed the stellar kinematics and populations of the boxy or peanut-shaped bulges in a sample of five edge-on galaxies. We placed slits along the major axis of each galaxy and at three offset but parallel positions to build up spatial coverage. The boxy bulge of NGC3390 rotates perfectly cylindrically within the spatial extent and uncertainties of the data. This is consistent with the metallicity and alpha-element enhancement of the bulge, which are the same as in the disk. This galaxy is thus a pure disc galaxy. The boxy bulge of ESO311-G012 also rotates very close to cylindrically. The boxy bulge of NGC1381 is neither clearly cylindrically nor non-cylindrically rotating, but it has a negative vertical metallicity gradient and is alpha-enhanced with respect to its disc, suggesting a composite bulge comprised of a classical bulge and bar (and possibly a discy pseudobulge) [abridged] Even this relatively small sample is sufficient to demonstrate that boxy bulges display a range of rotational and population properties, indicating that they do not form a homogeneous class of object.

To investigate the origins of S0 galaxies, we present a new method of analyzing their stellar kinematics from discrete tracers such as planetary nebulae. This method involves binning the data in the radial direction so as to extract the most general possible non-parametric kinematic profiles, and using a maximum likelihood fit within each bin in order to make full use of the information in the discrete kinematic tracers. Both disk and spheroid kinematic components are fitted, with a two-dimensional decomposition of imaging data used to attribute to each tracer a probability of membership in the separate components. Likelihood clipping also allows us to identify objects whose properties are not consistent with the adopted model, rendering the technique robust against contaminants and able to identify additional kinematic features. The method is first tested on an N-body simulated galaxy to assess possible sources of systematic error associated with the structural and kinematic decomposition, which are found to be small. It is then applied to the S0 system NGC~1023, for which a planetary nebula catalogue has already been released and analyzed by (Noordermeer et al., 2008). The correct inclusion of the spheroidal component allows us to show that, contrary to previous claims, the stellar kinematics of this galaxy are indistinguishable from those of a normal spiral galaxy, indicating that it may have evolved directly from such a system via gas stripping or secular evolution. The method also successfully identifies a population of outliers whose kinematics are different from those of the main galaxy; these objects can be identified with a stellar stream associated with the companion galaxy NGC~1023A.

We present an analysis of the properties of the lowest Halpha-luminosity galaxies (L_Halpha<4x10^32 W; SFR<0.02 Msun/yr) in the Galaxy And Mass Assembly (GAMA) survey. These galaxies make up the the rise above a Schechter function in the number density of systems seen at the faint end of the Halpha luminosity function. Above our flux limit we find that these galaxies are principally composed of intrinsically low stellar mass systems (median stellar mass =2.5x10^8 Msun) with only 5/90 having stellar masses M>10^10 Msun. The low SFR systems are found to exist predominantly in the lowest density environments (median density ~0.02 galaxy Mpc^-2 with none in environments more dense than ~1.5 galaxy Mpc^-2). Their current specific star formation rates (SSFR; -8.5 < log(SSFR[yr^-1])<-12.) are consistent with their having had a variety of star formation histories. The low density environments of these galaxies demonstrates that such low-mass, star-forming systems can only remain as low-mass and forming stars if they reside sufficiently far from other galaxies to avoid being accreted, dispersed through tidal effects or having their gas reservoirs rendered ineffective through external processes.

Current theories of structure formation predict specific density profiles of galaxy dark matter haloes, and with weak gravitational lensing we can probe these profiles on several scales. On small scales, higher-order shape distortions known as flexion add significant detail to the weak lensing measurements. We present here the first detection of a galaxy-galaxy flexion signal in space-based data, obtained using a new Shapelets pipeline introduced here. We combine this higher-order lensing signal with shear to constrain the average density profile of the galaxy lenses in the Hubble Space Telescope COSMOS survey. We also show that light from nearby bright objects can significantly affect flexion measurements. After correcting for the influence of lens light, we show that the inclusion of flexion provides tighter constraints on density profiles than does shear alone. Finally we find an average density profile consistent with an isothermal sphere.

As part of our current programme to test LCDM predictions for dark matter (DM) haloes using extended kinematical observations of early-type galaxies, we present a dynamical analysis of the bright elliptical galaxy NGC 4374 (M84) based on ~450 Planetary Nebulae (PNe) velocities from the PN.Spectrograph, along with extended long-slit stellar kinematics. This is the first such analysis of a galaxy from our survey with a radially constant velocity dispersion profile. We find that the spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The velocity kurtosis is consistent with zero at almost all radii. We construct a series of Jeans models, fitting both velocity dispersion and kurtosis to help break the mass-anisotropy degeneracy. Our mass models include DM halos either with shallow cores or with central cusps as predicted by cosmological simulations - along with the novel introduction in this context of adiabatic halo contraction from baryon infall. Both classes of models confirm a very massive dark halo around NGC 4374, demonstrating that PN kinematics data are well able to detect such haloes when present. Considering the default cosmological mass model, we confirm earlier suggestions that bright galaxies tend to have halo concentrations higher than LCDM predictions, but this is found to be solved if either a Salpeter IMF or adiabatic contraction with a Kroupa IMF is assumed. Thus for the first time a case is found where the PN dynamics may well be consistent with a standard dark matter halo. A cored halo can also fit the data, and prefers a stellar mass consistent with a Salpeter IMF. The less dramatic dark matter content found in lower-luminosity "ordinary" ellipticals suggests a bimodality in the halo properties which may be produced by divergent baryonic effects during their assembly histories.

Dark matter is an important ingredient of galaxies, as was recognised early on by Ken Freeman himself! Evidence for dark matter halos is still indirect, based on analysing motions of tracers such as gas and stars. In a sense the visible galaxy is the mask through which we can study the dark matter. Light rays are also sensitive to gravitational fields, and dark haloes cause observable gravitational lensing effects. There are three regimes: microlensing (which probes the clumpiness of dark matter haloes), strong lensing (sensitive to the inner mass distribution) and weak lensing (which can probe haloes out to 100s of kpc from the center). This review will concentrate on weak lensing, and describe a new survey, the Kilo-Degree Survey (KiDS) that is designed to study galaxy halo masses, extents and shapes as a function of environment, galaxy type and redshift.

We present self-consistent star formation rates derived through pan-spectral analysis of galaxies drawn from the Galaxy and Mass Assembly (GAMA) survey. We determine the most appropriate form of dust obscuration correction via application of a range of extinction laws drawn from the literature as applied to Halpha, [O{II}] and UV luminosities. These corrections are applied to a sample of 31,508 galaxies from the GAMA survey at z < 0.35. We consider several different obscuration curves, including those of Milky Way, Calzetti (2001) and Fischera and Dopita (2005) curves and their effects on the observed luminosities. At the core of this technique is the observed Balmer decrement, and we provide a prescription to apply optimal obscuration corrections using the Balmer decrement. We carry out an analysis of the star formation history (SFH) using stellar population synthesis tools to investigate the evolutionary history of our sample of galaxies as well as to understand the effects of variation in the Initial Mass Function (IMF) and the effects this has on the evolutionary history of galaxies. We find that the Fischera and Dopita (2005) obscuration curve with an R_{v} value of 4.5 gives the best agreement between the different SFR indicators. The 2200A feature needed to be removed from this curve to obtain complete consistency between all SFR indicators suggesting that this feature may not be common in the average integrated attenuation of galaxy emission. We also find that the UV dust obscuration is strongly dependent on the SFR.

In order to generate credible 0.1-2 {\mu}m SEDs, the GAMA project requires many Gigabytes of imaging data from a number of instruments to be re-processed into a standard format. In this paper we discuss the software infrastructure we use, and create self-consistent ugrizYJHK photometry for all sources within the GAMA sample. Using UKIDSS and SDSS archive data, we outline the pre-processing necessary to standardise all images to a common zeropoint, the steps taken to correct for seeing bias across the dataset, and the creation of Gigapixel-scale mosaics of the three 4x12 deg GAMA regions in each filter. From these mosaics, we extract source catalogues for the GAMA regions using elliptical Kron and Petrosian matched apertures. We also calculate S\'ersic magnitudes for all galaxies within the GAMA sample using SIGMA, a galaxy component modelling wrapper for GALFIT 3. We compare the resultant photometry directly, and also calculate the r band galaxy LF for all photometric datasets to highlight the uncertainty introduced by the photometric method. We find that (1) Changing the object detection threshold has a minor effect on the best-fitting Schechter parameters of the overall population (M* +/- 0.055mag, {\alpha} +/- 0.014, {\Phi}* +/- 0.0005 h^3 Mpc^{-3}). (2) An offset between datasets that use Kron or Petrosian photometry regardless of the filter. (3) The decision to use circular or elliptical apertures causes an offset in M* of 0.20mag. (4) The best-fitting Schechter parameters from total-magnitude photometric systems (such as SDSS modelmag or S\'ersic magnitudes) have a steeper faint-end slope than photometry dependent on Kron or Petrosian magnitudes. (5) Our Universe's total luminosity density, when calculated using Kron or Petrosian r-band photometry, is underestimated by at least 15%.

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)

Aims. The Large Binocular Cameras (LBC) are two twin wide field cameras (FOV ~ 23'x 25') mounted at the prime foci of the 8.4m Large Binocular Telescope (LBT). We performed a weak lensing analysis of the z=0.288 cluster Abell 611 on g-band data obtained by the blue-optimized Large Binocular Camera in order to estimate the cluster mass. Methods. Due to the complexity of the PSF of LBC, we decided to use two different approaches, KSB and Shapelets, to measure the shape of background galaxies and to derive the shear signal produced by the cluster. Then we estimated the cluster mass with both aperture densitometry and parametric model fits. Results. The combination of the large aperture of the telescope and the wide field of view allowed us to map a region well beyond the expected virial radius of the cluster and to get a high surface density of background galaxies (23 galaxies/arcmin^2). This made possible to estimate an accurate mass for Abell 611. We find that the mass within 1.5 Mpc is: $(8\pm3)\times 10^{14} M_\odot$ from the aperture mass technique and $(5\pm1)\times 10^{14} M_\odot$ using the model fitting by a NFW mass density profile, for both Shapelets and KSB methods. This analysis demonstrates that LBC is a powerful instrument for weak gravitational lensing studies.

A heuristic greedy algorithm is developed for efficiently tiling spatially dense redshift surveys. In its first application to the Galaxy and Mass Assembly (GAMA) redshift survey we find it rapidly improves the spatial uniformity of our data, and naturally corrects for any spatial bias introduced by the 2dF multi object spectrograph. We make conservative predictions for the final state of the GAMA redshift survey after our final allocation of time, and can be confident that even if worse than typical weather affects our observations, all of our main survey requirements will be met.

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

We present new results from BRAVA, a large scale radial velocity survey of the Galactic bulge, using M giant stars selected from the Two Micron All Sky Survey catalogue as targets for the Cerro Tololo Inter-American Observatory 4m Hydra multi-object spectrograph. The purpose of this survey is to construct a new generation of self-consistent bar models that conform to these observations. We report the dynamics for fields at the edge of the Galactic bulge at latitudes b=-8 deg. and compare to the dynamics at b=-4 deg. We find that the rotation curve V(r) is the same at b=-8 deg. as at b=-4 deg. That is, the Galactic boxy bulge rotates cylindrically, as do boxy bulges of other galaxies. The summed line of sight velocity distribution at b=-8 deg. is Gaussian, and the binned longitude-velocity plot shows no evidence for either a (disk) population with cold dynamics or for a (classical bulge) population with hot dynamics. The observed kinematics are well modeled by an edge-on N-body bar, in agreement with published structural evidence. Our kinematic observations indicate that the Galactic bulge is a prototypical product of secular evolution in galaxy disks, in contrast with stellar population results that are most easily understood if major mergers were the dominant formation process.

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses on a problem that is of crucial importance for future observations in cosmology. The shapes of distant galaxies can be used to determine the properties of dark energy and the nature of gravity, because light from those galaxies is bent by gravity from the intervening dark matter. The observed galaxy images appear distorted, although only slightly, and their shapes must be precisely disentangled from the effects of pixelisation, convolution and noise. The worldwide gravitational lensing community has made significant progress in techniques to measure these distortions via the Shear TEsting Program (STEP). Via STEP, we have run challenges within our own community, and come to recognise that this particular image analysis problem is ideally matched to experts in statistical inference, inverse problems and computational learning. Thus, in order to continue the progress seen in recent years, we are seeking an infusion of new ideas from these communities. This document details the GREAT08 Challenge for potential participants. Please visit http://www.great08challenge.info for the latest information.

We present new planetary nebulae (PNe) positions, radial velocities, and magnitudes for 6 early-type galaxies obtained with the Planetary Nebulae Spectrograph, their two-dimensional velocity and velocity dispersion fields. We extend this study to include an additional 10 early-type galaxies with PNe radial velocity measurements available from the literature, to obtain a broader description of the outer-halo kinematics in early-type galaxies. These data extend the information derived from stellar kinematics to typically up to ~8 Re. The combination of photometry, stellar and PNe kinematics shows: i) good agreement between the PNe number density and the stellar surface brightness in the region where the two data sets overlap; ii) good agreement between PNe and stellar kinematics; iii) that the mean rms velocity profiles fall into two groups: with of the galaxies characterized by slowly decreasing profiles and the remainder having steeply falling profiles; iv) a larger variety of velocity dispersion profiles; v) that twists and misalignments in the velocity fields are more frequent at large radii, including some fast rotators; vi) that outer haloes are characterised by more complex radial profiles of the specific angular momentum-related lambda_R parameter than observed within 1Re; vii) that many objects are more rotationally dominated at large radii than in their central parts; and viii) that the halo kinematics are correlated with other galaxy properties, such as total luminosity, isophotal shape, total stellar mass, V/sigma, and alpha parameter, with a clear separation between fast and slow rotators.

We present first results of a study of the halo kinematics for a sample of early type galaxies using planetary nebulae (PNe) as kinematical tracers. PNe allow to extend up to several effective radii (Re) the information from absorption line kinematics (confined to within 1 or 2 Re), providing valuable information and constraints for merger simulations and galaxy formation models. We find that the specific angular momentum per unit mass has a more complex radial dependence when the halo region is taken into account and that the halo velocity dispersion is related to the total galaxy luminosity, isophotal shape, and number of PNe per unit of luminosity

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 PNe out to 7 effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside 1 Re. The velocity dispersion profile declines with radius, though not very steeply, down to ~70 km/s at the last data point. We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component LCDM-motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fit solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio, and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model. Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration halos, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.

Recent results from the Planetary Nebula Spectrograph (PN.S) survey have revealed a rapidly falling velocity dispersion profile in the nearby elliptical galaxy NGC 3379, casting doubts on whether this intermediate-luminosity galaxy has the kind of dark matter halo expected in LambdaCDM cosmology. We present a detailed dynamical study of this galaxy, combining long-slit spectroscopy, SAURON integral-field data, and PN.S velocities, reaching to more than seven effective radii (R_e). We construct spherical and axisymmetric dynamical models for these data with the flexible made-to-measure NMAGIC code, in a sequence of gravitational potentials with varying dark halo mass. We find that the data are consistent both with near-isotropic spherical systems dominated by the stellar mass, and with models in massive halos with strongly radially anisotropic outer parts (beta >~ 0.8 at 7R_e). Formal likelihood limits would exclude (at 1 sigma) the model with stars only, as well as halo models with v_circ(7R_e) >~ 250 km/s. A sequence of more realistic axisymmetric models of different inclinations and a small number of triaxial tests confirm the spherical results. All valid models fitting all the data are dynamically stable over Gyrs, including the most anisotropic ones. NGC 3379 may well have a dark matter halo as predicted by recent merger models within LambdaCDM cosmology, provided its outer envelope is strongly radially anisotropic. (abridged)

Context: Difference imaging has proven to be a powerful technique for detecting and monitoring the variability of unresolved stellar sources in M 31. Using this technique in surveys of galaxies outside the Local Group could have many interesting applications. Aims: The goal of this paper is to test difference imaging photometry on Centaurus A, the nearest giant elliptical galaxy, at a distance of 4 Mpc. Methods: We obtained deep photometric data with the Wide Field Imager at the ESO/MPG 2.2m at La Silla spread over almost two months. Applying the difference imaging photometry package DIFIMPHOT, we produced high-quality difference images and detected variable sources. The sensitivity of the current observational setup was determined through artificial residual tests. Results: In the resulting high-quality difference images, we detect 271 variable stars. We find a difference flux detection limit corresponding to m_R~24.5. Based on a simple model of the halo of Centaurus A, we estimate that a ground-based microlensing survey would detect in the order of 4 microlensing events per year due to lenses in the halo. Conclusions: Difference imaging photometry works very well at the distance of Centaurus A and promises to be a useful tool for detecting and studying variable stars in galaxies outside the local group. For microlensing surveys, a higher sensitivity is needed than achieved here, which would be possible with a large ground-based telescope or space observatory with wide-field imaging capabilities.

We investigate the manner in which lenticular galaxies are formed by studying their stellar kinematics: an S0 formed from a fading spiral galaxy should display similar cold outer disc kinematics to its progenitor, while an S0 formed in a minor merger should be more dominated by random motions. In a pilot study to attempt to distinguish between these scenarios, we have measured the planetary nebula (PN) kinematics of the nearby S0 system NGC 1023. Using the Planetary Nebula Spectrograph, we have detected and measured the line-of-sight velocities of 204 candidate PNe in the field of this galaxy. Out to intermediate radii, the system displays the kinematics of a normal rotationally-supported disc system. After correction of its rotational velocities for asymmetric drift, the galaxy lies just below the spiral galaxy Tully-Fisher relation, as one would expect for a fading system. However, at larger radii the kinematics undergo a gradual but major transition to random motion with little rotation. This transition does not seem to reflect a change in the viewing geometry or the presence of a distinct halo component, since the number counts of PNe follow the same simple exponential decline as the stellar continuum with the same projected disc ellipticity out to large radii. The galaxy's small companion, NGC 1023A, does not seem to be large enough to have caused the observed modification either. This combination of properties would seem to indicate a complex evolutionary history in either the transition to form an S0 or in the past life of the spiral galaxy from which the S0 formed. More data sets of this type from both spirals and S0s are needed in order to definitively determine the relationship between these types of system.

We consider the signature and detectability of gravitational microlensing of distant quasars by cosmic strings. Because of the simple image configuration such events will have a characteristic light curve, in which a source would appear to brighten by exactly a factor of two, before reverting to its original apparent brightness. We calculate the optical depth and event rate, and conclude that current predictions and limits on the total length of strings on the sky imply optical depths of $\la 10^{-8}$ and event rates of fewer than one event per $10^9$ sources per year. Disregarding those predictions but replacing them with limits on the density of cosmic strings from the CMB fluctuation spectrum, leaves only a small region of parameter space (in which the sky contains about $3\times10^5$ strings with deficit angle of order 0.3 milli-arcseconds) for which a microlensing survey of exposure $10^7$ source-years, spanning a 20--40-year period, might reveal the presence of cosmic strings.

We examine the dark matter properties of nearby early-type galaxies using planetary nebulae (PNe) as mass probes. We have designed a specialised instrument, the Planetary Nebula Spectrograph (PN.S) operating at the William Herschel telescope, with the purpose of measuring PN velocities with best efficiency. The primary scientific objective of this custom-built instrument is the study of the PN kinematics in 12 ordinary round galaxies. Preliminary results showing a dearth of dark matter in ordinary galaxies (Romanowsky et al. 2003) are now confirmed by the first complete PN.S datasets. On the other hand early-type galaxies with a "regular" dark matter content are starting to be observed among the brighter PN.S target sample, thus confirming a correlation between the global dark-to-luminous mass virial ratio (f_DM=M_DM/M_star) and the galaxy luminosity and mass.

We present results from Planetary Nebula Spectrograph (PN.S) observations of the elliptical galaxy NGC 3379 and a description of the data reduction pipeline. We detected 214 planetary nebulae of which 191 are ascribed to NGC 3379, and 23 to the companion galaxy NGC 3384. Comparison with data from the literature show that the PN.S velocities have an internal error of <20km/s and a possible offset of similar magnitude. We present the results of kinematic modelling and show that the PN kinematics are consistent with absorption-line data in the region where they overlap. The resulting combined kinematic data set, running from the center of NGC 3379 out to more than seven effective radii (Reff), reveals a mean rotation velocity that is small compared to the random velocities, and a dispersion profile that declines rapidly with radius. From a series of Jeans dynamical models we find the B-band mass-to-light ratio inside 5 Reff to be 8 to 12 in solar units, and the dark matter fraction inside this radius to be less than 40%. We compare these and other results of dynamical analysis with those of dark-matter-dominated merger simulations, finding that significant discrepancies remain, reiterating the question of whether NGC 3379 has the kind of dark matter halo that the current LambdaCDM paradigm requires.

The recent explosion of recorded digital data and its processed derivatives threatens to overwhelm researchers when analysing their experimental data or when looking up data items in archives and file systems. While current hardware developments allow to acquire, process and store 100s of terabytes of data at the cost of a modern sports car, the software systems to handle these data are lagging behind. This general problem is recognized and addressed by various scientific communities, e.g., DATAGRID/EGEE federates compute and storage power over the high-energy physical community, while the astronomical community is building an Internet geared Virtual Observatory, connecting archival data. These large projects either focus on a specific distribution aspect or aim to connect many sub-communities and have a relatively long trajectory for setting standards and a common layer. Here, we report "first light" of a very different solution to the problem initiated by a smaller astronomical IT community. It provides the abstract "scientific information layer" which integrates distributed scientific analysis with distributed processing and federated archiving and publishing. By designing new abstractions and mixing in old ones, a Science Information System with fully scalable cornerstones has been achieved, transforming data systems into knowledge systems. This break-through is facilitated by the full end-to-end linking of all dependent data items, which allows full backward chaining from the observer/researcher to the experiment. Key is the notion that information is intrinsic in nature and thus is the data acquired by a scientific experiment. The new abstraction is that software systems guide the user to that intrinsic information by forcing full backward and forward chaining in the data modelling.

The OmegaCAM wide-field imager will start operations at the ESO VLT Survey Telescope at Paranal, Chile, in 2007. The photometric calibration of OmegaCAM data depends on standard star measurements that cover the complete 1 square degree FoV of OmegeaCAM. A catalog fullfilling this requirement for eight Landolt equatorial fields, denoted the OmegaCAM Secondary Standards Catalog, will be constructed from OmegaCAM observations during the first year of operations. Here we present the 'Preliminary Catalog' which will be used to bootstrap the construction of the OmegaCAM Secondary Standards Catalog. Thus the Preliminary Catalog will be used to assess the performance of OmegaCAM+VST early-on. The catalog is based on WFC data from the INT at La Palma.

We combine radial velocities, proper motions, and low resolution abundances for a sample of 315 K and M giants in the Baade's Window (l,b)=(0.9,-4) Galactic bulge field. The velocity ellipsoid of stars with [Fe/H]>-0.5 dex shows a vertex deviation in the plot of radial versus transverse velocity, consistent with that expected from a population with orbits supporting a bar. We demonstrate that the significance of this vertex deviation using non-parametric rank correlation statistic is >99%. The velocity ellipsoid for the metal poor ([FeH]<-0.5) part of the population shows no vertex deviation and is consistent with an isotropic, oblate rotating population. We find no evidence for kinematic subgroups, but there is a mild tendency for the vertical velocity dispersion sigma_b to decrease with increasing metallicity.

We describe a new technique for measuring accurate galaxy colours from images taken under different seeing conditions. The method involves two ingredients. First we define the Gaussian-aperture-and-PSF flux, which is the Gaussian-weighted flux a galaxy would have if it were observed with a Gaussian PSF. This theoretical aperture flux is independent of the PSF or pixel scale that the galaxy was observed with. Second we develop a procedure to measure such a `GaaP' flux from observed, pixellated images. This involves modelling source and PSF as a superposition of orthogonal shapelets. A correction scheme is also described, which approximately corrects for any residuals to the shapelet expansions. A series of tests on simulated images shows that with this method it is possible to reduce systematic errors in the matched-aperture fluxes to a percent, which makes it useful for deriving photometric redshifts from large imaging surveys.

The Shear TEsting Programme (STEP) is a collaborative project to improve the accuracy and reliability of weak lensing measurement, in preparation for the next generation of wide-field surveys. We review sixteen current and emerging shear measurement methods in a common language, and assess their performance by running them (blindly) on simulated images that contain a known shear signal. We determine the common features of algorithms that most successfully recover the input parameters. We achieve previously unattained discriminatory precision in our analysis, via a combination of more extensive simulations, and pairs of galaxy images that have been rotated with respect to each other, thus removing noise from their intrinsic ellipticities. The robustness of our simulation approach is also confirmed by testing the relative calibration of methods on real data. Weak lensing measurement has improved since the first STEP paper. Several methods now consistently achieve better than 2% precision, and are still being developed. However, the simulations can now distinguish all methods from perfect performance. Our main concern continues to be the potential for a multiplicative shear calibration bias: not least because this can not be internally calibrated with real data. We determine which galaxy populations are responsible and, by adjusting the simulated observing conditions, we also investigate the effects of instrumental and atmospheric parameters. We have isolated several previously unrecognised aspects of galaxy shape measurement, in which focussed development could provide further progress towards the sub-percent level of precision desired for future surveys. [ABRIDGED]

We derive the distribution of maser stars in the inner Milky Way (MW) based on an analysis of lV-diagrams (lVd) for two samples of maser stars: 771 OH/IR stars and 363 SiO-maser stars. They are all close to the plane of the MW and have long. from -45 to +45deg. The two lVds are compared and found to be very similar. They also compare well with the lVd of interstellar CO, but there are significant differences in detail between the stellar lVds and that of the ISM. Based on the qualitative discussion we divide the lVds into seven areas. In each area we compare the number of stars observed with those predicted by an assumed set of orbits in a galactic potential. This potential is axially symmetric but a weak rotating bar has been added. We conclude that the maser stars move on almost circular orbits outside of about 3.5 kpc, but that the orbits become more and more elongated when one goes deep inside our MW. We find a strong effect of the Corotation (CR) resonance (res) at 3.3 kpc, we see a small but noticeable effect of the Outer Lindblad res at 5 kpc and no effect of the Inner Lindblad res at r=0.8 kpc. We find a set of 6 groups of orbits that together predict counts in agreement with the counts of stars observed. We then calculate the trajectory of each orbit and so find the distribution of the maser stars in the plane of the MWG. This distribution has two new (but not unexpected) features. The first is a bar-like distribution within 2 kpc from the GC outlined. These orbits explain the high-vel stars near l=0deg in the forbidden and the permitted quadrants. The second feature are two "croissant"-like voids in the distribution close to the CR radius (3.3 kpc), which are the consequence of the presence of the CR res. We find excellent agreement with an earlier reconstruction by Sevenster (1999).

Accurate measurement of gravitational shear from images of distant galaxies is one of the most direct ways of studying the distribution of mass in the universe. We describe an implementation of a technique that is based on the shapelets formalism. The shapelets technique describes PSF and observed images in terms of Gauss-Hermite expansions (Gaussians times polynomials). It allows the various operations that a galaxy image undergoes before being registered in a camera (gravitational shear, PSF convolution, pixelation) to be modeled in a single formalism, so that intrinsic ellipticities can be derived in a single modeling step. The resulting algorithm, and tests of it on idealized data as well as more realistic simulated images from the STEP project, are described. Results are very promising, with attained calibration accuracy better than four percent (1 percent for round PSFs) and PSF ellipticity correction better than a factor of 20. Residual calibration problems are discussed.

We present radial velocities for a sample of 723 planetary nebulae (PNe) in the disk and bulge of M31, measured using the WYFFOS fibre spectrograph on the William Herschel telescope. Velocities are determined using the [OIII] 5007 Angstrom emission line. Rotation and velocity dispersion are measured to a radius of 50 arcminutes (11.5 kpc), the first stellar rotation curve and velocity dispersion profile for M31 to such a radius. Our kinematics are consistent with rotational support at radii well beyond the bulge effective radius of 1.4kpc, although our data beyond a radius of 5kpc are limited. We present tentative evidence for kinematic substructure in the bulge of M31 to be studied fully in a later work. This paper is part of an ongoing project to constrain the total mass, mass distribution and velocity anisotropy of the disk, bulge and halo of M31.

We present a catalogue of positions, magnitudes and velocities for 3300 emission-line objects found by the Planetary Nebula Spectrograph in a survey of the Andromeda Galaxy, M31. Of these objects, 2615 are found likely to be planetary nebulae (PNe) associated with M31. Initial results from this survey include: the likely non-existence of Andromeda VIII; a universal PN luminosity function, with the exception of a small amount of obscuration, and a small offset in normalization between bulge and disk components; very faint kinematically-selected photometry implying no cut-off in the disk to beyond 4 scalelengths and no halo population in excess of the bulge out to 10 effective bulge radii; disk kinematics that show significant dispersion and asymmetric drift out to large radii, consistent with a warm flaring disk; and no sign of any variation in kinematics with PN luminosity, suggesting that PNe arise from a fairly uniform population of old stars.

The Shear TEsting Programme, STEP, is a collaborative project to improve the accuracy and reliability of all weak lensing measurements in preparation for the next generation of wide-field surveys. In this first STEP paper we present the results of a blind analysis of simulated ground-based observations of relatively simple galaxy morphologies. The most successful methods are shown to achieve percent level accuracy. From the cosmic shear pipelines that have been used to constrain cosmology, we find weak lensing shear measured to an accuracy that is within the statistical errors of current weak lensing analyses, with shear measurements accurate to better than 7%. The dominant source of measurement error is shown to arise from calibration uncertainties where the measured shear is over or under-estimated by a constant multiplicative factor. This is of concern as calibration errors cannot be detected through standard diagnostic tests. The measured calibration errors appear to result from stellar contamination, false object detection, the shear measurement method itself, selection bias and/or the use of biased weights. Additive systematics (false detections of shear) resulting from residual point-spread function anisotropy are, in most cases, reduced to below an equivalent shear of 0.001, an order of magnitude below cosmic shear distortions on the scales probed by current surveys. Our results provide a snapshot view of the accuracy of current ground-based weak lensing methods and a benchmark upon which we can improve. To this end we provide descriptions of each method tested and include details of the eight different implementations of the commonly used Kaiser, Squires and Broadhurst (1995) method (KSB+) to aid the improvement of future KSB+ analyses.

We present results of a microlensing survey toward the Andromeda Galaxy (M31) carried out during four observing seasons at the Isaac Newton Telescope (INT). This survey is part of the larger microlensing survey toward M31 performed by the Microlensing Exploration of the Galaxy and Andromeda (MEGA) collaboration. Using a fully automated search algorithm, we indentify 14 candidate microlensing events, three of which are reported here for the first time. Observations obtained at the Mayall telescope are combined with the INT data to produce composite lightcurves for these candidates. The results from the survey are compared with theoretical predictions for the number and distribution of events. These predictions are based on a Monte Carlo calculation of the detection efficiency and disk-bulge-halo models for M31. The models provide the full phase-space distribution functions for the lens and source populations and are motivated by dynamical and observational considerations. They include differential extinction and span a wide range of parameter space characterised primarily by the mass-to-light ratios for the disk and bulge. For most models, the observed event rate is consistent with the rate predicted for self-lensing -- a MACHO halo fraction of 30% or higher can be ruled at the 95% confidence level. The event distribution does show a large near-far asymmetry hinting at a halo contribution to the microlensing signal. Two candidate events are located at particularly large projected radii on the far side of the disk. These events are difficult to explain by self lensing and only somewhat easier to explain by MACHO lensing. A possibility is that one of these is due to a lens in a giant stellar stream.

We present deep near-infrared Js, H, and Ks band imaging of a field around MS1054-03, a massive cluster at z=0.83. The observations were carried out with ISAAC at the ESO VLT as part of the Faint InfraRed Extragalactic Survey (FIRES). The total integration time amounts to 25.9h in Js, 24.4h in H, and 26.5h in Ks, divided nearly equally between four pointings covering 5.5'x5.3'. The 3-sigma total limiting AB magnitudes for point sources from the shallowest to deepest pointing are Js=26.0-26.2, H=25.5-25.8, and Ks=25.3-25.7. The effective spatial resolution of the coadded images has FWHM=0.48", 0.46", and 0.52" in Js, H, and Ks. We complemented the ISAAC data with deep optical imaging using existing HST WFPC2 mosaics in the F606W and F814W filters and new U, B and V band data from VLT FORS1. We constructed a Ks-band limited multicolour source catalogue to Ks(total,AB)=25 (about 5-sigma for point sources). The catalogue contains 1858 objects, of which 1663 have eight-band photometry. We describe the observations, data reduction, source detection and photometric measurements method. We present the number counts, colour distributions, and photometric redshifts z_ph of the catalogue sources. We find that our counts at the faint end 22<Ks(AB)<25, with slope dlog(N)/dm=0.20, lie at the flatter end of published counts in other deep fields and are consistent with those we derived in the HDF-South, the other FIRES field. Spectroscopic redshifts z_sp are available for about 330 sources in the MS1054-03 field; comparison between the z_ph and z_sp shows very good agreement, with <|z_sp-z_ph|/(1+z_sp)>=0.078. The MS1054-03 field observations complement our HDF-South data set with nearly five times larger area at about 0.7 brighter magnitudes. [ABRIDGED]

We investigate $HST$/ACS and WFPC2 images at the positions of five candidate microlensing events from a large survey of variability in M31 (MEGA). Three closely match unresolved sources, and two produce only flux upper limits. All are confined to regions of the color-magnitude diagram where stellar variability is unlikely to be easily confused with microlensing. Red variable stars cannot explain these events (although background supernova are possible for two). If these lenses arise in M31's halo, they are due to masses $0.15 < m / M_\odot < 0.49$ (95% certainty, for a $\delta$-function mass distribution), brown dwarfs for disk lenses, and stellar masses for bulge lenses.

We have obtained deep color-magnitude diagrams of the Sagittarius dwarf irregular galaxy with the HST/ACS. The new diagrams unveil for the first time the star formation history of SagDIG. We identify a small population of red horizontal branch stars, which sets the first epoch of star formation in SagDIG back to ~9-10 Gyr ago. This shows that SagDIG has been able to maintain a very low mean metallicity (-2.2<[Fe/H]<-1.9) over a 10 Gyr life time, and most importantly, that all Local Group dwarf galaxies share an ancient (~10 Gyr) initial episode of star formation. Thus, at the moment, IZw18 remains the only exception to this general trend, favoring the so-called ``young galaxy'' scenario. However, a re-analysis of ACS observations of IZw18 suggests that an older RGB population may be present also in IZw18. We also compare the distribution of atomic hydrogen with that of stellar populations of various ages. We find little correlation between stellar populations older than ~1 Gyr and the HI. In particular we fail to find any generation of stars that preferentially lie within the large HI hole. We therefore exclude the possibility that the SagDIG ring-like HI distribution is the result of multiple supernova events. Alternative scenarios have to be explored.

Since the near future should see a rapidly expanding set of probes of the halo masses of individual early-type galaxies, we introduce a convenient parameter for characterising the halo masses from both observational and theoretical results: \dML, the logarithmic radial gradient of the mass-to-light ratio. Using halo density profiles from LCDM simulations, we derive predictions for this gradient for various galaxy luminosities and star formation efficiencies $\epsilon_{SF}$. As a pilot study, we assemble the available \dML\ data from kinematics in early-type galaxies - representing the first unbiassed study of halo masses in a wide range of early-type galaxy luminosities - and find a correlation between luminosity and \dML, such that the brightest galaxies appear the most dark-matter dominated. We find that the gradients in most of the brightest galaxies may fit in well with the LCDM predictions, but that there is also a population of fainter galaxies whose gradients are so low as to imply an unreasonably high star formation efficiency $\epsilon_{SF} > 1$. This difficulty is eased if dark haloes are not assumed to have the standard LCDM profiles, but lower central concentrations.

Using the Planetary Nebula Spectrograph, we have observed and measured the velocities for some 2764 PNe in the disk and halo of the Andromeda galaxy. Preliminary analysis using a basic ring model shows a rotation curve in good agreement with that obtained from HI data out to ~20kpc. Some substructure has also been detected within the velocity field, which can be modeled as the continuation of the tidal-remnant known as the Southern Stream, as it passes through Andromeda's disk.

The Planetary Nebula Spectrograph is a dedicated instrument for measuring radial velocity of individual Planetary Nebulae (PNe) in galaxies. This new instrument is providing crucial data with which to probe the structure of dark halos in the outskirts of elliptical galaxies in particular, which are traditionally lacking of easy interpretable kinematical tracers at large distance from the center. Preliminary results on a sample of intermediate luminosity galaxies have shown little dark matter within 5 ~ R_eff implying halos either not as massive or not as centrally concentrated as CDM predicts (Romanowsky et al. 2003). We briefly discuss whether this is consistent with a systematic trend of the dark matter content with the luminosity as observed in an extended sample of early-type galaxies.

The luminosity-size and stellar mass-size distributions of galaxies out to z~3 is presented. We use very deep near-infrared images of the Hubble Deep Field South in the J_s, H, and K_s bands, taken as part FIRES at the VLT, to follow the evolution of the optical rest-frame sizes of galaxies. For a total of 168 galaxies with K_{s,AB}<23.5, we find that the rest-frame V-band sizes r_{e,V} of luminous galaxies (<L_V>~2x10^{10}h^{-2}L_sun) at 2<z<3 are 3 times smaller than for equally luminous galaxies today. In contrast, the mass-size relation has evolved relatively little: the size at mass <M_*>~2x10^{10}h^{-2}M_sun, has changed by 20(+-20)% since z~2.5. Both results can be reconciled by the fact that the stellar M/L ratio is lower in the luminous high z galaxies than in nearby ones because they have young stellar populations. The lower incidence of large galaxies at z~3 seems to reflect the rarity of galaxies with high stellar mass.

We present a possible orbit for the Southern Stream of stars in M31, which connects it to the Northern Spur. Support for this model comes from the dynamics of planetary nebulae (PNe) in the disk of M31: analysis of a new sample of 2611 PNe obtained using the Planetary Nebula Spectrograph reveals ~20 objects whose kinematics are inconsistent with the normal components of the galaxy, but which lie at the right positions and velocities to connect the two photometric features via this orbit. The satellite galaxy M32 is coincident with the stream both in position and velocity, adding weight to the hypothesis that the stream comprises its tidal debris.

A new instrument is providing crucial data with which to probe the structure of dark halos in elliptical galaxies.

I report on an ongoing programme of proper motion measurements in the Galactic Bulge, using HST images. As a first application, the data are used to derive the vertical gravitational acceleration on the minor axis of the Galaxy. Implications of the measurement for the flattening of the Galactic potential are discussed.

Planetary nebulae (PNe) may be the most promising tracers in the halos of early-type galaxies. We have used multi-object spectrographs on the WHT and the VLT, and the new Planetary Nebula Spectrograph on the WHT, to obtain hundreds of PN velocities in a small sample of nearby galaxies. These ellipticals show weak halo rotation, which may be consistent with ab initio models of galaxy formation, but not with more detailed major merger simulations. The galaxies near L* show evidence of a universal declining velocity dispersion profile, and dynamical models indicate the presence of little dark matter within 5 Reff -- implying halos either not as massive or not as centrally concentrated as CDM predicts.

We have found a correlation betwen the M/L global gradients and the structural parameters of the luminous components of a sample of 19 early-type galaxies. Such a correlation supports the hypothesis that there is a connection between the dark matter content and the evolution of the baryonic component in such systems.

We derive the luminosity function of high-redshift Lyman alpha emitting sources from a deep, blind, spectroscopic survey that utilized strong-lensing magnification by intermediate-redshift clusters of galaxies. We observed carefully selected regions near 9 clusters, consistent with magnification factors generally greater than 10 for the redshift range 4.5<z<6.7. Eleven emission-line candidates were located in the range 2.2<z<5.6 whose identification we justify as Lyman alpha, in most cases via further spectroscopic observations. The selection function we constructed for our survey takes into account our varying intrinsic Lyman alpha line sensitivity as a function of wavelength and sky position. By virtue of the strong magnification factor, we provide constraints on the Lyman alpha luminosity function to unprecedented limits of 10^40 erg/s, corresponding to a star-formation rate of 0.01 Msun/yr. Our cumulative z=5 Lyman alpha luminosity function is consistent with a power law form, n(>L) proportional to L^-1 over 10^41 to 10^42.5 erg/s. When combined with the results of other surveys, limited at higher luminosities, our results suggest evidence for the suppression of star formation in low-mass halos, as predicted in popular models of galaxy formation.

We present the first M31 candidate microlensing events from the Microlensing Exploration of the Galaxy and Andromeda (MEGA) survey. MEGA uses several telescopes to detect microlensing towards the nearby Andromeda galaxy, M31, in order to establish whether massive compact objects are a significant contribution to the mass budget of the dark halo of M31. The results presented here are based on observations with the Isaac Newton Telescope on La Palma, during the 1999/00 and 2000/01 observing seasons. In this data set, 14 variable sources consistent with microlensing have been detected, 12 of which are new and 2 have been reported previously by the POINT-AGAPE group. A preliminary analysis of the spatial and timescale distributions of the candidate events support their microlensing nature. We compare the spatial distributions of the candidate events and of long-period variable stars, assuming the chances of finding a long-period variable and a microlensing event are comparable. The spatial distribution of our candidate microlensing events is more far/near side asymmetric than expected from the detected long-period variable distribution. The current analysis is preliminary and the asymmetry not highly significant, but the spatial distribution of candidate microlenses is suggestive of the presence of a microlensing halo.

The kinematics of the outer parts of three intermediate-luminosity elliptical galaxies have been studied using the Planetary Nebula Spectrograph. The galaxies' velocity dispersion profiles are found to decline with radius; dynamical modeling of the data indicates the presence of little if any dark matter in these galaxies' halos. This surprising result conflicts with findings in other galaxy types, and poses a challenge to current galaxy formation theories.

The pattern speed in one of the fundamental parameters that determines the structure of barred galaxies. This quantity is usually derived from indirect methods or by employing model assumptions. The number of bar pattern speeds derived using the model independent Tremaine and Weinberg technique is still very limited. We present the results of model independent measurements of the bar pattern speed in four galaxies ranging in Hubble type from SB0 to SBbc. Three of the four galaxies in our sample are consistent with bars being fast rotators. The lack of slow bars is consistent with previous observations and suggests that barred galaxies do not have centrally concentrated dark matter halos. This contradicts simulations of cosmological structure formation and observations of the central mass concentration in nonbarred galaxies.

Using deep near-infrared imaging of the Hubble Deep Field South with ISAAC on the Very Large Telescope we find 6 large disk-like galaxies at redshifts z = 1.4-3.0. The galaxies, selected in K_s (2.2 micron), are regular and surprisingly large in the near-infrared (rest-frame optical), with face-on effective radii r_e = 0.65"-0.9" or 5.0-7.5 h_70^-1 kpc in a Lambda-CDM cosmology, comparable to the Milky Way. The surface brightness profiles are consistent with an exponential law over 2-3 effective radii. The WFPC2 morphologies in Hubble Space Telescope imaging (rest-frame UV) are irregular and show complex aggregates of star-forming regions ~2" (~15 h_70^-1 kpc) across, symmetrically distributed around the K_s-band centers. The spectral energy distributions show clear breaks in the rest-frame optical. The breaks are strongest in the central regions of the galaxies, and can be identified as the age-sensitive Balmer/4000 Angstrom break. The most straightforward interpretation is that these galaxies are large disk galaxies; deep NIR data are indispensable for this classification. The candidate disks constitute 50% of galaxies with L_V > 6 x 10^10 h_70^-2 L_sun at z = 1.4-3.0. This discovery was not expected on the basis of previously studied samples. In particular, the Hubble Deep Field North is deficient in large galaxies with the morphologies and profiles we report here.

The basic equations and geometry of gravitational lensing are described, as well as the most important contexts in which it is observed in astronomy: strong lensing, weak lensing and microlensing.

We present deep near-infrared (NIR) Js, H, and Ks-band ISAAC imaging of the WFPC2 field of the HDF-S. The 2.5'x 2.5' high Galactic latitude field was observed with the VLT under the best seeing conditions with integration times amounting to 33.6 hours in Js, 32.3 hours in H, and 35.6 hours in Ks. We reach total AB magnitudes for point sources of 26.8, 26.2, and 26.2 respectively (3 sigma), which make it the deepest ground-based NIR observations to date, and the deepest Ks-band data in any field. The effective seeing of the coadded images is ~0.45" in Js, ~0.48" in H, and ~0.46" in Ks. Using published WFPC2 optical data, we constructed a Ks-limited multicolor catalog containing 833 sources down to Ks,tot ~< 26 (AB), of which 624 have seven-band optical-to-NIR photometry. These data allow us to select normal galaxies from their rest-frame optical properties to high redshift (z ~< 4). The observations, data reduction and properties of the final images are discussed, and we address the detection and photometry procedures that were used in making the catalog. In addition, we present deep number counts, color distributions and photometric redshifts of the HDF-S galaxies. We find that our faint Ks-band number counts are flatter than published counts in other deep fields, which might reflect cosmic variations or different analysis techniques. Compared to the HDF-N, we find many galaxies with very red V-H colors at photometric redshifts 1.95 < z < 3.5. These galaxies are bright in Ks with infrared colors redder than Js-Ks > 2.3 (in Johnson magnitudes). Because they are extremely faint in the observed optical, they would be missed by ultraviolet-optical selection techniques, such as the U-dropout method.

We have observed two blank fields of approximately 30 by 23 arcminutes using the William Herschel Telescope. The fields have been studied as part of the CNOC2 survey, and spectroscopic redshifts are available for 1125 galaxies in the two fields. We measured the lensing signal caused by large scale structure, and found that the result is consistent with current, more accurate measurements. We study the galaxy-galaxy lensing signal of three overlapping samples of lenses, and detect a significant signal in all cases. The estimates for the velocity dispersion of an L* galaxy agree well for the various samples. The best fit singular isothermal sphere model to the ensemble averaged tangential distortion around the galaxies with redshifts yields a velocity dispersion of \sigma_*=130^{+15}_{-17} km/s, in good agreement with other studies. We use a maximum likelihood analysis, where a parameterized mass model is compared to the data, to study the extent of galaxy dark matter halos. Making use of all available data, we find \sigma_*=111\pm12 km/s (68.3% confidence) for a truncated isothermal sphere model in which all galaxies have the same mass-to-light ratio. The value of the truncation parameter is not constrained that well, and we find s_*=260^{+124}_{-73} h^{-1} kpc (68.3% confidence The galaxy-galaxy lensing analysis allows us to estimate the average mass-to-light ratio of the field, which can be used to estimate \Omega_m. The current result, however, depends strongly on the assumed scaling relation for s.

Planetary nebulae are now well established as probes of galaxy dynamics and as standard candles in distance determinations. Motivated by the need to improve the efficiency of planetary nebulae searches and the speed with which their radial velocities are determined, a dedicated instrument - the Planetary Nebulae Spectrograph or PN.S - has been designed and commissioned at the 4.2m William Herschel Telescope. The high optical efficiency of the spectrograph results in the detection of typically ~ 150 PN in galaxies at the distance of the Virgo cluster in one night of observations. In the same observation the radial velocities are obtained with an accuracy of ~ 20 km/s

We have observed planetary nebulae (PNe) in several early-type galaxies using new techniques on 4- to 8-meter-class telescopes. We obtain the first large data sets (~ 100 velocities each) of PN kinematics in galaxies at > 15 Mpc, and present some preliminary dynamical results.

We report here the successful commissioning of the PN.Spectrograph, the first special-purpose instrument for the measurement of galaxy kinematics through the PN population

We report proper motion dispersions for stars in the direction of two fields of the Galactic bulge, using HST/WFPC2 images taken six years apart. Our two fields are Baade's Window (l,b)=(1.13,-3.77) and Sgr I (l,b)=(1.25 ,-2.65). Our proper motion dispersions are in good agreement with prior ground- and space-based proper motion studies in bulge fields, but in contrast to some prior studies, we do not exclude any subset of stars from our studies. In Baade's Window, we find the $l$ and $b$ proper motion dispersions are 2.9 and 2.5 mas/yr, while in Sgr I, they are 3.3 and 2.7 mas/yr, respectively. For the first time, we can clearly separate the foreground disk stars out from the bulge because of their large mean apparent proper motion. The population with non-disk kinematics (which we conclude to be the bulge) has an old main sequence turnoff point, similar to those found in old, metal rich bulge globular clusters while those stars selected to have disk kinematics lie on a fully populated main sequence. Separating main sequence stars by luminosity, we find strong evidence that the bulge population is rotating, largely explaining observations of proper motion anisotropy in bulge fields. Because we have isolated such a pure sample of stars in the bulge, we have one of the clearest demonstrations that the old stellar population of the inner bulge/bar is in fact rotating.

We present 21-cm HI line and optical R-band observations for a sample of 26 edge-on galaxies. The HI observations were obtained with the Westerbork Synthesis Radio Telescope, and are part of the WHISP database (Westerbork HI Survey of Spiral and Irregular Galaxies). We present HI maps, optical images, and radial HI density profiles. We have also derived the rotation curves and studied the warping and lopsidedness of the HI disks. 20 out of the 26 galaxies of our sample are warped, confirming that warping of the HI disks is a very common phenomenon in disk galaxies. Indeed, we find that all galaxies that have an extended HI disk with respect to the optical are warped. The warping usually starts around the edge of the optical disk. The degree of warping varies considerably from galaxy to galaxy. Furthermore, many warps are asymmetric, as they show up in only one side of the disk or exhibit large differences in amplitude in the approaching and receding sides of the galaxy. These asymmetries are more pronounced in rich environments, which may indicate that tidal interactions are a source of warp asymmetry. A rich environment tends to produce larger warps as well. The presence of lopsidedness seems to be related to the presence of nearby companions.

This paper investigates the detailed dynamical properties of a relatively homogeneous sample of disc-dominated S0 galaxies, with a view to understanding their formation, evolution and structure. By using high signal-to-noise ratio long-slit spectra of edge-on systems, we have been able to reconstruct the complete line-of-sight velocity distributions of stars along the galaxies' major axes. From these data, we have derived both model distribution functions (the phase density of their stars) and the approximate form of their gravitational potentials. The derived distribution functions are all consistent with these galaxies being simple disc systems, with no evidence for a complex formation history. Essentially no correlation is found between the characteristic mass scale-lengths and the photometric scale-lengths in these galaxies, suggesting that they are dark-matter dominated even in their inner parts. Similarly, no correlation is found between the mass scale-lengths and asymptotic rotation speed, implying a wide range of dark matter halo properties. By comparing their asymptotic rotation speeds with their absolute magnitudes, we find that these S0 galaxies are systematically offset from the Tully-Fisher relation for later-type galaxies. The offset in luminosity is what one would expect if star formation had been suddenly switched off a few Gyrs ago, consistent with a simple picture in which these S0s were created from ordinary later-type spirals which were stripped of their star-forming ISM when they encountered a dense cluster environment.

We have detected the weak lensing signal induced by the cluster of galaxies MS 2053-04 (z=0.58) from a two-colour mosaic of 6 HST WFPC2 images. The best fit singular isothermal sphere model to the observed tangential distortion yields an Einstein radius r_E=6.2"+-1.8", which corresponds to a velocity dispersion of 886^{+121}_{-139} km/s. This result is in good agreement with the observed velocity dispersion of 817+-80 km/s from cluster members. MS 2053 is the third cluster we studied using mosaics of deep WFPC2 images. For all three clusters we find good agreement between dynamical and weak lensing velocity dispersions, in contrast to weak lensing studies based on single WFPC2 pointings on cluster cores. This result demonstrates the importance of wide field data. We have compared the ensemble averaged cluster profile to the predicted NFW profile, and find that a NFW profile can fit the observed lensing signal well. The best fit concentration parameter is found to be 0.79^{+0.44}_{-0.15} (68% confidence) times the predicted value from an open CDM model. The observed mass-to-light ratios of the clusters in our sample evolve with redshift, and are inconsistent with a constant, non-evolving, mass-to-light ratio at the 99% confidence level. The evolution is consistent with the results derived from the evolution of the fundamental plane of early type galaxies. The resulting average mass-to-light ratio for massive clusters at z=0 is found to be 239+-18+-9 M/L_B. (abridged)

We discuss the physical nature of a remarkably faint pair of Lyman alpha-emitting images discovered close to the giant cD galaxy in the lensing cluster Abell 2218 (z=0.18) during a systematic survey for highly-magnified star-forming galaxies beyond z=5. A well-constrained mass model suggests the pair arises via a gravitationally-lensed source viewed at high magnification. Keck spectroscopy confirms the lensing hypothesis and implies the unlensed source is a very faint (I~30) compact (<150 pc) and isolated object at z=5.576 whose optical emission is substantially contained within the Lyman alpha emission line; no stellar continuum is detectable. The available data suggest the source is a promising candidate for an isolated ~10^6 solar mass system seen producing its first generation of stars close to the epoch of reionization.

One of the possible astrophysical solutions to the galactic dark matter problem is the presence of a significant amount of ``dark'' compact objects (MACHOs) in galactic dark matter halos. MEGA (Microlensing Exploration of the Galaxy and Andromeda) tries to find proof for or against the presence of compact objects in the halo of the Andromeda galaxy (M31) by looking for the microlensing signature that would be induced by these objects. The lightcurves presented here are preliminary and based on observations of M31 with the Isaac Newton Telescope (INT) on La Palma during the second half of 1999.

We present the first results from the Faint Infra-Red Extragalactic Survey (FIRES) of the Hubble Deep Field South. Using a combination of deep near infrared data obtained with ISAAC at the VLT with the WFPC2 HST data, we construct a K-band selected sample of 136 galaxies with K(AB)<=23.5. We derive zphot's using a newly developed technique which models the observed spectral energy distribution with a linear combination of empirical galaxy templates. Testing our zphot technique against spectroscopic redshifts in the HDF-N yields dz/(1+z)~0.07 for z<6. We show that we can derive realistic error estimates in zphot by combining the systematic uncertainties derived from the HDF-N with errors in zphot which depend on the observed flux errors. The redshift histogram of galaxies in the HDF-S shows distinct structure with a sharp peak at z~0.5 and a broad enhancement at z~1-1.4. We find that 12% of our galaxies with K(vega)<21 lie at z>=2. While this is higher than the fraction predicted in Omega_M=1 hierarchical models of galaxy formation we find that published predictions using pure luminosity evolution models produce too many bright galaxies at redshifts greater than unity. Finally, we use our broad wavelength coverage to measure the rest-frame UBV luminosities Lrest for z<=3. There is a paucity of galaxies brighter than Lrest_V>=1.4e10 h^{-2}Lsun at z~1.5-2, however, at z>2 we find very luminous galaxies with Lrest_V>=5e10 h^{-2}Lsun (for Omega_M=0.3,Lambda=0.7, H_o=100 h km s^{-1}Mpc^{-1}). Local B-band luminosity functions predict 0.1 galaxies in the redshift range 2<=z<=3.5 and with Lrest_B>= 5e10 h^{-2}Lsun,B but we find 9. The discrepancy can be explained if L*_B increases by a factor of 2.4-3.2 with respect to locally determined values (abridged).

We report on the first measurement of the average mass and mass-to-light ratio of galaxy groups by analysing the weak lensing signal induced by these systems. The groups, which have velocity dispersions of 50-400 km/s, have been selected from the Canadian Network for Observational Cosmology Field Galaxy Redshift Survey (CNOC2). This survey allows the identification of a large number of groups with redshifts ranging from z=0.12-0.55, ideal for a weak lensing analysis of their mass distribution. For our analysis we use a sample of 50 groups which are selected on the basis of a careful dynamical analysis of group candidates. We detect a signal at the 99% confidence limit. The best fit singular isothermal sphere model yields an Einstein radius of 0.72+-0.29". This corresponds to a velocity dispersion of 274^{+48}_{-59} km/s (using photometric redshift distributions for the source galaxies), which is in good agreement with the dynamical estimate. Under the assumption that the light traces the mass, we find an average mass-to-light ratio of 191+-83 h Msun/Lsun in the restframe B band. Unlike dynamical estimates, this result is insensitive to problems associated with determining group membership. After correction of the observed mass-to-light ratio for luminosity evolution to z=0, we find 254+-110 h Msun/Lsun, lower than what is found for rich clusters. We use the observed mass-to-light ratio to estimate the matter density of the universe, for which we find Omega_m=0.19+-0.10 (Omega_Lambda=0), in good agreement with other recent estimates. For a closed universe, we obtain Omega_m=0.13+-0.07.

This review addresses recent developments in the field of disk galaxy warps. Both results from a new HI survey of edgeon disk galaxies, and of simulations of the interaction between a disk+halo and an orbiting satelite, will be discussed.

Recent results of microlensing surveys, show that 10-20% of the dark halo mass of the Milky Way consists of compact objects. The masses of these compact objects range from 0.1 to 1 solar mass. New theoretical cooling models for white dwarfs, and the detection of faint blue high proper motion objects, imply that ancient white dwarfs might make up part of this population. In this pilot project, using data from the C and D patches of the ESO Imaging Survey and data obtained at CTIO, we attempted to find such halo white dwarfs in the solar neighbourhood. With a time baseline of approximately one year between data sets and a limiting I-band magnitude of 23, we can find these high proper motion objects up to distances of 85 parsecs. In the 2.5 square degrees studied so far we have found three high proper motion candidates.

By making spectral absorption-line observations of disk galaxies at intermediate inclinations, we have been able to determine the amplitude of their constituent stars' random motions in three dimensions. This full measure of the shape of the velocity ellipsoid is a useful diagnostic for determining the "heating" mechanism responsible for creating the random motions. The analysis implies that the main heating process varies with galaxy type, with minor mergers dominating the heating in early-type disk galaxies, and spiral density waves the most important mechanism in late-type systems.

Now that weak lensing signals on the order of a percent are actively being searched for (cosmic shear, galaxy-galaxy lensing, large radii in clusters...) it is important to investigate how accurately weak shears can be determined. Many systematic effects are present, and need to be understood. I show that the Kaiser et al. technique can leave residual systematic errors at the percent level (through imperfect PSF anisotropy correction), and present an alternative technique which is able to recover shears a factor of ten weaker.

January 2000 completes the first season of intensive, wide-field observations of microlensing and stellar variability in M31 by MEGA (Microlensing Exploration of the Galaxy and Andromeda) at the Isaac Newton 2.5m Telescope, the KPNO 4m, and the 1.3m and 2.5m telescopes of MDM Observatory. In preliminary analysis, we detect ~50000 variable objects, including some consistent with microlensing events. We present the level of sensitivity to be reached in our planned three-year program to test for the presence of a significant halo microlensing population in M31, as well as its spatial distribution and mass-function. We also discuss our application of image subtraction to these wide fields and HST WFPC2 Snapshot followup observations to confirm candidates identified from previous years' surveys. We present intermediate results from our smaller-field survey, on the MDM 1.3m and Vatican Advanced Technology 1.8m Telescope, from 1994-1998, wherein we have discovered 8 additional probable microlensing events, over about one-half the time base of the project, in addition to confirming three of our original 6 microlensing candidates from 1995.

Various processes have been proposed to explain how galaxy discs acquire their thickness. A simple diagnostic for ascertaining this ``heating'' mechanism is provided by the ratio of the vertical to radial velocity dispersion components. In a previous paper we have developed a technique for measuring this ratio, and demonstrated its viability on the Sb system NGC 488. Here we present follow-up observations of the morphologically similar Sab galaxy NGC 2985, still only the second galaxy for which this ratio has been determined outside of the solar neighbourhood. The result is consistent with simple disc heating models which predict ratios of $\sigma_z / \sigma_R$ less than one

The planetary nebula populations of relatively nearby galaxies can be easily observed and provide both a distance estimate and a tool with which dynamical information can be obtained. Usually the requisite radial velocities are obtained by multi-object spectroscopy once the planetary nebulae have been located by direct imaging. Here we report on a technique for measuring planetary nebula kinematics using the double-beam ISIS spectrograph at the William Herschel Telescope in a novel slitless mode, which enables the detection and radial velocity measurements to be combined into a single step. The results on our first target, the Sab galaxy NGC 4736, allow the velocity dispersion of the stellar population in a disk galaxy to be traced out to four scale lengths for the first time and are consistent with a simple isothermal sheet model.

After studying the orientation of a warp generated by a companion satellite, we show that the Galactic Warp would be oriented in a different way if the Magellanic Clouds were its cause. We have treated the problem analytically, and complemented it with numerical N-Body simulations.

We present a gravitational lensing study of the massive galaxy cluster A2219 (redshift 0.22). This investigation is based on multicolour images from U through H, which allows photometric redshifts to be estimated for the background sources. The redshifts provide useful extra information for the lensing models: we show how they can be used to identify a new multiple-image system (and rule out an old one), how this information can be used to anchor the mass model for the cluster, and how the redshifts can be used to construct optimal samples of background galaxies for a weak lensing analysis. Combining all results, we obtain the mass distribution in this cluster from the inner, strong lensing region, out to a radius of 1.5 Mpc. The mass profile is consistent with a singular isothermal model over this radius range. Parametric and non-parametric reconstructions of the mass distribution in the cluster are compared. The main features (elongation, sub-clumps, radial mass profile) are in good agreement.

We have measured for the first time the weak lensing signal due to groups of galaxies. The groups are at intermediate redshifts, and have been identified in the CNOC2 field survey, which was kindly made available by Ray Carlberg and Howard Yee. The ensemble averaged group velocity dispersion, based on a preliminary selection of 59 groups, is found to be sigma= 320^{+46}_{-54} km/s, which is in fair agreement with the dynamical estimate from the spectroscopic redshifts. Under the assumption that mass traces light, we find an average mass-to-light ratio in restframe B of (256\pm84)h Msun/Lsun. Using this result we obtain Omega_m=0.22+-0.09 (\Omega_\Lambda=0).

We have measured the weak gravitational lensing signal of MS 1054-03, a rich and X-ray luminous cluster of galaxies at a redshift of z=0.83, using a two-colour mosaic of deep WFPC2 images. The small corrections for the size of the PSF and the high number density of background galaxies obtained in these observations result in an accurate and well calibrated measurement of the lensing induced distortion. The strength of the lensing signal depends on the redshift distribution of the background galaxies. We used photometric redshift distributions from the Northern and Southern Hubble Deep Fields to relate the lensing signal to the mass. The predicted variations of the signal as a function of apparent source magnitude and colour agrees well with the observed lensing signal. We determine a mass of (1.2+-0.2)x10^15 Msun within an aperture of radius 1 Mpc. Under the assumption of an isothermal mass distribution, the corresponding velocity dispersion is 1311^{+83}_{-89} km/s. For the mass-to-light ratio we find 269+-37 Msun/Lsun. The errors in the mass and mass-to-light ratio include the contribution from the random intrinsic ellipticities of the source galaxies, but not the (systematic) error due to the uncertainty in the redshift distribution. However, the estimates for the mass and mass-to-light ratio of MS 1054-03 agree well with other estimators, suggesting that the mass calibration works well. The reconstruction of the projected mass surface density shows a complex mass distribution, consistent with the light distribution. The results indicate that MS 1054-03 is a young system. The timescale for relaxation is estimated to be at least 1 Gyr. Averaging the tangential shear around the cluster galaxies, we find that the velocity dispersion of an Lstar galaxy is 203+-33 km/s.

We have developed a new technique for weak lensing analysis, with which the effect of the point spread function (PSF) on small galaxy images can be corrected for accurately. Rather than relying on weighted second moments of detected images, which we show can leave residuals at the level of a percent in the shear, we directly fit (stacked or individual) galaxy images as PSF-convolved, sheared circular sources. We show by means of simulations that this technique is able to recover shears well below the percent level for a variety of PSF shapes, and that its noise properties are similar to existing methods.

It has been suggested that the boxy and peanut-shaped bulges found in some edge-on galaxies are galactic bars viewed from the side. We investigate this hypothesis by presenting emission-line spectra for a sample of 10 edge-on galaxies that display a variety of bulge morphologies. To avoid potential biases in the classification of this morphology, we use an objective measure of bulge shape. Generally, bulges classified as more boxy show the more complicated kinematics characteristic of edge-on bars, confirming the intimate relation between the two phenomena.

The disks of galaxies are primarily stellar systems, and fundamentally dynamical entities. Thus, to fully understand galactic disks, we must study their stellar kinematics as well as their morphologies. Observational techniques have now advanced to a point where quite detailed stellar-kinematic information can be extracted from spectral observations. This review presents three illustrative examples of analyses that make use of such information to study the formation and evolution of these systems: the derivation of the pattern speed of the bar in NGC 936; the calculation of the complete velocity ellipsoid of random motions in NGC 488; and the strange phenomenon of counter-rotation seen in NGC 3593.

We report on the detection of weak gravitational lensing of faint, distant background galaxies by Cl 1358+62, a cluster of galaxies at a redshift of z=0.33. The observations were made using the HST. The measured shear is consistent with a velocity dispersion of 780+-50 km/s. The weak lensing mass is slightly lower than dynamical estimates and agrees well with X-ray mass estimates. The mass distribution is elongated similar to the light. The axis ratio of 0.30+-0.15 and position angle of -21+-8 degrees were measured directly from the observed shear and agree very well with a previous strong lensing analysis. We estimate the mass-to-light ratio to be 90+-13 h50 M_sun/L_Vsun. The HST point spread function is highly anisotropic at the edges of the individual chips. This systematically perturbs the shapes of objects and we present a method for applying the appropriate correction.

We show that the kinematics of the shells seen around some elliptical galaxies provide a new, independent means for measuring the gravitational potentials of elliptical galaxies out to large radii. A numerical simulation of a set of shells formed in the merger between an elliptical and a smaller galaxy reveals that the shells have a characteristic observable kinematic structure, with the maximum line-of-sight velocity increasing linearly as one moves inward from a shell edge. A simple analytic calculation shows that this structure provides a direct measure of the gradient of the gravitational potential at the shell radius. In order to extract this information from attainable data, we have also derived the complete distribution of line-of-sight velocities for material within a shell; comparing the observed spectra of a shell to a stellar template convolved with this distribution will enable us to measure the gradient of the potential at this radius. Repeating the analysis for a whole series of nested shells in a galaxy allows the complete form of the gravitational potential as a function of radius to be mapped out. The requisite observations lie within reach of the up-coming generation of large telescopes.

Secondary infall in galaxies could cause the angular momentum of the outer halo to change its orientation. The generation of warps due to this effect is studied using N-body simulations.

It has been known that recent microlensing results towards the bulge imply mass densities that are surprisingly high given dynamical constraints on the Milky Way mass distribution. We derive the maximum optical depth towards the bulge that may be generated by axisymmetric structures in the Milky Way, and show that observations are close to surpassing these limits. This result argues in favor of a bar as a source of significantly enhanced microlensing. Several of the bar models in the literature are discussed.

Theories of stellar orbit diffusion in disk galaxies predict different rates of increase of the velocity dispersions parallel and perpendicular to the disk plane, and it is therefore of interest to measure the different velocity dispersion components in galactic disks of different types. We show that it is possible to extract the three components of the velocity ellipsoid in an intermediate-inclination disk galaxy from measured line-of-sight velocity dispersions on the major and minor axes. On applying the method to observations of the Sb galaxy NGC 488, we find evidence for a higher ratio of vertical to radial dispersion in NGC 488 than in the solar neighbourhood of the Milky Way (the only other place where this quantity has ever been measured). The difference is qualitatively consistent with the notion that spiral structure has been relatively less important in the dynamical evolution of the disk of NGC 488 than molecular clouds.

The normal mode theory for warping of galaxy disks, in which disks are assumed to be tilted with respect to the equator of a massive, flattened dark halo, assumes a rigid, fixed halo. However, consideration of the back-reaction by a misaligned disk on a massive particle halo shows there to be strong coupling leading to efficient damping (or in some circumstances excitation) of the misalignment, and hence the warp. We therefore discuss possible alternative explanations of the warp phenomenon, with emphasis on the effect of a responsive, gravitationally live massive galactic halo.

We have obtained high signal-to-noise spectra along the major axes of 28 S0 galaxies in order to search for the presence of disk stars on retrograde orbits. Full line-of-sight velocity distributions were extracted from the data, and the velocity distributions were modelled as arising from the superposition of populations of stars on prograde and retrograde orbits. We find no new cases in which a significant fraction of disk stars lie on retrograde orbits; an identical analysis of NGC~4550 does reveal the previously-known counter-rotating stellar disk in this system. Upper limits determined for each object indicate that no more than about 5% of the observed disk star light could arise from counter-rotating stellar components. These results suggest that previously-discovered disk galaxies with counter-rotating stars are exceptional and that (at 95% confidence) at most 10% of S0 galaxies contain significant counter-rotating populations. The most likely value for the fraction of such S0 galaxies lies closer to 1%. This result contrasts with the prevalence of counter-rotating gas in these systems; combining our new observations with existing data, we find that 24+/-8% (1-sigma error) of the gas disks in S0 galaxies counter-rotate relative to their stellar components.

Evidence from a variety of sources points towards the existence of a bar in the central few kpc of the Galaxy. The measurements roughly agree on the direction of the bar major axis, but other parameters (axis ratio, size, pattern speed) are still poorly determined. Current dynamical models are limited by the quality of hydro simulations, the degeneracy of stellar orbit models, stellar-kinematic data and the significant lopsidedness of the central kpc. Microlensing promises new constraints on the mass distribution in the bulge/bar region.

Using a spectrum obtained under moderate (of order 1 arcsecond) seeing, we show that the double nucleus in M31 produces a strong kinematic signature even though the individual components are not spatially resolved. The signature consists of a significant asymmetric wing in the stellar velocity distribution close to the center of the system. The properties of the second nucleus derived from this analysis agree closely with those measured from high-spatial resolution Hubble Space Telescope images. Even Space Telescope only has sufficient resolution to study the structure of very nearby galactic nuclei photometrically; this spectroscopic approach offers a tool for detecting structure such as multiple nuclei in a wider sample of galaxy cores.

The recent observations of high-velocity maser emissions in NGC 4258 by Myoshi et al provide strong evidence for a central object of mass 3.6 10**7 solar masses in the nucleus of this galaxy based on a kinematic argument. However, they should allow an independent determination of the central object mass via gravitational redshift. We briefly discuss the size of the effect, and show that it should be detected in a more detailed analysis of the data.

We have used the Tremaine-Weinberg method to measure the angular speed of rotation for the bar in the SB0 galaxy NGC 936. With this technique, the bar's pattern speed, Omega_p, can be derived from the luminosity and stellar-kinematic information in long-slit spectral observations taken parallel to the major axis of the galaxy. The kinematic measurement required is the mean line-of-sight velocity of all stellar light entering the slit. This quantity can only be calculated reliably if any asymmetry in the shape of the broadening function of the spectral lines is also measured, and so we present a method which allows for such asymmetry. The technique also returns a true measure of the RMS uncertainty in the estimate. Application of the analysis to a set of long-slit spectra of NGC 936 returns four separate measures of Omega_p which are mutually consistent. Combining these data produces a best estimate for the bar pattern speed of Omega_p = 60 +/- 14 km/s/kpc (assuming a distance of 16.6 Mpc). This result refines the only previous attempt to make this measurement, which yielded an estimate for Omega_p in NGC 936 of 104 +/- 37 km/s/kpc (Kent 1987). The new measurement places the co-rotation radius just beyond the end of the bar, in agreement with theoretical calculations.

We describe methods for setting up self-consistent disk-bulge-halo galaxy models. The bulge and halo distribution functions (\df) are functions of $E$ and $L_z$ only. The halo's flattening and rotation can be specified. The disk \df\ is a function of $E$ and $L_z$ and a third ``integral'', $E_z$, the vertical energy, which is approximately conserved in a warm disk with vertical extent. The models also have finite extent making them suitable for N-body simulation. A simulation of a sample model shows that in practice the models are very close to equilibrium making them ideal for experiments on instabilities in galactic disks. We also present a sequence of models closely resembling the Milky Way mass distribution with 5 exponential scale radii and varying halo mass and radial extent.

It has been suggested that the peanut-shaped bulges seen in some edge-on disk galaxies are due to the presence of a central bar. Although bars cannot be detected photometrically in edge-on galaxies, we show that barred potentials produce a strong kinematic signature in the form of double-peaked line-of-sight velocity distributions with a characteristic ``figure-of-eight'' variation with radius. We have obtained spectroscopic observations of two edge-on galaxies with peanut-shaped bulges (NGC~5746 and NGC~5965), and they reveal exactly such line-of-sight velocity distributions in both their gaseous (emission line) and their stellar (absorption line) components. These observations provide strong observational evidence that peanut-shaped bulges are a by-product of bar formation.

When a galaxy forms, the disk may initially be tilted with respect to a flattened dark halo. The misalignment between the disk and the halo is a common explanation for galactic disk warps, since in this state, disks have precessing bending modes which resemble real warps. A dubious assumption in this theory is that the gravitational response of the halo is negligible. We therefore calculate the response of an oblate halo to a precessing inclined exponential disk using a variety of techniques. We construct models with a rigid exponential disk precessing in a particle halo, a particle disk precessing inside a static bulge/halo potential, and a self-consistent model with a particle disk, bulge and halo. When the disk:halo mass ratio is small ($\sim$10\%) within 5 exponential scale radii, the disk settles to the equatorial plane of the halo within 5 orbital times. When the disk and halo mass are comparable, the halo rapidly aligns with the disk within a few orbital times, while the disk inclination drops. The rapid response of the halo to a inclined disk suggests that the warps seen in galactic disks are not likely the bending modes of a precessing disk inside of triaxial halos. If a galaxy forms inclined to the principal plane of a dark halo, either the disk will settle to a principal plane or the inner halo will twist to align with the disk. The outer halo will remain misaligned for a much longer time and therefore may still exert a torque. Warped bending modes may still exist if the misalignment of the outer halo persists for a Hubble time.

We describe a method for parameterizing two-integral distribution functions, based on triangular tesselations of the integral plane. We apply the method to the axisymmetric isotropic rotator model for the Galactic bulge of Kent~(1992), and compare the results with observations of proper motions in Baade's Window, and with radial velocity surveys. In spite of mounting evidence from surface photometry and from study of the gas kinematics that the Galactic bulge is not axisymmetric, the stellar kinematics in Baade's Window are very similar to those of an isotropic oblate rotator. Another field at large radius does not fit this model, though. In any case, the edge-on kinematics of a hot stellar population are a poor handle on the existence or otherwise of a bar.

We describe an analytic distribution function of a finite, oblate stellar system that is useful for the practical modelling of dark halos. The function is determined by lowering Evans's (1993) distribution function of a flattened, cored isothermal system in analogy to the lowering of the singular, isothermal sphere in the definition of the King (1966) models. We derive analytic expressions of the density, maximal streaming velocity and velocity dispersion profiles as a function of the potential. As for the King models, the potential must be calculated numerically. We also present a recipe for generating N-body realizations of this distribution function and examine the stability in three models with dimensionless spins $\lambda=$0.0, 0.05 and 0.18 using N-body simulations with 50,000 particles. The $\lambda=0.18$ model is unstable to the formation of a triaxial bar within $\sim 5$ King radii while the other models appear stable. We conclude that the slowly rotating systems are useful for modelling flattened dark halos.