Gravitational microlensing is a key probe of the nature of dark matter and its distribution on the smallest scales. For many practical purposes, confronting theory to observation requires to model the probability that a light source is highly amplified by many-lens systems. This article reviews four simple analytic models of the amplification probability distribution, based on different approximations: (i) the strongest-lens model; (ii) the multiplicative model, where the total amplification is assumed to be the product of all the lenses' individual amplifications; (iii) a hybrid version of the previous two; and (iv) an empirical fitting function. In particular, a new derivation of the multiplicative amplification distribution is proposed, thereby correcting errors in the literature. Finally, the accuracy of these models is tested against ray-shooting simulations. They all produce excellent results as long as lenses are light and rare (low optical depth); however, for larger optical depths, none of them succeeds in capturing the relevant features of the amplification distribution. This conclusion emphasizes the crucial role of lens-lens coupling at large optical depths.

The weak equivalence principle is one of the cornerstone of general relativity. Its validity has been tested with impressive precision in the Solar System, with experiments involving baryonic matter and light. However, on cosmological scales and when dark matter is concerned, the validity of this principle is still unknown. In this paper we construct a null test that probes the validity of the equivalence principle for dark matter. Our test has the strong advantage that it can be applied on data without relying on any modelling of the theory of gravity. It involves a combination of redshift-space distortions and relativistic effects in the galaxy number-count fluctuation, that vanishes if and only if the equivalence principle holds. We show that the null test is very insensitive to typical uncertainties in other cosmological parameters, including the magnification bias parameter, and to non-linear effects, making this a robust null test for modified gravity.

Standard sirens have been proposed as probes of alternative theories of gravity, such as Horndeski models. Hitherto, all studies have been conducted on a homogeneous-isotropic cosmological background, which is unable to consistently account for realistic distributions of matter, and for inhomogeneities in the Horndeski scalar field. Yet, the latter are essential for screening mechanisms. In this article, we comprehensively analyze the propagation of Horndeski gravitational waves in an arbitrary background spacetime and scalar field. We restrict to the class of theories in which gravitational waves propagate at light speed, and work in the geometric-optics regime. We find that kinetic braiding only produces a nonphysical longitudinal mode, whereas conformal coupling affects the amplitude of the standard transverse modes but not their polarization. We confirm that any observable deviation from general relativity depends on the local value of the effective Planck mass at emission and reception of the wave. This result is interpreted as the conservation of the number of gravitons.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The standard theory of weak gravitational lensing relies on the approximation that light beams are infinitesimal. Our recent work showed that the finite size of sources, and the associated light beams, can cause nonperturbative corrections to the weak-lensing convergence and shear. This article thoroughly investigates these corrections in a realistic cosmological model. The continuous transition from infinitesimal to finite beams is understood, and reveals that the previous results overestimated finite-beam effects due to simplistic assumptions on the distribution of matter in the Universe. In a KiloDegree Survey-like setting, finite-beam corrections to the cosmic shear remain subpercent, while percent-level corrections are only reached on subarcmin scales. This article thus demonstrates the validity of the infinitesimal-beam approximation in the interpretation of current weak-lensing data.

When a luminous source is extended, its distortions by weak gravitational lensing are richer than a mere combination of magnification and shear. In a recent work, we proposed an elegant formalism based on complex analysis to describe and calculate such distortions. The present article further elaborates this finite-beam approach, and applies it to a realistic cosmological model. In particular, the cosmic correlations of image distortions beyond shear are predicted for the first time. These constitute new weak-lensing observables, sensitive to very-small-scale features of the distribution of matter in the Universe. While the major part of the analysis is performed in the approximation of circular sources, a general method for extending it to noncircular sources is presented and applied to the astrophysically relevant case of elliptic sources.

In view of interpreting the cosmological observations precisely, especially when they involve narrow light beams, it is crucial to understand how light propagates in our statistically homogeneous, clumpy, Universe. Among the various approaches to tackle this issue, Swiss-cheese models propose an inhomogeneous spacetime geometry which is an exact solution of Einstein's equation, while the Dyer-Roeder approximation deals with inhomogeneity in an effective way. In this article, we demonstrate that the distance-redshift relation of a certain class of Swiss-cheese models is the same as the one predicted by the Dyer-Roeder approach, at a well-controlled level of approximation. Both methods are therefore equivalent when applied to the interpretation of, e.g., supernova observations. The proof relies on completely analytical arguments, and is illustrated by numerical results.

The equivalence principle, that is one of the main pillars of general relativity, is very well tested in the Solar system; however, its validity is more uncertain on cosmological scales, or when dark matter is concerned. This article shows that relativistic effects in the large-scale structure can be used to directly test whether dark matter satisfies Euler's equation, i.e. whether its free fall is characterised by geodesic motion, just like baryons and light. After having proposed a general parametrisation for deviations from Euler's equation, we perform Fisher-matrix forecasts for future surveys like DESI and the SKA, and show that such deviations can be constrained with a precision of order 10%. Deviations from Euler's equation cannot be tested directly with standard methods like redshift-space distortions and gravitational lensing, since these observables are not sensitive to the time component of the metric. Our analysis shows therefore that relativistic effects bring new and complementary constraints to alternative theories of gravity.

It has been claimed in Ref.[arXiv:1712.02240] that massive primordial black holes (PBH) cannot constitute all of the dark matter (DM), because their gravitational-lensing imprint on the Hubble diagram of type Ia supernovae (SN) would be incompatible with present observations. In this paper, we critically review those constraints and find several caveats on the analysis. First of all, the constraints on the fraction $\alpha$ of PBH in matter seem to be driven by a very restrictive choice of priors on the cosmological parameters. In particular, the degeneracy between $\Omega_{\rm M}$ and $\alpha$ is ignored and thus, by fixing $\Omega_{\rm M}$, transferred the constraining power of SN magnitudes to $\alpha$. Furthermore, by considering more realistic physical sizes for the type-Ia supernovae, we find an effect on the SN lensing magnification distribution that leads to significantly looser constraints. Moreover, considering a wide mass spectrum of PBH, such as a lognormal distribution, further softens the constraints from SN lensing. Finally, we find that the fraction of PBH that could constitute DM today is bounded by $f_{\rm PBH} < 1.09\ (1.38)$, for JLA (Union 2.1) catalogs, and thus it is perfectly compatible with an all-PBH dark matter scenario in the LIGO band.

The standard theory of weak gravitational lensing relies on the infinitesimal light beam approximation. In this context, images are distorted by convergence and shear, the respective sources of which unphysically depend on the resolution of the distribution of matter---the so-called Ricci-Weyl problem. In this letter, we propose a strong-lensing-inspired formalism to describe the lensing of finite beams. We address the Ricci-Weyl problem by showing explicitly that convergence is caused by the matter enclosed by the beam, regardless of its distribution. Furthermore, shear turns out to be systematically enhanced by the finiteness of the beam. This implies, in particular, that the Kaiser-Squires relation between shear and convergence is violated, which could have profound consequences on the interpretation of weak lensing surveys.

On small scales the observable Universe is highly inhomogeneous, with galaxies and clusters forming a complex web of voids and filaments. The optical properties of such configurations can be quite different from the perfectly smooth Friedmann-Lema\^itre-Robertson-Walker (FLRW) solutions that are frequently used in cosmology, and must be well understood if we are to make precise inferences about fundamental physics from cosmological observations. We investigate this problem by calculating redshifts and luminosity distances within a class of cosmological models that are constructed explicitly in order to allow for large density contrasts on small scales. Our study of optics is then achieved by propagating one hundred thousand null geodesics through such space-times, with matter arranged in either compact opaque objects or diffuse transparent haloes. We find that in the absence of opaque objects, the mean of our ray tracing results faithfully reproduces the expectations from FLRW cosmology. When opaque objects with sizes similar to those of galactic bulges are introduced, however, we find that the mean of distance measures can be shifted up from FLRW predictions by as much as $10\%$. This bias is due to the viable photon trajectories being restricted by the presence of the opaque objects, which means that they cannot probe the regions of space-time with the highest curvature. It corresponds to a positive bias of order $10\%$ in the estimation of $\Omega_{\Lambda}$ and highlights the important consequences that astronomical selection effects can have on cosmological observables.

Cosmic backreaction refers to the general question of whether a homogeneous and isotropic cosmological model is able to predict the correct expansion dynamics of our inhomogeneous Universe. One aspect of this issue concerns the validity of the continuous approximation: does a system of point masses expand the same way as a fluid does? This article shows that it is not exactly the case in Newtonian gravity, although the associated corrections vanish in an infinite Universe. It turns out that Gauss's law is a key ingredient for such corrections to vanish. Backreaction, therefore, generically arises in alternative theories of gravitation, which threatens the trustworthiness of their cosmological tests. This phenomenon is illustrated with a toy model of massive gravity.

The Hubble diagram is one of the cornerstones of observational cosmology. It is usually analysed assuming that, on average, the underlying relation between magnitude and redshift matches the prediction of a Friedmann-Lema\^itre-Robertson-Walker model. However, the inhomogeneity of the Universe generically biases these observables, mainly due to peculiar velocities and gravitational lensing, in a way that depends on the notion of average used in theoretical calculations. In this article, we carefully derive the notion of average which corresponds to the observation of the Hubble diagram. We then calculate its bias at second-order in cosmological perturbations, and estimate the consequences on the inference of cosmological parameters, for various current and future surveys. We find that this bias deeply affects direct estimations of the evolution of the dark-energy equation of state. However, errors in the standard inference of cosmological parameters remain smaller than observational uncertainties, even though they reach percent level on some parameters; they reduce to sub-percent level if an optimal distance indicator is used.

The time delay between the receptions of ultra-relativistic particles emitted simultaneously is a useful observable for both fundamental physics and cosmology. The expression of the delay when the particles travel through an arbitrary spacetime has been derived recently by Fanizza et al., using a particular coordinate system and self-consistent assumptions. The present article shows that this formula enjoys a simple physical interpretation: the relative velocity between two ultra-relativistic particles is constant. This result reveals an interesting kinematical property of general relativity, namely that the tidal forces experienced by ultra-relativistic particles in the direction of their motion are much smaller than those experienced orthogonally to their motion.

The geodesic-light-cone (GLC) coordinates are a useful tool to analyse light propagation and observations in cosmological models. In this article, we propose a detailed, pedagogical, and rigorous introduction to this coordinate system, explore its gauge degrees of freedom, and emphasize its interest when geometric optics is at stake. We then apply the GLC formalism to the homogeneous and anisotropic Bianchi I cosmology. More than a simple illustration, this application (i) allows us to show that the Weinberg conjecture according to which gravitational lensing does not affect the proper area of constant-redshift surfaces is significantly violated in a globally anisotropic universe; and (ii) offers a glimpse into new ways to constrain cosmic isotropy from the Hubble diagram.

On the scale of the light beams subtended by small sources, e.g. supernovae, matter cannot be accurately described as a fluid, which questions the applicability of standard cosmic lensing to those cases. In this article, we propose a new formalism to deal with small-scale lensing as a diffusion process: the Sachs and Jacobi equations governing the propagation of narrow light beams are treated as Langevin equations. We derive the associated Fokker-Planck-Kolmogorov equations, and use them to deduce general analytical results on the mean and dispersion of the angular distance. This formalism is applied to random Einstein-Straus Swiss-cheese models, allowing us to: (1) show an explicit example of the involved calculations; (2) check the validity of the method against both ray-tracing simulations and direct numerical integrations of the Langevin equation. As a byproduct, we obtain a post-Kantowski-Dyer-Roeder approximation, accounting for the effect of tidal distortions on the angular distance, in excellent agreement with numerical results. Besides, the dispersion of the angular distance is correctly reproduced in some regimes.

The standard model of cosmology is based on the hypothesis that the Universe is spatially homogeneous and isotropic. When interpreting most observations, this cosmological principle is applied stricto sensu: the light emitted by distant sources is assumed to propagate through a Friedmann-Lema\^itre spacetime. The main goal of the present thesis was to evaluate how reliable this assumption is, especially when small scales are at stake. After having reviewed the laws of geometric optics in curved spacetime, and the standard interpretation of cosmological observables, the dissertation reports a comprehensive analysis of light propagation in Swiss-cheese models, designed to capture the clumpy character of the Universe. The resulting impact on the interpretation of the Hubble diagram is quantified, and shown to be relatively small, thanks to the cosmological constant. When applied to current supernova data, the associated corrections tend however to improve the agreement between the cosmological parameters inferred from the Hubble diagram and from the cosmic microwave background. This is a hint that the effect of small-scale structures on light propagation may become non-negligible in the era of precision cosmology. This motivated the development of a new theoretical framework, based on stochastic processes, which aims at describing small-scale gravitational lensing with a better accuracy. Regarding the isotropy side of the cosmological principle, this dissertation addresses, on the one hand, the potential effect of a large-scale anisotropy on light propagation, by solving all the equations of geometric optics in the Bianchi I spacetime. On the other hand, possible sources of such an anisotropy, namely scalar-vector models for inflation or dark energy, are analysed. Most of them turn out to be excluded as physically viable theories.

This article proposes a comprehensive analysis of light propagation in an anisotropic and spatially homogeneous Bianchi I universe. After recalling that null geodesics are easily determined in such a spacetime, we derive the expressions of the redshift and direction drifts of light sources; by solving analytically the Sachs equation, we then obtain an explicit expression of the Jacobi matrix describing the propagation of narrow light beams. As a byproduct, we recover the old formula by Saunders for the angular diameter distance in a Bianchi I spacetime, but our derivation goes further since it also provides the optical shear and rotation. These results pave the way to the analysis of both supernovae data and weak lensing by the large-scale structure in Bianchi universes.

Various extensions of standard inflationary models have been proposed recently by adding vector fields. Because they are generally motivated by large-scale anomalies, and the possibility of statistical anisotropy of primordial fluctuations, such models require to introduce non-standard couplings between vector fields on the one hand, and either gravity or scalar fields on the other hand. In this article, we study models involving a vector field coupled to a scalar field. We derive restrictive necessary conditions for these models to be both stable (Hamiltonian bounded by below) and causal (hyperbolic equations of motion).

The recent analysis of the Planck results reveals a tension between the best fits for ({\Omega}m0, H0) derived from the cosmic microwave background or baryonic acoustic oscillations on the one hand, and the Hubble diagram on the other hand. These observations probe the universe on very different scales since they involve light beams of very different angular sizes; hence the tension between them may indicate that they should not be interpreted the same way. More precisely, this Letter questions the accuracy of using only the (perturbed) Friedmann-Lema\^itre geometry to interpret all the cosmological observations, regardless of their angular or spatial resolution. We show that using an inhomogeneous "Swiss-cheese" model to interpret the Hubble diagram allows to reconcile it with the Planck results. Such an approach does not require us to invoke new physics nor to violate the Copernican principle.

In the standard cosmological framework, the Hubble diagram is interpreted by assuming that the light emitted by standard candles propagates in a spatially homogeneous and isotropic spacetime. However, the light from "point sources"--such as supernovae--probes the Universe on scales where the homogeneity principle is no longer valid. Inhomogeneities are expected to induce a bias and a dispersion of the Hubble diagram. This is investigated by considering a Swiss-cheese cosmological model, which (1) is an exact solution of the Einstein field equations, (2) is strongly inhomogeneous on small scales, but (3) has the same expansion history as a strictly homogeneous and isotropic universe. By simulating Hubble diagrams in such models, we quantify the influence of inhomogeneities on the measurement of the cosmological parameters. Though significant in general, the effects reduce drastically for a universe dominated by the cosmological constant.

We have converted the former solar electrical plant THEMIS (French Pyrenees) into an atmospheric Cherenkov detector called CELESTE, which records gamma rays above 30 GeV (7E24 Hz). Here we present the first sub-100 GeV detection by a ground based telescope of a gamma ray source, the Crab nebula, in the energy region between satellite measurements and imaging atmospheric Cherenkov telescopes. At our analysis threshold energy of 60 +/- 20 GeV we measure a gamma ray rate of 6.1 +/- 0.8 per minute. Allowing for 30% systematic uncertainties and a 30% error on the energy scale yields an integral gamma ray flux of I(E>60 GeV) = 6.2^{+5.3}_{-2.3} E-6 photons m^-2 s^-1. The analysis methods used to obtain the gamma ray signal from the raw data are detailed. In addition, we determine the upper limit for pulsed emission to be <12% of the Crab flux at the 99% confidence level, in the same energy range. Our result indicates that if the power law observed by EGRET is attenuated by a cutoff of form e^{-E/E_0} then E_0 < 26 GeV. This is the lowest energy probed by a Cherenkov detector and leaves only a narrow range unexplored beyond the energy range studied by EGRET.

The gamma-ray emission above 250 GeV from the BL Lac object Markarian 421 was observed by the CAT Cherenkov imaging telescope between December, 1996, and June, 2000. In 1998, the source produced a series of small flares, making it the second extragalactic source detected by CAT. The time-averaged differential spectrum has been measured from 0.3 to 5 TeV, which is well fitted with a power law with an index of -2.88+-0.12(stat)+-0.06(syst). In 2000, the source showed an unprecedented activity, with variability time-scales as short as one hour, as for instance observed during the night between 4 and 5 February. The 2000 time-averaged spectrum measured is compatible with that of 1998, but some indication of a spectral curvature is found between 0.3 and 5 TeV. The possibility of TeV spectral hardening during flares is also discussed, and the results are compared to those obtained on the other TeV BL Lac, Markarian 501.

The BL Lac object Mrk 501 went into a very high state of activity during 1997, both in VHE gamma-rays and X-rays. We present here results from observations at energies above 250 GeV carried out between March and October 1997 with the CAT Cerenkov imaging Telescope. The average differential spectrum between 30 GeV and 13 TeV shows significant curvature and is well represented by phi_0 * E_TeV^{-(alpha + beta*log10(E_TeV))}, with: phi_0 = 5.19 +/- 0.13 {stat} +/- 0.12 {sys-MC} +1.66/-1.04 {sys-atm} * 10^-11 /cm^2/s/TeV alpha = 2.24 +/- 0.04 {stat} +/- 0.05 {sys} beta = 0.50 +/- 0.07 {stat} (negligible systematics). The TeV spectral energy distribution of Mrk 501 clearly peaks in the range 500 GeV-1 TeV. Investigation of spectral variations shows a significant hardness-intensity correlation with no measurable effect on the curvature. This can be described as an increase of the peak TeV emission energy with intensity. Simultaneous and quasi-simultaneous CAT VHE gamma-ray and BeppoSAX hard X-ray detections for the highest recorded flare on 16th April and for lower-activity states of the same period show correlated variability with a higher luminosity in X-rays than in gamma-rays. The observed spectral energy distribution and the correlated variability between X-rays and gamma-rays, both in amplitude and in hardening of spectra, favour a two-component emission scheme where the low and high energy components are attributed to synchrotron and inverse Compton (IC) radiation, respectively.

A new method of shower-image analysis is presented which appears very powerful as applied to those Cherenkov Imaging Telescopes with very high definition imaging capability. It provides hadron rejection on the basis of a single cut on the image shape, and simultaneously determines the energy of the electromagnetic shower and the position of the shower axis with respect to the detector. The source location is also reconstructed for each individual gamma-ray shower, even with one single telescope, so for a point source the hadron rejection can be further improved. As an example, this new method is applied to data from the CAT (Cherenkov Array at Themis) imaging telescope, which has been operational since Autumn, 1996.

The CAT (Cherenkov Array at Themis) imaging telescope, equipped with a very-high-definition camera (546 fast phototubes with 0.12 degrees spacing surrounded by 54 larger tubes in two guard rings) started operation in Autumn 1996 on the site of the former solar plant Themis (France). Using the atmospheric Cherenkov technique, it detects and identifies very high energy gamma-rays in the range 250 GeV to a few tens of TeV. The instrument, which has detected three sources (Crab nebula, Mrk 421 and Mrk 501), is described in detail.

The CELESTE experiment will be an Atmospheric Cherenkov detector designed to bridge the gap in energy sensitivity between current satellite and ground-based gamma-ray telescopes, 20 to 300 GeV. We present test results made at the former solar power plant, Themis, in the French Pyrenees. The tests confirm the viability of using a central tower heliostat array for Cherenkov wavefront sampling.

The CAT imaging telesope, which uses the Atmospheric Cherenkov Technique to observe TeV gamma-ray sources, has been operating since September 1996. Located in southern France, it features a fine grained camera consisting of 546 PMT with a pixel size of 0.12 degrees. Strong gamma-ray signals from Mkn 501 above 200 GeV were detected during observations made in April 1997. The high signal/noise ratio for these data enables direct characterization of the gamma-ray images as recorded in the camera. These results are used to validate the Monte-Carlo simulations of the response of the telescope to gamma-ray showers and to improve its performances. Observational results on Mkn 180 and the Crab nebula spectrum are presented.

The CAT Imaging Telescope has observed the BL Lac object Markarian 501 between March and August 1997. We report here on the variability over this time including several large flares. We present also preliminary spectra for all these data, for the low emission state, and for the largest flare.

During remarkable flaring activity in 1997 Markarian 501 was the brightest source in the sky at TeV energies, outshining the Crab Nebula by a factor of up to 10. Periods of flaring activity each lasting a few days were observed simultaneously by several gamma ray telescopes. Contemporaneous multiwavelength observations in April 1997 show that a substantial fraction of the total AGN power is at TeV energies indicating that high energy processes dominate the energetics of this object. Rapid variability, on time scales of less than a day, and a flat spectrum extending up to at least 10 TeV characterize this object. Results of 1997 observations by 6 telescopes are summarized and some of the implications of these results are discussed.

The VHE gamma-ray imaging telescope CAT started taking data in October 1996. Located at the Themis solar site in southern France, it features a 17.7 m^2 Davies-Cotton mirror equipped with 600 PMT camera at the focal plane. The mechanics and optics, the PMTs and the electronics are presented. The performance, based on the first 7 months of operation, is discribed.

The CAT imaging telescope on the site on the former solar plant Themis has been observing gamma-rays from Mrk501 above 220 GeV in March and April 1997. This source is shown to be highly variable and the light curve is presented. The detected gamma-ray rate for the most intense flare is in excess of 10 per minute.

The CAT imaging telescope, recently built on the site of the former solar plant Themis (French Pyrenees), observed gamma-rays from the Crab nebula from October 1996 to March 1997. This steady source, often considered as the standard candle of very-high-energy gamma-ray astronomy, is used as a test-beam to probe the performances of the new telescope, particularly its energy threshold (220 GeV at 20 degrees zenith angle) and the stability of its response. Due to the fine-grain camera, an accurate analysis of the longitudinal profiles of shower images is performed, yielding the source position in two dimensions for each individual shower.

There is a growing interest in the possibility that dark matter could be formed of weakly interacting particles with a mass in the 100 GeV - 2 TeV range, and supersymmetric particles are favorite candidates. If they constitute the dark halo of our Galaxy, their mutual annihilations produce energetic gamma rays that could be detected using existing atmospheric \u{C}erenkov techniques.