Constraint on Brans-Dicke theory from Intermediate/Extreme Mass Ratio Inspirals with self-force method Tong Jiang,1, ∗ Ning Dai,1, y Yungui Gong,1, z Dicong Liang,1, x 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 y [email protected] two-year observation of a 1:4M NS on a quasicircular in- z 3 Corresponding author. [email protected] spiral into a 10 M BH by LISA could give the constraint x 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 p 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 ,α p 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].