We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicative M(l) systematics typically reported in current lensing measurements. We find that overall performance is driven by a product of a telescope/camera's *absolute performance*, and our *knowledge about its performance*. The second half of this paper propagates any residual shear measurement biases through to their effect on cosmological parameter constraints. Fully exploiting the statistical power of Stage IV weak lensing surveys will require additive biases A<1.8e-12 and multiplicative biases M<4.0e-3. These can be allocated between individual budgets in hardware, calibration data and software, using results from the first half of the paper. If instrumentation is stable and well-calibrated, we find extant shear measurement software from GREAT10 already meet requirements on galaxies detected at S/N=40. Averaging over a population of galaxies with a realistic distribution of sizes, it also meets requirements for a 2D cosmic shear analysis from space. If used on fainter galaxies or for 3D cosmic shear tomography, existing algorithms would need calibration on simulations to avoid introducing bias at a level similar to the statistical error. Requirements on hardware and calibration data are discussed in more detail in a companion paper. Our analysis is intentionally general, but is specifically being used to drive the hardware and ground segment performance budget for the design of the European Space Agency's recently-selected Euclid mission.

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.

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

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.

Searching for microlensing in M31 using automated superpixel surveys raises a number of difficulties which are not present in more conventional techniques. Here we focus on the problem that the list of microlensing candidates is sensitive to the selection criteria or "cuts" imposed and some subjectivity is involved in this. Weakening the cuts will generate a longer list of microlensing candidates but with a greater fraction of spurious ones; strengthening the cuts will produce a shorter list but may exclude some genuine events. We illustrate this by comparing three analyses of the same data-set obtained from a 3-year observing run on the INT in La Palma. The results of two of these analyses have been already reported: Belokurov et al. (2005) obtained between 3 and 22 candidates, depending on the strength of their cuts, while Calchi Novati et al. (2005) obtained 6 candidates. The third analysis is presented here for the first time and reports 10 microlensing candidates, 7 of which are new. Only two of the candidates are common to all three analyses. In order to understand why these analyses produce different candidate lists, a comparison is made of the cuts used by the three groups...

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.

Cosmic shear measurements rely on our ability to measure and correct the Point Spread Function (PSF) of the observations. This PSF is measured using stars in the field, which give a noisy measure at random points in the field. Using Wiener filtering, we show how errors in this PSF correction process propagate into shear power spectrum errors. This allows us to test future space-based missions, such as Euclid or JDEM, thereby allowing us to set clear engineering specifications on PSF variability. For ground-based surveys, where the variability of the PSF is dominated by the environment, we briefly discuss how our approach can also be used to study the potential of mitigation techniques such as correlating galaxy shapes in different exposures. To illustrate our approach we show that for a Euclid-like survey to be statistics limited, an initial pre-correction PSF ellipticity power spectrum, with a power-law slope of -3 must have an amplitude at l =1000 of less than 2 x 10^{-13}. This is 1500 times smaller than the typical lensing signal at this scale. We also find that the power spectrum of PSF size \dR^2) at this scale must be below 2 x 10^{-12}. Public code available as part of iCosmo at http://www.icosmo.org

We investigate the impact of point spread function (PSF) fitting errors on cosmic shear measurements using the concepts of complexity and sparsity. Complexity, introduced in a previous paper, characterizes the number of degrees of freedom of the PSF. For instance, fitting an underlying PSF with a model with low complexity will lead to small statistical errors on the model parameters, however these parameters could suffer from large biases. Alternatively, fitting with a large number of parameters will tend to reduce biases at the expense of statistical errors. We perform an optimisation of scatters and biases by studying the mean squared error of a PSF model. We also characterize a model sparsity, which describes how efficiently the model is able to represent the underlying PSF using a limited number of free parameters. We present the general case and illustrate it for a realistic example of PSF fitted with shapelet basis sets. We derive the relation between complexity and sparsity of the PSF model, signal-to-noise ratio of stars and systematic errors on cosmological parameters. With the constraint of maintaining the systematics below the statistical uncertainties, this lead to a relation between the required number of stars to calibrate the PSF and the sparsity. We discuss the impact of our results for current and future cosmic shear surveys. In the typical case where the biases can be represented as a power law of the complexity, we show that current weak lensing surveys can calibrate the PSF with few stars, while future surveys will require hard constraints on the sparsity in order to calibrate the PSF with 50 stars.

The control of systematic effects when measuring galaxy shapes is one of the main challenges for cosmic shear analyses. In this context, we study the fundamental limitations on shear accuracy due to the measurement of the Point Spread Function (PSF) from the finite number of stars. In order to do that, we translate the accuracy required for cosmological parameter estimation to the minimum number of stars over which the PSF must be calibrated. We first derive our results analytically in the case of infinitely small pixels (i.e. infinitely high resolution). Then image simulations are used to validate these results and investigate the effect of finite pixel size in the case of an elliptical gaussian PSF. Our results are expressed in terms of the minimum number of stars required to calibrate the PSF in order to ensure that systematic errors are smaller than statistical errors when estimating the cosmological parameters. On scales smaller than the area containing this minimum number of stars, there is not enough information to model the PSF. In the case of an elliptical gaussian PSF and in the absence of dithering, 2 pixels per PSF Full Width at Half Maximum (FWHM) implies a 20% increase of the minimum number of stars compared to the ideal case of infinitely small pixels; 0.9 pixels per PSF FWHM implies a factor 100 increase. In the case of a good resolution and a typical Signal-to-Noise Ratio distribution of stars, we find that current surveys need the PSF to be calibrated over a few stars, which may explain residual systematics on scales smaller than a few arcmins. Future all-sky cosmic shear surveys require the PSF to be calibrated over a region containing about 50 stars.

We present a joint weak lensing and X-ray analysis of 4 deg$^2$ from the CFHTLS and XMM-LSS surveys. Our weak lensing analysis is the first analysis of a real survey using shapelets, a new generation weak lensing analysis method. We create projected mass maps of the images, and extract 6 weak-lensing-detected clusters of galaxies. We show that their counts can be used to constrain the power spectrum normalisation $\sigma_8 =0.92_{-0.30}^{+0.26}$ for $\Omega_m=0.24$. We show that despite the large scatter generally observed in the M-T relation derived from lensing masses, tight constraints on both its slope and normalisation $M_*$ can be obtained with a moderate number of sources provided that the covered mass range is large enough. Adding clusters from Bardeau et al. (2007) to our sample, we measure $M_* = 2.71_{-0.61}^{+0.79} 10^{14} h^{-1} M_\odot$. Although they are dominated by shot noise and sample variance, our measurements are consistent with currently favoured values, and set the stage for future surveys. We thus investigate the dependence of those estimates on survey size, depth, and integration time, for joint weak lensing and X-ray surveys. We show that deep surveys should be dedicated to the study of the physics of clusters and groups of galaxies. For a given exposure time, wide surveys provide a larger number of detected clusters and are therefore preferred for the measurement of cosmological parameters such as $\sigma_8$ and $M_*$. We show that a wide survey of a few hundred square degrees is needed to improve upon current measurements of these parameters. More ambitious surveys covering 7000 deg$^2$ will provide the 1% accuracy in the estimation of the power spectrum and the M-T relation normalisations.

Weak lensing applied to deep optical images of clusters of galaxies provides a powerful tool to reconstruct the distribution of the gravitating mass associated to these structures. We use the shear signal extracted by an analysis of deep exposures of a region centered around the galaxy cluster Abell 209, at redshift z=0.2, to derive both a map of the projected mass distribution and an estimate of the total mass within a characteristic radius. We use a series of deep archival R-band images from CFHT-12k, covering an area of 0.3 deg^2. We determine the shear of background galaxy images using a new implementation of the modified Kaiser-Squires-Broadhurst pipeline for shear determination, which we has been tested against the ``Shear TEsting Program 1 and 2'' simulations. We use mass aperture statistics to produce maps of the 2 dimensional density distribution, and parametric fits using both Navarro-Frenk-White (NFW) and singular-isothermal-sphere profiles to constrain the total mass. The projected mass distribution shows a pronounced asymmetry, with an elongated structure extending from the SE to the NW. This is in general agreement with the optical distribution previously found by other authors. A similar elongation was previously detected in the X-ray emission map, and in the distribution of galaxy colours. The circular NFW mass profile fit gives a total mass of M_{200} = 7.7^{+4.3}_{-2.7} 10^{14} solar masses inside the virial radius r_{200} = 1.8\pm 0.3 Mpc. The weak lensing profile reinforces the evidence for an elongated structure of Abell 209, as previously suggested by studies of the galaxy distribution and velocities.

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]

The POINT-AGAPE collaboration is carrying out a search for gravitational microlensing toward M31 to reveal galactic dark matter in the form of MACHOs (Massive Astrophysical Compact Halo Objects) in the halos of the Milky Way and M31. A high-threshold analysis of 3 years of data yields 6 bright, short--duration microlensing events, which are confronted to a simulation of the observations and the analysis. The observed signal is much larger than expected from self lensing alone and we conclude, at the 95% confidence level, that at least 20% of the halo mass in the direction of M31 must be in the form of MACHOs if their average mass lies in the range 0.5-1 M$_\odot$. This lower bound drops to 8% for MACHOs with masses $\sim 0.01$ M$_\odot$. In addition, we discuss a likely binary microlensing candidate with caustic crossing. Its location, some 32' away from the centre of M31, supports our conclusion that we are detecting a MACHO signal in the direction of M31.

We report the results of our statistical study of the POINT-AGAPE CNe catalogue. The first Sloan $r'$ and $i'$ analysis of the MMRD relationship in M31 is performed using the resulting POINT-AGAPE CN catalogue. Good fits are produced to the MMRD in the two filters. However, we are unable to verify the $t_{15}$ relationship for either Sloan filter. The subsequent analysis of our automated CN detection pipeline has provided us with the most thorough knowledge of the completeness of a CN survey to-date. In addition, the large field of view of the survey has permitted us to probe the outburst rate well into the galactic disk, unlike previous CCD imaging surveys. Using this analysis we are able to probe the CN distribution of M31 and evaluate the global nova rate. Using models of the galactic surface brightness of M31, we show that the observed CN distribution consists of a separate bulge and disk population. We also show that the M31 bulge CN eruption rate per unit $r'$ flux is more than five times greater than that of the disk. We deduce a global M31 CN rate of $65^{+16}_{-15}$ year$^{-1}$, a value much higher than found by previous surveys. Using the global rate, we derive a M31 bulge rate of $38^{+15}_{-12}$ year$^{-1}$ and a disk rate of $27^{+19}_{-15}$ year$^{-1}$. Given our understanding of the completeness and an analysis of other sources of error, we conclude that the true global nova rate of M31 is at least 50% higher than was previously thought and this has consequent implications for the presumed CN rate in the Milky Way. We deduce a Galactic bulge rate of $14^{+6}_{-5}$ year$^{-1}$, a disk rate of $20^{+14}_{-11}$ year$^{-1}$ and a global Galactic rate of $34^{+15}_{-12}$ year$^{-1}$.

The POINT-AGAPE collaboration surveyed M31 with the primary goal of optical detection of microlensing events, yet its data catalogue is also a prime source of lightcurves of variable and transient objects, including classical novae (CNe). A reliable means of identification, combined with a thorough survey of the variable objects in M31, provides an excellent opportunity to locate and study an entire galactic population of CNe. This paper presents a set of 440 neural networks, working in 44 committees, designed specifically to identify fast CNe. The networks are developed using training sets consisting of simulated novae and POINT-AGAPE lightcurves, in a novel variation on K-fold cross-validation. They use the binned, normalised power spectra of the lightcurves as input units. The networks successfully identify 9 of the 13 previously identified M31 CNe within their optimal working range (and 11 out of 13 if the network error bars are taken into account). They provide a catalogue of 19 new candidate fast CNe, of which 4 are strongly favoured.

An automated search is carried out for microlensing events using a catalogue of 44554 variable superpixel lightcurves derived from our three-year monitoring program of M31. Each step of our candidate selection is objective and reproducible by a computer. Our search is unrestricted, in the sense that it has no explicit timescale cut. So, it must overcome the awkward problem of distinguishing long-timescale microlensing events from long-period stellar variables. The basis of the selection algorithm is the fitting of the superpixel lightcurves to two different theoretical models, using variable star and blended microlensing templates. Only if microlensing is preferred is an event retained as a possible candidate. Further cuts are made with regard to (i) sampling, (ii) goodness of fit of the peak to a Paczynski curve, (iii) consistency of the microlensing hypothesis with the absence of a resolved source, (iv) achromaticity, (v) position in the colour-magnitude diagram and (vi) signal-to-noise ratio. Our results are reported in terms of first-level candidates, which are the most trustworthy, and second-level candidates, which are possible microlensing but have lower signal-to-noise and are more questionable. The pipeline leaves just 3 first-level candidates, all of which have very short full-width half-maximum timescale (<5 days) and 3 second-level candidates, which have timescales of 31, 36 and 51 days respectively. We also show 16 third-level lightcurves, as an illustration of the events that just fail the threshold for designation as microlensing candidates. They are almost certainly mainly variable stars. Two of the 3 first-level candidates correspond to known events (PA 00-S3 and PA 00-S4) already reported by the POINT-AGAPE project. The remaining first-level candidate is new.

To reveal the galactic dark matter in the form of MACHOs ("Massive Astrophysical Compact Halo Objects"), the POINT-AGAPE collaboration is carrying out a search for gravitational microlensing towards M31. A clear microlensing signal is detected. The high-threshold analysis of 3-year data leads to 7 bright and short microlensing candidates. The preliminary estimation of the detection efficiency implies that less than 25% (60%) of standard halos can be composed of objects with masses between 10^{-4} and 10^{-1} (10^{-1} and 1) solar mass at the 95% C.L. This result is compatible with previous microlensing results towards the Magellanic Clouds and gives the first constraints on MACHOs for a distant spiral galaxy.

The POINT-AGAPE survey is an optical search for gravitational microlensing events towards the Andromeda Galaxy (M31). As well as microlensing, the survey is sensitive to many different classes of variable stars and transients. Here we describe the automated detection and selection pipeline used to identify M31 classical novae (CNe) and we present the resulting catalogue of 20 CN candidates observed over three seasons. CNe are observed both in the bulge region as well as over a wide area of the M31 disk. Nine of the CNe are caught during the final rise phase and all are well sampled in at least two colours. The excellent light-curve coverage has allowed us to detect and classify CNe over a wide range of speed class, from very fast to very slow. Among the light-curves is a moderately fast CN exhibiting entry into a deep transition minimum, followed by its final decline. We have also observed in detail a very slow CN which faded by only 0.01 mag day$^{-1}$ over a 150 day period. We detect other interesting variable objects, including one of the longest period and most luminous Mira variables. The CN catalogue constitutes a uniquely well-sampled and objectively-selected data set with which to study the statistical properties of classical novae in M31, such as the global nova rate, the reliability of novae as standard-candle distance indicators and the dependence of the nova population on stellar environment. The findings of this statistical study will be reported in a follow-up paper.

The POINT-AGAPE collaboration has been monitoring M31 for three seasons with the Wide Field Camera on the Isaac Newton Telescope. In each season, data are taken for one hour per night for roughly sixty nights during the six months that M31 is visible. The two fields of view straddle the central bulge, northwards and southwards. We have calculated the locations, periods and amplitudes of 35414 variable stars in M31 as a by-product of our microlensing search. The variables are classified according to their period and amplitude of variation. They are classified into population I and II Cepheids, Miras and semi-regular long-period variables. The population I Cepheids are associated with the spiral arms, while the central concentration of the Miras and long-period variables varies noticeably, the stars with brighter (and shorter) variations being much more centrally concentrated. A crucial role in the microlensing experiment is played by the asymmetry signal. It was initially assumed that the variable stars would not be a serious problem as their distributions would be symmetric. We demonstrate that this assumption is not correct. We find that differential extinction associated with the dust lanes causes the variable star distributions to be asymmetric. The size and direction of the asymmetry of the variable stars is measured as a function of period and amplitude of variation. The implications of this discovery for the successful completion of the microlensing experiments towards M31 are discussed. (Abridged)

The lightcurve of PA-99-N2, one of the recently announced microlensing candidates towards M31, shows small deviations from the standard Paczynski form. We explore a number of possible explanations, including correlations with the seeing, the parallax effect and a binary lens. We find that the observations are consistent with an unresolved RGB or AGB star in M31 being microlensed by a binary lens. We find that the best fit binary lens mass ratio is about one hundredth, which is one of most extreme values found for a binary lens so far. If both the source and lens lie in the M31 disk, then the standard M31 model predicts the probable mass range of the system to be 0.02-3.6 solar masses (95 % confidence limit). In this scenario, the mass of the secondary component is therefore likely to be below the hydrogen-burning limit. On the other hand, if a compact halo object in M31 is lensing a disk or spheroid source, then the total lens mass is likely to lie between 0.09-32 solar masses, which is consistent with the primary being a stellar remnant and the secondary a low mass star or brown dwarf. The optical depth (or alternatively the differential rate) along the line of sight toward the event indicates that a halo lens is more likely than a stellar lens provided that dark compact objects comprise no less than 15 per cent (or 5 per cent) of haloes.

The POINT-AGAPE collaboration is currently searching for massive compact halo objects (MACHOs) towards the Andromeda galaxy (M31). The survey aims to exploit the high inclination of the M31 disk, which causes an asymmetry in the spatial distribution of M31 MACHOs. Here, we investigate the effects of halo velocity anisotropy and flattening on the asymmetry signal using simple halo models. For a spherically symmetric and isotropic halo, we find that the underlying pixel-lensing rate in far-disk M31 MACHOs is more than 5 times the rate of near-disk events. We find that the asymmetry is increased further by about 30% if the MACHOs occupy radial orbits rather than tangential orbits, but is substantially reduced if the MACHOs lie in a flattened halo. However, even for haloes with a minor-to-major axis ratio q = 0.3, the numbers of M31 MACHOs in the far-side outnumber those in the near-side by a factor of ~2. We show that, if positional information is exploited in addition to number counts, then the number of candidate events required to confirm asymmetry for a range of flattened and anisotropic halo models is achievable, even with significant contamination by variable stars and foreground microlensing events. For pixel-lensing surveys which probe a representative portion of the M31 disk, a sample of around 50 candidates is likely to be sufficient to detect asymmetry within spherical haloes, even if half the sample is contaminated, or to detect asymmetry in haloes as flat as q = 0.3 provided less than a third of the sample comprises contaminants. We also argue that, provided its mass-to-light ratio is less than 100, the recently observed stellar stream around M31 is not problematic for the detection of asymmetry. (Abstract slightly abridged.)

We report the final analysis of a search for microlensing events in the direction of the Andromeda galaxy, which aimed to probe the MACHO composition of the M31 halo using data collected during the 1998-99 observational campaign at the MDM observatory. In a previous paper, we discussed the results from a first set of observations. Here, we deal with the complete data set, and we take advantage of some INT observations in the 1999-2000 seasons. This merging of data sets taken by different instruments turns out to be very useful, the study of the longer baseline available allowing us to test the uniqueness characteristic of microlensing events. As a result, all the candidate microlensing events previously reported turn out to be variable stars. We further discuss a selection based on different criteria, aimed at the detection of short--duration events. We find three candidates whose positions are consistent with self--lensing events, although the available data do not allow us to conclude unambiguously that they are due to microlensing.

We have carried out a survey of the Andromeda galaxy for unresolved microlensing (pixel lensing). We present a subset of four short timescale, high signal-to-noise microlensing candidates found by imposing severe selection criteria: the source flux variation exceeds the flux of an R=21 magnitude star and the full width at half maximum timescale is less than 25 days. Remarkably, in three out of four cases, we have been able to measure or strongly constrain the Einstein crossing time of the event. One event, which lies projected on the M31 bulge, is almost certainly due to a stellar lens in the bulge of M31. The other three candidates can be explained either by stars in M31 and M32 or by MACHOs.

We report the discovery of a microlensing candidate projected 2'54" from the center of M32, on the side closest to M31. The blue color (R-I= 0.00 +/- 0.14) of the source argues strongly that it lies in the disk of M31, while the proximity of the line of sight to M32 implies that this galaxy is the most likely host of the lens. If this interpretation is correct, it would confirm previous arguments that M32 lies in front of M31. We estimate that of order one such event or less should be present in the POINT-AGAPE data base. If more events are discovered in this direction in a dedicated experiment, they could be used to measure the mass function of M32 up to an unknown scale factor. By combining microlensing observations of a binary-lens event with a measurement of the M31-M32 relative proper motion using the astrometric satellites SIM or GAIA, it will be possible to measure the physical separation of M31 and M32, the last of the six phase-space coordinates needed to assign M32 an orbit.

We report the discovery of a short-duration microlensing candidate in the Northern field of the POINT-AGAPE pixel lensing survey towards M31. The full-width at half-maximum timescale is very short (tfwhm = 1.8 days) and the source star has almost certainly been identified on HST images, allowing to infer an Einstein crossing time : te=10.4 days and a magnification at maximum : Amax=18.

We report the discovery of a short-duration microlensing candidate in the northern field of the POINT-AGAPE pixel lensing survey towards M31. The full-width half-maximum timescale is very short, just 1.8 days. Almost certainly, the source star has been identified on Hubble Space Telescope archival images, allowing us to infer an Einstein crossing time of 10.4 days, a maximum magnification of about 18, and a lens-source proper motion greater than 0.3 microarcsec/day. The event lies projected at 8' from the center of M31, which is beyond the bulk of the stellar lens population. The lens is likely to reside in one of three locations. It may be a star in the M31 disk, or a massive compact halo object (Macho) in either M31 or the Milky Way. The most probable mass is 0.06 solar masses for an M31 Macho, 0.02 solar masses for a Milky Way Macho and 0.2 solar masses for an M31 stellar lens. Whilst the stellar interpretation is plausible, the Macho interpretation is the most probable for halo fractions above 20%.