Direction dependent calibration and imaging is a vital part of producing deep, high fidelity, high-dynamic range radio images with a wide-field low-frequency array like LOFAR. Currently, state-of-the-art facet-based direction dependent calibration algorithms rely on the assumption that the isoplanatic-patch size is much larger than the separation between bright in-field calibrators. This assumption is often violated due to the dynamic nature of the ionosphere, and as a result direction dependent errors affect image quality between calibrators. In this paper we propose a probabilistic physics-informed model for inferring ionospheric phase screens, providing a calibration for all sources in the field of view. We apply our method to a randomly selected observation from the LOFAR Two-Metre Sky Survey archive, and show that almost all direction dependent effects between bright calibrators are corrected and that the root-mean-squared residuals around bright sources is reduced by 32\% on average.

The largest uncertainty on measurements of dark energy using type Ia supernovae is presently due to systematics from photometry; specifically to the relative uncertainty on photometry as a function of wavelength in the optical spectrum. We show that a precise constraint on relative photometry between the visible and near-infrared can be achieved at upcoming survey telescopes (such as at the Vera Rubin Observatory [VRO]) via a mountaintop-located laser source tuned to the 342.78 nm vacuum excitation wavelength of neutral sodium atoms. Using a high-power (500 W) laser modified from laser guide star studies, this excitation will produce an artificial star (which we term a "laser photometric ratio star," or LPRS) of de-excitation light in the mesosphere that is observable from the ground at approximately 20 magnitude (i.e., well within the expected single-image magnitude limit of VRO) at wavelengths in vacuum of 589.16 nm, 589.76 nm, 818.55 nm, and 819.70 nm, with the sum of the numbers of 589.16 nm and 589.76 nm photons produced by this process equal to the sum of the numbers of 818.55 nm and 819.70 nm photons, establishing a precise calibration ratio between, for example, the VRO r and z filters. This technique can thus provide a novel mechanism for establishing a spectrophotometric calibration ratio of unprecedented precision, from above most of the Earth's atmosphere, for upcoming telescopic observations across astronomy and atmospheric physics. This article is the first in a pair of articles on this topic. The second article of the pair describes an alternative technique to achieve a similar, but brighter, LPRS than the technique described in this paper, by using two mountaintop-located lasers, at optical frequencies approximately 4 GHz away from resonances at wavelengths in vacuum of 589.16 nm and 819.71 nm respectively.

Calibrating for direction-dependent ionospheric distortions in visibility data is one of the main technical challenges that must be overcome to advance low-frequency radio astronomy. In this paper, we propose a novel probabilistic, tomographic approach that utilises Gaussian processes to calibrate direction-dependent ionospheric phase distortions in low-frequency interferometric data. We suggest that the ionospheric free electron density can be modelled to good approximation by a Gaussian process restricted to a thick single layer, and show that under this assumption the differential total electron content must also be a Gaussian process. We perform a comparison with a number of other widely successful Gaussian processes on simulated differential total electron contents over a wide range of experimental conditions, and find that, in all experimental conditions, our model is better able to represent observed data and generalise to unseen data. The mean equivalent source shift imposed by our predictive errors are half as large as the best competitor model's. We find that it is possible to partially constrain the ionosphere's hyperparameters from sparse-and-noisy observed data. Our model provides an alternative explanation for observed phase structure functions deviating from Kolmogorov's 5/3 turbulence, turnover at high baselines, and diffractive scale anisotropy. We show that our model implicitly cheaply performs tomography of the free electron density. Moreover, we find that even a fast, low-resolution approximation to our model yields better results than the best alternative Gaussian process, implying that the geometric coupling between directions and antennae is a powerful prior that should not be ignored.

Hypervelocity stars (HVSs) are amongst the fastest objects in our Milky Way. These stars are predicted to come from the Galactic center (GC) and travel along unbound orbits across the Galaxy. In the coming years, the ESA satellite Gaia will provide the most complete and accurate catalogue of the Milky Way, with full astrometric parameters for more than $1$ billion stars. In this paper, we present the expected sample size and properties (mass, magnitude, spatial, velocity distributions) of HVSs in the Gaia stellar catalogue. We build three Gaia mock catalogues of HVSs anchored to current observations, exploring different ejection mechanisms and GC stellar population properties. In all cases, we predict hundreds to thousands of HVSs with precise proper motion measurements within a few tens of kpc from us. For stars with a relative error in total proper motion below $10 \%$, the mass range extends to ~$10 M_{\odot}$ but peaks at ~$1$ $M_\odot$. The majority of Gaia HVSs will therefore probe a different mass and distance range compared to the current non-Gaia sample. In addition, a subset of a few hundreds to a few thousands of HVSs with $M$ ~ $3$ $M_\odot$ will be bright enough to have a precise measurement of the three-dimensional velocity from Gaia alone. Finally, we show that Gaia will provide more precise proper motion measurements for the current sample of HVS candidates. This will help identifying their birthplace narrowing down their ejection location, and confirming or rejecting their nature as HVSs. Overall, our forecasts are extremely encouraging in terms of quantity and quality of HVS data that can be exploited to constrain both the Milky Way potential and the GC properties.

We present radio observations of the galaxy cluster PLCK G004.5-19.5 ($z=0.52$) using the Giant Metrewave Radio Telescope at 150~MHz, 325~MHz, and 610~MHz. We find an unusual arrangement of diffuse radio emission in the center and periphery of the cluster, as well as several radio galaxies with head-tail emission. A patch of peripheral emission resembles a radio relic, and central emission resembles a radio halo. Reanalysis of archival XMM-Newton X-ray data shows that PLCK G004.5-19.5 is disturbed, which has a known correlation with the existence of radio relics and halos. Given that the number of known radio halos and radio relics at $z>0.5$ is very limited, PLCK G004.5-19.5 is an important addition to understanding merger-related particle acceleration at higher redshifts.

Procedures and results on hardware level detector calibration in Super-Kamiokande (SK) are presented in this paper. In particular, we report improvements made in our calibration methods for the experimental phase IV in which new readout electronics have been operating since 2008. The topics are separated into two parts. The first part describes the determination of constants needed to interpret the digitized output of our electronics so that we can obtain physical numbers such as photon counts and their arrival times for each photomultiplier tube (PMT). In this context, we developed an in-situ procedure to determine high-voltage settings for PMTs in large detectors like SK, as well as a new method for measuring PMT quantum efficiency and gain in such a detector. The second part describes the modeling of the detector in our Monte Carlo simulation, including in particular the optical properties of its water target and their variability over time. Detailed studies on the water quality are also presented. As a result of this work, we achieved a precision sufficient for physics analysis over a wide energy range (from a few MeV to above a TeV). For example, the charge determination was understood at the 1% level, and the timing resolution was 2.1 nsec at the one-photoelectron charge level and 0.5 nsec at the 100-photoelectron charge level.

The GMRT Epoch of Reionization (EoR) experiment is an ongoing effort to measure the power spectrum from neutral hydrogen at high redshift. We have previously reported an upper limit of (70 mK)^2 at wavenumbers of k=0.65 h/Mpc using a basic piecewise-linear foreground subtraction. In this paper we explore the use of a singular value decomposition to remove foregrounds with fewer assumptions about the foreground structure. Using this method we also quantify, for the first time, the signal loss due to the foreground filter and present new power spectra adjusted for this loss, providing a revised measurement of a 2-sigma upper limit at (248 mK)^2 for k=0.50 h/Mpc. While this revised limit is larger than previously reported, we believe it to be more robust and still represents the best current constraints on reionization at z=8.6.

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

Deeper understanding of the properties of dark energy via SNIa surveys, and to a large extent other methods as well, will require unprecedented photometric precision. Laboratory and solar photometry and radiometry regularly achieve precisions on the order of parts in ten thousand, but photometric calibration for non-solar astronomy presently remains stuck at the percent or greater level. We discuss our project to erase this discrepancy, and our steps toward achieving laboratory-level photometric precision for surveys late this decade. In particular, we show near-field observations of the balloon-borne light source we are presently testing, in addition to previous work with a calibrated laser source presently in low-Earth orbit. Our technique is additionally applicable to microwave astronomy. Observation of gravitational waves in the polarized CMB will similarly require unprecedented polarimetric and radiometric precision, and we briefly discuss our plans for a calibrated microwave source above the atmosphere as well.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton decays via the Callan-Rubakov process. In this scenario, protons, which initially decay into pions, will ultimately produce \nu_{e}, \nu_{\mu} and \bar{\nu}_{\mu}. After undergoing neutrino oscillation, all neutrino species appear when they arrive at the Earth, and can be detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). A search for low energy neutrinos in the electron total energy range from 19 to 55 MeV was carried out with SK and gives a monopole flux limit of F_M(\sigma_0/1 mb) < 6.3 \times 10^{-24} (\beta_M/10^{-3})^2 cm^{-2} s^{-1} sr^{-1} at 90% C.L., where \beta_M is the monopole velocity in units of the speed of light and \sigma_0 is the catalysis cross section at \beta_M=1. The obtained limit is more than eight orders of magnitude more stringent than the current best cosmic-ray supermassive monopole flux limit, F_M < 1 \times 10^{-15} cm^{-2} s^{-1} sr^{-1} for \beta_M < 10^{-3} and also two orders of magnitude lower than the result of the Kamiokande experiment, which used a similar detection method.

A significant and growing portion of systematic error on a number of fundamental parameters in astrophysics and cosmology is due to uncertainties from absolute photometric and flux standards. A path toward achieving major reduction in such uncertainties may be provided by satellite-mounted light sources, resulting in improvement in the ability to precisely characterize atmospheric extinction, and thus helping to usher in the coming generation of precision results in astronomy. Using a campaign of observations of the 532 nm pulsed laser aboard the CALIPSO satellite, collected using a portable network of cameras and photodiodes, we obtain initial measurements of atmospheric extinction, which can apparently be greatly improved by further data of this type. For a future satellite-mounted precision light source, a high-altitude balloon platform under development (together with colleagues) can provide testing as well as observational data for calibration of atmospheric uncertainties.

A significant and growing portion of systematic error on a number of fundamental parameters in astrophysics and cosmology is due to uncertainties from absolute photometric and flux standards. A path toward achieving major reduction in such uncertainties may be provided by satellite-mounted light sources, resulting in improvement in the ability to precisely characterize atmospheric extinction, and thus helping to usher in the coming generation of precision results in astronomy. Toward this end, we have performed a campaign of observations of the 532 nm pulsed laser aboard the CALIPSO satellite, using a portable network of cameras and photodiodes, to precisely measure atmospheric extinction.

R Coronae Borealis stars (RCB) are a rare type of evolved carbon-rich supergiant stars that are increasingly thought to result from the merger of two white dwarfs, called the Double degenerate scenario. This scenario is also studied as a source, at higher mass, of type Ia Supernovae (SnIa) explosions. Therefore a better understanding of RCBs composition would help to constrain simulations of such events. We searched for and studied RCB stars in the EROS Magellanic Clouds database. We also extended our research to DY Per type stars (DYPers) that are expected to be cooler RCBs (T~3500 K) and much more numerous than their hotter counterparts. The light curves of ~70 millions stars have been analysed to search for the main signature of RCBs and DYPers: a large drop in luminosity. Follow-up optical spectroscopy was used to confirm each photometric candidate found. We have discovered and confirmed 6 new Magellanic Cloud RCB stars and 7 new DYPers, but also listed new candidates: 3 RCBs and 14 DYPers. We estimated a range of Magellanic RCB shell temperatures between 360 and 600 K. We confirm the wide range of absolute luminosity known for RCB stars, M_V~-5.2 to -2.6. Our study further shows that mid-infrared surveys are ideal to search for RCB stars, since they have thinner and cooler circumstellar shells than classical post-AGB stars. In addition, by increasing the number of known DYPers by ~400%, we have been able to shed light on the similarities in the spectral energy distribution between DYPers and ordinary carbon stars. We also observed that DYPer circumstellar shells are fainter and hotter than those of RCBs. This suggests that DYPers may simply be ordinary carbon stars with ejection events, but more abundance analysis is necessary to give a status on a possible evolutionnary connexion between RCBs and DYPers.

The EROS-2 project has been designed to search for microlensing events towards any dense stellar field. The densest parts of the Galactic spiral arms have been monitored to maximize the microlensing signal expected from the stars of the Galactic disk and bulge. 12.9 million stars have been monitored during 7 seasons towards 4 directions in the Galactic plane, away from the Galactic center. A total of 27 microlensing event candidates have been found. Estimates of the optical depths from the 22 best events are provided. A first order interpretation shows that simple Galactic models with a standard disk and an elongated bulge are in agreement with our observations. We find that the average microlensing optical depth towards the complete EROS-cataloged stars of the spiral arms is $\bar{\tau} =0.51\pm .13\times 10^{-6}$, a number that is stable when the selection criteria are moderately varied. As the EROS catalog is almost complete up to $I_C=18.5$, the optical depth estimated for the sub-sample of bright target stars with $I_C<18.5$ ($\bar{\tau}=0.39\pm >.11\times 10^{-6}$) is easier to interpret. The set of microlensing events that we have observed is consistent with a simple Galactic model. A more precise interpretation would require either a better knowledge of the distance distribution of the target stars, or a simulation based on a Galactic model. For this purpose, we define and discuss the concept of optical depth for a given catalog or for a limiting magnitude.

A number of microlensing dark-matter surveys have produced tens of millions of light curves of individual background stars. These data provide an unprecedented opportunity for systematic studies of whole classes of variable stars and their host galaxies. We aim to use the EROS-2 survey of the Magellanic Clouds to detect and study the population of beat Cepheids (BCs) in both Clouds. BCs pulsating simultaneously in the first overtone and fundamental modes (FO/F) or in the second and first overtone modes (SO/FO) are of particular interest. Using special software designed to search for periodic variables, we have scanned the EROS-2 data base for variables in the typical period range of Cepheids. Metallicities of FO/F objects were then calculated from linear nonadiabatic convective stellar models. We identify 74 FO/F BCs in the LMC and 41 in the SMC, and 173 and 129 SO/FO pulsators in the LMC and SMC, respectively; 185 of these stars are new discoveries. For nearly all the FO/F objects we determine minimum, mean, and maximum values of the metallicity. The EROS data have expanded the samples of known BCs in the LMC by 31%, in the SMC by 110%. The FO/F objects provide independent measures of metallicities in these galaxies. The mean value of metallicity is 0.0045 in the LMC and 0.0018 in the SMC.

M87 is the only known non blazar radio galaxy to emit very high energy (VHE) gamma-rays. During a monitoring program of M87, a rapid flare in VHE gamma-rays was detected by the MAGIC telescope in early 2008. The flux was found to be variable above 350 GeV on a timescale as short as 1 day at a significance level of $5.6\sigma$. The highest measured flux reached 15% of the Crab Nebula flux. We observed several substantial changes of the flux level during the 13 day observing period. The flux at lower energies (150 -- 350 GeV), instead, is compatible with being constant. The energy spectrum can be described by a power law with a photon index of $2.30 \pm 0.11_\mathrm{stat} \pm 0.20_\mathrm{syst}$. The observed day-scale flux variability at VHE prefers the M87 core as source of the emission and implies that either the emission region is very compact (just a few Schwarzschild radii) or the Doppler factor of the emitting blob is rather large in the case of a non expanding emission region.

The atmospheric Cherenkov gamma-ray telescope MAGIC, designed for a low-energy threshold, has detected very-high-energy gamma rays from a giant flare of the distant Quasi-Stellar Radio Source (in short: radio quasar) 3C 279, at a distance of more than 5 billion light-years (a redshift of 0.536). No quasar has been observed previously in very-high-energy gamma radiation, and this is also the most distant object detected emitting gamma rays above 50 gigaelectron volts. Since high-energy gamma rays may be stopped by interacting with the diffuse background light in the universe, the observations by MAGIC imply a low amount for such light, consistent with that known from galaxy counts.

We report about very high energy (VHE) gamma-ray observations of the Crab Nebula with the MAGIC telescope. The gamma-ray flux from the nebula was measured between 60 GeV and 9 TeV. The energy spectrum can be described with a curved power law dF/dE=f0 (E/300 GeV)^(a+b*log10(E/300 GeV)) with a flux normalization f0 of(6.0+-0.2)*10^-10 1/(cm^2 s TeV), a=-2.31+-0.06 and b=-0.26+-0.07. The position of the IC-peak is determined at (77+-47) GeV. Within the observation time and the experimental resolution of the telescope, the gamma-ray emission is steady and pointlike. The emission's center of gravity coincides with the position of the pulsar. Pulsed gamma-ray emission from the pulsar could not be detected. We constrain the cutoff energy of the spectrum to be less than 27 GeV, assuming that the differential energy spectrum has an exponential cutoff. For a super-exponential shape, the cutoff energy can be as high as 60 GeV.

Until April 2007 the MAGIC telescope used a 300 MSamples/s FADC system to sample the shaped PMT signals produced by the captured Cherenkov photons of air showers. Different algorithms to reconstruct the signal from the read-out samples (extractors) have been implemented and are described and compared. Criteria based on the obtained charge and time resolution/bias are defined and used to judge the different extractors, by applying them to calibration, cosmic and pedestal signals. The achievable charge and time resolution have been derived as functions of the incident number of photo-electrons.

The MAGIC collaboration has recently reported the discovery of gamma-ray emission from the binary system LS I +61 303 in the TeV energy region. Here we present new observational results on this source in the energy range between 300 GeV and 3 TeV. In total 112 hours of data were taken between September and December 2006 covering 4 orbital cycles of this object. This large amount of data allowed us to produce an integral flux light curve covering for the first time all orbital phases of LS I +61 303. In addition, we also obtained a differential energy spectrum for two orbital phase bins covering the phase range 0.5<phi<0.6 and 0.6<phi<0.7. The photon index in the two phase bins is consistent within the errors with an average index Gamma=2.6+-0.2_{stat}+-0.2_{sys}. LS I +61 303 was found to be variable at TeV energies on timescales of days. These new MAGIC measurements allowed us to search for intra-night variability of the VHE emission; however, no evidence for flux variability on timescales down to 30 minutes was found. To test for possible periodic structures in the light curve, we apply the formalism developed by Lomb and Scargle to the LS I +61 303 data taken in 2005 and 2006. We found the LS I +61 303 data set to be periodic with a period of (26.8+-0.2) days (with a post-trial chance probability of 10^{-7}), close to the orbital period.

For seven years, the EROS-2 project obtained a mass of photometric data on variable stars. We present a peculiar Cepheid-like star, in the direction of the Small Magellanic Cloud, which demonstrates unusual photometric behaviour over a short time interval. We report on data of the photometry acquired by the MARLY telescope and spectroscopy from the EFOSC instrument for this star, called EROS2 J005135-714459(sm0060n13842), which resembles the unusual Cepheid HR 7308. The light curve of our target is analysed using the Analysis of Variance method to determine a pulsational period of 5.5675 days. A fit of time-dependent Fourier coefficients is performed and a search for proper motion is conducted. The light curve exhibits a previously unobserved and spectacular change in both mean magnitude and amplitude, which has no clear theoretical explanation. Our analysis of the spectrum implies a radial velocity of 104 km s$^{-1}$ and a metallicity of -0.4$\pm$0.2 dex. In the direction of right ascension, we measure a proper motion of 17.4$\pm$6.0 mas yr$^{-1}$ using EROS astrometry, which is compatible with data from the NOMAD catalogue. The nature of EROS2 J005135-714459(sm0060n13842) remains unclear. For this star, we may have detected a non-zero proper motion for this star, which would imply that it is a foreground object. Its radial velocity, pulsational characteristics, and photometric data, however, suggest that it is instead a Cepheid-like object located in the SMC. In such a case, it would present a challenge to conventional Cepheid models.

We present the results of the first multiwavelength observing campaign on the high-mass X-ray binary LS I +61 303 comprising observations at the TeV regime with the MAGIC telescope, along with X-ray observations with Chandra, and radio interferometric observations with the MERLIN, EVN and VLBA arrays, in October and November 2006. From our MERLIN observations, we can exclude the existence of large scale (~100 mas) persistent radio-jets. Our 5.0 GHz VLBA observations display morphological similarities to previous 8.4 GHz VLBA observations carried out at the same orbital phase, suggesting a high level of periodicity and stability of the processes behind the radio emission. This makes it unlikely that variability of the radio emission is due to the interaction of an outflow with variable wind clumps. If the radio emission is produced by a milliarcsecond scale jet, it should also show a stable, periodic behavior. It is then difficult to reconcile the absence of a large scale jet (~100 mas) in our observations with the evidence of a persistent relativistic jet reported previously. We find a possible hint of temporal correlation between the X-ray and TeV emissions and evidence for radio/TeV non-correlation, which points to the existence of one population of particles producing the radio emission and a different one producing the X-ray and TeV emissions. Finally, we present a quasi-simultaneous energy spectrum including radio, X-ray and TeV bands.

We observed the first known very high energy (VHE) gamma-ray emitting unidentified source, TeV J2032+4130, for 94 hours with the MAGIC telescope. The source was detected with a significance of 5.6 sigma. The flux, position, and angular extension are compatible with the previous ones measured by the HEGRA telescope system five years ago. The integral flux amounts to (4.5+-0.3stat+-0.35sys)x10^{-13} ph cm$^{-2}$ s$^{-1}$ above 1 TeV. The source energy spectrum, obtained with the lowest energy threshold to date, is compatible with a single power law with a hard photon index of Gamma=-2.0+-0.3stat+-0.2sys.

We present the results from a multiwavelength campaign on the TeV blazar 1ES 1959+650, performed in May, 2006. Data from the optical, UV, soft- and hard-X-ray and very high energy (VHE) gamma-ray (E > 100 GeV) bands were obtained with the SUZAKU and SWIFT satellites, with the MAGIC telescope and other ground based facilities. The source spectral energy distribution (SED), derived from SUZAKU and MAGIC observations at the end of May 2006, shows the usual double hump shape, with the synchrotron peak at a higher flux level than the Compton peak. With respect to historical values, during our campaign the source exhibited a relatively high state in X-rays and optical, while in the VHE band it was at one of the lowest level so far recorded. We also monitored the source for flux-spectral variability on a time window of 10 days in the optical-UV and X-ray bands and 7 days in the VHE band. The source varies more in the X-ray, than in the optical band, with the 2-10 keV X-ray flux varying by a factor of ~2. The synchrotron peak is located in the X-ray band and moves to higher energies as the source gets brighter, with the X-ray fluxes above it varying more rapidly than the X-ray fluxes at lower energies. The variability behaviour observed in the X-ray band cannot be produced by emitting regions varying independently, and suggests instead some sort of ``standing shock'' scenario. The overall SED is well represented by an homogeneous one-zone synchrotron inverse Compton emission model, from which we derive physical parameters that are typical of high energy peaked blazars.

The nearby dwarf spheroidal galaxy Draco with its high mass to light ratio is one of the most auspicious targets for indirect dark matter searches. Annihilation of hypothetical DM particles can result in high-energy gamma-rays, e.g. from neutralino annihilation in the supersymmetric framework. With the MAGIC telescope a search for a possible DM signal originating from Draco was performed during 2007. The analysis of the data results in a flux upper limit of 1.1x10^-11 photons cm^-2 sec^-1 for photon energies above 140 GeV, assuming a point like source. Furthermore, a comparison with predictions from supersymmetric models is given. While our results do not constrain the mSUGRA phase parameter space, a very high flux enhancement can be ruled out.

Rare types of variable star may give unique insight into short-lived stages of stellar evolution. The systematic monitoring of millions of stars and advanced light curve analysis techniques of microlensing surveys make them ideal for discovering also such rare variable stars. One example is the R Coronae Borealis (RCB) stars, a rare type of evolved carbon-rich supergiant. We have conducted a systematic search of the EROS-2 database for the Galactic catalogue Bulge and spiral arms to find Galactic RCB stars. The light curves of $\sim$100 million stars, monitored for 6.7 years (from July 1996 to February 2003), have been analysed to search for the main signature of RCB stars, large and rapid drops in luminosity. Follow-up spectroscopy has been used to confirm the photometric candidates. We have discovered 14 new RCB stars, all in the direction of the Galactic Bulge, bringing the total number of confirmed Galactic RCB stars to about 51. After reddening correction, the colours and absolute magnitudes of at least 9 of the stars are similar to those of Magellanic RCB stars. This suggests that these stars are in fact located in the Galactic Bulge, making them the first RCB stars discovered in the Bulge. The localisation of the 5 remaining RCBs is more uncertain: 4 are either located behind the Bulge at an estimated maximum distance of 14 kpc or have an unusual thick circumstellar shell; the other is a DY Per RCB which may be located in the Bulge, even if it is fainter than the known Magellanic DY Per. From the small scale height found using the 9 new Bulge RCBs, $61<h^{RCB}_{Bulge}<246$ pc (95% C.L.), we conclude that the RCB stars follow a disk-like distribution inside the Bulge.

The MAGIC Project: Contributions to ICRC 2007, Merida, Mexico. Contents pages for the Contribution on behalf of the MAGIC Collaboration to the 30th ICRC that took place in July 2007 in Merida, Mexico. The contents are in html form with clickable links to the papers that exist on the Astrophysics archive. We hope that this will make it easier to access the output of the conference in a systematic way. Comments on how useful this is/ how it could be improved should be sent to michela.demaria@iuav.it.

We report on very high energy gamma-observations with the MAGIC Telescope of the pulsar PSR B1951+32 and its associated nebula, CTB 80. Our data constrain the cutoff energy of the pulsar to be less than 32 GeV, assuming the pulsed gamma-ray emission to be exponentially cut off. The upper limit on the flux of pulsed gamma-ray emission above 75 GeV is 4.3*10^-11 photons cm^-2 sec^-1, and the upper limit on the flux of steady emission above 140 GeV is 1.5*10^-11 photons cm^-2 sec^-1. We discuss our results in the framework of recent model predictions and other studies.

We searched for very high energy (VHE) gamma-ray emission from the supernova remnant Cassiopeia A. The shell-type supernova remnant Cassiopeia A was observed with the 17 meter MAGIC telescope between July 2006 and January 2007 for a total time of 47 hours. The source was detected above an energy of 250 GeV with a significance of 5.2sigma and a photon flux above 1 TeV of (7.3+-0.7_stat+-2.2_sys)10^-13 cm^-2s^-1. The photon spectrum is compatible with a power law dN/dE = A E^-Gamma with a photon index Gamma=2.3+-0.2_stat+-0.2_sys. The source is point-like within the angular resolution of the telescope.

We report on the discovery of Very High Energy (VHE) gamma-ray emission from the BL Lacertae object 1ES1011+496. The observation was triggered by an optical outburst in March 2007 and the source was observed with the MAGIC telescope from March to May 2007. Observing for 18.7 hr we find an excess of 6.2 sigma with an integrated flux above 200 GeV of (1.58$\pm0.32) 10^{-11}$ photons cm$^{-2}$ s$^{-1}$. The VHE gamma-ray flux is >40% higher than in March-April 2006 (reported elsewhere), indicating that the VHE emission state may be related to the optical emission state. We have also determined the redshift of 1ES1011+496 based on an optical spectrum that reveals the absorption lines of the host galaxy. The redshift of z=0.212 makes 1ES1011+496 the most distant source observed to emit VHE gamma-rays up to date.

The paper describes an application of the tree classification method Random Forest (RF), as used in the analysis of data from the ground-based gamma telescope MAGIC. In such telescopes, cosmic gamma-rays are observed and have to be discriminated against a dominating background of hadronic cosmic-ray particles. We describe the application of RF for this gamma/hadron separation. The RF method often shows superior performance in comparison with traditional semi-empirical techniques. Critical issues of the method and its implementation are discussed. An application of the RF method for estimation of a continuous parameter from related variables, rather than discrete classes, is also discussed.

The MAGIC collaboration observed BL Lacertae for 22.2 hr during 2005 August to December and for 26 hr during 2006 July to September. The source is the historical prototype and eponym of a class of low-frequency peaked BL Lacertae (LBL) objects. A very high energy (VHE) gamma-ray signal was discovered with a 5.1 sigma excess in the 2005 data. Above 200 GeV, an integral flux of 0.6+-0.2 x 10^(-11) cm^(-2) s^(-1) was measured, corresponding to approximately 3% of the Crab flux. The differential spectrum between 150 and 900 GeV is rather steep with a photon index of -3.6+-0.5. The light curve shows no significant variability during the observations in 2005. For the first time a clear detection of VHE gamma-ray emission from an LBL object was obtained with a signal below previous upper limits. The 2006 data show no significant excess. This drop in flux follows the observed trend in optical activity.

The paper describes the different methods, used in the MAGIC experiment, to unfold experimental energy distributions of cosmic ray particles (gamma-rays). Questions and problems related to the unfolding are discussed. Various procedures are proposed which can help to make the unfolding robust and reliable. The different methods and procedures are implemented in the MAGIC software and are used in most of the analyses.

We report on the results from the observations in very high energy band (VHE, E_gamma > 100 GeV) of the black hole X-ray binary (BHXB) Cygnus X-1. The observations were performed with the MAGIC telescope, for a total of 40 hours during 26 nights, spanning the period between June and November 2006. Searches for steady gamma-ray signals yielded no positive result and upper limits to the integral flux ranging between 1 and 2% of the Crab nebula flux, depending on the energy, have been established. We also analyzed each observation night independently, obtaining evidence of gamma-ray signals at the 4.0 standard deviations (sigma) significance level (3.2 sigma after trial correction) for 154 minutes effective on-time (EOT) on September 24 between 20h58 and 23h41 UTC, coinciding with an X-ray flare seen by RXTE, Swift and INTEGRAL. A search for faster-varying signals within a night resulted in an excess with a significance of 4.9 sigma (4.1 sigma after trial correction) for 79 minutes EOT between 22h17 and 23h41 UTC. The measured excess is compatible with a point-like source at the position of Cygnus X-1, and excludes the nearby radio nebula powered by its relativistic jet. The differential energy spectrum is well fitted by an unbroken power-law described by dN/(dA dt dE) = (2.3+- 0.6) x 10^{-12} (E/1 TeV)^{-3.2 +- 0.6}. This is the first experimental evidence of VHE emission from a stellar-mass black hole, and therefore from a confirmed accreting X-ray binary.

The fast repositioning system of the MAGIC Telescope has allowed during its first data cycle, between 2005 and the beginning of year 2006, observing nine different GRBs as possible sources of very high energy gammas. These observations were triggered by alerts from Swift, HETE-II, and Integral; they started as fast as possible after the alerts and lasted for several minutes, with an energy threshold varying between 80 and 200 GeV, depending upon the zenith angle of the burst. No evidence for gamma signals was found, and upper limits for the flux were derived for all events, using the standard analysis chain of MAGIC. For the bursts with measured redshift, the upper limits are compatible with a power law extrapolation, when the intrinsic fluxes are evaluated taking into account the attenuation due to the scattering in the Metagalactic Radiation Field (MRF).

The EROS-2 project was designed to test the hypothesis that massive compact halo objects (the so-called ``machos'') could be a major component of the dark matter halo of the Milky Way galaxy. To this end, EROS-2 monitored over 6.7 years $33\times10^6$ stars in the Magellanic clouds for microlensing events caused by such objects. In this work, we use only a subsample of $7\times10^6$ bright stars spread over $84 \deg^2$ of the LMC and $9 \deg^2$ of the SMC. The strategy of using only bright stars helps to discriminate against background events due to variable stars and allows a simple determination of the effects of source confusion (blending). The use of a large solid angle makes the survey relatively insensitive to effects that could make the optical depth strongly direction dependent. Using this sample of bright stars, only one candidate event was found, whereas $\sim39$ events would have been expected if the Halo were entirely populated by objects of mass $M\sim0.4M_{\odot}$. Combined with the results of EROS-1, this implies that the optical depth toward the Large Magellanic Cloud (\object{LMC}) due to such lenses is $\tau<0.36\times10^{-7}$ (95%CL), corresponding to a fraction of the halo mass of less than 8%. This optical depth is considerably less than that measured by the MACHO collaboration in the central region of the LMC. More generally, machos in the mass range $0.6\times10^{-7}M_\odot<M<15M_{\odot}$ are ruled out as the primary occupants of the Milky Way Halo.

We study the capability of the MAGIC telescope to observe under moderate moonlight. TeV gamma-ray signals from the Crab nebula were detected with the MAGIC telescope during periods when the Moon was above the horizon and during twilight. This was accomplished by increasing the trigger discriminator thresholds. No change is necessary in the high voltage settings since the camera PMTs were especially designed to avoid high currents. We characterize the telescope performance by studying the effect of the moonlight on the gamma-ray detection efficiency and sensitivity, as well as on the energy threshold.

(abridged) We examine the spatial and temporal stability of the HST ACS Wide Field Camera (WFC) point spread function (PSF) using the two square degree COSMOS survey. We show that stochastic aliasing of the PSF necessarily occurs during `drizzling'. This aliasing is maximal if the output pixel scale is equal to the input pixel scale of 0.05''. We show that this source of PSF variation can be significantly reduced by choosing a Gaussian drizzle kernel and by setting the output pixel size to 0.03''. We show that the PSF is temporally unstable, most likely due to thermal fluctuations in the telescope's focus. We find that the primary manifestation of this thermal drift in COSMOS images is an overall slow periodic focus change. Using a modified version of TinyTim, we create undistorted stars in a 30x30 grid across the ACS WFC CCDs. These PSF models are created for telescope focus values in the range -10 microns to +5 microns, thus spanning the allowed range of telescope focus values. We then use the approximately ten well measured stars in each COSMOS field to pick the best-fit focus value for each field. The TinyTim model stars are then used to perform PSF corrections for weak lensing allowing systematics due to incorrectly modeled PSFs to be greatly reduced. We have made the software for PSF modeling using our modified version of TinyTim available to the astronomical community. We show the effects of Charge Transfer Efficiency (CTE) degradation, which distorts the object in the readout direction, mimicking a weak lensing signal. We derive a parametric correction for the effect of CTE on the shapes of objects in the COSMOS field as a function of observation date, position within the ACS WFC field, and object flux.

We present a three dimensional cosmic shear analysis of the Hubble Space Telescope COSMOS survey, the largest ever optical imaging program performed in space. We have measured the shapes of galaxies for the tell-tale distortions caused by weak gravitational lensing, and traced the growth of that signal as a function of redshift. Using both 2D and 3D analyses, we measure cosmological parameters Omega_m, the density of matter in the universe, and sigma_8, the normalization of the matter power spectrum. The introduction of redshift information tightens the constraints by a factor of three, and also reduces the relative sampling (or "cosmic") variance compared to recent surveys that may be larger but are only two dimensional. From the 3D analysis, we find sigma_8*(Omega_m/0.3)^0.44=0.866+^0.085_-0.068 at 68% confidence limits, including both statistical and potential systematic sources of error in the total budget. Indeed, the absolute calibration of shear measurement methods is now the dominant source of uncertainty. Assuming instead a baseline cosmology to fix the geometry of the universe, we have measured the growth of structure on both linear and non-linear physical scales. Our results thus demonstrate a proof of concept for tomographic analysis techniques that have been proposed for future weak lensing surveys by a dedicated wide-field telescope in space.

We present a new EROS-2 measurement of the microlensing optical depth toward the Galactic Bulge. Light curves of $5.6\times 10^{6}$ clump-giant stars distributed over $66 \deg^2$ of the Bulge were monitored during seven Bulge seasons. 120 events were found with apparent amplifications greater than 1.6 and Einstein radius crossing times in the range $5 {\rm d}<t_\e <400 {\rm d}$. This is the largest existing sample of clump-giant events and the first to include northern Galactic fields. In the Galactic latitude range $1.4\degr<|b|<7.0\degr$, we find $\tau/10^{-6}=(1.62 \pm 0.23)\exp[-a(|b|-3 {\rm deg})]$ with $a=(0.43 \pm0.16)\deg^{-1}$. These results are in good agreement with our previous measurement, with recent measurements of the MACHO and OGLE-II groups, and with predictions of Bulge models.

Microquasars are binary star systems with relativistic radio-emitting jets. They are potential sources of cosmic rays and laboratories for elucidating the physics of relativistic jets. Here we report the detection of variable gamma-ray emission above 100 gigaelectron volts from the microquasar LS I +61 303. Six orbital cycles were recorded. Several detections occur at a similar orbital phase, suggesting the emission is periodic. The strongest gamma-ray emission is not observed when the two stars are closest to one another, implying a strong orbital modulation of the emission or the absorption processes.

We propose a tunable laser-based satellite-mounted spectrophotometric and absolute flux calibration system, to be utilized by ground- and space-based telescopes. As spectrophotometric calibration may play a significant role in the accuracy of photometric redshift measurement, and photometric redshift accuracy is important for measuring dark energy using SNIa, weak gravitational lensing, and baryon oscillations, a method for reducing such uncertainties is needed. We propose to improve spectrophotometric calibration, currently obtained using standard stars, by placing a tunable laser and a wide-angle light source on a satellite by early next decade (perhaps included in the upgrade to the GPS satellite network) to improve absolute flux calibration and relative spectrophotometric calibration across the visible and near-infrared spectrum. As well as fundamental astrophysical applications, the system proposed here potentially has broad utility for defense and national security applications such as ground target illumination and space communication.

The MAGIC collaboration has studied the high peaked BL-Lac object 1ES1218+30.4 at a redshift z = 0.182, using the MAGIC imaging air Cherenkov telescope located on the Canary island of La Palma. A gamma-ray signal was observed with 6.4sigma significance. The differential energy spectrum for an energy threshold of 120GeV can be fitted by a simple power law yielding F_E(E) = (8.1+-2.1)*10^-7 (E/250GeV)^(-3.0+-0.4) TeV^-1 m^-2 s^-1. During the six days of observation in January 2005 no time variability on time scales of days was found within the statistical errors. The observed integral flux above 350GeV is nearly a factor two below the the upper limit reported by the Whipple Collaboration in 2003.

The ability to accurately measure the shapes of faint objects in images taken with the Advanced Camera for Surveys(ACS) on the Hubble Space Telescope (HST) depends upon detailed knowledge of the Point Spread Function (PSF). We show that thermal fluctuations cause the PSF of the ACS Wide Field Camera (WFC) to vary over time. We describe a modified version of the TinyTim PSF modeling software to create artificial grids of stars across the ACS field of view at a range of telescope focus values. These models closely resemble the stars in real ACS images. Using ~10 bright stars in a real image, we have been able to measure HST's apparent focus at the time of the exposure. TinyTim can then be used to model the PSF at any position on the ACS field of view. This obviates the need for images of dense stellar fields at different focus values, or interpolation between the few observed stars. We show that residual differences between our TinyTim models and real data are likely due to the effects of Charge Transfer Efficiency (CTE) degradation. Furthermore, we discuss stochastic noise that is added to the shape of point sources when distortion is removed, and we present Multidrizzle parameters that are optimal for weak lensing science. Specifically, we find that reducing the Multidrizzle output pixel scale and choosing a Gaussian kernel significantly stabilizes the resulting PSF after image combination, while still eliminating cosmic rays/bad pixels, and correcting the large geometric distortion in the ACS. We discuss future plans, which include more detailed study of the effects of CTE degradation on object shapes and releasing our TinyTim models to the astronomical community.

The long-duration GRB050713a was observed by the MAGIC Telescope, 40 seconds after the burst onset, and followed up for 37 minutes, until twilight. The observation, triggered by a SWIFT alert, covered energies above ~175 GeV. Using standard MAGIC analysis, no evidence for a gamma signal was found. As the redshift of the GRB was not measured directly, the flux upper limit, estimated by MAGIC, is still compatible with the assumption of an unbroken power-law spectrum extending from a few hundred keV to our energy range.

Recently, the Galactic Center has been reported to be a source of very high energy (VHE) gamma-rays by the VERITAS, CANGAROO and HESS experiments. The energy spectra as measured by these experiments show substantial differences. In this Letter we present MAGIC observations of the Galactic Center, resulting in the detection of a differential gamma-ray flux consistent with a steady, hard-slope power law, described as dN/(dA dt dE) = (2.9+/-0.6)*10^{-12} (E/TeV)^{-2.2+/-0.2} cm^{-2}s^{-1}TeV^{-1}. The gamma-ray source is centered at (Ra, Dec)=(17h 45m 20s, -29deg 2'. This result confirms the previous measurements by the HESS experiment and indicates a steady source of TeV gamma-rays. We briefly describe the observational technique used, the procedure implemented for the data analysis, and discuss the results in the perspective of different models proposed for the acceleration of the VHE gamma-rays.

Recently, the HESS collaboration has reported the detection of gamma-ray emission above a few hundred GeV from eight new sources located close to the Galactic Plane. The source HESS J1813-178 has sparked particular interest, as subsequent radio observations imply an association with SNR G12.82-0.02. Triggered by the detection in VHE gamma-rays, a positionally coincident source has also been found in INTEGRAL and ASCA data. In this Letter we present MAGIC observations of HESS J1813-178, resulting in the detection of a differential gamma-ray flux consistent with a hard-slope power law, described as dN/(dA dt dE) = (3.3+/-0.5)*10^{-12} (E/TeV)^{-2.1+/-0.2} cm^(-2)s^(-1)TeV^(-1). We briefly discuss the observational technique used, the procedure implemented for the data analysis, and put this detection in the perspective of multifrequency observations.

We present the type Ia rate measurement based on two EROS supernova search campaigns (in 1999 and 2000). Sixteen supernovae identified as type Ia were discovered. The measurement of the detection efficiency, using a Monte Carlo simulation, provides the type Ia supernova explosion rate at a redshift ~ 0.13. The result is $0.125^{+0.044+0.028}_{-0.034-0.028} h_{70}^2$ SNu where 1 SNu = 1 SN / $10^{10} L_{sun}^B$ / century. This value is compatible with the previous EROS measurement (Hardin et al. 2000), done with a much smaller sample, at a similar redshift. Comparison with other values at different redshifts suggests an evolution of the type Ia supernova rate.

We study the accuracy with which weak lensing measurements could be made from a future space-based survey, predicting the subsequent precisions of 3-dimensional dark matter maps, projected 2-dimensional dark matter maps, and mass-selected cluster catalogues. As a baseline, we use the instrumental specifications of the Supernova/Acceleration Probe (SNAP) satellite. We first compute its sensitivity to weak lensing shear as a function of survey depth. Our predictions are based on detailed image simulations created using `shapelets', a complete and orthogonal parameterization of galaxy morphologies. We incorporate a realistic redshift distribution of source galaxies, and calculate the average precision of photometric redshift recovery using the SNAP filter set to be Delta z=0.034. The high density of background galaxies resolved in a wide space-based survey allows projected dark matter maps with a rms sensitivity of 3% shear in 1 square arcminute cells. This will be further improved using a proposed deep space-based survey, which will be able to detect isolated clusters using a 3D lensing inversion techniques with a 1 sigma mass sensitivity of approximately 10^13 solar masses at z~0.25. Weak lensing measurements from space will thus be able to capture non-Gaussian features arising from gravitational instability and map out dark matter in the universe with unprecedented resolution.

Weak gravitational lensing provides a unique method to directly map the dark matter in the universe and measure cosmological parameters. Current weak lensing surveys are limited by the atmospheric seeing from the ground and by the small field of view of existing space telescopes. We study how a future wide-field space telescope can measure the lensing power spectrum and skewness, and set constraints on cosmological parameters. The lensing sensitivity was calculated using detailed image simulations and instrumental specifications studied in earlier papers in this series. For instance, the planned SuperNova/Acceleration Probe (SNAP) mission will be able to measure the matter density parameter Omega_m and the dark energy equation of state parameter w with precisions comparable and nearly orthogonal to those derived with SNAP from supernovae. The constraints degrade by a factor of about 2 if redshift tomography is not used, but are little affected if the skewness only is dropped. We also study how the constraints on these parameters depend upon the survey geometry and define an optimal observing strategy.

We present a measurement of the microlensing optical depth toward the Galactic bulge based on the analysis of 15 contiguous1 square degrees fields centered on (l=2.5 deg, b=-4.0 deg) and containing 1.42 million clump-giant stars (belonging to the extended clump area) monitored during almost three bulge seasons by EROS (Experience de Recherche d'Objets Sombres). We find a microlensing optical depth towards the bulge tau_bulge=0.94 +/- 0.29 10^-6 averaged over all fields, based on 16 microlensing events with clump giants as sources. This value is substantially below several other determinations by the MACHO and OGLE groups and is more in agreement with what is expected from axisymmetric and non-axisymmetric bulge models.

We are conducting a 377-square-degree proper motion survey in the ~V and I bands in order to determine the cool white dwarf contribution to the Galactic dark matter. Using the 250 square degrees for which we possess three epochs, and applying selection criteria designed to isolate halo-type objects, we find no candidates in a 5500 pc^3 effective volume for old, fast M_V=17 white dwarfs. We check the detection efficiency by cross-matching our catalogue with Luyten's NLTT catalogue. The halo white dwarf contribution cannot exceed 5% (95% C.L.) for objects with M_V=17 and 1<V-I<1.5. The same conclusion applies to a 14Gyr halo composed of white dwarfs with hydrogen atmosphere, as modeled by Chabrier (99).

We present the results of a massive variability search based on a photometric survey of a six square degree region along the Galactic plane at ($l = 305^\circ$, $b = -0.8^\circ$) and ($l = 330^\circ$, $b = -2.5^\circ$). This survey was performed in the framework of the EROS II (Exp\'erience de Recherche d'Objets Sombres) microlensing program. The variable stars were found among 1,913,576 stars that were monitored between April and June 1998 in two passbands, with an average of 60 measurements. A new period-search technique is proposed which makes use of a statistical variable that characterizes the overall regularity of the flux versus phase diagram. This method is well suited when the photometric data are unevenly distributed in time, as is our case. 1,362 objects whose luminosity varies were selected. Among them we identified 9 Cepheids, 19 RR Lyrae, 34 Miras, 176 eclipsing binaries and 266 Semi-Regular stars. Most of them are newly identified objects. The cross-identification with known catalogues has been performed. The mean distance of the RR Lyrae is estimated to be $\sim 4.9 \pm 0.3$ kpc undergoing an average absorption of $\sim 3.4 \pm 0.2$ magnitudes. This distance is in good agreement with the one of disc stars which contribute to the microlensing source star population.Our catalogue and light curves are available electronically from the CDS, Strasbourg and from our Web site http://eros.in2p3.fr.

We present the EROS nearby supernova ($z \sim 0.02 - 0.2$) search and the analysis of the first year of data (1997). A total of 80 square degrees were surveyed. Eight supernov{\ae} were detected, four of which were spectroscopically identified as type Ia supernov{\ae}. The search efficiency was determined with a Monte-Carlo simulation taking into account the efficiencies for both supernova detection and host galaxy identification. Assuming that for a given galaxy the supernova rate is proportional to the galactic luminosity, we compute a type Ia supernova explosion rate of: ${\cal R} = 0.44 {}_{-0.21}^{+0.35} {}_{-0.07}^{+0.13} h^2: / 10^{10} \lbsun / 100 {\rm yrs}$ at an average redshift of $\sim 0.1$ where the errors are respectively statistical and systematic (type misidentification included).

We present the analysis of the light curves of 8.5 million stars observed during two seasons by EROS (Experience de Recherche d'Objets Sombres), in the galactic plane away from the bulge. Three stars have been found that exhibit luminosity variations compatible with gravitational microlensing effects due to unseen objects. The corresponding optical depth, averaged over four directions, is 0.38 (+0.53, -0.15) 10^{-6}. All three candidates have long Einstein radius crossing times ($\sim$ 70 to 100 days). For one of them, the lack of evidence for a parallax or a source size effect enabled us to constrain the lens-source % geometric configuration. Another candidate displays a modulation of the magnification, which is compatible with the lensing of a binary source. The interpretation of the optical depths inferred from these observations is hindered by the imperfect knowledge of the distance to the target stars. Our measurements are compatible with expectations from simple galactic models under reasonable assumptions on the target distances.