This article describes the dynamics of a charge particle in a electromagnetic field in presence of a scalar fifth force. Focusing to the fifth force that would be induced by a chameleon field, the profile of which can be designed properly in the laboratory, it draws its physical effects on the cyclotron motion of a particle in a static and uniform magnetic field. The fifth force induces a drift of the trajectory that is estimated analytically and compared to numerical computations for profiles motivated by the ones of a chameleon field within two nested cylinders. The magnitude of the effect and the detectability of this drift are discussed to argue that this may offer a new experimental design to test small fifth force in the laboratory. More important, at the macroscopic level it induces a current that can in principle also be measured, and would even allow one to access the transverse profile of the scalar field within the cavity. In both cases, aligning the magnetic field with the local gravity field suppresses the effects of Newtonian gravity that would be several order larger than the ones of the fifth force otherwise and the Newtonian gravity of the cavity on the particle is also argued to be negligible. Given this insight, this experimental set-up, with its two effects -- on a single particle and at the macroscopic level -- may require attention to demonstrate its actual feasability in the laboratory.

This article investigates the properties of a scalar fifth force that arises in a scalar tensor-theory with a chameleon screening mechanism in the context of gravity space missions like the MICROSCOPE experiment. In such an experiment, the propagation of the chameleon field inside the nested cylinders of the experiment causes a fifth force when the cylinders are not perfectly co-axial. We propose a semi-analytic method to compute the field distribution and the induced fifth force and compare it to a full numerical simulation, in settings where the cylindrical symmetry is broken. The scaling of the fifth force with both the parameters of the model and the geometry of the experiment is discussed. We show that the fifth force is repulsive, hence adds a destabilizing stiffness that should be included in the force budget acting on the detector. This opens the way to a new method to constrain a scalar fifth force in screened models.

We present the scientific motivation for future space tests of the equivalence principle, and in particular the universality of free fall, at the $10^{-17}$ level or better. Two possible mission scenarios, one based on quantum technologies, the other on electrostatic accelerometers, that could reach that goal are briefly discussed.

The Weak Equivalence Principle (WEP), stating that two bodies of different compositions and/or mass fall at the same rate in a gravitational field (universality of free fall), is at the very foundation of General Relativity. The MICROSCOPE mission aims to test its validity to a precision of $10^{-15}$, two orders of magnitude better than current on-ground tests, by using two masses of different compositions (titanium and platinum alloys) on a quasi-circular trajectory around the Earth. This is realised by measuring the accelerations inferred from the forces required to maintain the two masses exactly in the same orbit. Any significant difference between the measured accelerations, occurring at a defined frequency, would correspond to the detection of a violation of the WEP, or to the discovery of a tiny new type of force added to gravity. MICROSCOPE's first results show no hint for such a difference, expressed in terms of E\"otv\"os parameter $\delta(Ti,Pt)=[-1\pm{}9{\rm (stat)}\pm{}9{\rm (syst)}] \times{}10^{-15}$ (both 1$\sigma$ uncertainties) for a titanium and platinum pair of materials. This result was obtained on a session with 120 orbital revolutions representing 7\% of the current available data acquired during the whole mission. The quadratic combination of 1$\sigma$ uncertainties leads to a current limit on $\delta$ of about $1.3\times{}10^{-14}$.

We speculate on the development and availability of new innovative propulsion techniques in the 2040s, that will allow us to fly a spacecraft outside the Solar System (at 150 AU and more) in a reasonable amount of time, in order to directly probe our (gravitational) Solar System neighborhood and answer pressing questions regarding the dark sector (dark energy and dark matter). We identify two closely related main science goals, as well as secondary objectives that could be fulfilled by a mission dedicated to probing the local dark sector: (i) begin the exploration of gravitation's low-acceleration regime with a man-made spacecraft and (ii) improve our knowledge of the local dark matter and baryon densities. Those questions can be answered by directly measuring the gravitational potential with an atomic clock on-board a spacecraft on an outbound Solar System orbit, and by comparing the spacecraft's trajectory with that predicted by General Relativity through the combination of ranging data and the in-situ measurement (and correction) of non-gravitational accelerations with an on-board accelerometer. Despite a wealth of new experiments getting online in the near future, that will bring new knowledge about the dark sector, it is very unlikely that those science questions will be closed in the next two decades. More importantly, it is likely that it will be even more urgent than currently to answer them. Tracking a spacecraft carrying a clock and an accelerometer as it leaves the Solar System may well be the easiest and fastest way to directly probe our dark environment.

This article investigates the profile of the scalar field of a scalar-tensor theory subject to the chameleon mechanism in the context of gravity space missions like the MICROSCOPE experiment. It analyses the experimental situations for models with an inverse power law potential that can in principle induce a fifth force inside the satellite, hence either be detected or constrained. As the mass of the scalar field depends on the local matter density, the screening of the scalar field depends crucially on both the parameters of the theory (potential and non-minimal coupling to matter) and on the geometry of the satellite. We calculate the profile of the scalar field in 1-, 2- and 3-dimensional satellite configurations without relying on the thick or thin shell approximations for the scalar field. In particular we consider the typical geometry with nested cylinders which is close to the MICROSCOPE design. In this case we evaluate the corresponding fifth force on a test body inside the satellite. This analysis clarifies previous claims on the detectability of the chameleon force by space-borne experiments.

We introduce one- and two-dimensional `exponential shapelets': orthonormal basis functions that efficiently model isolated features in data. They are built from eigenfunctions of the quantum mechanical hydrogen atom, and inherit mathematics with elegant properties under Fourier transform, and hence (de)convolution. For a wide variety of data, exponential shapelets compress information better than Gauss-Hermite/Gauss-Laguerre (`shapelet') decomposition, and generalise previous attempts that were limited to 1D or circularly symmetric basis functions. We discuss example applications in astronomy, fundamental physics and space geodesy.

The tests of the deviations from Newton's or Einstein's gravity in the Earth neighbourhood are tied to our knowledge of the shape and mass distribution of our planet. On the one hand estimators of these "modified" theories of gravity may be explicitly Earth-model-dependent whilst on the other hand the Earth gravitational field would act as a systematic error. We revisit deviations from Newtonian gravity described by a Yukawa interaction that can arise from the existence of a finite range fifth force. We show that the standard multipolar expansion of the Earth gravitational potential can be generalised. In particular, the multipolar coefficients depend on the distance to the centre of the Earth and are therefore not universal to the Earth system anymore. This offers new ways of constraining such Yukawa interactions and demonstrates explicitly the limits of the Newton-based interpretation of geodesy experiments. In turn, limitations from geodesy data restrict the possibility of testing gravity in space. The gravitational acceleration is described in terms of spin-weighted spherical harmonics allowing us to obtain the perturbing force entering the Lagrange-Gauss secular equations. This is then used to discuss the correlation between geodesy and modified gravity experiments and the possibility to break their degeneracy. Finally we show that, given the existing constraints, a Yukawa fifth force is expected to be sub-dominant in satellite dynamics and space geodesy experiments, as long as they are performed at altitudes greater than a few hundred kilometres. Gravity surveys will have to gain at least two orders of magnitude in instrumental precision before satellite geodesy could be used to improve the current constraints on modified gravity.

The existence of a light or massive scalar field with a coupling to matter weaker than gravitational strength is a possible source of violation of the weak equivalence principle. We use the first results on the E\"otv\"os parameter by the MICROSCOPE experiment to set new constraints on such scalar fields. For a massive scalar field of mass smaller than $10^{-12}$ eV (i.e. range larger than a few $10^5$ m) we improve existing constraints by one order of magnitude to $|\alpha|<10^{-11}$ if the scalar field couples to the baryon number and to $|\alpha|<10^{-12}$ if the scalar field couples to the difference between the baryon and the lepton numbers. We also consider a model describing the coupling of a generic dilaton to the standard matter fields with five parameters, for a light field: we find that for masses smaller than $10^{-12}$eV, the constraints on the dilaton coupling parameters are improved by one order of magnitude compared to previous equivalence principle tests.

According to the Weak Equivalence Principle, all bodies should fall at the same rate in a gravitational field. The MICROSCOPE satellite, launched in April 2016, aims to test its validity at the $10^{-15}$ precision level, by measuring the force required to maintain two test masses (of titanium and platinum alloys) exactly in the same orbit. A non-vanishing result would correspond to a violation of the Equivalence Principle, or to the discovery of a new long-range force. Analysis of the first data gives $\delta\rm{(Ti,Pt)}= [-1 \pm 9 (\mathrm{stat}) \pm 9 (\mathrm{syst})] \times 10^{-15}$ (1$\sigma$ statistical uncertainty) for the titanium-platinum E\"otv\"os parameter characterizing the relative difference in their free-fall accelerations.

The MICROSCOPE space mission, launched on April 25, 2016, aims to test the weak equivalence principle (WEP) with a 10^-15 precision. To reach this performance requires an accurate and robust data analysis method, especially since the possible WEP violation signal will be dominated by a strongly colored noise. An important complication is brought by the fact that some values will be missing -therefore, the measured time series will not be strictly regularly sampled. Those missing values induce a spectral leakage that significantly increases the noise in Fourier space, where the WEP violation signal is looked for, thereby complicating scientific returns. Recently, we developed an inpainting algorithm to correct the MICROSCOPE data for missing values. This code has been integrated in the official MICROSCOPE data processing pipeline because it enables us to significantly measure an equivalence principle violation (EPV) signal in a model-independent way, in the inertial satellite configuration. In this work, we present several improvements to the method that may allow us now to reach the MICROSCOPE requirements for both inertial and spin satellite configurations. The main improvement has been obtained using a prior on the power spectrum of the colored-noise that can be directly derived from the incomplete data. We show that after reconstructing missing values with this new algorithm, a least-squares fit may allow us to significantly measure an EPV signal with a 0.96x10^-15 precision in the inertial mode and 1.2x10^-15 precision in the spin mode. Although, the inpainting method presented in this paper has been optimized to the MICROSCOPE data, it remains sufficiently general to be used in the general context of missing data in time series dominated by an unknown colored-noise. The improved inpainting software, called ICON, is freely available at http://www.cosmostat.org/software/icon.

Missing data are a common problem in experimental and observational physics. They can be caused by various sources, either an instrument's saturation, or a contamination from an external event, or a data loss. In particular, they can have a disastrous effect when one is seeking to characterize a colored-noise-dominated signal in Fourier space, since they create a spectral leakage that can artificially increase the noise. It is therefore important to either take them into account or to correct for them prior to e.g. a Least-Square fit of the signal to be characterized. In this paper, we present an application of the {\it inpainting} algorithm to mock MICROSCOPE data; {\it inpainting} is based on a sparsity assumption, and has already been used in various astrophysical contexts; MICROSCOPE is a French Space Agency mission, whose launch is expected in 2016, that aims to test the Weak Equivalence Principle down to the $10^{-15}$ level. We then explore the {\it inpainting} dependence on the number of gaps and the total fraction of missing values. We show that, in a worst-case scenario, after reconstructing missing values with {\it inpainting}, a Least-Square fit may allow us to significantly measure a $1.1\times10^{-15}$ Equivalence Principle violation signal, which is sufficiently close to the MICROSCOPE requirements to implement {\it inpainting} in the official MICROSCOPE data processing and analysis pipeline. Together with the previously published KARMA method, {\it inpainting} will then allow us to independently characterize and cross-check an Equivalence Principle violation signal detection down to the $10^{-15}$ level.

The analysis of physical measurements often copes with highly correlated noises and interruptions caused by outliers, saturation events or transmission losses. We assess the impact of missing data on the performance of linear regression analysis involving the fit of modeled or measured time series. We show that data gaps can significantly alter the precision of the regression parameter estimation in the presence of colored noise, due to the frequency leakage of the noise power. We present a regression method which cancels this effect and estimates the parameters of interest with a precision comparable to the complete data case, even if the noise power spectral density (PSD) is not known a priori. The method is based on an autoregressive (AR) fit of the noise, which allows us to build an approximate generalized least squares estimator approaching the minimal variance bound. The method, which can be applied to any similar data processing, is tested on simulated measurements of the MICROSCOPE space mission, whose goal is to test the Weak Equivalence Principle (WEP) with a precision of $10^{-15}$. In this particular context the signal of interest is the WEP violation signal expected to be found around a well defined frequency. We test our method with different gap patterns and noise of known PSD and find that the results agree with the mission requirements, decreasing the uncertainty by a factor 60 with respect to ordinary least squares methods. We show that it also provides a test of significance to assess the uncertainty of the measurement.

MICROSCOPE is a French Space Agency mission that aims to test the Weak Equivalence Principle in space down to an accuracy of $10^{-15}$. This is two orders of magnitude better than the current constraints, which will allow us to test General Relativity as well as theories beyond General Relativity which predict a possible Weak Equivalence Principle violation below $10^{-13}$. In this communication, we describe the MICROSCOPE mission, its measurement principle and instrument, and we give an update on its status. After a successful instrument's commissioning, MICROSCOPE is on track for on-schedule launch, expected in 2016.

This work attains a threefold objective: first, we derived the richness-mass scaling in the local Universe from data of 53 clusters with individual measurements of mass. We found a 0.46+-0.12 slope and a 0.25+-0.03 dex scatter measuring richness with a previously developed method. Second, we showed on a real sample of 250 0.06<z<0.9 clusters, most of which are at z<0.3, with spectroscopic redshift that the colour of the red sequence allows us to measure the clusters' redshift to better than Delta z=0.02. Third, we computed the predicted prior of the richness-mass scaling to forecast the capabilities of future wide-field-area surveys of galaxy clusters to constrain cosmological parameters. We computed the uncertainty and the covariance matrix of the (evolving) richness-mass scaling of a PanStarrs 1+Euclid-like survey accounting for a large suite of sources of errors. We find that the richness-mass scaling parameters, which are the input ingredients of cosmological forecasts using cluster counts, can be determined 10^5 times better than estimated in previous works that did not use weak-lensing mass estimates. The better knowledge of the scaling parameters likely has a strong impact on the relative importance of the different probes used to constrain cosmological parameters. Richness-mass scaling parameters were recovered, but only if the cluster mass function and the weak-lensing redshift-dependent selection function were accounted for in the fitting of the mass-richness scaling. This emphasizes the limitations of often adopted simplifying assumptions, such as having a mass-complete redshift-independent sample. The fitting code used for computing the predicted prior, including the treatment of the mass function and of the weak-lensing selection function, is provided in the appendix. [Abridged]

Simulated wide-field images are becoming an important part of observational astronomy, either to prepare for new surveys or to test measurement methods. In order to efficiently explore vast parameter spaces, the computational speed of simulation codes is a central requirement to their implementation. We introduce the Ultra Fast Image Generator (UFig) which aims to bring wide-field imaging simulations to the current limits of computational capabilities. We achieve this goal through: (1) models of galaxies, stars and observational conditions, which, while simple, capture the key features necessary for realistic simulations, and (2) state-of-the-art computational and implementation optimizations. We present the performances of UFig and show that it is faster than existing public simulation codes by several orders of magnitude. It allows us to produce images more quickly than SExtractor needs to analyze them. For instance, it can simulate a typical 0.25 deg^2 Subaru SuprimeCam image (10k x 8k pixels) with a 5-sigma limiting magnitude of R=26 in 30 seconds on a laptop, yielding an average simulation time for a galaxy of 30 microseconds. This code is complementary to end-to-end simulation codes and can be used as a fast, central component of observational methods relying on simulations.

Point Spread Function (PSF) modeling is a central part of any astronomy data analysis relying on measuring the shapes of objects. It is especially crucial for weak gravitational lensing, in order to beat down systematics and allow one to reach the full potential of weak lensing in measuring dark energy. A PSF modeling pipeline is made of two main steps: the first one is to assess its shape on stars, and the second is to interpolate it at any desired position (usually galaxies). We focus on the second part, and compare different interpolation schemes, including polynomial interpolation, radial basis functions, Delaunay triangulation and Kriging. For that purpose, we develop simulations of PSF fields, in which stars are built from a set of basis functions defined from a Principal Components Analysis of a real ground-based image. We find that Kriging gives the most reliable interpolation, significantly better than the traditionally used polynomial interpolation. We also note that although a Kriging interpolation on individual images is enough to control systematics at the level necessary for current weak lensing surveys, more elaborate techniques will have to be developed to reach future ambitious surveys' requirements.

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