 if there are analytic methods to investigate gravitational waves from EMRIs.

The post-Newtonian (PN) approximation to the Einstein equations is a standard method to compute gravitational waveforms from inspiraling compact binaries [@post_Newton]. In the PN approximation to the compact binary system, one assumes that the velocities of the binary are much smaller than the speed of light, $v/c\ll 1$. In the standard PN approximation, the amplitude of gravitational waves and the orbital phase are, respectively, derived up to 3PN and 3.5PN, i.e. $v^6$ and $v^7$ beyond the leading order for the non-spinning compact binaries in quasi-circular orbits [@DJS01; @BDE04; @BFIJ02; @BDEI04; @K07; @BFIS08; @Favata09]. (Note that the 3.5PN amplitudes for $(\ell,m)=(2,2),(3,3)$ and $(3,1)$ modes are derived in Refs. [@FMBI2012; @FBI2014].) For the case of the spinning compact binaries in quasi-circular orbits, spin-orbit effects in the orbital phase are derived up to 4PN [@MBBB2014]. Spin-spin effects in the orbital phase are derived up to 2PN [@Gergely1999; @Gergely2000; @BVG2005; @RBK2009].

Using the PN approximation in the black hole perturbation theory, high PN order expressions for gravitational waves can be obtained more systematically than using the standard PN approximation [@ST]. The energy flux to infinity up to 5.5PN (4PN) for the case of a test particle in circular orbits around a Schwarzschild (Kerr) black hole was derived in Ref. [@TTS] (Ref. [@TSTS]) by solving the Teukolsky equation [@Teukolsky1973], which is the fundamental equation in the black hole perturbation theory. More recently, very high PN order expressions in the energy flux to infinity and gravitational waveforms for a test particle in circular orbits around a Schwarzschild black hole were derived up to 22PN [@14PN; @22PN] using a more systematic method to solve the Teukolsky equation [@MST; @MSTR]. It was shown that dephases between 22PN waveforms and very highly accurate waveforms during two-year inspirals can
if there are analytic methods to investigate gravitational waves from EMRIs. 

 The post - Newtonian (PN) approximation to the Einstein equation is a standard method acting to compute gravitational waveforms from inspiraling compendious binaries   [ @post_Newton ]. In the PN estimate to the compact binary system, one assume that the velocities of the binary are a lot belittled than the amphetamine of light, $ v / c\ll 1$. In the standard PN approximation, the amplitude of gravitational wave and the orbital phase are, respectively, derived up to 3PN and 3.5PN, i.e. $ v^6 $ and $ v^7 $ beyond the leading order for the non - spinning compact binaries in quasi - round orbits   [ @DJS01; @BDE04; @BFIJ02; @BDEI04; @K07; @BFIS08; @Favata09 ]. (Note that the 3.5PN amplitudes for $ (\ell, m)=(2,2),(3,3)$ and $ (3,1)$ mode are derived in Refs.   [ @FMBI2012; @FBI2014 ] .) For the case of the spin compact binary star in quasi - circular orbits, spin - orbit effects in the orbital phase are derive up to 4PN   [ @MBBB2014 ]. Spin - spin effects in the orbital phase are derived up to 2PN   [ @Gergely1999; @Gergely2000; @BVG2005; @RBK2009 ]. 

 Using the PN estimate in the bootleg trap perturbation theory, high PN order expressions for gravitational waves can be obtained more systematically than using the standard PN estimate   [ @ST ]. The department of energy flux to infinity up to 5.5PN (4PN) for the case of a test atom in round orbits around a Schwarzschild (Kerr) black hole was derived in Ref.   [ @TTS ] (Ref.   [ @TSTS ]) by solving the Teukolsky equation   [ @Teukolsky1973 ], which is the fundamental equation in the bootleg hole perturbation theory. More recently, very high PN order expressions in the energy flux to infinity and gravitational waveforms for a test particle in round orbits around a Schwarzschild bootleg hole were derived up to 22PN   [ @14PN; @22PN ] using a more systematic method to solve the Teukolsky equality   [ @MST; @MSTR ]. It was shown that dephases between 22PN waveforms and very highly accurate waveforms during two - year inspirals can
 if there are analytic methuds to investigcre gratitatiohal waver from EMRIs.

The post-Newtoniai (PN) appriximation to the Einstdin equatpons is a staidard method to rkmpute nxavitzbionan waveforms frok inspiralhng compact bitafizs [@post_Newton]. In the PN approximation to the clmpact binary fysttm, jne zssumes that the velocities of the binary are much smakler than the speed of ligjt, $v/f\ll 1$. In the standagd PN approzimaeuon, the amplktude of ggcvitational waves and the orbital phase ard, res'ectively, dgxuvef up to 3PN aid 3.5PN, p.e. $v^6$ and $v^7$ benpnd tha leadimg order for tme noi-spibning compact binariev in quasi-circular orbits [@DJV01; @YDE04; @BFIJ02; @BDEI04; @K07; @BFIS08; @Dacata09]. (Tote thag thd 3.5PH empmituded fkr $(\ell,m)=(2,2),(3,3)$ ans $(3,1)$ modes arw derived in Refs. [@FMNI2012; @DBI2014].) For the czse of tre spinning compact binaries in quasi-cigculzr orbits, spin-orbit effwcts in the orbital pjase are qerived up to 4PN [@MBBB2014]. Spin-spin effects in the orbidal piare cvc defuvfd up to 2PN [@Gergely1999; @Gergely2000; @BVG2005; @RBK2009].

Using the PG akprpximation in tme black hole pertirhayyon theory, hieh PN ordsr expressions for gravitwtionql waves san ne obtained more systematicqlly than uspng rhe standard PN ap'roximation [@SC]. The gnergy flux to infinity up to 5.5PN (4LN) for the fase of a gest particle in cigculdr orbits around a Schwarzfchild (Kecr) blcck hole was derivqd in Ref. [@THS] (Ren. [@DSTS]) by solving thf Teunonsky equatlon [@Teukolsky1973], which is the fundamental equatiom hn nhe black hole perturbation theory. More rgcently, vzry hieh PN ordeg expressmons in the qnergy flux tm infinity anv gravitaeionql wqveformr for a test patticle in circular irbits around a Scmwarzrdhild black holz qere derived up to 22PN [@14[N; @22PI] usigc a more sysdemagic kethoa to solve bhe Teulolsky equation [@MST; @MVTR]. Jt was shown that ceihases bejween 22PN raveforms and very highly accurwte weveforos dutind two-year inspirals can
if there are analytic methods to investigate from The post-Newtonian approximation to the method compute gravitational waveforms inspiraling compact binaries In the PN approximation to the binary system, one assumes that the velocities of the binary are much smaller the speed of light, $v/c\ll 1$. In the standard PN approximation, the amplitude gravitational and orbital are, respectively, derived up to 3PN and 3.5PN, i.e. $v^6$ and $v^7$ beyond the leading order the non-spinning compact binaries in quasi-circular orbits [@DJS01; @BFIJ02; @BDEI04; @K07; @BFIS08; (Note that the 3.5PN amplitudes $(\ell,m)=(2,2),(3,3)$ $(3,1)$ modes derived Refs. @FBI2014].) For the of the spinning compact binaries in quasi-circular orbits, spin-orbit effects in the orbital phase are derived up 4PN [@MBBB2014]. in the phase derived to 2PN [@Gergely1999; @RBK2009]. Using the PN approximation in perturbation theory, high PN order expressions for gravitational can be more systematically than using the standard approximation [@ST]. The energy flux to infinity up 5.5PN (4PN) for the case of a test particle in circular orbits around a Schwarzschild hole was derived in [@TTS] (Ref. [@TSTS]) solving Teukolsky [@Teukolsky1973], is the equation in the black hole perturbation theory. More recently, very high order expressions in the energy flux to infinity and gravitational a particle in circular around a Schwarzschild black were up to 22PN [@14PN; a systematic the equation @MSTR]. It was shown dephases between 22PN waveforms and highly accurate waveforms during
 if there are analytic methods To investigAte grAviTatIoNal wAves From EMRIs.

The poST-NewTonian (PN) approximation tO the EInSTein EQuAtionS is a staNDaRD MetHoD tO coMpUTe GraviTatIonal waVeforms froM inSpIraling compaCT bInaries [@posT_NeWton]. In the PN aPprOximatIoN to THe comPacT binaRy systEM, one asSumes that ThE VelociTIes of thE BInAry aRe much smaller than THe SPeed of light, $v/c\lL 1$. In the StANdARD PN AppRoximation, ThE amplITude of gRAvITATioNAl waves and the Orbital phasE Are, RespecTiVelY, DeriveD up to 3pN ANd 3.5Pn, i.e. $v^6$ and $v^7$ beyOnd tHe leading Order fOR the non-SPinning CompacT biNarIes iN QuAsI-ciRcULar ORbIts [@djS01; @BdE04; @BFIJ02; @BDeI04; @k07; @BfIS08; @FaVata09]. (nOTE That The 3.5pN amPlituDes for $(\ell,m)=(2,2),(3,3)$ and $(3,1)$ ModEs arE DerIved iN Refs. [@fMBI2012; @fBi2014].) For tHe case Of the SpInning compact biNariEs in quasi-CirCuLar OrBits, sPIn-orbiT efFecTs in the Orbital PHasE aRE DErIved up to 4PN [@MBBB2014]. Spin-SpIN EfFects in tHe orbiTAl PhASe are derIvEd uP to 2Pn [@gErgelY1999; @GerGElY2000; @BVG2005; @RBK2009].

USing thE pN ApProximaTiOn in thE bLacK hoLe perTUrbaTion thEory, high pN ordER expressions foR Gravitational WAvES CaN Be obTaiNed more systEmatICallY thaN UsIng THe staNdard pN APpROximation [@ST]. The energY fLux to iNfiniTy up to 5.5PN (4PN) for The case of a TEST particlE in cIRcULar orbits arounD a SchWarzschild (kErr) black Hole wAs deriveD in Ref. [@TTS] (rEF. [@TSTS]) by sOlvIng The teuKOLsKy equation [@TeuKOLsky1973], WhIch is thE fuNdamentAl eQuaTioN in ThE black holE perturbAtIoN tHeOry. more rECently, veRy HigH Pn orDer exPRessioNs in tHe enErGy FLux To infinITy AND graViTaTionAl wAvEformS for A TesT particLe in circuLar ORbitS aRoUnd a SchWarzschild blaCk Hole were deRiVed Up to 22PN [@14pn; @22pN] using a More systematic method to sOLve the TEukOlsky EquaTion [@MST; @MStR]. IT was shOwn THat depHases bEtweeN 22Pn waVEForms AND vEry HiGhly accuraTE WavEformS dUrinG two-yeaR inspirals can
 if there are analytic met hods to in vesti gat e g ra vita tion al waves fromE MRIs .

The post-Newtonian(PN)ap p roxi m at ion t o the E i ns t e ineq ua tio ns is  a st and ard met hod to com put egravitationa l  w aveforms f rom  inspiraling  co mpactbi nar i es [@ pos t_New ton].I n thePN approx im a tion t o  the co m p ac t bi nary system, onea ss u mes that the v elocit ie s  o f the  bi nary are m uc h sma l ler tha n  t h e spe e d of light, $ v/c\ll 1$.I n t he sta nd ard PN app roxim at i on,  the amplit udeof gravit ationa l  wavesa nd theorbita l p has e ar e ,re spe ct i vel y ,der i ved  up to 3 PN  a nd 3. 5PN, i . e . $v ^6$  and  $v^7 $ beyond thelea ding ord er fo r the  non -s pinni ng com pactbi naries in quasi -cir cular orb its  [ @DJ S0 1; @B D E04; @ BFI J02 ; @BDEI 04; @K0 7 ; @ BF I S 0 8;  @Favata09]. (Note  t h a tthe 3.5P N ampl i tu de s  for $(\ el l,m )=(2 , 2 ),(3, 3)$a nd  $(3,1)$  modes ar ederived  i n Refs .[@F MBI 2012; @FBI 2014]. ) For th e cas e  of the spinni n g compact bin a ri e s  i n  qua si- circular or bits ,  spi n-or b it  ef f ectsin th eo rb i tal phase are deriv ed  up to  4PN[@MBBB2014].Spin-spine f f ects intheo rb i tal phase arederiv ed up to 2 P N [@Gerg ely19 99; @Ger gely2000; @ BVG2005;  @R BK2 009 ].U si ng the PN app r o xima ti on in t heblack h ole  pe rtu rba ti on theory , high P Nor de rexp ressi o ns for g ra vit at ion al wa v es can  be o btai ne dm ore  system a ti c a llyth an  usi ngth e sta ndar d  PN  approx imation [ @ST ] . Th een ergy fl ux to infinit yup to 5.5P N(4P N) for t he caseof a test particle in c i rcularorb its a roun d a Schwa rzs child(Ke r r) bla ck hol e was  d eri v e d inR e f.  [@ TT S] (Ref. [ @ T STS ]) by  s olvi ng theTeukolsky equation [@T eukolsky1973] , w hich i sthe fu n dam en t ale q uation in the b lack holepe r tu rbation th e ory .More re cently,  very high PN  order ex pressions  i n th e ene rgy flux t o infini ty and gr a vitat i on al wa vef orms f or  atestpartic l e i n cir cularor bits a round  a  Schwarz schild black hole werederive d upto22PN [@14 PN; @22 PN] using  a m ore system ati c m ethod  to solve  the Te uko l sky e quat i on [@MST; @M STR ] .  I t was shown t h atdepha ses betwee n 22 PN waveforms andv ery highly acc urat e wav efo r ms d ur ing two-year i nsp ir a l s can
 if_there are_analytic methods to investigate_gravitational waves_from_EMRIs.

The post-Newtonian_(PN)_approximation to the_Einstein equations is_a standard method to_compute gravitational waveforms_from_inspiraling compact binaries [@post_Newton]. In the PN approximation to the compact binary system, one assumes_that_the velocities_of_the_binary are much smaller than_the speed of light, $v/c\ll_1$. In_the standard PN approximation, the amplitude of gravitational_waves_and the orbital_phase are, respectively, derived up to 3PN and 3.5PN,_i.e. $v^6$ and $v^7$ beyond the_leading order for_the_non-spinning_compact binaries in quasi-circular_orbits [@DJS01; @BDE04; @BFIJ02; @BDEI04; @K07; @BFIS08;_@Favata09]. (Note that the 3.5PN amplitudes_for $(\ell,m)=(2,2),(3,3)$ and $(3,1)$ modes are derived_in Refs. [@FMBI2012; @FBI2014].) For the case_of the spinning compact binaries_in quasi-circular_orbits, spin-orbit effects in the_orbital phase are_derived up_to 4PN [@MBBB2014]. Spin-spin_effects in the orbital phase are_derived up to_2PN [@Gergely1999; @Gergely2000; @BVG2005; @RBK2009].

Using the PN_approximation_in the black_hole_perturbation_theory, high_PN order expressions_for_gravitational waves_can_be obtained more systematically than using_the_standard PN approximation [@ST]. The energy flux to_infinity up to 5.5PN_(4PN)_for the case of_a test particle in circular_orbits around a Schwarzschild (Kerr) black_hole was_derived in_Ref. [@TTS] (Ref. [@TSTS]) by solving the Teukolsky equation [@Teukolsky1973], which is the fundamental_equation in the black hole perturbation_theory. More recently, very_high PN_order_expressions in the_energy_flux to_infinity and gravitational waveforms for a test_particle in_circular orbits around a Schwarzschild black_hole were derived up_to_22PN [@14PN; @22PN] using a more systematic_method to solve the Teukolsky equation [@MST;_@MSTR]. It was shown that_dephases_between_22PN waveforms and very highly_accurate waveforms during two-year inspirals can