Octupole と GR の共鳴

Perturber が伴星である場合と惑星である場合、e_max vs a1 を調べて、GR と octupole terms もつけたり消したりして、期待通りの結果が出た（惑星２つのケースでは、GR なしでも古在がきちんとかかっていないことがあるので、計算し直しの最中）。

GR は基本的に dg1/dt を足しているだけなのに、これだけの違いが出る。de1/dt の中の C3 の項との関係性を見てみる。

Q1: 何故 octupole terms としか共鳴しないのか（何故 quadrupole terms とは共鳴しないのか）?
Q2: 何故 inclination が低い時に、共鳴が顕著に見られるのか？
Q3: Kozai と perturber mass (m2) のリミット
　　すっかり見落としていたが、m2 があまりに小さいときちんと e_max が上がっていかないな。
Q3: Kozai と perturber mass (m2) のリミット
　　a1 は固定、a2 を動かして、e_max の位置を調べてみる。簡単な論文にはすぐなりそう。GJ436 b の文献を漁ってみる必要がある。