Search

With the rise of simulation-based inference (SBI) methods, simulations need to be fast as well as realistic. $\texttt{UFig v1}$ is a public Python package that generates simulated astronomical images with exceptional speed - taking approximately the same time as source extraction. This makes it particularly well-suited for simulation-based inference (SBI) methods where computational efficiency is crucial. To render an image, $\texttt{UFig}$ requires a galaxy catalog and a description of the point spread function (PSF). It can also add background noise, sample stars using the Besan\c{c}on model of the Milky Way, and run $\texttt{SExtractor}$ to extract sources from the rendered image. The extracted sources can be matched to the intrinsic catalog, flagged based on $\texttt{SExtractor}$ output and survey masks, and emulators can be used to bypass the image simulation and extraction steps. A first version of $\texttt{UFig}$ was presented in Berg\'e et al. (2013) and the software has since been used and further developed in a variety of forward modelling applications., Comment: 8 pages, 4 figures, submitted to JOSS

Scalar-tensor theories with screening mechanisms come with non-linearities that make it difficult to study setups of complex geometry without resorting to numerical simulations. In this article, we use the $\textit{femtoscope}$ code that we introduced in a previous work in order to compute the fifth force arising in the chameleon model in the Earth orbit. We go beyond published works by introducing a departure from spherical symmetry $\unicode{x2014}$ embodied by a mountain on an otherwise spherical Earth $\unicode{x2014}$ as well as by implementing several atmospheric models, and quantify their combined effect on the chameleon field. Building on the numerical results thus obtained, we address the question of the detectability of a putative chameleon fifth force by means of space geodesy techniques and, for the first time, quantitatively assess the back-reaction created by the screening of a satellite itself. We find that although the fifth force has a supposedly measurable effect on the dynamics of an orbiting spacecraft, the imprecise knowledge of the mass distribution inside the Earth greatly curtails the constraining power of such space missions. Finally, we show how this degeneracy can be lifted when several measurements are performed at different altitudes., Comment: 36 pages, 20 figures. Accepted for publication in PRD

The space mission MICROSCOPE dedicated to the test of the Equivalence Principle (EP) operated from April 25, 2016 until the deactivation of the satellite on October 16, 2018. In this analysis we compare the free-fall accelerations ($a_{\rm A}$ and $a_{\rm B}$) of two test masses in terms of the E\"otv\"os parameter $\eta({\rm{A, B}}) = 2 \frac{a_{\rm A}- a_{\rm B}}{a_{\rm A}+ a_{\rm B}}$. No EP violation has been detected for two test masses, made from platinum and titanium alloys, in a sequence of 19 segments lasting from 13 to 198 hours down to the limit of the statistical error which is smaller than $10^{-14}$ for $ \eta({\rm{Ti, Pt}})$. Accumulating data from all segments leads to $\eta({\rm{Ti, Pt}}) =[-1.5\pm{}2.3{\rm (stat)}\pm{}1.5{\rm (syst)}] \times{}10^{-15}$ showing no EP violation at the level of $2.7\times{}10^{-15}$ if we combine stochastic and systematic errors quadratically. This represents an improvement of almost two orders of magnitude with respect to the previous best such test performed by the E\"ot-Wash group. The reliability of this limit has been verified by comparing the free falls of two test masses of the same composition (platinum) leading to a null E\"otv\"os parameter with a statistical uncertainty of $1.1\times{}10^{-15}$., Comment: Class. Quantum Grav. 39 204009

The MICROSCOPE mission was designed to test the Weak Equivalence Principle (WEP), stating the equality between the inertial and the gravitational masses, with a precision of $10^{-15}$ in terms of the E\"otv\"os ratio $\eta$. Its experimental test consisted of comparing the accelerations undergone by two collocated test masses of different compositions as they orbited the Earth, by measuring the electrostatic forces required to keep them in equilibrium. This was done with ultra-sensitive differential electrostatic accelerometers onboard a drag-free satellite. The mission lasted two and a half years, cumulating five-months-worth of science free-fall data, two thirds with a pair of test masses of different compositions -- Titanium and Platinum alloys -- and the last third with a reference pair of test masses of the same composition -- Platinum. We summarize the data analysis, with an emphasis on the characterization of the systematic uncertainties due to thermal instabilities and on the correction of short-lived events which could mimic a WEP violation signal. We found no violation of the WEP, with the E\"otv\"os parameter of the Titanium and Platinum pair constrained to $\eta({\rm Ti, Pt})~=~ [-1.5 \pm 2.3~{\rm (stat)} \pm 1.5~{\rm (syst)}]~\times 10^{-15}$ at $1\sigma$ in statistical errors., Comment: Phys. Rev. Lett. 129, 121102

A large class of scalar-tensor theories of gravity exhibit a screening mechanism that dynamically suppresses fifth forces in the Solar system and local laboratory experiments. Technically, at the scalar field equation level, this usually translates into nonlinearities which strongly limit the scope of analytical approaches. This article presents $femtoscope$ $-$ a Python numerical tool based on the Finite Element Method (FEM) and Newton method for solving Klein-Gordon-like equations that arise in particular in the symmetron or chameleon models. Regarding the latter, the scalar field behavior is generally only known infinitely far away from the its sources. We thus investigate existing and new FEM-based techniques for dealing with asymptotic boundary conditions on finite-memory computers, whose convergence are assessed. Finally, $femtoscope$ is showcased with a study of the chameleon fifth force in Earth orbit., Comment: 35 pages, 20 figures, 2 tables. Accepted for publication in PRD

The objective of the proposed MAQRO mission is to harness space for achieving long free-fall times, extreme vacuum, nano-gravity, and cryogenic temperatures to test the foundations of physics in macroscopic quantum experiments. This will result in the development of novel quantum sensors and a means to probe the foundations of quantum physics at the interface with gravity. Earlier studies showed that the proposal is feasible but that several critical challenges remain, and key technologies need to be developed. These new technologies will open up the potential for achieving additional science objectives. The proposed research campaign aims to advance the state of the art and to perform the first macroscopic quantum experiments in space. Experiments on the ground, in micro-gravity, and in space will drive the proposed research campaign during the current decade to enable the implementation of MAQRO within the subsequent decade., Comment: 9 pages, 2 figures, submitted as a Research Campaign Whitepaper for the BPS2023 Decadal Survey

Testing the Weak Equivalence Principle (WEP) to a precision of $10^{-15}$ requires a quantity of data that give enough confidence on the final result: ideally, the longer the measurement the better the rejection of thestatistical noise. The science sessions had a duration of 120 orbits maximum and were regularly repeated and spaced out to accommodate operational constraints but also in order to repeat the experiment in different conditions and to allow time to calibrate the instrument. Several science sessions were performed over the 2.5 year duration of the experiment. This paper aims to describe how the data have been produced on the basis of a mission scenario and a data flow process, driven by a tradeoff between the science objectives and the operational constraints. The mission was led by the Centre National d'Etudes Spatiales (CNES) which provided the satellite, the launch and the ground operations. The ground segment was distributed between CNES and Office National d'Etudes et de Recherches A\'erospatiales (ONERA). CNES provided the raw data through the Centre d'Expertise de Compensation de Tra\^{i}n\'{e}e (CECT: Drag-free expertise centre). The science was led by the Observatoire de la C\^ote d{'}Azur (OCA) and ONERA was in charge of the data process. The latter also provided the instrument and the Science Mission Centre of MICROSCOPE (CMSM)., Comment: To be released in special MICROSCOPE edition of CQG

We summarize the discussions at a virtual Community Workshop on Cold Atoms in Space concerning the status of cold atom technologies, the prospective scientific and societal opportunities offered by their deployment in space, and the developments needed before cold atoms could be operated in space. The cold atom technologies discussed include atomic clocks, quantum gravimeters and accelerometers, and atom interferometers. Prospective applications include metrology, geodesy and measurement of terrestrial mass change due to, e.g., climate change, and fundamental science experiments such as tests of the equivalence principle, searches for dark matter, measurements of gravitational waves and tests of quantum mechanics. We review the current status of cold atom technologies and outline the requirements for their space qualification, including the development paths and the corresponding technical milestones, and identifying possible pathfinder missions to pave the way for missions to exploit the full potential of cold atoms in space. Finally, we present a first draft of a possible road-map for achieving these goals, that we propose for discussion by the interested cold atom, Earth Observation, fundamental physics and other prospective scientific user communities, together with ESA and national space and research funding agencies., Comment: Summary of the Community Workshop on Cold Atoms in Space and corresponding Road-map: https://indico.cern.ch/event/1064855/

The MICROSCOPE mission aims to test the Weak Equivalence Principle (WEP) in orbit with an unprecedented precision of 10$^{-15}$ on the E\"otv\"os parameter thanks to electrostatic accelerometers on board a drag-free micro-satellite. The precision of the test is determined by statistical errors, due to the environment and instrument noises, and by systematic errors to which this paper is devoted. Systematic error sources can be divided into three categories: external perturbations, such as the residual atmospheric drag or the gravity gradient at the satellite altitude, perturbations linked to the satellite design, such as thermal or magnetic perturbations, and perturbations from the instrument internal sources. Each systematic error is evaluated or bounded in order to set a reliable upper bound on the WEP parameter estimation uncertainty., Comment: To be released in CQG MICROSCOPE Special Edition

The MICROSCOPE experiment was designed to test the weak equivalence principle in space, by comparing the low-frequency dynamics of cylindrical "free-falling" test masses controlled by electrostatic forces. We use data taken during technical sessions aimed at estimating the electrostatic stiffness of MICROSCOPE's sensors to constrain a short-range Yukawa deviation from Newtonian gravity. We take advantage of the fact that in the limit of small displacements, the gravitational interaction (both Newtonian and Yukawa-like) between nested cylinders is linear, and thus simply characterised by a stiffness. By measuring the total stiffness of the forces acting on a test mass as it moves, and comparing it with the theoretical electrostatic stiffness (expected to dominate), it is a priori possible to infer constraints on the Yukawa potential parameters. However, we find that measurement uncertainties are dominated by the gold wires used to control the electric charge of the test masses, though their related stiffness is indeed smaller than the expected electrostatic stiffness. Moreover, we find a non-zero unaccounted for stiffness that depends on the instrument's electric configuration, hinting at the presence of patch-field effects. Added to significant uncertainties on the electrostatic model, they only allow for poor constraints on the Yukawa potential. This is not surprising, as MICROSCOPE was not designed for this measurement, but this analysis is the first step to new experimental searches for non-Newtonian gravity in space., Comment: Accepted for publication in CQG's MICROSCOPE special issue

This article is dedicated to the use the MICROSCOPE mission's data to test chameleon theory of gravity. We take advantage of the technical sessions aimed to characterize the electrostatic stiffness of MICROSCOPE's instrument intrinsic to its capacitive measurement system. Any discrepancy between the expected and measured stiffness may result from unaccounted-for contributors, i.e. extra-forces. This work considers the case of chameleon gravity as a possible contributor. It was previously shown that in situations similar to these measurement sessions, a chameleon fifth force appears and acts as a stiffness for small displacements. The magnitude of this new component of the stiffness is computed over the chameleon's parameter space. It allows us to derive constraints by excluding any force inconsistent with the MICROSCOPE data. As expected --since MICROSCOPE was not designed for the purpose of such an analysis--, these new bounds are not competitive with state-of-the-art constraints, but they could be improved by a better estimation of all effects at play in these sessions. Hence our work illustrates this novel technique as a new way of constraining fifth forces., Comment: Submitted to PRD

Dedicated accelerometers have been developed for the MICROSCOPE mission taking into account the specific range of acceleration to be measured on board the satellite. Considering one micro-g and even less as the full range of the instrument, leads to a customized concept and a high performance electronics for the sensing and servo-actuations of the accelerometer test-masses. In addition to a very accurate geometrical sensor core, a high performance electronics architecture provides the measurement of the weak electrostatic forces and torques applied to the test-masses. A set of capacitive sensors delivers the position and the attitude of the test-mass with respect to a very steady gold coated cage made in silica. The voltages applied on the electrodes surrounding each test-mass are finely controlled to generate the adequate electrical field and so the electrostatic pressures on the test-mass. This field maintains the test-mass motionless with respect to the instrument structure. Digital control laws are implemented in order to enable instrument operation flexibility and a weak position sensor noise. These electronics provide both the scientific data for MICROSCOPE's test of General Relativity and the data for the satellite drag-free and attitude control system (DFACS)., Comment: To be submitted to CQG's MICROSCOPE special issue

MICROSCOPE's space test of the weak equivalence principle (WEP) is based on the minute measurement of the difference of accelerations experienced by two test masses as they orbit the Earth. A detection of a violation of the WEP would appear at a well-known frequency $f_{\rm EP}$ depending on the satellite's orbital and spinning frequencies. Consequently, the experiment was optimised to miminise systematic errors at $f_{\rm EP}$. Glitches are short-lived events visible in the test masses' measured acceleration, most likely originating in cracks of the satellite's coating. In this paper, we characterise their shape and time distribution. Although intrinsically random, their time of arrival distribution is modulated by the orbital and spinning periods. They have an impact on the WEP test that must be quantified. However, the data available prevents us from unequivocally tackling this task. We show that glitches affect the test of the WEP, up to an a priori unknown level. Discarding the perturbed data is thus the best way to reduce their effect., Comment: References updated

After performing highly sensitive acceleration measurements during two years of drag-free flight around the Earth, MICROSCOPE provided the best constraint on the Weak Equivalence Principle (WEP) to date. Beside being a technological challenge, this experiment required a specialised data analysis pipeline to look for a potential small signal buried in the noise, possibly plagued by instrumental defects, missing data and glitches. This paper describes the frequency-domain iterative least-square technique that we developed for MICROSCOPE. In particular, using numerical simulations, we prove that our estimator is unbiased and provides correct error bars. This paper therefore justifies the robustness of the WEP measurements given by MICROSCOPE., Comment: References updated

This paper focus on the description of the design and performance of the MICROSCOPE satellite and its Drag-Free and Attitude Control System (DFACS). The satellite is derived from CNES' Myriade platform family, albeit with significant upgrades dictated by the unprecedented MICROSCOPE's mission requirements. The 300kg drag-free microsatellite has completed its 2-year flight with higher-than-expected performances. Its passive thermal concept allowed for variations smaller than 1 $\mu$K at the measurement frequency $f_{\rm{EP}}$. The propulsion system provided a 6 axis continuous and very low noise thrust from zero to some hundreds of micronewtons. Finally, the performance of its DFACS (aimed at compensating the disturbing forces and torques applied to the satellite) is the finest ever achieved in low Earth orbit, with residual accelerations along the three axes are lower than $10^{-12} {\rm m/s}^2$ at $f_{\rm{EP}}$ over 8 days., Comment: References updated

The MICROSCOPE mission aimed to test the Weak Equivalence Principle (WEP) to a precision of $10^{-15}$. The WEP states that two bodies fall at the same rate on a gravitational field independently of their mass or composition. In MICROSCOPE, two masses of different compositions (titanium and platinum alloys) are placed on a quasi-circular trajectory around the Earth. They are the test-masses of a double accelerometer. The measurement of their accelerations is used to extract a potential WEP violation that would occur at a frequency defined by the motion and attitude of the satellite around the Earth. This paper details the major drivers of the mission leading to the specification of the major subsystems (satellite, ground segment, instrument, orbit...). Building upon the measurement equation, we derive the objective of the test in statistical and systematic error allocation and provide the mission's expected error budget., Comment: References updated

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., Comment: 16 pages, 21 figures

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., Comment: Submitted to Phys. Rev. D

We use data from the T-SAGE instrument on board the MICROSCOPE space mission to search for Lorentz violation in matter-gravity couplings as described by the Lorentz violating Standard-Model Extension (SME) coefficients $(\bar{a}_\text{eff})_\mu^w$, where ($\mu = T,X,Y,Z$) and ($w = e,p,n$) for the electron, proton and neutron. One of the phenomenological consequences of a non-zero value of those coefficients is that test bodies of different composition fall differently in an external gravitational field. This is similar to "standard" tests of the universality of free fall, but with a specific signature that depends on the orbital velocity and rotation of the Earth. We analyze data from five measurement sessions of MICROSCOPE spread over a year finding no evidence for such a signature, but setting constraints on linear combinations of the SME coefficients that improve on best previous results by one to two orders of magnitude. Additionally, our independent linear combinations are different from previous ones, which increases the diversity of available constraints, paving the way towards a full decorrelation of the individual coefficients.

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}$., Comment: To appear in CQG

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., Comment: White Paper submitted to ESA's Voyage 2050 call for papers. Reference updated

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., Comment: White paper submitted to the ESA Voyage 2050 long term plan. Partly derived from Altschul et al., Advances in Space Research, 2015, 55, 501 - 524, arXiv:1404.4307

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., Comment: Submitted to Phys. Rev. D

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., Comment: Accepted for publication in MNRAS. 13+4 pages

The ISLAND (Inverse Square Law And Newtonian Dynamics) Space Explorer is a new concept to test the gravitational Inverse Square Law at: (1) submillimeter scale and (2) at the largest Solar System scales (dozens of Astronomical Units --AU). The main idea is to embark a torsion pendulum at the center of gravity of a dedicated, possibly drag-free and attitude-controlled, interplanetary probe whose gravitational environment is accurately probed by, and corrected for thanks to six ultrasensitive accelerometers arranged as a cross around the torsion pendulum., Comment: Moriond Gravitation session 2017 proceeding + COSPAR meeting 2018

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., Comment: submitted to CQG

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., Comment: Typos corrected

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., Comment: Matches publlished version. Supplemental material added

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., Comment: 10 pages, 8 figures, 1 table, Accepted for publication in Physical Review D

We present a Gaussian regression method for time series with missing data and stationary residuals of unknown power spectral density (PSD). The missing data are efficiently estimated by their conditional expectation as in universal Kriging, based on the circulant approximation of the complete data covariance. After initialization with an autoregessive fit of the noise, a few iterations of estimation/reconstruction steps are performed until convergence of the regression and PSD estimates, in a way similar to the expectation-conditional-maximization algorithm. The estimation can be performed for an arbitrary PSD provided that it is sufficiently smooth. The algorithm is developed in the framework of the MICROSCOPE space mission whose goal is to test the weak equivalence principle (WEP) with a precision of $10^{-15}$. We show by numerical simulations that the developed method allows us to meet three major requirements: to maintain the targeted precision of the WEP test in spite of the loss of data, to calculate a reliable estimate of this precision and of the noise level, and finally to provide consistent and faithful reconstructed data to the scientific community., Comment: 17 pages, 6 figures

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., Comment: Accepted for publication in Physical Review D. 12 pages, 6 figures

Weak lensing by large-scale structure is a powerful technique to probe the dark components of the universe. To understand the measurement process of weak lensing and the associated systematic effects, image simulations are becoming increasingly important. For this purpose we present a first implementation of the $\textit{Monte Carlo Control Loops}$ ($\textit{MCCL}$; Refregier & Amara 2014), a coherent framework for studying systematic effects in weak lensing. It allows us to model and calibrate the shear measurement process using image simulations from the Ultra Fast Image Generator (UFig; Berge et al. 2013). We apply this framework to a subset of the data taken during the Science Verification period (SV) of the Dark Energy Survey (DES). We calibrate the UFig simulations to be statistically consistent with DES images. We then perform tolerance analyses by perturbing the simulation parameters and study their impact on the shear measurement at the one-point level. This allows us to determine the relative importance of different input parameters to the simulations. For spatially constant systematic errors and six simulation parameters, the calibration of the simulation reaches the weak lensing precision needed for the DES SV survey area. Furthermore, we find a sensitivity of the shear measurement to the intrinsic ellipticity distribution, and an interplay between the magnitude-size and the pixel value diagnostics in constraining the noise model. This work is the first application of the $\textit{MCCL}$ framework to data and shows how it can be used to methodically study the impact of systematics on the cosmic shear measurement., Comment: 14 pages, 9 Figures, submitted to ApJ

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., Comment: Proceedings of LISA Symposium X; accepted for publication in Journal of Physics: Conference Series

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., Comment: Submitted to Astronomy and Computing. 13 pages, 9 figures

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed., Comment: 26 pages and 7 figures; to appear in Astroparticle Physics

In the theory of structure formation, galaxies are biased tracers of the underlying matter density field. The statistical relation between galaxy and matter density field is commonly referred as galaxy bias. In this paper, we test the linear bias model with weak-lensing and galaxy clustering measurements in the 2 square degrees COSMOS field (Scoville et al. 2007). We estimate the bias of galaxies between redshifts z=0.2 and z=1, and over correlation scales between R=0.2 h^-1 Mpc and R=15 h^-1 Mpc. We focus on three galaxy samples, selected in flux (simultaneous cuts I_814W < 26.5 and K_s < 24), and in stellar-mass (10^9 < M_* < 10^10 h^-2 Msun and 10^10 < M^*< 10^11 h^-2 Msun). At scales R > 2 h^-1 Mpc, our measurements support a model of bias increasing with redshift. The Tinker et al. (2010) fitting function provides a good fit to the data. We find the best fit mass of the galaxy halos to be log(M_200 h^-1 Msun) = 11.7^+0.6_-1.3 and log(M_200 h^-1 Msun) = 12.4^+0.2_-2.9 respectively for the low and high stellar-mass samples. In the halo model framework, bias is scale-dependent with a change of slope at the transition scale between the one and the two halo terms. We detect a scale-dependence of bias with a turn-down at scale R=2.3\pm1.5 h^-1 Mpc, in agreement with previous galaxy clustering studies. We find no significant amount of stochasticity, suggesting that a linear bias model is sufficient to describe our data. We use N-body simulations to quantify both the amount of cosmic variance and systematic errors in the measurement., Comment: 20 pages, 18 figures, accepted for publication in ApJ

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., Comment: Accepted for publication in MNRAS

We compare the efficiency of weak lensing-selected galaxy clusters counts and of the weak lensing bispectrum at capturing non-Gaussian features in the dark matter distribution. We use the halo model to compute the weak lensing power spectrum, the bispectrum and the expected number of detected clusters, and derive constraints on cosmological parameters for a large, low systematic weak lensing survey, by focusing on the $\Omega_m$-$\sigma_8$ plane and on the dark energy equation of state. We separate the power spectrum into the resolved and the unresolved parts of the data, the resolved part being defined as detected clusters, and the unresolved part as the rest of the field. We consider four kinds of clusters counts, taking into account different amount of information : signal-to-noise ratio peak counts; counts as a function of clusters' mass; counts as a function of clusters' redshift; and counts as a function of clusters' mass and redshift. We show that when combined with the power spectrum, those four kinds of counts provide similar constraints, thus allowing one to perform the most direct counts, signal-to-noise peaks counts, and get percent level constraints on cosmological parameters. We show that the weak lensing bispectrum gives constraints comparable to those given by the power spectrum and captures non-Gaussian features as well as clusters counts, its combination with the power spectrum giving errors on cosmological parameters that are similar to, if not marginally smaller than, those obtained when combining the power spectrum with cluster counts. We finally note that in order to reach its potential, the weak lensing bispectrum must be computed using all triangle configurations, as equilateral triangles alone do not provide useful information., Comment: Matches ApJ-accepted version

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses on a problem that is of crucial importance for future observations in cosmology. The shapes of distant galaxies can be used to determine the properties of dark energy and the nature of gravity, because light from those galaxies is bent by gravity from the intervening dark matter. The observed galaxy images appear distorted, although only slightly, and their shapes must be precisely disentangled from the effects of pixelisation, convolution and noise. The worldwide gravitational lensing community has made significant progress in techniques to measure these distortions via the Shear TEsting Program (STEP). Via STEP, we have run challenges within our own community, and come to recognise that this particular image analysis problem is ideally matched to experts in statistical inference, inverse problems and computational learning. Thus, in order to continue the progress seen in recent years, we are seeking an infusion of new ideas from these communities. This document details the GREAT08 Challenge for potential participants. Please visit http://www.great08challenge.info for the latest information., Comment: Published in at http://dx.doi.org/10.1214/08-AOAS222 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

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., Comment: MNRAS in press Matches accepted version. References updated

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

We derive expressions, in terms of "polar shapelets", for the image distortion operations associated with weak gravitational lensing. Shear causes galaxy shapes to become elongated, and is sensitive to the second derivative of the projected gravitational potential along their line of sight; flexion bends galaxy shapes into arcs, and is sensitive to the third derivative. Polar shapelets provide a natural representation, in which both shear and flexion transformations are compact. Through this tool, we understand progress in several weak lensing methods. We then exploit various symmetries of shapelets to construct a range of shear estimators with useful properties. Through an analogous investigation, we also explore several flexion estimators. In particular, some of the estimators can be measured simultaneously and independently for every galaxy, and will provide unique checks for systematics in future weak lensing analyses. Using simulated images from the Shear TEsting Programme (STEP), we show that we can recover input shears with no significant bias. A complete software package to parametrize astronomical images in terms of polar shapelets, and to perform a full weak lensing analysis, is available at http://www.astro.caltech.edu/~rjm/shapelets/ ., Comment: 17 pages, 3 figures. Matches published version

The Shear TEsting Programme (STEP) is a collaborative project to improve the accuracy and reliability of weak lensing measurement, in preparation for the next generation of wide-field surveys. We review sixteen current and emerging shear measurement methods in a common language, and assess their performance by running them (blindly) on simulated images that contain a known shear signal. We determine the common features of algorithms that most successfully recover the input parameters. We achieve previously unattained discriminatory precision in our analysis, via a combination of more extensive simulations, and pairs of galaxy images that have been rotated with respect to each other, thus removing noise from their intrinsic ellipticities. The robustness of our simulation approach is also confirmed by testing the relative calibration of methods on real data. Weak lensing measurement has improved since the first STEP paper. Several methods now consistently achieve better than 2% precision, and are still being developed. However, the simulations can now distinguish all methods from perfect performance. Our main concern continues to be the potential for a multiplicative shear calibration bias: not least because this can not be internally calibrated with real data. We determine which galaxy populations are responsible and, by adjusting the simulated observing conditions, we also investigate the effects of instrumental and atmospheric parameters. We have isolated several previously unrecognised aspects of galaxy shape measurement, in which focussed development could provide further progress towards the sub-percent level of precision desired for future surveys. [ABRIDGED], Comment: 27 pages, 10 figures, MNRAS submitted

The Shear TEsting Programme, STEP, is a collaborative project to improve the accuracy and reliability of all weak lensing measurements in preparation for the next generation of wide-field surveys. In this first STEP paper we present the results of a blind analysis of simulated ground-based observations of relatively simple galaxy morphologies. The most successful methods are shown to achieve percent level accuracy. From the cosmic shear pipelines that have been used to constrain cosmology, we find weak lensing shear measured to an accuracy that is within the statistical errors of current weak lensing analyses, with shear measurements accurate to better than 7%. The dominant source of measurement error is shown to arise from calibration uncertainties where the measured shear is over or under-estimated by a constant multiplicative factor. This is of concern as calibration errors cannot be detected through standard diagnostic tests. The measured calibration errors appear to result from stellar contamination, false object detection, the shear measurement method itself, selection bias and/or the use of biased weights. Additive systematics (false detections of shear) resulting from residual point-spread function anisotropy are, in most cases, reduced to below an equivalent shear of 0.001, an order of magnitude below cosmic shear distortions on the scales probed by current surveys. Our results provide a snapshot view of the accuracy of current ground-based weak lensing methods and a benchmark upon which we can improve. To this end we provide descriptions of each method tested and include details of the eight different implementations of the commonly used Kaiser, Squires and Broadhurst (1995) method (KSB+) to aid the improvement of future KSB+ analyses., Comment: 18 pages, 5 figures. Version accepted by MNRAS includes 2 extra explanatory figures and updated results for the Kuijken analysis (see astroph/0601011)

We use the non-linear spherical model in cold dark matter (CDM) cosmologies with dark energy to investigate the effects of dark energy on the growth of structure and the formation of virialised structures. We consider dark energy models with a constant equation of state parameter w. For -1