Constraint on Brans-Dicke theory from Intermediate/Extreme Mass Ratio Inspirals with self-force method
Tong Jiang,1, ∗ Ning Dai,1, † Yungui Gong,1, ‡ Dicong Liang,1, § and Chao Zhang1, ¶ 1School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China Intermediate/Extreme mass ratio inspiral (IMRI/EMRI) system provides a good tool to test the nature of gravity in strong field. We construct the self-force and use the self-force method to generate accurate waveform templates for IMRIS/EMRIs on quasi-elliptical orbits in Brans-Dicke theory. The extra monopole and dipole emissions in Brans-Dicke theory accelerate the orbital decay, so the observations of gravitational waves may place stronger constraint on Brans-Dicke theory. With a two-year observations of gravitational waves emitted from IMRIs/EMRIs with LISA, we can get the 5 most stringent constraint on the Brans-Dicke coupling parameter ω0 > 10 .
Since the first direct detection of gravitational wave orbital dynamics and the GW waveform of a compact (GW) event GW150914 [1,2] by the Laser Interferom- binary system different from those in GR [46–52], so BD eter Gravitational-Wave Observatory (LIGO) Scientific theory can also be probed by the observations of GWs Collaboration [3,4] and Virgo Collaboration [5], there [53–63]. However, the orbital evolution and gravitational have been tens of confirmed GW detection [6,7]. It is radiation from binary black holes (BBHs) are identical well known that in general relativity (GR) GWs prop- in GR and BD theory, so GW observations of BBHs are agate with the speed of light and have only two tensor unable to distinguish BD theory from GR [38, 46, 61]. By polarizations. However, six possible polarization states in using the simple approximate waveform template with GWs are allowed in general metric theory of gravity [8,9] the dipolar correction in the phase and the Fisher matrix and the number of polarization states depends on the par- method, it was estimated that the observation of a 0.7M ticular theory of gravity [10–14]. The detected GWs are neutron star (NS) on a quasicircular inspiral into a 3M useful to understand the nature of gravity and test Ein- black hole (BH) with a signal-to-noise (SNR) of 10 by stein’s GR in strong-field and nonlinear regions [15–18]. LIGO/Virgo detectors could give the constraint ωBD & The observation of GW170817 and its electromagnetic 2000 [53]. Using the Bayesian inference method, the GW counterpart GRB170817A constrained the speed of GWs event GW190426 152155 of a possible 1.5M NS/5.7M −15 −16 as −3 × 10 < cgw/c − 1 ≤ 7 × 10 [19] and this BH binary gave the constraint ωBD & 10 [64]. measurement on the propagation speed of GWs was al- If the BHs are more massive, then the frequency of ready used to exclude some alternative theories of gravity emitted GWs is in the mHz band and the GWs should [20–31]. be measured by space-based GW detectors like the Laser One of the simple alternative theories of gravity is Interferometer Space Antenna (LISA) [65, 66], TianQin Brans-Dicke (BD) theory of gravity [32, 33]. In BD the- [67] and Taiji [68]. One particular interest in space-based ory, the BD scalar field ϕ not only takes the role of 1/G GW detectors is the stellar-mass BHs or NS captured into but also mediates gravity and excites the scalar breath- inspiral orbits around massive BHs (MBHs), the Inter- ing mode in GWs. The most stringent constraint on mediate/Extreme mass ratio inspiral (IMRI/EMRI). The BD theory comes from the Cassini measurements of the mass ratio between the MBHs and the compact stellar ob- Shapiro time delay in the solar system [34] and the re- ject is about 102-104 : 1 for IMRIs and 104 : 1 for EM- sult is ω > 40000 [35]. For a binary system, the or- & BD RIs. In EMRIs, the timescale on the modification of the bital period of the system will decrease due to the loss orbit due to the back-reaction from gravitational radia- of energy by the emission of GWs. In BD theory [36], tion is much larger than the orbital period, so it takes the arXiv:2107.02700v2 [gr-qc] 7 Jul 2021 the extra dipolar emission channel of GWs can further compact object (CO) the last few years to inspiral deep decrease the orbital period of a binary system [36, 37], inside the strong field region of the MBH with a speed of so the measurement of the change of the orbital period a significant fraction of the speed of light and there are of a binary can be used to constrain BD theory [36–44]. 105-106 GW cycles in the detector band [69]. The emit- By using the measurement of the orbital decay from the ted GWs from (IMRIs/EMRIs) encode rich information pulsar-white dwarf binary PSR J1738+0333, the BD pa- about the spacetime geometry around the MBH and they rameter ω was constrained to be ω > 25000 [45]. 0 0 can be used to confirm whether the MBH is a Kerr BH The extra energy loss in BD theory makes both the predicated by GR. Using the Fisher matrix method, it was shown that the observation of a 1.4M NS on a qua- 3 sicircular inspiral into a 10 M BH with a SNR of 10 by ∗ [email protected] LISA could give the constraint ωBD & 24000 [58], and a † [email protected] two-year observation of a 1.4M NS on a quasicircular in- ‡ 3 Corresponding author. [email protected] spiral into a 10 M BH by LISA could give the constraint § 5 [email protected] ωBD > 3×10 [59]. Including the spin-spin coulpling and ¶ chao [email protected] small eccentricity, the constraint from a 1.4M NS inspi- 2 √ 3 ralling into a 10 M BH with a SNR of 200 by LISA elliptical orbits in the framework of GR and BD theory. became ωBD & 6944 [60]. Because BD theory modifies We then discuss the constraint on ωBD with LISA. With gravity in the weak field and the deviation in the en- two-year observation of IMRI/EMRI, we derive the con- 5 ergy flux are largest at small velocities, the constraints straint ω0 > 10 , which is more stringent than current on ωBD from EMRIs on quasicircular orbits are worse solar system tests. The constraint obtained in this paper than those derived from comparable-mass binaries [61]. should be more reliable than those derived in [58–60] with For EMRIs, half the total energy is radiated earlier than the approximate analytic waveform in frequency domain 10 years before the final plunge [70] and monopole radi- and the Fisher matrix approximation. In this paper we ation appears when the eccentricity is nonzero [38, 57], adopt natural units c = G = 1. so the early inspiral of EMRIs on quasi-elliptical orbits For EMRI systems, at the zeroth order of approxima- may place stronger constraint on BD theory. Therefore, tion the small compact object moves along the geodesic it is interesting to consider EMRIs when the small CO of the central MBH. The bound geodesic is a Keplerian moves slowly in a quasi-elliptical orbit. orbit parameterized by χ as The long inspiral time of EMRIs makes the genera- tion of accurate template waveforms for matched filter- pm r(χ) = 1 , (1) ing with the numerical relativity method computation- 1 + e cos (χ − w) ally impossible. However, the problem can be approached based on the expansion in mass ratio. To the lowest or- der, the small CO can be treated as a point like test parti- where m1 is the mass of the central MBH, the parameter cle moving in the geodesics of the central MBH. To higher χ runs from 0 to 2π over one radial cycle, w is the value orders, the gravitational field, the internal structure of of χ at periapsis, e is the orbital eccentricity and p is the small CO and the back-reaction of gravitational ra- the semilatus rectum. The self-force including the back diation are treated as perturbations. The waveforms are reaction of radiative GW and the interaction between calculated by solving the Teukolsky equation [71] and the small compact object and its own gravitational per- summing all the multipole modes. The Teukolsky-based turbation, is considered as a perturbation acting on the waveforms are computationally expensive. In order to geodesic of the central MBH. With the self-force, the or- quickly derive the equation of motion (EoM) and the bital parameters p, e and w evolve with time and we can waveform template for EMRIs, the kludge models in- get more accurate template waveforms. Here we mainly cluding the analytical kludge (AK) model [70] and the discuss the method of Pound and Poisson which uses the numerical kludge (NK) model [72] have been proposed. PN theory to construct the SF up to 2.5PN [88]. The AK model assumes that the small CO moves on a Now we apply the above SF method to the BD theory Keplerian orbit with relativistic corrections such as peri- with the action [32, 36] apsis precession, Lense-Thirring precession and inspiral from radiation reaction given by analytic post-Newtonian Z −1 −1 ,α √ 4 (PN) evolution equations. It is extremely quick to cal- S =(16π) ϕR − ϕ ω(ϕ)ϕ ϕ,α −gd x culate, but it dephases relative to true waveforms within (2) hours. The augmented AK model improved the accuracy + Sm(Ψ, gαβ), of waveform template [73, 74]. The NK model combines Kerr geodesic with PN orbital evolution which is caused where ϕ is the BD scalar field, the coupling parameter ω by radiation reaction of GWs, and numerically integrates is a function of ϕ and Sm is the matter action in which the Kerr geodesic equations along the inspiral trajectory. the matter field Ψ couples to the metric gαβ only but It is more accurate and computationally expensive than the mass of a self-gravitating body depends on the BD the AK model. scalar field. We denote ϕ0 as the vacuum expectation On the other hand, the interaction between the value of ϕ and ω0 = ω(ϕ0). To parametrize the sensi- small CO and its own gravitational perturbation can tivity of the body’s binding energy on the background be thought as an effective gravitational self-force driving scalar field, we introduce the the sensitivity of body A as the radiative evolution of the geodesic orbit of the cen- sA = d ln mA(ϕ)/d ln ϕ0. In the original BD theory [32], tral MBH. The perturbative force includes the radiation- ω(φ) = ωBD, the effective Newtonian gravitational cou- reaction of GWs and the gravitational effect caused by pling constant measured by Cavendish-like experiments the small CO [75, 76]. There were some developments is G = (4 + 2ωBD)/[φ0(3 + 2ωBD)], and the stationary, on the construction of self-force [77–87]. To integrate asymptotically flat BHs in vacuum in BD theory are the the equations of motion that govern accelerated motion BHs of GR. due to self-force in Schwarzschild spacetime, the method For a compact binary, the relative acceleration between of osculating orbits was proposed [88]. The self-force the two COs up to the 2.5PN in BD theory is [46] method can be applied in strong field regions and can generate waveform templates for EMRIs with a good accuracy [88–90]. In this paper, we use the self-force d2~x αm 2 = − 2 (ABDn − BBDv), (3) method to generate waveforms for EMRIs with quasi- dt rh 3 where the orbits in GR and BD theory are almost identical. Af- 8 αm ter 18 months of inspirals, the orbits in BD theory with ABD = 1 − A1PN − A2PN − η r˙h(A1.5PN + A2.5PN), ω0 = 50000 start to deviate from those in GR and the or- 5 rh bits in BD theory with ω = 2000 are different from those 8 αm 0 BBD = (B1PN + B2PN)r ˙h − η (B1.5PN + B2.5PN), in GR. The main difference in the orbital evolution can 5 rh be manifested by the orbital phase accumulation. In Fig. (4) 2, we show the evolutions of the orbital phase difference ∆φ = φ − φ between GR and BD theory with dif- r ≡ |~x| = |~x − ~x |, the subscript h means that we BD GR 1 2 ferent ω . As expected, the phase difference increases as use the harmonic gauge in the PN expansion, the har- 0 ω becomes smaller which is consistent with the results monic coordinate is related to Schwarzschild coordinate 0 of the orbital evolution. From Fig.2, we see that ∆ φ by the transformation r = r − m , n = ~x/r , m = h 1 h oscillates with time. After around one-year inspiral, the m + m , η = m m /m2, v = v − v ,r ˙ = dr /dt, 1 2 1 2 1 2 h h phase difference is larger than 40 rad for ω = 2000; the α = 1 − ζ + ζ(1 − 2s )(1 − 2s ), ζ = 1/(4 + 2ω ), the 0 1 2 0 phase difference oscillates between 6 rad and 12 rad for variables A and B through 2.5PN are given in the ap- ω = 10000; and the maximum phase difference reaches pendix. For the equations of motion through 2PN and 0 5 rad for ω = 50000. The accumulated phase difference 3PN, please refer [47, 48, 50]. From Eq. (3) we construct 0 will be manifested in the GW waveform which provides the SF and then we solve the osculating orbits for the us tool to constrain BD. orbital parameters to see whether the orbital evolutions in GR and BD theory are distinguishable. Since the or- 50 ω = bital evolution of BBHs cannot distinguish GR from BD 0 2000 ω0=10000 ω = theory, so we consider a binary with a NS and a BH. 40 0 50000 The mass and the sensitivity of the NS are m2 = 1.3M and s2 = 0.2 respectively, and the mass ratio between −4 30 the NS and the BH is q = m2/m1 = 10 . The initial rad conditions are as follows: at the initial time t0 = 0, the Δϕ / 20 dimensionless semilatus rectum p0 = 80, the eccentricity e0 = 0.8, the orbital parameter at periapsis w0 = 10, and the phase at periapsis φ0 = 0. The results for the orbital 10 evolution are shown in Fig.1.
0 0 200 400 600 800 1000 1200