 into the values around which the angles of a resonant system should librate. The method described in the following is a generalization of the method presented in @delisle_dissipation_2012 which focuses on finding the fixed points for two-planet MMR.

I denote by $\lambda_i$ and $\varpi_i$ the mean longitude and longitude of periastron of planet $i$ (in astrocentric coordinates), respectively. The actions canonically conjugated to the angles $\lambda_i$ and $-\varpi_i$ are the circular angular momentum $\Lambda_i$ and the angular momentum deficit [AMD, see @laskar_spacing_2000] $D_i$, respectively. These actions are defined as follows $$\begin{aligned}
  \Lambda_i &=& \beta_i\sqrt{\mu_i a_i},\\
  D_i &=& \Lambda_i-G_i = \Lambda_i \left(1-\sqrt{1-e_i^2}\right),\end{aligned}$$ where $G_i=\Lambda_i\sqrt{1-e_i^2}$ is the angular momentum of planet $i$, $\beta_i = m_i m_0/(m_0+m_i)$, $\mu_i = \G (m_0+m_i)$. At low eccentricities the deficit of angular momentum $D_i$ is proportional to $e_i^2$. The Hamiltonian (Eq. (\[eq:hamposvel\])) can be expressed using these action-angle coordinates $$\label{eq:hamLaD}
  \H = - \sum_{i=1}^n \frac{\mu_i^2\beta_i^3}{2\Lambda_i^2}
  + \sum_{1\leq i<j\leq n} \H_{i,j}(\Lambda_i,\Lambda_j,D_i,D_j,\lambda_i,\lambda_j,\varpi_i,\varpi_j),$$ where the first sum is the Keplerian part, which depends only on $\Lambda_i$ (or equivalently $a_i$), and $\H_{i,j}$ is the perturbation between the planets $i$ and $j$, which depends on the eight coordinates associated to $i$ and $j$. I follow the method described in @laskar_stability_1995 to compute $H
into the values around which the angles of a resonant system should weigh. The method acting identify in the following is a abstraction of the method deliver in @delisle_dissipation_2012 which focuses on finding the pay back points for two - planet MMR. 

 I announce by $ \lambda_i$ and $ \varpi_i$ the base longitude and longitude of periastron of planet $ i$ (in astrocentric coordinate), respectively. The actions canonically conjugated to the angles $ \lambda_i$ and $ -\varpi_i$ are the circular angular momentum $ \Lambda_i$ and the angular momentum deficit [ AMD, determine @laskar_spacing_2000 ] $ D_i$, respectively. These actions are defined as follows $ $ \begin{aligned } 
   \Lambda_i & = & \beta_i\sqrt{\mu_i a_i},\\ 
   D_i & = & \Lambda_i - G_i = \Lambda_i \left(1-\sqrt{1 - e_i^2}\right),\end{aligned}$$ where $ G_i=\Lambda_i\sqrt{1 - e_i^2}$ is the angular momentum of planet $ i$, $ \beta_i = m_i m_0/(m_0+m_i)$, $ \mu_i = \G (m_0+m_i)$. At low eccentricity the deficit of angular momentum $ D_i$ is proportional to $ e_i^2$. The Hamiltonian (Eq.   (\[eq: hamposvel\ ]) ) can be expressed using these action - slant coordinates $ $ \label{eq: hamLaD } 
   \H = - \sum_{i=1}^n \frac{\mu_i^2\beta_i^3}{2\Lambda_i^2 } 
   + \sum_{1\leq i < j\leq n } \H_{i, j}(\Lambda_i,\Lambda_j, D_i, D_j,\lambda_i,\lambda_j,\varpi_i,\varpi_j),$$ where the first kernel is the Keplerian part, which depend only on $ \Lambda_i$ (or equivalently $ a_i$), and $ \H_{i, j}$ is the perturbation between the planets $ i$ and $ j$, which depends on the eight coordinates consociate to $ i$ and $ j$. I follow the method described in @laskar_stability_1995 to compute $ H
 inho the values around whigh the angles of a resoiant syatem shojld librate. The method descrmbed in tye following is a genefalizatioj of the metiod presented in @delisle_dissipzbion_2012 chmch focuses on ninding the fixed points xof cwo-planet MMR.

I denote by $\lambda_i$ and $\varpi_i$ tje mean longitode amq lohgitude of periastron of planet $i$ (jn astrmcentric coorcinates), respectively. The aftiojs canonically conuugated to jge wbgles $\lambda_k$ and $-\varpp_n$ are the cjrcular angular momentum $\Lambda_k$ and the angulqr molgntum deficiv [AMD, fee @laskar_spacing_2000] $D_i$, respecyively. These agtionv aee defined as follows $$\begin{aligned}
  \Lamfda_i &=& \betd_i\aqrt{\mu_i a_i},\\
  D_i &=& \Lqmvda_i-G_h = \Ldmbdx_u \ldft(1-\aqct{1-e_j^2}\right),\fnd{eligned}$$ whede $G_i=\Lambda_u\sqrt{1-e_i^2}$ is the angukaw momentum of llanet $i$, $\beta_i = m_i m_0/(m_0+m_i)$, $\mu_i = \G (m_0+m_i)$. At low eccenuricifies the deficit of angylar momentum $D_i$ is ptoportionaj to $e_i^2$. The Hamiltonian (Eq. (\[eq:hamposvel\])) can be exprassed jsiun gyede action-angle coordinates $$\label{eq:hamLaD}
  \H = - \auk_{i=1}^k \frac{\mu_i^2\beta_i^3}{2\Laibda_i^2}
  + \sum_{1\kee o<t\leq n} \H_{i,j}(\Lamcda_i,\Lambdz_j,D_i,D_j,\lambda_i,\lambdw_j,\varpi_y,\varpu_j),$$ where uhe forst sum is the Keplerian pqrt, which deiendw only on $\Lambda_i$ (lr equivaleutly $a_o$), and $\H_{i,j}$ is the perturbatiou betwsen the plajets $i$ ans $j$, which depends on tve eight coordinates assocyated to $m$ and $j$. I foulow the mqthod descgibed in @laskar_stability_1995 to clm[ute $H
into the values around which the angles resonant should librate. method described in of method presented in which focuses on the fixed points for two-planet MMR. denote by $\lambda_i$ and $\varpi_i$ the mean longitude and longitude of periastron of $i$ (in astrocentric coordinates), respectively. The actions canonically conjugated to the angles $\lambda_i$ $-\varpi_i$ the angular $\Lambda_i$ and the angular momentum deficit [AMD, see @laskar_spacing_2000] $D_i$, respectively. These actions are defined as $$\begin{aligned} \Lambda_i &=& \beta_i\sqrt{\mu_i a_i},\\ D_i &=& \Lambda_i-G_i \Lambda_i \left(1-\sqrt{1-e_i^2}\right),\end{aligned}$$ where $G_i=\Lambda_i\sqrt{1-e_i^2}$ the angular momentum of planet $\beta_i m_i m_0/(m_0+m_i)$, = (m_0+m_i)$. low eccentricities the of angular momentum $D_i$ is proportional to $e_i^2$. The Hamiltonian (Eq. (\[eq:hamposvel\])) can be expressed using these coordinates $$\label{eq:hamLaD} - \sum_{i=1}^n + i<j\leq \H_{i,j}(\Lambda_i,\Lambda_j,D_i,D_j,\lambda_i,\lambda_j,\varpi_i,\varpi_j),$$ where the is the Keplerian part, which depends (or equivalently $a_i$), and $\H_{i,j}$ is the perturbation the planets and $j$, which depends on the coordinates associated to $i$ and $j$. I follow method described in @laskar_stability_1995 to compute $H
 into the values around which tHe angles of A resoNanT syStEm shOuld Librate. The methOD desCribed in the following is A geneRaLIzatIOn Of the Method pREsENTed In @DeLisLe_DIsSipatIon_2012 Which foCuses on finDinG tHe fixed pointS FoR two-planet mMR.

i denote by $\lamBda_I$ and $\vaRpI_i$ tHE mean LonGitudE and loNGitude Of periastRoN Of planET $i$ (in astROCeNtriC coordinates), respeCTiVEly. The actions cAnonicAlLY cONJugAteD to the anglEs $\LambdA_I$ and $-\varPI_i$ ARE The CIrcular angulaR momentum $\LaMBda_I$ and thE aNguLAr momeNtum dEfICit [aMD, see @laskaR_spaCing_2000] $D_i$, resPectivELy. These ACtions aRe defiNed As fOlloWS $$\bEgIn{aLiGNed}
  \lAmBda_I &=& \BetA_i\sqrt{\mu_I a_I},\\
  D_I &=& \LambDa_i-G_I = \lAMBda_i \LefT(1-\sqrT{1-e_i^2}\riGht),\end{aligned}$$ WheRe $G_i=\lAmbDa_i\sqRt{1-e_i^2}$ iS the AnGular MomentUm of pLaNet $i$, $\beta_i = m_i m_0/(m_0+m_i)$, $\Mu_i = \G (M_0+m_i)$. At low eCceNtRicItIes thE DeficiT of AngUlar momEntum $D_i$ IS prOpORTIoNal to $e_i^2$. The HamiltonIaN (eQ. (\[eQ:hamposvEl\])) can bE ExPrESsed usinG tHesE actION-anglE cooRDiNates $$\labEl{eq:haMlad}
  \H = - \Sum_{i=1}^n \frAc{\Mu_i^2\betA_i^3}{2\lamBda_I^2}
  + \sum_{1\lEQ i<j\lEq n} \H_{i,j}(\lambda_i,\LAmbda_J,d_i,D_j,\lambda_i,\lamBDa_j,\varpi_i,\varpI_J),$$ wHERe THe fiRst Sum is the KepLeriAN parT, whiCH dEpeNDs onlY on $\LaMbDA_i$ (OR equivalently $a_i$), and $\H_{I,j}$ Is the pErturBation between The planets $I$ AND $j$, which dEpenDS oN The eight coordiNates Associated TO $i$ and $j$. I fOllow The methoD describeD IN @laskar_sTabIliTy_1995 tO coMPUtE $H
 into the values around wh ich the an glesofa r es onan t sy stem should li b rate . The method described  in t he foll o wi ng is  a gene r al i z ati on  o f t he me thodpre sentedin @delisl e_d is sipation_201 2  w hich focus eson finding t hefixedpo int s  fortwo -plan et MMR . 

I de note by $ \l a mbda_i $  and $\ v a rp i_i$  the mean longitu d ea nd longitude o f peri as t ro n ofpla net $i$ (i nastro c entricc oo r d i nat e s), respectiv ely. The ac t ion s cano ni cal l y conj ugate dt o t he angles $ \lam bda_i$ an d $-\v a rpi_i$a re thecircul arang ular mo me ntu m$ \La m bd a_i $  an d the an gu la r mom entu m d e fici t [ AMD,  see@laskar_spaci ng_ 2000 ]  $D _i$,respe ctiv el y. Th ese ac tions  a re defined as f ollo ws $$\beg in{ al ign ed }
  \ L ambda_ i & =&\beta_i \sqrt{\ m u_i  a _ i } ,\ \
  D_i &=& \Lambd a_ i - G_ i = \Lam bda_i\ le ft ( 1-\sqrt{ 1- e_i ^2}\ r i ght), \end { al igned}$$  where $G _i =\Lambd a_ i\sqrt {1 -e_ i^2 }$ is theangula r moment um of planet $i$, $\ b eta_i = m_i m _ 0/ ( m _0 + m_i) $,$\mu_i = \G  (m_ 0 +m_i )$.A tlow eccen trici ti e st he deficit of angul ar  momen tum $ D_i$ is propo rtional to $ e _i^2$. T he H a mi l tonian (Eq. (\ [eq:h amposvel\] ) ) can be  expr essed us ing these a ction-an gle  co ord ina t e s$$\label{eq:h a m LaD} 
 \H = -  \s um_{i=1 }^n  \f rac {\m u_ i^2\beta_ i^3}{2\L am bd a_ i^ 2}  + \ s um_{1\le qi<j \l eqn} \H _ {i,j}( \Lamb da_i ,\ La m bda _j,D_i, D _j , \ lamb da _i ,\la mbd a_ j,\va rpi_ i ,\v arpi_j) ,$$ where  th e  fir st  s um is t he Keplerianpa rt, whichde pen ds onl y on $\Lam bda_i$ (or equivalently $a_i$),  an d $\H _{i, j}$ is th e p erturb ati o n betw een th e pla ne ts$ i $ and $ j$ , w hi ch depends o n t he ei gh t co ordinat es associated to $ i $ a nd $j$. I fol low  the m et hod de s cri be d  in @ laskar_stabilit y_1995 toco m pu te $H
 into_the values_around which the angles_of a_resonant_system should_librate._The method described_in the following_is a generalization of_the method presented_in_@delisle_dissipation_2012 which focuses on finding the fixed points for two-planet MMR.

I denote by $\lambda_i$_and_$\varpi_i$ the_mean_longitude_and longitude of periastron of_planet $i$ (in astrocentric coordinates),_respectively. The_actions canonically conjugated to the angles $\lambda_i$ and_$-\varpi_i$_are the circular_angular momentum $\Lambda_i$ and the angular momentum deficit [AMD,_see @laskar_spacing_2000] $D_i$, respectively. These actions_are defined as_follows_$$\begin{aligned}
_ \Lambda_i &=& \beta_i\sqrt{\mu_i_a_i},\\
  D_i &=& \Lambda_i-G_i =_\Lambda_i \left(1-\sqrt{1-e_i^2}\right),\end{aligned}$$ where $G_i=\Lambda_i\sqrt{1-e_i^2}$ is the_angular momentum of planet $i$, $\beta_i =_m_i m_0/(m_0+m_i)$, $\mu_i = \G (m_0+m_i)$._At low eccentricities the deficit_of angular_momentum $D_i$ is proportional to_$e_i^2$. The Hamiltonian_(Eq. (\[eq:hamposvel\])) can_be expressed using_these action-angle coordinates $$\label{eq:hamLaD}
  \H_= - \sum_{i=1}^n_\frac{\mu_i^2\beta_i^3}{2\Lambda_i^2}
  + \sum_{1\leq i<j\leq n}_\H_{i,j}(\Lambda_i,\Lambda_j,D_i,D_j,\lambda_i,\lambda_j,\varpi_i,\varpi_j),$$_where the first_sum_is_the Keplerian_part, which depends_only_on $\Lambda_i$_(or_equivalently $a_i$), and $\H_{i,j}$ is the_perturbation_between the planets $i$ and $j$, which_depends on the eight_coordinates_associated to $i$ and_$j$. I follow the method_described in @laskar_stability_1995 to compute $H