We present the analysis of high-resolution images of MOA-2013-BLG-220, taken with the Keck adaptive optics system 6 years after the initial observation, and identify the lens as a solar type star hosting a super-Jupiter mass planet. The masses of planets and host stars discovered by microlensing are often not determined from the light curve data, while the star-planet mass ratio and projected separation in units of Einstein ring radius are well measured. High resolution follow-up observations after the lensing event is complete, can resolve the source and lens, as well as measure their fluxes, amplitude and direction of relative proper motion, giving us strong constraints on the observed system parameters. Due to the high relative proper motion of this event, $\mu_{rel,G} = 12.6 \pm 0.3$ mas/yr, we were able to resolve the source and lens with a separation of $77.6 \pm 0.4$ mas. As a consequence we constrain the lens NIR flux to $K_L = 17.92\pm 0.05$. By combining constraints from the initial model and the AO flux, we find the lens star to have a mass of $M_L = 0.89 \pm 0.05\ \rm M_\odot$ located at $D_L = 6.96 \pm 0.60$ kpc. With a mass ratio of $q = (3.26 \pm 0.04) \times 10^{-3}$ from reanalysis of the light-curve, the planet is found to have a mass of $M_P = 3.06 \pm 0.18\ M_{J}$ and a separation of $r_\perp = 2.46 \pm 0.23$ AU. This mass is much higher than the prediction of a Bayesian analysis that assumes that all stars have an equal probability to host a planet of the measured mass ratio, and it suggests that planets with a mass ratio of a few $\times 10^{-3}$ are more common orbiting massive stars. This demonstrates the importance of high angular resolution follow-up observations.

We propose a new semi-Lagrangian Vlasov-Poisson solver. It employs elements of metric to follow locally the flow and its deformation, allowing one to find quickly and accurately the initial phase-space position $Q(P)$ of any test particle $P$, by expanding at second order the geometry of the motion in the vicinity of the closest element. It is thus possible to reconstruct accurately the phase-space distribution function at any time $t$ and position $P$ by proper interpolation of initial conditions, following Liouville theorem. When distorsion of the elements of metric becomes too large, it is necessary to create new initial conditions along with isotropic elements and repeat the procedure again until next resampling. To speed up the process, interpolation of the phase-space distribution is performed at second order during the transport phase, while third order splines are used at the moments of remapping. We also show how to compute accurately the region of influence of each element of metric with the proper percolation scheme. The algorithm is tested here in the framework of one-dimensional gravitational dynamics but is implemented in such a way that it can be extended easily to four or six-dimensional phase-space. It can also be trivially generalised to plasmas.

The cosmic horseshoe gravitational lens is analyzed using the perturbative approach. The two first order perturbative fields are expanded in Fourier series. The source is reconstructed using a fine adaptive grid. The expansion of the fields at order 2 produces a higher value of the chi-square. Expanding at order 3 provides a very significant improvement, while order 4 does not bring a significant improvement over order 3. The presence of the order 3 terms is not a consequence of limiting the perturbative expansion to the first order. The amplitude and signs of the third order terms are recovered by including the contribution of the other group members. This analysis demonstrates that the fine details of the potential of the lens could be recovered independently of any assumptions by using the perturbative approach.

The extension of the singular perturbative approach to the second order is presented in this paper. The general expansion to the second order is derived. The second order expansion is considered as a small correction to the first order expansion. Using this approach it is demonstrated that the second order expansion is reducible to a first order expansion via a re-definition of the first order pertubative fields. Even if in practice the second order correction is small the reducibility of the second order expansion to the first order expansion indicates a degeneracy problem. In general this degeneracy is hard to break. A useful and simple second order approximation is the thin source approximation which offers a direct estimation of the correction. The practical application of the corrections derived in this paper are illustrated by using an elliptical NFW lens model. The second order pertubative expansion provides a noticeable improvement, even for the simplest case of thin source approximation. To conclude it is clear that for accurate modelisation of gravitational lenses using the perturbative method the second order perturbative expansion should be considered. In particular an evaluation of the degeneracy due to the second order term should be performed, for which the thin source approximation is particularly useful.

We report periods and JHKL observations for 648 oxygen-rich Mira variables found in two outer bulge fields at b=-7 degrees and l=+/-8 degrees and combine these with data on 8057 inner bulge Miras from the OGLE, Macho and 2MASS surveys, which are concentrated closer to the Galactic centre. Distance moduli are estimated for all these stars. Evidence is given showing that the bulge structure is a function of age. The longer period Miras (log P > 2.6, age about 5 Gyr and younger) show clear evidence of a bar structure inclined to the line of sight in both the inner and outer regions. The distribution of the shorter period (metal-rich globular cluster age) Miras, appears spheroidal in the outer bulge. In the inner region these old stars are also distributed differently from the younger ones and possibly suggest a more complex structure. These data suggest a distance to the Galactic centre, R0, of 8.9 kpc with an estimated uncertainty of 0.4 kpc. The possible effect of helium enrichment on our conclusions is discussed.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

The cosmological simulations indicates that the dark matter haloes have specific self similar properties. However the halo similarity is affected by the baryonic feedback, the momentum injected by the supernovae re-shape the dark matter core and transform it to a flat density core, with a scale length imposed by the baryonic feedback. Additionally the baryon feedback impose also an equilibrium condition, which when coupled with the imposed baryonic scale length induce a new type of similarity. The new self similar solution implies that the acceleration generated by dark matter is scale free, which in turns implies that the baryonic acceleration at a reference radius is also scale free. Constant dark matter and baryonic accelerations at a reference radius have effectively been observed for a large class of different galaxies, which is in support of this approach. The new self similar properties implies that the total acceleration at larger distances is scale free, the transition between the dark matter and baryons dominated regime occurs at a constant acceleration, and the maximum of the velocity curve which defines the amplitude of the velocity curve at larger distances is proportional to $M^{\frac{1}{4}}$. These results demonstrates that in this self similar model, cold dark matter is consistent with the basics of MOND phenomenology for the galaxies. In agreement with the observation the coincidence between the self similar model and MOND is expected to break at the scale of clusters of galaxies. Some numerical experiments shows that the behavior of the density near the origin is closely approximated by a Einasto profile.

This papers explores the self similar solutions of the Vlasov-Poisson system and their relation to the gravitational collapse of dynamically cold systems. Analytic solutions are derived for power law potential in one dimension, and extensions of these solutions in three dimensions are proposed. Next the self similarity of the collapse of cold dynamical systems is investigated numerically. The fold system in phase space is consistent with analytic self similar solutions, the solutions present all the proper self-similar scalings. An additional point is the appearance of an $x^{-(1/2)}$ law at the center of the system for initial conditions with power law index larger than $-(1/2)$. It is found that the first appearance of the $x^{-(1/2)}$ law corresponds to the formation of a singularity very close to the center. Finally the general properties of self similar multi dimensional solutions near equilibrium are investigated. Smooth and continuous self similar solutions have power law behavior at equilibrium. However cold initial conditions result in discontinuous phase space solutions, and the smoothed phase space density looses its auto similar properties. This problem is easily solved by observing that the probability distribution of the phase space density $P$ is identical except for scaling parameters to the probability distribution of the smoothed phase space density $P_S$. As a consequence $P_S$ inherit the self similar properties of $P$. This particular property is at the origin of the universal power law observed in numerical simulation for ${\rho}/{\sigma^3}$. The self similar properties of $P_S$ implies that other quantities should have also an universal power law behavior with predictable exponents. This hypothesis is tested using a numerical model of the phase space density of cold dark matter halo's, an excellent agreement is obtained.

We present the Strong Lensing Legacy Survey - ARCS (SARCS) sample compiled from the final T0006 data release of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) covering a total non-overlapping area of 159 sq.deg. We adopt a semi-automatic method to find gravitational arcs in the survey that makes use of an arc-finding algorithm. The candidate list is pruned by visual inspection and ranking to form the SARCS sample. This list also includes some serendipitously discovered lens candidates. The SARCS sample consists of 127 lens candidates which span arc radii~2"-18" within the unmasked area of ~ 150 sq. deg. Within the sample, 54 systems are promising lenses amongst which, we find 12 giant arcs and 2 radial arc candidates. From our sample, we detect a systematic alignment of giant arcs with the ellipticity of the baryonic component of the lens in concordance with previous studies. The lens redshift distribution corresponding to both the giant arcs and all arcs, estimated from photometric catalogs, peaks at z~0.6. Owing to the large area and depth of the CFHTLS, we find the largest lens sample probing group-scales for the first time. We compare the observed image separation distribution (ISD) of our arcs with theoretical models. A two-component density profile for the lenses which accounts for both the central galaxy and dark matter component is required by the data to explain the observed ISD. Unfortunately, current levels of uncertainties and degeneracies accommodate models both with and without adiabatic contraction. We also show the effects of changing parameters of the model that predict the ISD and that a larger lens sample might constrain relations such as the concentration-mass relation, mass-luminosity relation and slope of the luminosity function. (abridged)

This paper presents a principal components analysis of rotation curves from a sample of low surface brightness galaxies. The physical meaning of the principal components is investigated, and related to the intrinsic properties of the galaxies. The rotation curves are re-scaled using the optical disk scale, the resulting principal component decomposition demonstrates that the whole sample is properly approximated using two components. The ratio of the second to the first component is related to the halo steepness in the central region, is correlated to the gas fraction in the galaxy, and is un-correlated to other parameters. As a consequence the gas fraction appear as a fundamental variable with respect to the galaxies rotation curves, and its correlation with the halo steepness is especially important. Since the gas fraction is related to the degree of galaxy evolution, it is very likely that the steepness of the halo at the center is a consequence of galaxy evolution. More evolved galaxies have shallower central profile and statistically less gas, most likely as a consequence of more star formation and supernovae. The differences in evolution, gas fractions and halo central steepness of the galaxies could be due to the influence of different environments.

In gravitational lensing the average colors of the images are not identical to the average color of the source. The highly non-linear mapping of gravitational lensing does not preserve the color balance of the source, and this mapping is different for each image. The color distortion of the images is illustrated using HST images of the lens SL2SJ02140. It is shown that in this lens the color of the images is variable, reflecting the variable color of the source. The average color of the images in SL2SJ02140 are interpreted as a variable amplification of different sources regions with different colors. The variation of the average image colors affects the measurements of the photometric redshift of the images. This is especially true for SL2SJ02140 where the color variations due to the non-linear mapping of the lens simulates pseudo redshifts variations.

This paper presents a reconstruction of the gravitational lens SL2S02176-0513 using the singular perturbative method presented in Alard 2007, MNRAS Letters, 382, 58 and Alard, C., 2008, MNRAS, 388, 375. The ability of the perturbative method to separate the inner and outer contributions of the potential in gravitational lenses is tested using SL2S02176-0513. In this lens, the gravitational field of the central galaxy is dominated by a nearby group of galaxies located at a distance of a few critical radius. The perturbative functionals are re-constructed using local polynomials. The polynomial interpolation is smoothed using Fourier series, and numerically fitted to HST data using a non-linear minimization procedure. The potential inside and outside the critical circle is derived from the reconstruction of the perturbative fields. The inner and outer potential contours are very different.The inner contours are consistent with the central galaxy, while the outer contours are fully consistent with the perturbation introduced by the group of galaxies. The ability of the perturbative method to separate the inner and outer contribution is confirmed, and indicates that in the perturbative approach the field of the central deflector can be separated from outer perturbations. The separation of the inner and outer contribution is especially important for the study of the shape of dark matter halo's as well as for the statistical analysis of the effect of dark matter substructures.

The structure and potential of a complex gravitational lens is reconstructed using the perturbative method presented in Alard 2007, MNRAS, 382L, 58; Alard 2008, MNRAS, 388, 375. This lens is composed of 6 galaxies belonging to a small group. The lens inversion is reduced to the problem of reconstructing non-degenerate quantities: the 2 fields of the perturbative theory of strong gravitational lenses. Since in the perturbative theory the circular source solution is analytical, the general properties of the perturbative solution can be inferred directly from the data. As a consequence, the reconstruction of the perturbative fields is not affected by degeneracy, and finding the best solution is only a matter of numerical refinement. The local shape of the potential and density of the lens are inferred from the perturbative solution, revealing the existence of an independent dark component that does not follow light. The most likely explanation is that the particular shape of the dark halo is due to the merging of cold dark matter halos. This is a new result illustrating the structure of dark halos at the scale of galaxies.

A new multi-dimensional Hierarchical Structure Finder (HSF) to study the phase-space structure of dark matter in N-body cosmological simulations is presented. The algorithm depends mainly on two parameters, which control the level of connectivity of the detected structures and their significance compared to Poisson noise. By working in 6D phase-space, where contrasts are much more pronounced than in 3D position space, our HSF algorithm is capable of detecting subhaloes including their tidal tails, and can recognise other phase-space structures such as pure streams and candidate caustics. If an additional unbinding criterion is added, the algorithm can be used as a self-consistent halo and subhalo finder. As a test, we apply it to a large halo of the Millennium Simulation, where 19 % of the halo mass are found to belong to bound substructures, which is more than what is detected with conventional 3D substructure finders, and an additional 23-36 % of the total mass belongs to unbound HSF structures. The distribution of identified phase-space density peaks is clearly bimodal: high peaks are dominated by the bound structures and low peaks belong mostly to tidal streams. In order to better understand what HSF provides, we examine the time evolution of structures, based on the merger tree history. Bound structures typically make only up to 6 orbits inside the main halo. Still, HSF can identify at the present time at least 80 % of the original content of structures with a redshift of infall as high as z <= 0.3, which illustrates the significant power of this tool to perform dynamical analyses in phase-space.

This paper reviews and analyses methods used to identify neighbours in 6D space and estimate the corresponding phase-space density. It compares SPH methods to 6D Delaunay tessellation on statical and dynamical realisation of single halo profiles, paying attention to the unknown scaling, S_G, used to relate the spatial dimensions to the velocity dimensions. The methods with local adaptive metric provide the best phase-space estimators. They make use of a Shannon entropy criterion combined with a binary tree partitioning and with SPH interpolation using 10-40 neighbours. Local scaling implemented by such methods, which enforces local isotropy of the distribution function, can vary by about one order of magnitude in different regions within the system. It presents a bimodal distribution, in which one component is dominated by the main part of the halo and the other one is dominated by the substructures. While potentially better than SPH techniques, since it yields an optimal estimate of the local softening volume (and the local number of neighbours required to perform the interpolation), the Delaunay tessellation in fact poorly estimates the phase-space distribution function. Indeed, it requires, the choice of a global scaling S_G. We propose two methods to estimate S_G that yield a good global compromise. However, the Delaunay interpolation still remains quite sensitive to local anisotropies in the distribution. We also compare 6D phase-space density estimation with the proxy, Q=rho/sigma^3, where rho is the local density and 3 sigma^2 is the local 3D velocity dispersion. We show that Q only corresponds to a rough approximation of the true phase-space density, and is not able to capture all the details of the distribution in phase-space, ignoring, in particular, filamentation and tidal streams.

We present a statistical study of the effects induced by substructures on the deflection potential of dark matter halos in the strong lensing regime. This investigation is based on the pertubative solution around the Einstein radius (Alard 2007) in which all the information on the deflection potential is specified by only a pair of one-dimensional functions on this ring. Using direct comparison with ray-tracing solutions, we found that the iso-contours of lensed images predicted by the pertubative solution is reproduced with a mean error on their radial extension of less than 1% - in units of the Einstein radius, for reasonable substructure masses. It demonstrates the efficiency of the approximation to track possible signatures of substructures. We have evaluated these two fields and studied their properties for different lens configurations modelled either through massive dark matter halos from a cosmological N-body simulation, or via toy models of Monte Carlo distribution of substructures embedded in a triaxial Hernquist potential. As expected, the angular power spectra of these two fields tend to have larger values for larger harmonic numbers when substructures are accounted for and they can be approximated by power-laws, whose values are fitted as a function of the profile and the distribution of the substructures.

In the perturbative approach, substructures in the lens can be reduced to their effect on the two perturbative fields $f_1$ and $\frac{d f_0}{d\theta}$. A simple generic model of elliptical lens with a substructure situated near the critical radius is investigated in details. Analytical expressions are derived for each perturbative field, and basic properties are analyzed. The power spectrum of the fields is well approximated by a power-law, resulting in significant tails at high frequencies. Another feature of the perturbation by a substructure is that the ratio of the power spectrum at order $n$ of the 2 fields $R_n$ is nearly 1. The ratio $R_n \simeq 1$ is specific to substructures, for instance an higher order distortion ($n>2$) but with auto-similar isophotes will result in $R_n \propto \frac{1}{n^2}$. Finally, the problem of reconstructing the perturbative field is investigated. Local field model are implemented and fitted to maximize image similarity in the source plane. The non-linear optimization is greatly facilitated, since in the perturbative approach the circular source solution is always known. Examples of images distortions in the subcritical regime due to substructures are presented, and analyzed for different source shapes. Provided enough images and signal is available, the substructure field can be identified confidently. These results suggests that the perturbative method is an efficient tool to estimate the contribution of substructures to the mass distribution of lenses.

The image of a point situated at the center of a circularly symmetric potential is a perfect circle. The perturbative effect of non-symmetrical potential terms is to displace and break the perfect circle. These 2 effects, displacement and breaking are directly related to the Taylor expansion of the perturbation at first order on the circle. The numerical accuracy of this perturbative approach is tested in the case of an elliptical potential with a core radius. The contour of the images and the caustics lines are well re-produced by the perturbative approach. These results suggests that the modeling of arcs, and in particular of tangential arcs may be simplified by using a general perturbative representation on the circle. An interesting feature of the perturbative approach, is that the equation of the caustic line depends only on the values on the circle of the lens displacement field along the $\theta$ direction.

AIMS: We present data from the CFHTLS Strong Lensing Legacy Survey (SL2S). Due to the unsurpassed combined depth, area and image quality of the Canada-France-Hawaii Legacy Survey it is becoming possible to uncover a large, statistically well-defined sample of strong gravitational lenses which spans the dark halo mass spectrum predicted by the concordance model from galaxy to cluster haloes. METHODS: We describe the development of several automated procedures to find strong lenses of various mass regimes in CFHTLS images. RESULTS: The preliminary sample of about 40 strong lensing candidates discovered in the CFHTLS T0002 release, covering an effective field of view of 28 deg$^2$ is presented. These strong lensing systems were discovered using an automated search and consist mainly of gravitational arc systems with splitting angles between 2 and 15 arcsec. This sample shows for the first time that it is possible to uncover a large population of strong lenses from galaxy groups with typical halo masses of about $10^{13}h^{-1}M_\odot$. We discuss the future evolution of the SL2S project and its main scientific aims for the next 3 years, in particular our observational strategy to extract the hundreds of gravitational rings also present in these fields.

This paper presents a method to identify gravitational arcs or more generally elongated structures in a given image. The method is based on the computation of a local estimator of the elongation. The estimation of the local elongation proceed in two steps: first the local orientation of the structure is computed, then in the next step, a rotation is performed, and the marginal distributions are used to compute the elongation. This procedure allows the computation of the local elongation at each point in the image. Then, using a threshold on the elongation map the elongated structures are identified and re-constructed using connectivity criteria. Finally a catalogue of elongated structures is produced, and the properties of each object are computed, allowing the selection of potential arc candidates. The final selection of the arc candidates is performed by visual inspection of multi-color images of a small number of objects. This method is a general tool that may be applied not only to gravitational arcs, but to all problems related to the mapping and measurement of elongated structures, in an image, or a volume.

We propose to integrate the Vlasov-Poisson equations giving the evolution of a dynamical system in phase-space using a continuous set of local basis functions. In practice, the method decomposes the density in phase-space into small smooth units having compact support. We call these small units ``clouds'' and choose them to be Gaussians of elliptical support. Fortunately, the evolution of these clouds in the local potential has an analytical solution, that can be used to evolve the whole system during a significant fraction of dynamical time. In the process, the clouds, initially round, change shape and get elongated. At some point, the system needs to be remapped on round clouds once again. This remapping can be performed optimally using a small number of Lucy iterations. The remapped solution can be evolved again with the cloud method, and the process can be iterated a large number of times without showing significant diffusion. Our numerical experiments show that it is possible to follow the 2 dimensional phase space distribution during a large number of dynamical times with excellent accuracy. The main limitation to this accuracy is the finite size of the clouds, which results in coarse graining the structures smaller than the clouds and induces small aliasing effects at these scales. However, it is shown in this paper that this method is consistent with an adaptive refinement algorithm which allows one to track the evolution of the finer structure in phase space. It is also shown that the generalization of the cloud method to the 4 dimensional and the 6 dimensional phase space is quite natural.

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.

The image subtraction method is a powerful tool to analyze the light variations in crowded fields. This method is able to achieve a nearly optimal differential photometry, even in very dense regions. However, image subtraction is not limited to photometry, and it is shown in this paper that the method can be generalized to the measurement of the differential proper motions. It is important to emphasize that image subtraction can re-construct an un-biased frame to frame astrometric transform. A nearly optimal determination of this astrometric transform can be performed by expanding the spatial variations of the kernel using the derivatives of the constant kernel solution. It is demonstrated that an expansion using first and second derivatives is optimal. Differential refraction can be corrected easily using an artificial image and the first derivatives of the kernel. To illustrate the ability of the method to measure proper motions, a small sub-area near the center of a Galactic Bulge field covered 265 times by the OGLE II experiment has been selected. The resulting astrometric accuracy is very close to the photon noise for the faint objects, while for the brighter ones it approaches closely the limit set by the residual seeing fluctuations on the smaller scales. The accuracy obtained with this new method is compared to the accuracy achieved using classical methods. The improvement obtained varies from a factor of $\simeq$ 2 for the brighter objects, to more than a factor of 3 for the fainter ones. This improvement is typical of the improvement that was obtained for the measurement of the photometric variations.

The ISOGAL project is an infrared survey of specific regions sampling the Galactic Plane selected to provide information on Galactic structure,stellar populations,stellar mass-loss and the recent star formation history of the inner disk and Bulge of the Galaxy. ISOGAL combines 7 and 15 micron ISOCAM observations - with a resolution of 6'' at worst - with DENIS IJKs data to determine the nature of the sources and theinterstellar extinction. We have observed about 16 square degrees with a sensitivity approaching 10-20mJy, detecting ~10^5 sources,mostly AGB stars,red giants and young stars. The main features of the ISOGAL survey and the observations are summarized in this paper,together with a brief discussion of data processing and quality. The primary ISOGAL products are described briefly (a full description is given in Schuller et al. 2003, astro-ph/0304309): viz. the images and theISOGAL-DENIS five-wavelength point source catalogue. The main scientific results already derived or in progress are summarized. These include astrometrically calibrated 7 and 15um images,determining structures of resolved sources; identification and properties of interstellar dark clouds; quantification of the infrared extinction law and source dereddening; analysis of red giant and (especially) AGB stellar populations in the central Bulge,determining luminosity,presence of circumstellar dust and mass--loss rate,and source classification,supplemented in some cases by ISO/CVF spectroscopy; detection of young stellar objects of diverse types,especially in the inner Bulge with information about the present and recent star formation rate; identification of foreground sources with mid-IR excess. These results are the subject of about 25 refereed papers published or in preparation.

We present version 1.0 of the ISOGAL-DENIS Point Source Catalogue (PSC), containing more than 100,000 point sources detected at 7 and/or 15 micron in the ISOGAL survey of the inner Galaxy with the ISOCAM instrument on board the Infrared Space Observatory (ISO). These sources are cross-identified, wherever possible, with near-infrared (0.8-2.2 micron) data from the DENIS survey. The overall surface covered by the ISOGAL survey is about 16 square degrees, mostly (95%) distributed near the Galactic plane (|b| < 1 deg), where the source extraction can become confusion limited and perturbed by the high background emission. Therefore, special care has been taken aimed at limiting the photometric error to ~0.2 magnitude down to a sensitivity limit of typically 10 mJy. The present paper gives a complete description of the entries and the information which can be found in this catalogue, as well as a detailed discussion of the data processing and the quality checks which have been completed. The catalogue is available via the VizieR Service at the Centre de Donn\'ees Astronomiques de Strasbourg (CDS, http://vizier.u-strasbg.fr/viz-bin/VizieR/) and also via the server at the Institut d'Astrophysique de Paris (http://www-isogal.iap.fr/). A more complete version of this paper, including a detailed description of the data processing, is available in electronic form through the ADS service.

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.

A map of the projected density of the old stellar population of the Galactic Bulge region is reconstructed using 2MASS data. By making a combination of the H and K photometric bands, it is possible to overcome the effect of reddening, and thus penetrate the inner structure of the Galactic Bulge. The main structure in the map corresponds to the well documented peanut shaped bar which is formed by the inner parts of the Galactic disk as a result of dynamical instabilities. As suggested by numerical simulations, the projected Z profile of the bar, has an almost exponential shape. After subtracting the exponential profile associated with the bar, a large residual appear near the Galactic Center. This residual is elongated and asymmetrical, which suggest a bar structure. Thus we arrive at the conclusion that in addition to the main bar a smaller bar with a different orientation may exist in the central region of the Milky Way. This finding makes the Milky Way very similar to a large number of barred spiral Galaxies which show as well a smaller bar in their central regions.

Metallicity gradients in the Sagittarius dwarf Galaxy (Sgr) are investigated by using infrared photometric data from the 2MASS survey. To search for metallicity effects, the giant branch in a field situated near the Center of the Sgr is compared to the giant branch in a field situated near its southern edge. The contamination of Sgr giant branch by foreground Galactic stars is canceled by statistical subtraction of diagrams symmetrical in Galactic latitude. After subtraction it is possible to reconstruct the Sgr giant branch with excellent accuracy. The giant branch in the two fields have similar slopes but are shifted in color. Even after correction for the differential reddening between the fields, the shift in color between the branch remains, and is very significant. This variation in the color of the giant branch corresponds to a metallicity variation of about -0.25 Dex. The existence of a metallicity gradient in Sgr may indicate that there are two different stellar population in Sgr. One has low metallicity, and another one of higher metallicity has a smaller spatial extension.

We estimate the optical depth, time-scale distribution and fraction of microlensing events originating from sources in the Sagittarius dwarf galaxy (Sgr) lensed by deflectors in the Milky Way. These events have a time-scale longer by a factor ~1.3 than the MW/MW events and occur mainly on sources fainter than V~21 mag below Sgr's turn off. The fraction of events involving a source in Sgr depends on the location and extinction of the field and on the limiting magnitude of the survey. The contribution of the MW/Sgr events is negligible (<1%) at very low latitudes (|b|<2 deg.) but increases continuously towards higher |b| and becomes dominant near the highest density region of the dwarf galaxy. Sgr is present within the fields of current microlensing surveys and any optical depth map inferred from observations will become biased by the presence of Sgr towards higher |b| where the contribution of MW/Sgr events is significant. Systematic spectroscopic measurements on the sources of all the microlensing events may allow detection of this kind of event for which the degeneracy on the lens mass can be significantly reduced.

By cross-correlating the results of two recent large-scale surveys, the general properties of a well defined sample of semi-regular variable stars have been determined. ISOGAL mid-infrared photometry and MACHO lightcurves are assembled for approximately 300 stars in the Baade's Windows of low extinction towards the Galactic bulge. These stars are mainly giants of late M spectral type, evolving along the asymptotic giant branch (AGB). They are found to possess a wide and continuous distribution of pulsation periods and to obey an approximate log~period -- bolometric magnitude relation or set of such relations. Approximate mass-loss rates in the range of 1e-8 to 5e-7 M_sun per year are derived from ISOGAL mid-infrared photometry and models of stellar spectra adjusted for the presence of optically-thin circumstellar silicate dust. Mass-loss rates depend on luminosity and pulsation period. Some stars lose mass as rapidly as short-period Miras but do not show Mira-like amplitudes. A period of 70 days or longer is a necessary but not a sufficient condition for mass loss to occur. For AGB stars in the mass-loss ranges that we observe, the functional dependence of mass-loss rate on temperature and luminosity is found to be in agreement with recent theoretical predictions. If we include our mass-loss rates with a sample of extreme mass-losing AGB stars in the Large Magellanic Cloud, we get the general result for AGB stars that mass-loss rate is proportional to luminosity^{2.7}, valid for AGB stars with 10^{-8} to 10^{-4} M_sun per year (Abridged).

We present the discovery of a previously unknown member of the immediate solar neighbourhood, DENIS-P J104814.7-395606.1 (hereafter DENIS 1048-39), identified while mining the DENIS database for new nearby stars. A HIRES echelle spectrum obtained with the 10-m Keck telescope shows that it is an M9 dwarf. DENIS 1048-39 has a very bright apparent magnitude (I=12.67) for its spectral type and colour (I-J=3.07), and is therefore very nearby. If it is single its distance is only 4.1 +- 0.6pc, ranking it as between our twelfth and fortyth closest neighbour. It is also the closest star or brown dwarf with a spectral type later than M7V. Its proper motion was determined through comparison of Sky atlas Schmidt plates, scanned by the MAMA microdensitometer, with the DENIS images. At 1.53"/yr it further attests the closeness of DENIS 1048-39, and hence its dwarf status. These characteristics make it an obvious target for further detailed studies.

We have compiled the near infrared Point Source Catalogue (PSC) towards the Magellanic Clouds (MCs) extracted from the data obtained with the Deep Near Infrared Survey of the Southern Sky - DENIS (Epchtein et al. 1997). The catalogue covers an area of of 19.87* 16 square degrees centered on (RA, DEC)=(5h27m20s, -69o00'00'') for the Large Magellanic Cloud (LMC) and 14.7* 10 square degrees centered on (RA, DEC)=(h02m40s, -73o00'00'' for the Small Magellanic Cloud (SMC) at the epoch J2000. It contains about 1300000 sources towards the LMC and 300000 sources towards the SMC each detected in at least 2 of the 3 photometric bands involved in the survey (I, J, Ks). 70% of the detected sources are true members of the Magellanic Clouds, respectively and consist mainly of red giants, asymptotic giant branch stars and super-giants. The observations have all been made with the same instrument and the data have been calibrated and reduced uniformly. The catalogue provides a homogeneous set of photometric data.

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

This paper presents an investigation of the outer disk structure by using data from a recent release of the 2 micron sky survey (2MASS). This 2MASS data show unambiguously that the stellar disk thickens with increasing distance from the Sun. In one of the field (longitude l$=$240) there is also strong evidence of an assymetry associated with the Galactic warp. This flaring and warping of the stellar disk is very similar to the features observed in the HI disk. The thickening of the stellar disk explains the drop in density observed near the Galactic plane: stars located at lower Galactic latitudes are re-distributed to higher latitudes. It is no longer necessary to introduce a disk cut-off to explain the drop in density. It is also not clear how this flaring disk is distinct from the thick disk in the outer disk region. At least, for lines of sight in the direction of the outer disk, the thickening of the disk is sufficient to account for the excess in star counts attributed by some models to the thick disk.

The large number of microlensing events discovered towards the Galactic Bulge bears the promise to reconstruct the stellar mass function. The more interesting issue concerning the mass function is certainly to probe its low mass end. However due to the source confusion, even if the distribution and the kinematics of the lenses are known, the estimation of the mass function is extremely biased at low masses. The blending due to the source confusion biases the duration of the event, which in turn dramatically biases the estimation of the mass of the lens. To overcome this problem we propose to use differential photometry of the microlensing events obtained using the image subtraction method. Differential photometry is free of any bias due to blending, however the drawback of differential photometry is that the baseline flux is unknown. In this paper we will show that even without knowing the baseline flux, purely differential photometry allow to estimate the mass function without any biases. The basis of the method is that taking the scalar product of the microlensing light curves with a given function and taking its sum over all the microlensing events is equivalent to project the mass function on another function. This method demonstrates that there is a direct correspondancy between the space of the observations and the space of the mass function. Concerning the function to use in order to project the observations, we show that the principal components of the light curves are an optimal set. To illustrate the method we simulate sets consistent with the microlensing experiments. By using 1000 of these simulations, we show that for instance the exponent of the mass function can be reconstructed without any biases.

We report the detection of $\sim$ 1500 RR Lyrae of Bailey type ab located in the Sagittarius dwarf galaxy (Sgr). These variables have been detected on two ESO Schmidt fields centred on (l,b)=(3.1 deg,-7.1 deg) and (6.6 deg, -10.8 deg), covering an area of $\sim$50 square degree. We present a surface density map of Sgr based on the spatial distribution of these RRab, allowing us to trace its structure in a region that was still almost unexplored between b=-14 deg and b=-4 deg. We present the results of the fit of different models to the density profile of Sgr. The best fit to the core of Sgr is an exponential with a scale length of 4.1 deg along the major axis. When we look at the extension of Sgr we find a break (significant at the $\sim2\sigma$ level) in the slope of the surface density along the main axis of Sgr. The nearly flat (or at least very slowly decreasing) profile in the outer region of Sgr shows that this dwarf galaxy is probably extending even further out our fields.

DENIS observations in the J (1.2 micron) and K_S (2.15 micron) bands together with isochrones calculated for the RGB and AGB phase are used to draw an extinction map of the inner Galactic Bulge. The uncertainty in this method is mainly limited by the optical depth of the Bulge itself. A comparison with fields of known extinction shows a very good agreement. We present an extinction map for the inner Galactic Bulge (approx. 20 sq. deg.)

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.

The ISOGAL project is a survey of the stellar populations, structure, and recent star formation history of the inner disk and bulge of the Galaxy. ISOGAL combines 15 and 7micron ISOCAM observations with DENIS IJKs data to determine the nature of a source and the interstellar extinction. In this paper we report an ISOGAL study of a small field in the inner Galactic Bulge (l=0deg, b=1.0deg, area=0.035 sq. deg) as a prototype of the larger area ISOGAL survey of the inner Galaxy. The five wavelengths of ISOGAL+DENIS, together with the relatively low and constant extinction in front of this specific field, allow reliable determination of the nature of the sources. The primary scientific result of this paper is evidence that the most numerous class of ISOGAL 15micron sources are Red Giants in the Galactic bulge and central disk, with luminosities just above or close to the RGB tip and weak mass-loss rates. They form loose sequences in the magnitude-colour diagrams [15]/Ks-[15] and [15]/[7]-[15]. Their large excesses at 15micron with respect to 2micron and 7micron is due to circumstellar dust produced by mass-loss at low rates. These ISOGAL results are the first systematic evidence and study of dust emission at this early stage (''Intermediate'' AGB), before the onset of the large mass-loss phase. It is thus well established that efficient dust formation is already associated with such low mass-loss rates during this early phase.

We present $V$-band light curves of 61 variables from the core of the globular cluster M5 obtained using a newly developed image subtraction method (ISM). Four of these variables were previously unknown. Only 26 variables were found in the same field using photometry obtained with DoPHOT software. Fourier parameters of the ISM light curves have relative errors up to 20 times smaller than parameters measured from DoPHOT photometry. We conclude that the new method is very promising for searching for variable stars in the cores of the globular clusters and gives very accurate relative photometry with quality comparable to photometry obtained by HST. We also show that the variable V104 is not an eclipsing star as has been suggested, but is an RRc star showing non-radial pulsations.

We present results from 2 years of monitoring of Huchra's lens (QSO 2237+0305) with the 1.3 m Warsaw telescope on Las Campanas, Chile. Photometry in the V band was done using a newly developed method for image subtraction. Reliable subtraction without Fourier division removes all complexities associated with the presence of a bright lensing galaxy. With positions of lensed images adopted from HST measurements it is relatively easy to fit the variable part of the flux in this system, as opposed to modeling of the underlying galaxy. For the first time we observed smooth light variation over a period of a few months, which can be naturally attributed to microlensing. We also describe automated software capable of real time analysis of the images of QSO 2237+0305. It is expected that starting from the next observing season in 1999 an alert system will be implemented for high amplification events (HAE) in this object. Time sampling and photometric accuracy achieved should be sufficient for early detection of caustic crossings.

The ISOGAL mid-infrared survey of areas close to the Galactic Plane aims to determine their stellar content and its possible bearing on the history of the Galaxy. The NGC6522 and Sgr I Baade's Windows of low obscuration towards the inner parts of the Bulge represent ideal places in which to calibrate and understand the ISOGAL colour-magnitude diagrams. The survey observations were made with the ISOCAM instrument of the ISO satellite. The filter bands chosen were LW2(~7 microns) and LW3(~15 microns). The results presented here show that most of the detected objects are late M-type giants on the AGB, with a cut-off for those earlier than M3-M4. The most luminous members of these two Bulge fields at 7 microns are the Mira variables. However, it is evident that they represent the end of a sequence of increasing 15 micron dust emission which commences with M giants of earlier sub-type. In observations of late-type giants the ISOCAM 15 micron band is mainly sensitive to the cool silicate or aluminate dust shells which overwhelm the photospheric emission. However, in ordinary M-giant stars, the 7 micron band is not strongly affected by dust emission and may be influenced instead by absorption. The nu2 band of water at 6.25 microns and the SiO fundamental at 7.9 microns are likely contribitors to this effect. A group of late M stars has been found which vary little or not at all but have infrared colours typical of well-developed dust shells. Their luminosities are similar to those of 200-300 day Miras but they have slightly redder [7]-[15] colours which form an extension of the ordinary M giant sequence. The Mira dust shells show a mid-infrared [7]-[15] colour-period relation. ca 700 days.

It is demonstrated that non-constant kernel solution, that can fit the spatial variations of the kernel can be obtained with minimum computing time. The CPU cost required with this new extension of the image subtraction method is almost the same as for fitting constant kernel solution. The method is demonstrated with a serie of Monte-Carlo images. Differential PSF variations and differential rotation between the images are simulated. It is shown that this new method is able to achieve an optimal results even in these difficult cases. It is clear that the most frequents instrumental problems are automatically corrected by the method. It is also demonstrated that the method does not suffer problems with under-sampling of the images. To conclude the method is applied to real images of crowded field. It is proved that much larger sub-areas of the images can be used for the fit, while keeping the same accuracy in the subtracted image. This result is especially important in case of variables located in low density fields, like the Huchra lens. Many other useful applications of the method are possible for major astrophysical problems. For instance supernovas search should benefit the method, but other applications like Cepheids surveys in other galaxies should also find some advantage in using the method. Many other applications will certainly show-up, since variability search is a major issue in astronomy.

The light curves of the OGLE microlensing candidates have been reconstructed using the image subtraction method. A large improvement of the photometric accuracy has been found in comparison with previous processing of the data with DoPHOT. On the mean, the residuals to the fit of a microlensing light curve are improved by a factor of 2 for baseline data points, and by a factor of 2.5 during magnification. The largest improvement was found for the OGLE #5 event, where we get an accuracy 7.5 times better than with DoPHOT. Despite some defects in the old CCD used during the OGLE I experiment we obtain most of the time errors that are only 30 % to 40 % in excess of the photon noise. Previous experiment showed that with modern CCD chips (OGLE II), residuals much closer to the photon noise were obtained. The better photometric quality enabled us to find a low amplitude, long term variability in the OGLE #12 and OGLE #11 baseline magnitude. We also found that the shape of the OGLE #14 candidate light curve is fairly inconsistent with microlensing of a point source by a point lens. A dramatic change in the light curve of the OGLE #9 candidate was also found, which indicates that very large biases can be present in data processed with DoPHOT. To conclude we made a detailed analysis of the blending issue. It is found that OGLE #5 and OGLE #6 are very likely highly blended microlensing events. These events result from the magnification of a faint star that would have been undetectable without microlensing.

We present a new method designed for optimal subtraction of two images with different seeing. Using image subtraction appears to be essential for the full analysis of the microlensing survey images, however a perfect subtraction of two images is not easy as it requires the derivation of an extremely accurate convolution kernel. Some empirical attempts to find the kernel have used the Fourier transform of bright stars, but solving the statistical problem of finding the best kernel solution has never really been tackled. We demonstrate that it is possible to derive an optimal kernel solution from a simple least square analysis using all the pixels of both images, and also show that it is possible to fit the differential background variation at the same time. We also show that PSF variations can also be easily handled by the method. To demonstrate the practical efficiency of the method, we analyzed some images from a Galactic Bulge field monitored by the OGLE II project. We find that the residuals in the subtracted images are very close to the photon noise expectations. We also present some light curves of variable stars, and show that, despite high crowding levels, we get an error distribution close to that expected from photon noise alone. We thus demonstrate that nearly optimal differential photometry can be achieved even in very crowded fields. We suggest that this algorithm might be particularly important for microlensing surveys, where the photometric accuracy and completeness levels could be very significantly improved by using this method.

Previous calculations of the rates and optical depths due to microlensing only considered resolved stars. However, if a faint unresolved star lens is close enough to a resolved star, the event will be seen by the microlensing experiments and attributed to the bighter star. The blending biases the duration, making the contribution of the unresolved stars very significant for short events. This contribution is confused with lensing by brown dwarfs. The exact rates of these blended events are extremly sensitive to the limiting magnitude achieved in the microlensing search. Appropriate calculations of the optical depth and rates are provided here, and illustrated in the case of the DUO and OGLE experiments. The additional contribution of unresolved stars is very significant and probably explains the high optical depth and rates observed towards the Galactic Bulge. The blended unresolved event can be identified using either the color shift or the light curve shape. However, neither of these two methods is apropriate to identify a large number of blended events towards the Bulge. In some cases of good photometry and small impact parameter, an identification is possible, as for the OGLE 5 event, which clearly appears as a case of lensing of an unresolved star. The recent results obtained by the PLANET collaboration indicate that a high resolution and dense sampling of the light curve is possible, and will probably provide a very interesting possibility to correct the blending bias, as demonstrated for OGLE 5. This possibility, is certainly better than a statistical estimation of the lensing rates, which are always prone to some uncertainty. But, at this time, the contribution of unresolved stars must be included in the analyses of microlensing experiments.

We present the light curve of an unusual variable object, DUO 2, detected during the search for microlensing events by the DUO project. The star remained stable for more than 150 days before it brightened by more than two magnitudes in 6 days in the B and R bands. The light curves are achromatic during the variability. We consider possible explanations of the photometric behavior, with particular emphasis on the binary lens interpretation of the event. The masses of the lenses are quite small, with the companion possibly in the range of a brown dwarf or even a few times of Jupiter. We report evidence of blending of the source by a companion through the first detection of shift in the light centroid among all the microlensing experiments. This shift sets a lower limit of $0.3^{\prime\prime}$ on the separation between the stars. The best lens model obtained requires moderate blending, which was what motivated us to check the centroid shift that was subsequently found. The best lens model predicts a separation of $1^{\prime\prime}$ between the two blended stars. This prediction was recently tested using two CCD images taken under good seeing conditions. Both images show two components. Their separation and position angle are in good agreement with our model.