MICROSCOPE mission: first constraints on the violation of the weak equivalence principle by a light scalar dilaton

Abstract

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ötvös 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.

Scalar-tensor theories are a wide class of gravity theories that contain general relativity [1]. In the Newtonian limit, they imply the existence of a fifth force, that can be well-described by a Yukawa deviation to Newtonian gravity. Its range depends mostly on the mass of the scalar field and can vary from sub-millimetric to cosmological scales [2, 3]. It has so far been constrained on all scales from a few microns to the largest scales of the Universe (see e.g. Refs.[1, 4, 5]).

This new force may or may not be composition-dependent. A non-universal coupling implies both a violation of the weak equivalence principle (WEP) and a variation of the fundamental constants [6, 7]. The former effect has already been exploited by the Eöt-Wash group to bring the current best constraints on Yukawa-type interactions and on light dilaton interactions [8, 9, 10], while the latter allows one to set constraints on cosmological to local scales [11].

The MICROSCOPE satellite aims to constrain the WEP in space [12, 13] by measuring the Eötvös parameter, defined as the normalized difference of acceleration between two bodies $i$ and $j$ in the same gravity field, $\eta=\left(\Delta a/a\right)_{ij}=2|\mathchoice{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\displaystyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\displaystyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\displaystyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\textstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\textstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\textstyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptstyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptscriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptscriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptscriptstyle a\hfil$\crcr}}}_{i}-\mathchoice{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\displaystyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\displaystyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\displaystyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\textstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\textstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\textstyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptstyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptscriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptscriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptscriptstyle a\hfil$\crcr}}}_{j}|/|\mathchoice{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\displaystyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\displaystyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\displaystyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\textstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\textstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\textstyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptstyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptscriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptscriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptscriptstyle a\hfil$\crcr}}}_{i}+\mathchoice{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\displaystyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\displaystyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\displaystyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\textstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\textstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\textstyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptstyle a\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptscriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptscriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptscriptstyle a\hfil$\crcr}}}_{j}|$. First results [14] give

$$\eta=\left(-1\pm 27\right)\times 10^{-15}$$ (1)

at a 2-$\sigma$ confidence level. MICROSCOPE tests the WEP by finely monitoring the difference of acceleration of freely-falling test masses of different composition (Platinum and Titanium) as they orbit the Earth, measured along the principal axis of the (cylindrical) test masses. The measurement equation is given e.g. in [14] as $a_{\rm Pt}-a_{\rm Ti}=g_{x}\eta+f(\mathchoice{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\displaystyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\displaystyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\displaystyle p\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\textstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\textstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\textstyle p\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptstyle p\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptscriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptscriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptscriptstyle p\hfil$\crcr}}},n)$, where $g_{x}$ is the projection of the Earth gravity field onto the axis of the test and $f(\mathchoice{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\displaystyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\displaystyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\displaystyle p\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\textstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\textstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\textstyle p\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptstyle p\hfil$\crcr}}}{\vbox{\halign{#\cr\kern-0.7pt\cr$\mkern 2.0mu\scriptscriptstyle\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraitd}$}}{{}\hbox{$\textstyle{\montraitd}$}}{{}\hbox{$\scriptstyle{\montraitd}$}}{{}\hbox{$\scriptscriptstyle{\montraitd}$}}}\mkern-1.5mu\leaders{\hbox{$\scriptscriptstyle\mkern 0.0mu\mathrel{\mathchoice{{}\hbox{$\displaystyle{\montraita}$}}{{}\hbox{$\textstyle{\montraita}$}}{{}\hbox{$\scriptstyle{\montraita}$}}{{}\hbox{$\scriptscriptstyle{\montraita}$}}}\mkern 0.0mu$}}{\hfill}\mkern-1.5mu\fldr$\crcr\kern-0.3pt\cr$\hfil\scriptscriptstyle p\hfil$\crcr}}},n)$ is a function of the instrumental and environmental parameters and measurement noise.

The constraint (1) was obtained after analyzing only one measurement session; therefore, the error bars should be considered as the largest that can be expected from the whole MICROSCOPE mission. The statistical error is expected to decrease with increasing data and with the refinement of the data analysis by the end of the mission in 2018. In the meantime, this new constraint of the WEP can already be used to set new bounds on fifth force characteristics. This letter focuses on the implications of the first results of MICROSCOPE for an interaction between matter and a light dilaton.

Scalar fifth force. The existence of a light scalar field $\phi$ modifies the Newtonian interaction between two bodies $i$ and $j$ of masses $m_{i}$ and $m_{j}$ by a Yukawa coupling [4, 15]

$$V_{ij}(r)=-\frac{Gm_{i}m_{j}}{r}\left(1+\alpha_{ij}\hbox{e}^{-r/\lambda}\right).$$ (2)

The scalar coupling to matter $\alpha_{ij}$ can be decomposed as the product $\alpha_{i}\alpha_{j}$ of the scalar couplings to matter for each test-body measured by the dimensionless factors (e.g. [17])

$$\alpha_{i}\equiv\frac{\partial\ln m_{i}/M_{P}}{\partial\phi/M_{P}}$$ (3)

with $M_{P}^{-1}=\sqrt{4\pi G}$ the Planck mass. The range $\lambda$ of the Yukawa interaction is related to the mass of the field by $\lambda=\hbar/m_{\phi}c$. The amplitude of the WEP violation is related to the presence of a scalar field that does not couple universally to all forms of energy, contrary to general relativity. The magnitude of the scalar force varies from element to element and is characterized by $\alpha_{i}(\phi)$ which requires the determination of $m_{i}(\phi)$ and thus the specification of the couplings of the scalar field to the standard model fields. Any dynamics or gradient of this scalar fields thus induce a spatial dependence of the fundamental constants [6, 7]. For two test masses in the external field of a body $E$, the Eötvös parameter reduces to

$$\eta=\frac{(\alpha_{i}-\alpha_{j})\alpha_{E}}{1+\frac{1}{2}(\alpha_{i}+\alpha_{j})\alpha_{E}}\simeq(\alpha_{i}-\alpha_{j})\alpha_{E}.$$ (4)

In order to set constraints, we need to specify the couplings of the field to matter as well as its masses.

Figure 1: Constraints on the Yukawa potential parameters ($\alpha,\lambda$) with $q=B$. The excluded region is shown in yellow and compared to earlier constraints from Ref. [18] (dotted), Ref. [8] (dashed) and Refs. [19, 20] (dot-dashed). MICROSCOPE (solid line) improves on the Eöt-Wash contraints by one order of magnitude for $\lambda>{\rm a\,few\,}10^{5}$ m.

Figure 2: Same as Fig. 1, but with $q=B-L$, compared to the earlier constraints from Ref. [18] (dotted), Ref. [9] (dashed) and Refs. [19, 20] (dot-dashed).

Baryonic/Leptonic charges. The simplest analysis consists in assuming that the composition-dependent coupling $\alpha_{ij}$ depends on a scalar dimensionless “Yukawa charge” $q$, characteristic of each material as [8, 9]

$$\alpha_{ij}=\alpha\left(\frac{q}{\mu}\right)_{i}\left(\frac{q}{\mu}\right)_{j},$$ (5)

where $\alpha$ is a universal dimensionless coupling constant which quantifies the strength of the interaction with respect to gravity and $\mu$ is the atomic mass in atomic units (e.g. $\mu=12$ for carbon-12, or $\mu=47.948$ for titanium). Different definitions of the charge $q$ are possible depending on the detailed microscopic coupling of the scalar field to the standard model fields. At the atomic levels, taking into account the electromagnetic and nuclear binding energies, the charge are usually reduced to the materials’s baryon and/or lepton numbers ($B$ and $L$) (see e.g. Refs. [21, 22]). Hence, for a macroscopic body, we must consider its isotopic composition. Hereafter, we shall set constraints on such interactions with either $q=B$ or $q=B-L$.

Following Ref. [14] and their approximations, it is straightforward to show (using Eqs. (2) and (4)) that for MICROSCOPE, the Eötvös parameter due to a Yukawa potential is

$$\eta=\alpha\left[\left(\frac{q}{\mu}\right)_{\rm Pt}\!\!\!\!\!-\left(\frac{q}{\mu}\right)_{\rm Ti}\right]\!\!\left(\frac{q}{\mu}\right)_{E}\left(1+\frac{r}{\lambda}\right)\hbox{e}^{-\frac{r}{\lambda}}$$ (6)

where $r$ is the mean distance from the satellite to the center of the Earth . The Earth charge takes into account the Earth differentiation between core and mantle

$$\left(\frac{q}{\mu}\right)_{E}=\left(\frac{q}{\mu}\right)_{\rm core}\Phi\left(\frac{R_{c}}{\lambda}\right)+\left(\frac{q}{\mu}\right)_{\rm mantle}\left[\Phi\left(\frac{R_{E}}{\lambda}\right)-\Phi\left(\frac{R_{c}}{\lambda}\right)\right],$$ (7)

where $R_{E}$ is the Earth mean radius and $R_{c}$ the Earth core radius. The function $\Phi(x)\equiv 3(x\cosh x-\sinh x)/x^{3}$ [4] takes into account the fact that all Earth elements do not contribute similarly to the Yukawa interaction at the satellite’s altitude ($\Phi=1$ for the test masses since their sizes are much smaller than the ranges $\lambda$ that can be probed in orbit). We assume that the core of the Earth is composed of iron and that the mantle is composed of silica (SiO₂) [25]. The baryonic and lepton charges for the MICROSCOPE experiment are summarized in Table 1.

At the 2-$\sigma$ level, MICROSCOPE’s constraints on the Eötvös parameter are given by Eq. (1), and can readily be transformed into constraints on Yukawa’s ($\alpha$, $\lambda$). Figs 1 and 2 depict the corresponding exclusion regions respectively for $q=B$ and $q=B-L$. In both analyses, we compare our new constraint to the bounds from Eöt-Wash’s torsion pendulum experiments [18, 8, 9] and the constraints from the Lunar-Laser Ranging (LLR) experiment [19, 20]. Note that while we plot only the latest, most competitive constraints, several other experimental constraints are available (e.g. [4, 26, 27, 28, 29, 30, 31, 32, 33, 34]). Moreover, the LLR constraint could be slightly strengthened in the near future [35]. This shows that MICROSCOPE’s first results allow us to gain one order of magnitude compared to previous analyses for $\lambda>{\rm a\,few\,}10^{5}$m. As MICROSCOPE orbits Earth at about 7000 km from its center, one would naively expect that it can only probe interactions with $\lambda>{\rm a\,few\,}10^{6}{\rm m}$; smaller ranges could not be probed as they imply too much of a damping at MICROSCOPE’s altitude. However, if a fifth force with $\lambda\approx{\rm a\,few\,}10^{5}{\rm m}$ was strong enough to affect MICROSCOPE, the contribution from the nearest point of the Earth (as seen from MICROSCOPE) would be higher than that of the farthest point of the Earth, implying an asymmetric behavior that can be probed by MICROSCOPE (as captured by the function $\Phi(x)$ above). Hence, MICROSCOPE is sensitive to scalar interactions with ranges as low as a few hundreds of kilometers.

Dilaton models. We now consider the characteristics of a generic dilaton with couplings described in Refs [36, 25, 37]. The mass of an atom (atomic number $Z$ and mass number $A$) can be decomposed as $m(A,Z)=Zm_{\rm p}+(A-Z)m_{\rm n}+Zm_{\rm e}+E_{1}+E_{3}$ where $m_{\rm n,p}$ is the mass of the neutron or proton and $E_{1}$ and $E_{3}$ are the electromagnetic and strong interaction binding energies. Following Ref. [25], we consider that the coupling coefficients of the dilaton to the electromagnetic and gluonic fields are $d_{e}$ and $d_{g}$ while $d_{m_{e}}$, $d_{m_{u}}$ and $d_{m_{d}}$ are its coupling to the electron, $u$ and $d$ quarks mass terms. The latter two can be replaced by the couplings $d_{\delta m}$ and $d_{\tilde{m}}$ to the symmetric and antisymmetric linear combination of $u$ and $d$. Assuming a linear coupling, one deduces that the variation of the fine structure constants and masses of the quarks are given by $\Delta\alpha_{\rm EM}/\alpha_{\rm EM}=d_{e}\phi/M_{p}$ and $\Delta m_{u,d}/m_{u,d}=d_{u,d}\phi/M_{p}$.

Table 1: Baryonic, leptonic and dilaton charges for MICROSCOPE’s test masses.

Material	$B/\mu$	$(B-L)/\mu$	$Q^{\prime}_{\tilde{m}}$	$Q^{\prime}_{e}$
Pt/Rh	1.00026	0.59668	0.0859	0.0038
Ti/Al/V	1.00105	0.54044	0.0826	0.0019

Figure 3: Constraints on the couplings of a massless dilaton $(D_{\tilde{m}},D_{e}$). The region allowed by the MICROSCOPE measurement (black band) is compared to earlier constraints by torsion pendulum experiments from Ref. [38] (green) and Ref. [9] (yellow, cyan). The difference of slopes arises from the difference of material used in these 3 experiments. MICROSCOPE allows us to shrink the allowed region by one order of magnitude.

First, we consider a massless dilaton ($m_{\phi}=0$), whose range $\lambda_{\phi}$ is infinite, as was done by the Eöt-Wash group [9]. The dilaton coupling to matter, and hence the fifth force, is parametrized by the 5 numbers ($d_{g}$, $d_{e}$, $d_{\tilde{m}}$, $d_{\delta m}$, $d_{m_{e}}$) so that the coupling to matter (3) takes the form

$$\alpha_{i}\approx d_{g}^{*}+\left[\left(d_{\tilde{m}}-d_{g}\right)Q^{\prime}_{\tilde{m}}+d_{e}Q^{\prime}_{e}\right]_{i},$$ (8)

where $d_{g}^{*}=d_{g}+0.093(d_{\tilde{m}}-d_{g})+0.00027d_{e}$. The dilaton charges depend on the chemical composition of the test masses and on the local value of the dilaton. Following Ref. [25], they are well-approximated by

$$Q^{\prime}_{\tilde{m}}=0.093-\frac{0.036}{A^{1/3}}-1.4\times 10^{-4}\frac{Z(Z-1)}{A^{4/3}}$$ (9)

and

$$Q^{\prime}_{e}=-1.4\times 10^{-4}+7.7\times 10^{-4}\frac{Z(Z-1)}{A^{4/3}}.$$ (10)

In the limit where $\lambda$ is much larger than any other spatial scales, the Eötvös parameter reduces to Eq. (4) so that (at first order in dilaton charges $Q^{\prime}_{j}$ –given that $|Q^{\prime}_{j}|\ll 1$)

$$\eta_{\rm massless}=D_{\tilde{m}}\left([Q^{\prime}_{\tilde{m}}]_{\rm Pt}-[Q^{\prime}_{\tilde{m}}]_{\rm Ti}\right)+D_{e}\left([Q^{\prime}_{e}]_{\rm Pt}-[{Q^{\prime}_{e}}]_{\rm Ti}\right),$$ (11)

where the coefficients $D_{\tilde{m}}=d_{g}^{*}(d_{\tilde{m}}-d_{g})$ and $D_{e}=d_{g}^{*}d_{e}$ are to be estimated. The values for $Q^{\prime}_{\tilde{m}}$ and $Q^{\prime}_{e}$ in the MICROSCOPE case are given in Table 1.

Fig. 3 summarizes our new constraints and compare them to the earlier ones from the Eöt-Wash [9] and the Moscow groups [38]. The different slopes of the allowed regions are due to the different pairs of materials used by each experiment.

Massive dilaton. The mass of the dilaton modifies the range of its interaction so that Eq. (11) is modified as

$$\eta=\eta_{\rm massless}\times\Phi\left(\frac{R_{E}}{\lambda_{\phi}}\right)\left(1+\frac{r}{\lambda_{\phi}}\right)\hbox{e}^{-r/\lambda_{\phi}}.$$ (12)

Note that this equation is simpler than Eq. (7) because Eq. (11) does not depend on the Earth dilaton charge, and it is therefore independent of the exact Earth model used.

From Figs. 1 and 2, we expect that MICROSCOPE shall mainly be sensitive to masses in the range $10^{-14}-10^{-12}$ eV. Lower masses will result in constraints similar to those for a massless dilaton (see Fig. 3) while larger masses cannot be constrained, as they correspond to ranges that MICROSCOPE cannot probe. This is indeed what we conclude from our analysis summarized in Fig. 4. Constraints in the $(D_{\tilde{m}},D_{e}$) plane are rather loose for high-enough masses, $m_{\phi}>10^{-12}$eV, and converge to those of a long-range dilaton for $m_{\phi}<10^{-14}$eV.

Figure 4: Constraints on the couplings of a massive dilaton for various values of its mass. Each color shows the allowed $(D_{\tilde{m}},D_{e})$ for a given mass of the scalar field. The inset is a zoom on smaller $(D_{\tilde{m}},D_{e})$. Constraints saturate for light fields $m_{\phi}<10^{-14}$ eV. MICROSCOPE is not sensitive to masses larger than a few $10^{-12}$ eV.

Finally, we assume that the dilaton field couples only to the electromagnetic field, i.e. the only non-vanishing coupling is $d_{e}$. The coupling to proton and neutron is then induced from their binding energy [39]. Several groups set constraints on such a dilaton from the fine structure constant oscillations in atomic frequency comparisons [40, 41, 42]. These results are based on the time evolution of the scalar field that oscillates within its self-potential. It has been argued that these oscillations may lead to oscillations of the Newtonian potentials if the scalar field behaves like cold dark matter [43] (thereby affecting MICROSCOPE in an unexpected way), or even break the Yukawa approximation [44]. Here, we do not tie our scalar field to describe dark matter and we restrict our analysis to linear couplings, thence avoiding those possible pitfalls . The MICROSCOPE constraints are obtained by considering the $D_{\tilde{m}}=0$-subspace of the parameter space $(D_{\tilde{m}},D_{e},m_{\phi})$ of Fig. 4, and recognizing that $D_{e}=d_{g}^{*}d_{e}=0.00027d_{e}^{2}$. Fig. 5 shows our constraints, compared with those from the Eöt-Wash test of the WEP and with atomic spectroscopy [40, 41]. MICROSCOPE allows us to exclude a new region above $|d_{e}|=10^{-4}$, for a field of mass $10^{-18}<m_{\phi}/{\rm eV}<10^{-11}$. Atomic spectroscopy stays more competitive for lighter fields.

Figure 5: Constraints on $d_{e}$, for a dilaton coupled only to the electromagnetic sector, compared with constraints from atomic spectroscopy (dot-dashed [40, 41]) and Eöt-Wash WEP test (dashed [8]).

Conclusion. This letter gave the first constraints on a composition-dependent scalar fifth force from MICROSCOPE’s first measurement of the WEP [14]. We first considered the case of a massive scalar field coupled to either $B$ or $B-L$ to conclude that MICROSCOPE is particularly competitive for a Yukawa potential of range larger than $10^{5}$m (corresponding to a field of mass smaller than $10^{-12}$eV). In that case, we improved existing constraints on the strength of the field by one order of magnitude. Below that range, torsion pendulum experiments remain unbeaten. Then, we considered a model describing the coupling of a generic dilaton to the standard matter field with 5 parameters, both for a massless and massive field. For $m_{\phi}<10^{-14}$ eV, our constraints are similar to those for a massless field and better by one order of magnitude than the previously published ones.

From a theoretical perspective, a scalar long-range interaction is severely constrained by its effects on planetary motion. Since general relativity passes all tests on Solar-System scales many mechanisms have been designed to hide this scalar field in dense regions (e.g. chameleons [46, 47], symmetron [48], K-mouflage [49, 50] or Vainshtein [51]). The generic dilaton model considered in this letter corresponds to another type of screening (the least coupling principle [36]) and can incorporate the behavior of many theories, such as string theory. The local prediction of the violation of the WEP can be compared to the variation of the fundamental constants on local and astrophysical scales (e.g. [52, 53, 54, 55]). Better constraints can be obtained from modeling the profile (and time variation) of the scalar field along MICROSCOPE’s orbit, as well as its propagation inside the satellite up to the test masses; this is non-trivial, requires some care, and will be done in a further work. Constraints on the violation of the WEP will also have strong consequences for bigravity models [56].

From an experimental perspective, these new constraints were obtained from only two MICROSCOPE’s measurement sessions of the Eötvös parameter [14]. As the mission is scheduled to continue until 2018, new data are currently coming in, thereby offering the possibility of decreasing the statistical errors. We are also refining our data analysis procedures to optimize the measurement of the WEP. We therefore expect to improve on MICROSCOPE’s constraint on the Eötvös parameter by the end of the mission: ten times as many data will be available than were used in Ref.[14]; furthermore, although we expect the data to become systematic-dominated, the control on systematics will be improved compared to Ref. [14], since calibration sessions have been performed, whose results will be used in the next data analysis. Therefore, we could improve the constraints reported in that letter by up to another order of magnitude (unless a WEP violation becomes apparent). But this forecast is valid only for $\lambda>{\rm a\,few\,}10^{5}$ m ($m_{\phi}<10^{-12}$ eV). Probing lower-range (more massive) scalar fields can be done only using small scale experiments. Torsion pendulum and atomic interferometry experiments represent our best hopes to look for such extra-fields. New, improved torsion pendulum will then be required to probe laboratory and smaller scale gravity, either through the measurement of the WEP or of the gravitational inverse square law. A torsion pendulum experiment in space seems the way forward to beat the current on-ground limits [57].

Acknowledgements. We thank the members of the MICROSCOPE Science Performance Group for useful discussions. We also thank Thibault Damour for useful comments on a first version of the manuscript. This work makes use of technical data from the CNES-ESA-ONERA-CNRS-OCA Microscope mission, and has received financial support from ONERA and CNES. We acknowledge the financial support of the UnivEarthS Labex program at Sorbonne Paris Cité (ANR-10-LABX-0023 and ANR-11-IDEX-0005-02). The work of JPU is made in the ILP LABEX (under reference ANR-10-LABX-63) was supported by French state funds managed by the ANR within the Investissements d’Avenir programme under reference ANR-11-IDEX-0004-02. This work is supported in part by the EU Horizon 2020 research and innovation programme under the Marie-Sklodowska grant No. 690575. This article is based upon work related to the COST Action CA15117 (CANTATA) supported by COST (European Cooperation in Science and Technology).

We are interested in scalar-tensor theories of the type

$$S=\int d^{4}x\sqrt{-g}\left(\frac{R}{16\pi G_{N}}-\frac{1}{2}(\partial\phi)^{2}-\frac{m^{2}}{2}\phi^{2}\right)+\sum_{i}S_{m}(A^{2}_{i}(\phi)g_{\mu\nu})$$ (S1)

where different particle species $i$ couple to matter via the coupling function $i$. In the core of the paper we also consider the coupling to gauge fields in the form of

$$S_{F}=-\int d^{4}x\sqrt{-g}\frac{B_{i}(\phi)}{4}F^{\mu\nu}_{i}F_{\mu\nu}^{i}$$ (S2)

which plays a role to give the gluon and electromagnetic parts of the atomic masses. In this section we focus on the role of the coupling $A_{i}$ whilst the generalisation which includes the $B_{i}$’s is taken into account in the text following Damour and Donoghue [37], and Damour and Polyakov [36].

The scalar field is a canonically normalised scalar field with a mass $m$. We focus on linearly coupled scalars where

$$A_{i}(\phi)=e^{\alpha_{i}\phi/\sqrt{2}m_{\rm Pl}}$$ (S3)

where here and contrarily to the main text $m_{\rm pl}^{2}=1/8\pi G_{N}$. This choice is crucial as if the potential were non-linear such as the inverse power law $V(\phi)\sim 1/\phi^{n}$ for instance, the effects of the scalar field on matter could be screened, see the review [58] and references therein. If the coupling function were universal $A_{i}\equiv A$ and of the scalarisation type $A(\phi)=e^{\beta\phi^{2}/2m_{\rm Pl}^{2}}$ as in [44], the solutions to the Klein-Gordon equation in the presence of matter would be destabilised due to the negative coupling $\beta<0$ and lead to a completely different phenomenology. Here we restrict ourselves to the linear couplings (S3) leading to a Yukawa interaction between massive objects as we recall below.

Notice too that we do not tie our scalar field to describe dark matter. Indeed, a massive scalar oscillating around the minimum of its potential behaves like Cold Dark Matter over time scales larger than the oscillation time. This is what is assumed in [43] where the oscillations of the coupled scalar lead to oscillations of the Newtonian potentials, which could become observable in the change of periods of binary pulsars. In our case, the field can be assumed to store an energy in the form of oscillations which is negligible compared to the dark matter energy density. This can be easily realised if for instance the scalar field is present during inflation and its amplitude decays exponentially fast in $a^{-3/2}$ towards zero, in such a way that the field is static at the end of inflation.

Let us recall how Eq. (2) of the main text is derived, see [59] that we adapt to the present situation. First of all let us recall that the Einstein equation reads

$$G_{\mu\nu}=8\pi G_{N}(T^{m}_{\mu\nu}+T^{\phi}_{\mu\nu})$$ (S4)

where the energy momentum of matter is given by

$$T^{m}_{\mu\nu}=A(\phi)\rho u_{\mu}u_{\nu}$$ (S5)

where $\rho$ is the conserved matter density and $u_{\mu}$ the velocity 4-vector of matter. Working in the Newtonian gauge for the metric

$$ds^{2}=-(1+2\Phi_{N})dt^{2}+(1-2\Phi_{N})dx^{2}$$ (S6)

we can expand both the scalar field and the Newtonian potential as

$$\Phi_{N}=\bar{\Phi}+\delta\Phi,\ \phi=\bar{\phi}+\delta\phi$$ (S7)

where a bar denotes a background quantity and “$\delta$” the fields sourced by objects. In the vicinity of an object, the background field can be expanded to linear order

$$\bar{\Phi}(\vec{x})=\bar{\Phi}_{0}+\partial_{i}\bar{\Phi}(\vec{x})x^{i},\ \bar{\phi}(\vec{x})=\bar{\phi}_{0}+\partial_{i}\bar{\phi}(\vec{x})x^{i}.$$ (S8)

The field outside a given object $B$ and created by this object satisfies the static Klein-Gordon equation, where here we assume that the size of the objects is much smaller than the range $1/m$ of the scalar interaction,

$$\Delta\delta\phi-m^{2}\delta\phi=\frac{\alpha_{B}m_{B}}{\sqrt{2}m_{\rm Pl}}\delta^{(3)}(\vec{x})$$ (S9)

and the matter density has been approximated by a point-like source. The solution is the Yukawa potential

$$\delta\phi=-\frac{\alpha_{B}m_{B}}{4\pi\sqrt{2}m_{\rm Pl}r}e^{-mr}.$$ (S10)

where $r=\sqrt{\vec{x}^{2}}$. For extended objects, we need to integrate this point-like solution, which leads to the shape function of the main text’s Eq. (7). We also assume that the object creates the Newtonian potential

$$\delta\Phi(r)=-\frac{G_{N}m_{B}}{r}.$$ (S11)

We assume that matter is responsible for the Newtonian potential, and that the scalar field energy scale is negligible compared to matter inside the object and very small outside. This is the case for instance if the scalar field sits at the minimum of the potential $\bar{\phi}_{0}=0$.

The Einstein equation can be rewritten as

$$G^{(1)}_{\mu\nu}=8\pi G_{N}(T^{m}_{\mu\nu}+T^{\phi}_{\mu\nu}+t_{\mu\nu}).$$ (S12)

where $G^{(1)}_{\mu\nu}$ is linear in the Newtonian potential and the pseudo-tensor is given by

$$t_{\mu\nu}=-\frac{1}{8\pi G_{N}}G^{(2)}_{\mu\nu}$$ (S13)

where $G^{(2)}_{\mu\nu}$ contains all the higher order terms in the Newtonian potential. The mass $m_{B}$ of the object becomes

$$m_{B}=-\int_{V}d^{3}x\tilde{T}_{0}^{0}.$$ (S14)

where we draw a sphere of volume $V$ encircling the object and

$$\tilde{T}_{\mu\nu}=T^{m}_{\mu\nu}+T^{\phi}_{\mu\nu}+t_{\mu\nu}.$$ (S15)

Neglecting the scalar contribution to the energy density, the mass is given by the integral over the object

$$m_{B}=\int_{V}A(\phi)\rho d^{3}x.$$ (S16)

The momentum of the object is then

$$P^{B}_{i}=\int_{V}d^{3}x\tilde{T}_{i}^{0}.$$ (S17)

Using the non-covariant Bianchi identity which implies that $\partial^{\mu}\tilde{T}_{\mu\nu}=0$, we find that

$$\dot{P}^{B}_{i}=-\int_{\partial V}dS_{j}\tilde{T}^{j}_{i}$$ (S18)

where the surface integral is on the surface of the outer sphere. There the matter energy momentum tensor is negligible, and similarly for the contribution from the scalar field energy density. Only two terms have a relevant flux: the scalar and gravitational ones. The gravitational flux has been computed in [59] and yields

$$\int_{\partial V}dS_{j}t^{j}_{i}=m_{B}\partial_{i}\bar{\Phi}.$$ (S19)

The new contribution from the scalar field is simply dominated by the large gradient of the scalar field $\delta\phi(r)$

$$-\int_{\partial V}dS_{j}T^{\phi j}_{i}=-\frac{\alpha_{B}m_{B}}{\sqrt{2}m_{\rm Pl}}\partial_{i}\bar{\phi}\,.$$ (S20)

As a result we obtain that

$$\dot{P}^{B}_{i}=-m_{B}\partial_{i}\bar{\Phi}-\frac{\alpha_{B}m_{B}}{\sqrt{2}m_{\rm Pl}}\partial_{i}\bar{\phi}.$$ (S21)

Similarly the centre of mass coordinates satisfies the modified Newton equation

$$\ddot{X}_{B}^{i}=-\partial^{i}\bar{\Phi}-\alpha_{B}\frac{\partial^{i}\bar{\phi}}{\sqrt{2}m_{\rm Pl}}.$$ (S22)

When the external fields $\bar{\Phi}$ and $\bar{\phi}$ are due to another extended object, we have

$$\bar{\Phi}=-\frac{G_{N}m_{A}}{|\vec{x}-\vec{x}_{A}|},\ \ \bar{\phi}=-\frac{\alpha_{A}m_{A}}{4\pi\sqrt{2}m_{\rm Pl}|\vec{x}-\vec{x}_{A}|}e^{-m|\vec{x}-\vec{x}_{A}|}.$$ (S23)

The motion of the object $B$ is due to the total potential

$$\Phi_{AB}=(1+\alpha_{A}\alpha_{B}e^{-m|\vec{x}-\vec{x}_{A}|})\Phi_{A}$$ (S24)

where $\Phi_{A}$ is the Newtonian potential due to the second object $A$ and we have

$$\ddot{X}^{i}_{A}=-\partial^{i}\Phi_{AB}.$$ (S25)