This document is part of https : // dcc.ligo.org/LIGO - T1500554/public

Load FeynCalc

<<HighEnergyPhysics`fc`

Loading FeynCalc from /home/riccardo/.Mathematica/Applications/HighEnergyPhysics

FeynCalc 8.2.0 For help, type ?FeynCalc, open FeynCalcRef8.nb or visit www.feyncalc.org

Loading FeynArts, see www.feynarts.de for documentation

FeynArts 3.10 patched for use with FeynCalc

Tensor definitions

V3[a_,b_]:=FourVector[a,b,Dimension->3];
δ[i_,j_]:=MetricTensor[i,j,Dimension->3];
ENS[x_]:=Expand[Normal[Simplify[x]]];

General definitions

η=m1 m2/(m1+m2)^2;
M=m1+m2;
(*Efac=-1/2η M;*)
δe=(m1-m2)/(m1+m2);
κp=κ1+κ2;
κm=κ1-κ2;
SS=SP[S1,S1]+SP[S2,S2]+2SP[S1,S2];
ΣΣ=(m1+m2)^2(SP[S2,S2]/m2^2+SP[S1,S1]/m1^2-2SP[S1,S2]/(m1 m2));
SΣ=(m1+m2)(SP[S2,S2]/m2-SP[S1,S1]/m1+SP[S1,S2](1/m2-1/m1));
projSS=(SS-(SP[S1,L]SP[S1,L]+SP[S2,L]SP[S2,L]+2SP[S1,L]SP[S2,L]))/2;
projΣΣ=(ΣΣ-(m1+m2)^2(SP[S2,L]SP[S2,L]/m2^2+SP[S1,L]SP[S1,L]/m1^2-2SP[S1,L]SP[S2,L]/(m1 m2)))/2;
projSΣ=(SΣ-(m1+m2)(SP[S2,L]SP[S2,L]/m2-SP[S1,L]SP[S1,L]/m1+SP[S1,L]SP[S2,L](1/m2-1/ m1)))/2;
projSSnn=(SP[S1,n]SP[S1,n]+SP[S2,n]SP[S2,n]+2SP[S1,n]SP[S2,n]);
projΣΣnn=(m1+m2)^2(SP[S2,n]SP[S2,n]/m2^2+SP[S1,n]SP[S1,n]/m1^2-2SP[S1,n]SP[S2,n]/(m1 m2));
projSΣnn=(m1+m2)(SP[S2,n]SP[S2,n]/m2-SP[S1,n]SP[S1,n]/m1+SP[S1,n]SP[S2,n](1/m2-1/ m1));
nS=SP[n,S1]+SP[n,S2];
nΣ=(m1+m2)(SP[n,S2]/m2-SP[n,S1]/m1);

ω - γ relationship

Eq. (3) of arXiv : astro - ph/0504538
γ=x(1+x γ1PN +x^(3/2)((M/m1)-1/3)SP[S1,L]+...+x^2/2(SP[S1,S2]-3SP[S1,L]SP[S2,L]))

From the equation a^i = -(ω^2 r) n^i
ω^2 r = M/r^2 (1 + (Mω)^2/(M/r)^2 g1PNa + M/r g1PNb + (Mω) gSO1p5PN + M^2/r^2  + (M/r)^3 gSS3PN +...) -> x^3 = γ^3 (1 + x^3/(M/r)^2 g1PNa + M/r g1PNb + x^(3/2) gSO1p5PN + γ^2 gSS2PN + 
gSS is given by (-1 times) the rhs of (A3).n of arXiv:1501.01529

2 PN

γ1PN=(1-η/3);

From eq.(A3) of arXiv : 1501.01529

gSS2PN=-1/4 (SS (-6 κp-12)+SΣ (-6 δe κp-12 δe+6 κm)+ΣΣ (3 δe κm-3 κp+η (6 κp+12))+projSS (30 κp+60)+projSΣ (30 δe κp+60 δe-30 κm)+projΣΣ (15 κp-15 δe κm+η (-30 κp-60))+projSS (-12 κp-24)+projSΣ 2 (-6 δe κp-12 δe+6 κm)+projΣΣ (6 δe κm-6 κp+η (12 κp+24)));
gSS2PNinst=-1/4 (SS (-6 κp-12)+SΣ (-6 δe κp-12 δe+6 κm)+ΣΣ (3 δe κm-3 κp+η (6 κp+12))+nS^2 (30 κp+60)+nS nΣ (30 δe κp+60 δe-30 κm)+nΣ^2 (15 κp-15 δe κm+η (-30 κp-60))+nS^2 (-12 κp-24)+2 nS nΣ (-6 δe κp-12 δe+6 κm)+nΣ^2 (6 δe κm-6 κp+η (12 κp+24)));
(*Now use that M^2ω^2 = M^3/r^3(1-M^2/(4r^2)(χ^2+..))*)
γSS2PN=-1/3gSS2PN;
γSS2PNinst=-1/3gSS2PNinst;

Checks

From (3.32) of arXiv : 1501.01529

γSS2PNcheck=SS(-κp/2-1)+SΣ(-δe κp/2-δe +κm/2)+ΣΣ(δe κm/4-κp/4+η(κp/2+1));
ENS[(γSS2PN-γSS2PNcheck)/.SP[L,S1]^2->SP[S1,S1]/.SP[L,S2]^2->SP[S2,S2]/.SP[L,S1]SP[L,S2]->SP[S1,S2]]
ENS[(γSS2PNinst-γSS2PNcheck)/.SP[S1,n]->0/.SP[S2,n]->0]

3 PN

g1PNa=1+3η(*v^2*);
g1PNb=-2(2+η)(*γ*);
g1p5PNSO=-SP[L,S1] (3M/m1-1)-SP[L,S2](3M/m2-1);
g2PNa=3/4(12+29η)(*γ^2*);
g2PNb=η(3-4η)(*v^4*);
g2PNc=-η/2(13-4η)(*γ v^2*);
γ2PN=(1-65η/12);

γx3[x_,γ_,ϵ_]:=(1+ϵ^2(x^3/γ^2g1PNax+γ g1PNbx)+ϵ^3 x^(3/2)g1p5PNSOx+ϵ^4γ^2gSS2PNx+ϵ^6γ^2(x gSS3PNax+γ gSS3PNbx))^(-1/3);
TotalN[list_,n_]:=Module[{max,out,idx},(max=Min[n,Length[list]];out=0;For[idx=1,idx≤max,(out+=list[[idx]];idx++;)];Return[out];)]

γL={1};
pNMaxOrder=3;
(*Print["γL start: ",γL];*)
Module[{tmp,idx,jdx,add,γxList},(γxList={1};
For[idx=1,idx≤2*pNMaxOrder,
(AppendTo[γxList,SeriesCoefficient[γx3[x,γ,ϵ],{ϵ,0,idx}]];
tmp=0;
For[jdx=0,jdx<idx,(
add=0;If[γL[[jdx+1]]=!=0,
(*Print["Inizio jdx ",jdx," idx ",idx," trova γL[",jdx+1,"]=",γL[[jdx+1]]];*)
add=ENS[SeriesCoefficient[γxList[[idx+1-jdx]]/.γ->(x TotalN[γL,jdx+1])/.x->ϵ^2x,{ϵ,0,idx}]];
(*Print["  finito jdx ",jdx," idx ",idx," add->",add];*),(*Print["  Salto jdx ",jdx];*)];
tmp+=add;jdx++;)];
AppendTo[γL,ENS[tmp]];idx++;)])];
Print["γL: ",γL]

ENS[SeriesCoefficient[1/(Total[γL])^2,{x,0,5/2}]]

γSO2p5PNm2=ENS[2/9g1p5PNSO(g1PNa-2g1PNb)];
γSO2p5PN=ENS[1/9g1p5PNSO(2g1PNa+5g1PNb)];

ENS[SeriesCoefficient[1/(Total[γL])^2,{x,0,3}]]

ENS[SeriesCoefficient[(Total[γL])^2,{x,0,3}]]

γSS3PNm2=ENS[2/3(g3PNSSa+g3PNSSb)-1/9g1p5PNSO^2-2/9(g1PNa+4g1PNb) gSS2PN];
γSS3PN2=ENS[-2/3(g3PNSSa+g3PNSSb)+5/9g1p5PNSO^2+2/9(5g1PNa+8g1PNb) gSS2PN];

g3PNSSa=-1/8(projSS(12δe κm-36κp-48+η(-72 κp-144)+6δe κm-6κp+84+60κp+120+η(180κp+360))+projSΣ(24κm -24δe κp+η(-36δe κp-72δe+12κm)-24δe κp-24δe +24κm+η(-36δe κp-72δe +12κm)-12δe κp+132δe+12κm+η(-24κm)+60δe κp+120δe-60κm+η(180 δe κp+360δe-180κm))+projΣΣ(24δe κm-24κp+24+η(24δe κm+24κp)+η^2(72κp+144)+6δe κm-6κp+48+η(-6δe κm+18κp-180)+30κp-30δe κm+η(-90δe κm+30κp-120)+η^2(-180κp-360))+SS(-6δe κm-6κp-48+η(-36κp-72))+SΣ(-72δe +η(-36δe κp-72δe+60κm))+ΣΣ(-24+η(24δe κm-24κp+96)+η^2(36κp+72)));
g3PNSSb=-1/4(projSS(-24δe κm+72κp+164+η(36κp+72)+48δe κm-192κp-420+η(-96κp-192))+projSΣ(48δe κp+72δe-48κm+η(18δe κp+36δe+30κm)+48δe κp+84δe-48κm+η(18δe κp+36δe+30κm)-240δe κp-396δe+240κm+η(-96δe κp-192 δe-96κm))+projΣΣ(48κp-48δe κm+η(6δe κm-102κp-148)+η^2(-36κp-72)+120δe κm-120κp+η(240κp+372)+η^2(96κp+192))+SS(-8δe κm+40κp+72+η(20κp+40))+SΣ(48δe κp+72δe-48κm+η(20δe κp+40δe+12κm))+ΣΣ(24κp-24δe κm+η(-2δe κm-46κp-72)+η^2(-20κp-40)));
γSS3PN=ENS[-1/3(g3PNSSa+g3PNSSb)+2/9g1p5PNSO^2+1/9(4 g1PNa+7g1PNb) gSS2PN];

Check

ENS[γ1PN-(-1/3(g1PNa+g1PNb))]

Kidder 2.2 b 9506022, minus from collecting Newtonian term in front

g1PNcheck=(1+3 η)-2 (2+η);
ENS[g1PNcheck-(g1PNa+g1PNb)]
(*Kidder 2.2c 9506022, minus from collecting Newtonian term in front*)
g1p5PNSOcheck=-5 (SP[S1,L]+SP[S2,L])-3 (m1-m2)(SP[S2,L]/m2-SP[S1,L]/m1);
(*A4 of arXiv:1501.01529*)
ENS[g1p5PNSOcheck-g1p5PNSO]

ENS[γ2PN-(-1/3(g2PNa+g2PNb+g2PNc)+1/3g1PNb(g1PNa+g1PNb))]

From eq.(3) of arXiv : astro - ph/0504538

γSO1p5PN=((M/m1)-1/3)SP[S1,L]+((M/m2)-1/3)SP[S2,L];
ENS[γSO1p5PN-(-1/3g1p5PNSO)]

γSS3PNcheck=SS(-11/12δe κm-11κp/12+14/9+η(-κp/6-1/3))+ΣΣ(1+η(δe κm-κp-2)+η^2(κp/6+1/3))+SΣ(5δe/3+η(-δe κp/6-δe/3+23κm/6));
ENS[(γSS3PN-γSS3PNcheck)/.SP[L,S1]^2->SP[S1,S1]/.SP[L,S1]SP[L,S2]->SP[S1,S2]/.SP[L,S2]^2->SP[S2,S2]]

Energy S1^2

1 PN

ENa=-x^2/γ^2+2γ/x;
E1PNa=-3/4(1-3η)(*x^4/γ^4*);
E1PNb=-(3+η)(*x/γ*);
E1PNc=-1(*γ^2/x^2*);
E1PN=E1PNa+E1PNb+E1PNc+4γ1PN;

1.5 PN

E1p5PN0=-2SP[L,S1](-1+M/m1)-2SP[L,S2](-1+M/m2);
E1p5PN=E1p5PN0+4γSO1p5PN;
ENS[SeriesCoefficient[E1p5PN,{SP[L,S1],0,1}]/.m2->MT-m1]
ENS[SeriesCoefficient[E1p5PN,{SP[L,S2],0,1}]/.m1->MT-m2]

2 PN

3.30 of arXiv : 1501.01529 v1

Energy2PN0=-1/2(SS(-2κp-4)+SΣ(-2δe κp-4δe+2κm)+ΣΣ(δe κm-κp+η(2κp+4))+projSS(6κp +12)+projΣΣ(3κp-3δe κm+η(-6κp-12))+projSΣ(6δe κp+12δe -6κm));
Energy2PN[S1_,S2_,m1_,m2_]:=Energy2PN0+4γSS2PN
(*The Newtonian Energy is -1/2 η M x-2γss=-1/2ηM(x+4 γSS)*)
Energy2PNnoquad[S1_,S2_,m1_,m2_]:=(-3/η SP[L,S1]SP[L,S2]+1/η SP[S1,S2]);
Energy2PNk1=ENS[SeriesCoefficient[Energy2PN[S1,S2,m1,m2],{κ1,0,1}]];
Energy2PNk2=ENS[SeriesCoefficient[Energy2PN[S1,S2,m1,m2],{κ2,0,1}]];
ENS[(Energy2PN[S1,S2,m1,m2]/.κ1->0/.κ2->0)-Energy2PNnoquad[S1,S2,m1,m2]]
ENS[Energy2PNk1-((M/m1)^2(-3/2SP[L,S1]^2+1/2SP[S1,S1]))]
ENS[Energy2PNk2-((M/m2)^2(-3/2SP[L,S2]^2+1/2SP[S2,S2]))]

Energy2PNinst=-1/2(SS(-2κp-4)+SΣ(-2δe κp-4δe+2κm)+ΣΣ(δe κm-κp+η(2κp+4))+(nS)^2(6κp +12)+(nΣ)^2(3κp-3δe κm+η(-6κp-12))+nS nΣ(6δe κp+12δe -6κm))+4γSS2PNinst;
Energy2PNinstnoquad[S1_,S2_,m1_,m2_]:=(6/η SP[n,S1]SP[n,S2]-2/η SP[S1,S2]);
ENS[(Energy2PNinst/.κ1->0/.κ2->0)-Energy2PNinstnoquad[S1,S2,m1,m2]]

Checks

ENS[(Energy2PN[S1,S2,m1,m2]-(κ1 Energy2PNk1+κ2 Energy2PNk2+Energy2PNnoquad[S1,S2,m1,m2]))]

(3.33) of arXiv : 1501.01529v1

E2PNcheck[S1_,S2_,m1_,m2_]:=(SS(-κp-2)+SΣ(-δe κp-2δe+κm)+ΣΣ(δe κm/2-κp/2+η(κp+2)))
ENS[(Energy2PN[S1,S2,m1,m2]-E2PNcheck[S1,S2,m1,m2])/.SP[L,S1]^2->SP[S1,S1]/.SP[L,S2]^2->SP[S2,S2]/.SP[L,S1]SP[L,S2]->SP[S1,S2]]

E2PNA=(-2/η)SP[S1,S2]+κ1(-(M/m1)^2)SP[S1,S1]+κ2(-(M/m2)^2)SP[S2,S2];
ENS[(E2PNcheck[S1,S2,m1,m2])-E2PNA]

ENS[(Energy2PNinst-E2PNcheck[S1,S2,m1,m2])/.SP[L,S1]^2->SP[S1,S1]/.SP[L,S2]^2->SP[S2,S2]/.SP[L,S1]SP[L,S2]->SP[S1,S2]/.SP[n,S1]->0/.SP[n,S2]->0/.κ1->0/.κ2->0]

Eq. 5.14e of 0812.4413 (this paper consider BHs, so κ1=κ2=1)

S0S0=SP[S1,S1](1+m2/m1)^2+SP[S2,S2](1+m1/m2)^2+2SP[S1,S2](1+m2/m1)(1+m1/m2);
projS0S0=(S0S0-((1+m2/m1)^2SP[S1,L]SP[S1,L]+(1+m1/m2)^2SP[S2,L]SP[S2,L]+2(1+m2/m1)(1+m1/m2)SP[S1,L]SP[S2,L]))/2;
Energy2PNRacine=ENS[S0S0-3projS0S0+4γSS2PN];
ENS[(Energy2PN[S1,S2,m1,m2]-Energy2PNRacine)/.κ1->1/.κ2->1]

Eq.(21) of arXiv : gr - qc/9709032 v1+eq.(6) of astro-ph/0504538v2

Energy2PNPoi=SP[S1,S2]/η-3/η SP[S1,L]SP[S2,L]+κ1(M/m1)^2(-3/2SP[L,S1]^2+1/2SP[S1,S1])+κ2(M/m2)^2(-3/2SP[L,S2]^2+1/2SP[S2,S2]);
ENS[Energy2PN[S1,S2,m1,m2]-Energy2PNPoi]

E2p5PNcheck=ENS[(11-61/9η)(SP[S1,L]+SP[S2,L])+(m1-m2)(3-10/3η)(SP[S2,L]/m2-SP[S1,L]/m1)];
ENS[SeriesCoefficient[E2p5PNcheck,{SP[S1,L],0,1}]/.m2->MT-m1]

E3p5PNcheck=ENS[(135/4-367/4η+29/12η^2)(SP[S1,L]+SP[S2,L])+(m1-m2)(27/4-39η+5/4η^2)(SP[S2,L]/m2-SP[S1,L]/m1)];
ENS[SeriesCoefficient[E3p5PNcheck,{SP[S1,L],0,1}]/.m2->MT-m1]

3 PN

(*-1/2η M x factored out*)
Energy3PN[S1_,S2_,m1_,m2_]:= -γSS3PNm2+2γSS3PN-4E1PNa γSS2PN-E1PNb γSS2PN+2E1PNc γSS2PN+E1p5PN0*γSO1p5PN+Energy2PN0 3γ1PN+(-1/4(SS((-2δe κm+8κp-28)+η(2κp+4))+SΣ((10δe κp-52δe-10κm)+η(2δe κp+4δe+6κm))+ΣΣ((-5δe κm+5κp-24)+η(δe κm-11κp+76)+η^2(-2κp-4))+projSS(60-18κp+η(-6κp-12)+6δe κm-6κp+28)+projΣΣ(9δe κm-9κp+48+η(3δe κm+15κp-156)+η^2(6κp+12)+6δe κm-6κp+24+η(-6δe κm+18κp-76))+projSΣ(-18δe κp+108δe+18κm+η(-6δe κp-12δe+6κm)-12δe κp+52 δe+12κm+η(-24κm)))-1/2(SS(-3δe κm+5κp+8)+SΣ(8δe κp+8δe-8κm+η(12κm))+ΣΣ(4κp-4δe κm+η(3δe κm-11κp-10))+projSS(9δe κm-15κp-36)+projΣΣ(12δe κm-12κp+η(-9δe κm+33κp+30))+projSΣ(-24δe κp-32δe+24κm+η(-36κm))));
Energy3PNnoquad[S1_,S2_,m1_,m2_]:=((6(M/m1)^2-3/2(M/m1)-11/18)SP[L,S1]^2+(-(M/m1)^2-1/6M/m1-1/2)SP[S1,S1]+(6(M/m2)^2-3/2(M/m2)-11/18)SP[L,S2]^2+(-(M/m2)^2-1/6M/m2-1/2)SP[S2,S2]+(-11/3/η+23/18)SP[L,S1]SP[L,S2]+(2/η-11/6)SP[S1,S2]);
Energy3PNk1=ENS[SeriesCoefficient[Energy3PN[S1,S2,m1,m2],{κ1,0,1}]];
Energy3PNk2=ENS[SeriesCoefficient[Energy3PN[S1,S2,m1,m2],{κ2,0,1}]];

Check

Check E1 .5 PN

ENS[E1p5PN-((2/3+2M/m1)SP[L,S1]+(2/3+2M/m2)SP[L,S2])]

Energy3PNcheck=SS(-5δe κm/3-25κp/6+50/9+η(5κp/6+5/3))+SΣ(-5δe κp/2+25δe/3+5/2κm+η(5δe κp/6+5δe/3+35κm/6))+ΣΣ(5δe κm/4-5κp/4+5+η(5δe κm/4+5κp/4-10)+η^2(-5κp/6-5/3));
ENS[ENS[Energy3PN[S1,S2,m1,m2]-Energy3PNcheck]/.SP[L,S1]^2->SP[S1,S1]/.SP[L,S2]^2->SP[S2,S2]/.SP[L,S1]SP[L,S2]->SP[S1,S2]]

ENS[(Energy3PNcheck)-((5(M/m1)^2-5/3M/m1-10/9)SP[S1,S1]+(5(M/m2)^2-5/3M/m2-10/9)SP[S2,S2]+(-5/3/η-5/9)SP[S1,S2]+κ1(-5/2(M/m1)^2-5/2M/m1-5/6)SP[S1,S1]+κ2(-5/2(M/m2)^2-5/2M/m2-5/6)SP[S2,S2])]

ENS[(Energy3PN[S1,S2,m1,m2]-(κ1 Energy3PNk1+κ2 Energy3PNk2+Energy3PNnoquad[S1,S2,m1,m2]))]

ENS[Energy3PNk1-((-3(M/m1)^2-3M/m1-1)5/4SP[L,S1]^2+((M/m1)^2+M/m1+1/3)5/4SP[S1,S1])]

ENS[Energy3PNk2-((-3(M/m2)^2-3M/m2-1)5/4SP[L,S2]^2+((M/m2)^2+M/m2+1/3)5/4SP[S2,S2])]

L

Energy SO 3.5 PN -- Eq.(4.6) of arXiv:12125520

E3p5PNSO=ENS[(135/4-367/4η+29/12η^2)(SP[S1,LN]+SP[S2,LN])+(27/4-39η+5/4η^2)(m1-m2)(SP[S2,LN]/m2-SP[S1,LN]/m1)]

ENS[(SeriesCoefficient[E3p5PNSO,{SP[S1,LN],0,1}]-(-75/4+27/(4m1/MTot)+53 m1/MTot/2+67((m1/MTot)^2)/6+17(m1/MTot)^3/12-(m1/MTot)^4/12))/.MTot->m1+m2]

Flux S1^2

2 PN

Qtriple[i_,j_]:=(V3[b,i]V3[x,j]+3V3[a,i]V3[v,j]+3V3[v,i]V3[a,j]+V3[x,i]V3[b,j]);
ENS[Contract[Module[{i,j},Qtriple[i,j]Qtriple[i,j]]]/.Pair[Momentum[x,3],Momentum[v,3]]->0/.Pair[Momentum[x,3],Momentum[b,3]]->0/.Pair[Momentum[v,3],Momentum[a,3]]->0/.Pair[Momentum[x,3],Momentum[x,3]]->r^2/.Pair[Momentum[x,3],Momentum[a,3]]->γ/.Pair[Momentum[v,3],Momentum[v,3]]->ω^2r^2/.Pair[Momentum[a,3],Momentum[a,3]]->γ^2/r^2/.Pair[Momentum[v,3],Momentum[b,3]]->ω^2γ/.Pair[Momentum[b,3],Momentum[b,3]]->ω^2γ^2/r^2]

Eq.(4.12) of arXiv : 1501529 v1

Ffac=1/96;(*The leading order is assumed 1 as 32/5 η^2 is factored out*)
Flux2PN0=Ffac(SS(288κp+576)+SΣ(288δe κp+576δe-288κm)+ΣΣ((-144δe κm+144κp+6)+η(-288κp-576))+
projSS((-1008κp-2016)+(144κp+288))+
projΣΣ((504δe κm-504κp)+η(1008 κp+2016)+(-72δe κm+72κp+2+η(-144κp-288)))+
projSΣ((-1008δe κp-2016δe+1008κm)+(144δe κp+288δe-144κm)));
Flux2PN[S1_,S2_,m1_,m2_]:=Flux2PN0+2γSS2PN;
Flux2PNnoquad[S1_,S2_,m1_,m2_]:=Ffac((M/m1)^2(7 SP[S1,S1]-SP[S1,L]SP[S1,L])+(M/m2)^2(7SP[S2,S2]-1SP[S2,L]SP[S2,L])-(*494*)206/η SP[S1,S2]+(*1442*)578/η SP[L,S1]SP[L,S2])
Flux2PNk1=ENS[SeriesCoefficient[Flux2PN[S1, S2, m1, m2], {κ1, 0, 1}]];
Flux2PNk2=ENS[SeriesCoefficient[Flux2PN[S1, S2, m1, m2], {κ2, 0, 1}]];

Checks

ENS[Flux2PN[S1, S2, m1, m2] - (κ1 Flux2PNk1+κ2 Flux2PNk2+Flux2PNnoquad[S1,S2,m1,m2])]

Eq. 5.4 of 0812.4413

Flux2PNRacine[S1_,S2_,m1_,m2_]:=2/5/(32/5)(48S0S0-8*21projS0S0+24projS0S0+ΣΣ+1/3projΣΣ)+(2γSS2PN/.κ1->1/.κ2->1)
ENS[(Flux2PN[S1, S2, m1, m2]/.κ1->1/.κ2->1)-Flux2PNRacine[S1,S2,m1,m2]]

Eq. 4.14 of arXiv:1501.01529v1

F2PNA=((1/16(M/m1)^2)SP[S1,S1]+(1/16(M/m2)^2)SP[S2,S2]+(31/8/η)SP[S1,S2]+κ1(2(M/m1)^2)SP[S1,S1]+κ2(2(M/m2)^2)SP[S2,S2]);
ENS[Flux2PNcheck[S1,S2,m1,m2]-F2PNA]

Flux2PNcheck[S1_,S2_,m1_,m2_]:=(SS(2κp+4)+SΣ(2δe κp+4δe-2κm)+ΣΣ(-δe κm+κp+1/16+η(-2κp-4)))
ENS[(Flux2PN[S1,S2,m1,m2]-Flux2PNcheck[S1,S2,m1,m2])/.SP[L,S1]^2->SP[S1,S1]/.SP[L,S2]^2->SP[S2,S2]/.SP[L,S1]SP[L,S2]->SP[S1,S2]]

Eq.(33) of arXiv : gr - qc/9709032

Flux2PNk1Poi[m_,S_]:=((M/m)^2(3SP[L,S]^2-SP[S,S]))
ENS[Flux2PNk1-Flux2PNk1Poi[m1,S1]]
ENS[Flux2PNk2-Flux2PNk1Poi[m2,S2]]

3 PN

fac1PN=1/336;
Flux1PNa=fac1PN(785-852η)(*x^6*);
Flux1PNb=fac1PN(-160(17-η))(*x^3γ^3*);
Flux1PNc=fac1PN(16(1-4η))(*γ^6*);
Flux1PN=Flux1PNa+Flux1PNb+Flux1PNc+2γ1PN;
fac1p5PN=1/12;
Flux1p5PN0a=fac1p5PN((-43M/m1+6)SP[L,S1]+(-43M/m2+6)SP[L,S2])(*γ^2*);
Flux1p5PN0b=fac1p5PN((4M/m1-16)SP[L,S1]+(4M/m2-16)SP[L,S2])(*γ^5*);
Flux1p5PN=Flux1p5PN0a+Flux1p5PN0b+2γSO1p5PN;
Flux3PN0=1/(32*21)(SS(366δe κm+5034κp+18276+η(-4584κp-9168))+SΣ(4668δe κp+20812δe-4668κm+η(-4584δe κp-9168δe+3120κm))+ΣΣ(-2334δe κm+2334κp+3796+η(1926δe κm-6594κp-23336)+η^2(4584κp+9168))+projSS(570δe κm-14778κp-6546+η(13632κp+27264)-1668δe κm-324κp-6478+η(120κp+240))+projΣΣ(7674δe κm-7674κp+3090+η(-7386δe κm+22734κp+10884)+η^2(-13632κp-27264)-672δe κm+672κp-2242+η(1608δe κm-2952κp+9788)+η^2(-120κp-240))+projSΣ(-15348δe κp -8664δe+15348κm+η(13632δe κp+27264δe-15912κm)+1344δe κp-8188δe-1344κm+η(120δe κp+240δe+6552κm))+SS(2572δe κm-21028κp-41832+η(720κp+1440))+SΣ(-23600δe κp-48872δe+23600κm+η(720δe κp+1440δe-11008κm))+ΣΣ(11800δe κm-11800κp-4408+η(-2932δe κm+26532κp+56288)+η^2(-720κp-1440))+projSS(-12380δe κm+77736κp+182752+η(-1992κp-3984)+4664δe κm-14652κp-27240+η(-168κp-336))+projΣΣ(-45058δe κm+45058κp+9700+η(13376δe κm-103492κp-185532)+η^2(1992κp+3984)+9658δe κm-9658κp+4784+η(-4580δe κm+23896κp+9808)+η^2(168κp+336))+projSΣ(90116δe κp+183716δe-90116κm+η(-1992δe κp-3984δe+51512κm)-19316δe κp-18480δe+19316κm+η(-168δe κp-336δe-18488κm))+SS(16δe κm+144κp+368+η(-576κp-1152))+SΣ(128δe κp+224δe-128κm+η(-576δe κp-1152δe+512κm))+ΣΣ(-64δe κm+64κp+24+η(272δe κm-400κp-384)+η^2(576κp+1152))+projSS(-48δe κm-432κp-936+η(1728κp+3456))+projΣΣ(192δe κm-192κp-16+η(-816δe κm+1200κp+928)+η^2(-1728κp-3456))+projSΣ(-384δe κp-784δe+384κm+η(1728δe κp+3456δe-1536κm)));
Flux3PN[S1_,S2_,m1_,m2_]:=γSS3PN2+(Flux1PNb 3+6Flux1PNc) γSS2PN+Flux1p5PN0a*2γSO1p5PN+Flux1p5PN0b*5γSO1p5PN+Flux2PN0*4 γ1PN+Flux3PN0;
Flux3PNnoquad=(189/16(M/m1)^2-35/144(M/m1)+47/144)SP[S1,S1]+(189/16(M/m2)^2-35/144(M/m2)+47/144)SP[S2,S2]+(-5647/168/η-2023/72) SP[L,S1]SP[L,S2]+(2123/84/η+821/72)SP[S1,S2]+(-239/16 (M/m1)^2+293/144M/m1+299/144)SP[L,S1]^2+(-239/16 (M/m2)^2+293/144M/m2+299/144)SP[L,S2]^2;
Flux3PNk1=(-837/112(M/m1)^2-135/16(M/m1)+129/8)SP[L,S1]^2+(279/112(M/m1)^2+45/16(M/m1)-43/8)SP[S1,S1];
Flux3PNk2=(-837/112(M/m2)^2-135/16(M/m2)+129/8)SP[L,S2]^2+(279/112(M/m2)^2+45/16(M/m2)-43/8)SP[S2,S2];

Checks

ENS[Flux1PN-(-(1247/336+η 35/12))]

ENS[Flux1p5PN-((-5/4M/m1-3/2)SP[L,S1]+(-5/4M/m2-3/2)SP[L,S2])]

Flux2p5PNSO=ENS[(-9/2+272/9η)(SP[S1,L]+SP[S2,L])+(-13/16 +43/4η)(m1-m2)(SP[S2,L]/m2-SP[S1,L]/m1)];
ENS[SeriesCoefficient[Flux2p5PNSO,{SP[S1,L],0,1}]/.m2->MT-m1]

Flux3PNSO=ENS[-16 Pi(SP[S1,L]+SP[S2,L])+(-31/6 Pi)(m1-m2)(SP[S2,L]/m2-SP[S1,L]/m1)];
ENS[SeriesCoefficient[Flux3PNSO,{SP[S1,L],0,1}]/.m2->MT-m1]

Flux3p5PNSO=ENS[(476645/6804+6172/189η-2810/27η^2)(SP[S1,L]+SP[S2,L])+(9535/336+1849/126η-1501/36η^2)(m1-m2)(SP[S2,L]/m2-SP[S1,L]/m1)];
ENS[SeriesCoefficient[Flux3p5PNSO,{SP[S1,L],0,1}]/.m2->MT-m1]

Flux4PNSO=ENS[Pi((-3485/96+13789/72η)(SP[S1,L]+SP[S2,L])+(-7163/672+130583/2016η)(m1-m2)(SP[S2,L]/m2-SP[S1,L]/m1))];
ENS[SeriesCoefficient[Flux4PNSO,{SP[S1,L],0,1}]/.m2->MT-m1]

Flux3PNcheck=SS(41/16δe κm-271/112κp-5239/504+η(-43/4κp-43/2))+SΣ(-279/56δe κp-817/56δe+279/56κm+η(-43/4δe κp-43/2δe+κm/2))+ΣΣ(279/112δe κm-279/112κp-25/8+η(45/16δe κm+243/112κp+344/21)+η^2(43κp/4+43/2));
ENS[ENS[Flux3PN[S1,S2,m1,m2]-Flux3PNcheck]/.SP[S1,L]^2->SP[S1,S1]/.SP[S2,L]^2->SP[S2,S2]/.SP[S1,L]SP[S2,L]->SP[S1,S2]]

ENS[(Flux3PNcheck)-((-25/8(M/m1)^2+43/24M/m1+173/72)SP[S1,S1]+(-25/8(M/m2)^2+43/24M/m2+173/72)SP[S2,S2]+(-467/56/η-601/36)SP[S1,S2]+κ1(-279/56(M/m1)^2-45/8M/m1+43/4)SP[S1,S1]+κ2(-279/56(M/m2)^2-45/8M/m2+43/4)SP[S2,S2])]

ENS[SeriesCoefficient[Flux3PN[S1,S2,m1,m2],{κ1,0,1}]-Flux3PNk1]
ENS[SeriesCoefficient[Flux3PN[S1,S2,m1,m2],{κ2,0,1}]-Flux3PNk2]
ENS[Flux3PN[S1,S2,m1,m2]-(κ1 Flux3PNk1+κ2 Flux3PNk2+Flux3PNnoquad)]

wdotTT5

En1PN[η_]:=-(3/4+η/12);
En1p5PNSO[mM_]:=(2/3+2/mM)SpOr;
En2PN[η_]:=-(27/8-19/8η+1/24η^2);
En2p5PNSO[mM_]:=(5/3+3/mM+29mM/9+mM^2/9)SpOr;
En3PN[η_]:=-(675/64-(34445/576-205/96*Pi^2)*η+155/96*η^2+35/5184*η^3);
En3p5PNSO[mM_]:=(-75/4+27/(4mM)+53mM/2+67mM^2/6+17mM^3/12-mM^4/12)SpOr;
En5PNTidal[mM_]:=-9mM^4(1-mM)Td;
En6PNTidal[mM_]:=(-33/2+11/2mM-11/2mM^2+33/2mM^3)mM^4 Td;

Fl1PN[η_]:=-(1247/336+35/12η);
Fl1p5PN=4Pi;
Fl1p5PNSO[mM_]:=(-3/2-5/4/mM )SpOr;
Fl2PN[η_]:=-(44711/9072-9271/504η-65/18η^2);
Fl2p5PNSO[mM_]:=(63/8-13/(16mM)-73mM/36-157mM^2/18)SpOr;
Fl2p5PN[η_]:=-(8191/672+583/24η)Pi;
Fl3PN[η_]:=
(6643739519/69854400+16/3Pi^2-1712/105*EulerGamma-1712/105Log[4]+(41/48Pi^2-134543/7776)η-94403/3024η^2-775/324η^3);
Fl3PNSO[mM_]:=Pi(-17/3-31/(6mM))SpOr;
Fl3p5PN[η_]:=-(16285/504-214745/1728η-193385/3024η^2)Pi;
Fl3p5PNSO[mM_]:=(380647/13608+9535/(336mM)-40115mM/756+3742mM^2/63-35mM^3/108-1117mM^4/54)SpOr;
Fl5PNTidal[mM_]:=6(3-2mM)mM^4 Td;
Fl6PNTidal[mM_]:=(-176/7-1803/28mM+643/4mM^2-155/2mM^3)mM^4 Td;

wdotm1=ENS[(1+2v^2 En1PN[η]+5/2 v^3En1p5PNSO[mM]+3v^4 (En2PN[η]+En2PNSS SpOr^2)+7/2v^5En2p5PNSO[mM]+4v^6En3PN[η]+9/2v^7En3p5PNSO[mM]+6v^10En5PNTidal[mM]+7v^12En6PNTidal[mM])/(1+v^2 Fl1PN[η]+v^3(Fl1p5PN+Fl1p5PNSO[mM])+v^4(Fl2PN[η]+Fl2PNSS SpOr^2)+v^5(Fl2p5PN[η]+Fl2p5PNSO[mM])+v^6(Fl3PN[η]+Fl3PNSO[mM])+v^7(Fl3p5PN[η]+Fl3p5PNSO[mM])+v^10Fl5PNTidal[mM]+v^12Fl6PNTidal[mM])];

1 PN

wdot1PNm1=ENS[SeriesCoefficient[wdotm1,{v,0,2}]]
wdot1PNm1Check=ENS[2En1PN[η]-Fl1PN[η]];
ENS[wdot1PNm1-wdot1PNm1Check]

1.5 PN

wdot1p5PNm1=ENS[SeriesCoefficient[wdotm1,{v,0,3}]]

ENS[113/12((SpOr1+SpOr2)/2(1-76 η/113)+(m1-m2)/M(SpOr1-SpOr2)/2)/.SpOr2->0/.η->m1/M(1-m1/M)/.m2->M-m1]

2 PN

wdot2PNm1=ENS[SeriesCoefficient[wdotm1,{v,0,4}]]

2.5 PN

wdot2p5PNm1=ENS[SeriesCoefficient[wdotm1,{v,0,5}]];
wdot2p5PNm1/.SpOr->0
ENS[SeriesCoefficient[wdot2p5PNm1,{SpOr,0,1}]/.η->mM(1-mM)]

ENS[1/2(SpOr1- SpOr2)(m1/M-m2/M)(7η/2+146597/2016)+(SpOr1+SpOr2)/2(-17η^2/2-1213η/18+146597/2016)/.m2->M-m1/.η->m1/M(1-m1/M)/.SpOr2->0]

3 PN

wdot3PNm1=ENS[SeriesCoefficient[wdotm1,{v,0,6}]];
wdot3PNm1/.SpOr->0
ENS[SeriesCoefficient[wdot3PNm1,{SpOr,0,1}]/.η->mM(1-mM)]
ENS[SeriesCoefficient[wdot3PNm1,{SpOr,0,2}]]

ENS[Fl1p5PNSO[mM](-5/2 En1p5PNSO[mM]+Fl1p5PNSO[mM])]

3.5 PN

wdot3p5PNm1=ENS[SeriesCoefficient[wdotm1,{v,0,7}]];
wdot3p5PNm1/.SpOr->0
ENS[SeriesCoefficient[wdot3p5PNm1,{SpOr,0,1}]/.η->mM(1-mM)]
ENS[SeriesCoefficient[wdot3p5PNm1,{SpOr,0,2}]]
ENS[SeriesCoefficient[wdot3p5PNm1,{SpOr,0,3}]]
ENS[SeriesCoefficient[wdot3p5PNm1,{SpOr,0,4}]]
ENS[-Fl1p5PN(3En2PNSS-2Fl2PNSS)]
ENS[-Fl1p5PNSO[mM](3 En2PNSS-2Fl2PNSS)-5/2En1p5PNSO[mM] Fl2PNSS]

linfit=a*x+b;
x={0.1,0.2,0.3,0.4,0.5,0.6,0.7};
yw={-44.53594134653,-32.52717259115,-20.51840383577,-8.50963508039,3.49913367499,15.50790243037,15.50790243037};
yef={-44.75649466685,-33.39624524743,-22.44425090614,-11.88080070548,-1.68618370797,8.15931102391,17.67539442765};

(yw[[1]]-yw[[2]])/0.1
(yw[[2]]-yw[[3]])/0.1
(yw[[3]]-yw[[4]])/0.1
(yef[[1]]-yef[[2]])/0.1
(yef[[2]]-yef[[3]])/0.1
(yef[[3]]-yef[[4]])/0.1

Tidal

wdotTd5PNm1=ENS[SeriesCoefficient[SeriesCoefficient[wdotm1,{v,0,10}],{Td,0,1}]]
wdotTd6PNm1=ENS[SeriesCoefficient[SeriesCoefficient[wdotm1,{v,0,12}],{Td,0,1}]/.η->mM(1-mM)]

wdotTT4

2 PN

wdot2PN=ENS[SeriesCoefficient[(1+x*x Flux2PN[S1,S2,m1,m2])/(1+3x*x Energy2PN[S1,S2,m1,m2]),{x,0,2}]];
wdot2PNcheck=ENS[Flux2PN[S1,S2,m1,m2]-3Energy2PN[S1,S2,m1,m2]];
wdot2PNnoquad=1/96((M/m1)^2(-SP[L,S1]^2+7SP[S1,S1])+(M/m2)^2(-SP[L,S2]^2+7SP[S2,S2])+721*2/η SP[L,S1]SP[L,S2]-247*2/η SP[S1,S2]);
wdot2PNk1=ENS[SeriesCoefficient[wdot2PN,{κ1,0,1}]];
wdot2PNk2=ENS[SeriesCoefficient[wdot2PN,{κ2,0,1}]];
ENS[(wdot2PN/.κ1->0/.κ2->0)-wdot2PNnoquad]
ENS[wdot2PNk1-5/2(M/m1)^2(3SP[L,S1]^2-SP[S1,S1])]
ENS[wdot2PNk2-5/2(M/m2)^2(3SP[L,S2]^2-SP[S2,S2])]

Checks

ENS[wdot2PN-wdot2PNcheck]

Eqs. 7 - 9 of arXiv : gr - qc/0504538 v2

wdot2PNHun=(721/48/η SP[L,S1]SP[L,S2]-247/48/η SP[S1,S2]-1/96(M/m1)^2SP[L,S1]^2-1/96(M/m2)^2SP[L,S2]^2+7/96(M/m1)^2SP[S1,S1]+7/96(M/m2)^2SP[S2,S2]+κ1(M/m1)^2(15/2SP[L,S1]^2-5/2SP[S1,S1])+κ2(M/m2)^2(15/2SP[L,S2]^2-5/2SP[S2,S2]));
ENS[(wdot2PNcheck-wdot2PNHun)]

Eqs 5.15 - 5.17 of arXiv : 0812 : 4413 v2

wdot2PNRacine=(721/48/η SP[L,S1]SP[L,S2]-247/48/η SP[S1,S2]+1/96(M/m1)^2(719SP[L,S1]^2-233SP[S1,S1])+1/96(M/m2)^2(719SP[L,S2]^2-233SP[S2,S2]));
ENS[(wdot2PNHun/.κ1->1/.κ2->1)-wdot2PNRacine]

3 PN

wdot3PN=ENS[SeriesCoefficient[(1+x Flux1PN+x^(3/2) Flux1p5PN+x*x Flux2PN[S1,S2,m1,m2]+x*x*x Flux3PN[S1,S2,m1,m2])/(1+2x E1PN+5/2x^(3/2)E1p5PN+3x*x Energy2PN[S1,S2,m1,m2]+4x*x*x Energy3PN[S1,S2,m1,m2]),{x,0,3}]];
wdot3PNcheck=ENS[Flux3PN[S1,S2,m1,m2]-4Energy3PN[S1,S2,m1,m2]+(-3Flux1PN Energy2PN[S1,S2,m1,m2]-2Flux2PN[S1,S2,m1,m2]E1PN+12E1PN Energy2PN[S1,S2,m1,m2]+25/4E1p5PN^2-5/2E1p5PN Flux1p5PN)];
ENS[wdot3PN-wdot3PNcheck]

wdot3PNnoquad=(-493/64(M/m1)^2+19747/576M/m1+5645/576)SP[L,S1]^2+(-493/64(M/m2)^2+19747/576M/m2+5645/576)SP[L,S2]^2+(10879/672/η+7525/288)SP[S1,S2]+(1019/64(M/m1)^2+251/576M/m1+1333/576)SP[S1,S1]+(1019/64(M/m2)^2+251/576M/m2+1333/576)SP[S2,S2]+(16225/224/η-12931/288)SP[L,S1]SP[L,S2];
ENS[(wdot3PNcheck/.κ1->0/.κ2->0)-wdot3PNnoquad]
wdot3PNk1=ENS[SeriesCoefficient[wdot3PN,{κ1,0,1}]];
wdot3PNk2=ENS[SeriesCoefficient[wdot3PN,{κ2,0,1}]];
ENS[wdot3PNk1-((1977/224(M/m1)^2-73/16M/m1+129/4)SP[L,S1]^2+(-659/224 (M/m1)^2+73/48M/m1-43/4)SP[S1,S1])]
ENS[wdot3PNk2-((1977/224(M/m2)^2-73/16M/m2+129/4)SP[L,S2]^2+(-659/224 (M/m2)^2+73/48M/m2-43/4)SP[S2,S2])]

wdot4PNSO=ENS[SeriesCoefficient[SeriesCoefficient[(1+x F1PN+x^(3/2) (F1p5PN+λ F1p5PNSO)+x^2F2PN+x^(5/2)(F2p5PN+λ F2p5PNSO)+x^3λ F3PNSO+x^4λ F4PNSO)/(1+2x E1PN+5/2x^(3/2)λ E1p5PNSO+3x*x E2PN+7/2x^(5/2)λ E2p5PNSO+4x^3 (E3PN+λ E3PNSO)+5x^4λ E4PNSO),{x,0,4}],{λ,0,1}]]

wdotTT2

2 PN

dtdv2PN=ENS[SeriesCoefficient[(1+3x*x Energy2PN[S1,S2,m1,m2])/(1+x*x Flux2PN[S1,S2,m1,m2]),{x,0,2}]];
dtdv2PNcheck=ENS[3Energy2PN[S1,S2,m1,m2]-Flux2PN[S1,S2,m1,m2]];
dtdv2PNnoquad=(1/96((M/m1)^2(SP[L,S1]^2-7SP[S1,S1])+(M/m2)^2(SP[L,S2]^2-7SP[S2,S2])-721*2/η SP[L,S1]SP[L,S2]+247*2/η SP[S1,S2]));
dtdv2PNk1=ENS[SeriesCoefficient[dtdv2PN,{κ1,0,1}]];
dtdv2PNk2=ENS[SeriesCoefficient[dtdv2PN,{κ2,0,1}]];

Check

ENS[dtdv2PN-dtdv2PNcheck]
ENS[(dtdv2PN/.κ1->0/.κ2->0)-dtdv2PNnoquad]
ENS[dtdv2PNk1-5/2(M/m1)^2(-3SP[L,S1]^2+SP[S1,S1])]
ENS[dtdv2PNk2-5/2(M/m2)^2(-3SP[L,S2]^2+SP[S2,S2])]

dtdv2PNA=ENS[SeriesCoefficient[(1+3x^2 E2PNA)/(1+x^2 F2PNA),{x,0,2}]];
dtdv2PNAcheck=ENS[3E2PNA-F2PNA];
dtdv2PNAnoquad=(-1/16(M/m1)^2SP[S1,S1]-1/16(M/m2)^2SP[S2,S2]-79/8/η SP[S1,S2]);
dtdv2PNk1A=ENS[SeriesCoefficient[dtdv2PNA,{κ1,0,1}]];
dtdv2PNk2A=ENS[SeriesCoefficient[dtdv2PNA,{κ2,0,1}]];
ENS[dtdv2PNA-dtdv2PNAcheck]
ENS[(dtdv2PNA/.κ1->0/.κ2->0)-dtdv2PNAnoquad]
ENS[dtdv2PNk1A-(-5(M/m1)^2SP[S1,S1])]
ENS[dtdv2PNk2A-(-5(M/m2)^2SP[S2,S2])]

3 PN

dtdv3PN=ENS[SeriesCoefficient[(1+2x E1PN+5/2x^(3/2)E1p5PN+3x*x Energy2PN[S1,S2,m1,m2]+4x *x*x Energy3PN[S1,S2,m1,m2])/(1+x Flux1PN+x^(3/2) Flux1p5PN+x*x Flux2PN[S1,S2,m1,m2]+x*x*x Flux3PN[S1,S2,m1,m2]),{x,0,3}]];
dtdv3PNcheck=ENS[4Energy3PN[S1,S2,m1,m2]-Flux3PN[S1,S2,m1,m2]+(-3Flux1PN Energy2PN[S1,S2,m1,m2]-2Flux2PN[S1,S2,m1,m2]E1PN+2Flux1PN Flux2PN[S1,S2,m1,m2]+Flux1p5PN^2-5/2E1p5PN Flux1p5PN)+(2E1PN*Flux1PN^2-Flux1PN^3)];
ENS[dtdv3PN-dtdv3PNcheck]

dtdv3PNnoquad=((754979/16128(M/m1)^2+1543/288M/m1+49/288)SP[L,S1]^2+(754979/16128(M/m2)^2+1543/288M/m2+49/288)SP[L,S2]^2+(52973/8064/η+313/144)SP[S1,S2]+(-37427/2304(M/m1)^2-241/288M/m1-551/288)SP[S1,S1]+(-37427/2304(M/m2)^2-241/288M/m2-551/288)SP[S2,S2]+(-170603/8064/η-2543/144)SP[L,S1]SP[L,S2]);
ENS[(dtdv3PN/.κ1->0/.κ2->0)-dtdv3PNnoquad-(2E1PN*Flux1PN^2-Flux1PN^3)]
dtdv3PNk1=ENS[SeriesCoefficient[dtdv3PN,{κ1,0,1}]];
dtdv3PNk2=ENS[SeriesCoefficient[dtdv3PN,{κ2,0,1}]];
ENS[dtdv3PNk1-((-9407/224(M/m1)^2-587/16M/m1+9)SP[L,S1]^2+(9407/672(M/m1)^2+587/48M/m1-3)SP[S1,S1])]
ENS[dtdv3PNk2-((-9407/224(M/m2)^2-587/16M/m2+9)SP[L,S2]^2+(9407/672(M/m2)^2+587/48M/m2-3)SP[S2,S2])]