No Arabic abstract
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.
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 missions expected error budget.
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.
MICROSCOPEs 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 satellites 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 satellites 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.
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.
Gravitational waves are radiative solutions of space-time dynamics predicted by Einsteins theory of General Relativity. A world-wide array of large-scale and highly sensitive interferometric detectors constantly scrutinizes the geometry of the local space-time with the hope to detect deviations that would signal an impinging gravitational wave from a remote astrophysical source. Finding the rare and weak signature of gravitational waves buried in non-stationary and non-Gaussian instrument noise is a particularly challenging problem. We will give an overview of the data-analysis techniques and associated observational results obtained so far by Virgo (in Europe) and LIGO (in the US), along with the prospects offered by the up-coming advance