Secular gravity gradients in non-dynamical Chern–Simons modified gravity for satellite gradiometry measurements

With continuous advances in related technologies, precision tests of modern gravitational theories with orbiting gradiometers becomes feasible, which may naturally be incorporated into future satellite gravity missions. In this work, we derive, at the post-Newtonian level, the new secular gravity gradient signals from the non-dynamical Chern–Simons modified gravity for satellite gradiometry measurements, which may be exploited to improve the constraints on the mass scale MCS\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{CS}$$\end{document} or the corresponding length scale θ˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\dot{\theta }}$$\end{document} of the theory with future missions. For orbiting superconducting gradiometers, a bound MCS≥10-7eV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{CS}\ge 10^{-7}\ \mathrm{eV}$$\end{document} and θ˙≤1m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\dot{\theta }} \le 1\ \mathrm{m}$$\end{document} could in principle be obtained, and for gradiometers with optical readout based on the similar technologies established in the LISA PathFinder mission, an even stronger bound MCS≥10-6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{CS}\ge 10^{-6}$$\end{document}–10-5eV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-5}\ \mathrm{eV}$$\end{document} and θ˙≤10-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\dot{\theta }} \le 10^{-1}$$\end{document}–10-2m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-2} \ \mathrm{m}$$\end{document} might be expected.


Introduction
Among modifications of Einstein's general relativity (GR), extensions to the Einstein-Hilbert action with second order curvature terms are of particular interest, which may arise from the full, but still lacking, quantum theory of gravity [1]. The Chern-Simons (CS) modified gravity [2][3][4][5][6] belongs to such extensions of GR, which has physical roots in particle physics and string theory. In particle physics, the CS modification or the Pontryagin term R R is known to be related to the chiral current anomaly that caused by spacetime curvature [7,8], which may also serve as a possible candidate that sources the baryon asymmetry of the Universe [9] through gravi-leptogenesis [10]. In string theory, the CS modification a e-mail: xupeng@amss.ac.cn b e-mail: qqllee815@chd.edu.cn emerges as an anomaly-canceling term through the Green-Schwarz mechanism [11]. More interestingly, the CS extension to GR may provide us insights into the physics of possible parity-violations in gravitation, that includes effects like amplitude birefringent gravitational waves [5,12,13], different gravito-magnetic (GM) sectors [12][13][14], and etc.. Therefore, experimental tests of the CS modified gravity and the resulted constraints are of importance.
The CS modified gravity is considered as an effective theory, that the ultra-violet modifications to gravitation and their possible observable effects are to be studied in more fundamental and sophisticated theories such as string theory or loop quantum gravity. Theoretical studies had shown that in CS modified gravity the Lorentz symmetry can be satisfied [5,15,16], and, up to now, experimental constraints on CS gravity are mainly from astrophysical observations and Solar system tests. Of course, it is natural to expect that stronger bounds on CS gravity might come from particle physics because of the connections between these theories. In this work, we will focus on the tests, with orbiting gradiometers, of the non-dynamical formulation of the CS modified gravity [5], where the coupling field or the deformation parameter θ is externally prescribed (one notices that the arbitrariness in θ could not be completely removed in the dynamical formulation due to the different choices of the potential V (θ )). The first constraint [14] on the time derivative of the coupling scalarθ and the corresponding mass scale M C S ∼ 1/θ of the non-dynamical CS gravity was obtained based on the observations from the LAGEOS I, II [17,18] and the Gravity Probe-B [19] missions, which had set M C S ≥ 10 −13 eV anḋ θ ≤ 10 6 m, and a stronger bounds M C S ≥ 4.7 × 10 −10 eV andθ ≤ 0.4 × 10 3 m (been revised in [20]) was obtained based on the data from double binary pulsars [21].
Theoretical studies of testing relativistic gravitational theories with orbiting gravity gradiometers in space were firstly carried out in 1980s [22][23][24], and such measurement schemes could be naturally incorporated into future satellite gradiometry missions or missions that carrying high sensitive gradiometer as one of the key payloads. For the baseline design of high sensitive gravity gradiometers in micro-gravity or zero-g environment in space, such as electrostatic or superconducting ones, one generally has pairs of proof masses aligned along each of the measurement axes, and a combinations of strategies of proof mass disturbances isolation, proof mass position sensing and control is employed, see [25][26][27][28] for reviews. The proof masses are generally enclosed within sensor cages or housings, vacuum maintenances and other shielding devices, and, with such setup, fluctuations forces subjected to proof masses are to be reduced or isolated as much as possible. The relative motions or accelerations between the "free-falling" proof masses (with respect to certain noise level) in space will give rise to measurements of the tidal matrix from spacetime curvature R j 0i0 along certain orbits (for Newtonian limits, R j 0i0 reduces to ∂ i ∂ j U with U the Newtonian gravitational potential). For electrostatic and superconducting gradiometers, the difference between the compensating forces that restoring the proof masses to their nominal positions can be used as the direct readouts of the tidal accelerations. The GOCE satellite, launched in March 2009, carried an electrostatic gravity gradiometer containing six proof masses in three pairs to map out the geopotential of Earth, whose sensitive had reached 10 mE/Hz 1/2 in the frequency band of 5 ∼ 100 mHz [27]. With the continuous advances, the multi-axis superconducting gravity gradiometer under the development could reach the sensitivity about 10 −2 mE/Hz 1/2 in the band around 1 mHz in space [25,29]. As an alternative optical readout method, the relative motions between proof masses as integrations of tidal accelerations can also be precisely measured by onboard laser interferometers 1 . The LISA PathFinder (LPF) mission [30,31], which can be view as a demonstration of an one dimensional optical gradiometer with the resolution of the onboard laser interferometer better than 9 pm/ √ Hz in the m H z band, had even reached the noise floor of 10 −3 -10 −4 mE/Hz 1/2 . Tests of relativistic gravitational theories including the CS modified gravity with satellite gradiometry now becomes more and more feasible, and for CS gravity a preliminary measurement scheme had been studied in [32,33]. It was firstly noticed by Mashhoon and Theiss [22,[34][35][36] that along orbit motions relativistic secular tidal effects that growing with time may exist (known as the Mashhoon-Theiss anomaly), which would greatly improve the measurement accuracy of relativistic tidal components. Recently, the physical mechanism behind such secular tidal effects had been studied and explained in [37,38]. From the post-Newtonian (PN) point of view, the difference between the relativistic precession of the local free-falling frame (or the parallel transported measurement axes) and the orbit plane with respect to the sidereal frame will produce modulations of Newtonian tidal forces along certain axes and then gives rise to periodic secular tidal signals [37,39]. Back to the tests of the CS modified gravity, in the far field expansion of the non-dynamical theory, the CS gravity will add modifications to the GM sector of GR [12,13], and therefore will give rise to new secular tidal effects that could be read out precisely along certain measurement axes of an orbiting gradiometer. In this work, we derive, at the PN level, the new secular tidal tensor from the non-dynamical CS modified gravity under the local Earth pointing frame along a relativistic polar and nearly circular orbit. For (possible) future experiments, we give the estimations of the bound on the characteristic mass scale M C S that could be drawn from such a measurement scheme.

Models and settings
The action of the non-dynamical CS modified gravity is given by for clarity the Natural units c = h = 1 are adopted hereafter and in the end the SI units will be recovered. In the nondynamical theory with the canonical coupling [5], the scalar θ is an externally prescribed and spatially isotropic function that proportional to the coordinate time. The field equation reads where The introduction of the new scalar degree of freedom gives rise to a new constraint If the constraint is satisfied, the Bianchi identities and the equations of motion for matter fields ∇ ν T νμ = 0 are recovered, which ranks the non-dynamical CS modified gravity a metric theory [40,41].
In this work, we model Earth as an ideal uniform and rotating spherical body with total mass M and angular momentum J. The geocentric inertial coordinates system {t, x i } is defined as follows, that one of its bases ∂ ∂ x 3 is parallel to the direction of J and the coordinate time t is measured in asymptotically flat regions. For an orbiting proof mass or satellite, we have the PN order relations where v is the 3-velocity, r = 3 i=1 (x i ) 2 and for low (with altitude below 2000 km) and medium (altitude between 2000 km and a geostationary orbit) Earth orbits = M r is about 10 −5 -10 −6 . Up to the required order, the metric field outside the ideal Earth model can be expanded as where the mass scale M C S of the CS modified gravity reads As mentioned, the non-dynamical CS modified gravity differs from GR only in the GM sector.

Reference orbit and local tetrad
Being a metric theory, motions of free-falling masses or satellites in CS modified gravity satisfy the geodesic equation.
According to the general choices of orbits for satellite gradiometry missions (like GOCE [27]), we choose the reference orbit followed by the mass center of the gradiometer to be a polar and circular one with the relativistic precession caused by the GM effect x 2 = a sin Ψ.
Here a denotes the orbit radius, Ψ = ωτ the true anomaly, ω is the mean angular frequency with respect to the proper time τ along the orbit and Ω is the longitude of ascending node with initial value Ω(0) = 0. The precession rate of the nodė Ω =Ω G R +Ω C S , where the Lense-Thirring precession ratė Ω G R = 2G J a 3 [42] and the correction from the non-dynamical CS modified gravity had been worked out in [14] aṡ where R is the averaged radius of Earth, j l (x) and y l (x) are spherical Bessel functions of the first and second kind, and χ was a new PN parameter introduced in [12,13] which can be related to the mass scale as χ = 4 a M C S . For polar circular orbit, the GM force in CS modified gravity generated by spherical sources [14] will also change slightly the orbital eccentricity to e ∼ χ O( 2 ), which is too small to be relevant to secular tidal effects at the PN level. Therefore, in this work the small eccentricity is ignored, and its effects together with other orbital perturbations in satellite gradiometry, such as those from geopotential harmonics, can be found in [27,43].
For satellite gradiometry missions, spacecraft attitudes are generally chosen to follow the Earth pointing orientation (like GOCE [27]). Then, we define the local free-falling Earth pointing frame by the tetrad {E  [44], we solve for the spatial bases {E μ (i) } in the following three steps. First, in the geocentric coordinates system, we solve for the precession of the local inertial frame (Fermishifted frame) along the orbit given in Eqs. (7), (8). Second, with respect to the local inertial frame, we rotate {E μ (i) } with an initial angular velocity to make it an Earth pointing triad. At last, since the local frame is moving along the orbit, we need to perform the boost Lorentz transformations of the bases {E μ (i) } with respect to the 4-velocity τ μ . The general time scales or periods of frame-dragging precessions in Earth orbit are about c 2 a 3 G J ∼ 10 7 years, which is extremely long compared with general mission lifetimes. Then, following the above three steps and within the short time limit τ a a 2 G J , the tetrad can be worked out up to the PN level as where the correction to the precessions of the bases from the non-dynamical CS modified gravity reads [14] Δ C S = 15

Secular gradient observables
We introduce the position difference vector Z μ between the pair of two adjacent free-falling proof masses of certain measurement axis in the orbiting gradiometer. Generally |Z | ∼ 10 −1 m, which is much shorter compared with the orbital radius a ∼ 10 7 m, therefore the relative motion between the test masses can be obtained by integrating the geodesic deviation equation along the reference orbit In the local frame {E μ (a) }, the above geodesic deviation equation can be expanded as where are the Ricci rotation coefficients [45]. The first line of the right hand side of the above equation is the relativistic analogue of the Coriolis force, the second line contains the inertial tidal forces and the last line is the tidal force from the spacetime curvature, where the tidal matrix from curvature is defined by For electrostatic and superconducting gradiometers, the motions of test masses are suppressed by compensating forces. Then the total tidal tensor T (a)(b) affecting the gradiometer will be After straightforward but tedious algebraic manipulations and leaving out all the terms beyond 1 a 2 O( 4 ) and 1 a 2 Ψ O( 4 ), we work out, to the PN level, the tidal tensors in the local free-falling Earth pointing frame along the reference orbit. For the tidal tensor K (a)(b) , as expected we have K (a)(0) = 0, and the Newtonian part which agrees exactly with the classical Newtonian tidal tensor ∂ i ∂ j G M r evaluated in such Earth-pointing frame. The PN part may be divided into the tidal tensor K G R (i)( j) from GR and the new tensor K C S (i)( j) from the CS modification, which, within the short time limit τ a a 2 G J , can be worked out as here K G R (i)( j) agrees exactly with the former result derived in [37]. Due to the relativistic precessions of the free-falling local frame and the orbit plane, the modulations of Newtonian tidal tensor given in Eq. (19) produces secular terms in the K (2)(3) and K (3) (2) components, which are the expected secular gradient observables appeared along polar and nearly circular orbits. Finally, with the tidal tensor from inertial forces in Eq. (18) been worked out, the total tidal tensor T (i)( j) turns out to be

Conclusions
In conclusion, how the new secular gradient components given in Eq. (21) can be measured by an orbiting 3-axis gradiometer is discussed, and the expected constraints on the mass scale of the non-dynamical CS gravity set by the proposed measurement scheme is given. Following [23,32,46], we orient two of the three gradiometer axes 45 degrees above and below the orbital plane and difference their outputs to reject the Newtonian and PN gravitoelectric terms and therefore measure only the GM gradient terms. In the local frame {E μ (a) }, the three axes of the gradiometer are oriented aŝ Fig. (1) for the illustration. The difference between the readouts in thep andq axes turns out to be The boxed term is the secular gradient signal s C S to be measured, which is from the non-dynamical CS modifications and grows linearly with time within the mission lifetime. Such combinations of readouts can be obtained without really re-orientating the gradiometer axes according to Fig. 1 in mission operations, but can be derived in the post data procession by combining the cross-track readouts of multi-axis gradiometers (like in GOCE [27]). As M C S → ∞ or χ → 0 we have Π C S → 0 and Δ C S → 0, the above combined readout will then reduce to the GM gradient terms from GR as expected [37]. The magnitudes of these GR gradient terms, especially for the secular one to the right of the boxed signal, are suppose to be very large compared with that of the CS secular signal s C S . In the low energy regime or the weak field limit, GR had passed many stringent tests [41]. From the PN point of view, the weak field limits of general metric theories of gravity [40] can be described by ten standard PN parameters (CS gravity and the related parameter not included), which, confronted with modern experiments, have been constrained to agree precisely with the values predicted by GR, see [41] for details. Therefore the relative large tidal gradient terms from GR in Eq. (23) can be treated as from a well-test theory, and are to be modeled by GR in post data analysis. The CS gradient signal s C S and the new parameter χ or M C S is then to be fitted out in the residuals or measured as the deviations from the prediction of GR. Similar strategy has been used in fitting the signals from the non-dynamical CS gravity in the data of the LAGEOS I, II missions and the Gravity Probe B mission [14]. Errors in such combination 1 2 (Tpp − Tqq) may also arise from misalignments and mispointings of the gradiometer axes, and the related analysis and possible solutions are discussed in [46].
Recovering the SI units, we have Form previous experiments [14,20,21], the new PN parameter χ = 4hc/a M C S had already been constrained to be a small quantity χ ≤ 10 −4 . Π C S and Δ C S can then be expanded as Therefore, one has s C S = − 45G 2 cos 4 χ sin 4R aχ χ J MΨ sin φ cos Ψ 8c 2 a 3 R 3 ω +O(χ 2 ).
To give the estimation, we assume the orbital altitude to be 500 km and the mission life time about one year. After one year's accumulation, the total orbital cycle N in s C S will be 3.5 × 10 4 , and the secular signal will reach about 4.3χ mE. With proper data analysis methods employed, the total signal-to-noise ratio can be further amplified by a factor of the square root of the total cycles √ N . Therefore, for superconducting gradiometers with potential sensitivity better than 10 −2 mE/ √ Hz in low frequency band near 0.1 mHz [29], a rather strong constraint on the CS mass scale of the non-dynamical theory may in principle be obtained as M C S ≥ 10 −7 eV andθ ≤ 1 m.
For future gradiometers with optical readout 1 based on similar measurement schemes and techniques from the LPF mission, an even stronger bound may be expected M C S ≥ 10 −6 −10 −5 eV andθ ≤ 10 −1 −10 −2 m. Therefore, to conclude, the proposed experiment could naturally be incorporated into future missions that carrying highsensitive gradiometers as key payloads, which in principle would improve the current constraints of the CS modified gravity and add potential scientific objectives to such satellite gradiometry missions.