_{i,j}$ as a power series of the eccentricities (or equivalently of $\sqrt{D_i}$ and $\sqrt{D_j}$), and a Fourrier series of the angles, where the coefficients are functions of $\Lambda_i$ and $\Lambda_j$ (i.e., of the semi-major axes). For a system that is close to the resonance or resonant, the semi-major axes remain close to the nominal resonant values (Kepler’s third law) $$\label{eq:nominalsma}
  \frac{a_i}{a_j} \approx \frac{a_{i,0}}{a_{j,0}} = \left(\frac{k_{i,j}}{k_{j,i}}\right)^{2/3} \left(\frac{\mu_i}{\mu_j}\right)^{1/3}.$$ I introduce $$\Delta\Lambda_i = \Lambda_i - \Lambda_{i,0},$$ where $$\Lambda_{i,0} = \beta_i\sqrt{\mu_i a_{i,0}},$$ and expand the Keplerian part at degree 2, and the perturbative part at degree 0 in $\Delta\Lambda_i$ $$\begin{aligned}
\label{eq:hamdLaD}
  \H &=& \sum_{i=1}^n n_{i,0} \Delta\Lambda_i
  - \frac{3}{2}\frac{n_{i,0}}{\Lambda_{i,0}}\Delta\Lambda_i^2\nonumber\\
  &+& \sum_{1\leq i<j\leq n} \H_{i,j}(D_i,D_j,\lambda_i,\lambda_j,\varpi_i,\varpi_j),\end{aligned}$$ where $n_{i,0}$ is the nominal mean motion of planet $i$, such that $$\frac{n_{i,0}}{n_{j,0}} = \frac{k_{j,i}}{k_{i,j}}.$$ The perturbative part does not depend on $\Lambda_i$ anymore, but is simply evaluated at $\Lambda_{i,0}$.

In order to perform the reductions associated to the conservation of angular momentum and to the averaging, I first change the system of coordinates. For the sake of readability, I present the general case (with any number of planets, in any resonance of any
_ { i, j}$ as a power series of the eccentricities (or equivalently of $ \sqrt{D_i}$ and $ \sqrt{D_j}$), and a Fourrier series of the angle, where the coefficient are functions of $ \Lambda_i$ and $ \Lambda_j$ (i.e., of the semi - major axes). For a arrangement that is airless to the resonance or resonant, the semi - major axes stay close to the nominal evocative value (Kepler ’s third law) $ $ \label{eq: nominalsma } 
   \frac{a_i}{a_j } \approx \frac{a_{i,0}}{a_{j,0 } } = \left(\frac{k_{i, j}}{k_{j, i}}\right)^{2/3 } \left(\frac{\mu_i}{\mu_j}\right)^{1/3}.$$ I precede $ $ \Delta\Lambda_i = \Lambda_i - \Lambda_{i,0},$$ where $ $ \Lambda_{i,0 } = \beta_i\sqrt{\mu_i a_{i,0}},$$ and expand the Keplerian part at academic degree 2, and the perturbative part at degree 0 in $ \Delta\Lambda_i$ $ $ \begin{aligned } 
 \label{eq: hamdLaD } 
   \H & = & \sum_{i=1}^n n_{i,0 } \Delta\Lambda_i 
   - \frac{3}{2}\frac{n_{i,0}}{\Lambda_{i,0}}\Delta\Lambda_i^2\nonumber\\ 
   & + & \sum_{1\leq i < j\leq n } \H_{i, j}(D_i, D_j,\lambda_i,\lambda_j,\varpi_i,\varpi_j),\end{aligned}$$ where $ n_{i,0}$ is the nominal mean motion of planet $ i$, such that $ $ \frac{n_{i,0}}{n_{j,0 } } = \frac{k_{j, i}}{k_{i, j}}.$$ The perturbative part does not count on $ \Lambda_i$ anymore, but is simply evaluated at $ \Lambda_{i,0}$. 

 In order to do the reductions associated to the conservation of angular momentum and to the averaging, I first change the system of coordinate. For the sake of readability, I give the general case (with any number of planet, in any resonance of any
_{i,j}$ ws a power series of the eccentricities (or equmvalentmy of $\sqft{D_i}$ and $\sqrt{D_j}$), and a Fourrixr swries of the angles, where tfe coeffibients arw fuictions of $\Lambde_j$ and $\Lambda_j$ (l.e., of vhe semi-major awes). For a sfstem that is wlusz to the resonance or resonant, the sqmi-majot wxes remain cljse uo ehe hominal resonant values (Kepler’s thjrd law) $$\label{eq:nominslsma}
  \frac{a_i}{a_j} \approx \fraf{a_{i,0}}{a_{u,0}} = \left(\frac{k_{i,j}}{k_{j,i}}\rlght)^{2/3} \left(\frqc{\mu_y}{\nu_j}\right)^{1/3}.$$ I ivtroduce $$\Dtlca\Lambda_i = \Mambda_i - \Lambda_{i,0},$$ where $$\Lambda_{i,0} = \beta_n\sqrt{\mu_i a_{i,0}},$$ abd ffpand the Kxpleriwn part at dcbree 2, dnd the perturbative iart et dwgree 0 in $\Delta\Lambda_m$ $$\begin{aligned}
\label{ez:hamdLaD}
  \H &=& \sum_{i=1}^n n_{i,0} \Delta\Lqmvda_i
  - \frdc{3}{2}\frxx{n_{i,0}}{\Uamgde_{i,0}}\Dslta\Lalbde_i^2\nonumber\\
  &+& \sum_{1\leq i<j\oeq n} \H_{i,j}(D_i,D_j,\lambda_o,\lwnbda_j,\varpi_i,\vadpi_j),\enq{ajigned}$$ where $n_{i,0}$ is the nominal mean motpon kf planet $i$, such that $$\feac{n_{i,0}}{n_{j,0}} = \frac{k_{j,i}}{k_{i,j}}.$$ Tje perturfative part does not depend on $\Lambda_i$ anymore, bud is akmpoy dcapuated at $\Lambda_{i,0}$.

In order to perform the redusfipnx associated tj the consetvwtojn of angular momenchm and to the averaglng, I fyrst xhange tht sysyem of coordinates. For the wake of readcbioity, I present the general care (woth amy number of planets, in any resonance lf any
_{i,j}$ as a power series of the equivalently $\sqrt{D_i}$ and and a Fourrier the are functions of and $\Lambda_j$ (i.e., the semi-major axes). For a system is close to the resonance or resonant, the semi-major axes remain close to nominal resonant values (Kepler’s third law) $$\label{eq:nominalsma} \frac{a_i}{a_j} \approx \frac{a_{i,0}}{a_{j,0}} = \left(\frac{k_{i,j}}{k_{j,i}}\right)^{2/3} \left(\frac{\mu_i}{\mu_j}\right)^{1/3}.$$ introduce = - where $$\Lambda_{i,0} = \beta_i\sqrt{\mu_i a_{i,0}},$$ and expand the Keplerian part at degree 2, and the perturbative at degree 0 in $\Delta\Lambda_i$ $$\begin{aligned} \label{eq:hamdLaD} \H \sum_{i=1}^n n_{i,0} \Delta\Lambda_i - &+& \sum_{1\leq i<j\leq n} \H_{i,j}(D_i,D_j,\lambda_i,\lambda_j,\varpi_i,\varpi_j),\end{aligned}$$ $n_{i,0}$ the nominal motion planet such that $$\frac{n_{i,0}}{n_{j,0}} \frac{k_{j,i}}{k_{i,j}}.$$ The perturbative part does not depend on $\Lambda_i$ anymore, but is simply evaluated at $\Lambda_{i,0}$. In to perform associated to conservation angular and to the first change the system of coordinates. of readability, I present the general case (with number of in any resonance of any
_{i,j}$ as a power series of the ecceNtricities (Or equIvaLenTlY of $\sQrt{D_I}$ and $\sqrt{D_j}$), and a fOurrIer series of the angles, whEre thE cOEffiCIeNts arE functiONs OF $\lamBdA_i$ And $\laMBdA_j$ (i.e., oF thE semi-maJor axes). For A syStEm that is closE To The resonanCe oR resonant, the SemI-major AxEs rEMain cLosE to thE nominAL resonAnt values (kePLer’s thIRd law) $$\laBEL{eQ:nomInalsma}
  \frac{a_i}{a_j} \apPRoX \Frac{a_{i,0}}{a_{j,0}} = \left(\frAc{k_{i,j}}{k_{J,i}}\RIgHT)^{2/3} \LefT(\frAc{\mu_i}{\mu_j}\riGhT)^{1/3}.$$ I intROduce $$\DeLTa\lAMBda_I = \lambda_i - \Lambda_{I,0},$$ where $$\LambdA_{I,0} = \beTa_i\sqrT{\mU_i a_{I,0}},$$ And expAnd thE KEPleRian part at dEgreE 2, and the peRturbaTIve part AT degree 0 In $\DeltA\LaMbdA_i$ $$\beGIn{AlIgnEd}
\LAbeL{Eq:HamDlaD}
  \h &=& \sum_{i=1}^n n_{i,0} \deLtA\LambDa_i
  - \fRAC{3}{2}\FRac{n_{I,0}}{\LaMbda_{I,0}}\DeltA\Lambda_i^2\nonumBer\\
  &+& \Sum_{1\lEQ i<j\Leq n} \H_{I,j}(D_i,D_J,\lamBdA_i,\lamBda_j,\vaRpi_i,\vArPi_j),\end{aligned}$$ whEre $n_{I,0}$ is the nomInaL mEan MoTion oF Planet $I$, suCh tHat $$\frac{N_{i,0}}{n_{j,0}} = \fraC{K_{j,i}}{K_{i,J}}.$$ tHE pErturbative part doeS nOT DePend on $\LaMbda_i$ aNYmOrE, But is simPlY evAluaTED at $\LaMbda_{I,0}$.

in Order to pErform THe ReDuctionS aSsociaTeD to The ConseRVatiOn of anGular momEntum ANd to the averagiNG, I first change THe SYStEM of cOorDinates. For tHe saKE of rEadaBIlIty, i PreseNt the GeNErAL case (with any number oF pLanets, In any Resonance of anY
_{i,j}$ as a power seriesof the ecc entri cit ies  ( or e quiv alently of $\s q rt{D _i}$ and $\sqrt{D_j}$) , and  a Four r ie r ser ies oft he a ngl es ,whe re th e coe ffi cientsare functi ons  o f $\Lambda_i $  a nd $\Lambd a_j $ (i.e., ofthe  semi- ma jor axes) . F or asystem that i s close t ot he res o nance o r re sona nt, the semi-majo r  a x es remain clos e to t he no m i nal  re sonant val ue s (Ke p ler’s t h ir d l aw) $$\label{eq:n ominalsma} \f rac{a_ i} {a_ j } \app rox \ fr a c{a _{i,0}}{a_{ j,0} } = \left (\frac { k_{i,j} } {k_{j,i }}\rig ht) ^{2 /3}\ le ft (\f ra c {\m u _i }{\ m u_j }\right) ^{ 1/ 3}.$$  I i n t r o duce  $$ \Del ta\La mbda_i = \Lam bda _i - \La mbda_ {i,0} ,$$wh ere $ $\Lamb da_{i ,0 } = \beta_i\sqr t{\m u_i a_{i, 0}} ,$ $ a nd  expa n d theKep ler ian par t at de g ree  2 , a nd  the perturbativepa r t  a t degree  0 in$ \D el t a\Lambda _i $ $ $\be g i n{ali gned } 
\ label{eq :hamdL a D} 
 \H &=&  \ sum_{i =1 }^n  n_ {i,0} \Del ta\Lam bda_i
 - \fr a c{3}{2}\frac{n _ {i,0}}{\Lambd a _{ i , 0} } \Del ta\ Lambda_i^2\ nonu m ber\ \
 & +&  \s u m_{1\ leq i <j \ le q  n} \H_{i,j}(D_i,D_ j, \lambd a_i,\ lambda_j,\var pi_i,\varp i _ j ),\end{a lign e d} $ $ where $n_{i, 0}$ i s the nomi n al meanmotio n of pla net $i$,s u ch that$$\ fra c{n _{i , 0 }} {n_{j,0}} = \ f r ac{k _{ j,i}}{k _{i ,j}}.$$  Th e p ert urb at ive partdoes not  d ep en don$\Lam b da_i$ an ym ore ,but  is s i mply e valua tedat  $ \ Lam bda_{i, 0 }$ . 
Inor de r to  pe rf orm t he r e duc tions a ssociated  to theco ns ervatio n of angularmo mentum and  t o t he ave r a ging, Ifirst change the system of coor din ates.  For  the sake  of  reada bil i ty, Ipresen t the  g ene r a l cas e (w ith  a ny numbero f  pl anets ,in a ny reso nance of any
_{i,j}$ as_a power_series of the eccentricities_(or equivalently_of_$\sqrt{D_i}$ and_$\sqrt{D_j}$),_and a Fourrier_series of the_angles, where the coefficients_are functions of_$\Lambda_i$_and $\Lambda_j$ (i.e., of the semi-major axes). For a system that is close to_the_resonance or_resonant,_the_semi-major axes remain close to_the nominal resonant values (Kepler’s_third law)_$$\label{eq:nominalsma}
  \frac{a_i}{a_j} \approx \frac{a_{i,0}}{a_{j,0}} = \left(\frac{k_{i,j}}{k_{j,i}}\right)^{2/3} \left(\frac{\mu_i}{\mu_j}\right)^{1/3}.$$_I_introduce $$\Delta\Lambda_i =_\Lambda_i - \Lambda_{i,0},$$ where $$\Lambda_{i,0} = \beta_i\sqrt{\mu_i a_{i,0}},$$ and_expand the Keplerian part at degree_2, and the_perturbative_part_at degree 0 in_$\Delta\Lambda_i$ $$\begin{aligned}
\label{eq:hamdLaD}
  \H &=& \sum_{i=1}^n_n_{i,0} \Delta\Lambda_i
  - \frac{3}{2}\frac{n_{i,0}}{\Lambda_{i,0}}\Delta\Lambda_i^2\nonumber\\
 _&+& \sum_{1\leq i<j\leq n} \H_{i,j}(D_i,D_j,\lambda_i,\lambda_j,\varpi_i,\varpi_j),\end{aligned}$$ where $n_{i,0}$_is the nominal mean motion of_planet $i$, such that $$\frac{n_{i,0}}{n_{j,0}}_= \frac{k_{j,i}}{k_{i,j}}.$$_The perturbative part does not_depend on $\Lambda_i$_anymore, but_is simply evaluated_at $\Lambda_{i,0}$.

In order to perform the_reductions associated to_the conservation of angular momentum and_to_the averaging, I_first_change_the system_of coordinates. For_the_sake of_readability,_I present the general case (with_any_number of planets, in any resonance of_any