A&A 400, 951-956 (2003)
DOI: 10.1051/0004-6361:20030087
C. Afonso^{1,4,8} - J. N. Albert^{2} - J. Andersen^{6} - R. Ansari^{2} - É. Aubourg^{1} - P. Bareyre^{1,4} - J. P. Beaulieu^{5} - G. Blanc^{1} - X. Charlot^{1} - F. Couchot^{2} - C. Coutures^{1} - R. Ferlet^{5} - P. Fouqué^{9,10} - J. F. Glicenstein^{1} - B. Goldman^{1,4,8} - A. Gould^{1,7} - D. Graff^{7} - M. Gros^{1} - J. Haissinski^{2} - C. Hamadache^{1} - J. de Kat^{1} - T. Lasserre^{1} - L. Le Guillou^{1} - É. Lesquoy^{1,5} - C. Loup^{5} - C. Magneville ^{1} - J. B. Marquette^{5} - É. Maurice^{3} - A. Maury^{9} - A. Milsztajn ^{1} - M. Moniez^{2} - N. Palanque-Delabrouille^{1} - O. Perdereau^{2} - L. Prévot^{3} - Y. R. Rahal^{2} - J. Rich^{1} - M. Spiro^{1} - P. Tisserand^{1} - A. Vidal-Madjar^{5} - L. Vigroux^{1} - S. Zylberajch^{1}
1 - DSM/DAPNIA, CEA/Saclay, 91191 Gif-sur-Yvette Cedex, France
2 -
Laboratoire de l'Accélérateur Linéaire,
IN2P3 CNRS, Université Paris-Sud, 91405 Orsay Cedex, France
3 -
Observatoire de Marseille,
2 Pl. Le Verrier, 13248 Marseille Cedex 04, France
4 -
Collège de France, Physique Corpusculaire et Cosmologie, IN2P3 CNRS,
11 Pl. M. Berthelot, 75231 Paris Cedex, France
5 -
Institut d'Astrophysique de Paris, INSU CNRS,
98bis boulevard Arago, 75014 Paris, France
6 -
Astronomical Observatory, Copenhagen University, Juliane Maries Vej 30,
2100 Copenhagen, Denmark
7 -
Departments of Astronomy and Physics, Ohio State University, Columbus,
OH 43210, USA
8 -
Department of Astronomy, New Mexico State University, Las Cruces, NM 88003-8001, USA
9 -
European Southern Observatory (ESO), Casilla 19001, Santiago 19, Chile
10 -
Observatoire de Paris, LESIA, 92195 Meudon Cedex, France
Received 16 December 2002 / Accepted 17 January 2003
Abstract
Five years of EROS data towards the Small Magellanic Cloud have been
searched for gravitational microlensing events, using a new, more
accurate method to assess the impact of stellar blending on the
efficiency. Four long-duration candidates have been found which, if
they are microlensing events, hint at a non-halo population of lenses.
Combined with results from other
EROS observation programs, this analysis yields strong limits on the amount of
Galactic dark matter made of compact objects. Less than 25% of a
standard halo can be composed of objects with a mass between
and 1
at the 95% C.L.
Key words: Galaxy: halo - Galaxy: kinematics and dynamics - Galaxy: stellar content - cosmology: dark matter - gravitational lensing
The idea of using gravitational microlensing as a tool to probe the dark matter in the Galactic halo goes back to 1986 with a paper by B. Paczynski (Paczynski 1986). Since then, several experiments have been monitoring millions of stars towards the Large and Small Magellanic Clouds ( LMC and SMC) and candidates have been observed towards the two targets (Alcock et al. 1993; Aubourg et al. 1993; Ansari et al. 1996; Alcock et al. 1997a; Afonso et al. 1999; Lasserre et al. 2000).
The detection of candidates towards the LMC suggested a population of lenses accounting for 20% of a standard halo, with a most probable mass between 0.2 and 0.9 (Alcock et al. 2000). Results by Lasserre et al. (2000) show that objects up to 1 cannot account for more than 40% of a standard halo. These lenses cannot be ordinary Galactic stars, whose density is much too low to account for the observed optical depth, and Goldman et al. (2002) demonstrated they cannot be halo white dwarfs with a hydrogen atmosphere. Much debate has occurred on the nature and the location of the lenses, the main two possibilities being 0.4 dark objects in the halo or stars at the low mass end of the main sequence in the clouds themselves (Sahu 1994; Wu 1994).
This makes the SMC a very valuable target. While the characteristics (optical depth and duration) of events towards the two clouds should be similar for halo lenses (Sackett & Gould 1993), the different dynamical properties of the clouds could account for differences if the events are due to self-lensing.
The EROS 2 experiment has been surveying the SMC since 1996. We present here an analysis of five years of data accumulated towards this target.
The telescope, camera, telescope operations and data reduction are as described in Palanque-Delabrouille et al. (1998), hereafter SM98, and references therein. Since July 1996, 10 one-square-degree fields have been monitored towards the SMC. The first 5 years of data from 8.6 square degrees spread over these 10 fields have been analyzed. The data set contains about 5.2 million pairs of light curves in two wide pass bands called thereafter "red'' and "blue'', sampled on average once every 2.5 days from end-April to mid-March when the SMC is visible. About 400-500 images of each field were taken, with exposure times ranging from 5 min in the center of the SMC to 15 min in the outermost region.
The analysis of the five-year data set is similar to those of the first year ( SM98) and of the first two years (Afonso et al. 1999) of data where details can be found.
A first set of cuts is designed to select light curves that exhibit a single significant fluctuation. This already excludes 99.8% of the stars, most of which exhibit flat light curves.
A second set of cuts rejects stars that are a priori likely to be variable, based upon their position in a color-magnitude diagram: bright blue stars in the upper main sequence with low-amplitude variations and stars much brighter and redder than those of the red clump are rejected. These stars are contained in two thinly populated regions of the color-magnitude diagram, so rejecting them does not reduce much the number of light curves on which we can search for microlensing. Stars exhibiting a strong correlation between the red and the blue light curves outside the period containing the main fluctuation are also rejected. Supernovae are rejected by the comparison of the rising and falling times of the variation. Half of the light curves surviving the first set of cuts are rejected at this stage, leaving about 5000 light curves.
A third set of cuts improves the signal-to-noise ratio of the set of selected candidates by comparing the measurements with the best-fit point-source point-lens microlensing light curve. These criteria are sufficiently loose not to reject light curves affected by blending, parallax and most cases of binary lenses or sources. The remaining sample consists of about 30 light curves exhibiting a clean and unique fluctuation with a smooth time variation.
A fourth set of cuts constrains the time coverage of the event, excluding in particular all events with Einstein radius crossing times exceeding 1200 days.
The final criterion is as follows: if blending significantly improves the fit, then the fitted blending should be physical, i.e. indeed correspond to the amplification of only a fraction of the total flux recovered.
To limit contamination of the set of selected candidates by spurious events, the levels of the above cuts are set such that most variable stars are rejected by at least two of them. The tuning of each cut and the estimation of the efficiency of the analysis is done with Monte Carlo simulated light curves, as described in SM98. Microlensing parameters are drawn uniformly in the following intervals: time of maximum magnification t_{0} within the observing period 300 days, impact parameter normalized to the Einstein radius and time-scale days.
Four candidates ( SMC-1 to SMC-4) pass all the cuts. Their light curves are shown in Fig. 1, their coordinates are given in Table 1 and the microlensing fit parameters (point-source, point-lens and no-blending) in Table 2.
Figure 1: Red light curves of the microlensing candidates detected towards the SMC. | |
Open with DEXTER |
V_{J} | I_{C} | |||
SMC-1 | 01:00:05.73 | -72:15:02.33 | 17.9 | 17.9 |
SMC-2 | 00:48:20.16 | -74:12:35.31 | 20.4 | 19.5 |
SMC-3 | 00:49:38.68 | -74:14:06.77 | 20.2 | 19.3 |
SMC-4 | 01:05:36.56 | -72:15:02.39 | 19.9 | 19.4 |
u_{0} | |||
SMC-1 | 0.52 | 101 | 1026/706 |
SMC-2 | 0.82 | 390 | 1526/990 |
SMC-3 | 0.66 | 612 | 2316/968 |
SMC-4 | 0.80 | 243 | 2923/951 |
5 | 15 | 50 | 125 | 325 | 500 | 900 | 1100 | 1300 | |
2.1 | 6.3 | 11.1 | 13.8 | 11.8 | 10.1 | 7.5 | 5.9 | 2.7 |
Candidate SMC-1 is the one already present in the previous analyses
(one year and two years of SMC data) done by EROS. Further analysis
of its light curve indicated that it was probably a self-lensing event
( SM98). It was also detected as a possible microlensing
candidate by the online trigger system of the MACHO experiment
(Alcock et al. 1997b). Candidates SMC-2, SMC-3 and SMC-4 have very long
durations, with few points outside the amplified region of the light curve. The
of candidate SMC-4 is large, yet far
from the cut level (which is set at 5.0 in the peak region). All three candidates, however, exhibit features in their baseline that are reminiscent of variable stars. These features appear clearly on the light curves plotted in 25-day bins (see Fig. 2). To allow a direct comparison of the red and blue fluxes, all fluxes have been normalized so that they would have a mean of zero and a variance of unity. Candidates SMC-2 and SMC-4 clearly look like irregular variables, while SMC-3 has an irregular light curve and seems to be rising at the end.
These candidates will need to be monitored for several more years to be either
confirmed but most probably ruled out as microlensing candidates. They are not in the 3
square degrees of the SMC on which the MACHO project took
data.
Figure 2: Red and blue light curves of the microlensing candidates detected towards the SMC, in 25-day bins. | |
Open with DEXTER |
The efficiency of the analysis for events with (Monte-Carlo) u_{0}<1.5 normalized to an observing period of five years is summarized in Table 3.
The above efficiencies have been derived without taking account of stellar blending.
We have performed a study to quantify precisely the influence of blending on our results. We simulate an experiment with synthetic images, regularly spaced in time, in which microlensing events are embedded. The images are built from the HST colour-magnitude diagram (CMD) and luminosity function (LF) of the Magellanic Clouds, which extend well beyond the EROS detection limit (down to V=23.5 in the HST data). We have checked that the CMD and LF that we construct from these simulated images are in good agreement with those obtained from our data. All stars are lensed for a timescale of 5 images. The impact parameter u_{0} is drawn uniformly between 0 and 1.5, and t_{0} is set so that no two stars less than 20 pixels apart are magnified simultaneously (typical EROS images have a seeing of 2.1 arcsec corresponding to a Gaussian with rms of 1.5 pixel).
The standard photometry chain of EROS is applied to these
images, leading to the detection of 5599 objects (in this
terminology an object usually encompasses several simulated stars) with their red and blue light curves. A total of 6304
generated microlensing events were found with u_{0}<1 and an average
timescale
images. The effect of
blending is readily seen from the fact that more events than objects
are detected. Of these events, 60% are due to the brightest
star in the two pixels around the object, with recovered
images and
,
close to
the generated values.
The remaining events are due to a fainter, blended component of the
object, with recovered
images, well
underestimated, and
,
clearly overestimated,
as expected from the impact of a significant blend. The
optical depth being proportional to the number of events and to their
mean ,
an estimate of the ratio R of the recovered
optical depth to the generated one is:
The light curves from the above blending simulation are then used as templates in our Monte Carlo chain designed to compute the efficiency of our analysis (Sect. 3) as follows: for each object detected in the actual data, an object is randomly chosen from the blending simulation in the same region of the CMD. The light curve of the latter is shifted and stretched to account for the randomly selected duration and time of maximum amplification, but the amplitudes in each colour are of course kept as found in the blending simulation to reflect the impact of the underlying stellar companions of the microlensed star. This template is used (instead of the standard shape derived by Paczynski in the absence of blending and used for the efficiency computation in the previous section) to add a microlensing event to the real (flat in most cases) light curve. To ensure similar photometric dispersion in the Monte Carlo and in the data, the rms flux deviation of each simulated data point from the blended microlensing light curve is taken to be the same as the deviation of the real light curve from constant brightness. The modified light curve goes through the same analysis cuts as in the standard computation of the efficiency.
The efficiency with the blending simulation described above is within 10% of the efficiency computed from the standard procedure, in agreement with estimates in previous works (Palanque-Delabrouille 1997). The loss due to the under-estimation of the magnification of a blended event (the observed magnification on a blended object is lower than the physical magnification of the underlying star) is almost exactly compensated by the fact that an object includes more than one star subject to lensing. In addition, as expected, the cuts of the analysis are sufficiently loose not to reject blended microlensing events. In particular, there are no cuts requiring a strict achromaticity of the event. For the EROS experiment, blending has a small impact on the efficiency and can therefore be neglected.
In order to set limits on the contribution of dark objects to the halo, we use the so-called "standard'' halo model described in SM98 as model 1, and take into account the efficiency of the analysis given in Sect. 3.
For a given experiment, assuming that the halo is made of compact
objects having a single mass M, we define
as the
probability distribution of expected event durations, taking efficiencies into
account, and
as the expected number of
events, where f is the halo mass fraction with respect to a full standard
halo. We construct a frequentist confidence level 1-p_{n} considering
only the number n of observed events as a Poisson process:
A well known prescription for the confidence level of the combined test is
We then set a frequentist 95% CL limit on f, taking into account both the number of candidates and their duration, by finding the value of f that yields .
Although this prescription could be used to combine several experiments, it gives the same weight to all and does not yield useful results if the sensitivities of some of the experiments differ significantly. Therefore, to combine various experiments, we use instead a Liptak-Stouffer (Liptak 1958) prescription where the flat p-values are converted to unit-variance, zero-mean Gaussian distributed variables. A weighted sum of these variables is also Gaussian distributed, and can be converted back to a p-value. We take as weights the expected number of events of each experiment for each lens mass^{}.
The results obtained by the various phases of EROS (which are independent experiments) are summarized in Table 2 for microlensing candidates towards the SMC - no SMC candidate in EROS1 (Renault et al. 1998), 4 candidates in EROS2 ( SM98and present work) - and in table 4 for microlensing candidates towards the LMC - no planetary mass candidate (Renault et al. 1997, 1998), LMC-1 from EROS1 photographic plates (Ansari et al. 1996) and the others from EROS2 (Lasserre et al. 2000; Lasserre 2000; Milsztajn et al. 2001).
u_{0} | ||
LMC-1 | 0.44 | 23 |
LMC-3 | 0.21 | 44 |
LMC-5 | 0.59 | 24 |
LMC-6 | 0.41 | 35 |
LMC-7 | 0.30 | 30 |
Figure 3 shows the 95% exclusion limit derived from this method on f, the halo mass fraction, at any given mass M - i.e. assuming all deflectors in the halo have mass M - for the EROS1 CCD LMC and EROS1 CCD SMC, EROS1 photographic plates, EROS2 3-year LMC and EROS2 5-year SMC experiments, and for the combination of all. We also show in the figure the limits that would be obtained with a single SMC event ( SMC-1), then with no event at all in any of the EROS experiments, which indicates the overall sensibility of the EROS project, considering presently analysed data.
Figure 3: Exclusion diagram at 95% C.L. for the standard halo model ( inside 50 kpc). The dashed lines are the limits towards the LMC by EROS 1 and EROS 2, the thin plain lines are limits towards the SMC and the thick line is the combined limit from the five EROS sub-experiments. The dotted lines are the limits that would be obtained considering no observed events: they indicate the overall sensibility of the EROS dataset. We also indicate the limit that would be obtained if the three very long duration events on the SMC were considered as background (label " SMC-1 only''): the limit is lowered in a negligible way in the 2-3 mass range. | |
Open with DEXTER |
The "dent'' in the EROS1 plate limit and in the EROS2 LMC limit at a mass near is the impact of the 30 day candidates observed towards the LMC. For any mass between and , we exclude at 95% C.L. that more than 20% of the mass of a standard halo be made of compact objects. It can be seen that the combined limit is above the best limit for some values of the mass. This occurs quite naturally since the observed SMC and LMC candidates have quite different characteristics. At high mass, for instance, our method will consider SMC candidates as signal and LMC candidates as background, thus weakening the limit obtained by the LMC alone. It illustrates the (marginally significant) incompatibility between the candidates observed by EROS towards the SMC and the LMC.
Another way to gauge the present results on microlensing towards the SMC is to compare them directly to the result of the MACHO collaboration towards the LMC (Alcock et al. 2000). Such a comparison is made easier by the fact that the expected distributions of halo microlenses towards the SMC and LMC are very similar: their averages and widths differ by only a few percent, and these differences are much smaller than the widths themselves. The ratio of the optical depths towards the SMC and LMC is more uncertain; for a spherical isothermal halo, it is expected to be about , but this value can be lower for flattened halos. Actually, it was proposed to use observations towards the SMC to evaluate the galactic dark halo flattening (Sackett & Gould 1993).
In their analysis A, the MACHO group has listed 13 microlens
candidates with an average duration
days and a width
of the
distribution
compatible with all microlenses having the same mass. From this,
they compute an optical depth
,
of which about
is attributed to halo lenses. This translates into an
expected number of events for the present EROS SMC data set
of
There are, however, no EROS SMC microlensing candidates that are compatible with the distribution observed by the MACHO group towards the LMC: candidate SMC-1, the shortest of all, is very likely due to a lens in the SMC (see parallax analysis and the discussion that follows in SM98); the other 3 candidates all have durations longer than 240 days, clearly incompatible with an average of 36 days. From this absence of microlensing candidates with the required duration, EROS can exclude for events similar to those in MACHO sample A, at better than 90% C.L. The value corresponding to the spherical isothermal halo model, is excluded at better than 96% C.L.
The analysis of five years of microlensing data towards the SMC yielded four candidates, and allowed us to put stringent limits on the amount of galactic dark matter made of compact objects. Objects with a mass between and 1 cannot account for more than 25% of the mass of a standard spherical, isothermal and isotropic Galactic halo of out to 50 kpc.
The new method described in this paper to combine the results and use all the information available on the events has one drawback compared to methods used in previous works: it could be more sensitive to the halo model, adding an extra dependence on the velocity distribution. Previously, only the influence of the halo model on the number of expected events has been studied, as in SM98. A change in the velocity distribution, however, would mainly shift the "dent'' in the EROS2 LMC limit. In addition, since the final exclusion limit is quite flat over a large range of masses, the final result should not be influenced by reasonable changes in the velocity distribution assumptions.
This limit excludes a sizeable fraction of the allowed domain shown in Fig. 12 of Alcock et al. (2000).
All the candidates that have been detected so far towards the SMC have long durations, and seem more compatible with unidentified variable stars or self-lensing within the cloud than with halo objects. These would need to have supersolar masses to account for such durations.
The statistics are still very low for a quantitative comparison of the SMC and LMC samples of microlensing candidates. They have, however, quite different characteristics which is poorly compatible with a unique population of lenses.
Acknowledgements
We thank J. Bouchez for very useful discussions. We are grateful to D. Lacroix and the technical staff at the Observatoire de Haute Provence and to A. Baranne for their help in refurbishing the MARLY telescope and remounting it in La Silla. We are also grateful for the support given to our project by the technical staff at ESO, La Silla. We thank J. F. Lecointe for assistance with the online computing. AG was supported by grant AST 02-01266 from the NSF.