Explicit expansion of the three-body disturbing function for arbitrary eccentricities and inclinations

J. Laskar; G. Boué

doi:10.1051/0004-6361/201014496

Home

All issues

Volume 522 (November 2010)

A&A, 522 (2010) A60

Full HTML

Free Access

Issue		A&A Volume 522, November 2010


Article Number		A60
Number of page(s)		11
Section		Celestial mechanics and astrometry
DOI		https://doi.org/10.1051/0004-6361/201014496
Published online		03 November 2010

A&A 522, A60 (2010)

Explicit expansion of the three-body disturbing function for arbitrary eccentricities and inclinations

J. Laskar and G. Boué

ASD, IMCCE-CNRS UMR8028, Observatoire de Paris, UPMC, 77 avenue Denfert-Rochereau, 75014 Paris, France
e-mail: This email address is being protected from spambots. You need JavaScript enabled to view it.

Received: 24 March 2010
Accepted: 27 June 2010

Abstract

Since the original work of Hansen and Tisserand in the XIXth century, there have been many variations in the analytical expansion of the three-body disturbing function in series of the semi-major axis ratio. With the increasing number of planetary systems of large eccentricity, these expansions are even more interesting as they allow us to obtain for the secular systems finite expressions that are valid for all eccentricities and inclinations. We revisited the derivation of the disturbing function in Legendre polynomial, with a special focus on the secular system. We provide here expressions of the disturbing function for the planar and spatial case at any order with respect to the ratio of the semi-major axes. Moreover, for orders in the ratio of semi-major axis up to ten in the planar case and five in the spatial case, we provide explicit expansions of the secular system, and simple algorithms with minimal computation to extend this to higher order, as well as the algorithms for the computation of non secular terms.

Key words: celestial mechanics / planets and satellites: dynamical evolution and stability / planetary systems / methods: analytical

© ESO, 2010

1. Introduction

In the three-body problem, there are two classical ways to compute the principal part of the disturbing function −Gmm′/Δ. The first approach is to expand it in powers of eccentricity and inclination, with coefficients that are expressed in term of Laplace coefficients (e.g. Laplace 1785; Abu-El-Ata & Chapront 1975; Laskar & Robutel 1995), but this approach, which is well suited to the study of the Solar System, has some limitations for some extra-solar planetary systems, where the eccentricity can reach very high values. Another drawback, is that without some special truncation corrections, the angular momentum will not be conserved exactly in the truncated system.

In the second approach, the expansion is made with respect to the ratio of the semi-major axes of the two bodies α = a/a′, where a′ is related to the external body. Although it may be less efficient for low eccentricity and large values of α, as for the inner planets in the Solar system, the advantage is that this expansion allows us to obtain finite expressions for arbitrary eccentricities in the secular system, while both approaches allow expansions for arbitrary inclinations. The most important contribution to this problem was made in the XIXth century (Hansen 1855; Hill 1875; Tisserand 1899). With the development of space technology, expansions of the disturbing function in term of Legendre polynomials have been rejuvenated for the construction of satellite theories in the vicinity of the Earth (Kozai 1959; Kaula 1962; Brumberg 1967; Brumberg et al. 1971; Giacaglia 1974; Abu-El-Ata & Chapront 1975).

The discovery of new planetary systems, starting with 51Pegb (Mayor & Queloz 1995), has raised the need to revise these methods as many systems have planets with very high eccentricity, as GL 581, HD 217107, HD 69830, HD 74156, HD 168443, HD 102272, HD 169830, HD 202206, HD 183263, or even HD 80606, where the eccentricity reaches 0.93¹. With the discovery of these numerous new planetary systems, the previous analytical expansions in Laplace coefficients developed for the Solar System are no longer the most appropriate, and we are faced with the need to understand more globally the dynamics of these systems. On the other hand, as the parameters of these extrasolar planetary systems remain not very well known, there is no necessity for very precise analytical approximations. This has led several authors to use the Legendre expansion of the potential for the dynamical study of the secular planetary system, following the previous studies of stellar systems (Krymolowski & Mazeh 1999; Ford et al. 2000; Blaes et al. 2002), where the secular spatial three-body system was expanded up to the octupolar order (α³). In particular, Lee & Peale (2003) used the secular planar system at octupolar order to study the secular dynamics of the HD 168443 and HD 12661 systems. Migaszewski & Goździewski (2008) computed the planar secular system to high order using computer algebra to average over the mean anomaly.

As it appears that there is a growing interest in analytical studies of extra solar planetary systems, we considered it interesting to present a derivation of the planetary (or stellar) disturbing function in a very simple and explicit way that does not already appear in the existing literature. We aimed to write a self-contained paper that allows one to construct explicitly the planetary Hamiltonian to high order, with minimal additional computation. In particular, in Sect. 3 we show that the planar secular system can be obtained explicitly at any order without the need for computer algebra. The spatial case is treated in Sect. 4 with minimal computation when expressed with respect to the mutual inclination J. In the presence of more than two planets, it may be more convenient to use a fixed reference frame, and the derivation of the spatial Hamiltonian is provided in this case in Sect. 4.2. Our derivation is close to the original formulations of Hansen (1855) and Tisserand (1899), but with emphasis on the secular system and the direct derivation of explicit expressions. We also present in Appendix A a new derivation of the Hansen coefficients for the secular terms.

2. Expansion of the disturbing function

To simplify the notations, and although everything can be generalized to an N-body system, we consider here a three-body problem with a central body of mass M and two other celestial bodies of masses m,m′. If we write the Hamiltonian of the Newtonian interactions among these three bodies in Poincaré canonical heliocentric variables, we obtain (Laskar & Robutel 1995, Eq. (15))

$ℋ = ℋ_{0} + 𝒯_{1} + 𝒰_{1},$ $Mathematical equation: \begin{equation} {\cal H} = {\cal H}_0 + {\cal T}_1 + {\cal U}_1 , \end{equation}$ (1)

where

$ℋ_{0} = - \frac{G (M + m)}{2 a} - \frac{G (M + m^{'})}{2 a^{'}}$ $Mathematical equation: \begin{equation} {\cal H}_0 = -\frac{G(M+m)}{2a} -\frac{G(M+m')}{2a'} \end{equation}$ (2)

is the Keplerian interaction with elliptical elements a,e,i,M,ω, and Ω that denote, respectively, the semi-major axis, eccentricity, inclination, mean anomaly, argument of perihelion, and longitude of the node (with primes for the external body of mass m′). The principal part U₁ part of the perturbation and indirect part T₁ are then

$𝒰_{1} = - G \frac{m m^{'}}{∥ r - r^{'} ∥}; 𝒯_{1} = \frac{\begin{matrix} ˜ \\ r \end{matrix} \cdot {\begin{matrix} ˜ \\ r \end{matrix}}^{'}}{M},$ $Mathematical equation: \begin{equation} {\cal U}_1 = - G \frac{mm'}{\norm{\br-\br'}} ; \quad {\cal T}_1 = \frac{\tilde\br\cdot \tilde\br'}{M} , \label{eq.hamt} \end{equation}$ (3)

where r,r′ are the radius vectors of the inner and outer planets, with norms r,r′, unit vectors u = r/r and u′ = r′/r′, and conjugate momentum $Mathematical equation: $\tilde\br, \tilde\br'$$ . We focus first on the principal part of the Hamiltonian U₁, which is the most difficult part to compute, while the computation of the indirect part will be made in Sect. 5. With ρ = r/r′, and $Mathematical equation: $\cF = a'/\norm{\br-\br'}$$ , we have

$ℱ = \frac{a^{'}}{r^{'}} {(1 + ρ^{2} - 2 ρ u \cdot u^{'})}^{- 1 / 2} = \frac{a^{'}}{r^{'}} \sum_{n = 0}^{\infty} P_{n} (u \cdot u^{'}) ρ^{n},$ $Mathematical equation: \begin{equation} \cF= \frac{a'}{r'}\left(1+\rho^2 -2\rho \uu\dpp\uu'\right)^{-1/2} = \frac{a' }{r' }\sum_{n=0}^\infty P_n(\uu\dpp\uu') \rho^n, \end{equation}$ (4)

where P_n(x) are the Legendre polynomials that can be written as

$P_{n} (z) = \sum_{k = 0}^{[n / 2]} p_{n,k} z^{n - 2 k}$ $Mathematical equation: \begin{equation} P_n(z) = \sum_{k=0}^{[n/2]} p_{n,k}\, z^{n-2k} \end{equation}$ (5)

with

$p_{n,k} = \frac{(- 1)^{k}}{2^{n}} \frac{(2 n - 2 k)!}{k! (n - k)! (n - 2 k)!} \cdot$ $Mathematical equation: \begin{equation} p_{n,k} = \frac{(-1)^k}{2^n} \frac{(2n-2k)!}{k! (n-k)! (n-2k)!} \cdot \label{eq.pnk} \end{equation}$ (6)

With α = a/a′,γ = r/a,γ′ = r′/a′, we have ρ = αγ/γ′, and thus, as γ,γ′,u,u′ do not depend on a,a′, the expansion of ℱ in powers of α is

$ℱ = \sum_{n = 0}^{\infty} ℱ_{n} α^{n}; with ℱ_{n} = \frac{γ^{n}}{γ^{' n + 1}} F_{n}$ $Mathematical equation: \begin{equation} \cF = \sum_{n=0}^\infty \cF_n \, \al^n; \quad \hbox{with}\ \cF_n= \frac{\ga^n}{\ga'^{n+1}} F_n \end{equation}$ (7)

and

$F_{n} = P_{n} (u \cdot u^{'}) = \sum_{k = 0}^{[n / 2]} p_{n,k} (u \cdot u^{'})^{n - 2 k} .$ $Mathematical equation: \begin{equation} F_n=P_n(\uu\dpp\uu')=\sum_{k=0}^{[n/2]} p_{n,k}\, (\uu\dpp\uu')^{n-2k}. \label{eq.fn} \end{equation}$ (8)

3. Planar case

Table 1

Tisserand functions for the planar case (Eq. (8)).

We first study the planar case that leads to some simplifications in the expansions. We define v,v′ to be the true anomalies of u,u′, and ω,ω′ their argument of perihelion, u = v + ω, u′ = v′ + ω′, and x = u − u′. We have then u·u′ = cosx and, following a classical computation (e.g. Whittaker & Watson 1927, p. 303), we have for all z ∈ [ 0,1 [ ,

$\begin{matrix} (1 - 2 z \cos x + z^{2})^{- 1 / 2} & = & (1 - z e^{i x})^{- 1 / 2} (1 - z e^{- i x})^{- 1 / 2} \\ = & \sum_{p = 0}^{\infty} \sum_{q = 0}^{\infty} \frac{(2 p)! (2 q)!}{2^{2 p + 2 q} (p!)^{2} (q!)^{2}} e^{i (p - q) x} z^{p + q} \end{matrix}$ $Mathematical equation: \begin{eqnarray} % (1-2z\cos x \! +\! z^2)^{-1/2} & =& (1-z \,\e^{\i x})^{-1/2}(1-z \,\e^{-\i x})^{-1/2} \nonumber \\ &=& \!\!\sum_{p=0}^{\infty}\sum_{q=0}^{\infty}\frac{(2p)!(2q)!}{2^{2p+2q}(p!)^2(q!)^2} \,\e^{\i(p-q)x} z^{p+q} % % \end{eqnarray}$ (9)

and, with n = p + q, and after changing q to n − q,

$= \sum_{n = 0}^{\infty} (\sum_{q = 0}^{n} \frac{(2 q)! (2 n - 2 q)!}{2^{2 n} (q!)^{2} ((n - q)!)^{2}} e^{i (2 q - n) x}) z^{n} .$ $Mathematical equation: \begin{equation} \qquad\qquad \quad\!\! \!= \sum_{n=0}^{\infty} \left(\sum_{q=0}^n \frac{(2q)!(2n-2q)!}{2^{2n}(q!)^2((n-q)!)^2} \,\e^{\i(2q-n)x}\right)z^n. \end{equation}$ (10)

Thus

$F_{n} = \sum_{q = 0}^{n} f_{n,q} e^{i (2 q - n) x}$ $Mathematical equation: \begin{equation} F_n= \sum_{q=0}^{n} f_{n,q} \e^{\i(2q-n)x} \label{eq.fnx} \end{equation}$ (11)

with

$f_{n,q} = \frac{(2 q)! (2 n - 2 q)!}{2^{2 n} (q!)^{2} ((n - q)!)^{2}}$ $Mathematical equation: \begin{equation} f_{n,q}= \frac{(2q)!(2n-2q)!}{2^{2n}(q!)^2((n-q)!)^2} \label{eq.fnq} \end{equation}$ (12)

for 0 ≤ q ≤ n. If we write x = v − v′ + ω − ω′, ℱ_n can now be expressed in the form

$ℱ_{n} = \sum_{q = 0}^{n} ℱ_{n,q} e^{i (2 q - n) (ω - ω^{'})}$ $Mathematical equation: \begin{equation} \cF_n= \sum_{q=0}^{n} \cF_{n,q} \, \e^{\i(2q-n)(\om-\om')} \end{equation}$ (13)

with

$ℱ_{n,q} = f_{n,q} \frac{γ^{n}}{γ^{' n + 1}} e^{i (2 q - n) (v - v^{'})} .$ $Mathematical equation: \begin{equation} \cF_{n,q} = f_{n,q} \frac{\ga^n}{\ga'^{n+1}} \e^{\i(2q-n)(v-v')}. \end{equation}$ (14)

The quantities ℱ_n,q can then be expressed in term of Hansen coefficients $Mathematical equation: $X_k^{n,m}(e)$$ defined for n,m ∈ Z as

${(\frac{r}{a})}^{n} e^{i mv} = \sum_{k = - \infty}^{+ \infty} X_{k}^{n,m} (e) e^{i kM} .$ $Mathematical equation: \begin{equation} \left(\frac{r}{a}\right)^n \e^{\i mv} = \sum_{k=-\infty}^{+\infty} X_k^{n,m}(e) \, \e^{\i kM}. \end{equation}$ (15)

For convenience, we denote $Mathematical equation: $X_k^{n,m}=X_k^{n,m}(e)$$ and $Mathematical equation: ${X'}_k^{n,m}=X_k^{n,m}(e')$$ . We have thus

$ℱ_{n,q} = f_{n,q} \sum_{k, k^{'} = - \infty}^{+ \infty} X_{k}^{n, 2 q - n} {X^{'}}_{k^{'}}^{- (n + 1), - 2 q + n} e^{i (kM + k^{'} M^{'})} .$ $Mathematical equation: \begin{equation} \cF_{n,q} = f_{n,q} \sum_{k,k'=-\infty}^{+\infty} X_k^{n,2q-n} {X'}_{k'}^{-(n+1),-2q+n}\e^{\i(kM+k'M')} . \end{equation}$ (16)

For arbitrary values of k ∈ Z, the Hansen coefficients $Mathematical equation: $X_k^{n,m}$$ can be expressed in an explicit manner as an infinite series involving Bessel functions and hypergeometric functions (Hansen 1855; Tisserand 1899), or in terms of generalized Laplace coefficients (Laskar 2005), but for k = 0, $Mathematical equation: $X_k^{n,m}$$ reduces to a finite polynomial expression in e, 1/e, $Mathematical equation: $\sqrt{1-e^2}$$ , and $Mathematical equation: $1/\sqrt{1-e^2}$$ (see Appendix A). We have thus a very compact expression for the coefficient of any argument kM + k′M′ in explicit form at all orders n in α. Indeed, if we denote this coefficient $Mathematical equation: $\cF_n^{(k,k')}$$ , that is

$ℱ_{n} = \sum_{k, k^{'} = - \infty}^{+ \infty} ℱ_{n}^{(k, k^{'})} e^{i (kM + k^{'} M^{'})},$ $Mathematical equation: \begin{equation} \cF_n = \sum_{k,k'=-\infty}^{+\infty} \cF_n^{(k,k')} \e^{\i(kM+k'M')} , \label{eq.fkk} \end{equation}$ (17)

we have

$ℱ_{n}^{(k, k^{'})} = \sum_{q = 0}^{n} f_{n,q} X_{k}^{n, 2 q - n} {X^{'}}_{k^{'}}^{- (n + 1), - 2 q + n} e^{i (2 q - n) (ω - ω^{'})} .$ $Mathematical equation: \begin{equation} \cF_n^{(k,k')} = \sum_{q=0}^{n} f_{n,q} X_k^{n,2q-n} {X'}_{k'}^{-(n+1), -2q+n} \e^{\i(2q-n)(\om-\om')}. \end{equation}$ (18)

In particular, for all n ∈ N, the secular part $Mathematical equation: $\cF_n^{(0,0)}$$ can even be simplified, using the classical relation $Mathematical equation: $X_{-k}^{n,-m}= X_k^{n,m}$$ among Hansen coefficients and the relation f_n,q = f_n,n − q from Eq. (12)

$\begin{matrix} ℱ_{n}^{(0, 0)} & = & ϵ_{n} f_{n, \frac{n}{2}} X_{0}^{n, 0} {X^{'}}_{0}^{- (n + 1), 0} \\ + \sum_{q = 0}^{[(n - 1) / 2]} 2 f_{n,q} X_{0}^{n,n - 2 q} {X^{'}}_{0}^{- (n + 1),n - 2 q} \cos ((n - 2 q) (ω - ω^{'})), \end{matrix}$ $Mathematical equation: \begin{eqnarray} % \cF_n^{(0,0)} &=& \epsilon_n f_{n,\frac{n}{2}} X_0^{n,0} {X'}_{0}^{-(n+1), 0} \nonumber \\ \label{eq.perplan} &&\hspace*{-1.2cm} +\sum_{q=0}^{[(n-1)/2]} 2\,f_{n,q} X_0^{n,n-2q} {X'}_{0}^{-(n+1), n-2q} \cos((n-2q)(\om-\om') ), \end{eqnarray}$ (19)

where ϵ_n = 0 if n is odd, and ϵ_n = 1 if n is even.

3.1. Practical algorithm

Equations (12) and (19) provide an explicit algorithm for the computation of the principal part of the disturbing function ℱ, and in particular for the computation of the secular Hamiltonian at any order N in α

$⟨ ℱ ⟩ = \sum_{n = 0}^{N} ℱ_{n}^{(0, 0)} α^{n},$ $Mathematical equation: \begin{equation} \langle \cF \rangle = \sum_{n=0}^N \cF_n^{(0,0)}\alpha^n , \end{equation}$ (20)

for which the involved Hansen coefficients $Mathematical equation: $ X_0^{n,m}, {X'}_{0}^{-(n+1), m}$$ reduce to finite expressions (see Appendix A). As there are no secular terms in the indirect part of the Hamiltonian (see Laskar & Robutel 1995), the computation of the secular Hamiltonian reduces to the computation of the rational constants f_n,k given by Eq. (12). In practice, for finite order n, it is even easier to use the expression of F_n (11) (Table 1) and to translate it into ℱ_n (19) by the simple transformation

$\begin{matrix} \frac{γ^{n}}{γ^{' n + 1}} \cos (mx) = \sum_{k, k^{'} = - \infty}^{+ \infty} X_{k}^{n,m} {X^{'}}_{k^{'}}^{- (n + 1), - m} \cos (kM + k^{'} M^{'} + m (ω - ω^{'})), \end{matrix}$ $Mathematical equation: \begin{eqnarray} \frac{\ga^n}{\ga'^{n+1}} \cos (mx) = \sum_{{k,k'=-\infty}}^{+\infty} X_k^{n,m} {X'}_{k'}^{-(n+1), -m} \cos(kM+k'M'+m(\om-\om') ) , \label{eq.transf} \end{eqnarray}$ (21)

which allows an immediate computation of any argument $Mathematical equation: $\cF_n^{(k,k')}$$ in terms of Hansen coefficients.

3.2. Computation of the secular part $Mathematical equation: $\cF_{\sf n}^{(0,0)}$$

The computation of the secular part of the Hamiltonian $Mathematical equation: $\cF_n^{(0,0)}$$ is quite simple, as we will just have to make the transformation

$\cos (mx) - \to X_{0}^{n,m} {X^{'}}_{0}^{- (n + 1),m} \cos (m (ω - ω^{'}))$ $Mathematical equation: \begin{equation} \cos (mx) \longrightarrow X_0^{n,m} {X'}_{0}^{-(n+1), m} \cos(m(\om-\om') ) \end{equation}$ (22)

in the expression of F_n. Moreover, we have $Mathematical equation: ${X'}_{0}^{-(n+1), m}= 0$$ for m ≥ n ≥ 1 (see Appendix A). The term in cosnx can thus be discarded in the expression of F_n (Table 1) which simplifies the expression of the Hamiltonian. The secular Hamiltonian is thus expressed in finite form, using the values of the Hansen coefficients given in Appendix A. We have

$\begin{matrix} ℱ_{0}^{(0, 0)} & = & 1 \\ ℱ_{1}^{(0, 0)} & = & 0 \\ ℱ_{2}^{(0, 0)} & = & \frac{1}{4} X_{0}^{2, 0} {X^{'}}_{0}^{- 3, 0} \\ ℱ_{3}^{(0, 0)} & = & \frac{3}{8} X_{0}^{3, 1} {X^{'}}_{0}^{- 4, 1} \cos (ω - ω^{'}) \\ ℱ_{4}^{(0, 0)} & = & \frac{9}{64} X_{0}^{4, 0} {X^{'}}_{0}^{- 5, 0} + \frac{5}{16} X_{0}^{4, 2} {X^{'}}_{0}^{- 5, 2} \cos (2 (ω - ω^{'})) \\ ℱ_{5}^{(0, 0)} & = & \frac{15}{64} X_{0}^{5, 1} {X^{'}}_{0}^{- 6, 1} \cos (ω - ω^{'}) \\ + \frac{35}{128} X_{0}^{5, 3} {X^{'}}_{0}^{- 6, 3} \cos (3 (ω - ω^{'})) \\ ℱ_{6}^{(0, 0)} & = & \frac{25}{256} X_{0}^{6, 0} {X^{'}}_{0}^{- 7, 0} + \frac{105}{512} X_{0}^{6, 2} {X^{'}}_{0}^{- 7, 2} \cos (2 (ω - ω^{'})) \\ + \frac{63}{256} X_{0}^{6, 4} {X^{'}}_{0}^{- 7, 4} \cos (4 (ω - ω^{'})) \\ ℱ_{7}^{(0, 0)} & = & \frac{175}{1024} X_{0}^{7, 1} {X^{'}}_{0}^{- 8, 1} \cos (ω - ω^{'}) \\ + \frac{189}{1024} X_{0}^{7, 3} {X^{'}}_{0}^{- 8, 3} \cos (3 (ω - ω^{'})) \\ + \frac{231}{1024} X_{0}^{7, 5} {X^{'}}_{0}^{- 8, 5} \cos (5 (ω - ω^{'})) . \end{matrix}$ $Mathematical equation: \begin{eqnarray*} \cF_0^{(0,0)} &=& 1\\ \cF_1^{(0,0)} &=& 0\\ \cF_2^{(0,0)} &=& \frac{1}{4} X_0^{2,0}{X'}_{0}^{-3, 0} \\ \cF_3^{(0,0)} &=& \frac{3}{8} X_0^{3,1} {X'}_{0}^{-4, 1} \cos(\om-\om' ) \\ \cF_4^{(0,0)} &=& \frac{9}{64} X_0^{4,0} {X'}_{0}^{-5, 0} + \frac{5}{16} X_0^{4,2} {X'}_{0}^{-5, 2} \cos(2(\om-\om') )\\ \cF_5^{(0,0)} &=& \frac{15}{64} X_0^{5,1} {X'}_{0}^{-6, 1} \cos( \om-\om') \\ &&+ \frac{35}{128} X_0^{5,3} {X'}_{0}^{-6, 3} \cos(3(\om-\om') ) \\ \cF_6^{(0,0)} &=& \frac{25}{256}X_0^{6,0} {X'}_{0}^{-7, 0} +\frac{105}{512} X_0^{6,2} {X'}_{0}^{-7, 2} \cos( 2(\om-\om'))\\ &&+ \frac{63}{256} X_0^{6,4} {X'}_{0}^{-7, 4} \cos(4(\om-\om') ) \\ \cF_7^{(0,0)} &=& \frac{175}{1024} X_0^{7,1} {X'}_{0}^{-8, 1} \cos(\om-\om') \\ &&+ \frac{189}{1024} X_0^{7,3} {X'}_{0}^{-8, 3} \cos(3(\om-\om')) \\ &&+ \frac{231}{1024} X_0^{7,5} {X'}_{0}^{-8, 5} \cos(5(\om-\om')) . \end{eqnarray*}$

These results are equivalent to the expression obtained by Migaszewski & Goździewski (2008) using computer algebra.

4. Spatial case

The spatial case is more complicated as it involves additional variables. Our goal is to derive explicit formulae that are as compact as possible. We thus expand in terms of the mutual inclination J. For each orbit, we use a reference frame (i,j,k) associated with the orbit, with first vector i in the direction of the ascending node of r′ over r. With u = v + ω and u′ = v′ + ω′, we have

$\begin{matrix} u \cdot u^{'} = \cos u \cos u^{'} + \cos J \sin u \sin u^{'} = μ \cos x + ν \cos y \end{matrix}$ $Mathematical equation: \begin{equation} \EQM{ \uu\dpp\uu' &=& \cos u \cos u' + \cos J \sin u \sin u' \\ &=& \mu \cos x + \nu \cos y \label{eq.uupJ} } \end{equation}$ (23)

with the same notations as Tisserand (1885)

$x = u - u^{'}; y = u + u^{'}; μ = \cos^{2} \frac{J}{2}; ν = \sin^{2} \frac{J}{2} \cdot$ $Mathematical equation: \begin{equation} x = u-u' ;\quad y= u+u' ;\quad \mu = \cos^2 \frac{J}{2} ;\quad \nu = \sin^2 \frac{J}{2} \cdot \end{equation}$ (24)

As in the planar case, we have F_n = P_n(u·u′), but now u·u′ is given by the slightly more complex expression (23). For all n, we have

$F_{n} = \sum_{s = 0}^{n} \sum_{q = 0}^{n} 𝒬_{s,q}^{(n)} (μ,ν) e^{i u (n - 2 s)} e^{i u^{'} (n - 2 q)},$ $Mathematical equation: \begin{equation} F_n = \sum_{s=0}^{n} \sum_{q= 0}^{n } \cQ^{(n)}_{s,q} (\mu,\nu) \, \e^{\i u\,(n-2s)}\e^{\i u'\,(n-2q)}, \label{eq.fnuu} \end{equation}$ (25)

where the $Mathematical equation: $ \cQ^{(n)}_{s,q} (\mu,\nu)$$ are polynomials in μ,ν of degree n that are called the Tisserand polynomials² as a recognition of the work of Tisserand (1885), although these expressions are already present in (Hansen 1855). As (u,u′) − → (− u, − u′) leaves u·u′ unchanged, we have $Mathematical equation: $ \cQ^{(n)}_{n-s,n-q} = \cQ^{(n)}_{s,q} $$ , thus F_n can be expressed as a trigonometric polynomial in cos(mu + m′u′). Although they can be computed explicitly for all n (see Appendix B), for a given n, it is often more efficient to make a direct computation of F_n on a computer algebra system. For example, F₂₀ is computed in less than 1.5 s in exact rational arithmetics with TRIP (Gastineau & Laskar 2009) on an average laptop computer using the simple expression (8). We can then express ℱ_n = (γⁿ/γ′^n + 1)F_n in term of Hansen coefficients. We thus have

$ℱ_{n} = \sum_{k = - \infty}^{+ \infty} \sum_{k^{'} = - \infty}^{+ \infty} ℱ_{n}^{(k, k^{'})} e^{i (kM + k^{'} M^{'})},$ $Mathematical equation: \begin{equation} \cF_n = \sum_{k=-\infty}^{+\infty}\sum_{k'=-\infty}^{+\infty} \cF_n^{(k,k')} \e^{\i(kM+k'M')}, \end{equation}$ (26)

with

$ℱ_{n}^{(k, k^{'})} = \sum_{s,q = 0}^{n} 𝒬_{s,q}^{(n)} X_{k}^{n,n - 2 s} {X^{'}}_{k^{'}}^{- (n + 1),n - 2 q} e^{i (n - 2 s) ω} e^{i (n - 2 q) ω^{'}} .$ $Mathematical equation: \begin{equation} \cF_n^{(k,k')} = \sum_{s,q=0}^{n} \cQ^{(n)}_{s,q} \, X_k^{n,n-2s}\, {X'}_{k'}^{-(n+1),n-2q}\,\e^{\i\,(n-2s)\om}\e^{\i\,(n-2q)\om'} . \end{equation}$ (27)

More practically, starting from the finite expressions in cosine polynomials (Table 2) of the Tisserand functions F_n, the expression of ℱ_n is obtained as in the planar case by means of the more general transformation

$\begin{matrix} \frac{γ^{n}}{γ^{' n + 1}} \cos (mu + m^{'} u^{'}) & = & \sum_{k = - \infty}^{+ \infty} \sum_{k^{'} = - \infty}^{+ \infty} X_{k}^{n,m} {X^{'}}_{k^{'}}^{- (n + 1), m^{'}} \\ \times \cos (kM + k^{'} M^{'} + mω + m^{'} ω^{'}) . \end{matrix}$ $Mathematical equation: \begin{eqnarray} % \frac{\ga^n}{\ga'^{n+1}}\cos (mu +m'u')&=&\sum_{k=-\infty}^{+\infty}\sum_{{k'=-\infty}}^{+\infty} X_k^{n,m} {X'}_{k'}^{-(n+1), m'} \nonumber \\ \label{eq.transfiiii} & &\times \cos(kM+k'M'+m\om+m'\om') . \end{eqnarray}$ (28)

4.1. Computation of the secular part $Mathematical equation: $\cF_n^{(0,0)}$$

As in the planar case, the computation of the secular part of the Hamiltonian $Mathematical equation: $\cF_n^{(0,0)}$$ is just a straightforward translation of F_n in Table 2, with the transformation

$\cos (mu + m^{'} u^{'}) - \to X_{0}^{n,m} {X^{'}}_{0}^{- (n + 1), m^{'}} \cos (mω + m^{'} ω^{'}) .$ $Mathematical equation: \begin{equation} \cos (mu+m'u') \longrightarrow X_0^{n,m} {X'}_{0}^{-(n+1), m'} \cos(m\om+m'\om') \ . \end{equation}$ (29)

Moreover, as $Mathematical equation: ${X'}_0^{-(n+1),n} = 0$$ for n ≥ 1, all terms in cos(mu ± nu′) can be discarded in F_n. We thus have,

$\begin{matrix} ℱ_{0}^{(0, 0)} & = & 1 \\ ℱ_{1}^{(0, 0)} & = & 0 \\ ℱ_{2}^{(0, 0)} & = & (- \frac{1}{2} + \frac{3}{4} ν^{2} + \frac{3}{4} μ^{2}) X_{0}^{2, 0} {X^{'}}_{0}^{- 3, 0} + \frac{3}{2} νμ X_{0}^{2, 2} {X^{'}}_{0}^{- 3, 0} \cos (2 ω) \\ ℱ_{3}^{(0, 0)} & = & X_{0}^{3, 1} {X^{'}}_{0}^{- 4, 1} [(- \frac{3}{2} μ + \frac{15}{4} ν^{2} μ + \frac{15}{8} μ^{3}) \cos (ω - ω^{'}) \\ + (- \frac{3}{2} ν + \frac{15}{8} ν^{3} + \frac{15}{4} ν μ^{2}) \cos (ω + ω^{'})] \\ + \frac{15}{8} X_{0}^{3, 3} {X^{'}}_{0}^{- 4, 1} [ν^{2} μ \cos (3 ω + ω^{'}) + ν μ^{2} \cos (3 ω - ω^{'})^{]} \\ ℱ_{4}^{(0, 0)} & = & X_{0}^{4, 0} {X^{'}}_{0}^{- 5, 0} (+ \frac{3}{8} - \frac{15}{8} ν^{2} \\ + \frac{105}{64} ν^{4} - \frac{15}{8} μ^{2} + \frac{105}{16} ν^{2} μ^{2} + \frac{105}{64} μ^{4}) \\ + X_{0}^{4, 2} {X^{'}}_{0}^{- 5, 2} [(- \frac{15}{8} μ^{2} \\ + \frac{105}{16} ν^{2} μ^{2} + \frac{35}{16} μ^{4}) \cos (2 ω - 2 ω^{'}) \\ + (- \frac{15}{8} ν^{2} + \frac{35}{16} ν^{4} + \frac{105}{16} ν^{2} μ^{2}) \cos (2 ω + 2 ω^{'})] \\ + (- \frac{15}{4} νμ + \frac{105}{16} ν^{3} μ + \frac{105}{16} ν μ^{3}) X_{0}^{4, 2} {X^{'}}_{0}^{- 5, 0} \cos (2 ω) \\ + (- \frac{15}{4} νμ + \frac{105}{16} ν^{3} μ + \frac{105}{16} ν μ^{3}) X_{0}^{4, 0} {X^{'}}_{0}^{- 5, 2} \cos (2 ω^{'}) \\ + X_{0}^{4, 4} {X^{'}}_{0}^{- 5, 2} [+ \frac{35}{16} ν μ^{3} \cos (4 ω - 2 ω^{'}) \\ + \frac{35}{16} ν^{3} μ \cos (4 ω + 2 ω^{'})] \\ + \frac{105}{32} X_{0}^{4, 4} {X^{'}}_{0}^{- 5, 0} ν^{2} μ^{2} \cos (4 ω) . \end{matrix}$ $Mathematical equation: \begin{eqnarray*} % \cF_{0}^{(0,0)} \!&=&\!1\nonumber\\ \cF_{1}^{(0,0)} \!&=&\!0\nonumber\\ \cF_{2}^{(0,0)} &=& \left(- \frac{1}{2} + \frac{3}{4} \nu^2 + \frac{3}{4} \mu^2 \right)X_0^{2,0}{X'}_{0}^{-3, 0} \!+ \! \frac{3}{2} \nu \mu X_0^{2,2}{X'}_{0}^{-3, 0} \cos(2\om) \nonumber\\[1mm] \cF_{3}^{(0,0)} &=&X_0^{3,1} {X'}_{0}^{-4,1} \left[ \left(- \frac{3}{2} \mu + \frac{15}{4} \nu^2 \mu + \frac{15}{8} \mu^3 \right) \cos(\om-\om') \right. \nonumber\\[1mm] &&+\left.\left(- \frac{3}{2} \nu + \frac{15}{8} \nu^3 + \frac{15}{4} \nu \mu^2\right) \cos(\om+\om')\right] \nonumber\\[1mm] && +\frac{15}{8}X_0^{3,3} {X'}_{0}^{-4,1} \left[ \nu^2 \mu \cos(3\om+\om') + \nu \mu^2 \cos(3\om-\om')\right]\nonumber\\ % % % \cF_{4}^{(0,0)} &=&X_0^{4,0}{X'}_{0}^{-5, 0}\left( + \frac{3}{8} - \frac{15}{8} \nu^2 \right.\nonumber \\ && \left. + \frac{105}{64} \nu^4 - \frac{15}{8} \mu^2 + \frac{105}{16} \nu^2 \mu^2 + \frac{105}{64} \mu^4\right) \nonumber\\ && + X_0^{4,2}{X'}_{0}^{-5, 2}\left[\left( - \frac{15}{8} \mu^2 \right. \right. \nonumber\\ && \left. + \frac{105}{16} \nu^2 \mu^2 + \frac{35}{16} \mu^4\right) \cos(2\om-2\om') \nonumber \\ &&+\left.\left( - \frac{15}{8} \nu^2 + \frac{35}{16} \nu^4 + \frac{105}{16} \nu^2 \mu^2 \right)\cos(2\om+2\om')\right]\nonumber\\[1mm] && +\left( - \frac{15}{4} \nu \mu + \frac{105}{16} \nu^3 \mu + \frac{105}{16} \nu \mu^3\right) X_0^{4,2}{X'}_{0}^{-5, 0}\cos(2\om) \nonumber\\[1mm] && +\left( - \frac{15}{4} \nu \mu + \frac{105}{16} \nu^3 \mu + \frac{105}{16} \nu \mu^3\right) X_0^{4,0}{X'}_{0}^{-5, 2} \cos(2\om')\nonumber\\[1mm] && +X_0^{4,4}{X'}_{0}^{-5,2}\left[ + \frac{35}{16} \nu \mu^3 \cos(4\om-2\om') \right.\nonumber\\ &&\left. + \frac{35}{16} \nu^3 \mu \cos(4\om+2\om')\right]\nonumber\\ &&+ \frac{105}{32}X_0^{4,4}{X'}_{0}^{-5, 0} \nu^2 \mu^2 \cos(4\om).\nonumber \end{eqnarray*}$

We thus observe here an important simplification in the quadrupolar secular Hamiltonian $Mathematical equation: $\cF_{2}^{(0,0)}$$ , as all terms involving the external planet longitude of perihelion ω′ vanish and we are left with an integrable Hamiltonian. This is what Lidov & Ziglin (1976) called a happy coincidence (see also Farago & Laskar 2010). This is no longer the case at higher orders.

4.2. Expression in a fixed reference frame

Table 2

Tisserand functions for the spatial case.

Table 3

Tisserand functions for the spatial case in a fixed reference frame.

The above expressions are given with respect to the mutual inclination to shorten the algebraic expansions. This is especially useful in a three-body problem, but to obtain expressions valid in a fixed reference frame, then one needs to substitute into u·u′ its expression in terms of the elliptical elements of the two bodies. We can then generalize the expression (23) and write it now as (Brumberg 1967; Abu-El-Ata & Chapront 1975; Laskar & Robutel 1995)

$u \cdot u^{'} = ℜ (μ_{*} e^{i x} + ν_{*} e^{i y})$ $Mathematical equation: \begin{equation} \uu\dpp\uu' = \Re \left( \tmu \,\e^{\i x} + \tnu \,\e^{\i y}\right) \end{equation}$ (30)

with as before x = u − u′, y = u + u′, and

$\begin{matrix} μ_{*} & = {(c c^{'} e^{i \frac{Ω - Ω^{'}}{2}} + s s^{'} e^{- i \frac{Ω - Ω^{'}}{2}})}^{2} \\ ν_{*} & = {(c s^{'} e^{i \frac{Ω - Ω^{'}}{2}} - s c^{'} e^{- i \frac{Ω - Ω^{'}}{2}})}^{2} \end{matrix}$ $Mathematical equation: \begin{equation} \begin{array}{ll} \tmu &= \left( \co \co'\e^{\i\frac{\Om-\Om'}{2}} + \si\si' \e^{-\i\frac{\Om-\Om'}{2}} \right)^2 \\ \tnu &= \left( \co \si'\e^{\i\frac{\Om-\Om'}{2}} - \si\co' \e^{-\i\frac{\Om-\Om'}{2}} \right)^2 \end{array} \end{equation}$ (31)

with c = cos(i/2), s = sin(i/2), and the same for the primes. We write

$μ_{*} = a + i a^{'}; ν_{*} = b + i b^{'}$ $Mathematical equation: \begin{equation} \tmu = a + {\rm i} a' \ ; \quad \tnu = b+{\rm i}b' \end{equation}$ (32)

with

$\begin{matrix} a & = (c^{2} {c^{'}}^{2} + s^{2} {s^{'}}^{2}) \cos (Ω - Ω^{'}) + 2 c c^{'} s s^{'} \\ a^{'} & = (c^{2} {c^{'}}^{2} - s^{2} {s^{'}}^{2}) \sin (Ω - Ω^{'}) \\ b & = (c^{2} {s^{'}}^{2} + s^{2} {c^{'}}^{2}) \cos (Ω - Ω^{'}) - 2 c c^{'} s s^{'} \\ b^{'} & = (c^{2} {s^{'}}^{2} - s^{2} {c^{'}}^{2}) \sin (Ω - Ω^{'}) . \end{matrix}$ $Mathematical equation: \begin{equation} \begin{array}{ll} a &= (\co^2{\co'}^2+\si^2{\si'}^2) \cos (\Om-\Om') + 2 \co\co'\si\si' \\ a' &= (\co^2{\co'}^2-\si^2{\si'}^2) \sin (\Om-\Om') \\ b &= (\co^2{\si'}^2+\si^2{\co'}^2) \cos (\Om-\Om') - 2 \co\co'\si\si' \\ b' &= (\co^2{\si'}^2-\si^2{\co'}^2) \sin (\Om-\Om'). \end{array} \label{eq.aabb} \end{equation}$ (33)

With these notations, we have

$| ν_{*} | = \sin^{2} (J / 2), | μ_{*} | = \cos^{2} (J / 2)$ $Mathematical equation: \begin{equation} \abs{\tnu} = \sin^2(J/2), \qquad \abs{\tmu} = \cos^2(J/2) \label{eq.tnutmu} \end{equation}$ (34)

and

$\begin{matrix} u \cdot u^{'} = a \cos (u - u^{'}) - a^{'} \sin (u - u^{'}) \\ + b \cos (u + u^{'}) - b^{'} \sin (u + u^{'}) . \end{matrix}$ $Mathematical equation: \begin{equation} \EQM{ \uu\dpp\uu' &= &a\cos(u-u') -a'\sin(u-u')\crm &&+b\cos(u+u')-b'\sin(u+u') . \label{eq.uuii} } \end{equation}$ (35)

The remaining part is identical to the previous case. The expressions of the Tisserand functions F_n are given for 0 ≤ n ≤ 3 in Table 3 and for all n in Appendix B. In practice, for small values of n we can use explicit expressions of F_n, using the straightforward approach consisting of computing F_n directly with computer algebra, using the relation (8), and then to translate it as previously in order to obtain the expression of any argument kM + k′M′ by means of the relations

$\begin{matrix} \frac{γ^{n}}{γ^{' n + 1}} & {\begin{matrix} \cos \\ \sin \end{matrix} (mu + m^{'} u^{'}) = \sum_{k = - \infty}^{+ \infty} \sum_{k^{'} = - \infty}^{+ \infty} X_{k}^{n,m} \end{matrix}$ $Mathematical equation: \begin{eqnarray} \frac{\ga^n}{\ga'^{n+1}}&&\!\!\!\!\! \left\{\EQM{\cos\cr\sin}\right. ( m u+m'u')=\sum_{k=-\infty}^{+\infty}\sum_{{k'=-\infty}}^{+\infty} X_k^{n,m} \crm && \times\, {X'}_{k'}^{-(n+1), m'} \left\{\EQM{\cos\cr\sin}\right.(kM+k'M'+m\om+m'\om'). \label{eq.transfii} \end{eqnarray}$ (36)

We thus have for the secular part $Mathematical equation: $\cF_n^{(0,0)}$$ ,

$\begin{matrix} ℱ_{0}^{(0, 0)} & = & 1 \\ ℱ_{1}^{(0, 0)} & = & 0 \\ ℱ_{2}^{(0, 0)} & = & X_{0}^{2, 0} {X^{'}}_{0}^{- 3, 0} [- \frac{1}{2} + \frac{3}{4} (b^{' 2} + b^{2} + a^{' 2} + a^{2})] \\ + X_{0}^{2, 2} {X^{'}}_{0}^{- 3, 0} \frac{3}{2} [(ba - b^{'} a^{'}) \cos (2 ω) - (b a^{'} + b^{'} a) \sin (2 ω)] \end{matrix}$ $Mathematical equation: \begin{eqnarray*} % \cF_{0}^{(0,0)} &=&1\nonumber\\ \cF_{1}^{(0,0)} &=&0\nonumber\\ \cF_{2}^{(0,0)} &=& X_0^{2,0}{X'}_{0}^{-3, 0}\left[ -\frac{1}{2} +\frac{3}{4}( b'^2 + b^2 + a'^2 + a^2) \right]\nonumber\\ &&\!\!\!\! \!\!\! + \, X_0^{2,2}{X'}_{0}^{-3, 0}\, \frac{3}{2}[ ( b a - b' a') \cos(2 \om) -( b a' + b' a )\sin(2 \om) ]\nonumber % \end{eqnarray*}$

and

$\begin{matrix} ℱ_{3}^{(0, 0)} & = & \frac{15}{8} \\ \times X_{0}^{3, 1} {X^{'}}_{0}^{- 4, 1} [+ (2 b^{' 2} a + 2 b^{2} a + a^{' 2} a + a^{3} - \frac{4}{5} a) \cos (ω - ω^{'}) \\ + (2 b a^{' 2} + 2 b a^{2} - \frac{4}{5} b + b^{' 2} b + b^{3}) \cos (ω + ω^{'}) \\ + (\frac{4}{5} b^{'} - b^{' 3} - b^{'} b^{2} - 2 b^{'} a^{' 2} - 2 b^{'} a^{2}) \sin (ω + ω^{'}) \\ - (a^{'} a^{2} - \frac{4}{5} a^{'} + 2 b^{' 2} a^{'} + 2 b^{2} a^{'} + a^{' 3}) \sin (ω - ω^{'})] \\ + \frac{15}{8} X_{0}^{3, 3} {X^{'}}_{0}^{- 4, 1} [(b^{2} a - 2 b^{'} b a^{'} - b^{' 2} a) \cos (3 ω + ω^{'}) \\ + (b a^{2} - b a^{' 2} - 2 b^{'} a^{'} a) \cos (3 ω - ω^{'}) \\ + (b^{' 2} a^{'} - b^{2} a^{'} - 2 b^{'} ba) \sin (3 ω + ω^{'}) \\ - (2 b a^{'} a + b^{'} a^{2} - b^{'} a^{' 2}) \sin (3 ω - ω^{'})] . \end{matrix}$ $Mathematical equation: \begin{eqnarray*} \cF_{3}^{(0,0)}&=& \frac{15}{8}\nonumber\\ &&\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\times\, X_0^{3,1}{X'}_{0} ^{-4, 1}\left[+ \left(2 b'^2 a +2 b^2 a + a'^2 a + a^3 -\frac{4}{5} a \right) \cos(\om-\om') \right.\nonumber\\ &&+\, \left(2 b a'^2 +2 b a^2 -\frac{4}{5} b + b'^2 b + b^3 \right)\cos(\om+\om')\nonumber\\ &&+\, \left(\frac{4}{5} b' - b'^3 - b' b^2 -2 b' a'^2 -2b' a^2 \right)\sin(\om+\om')\nonumber\\ &&-\, \left.\left( a' a^2 -\frac{4}{5} a' +2 b'^2 a' +2 b^2 a' + a'^3 \right) \sin(\om-\om')\right]\nonumber\\ &&+\, \frac{15}{8}X_0^{3,3}{X'}_{0}^{-4, 1}[ ( b^2 a -2 b' b a' - b'^2 a )\cos(3 \om+\om')\nonumber\\ &&+\, (b a^2 - b a'^2 -2 b' a' a )\cos(3 \om-\om')\nonumber\\[1mm] &&+\, ( b'^2 a' - b^2 a' -2 b' b a )\sin(3 \om+\om')\nonumber\\[1mm] &&-\, (2 b a' a + b' a^2 - b' a'^2 )\sin(3 \om-\om')] .\nonumber \end{eqnarray*}$

5. Indirect part

Until now, we considered only the principal part of the perturbing Hamiltonian. The computation of the indirect part T₁ (3) is more straightforward, as it compares with the computation of r·r′. Indeed, we have (Laskar & Robutel 1995, Eq. (24)),

$𝒯_{1} = \frac{m m^{'}}{M} V with V = \dot{r} \cdot \dot{r}',$ $Mathematical equation: \begin{equation} {\cal T}_1 = \frac{mm'}{M} V \quad \hbox{with} \quad V=\dot\br\dpp\dot\br' , \end{equation}$ (37)

where $Mathematical equation: $\dot\br,\dot\br'$$ are the velocities in the corresponding Keplerian problem. With classical computation, we obtain

$\begin{matrix} \frac{a}{r} \cos E = \frac{2}{e} \sum_{k = 1}^{+ \infty} J_{k} (ke) \cos kM \frac{a}{r} \sin E = 2 \sum_{k = 1}^{+ \infty} J_{k}^{'} (ke) \sin kM \end{matrix}$ $Mathematical equation: \begin{equation} \EQM{ \frac{a}{r} \cos E &=& \frac{2}{e} \sum_{k=1}^{+\infty} J_{k}(ke) \cos kM \\ \frac{a}{r} \sin E &=& 2 \sum_{k=1}^{+\infty} J'_{k}(ke) \sin kM } \end{equation}$ (38)

where E is the eccentric anomaly and J_k(x) are the Bessel functions. The coordinates (Ẋ,Ẏ) of the velocity $Mathematical equation: $\dot\br$$ in the reference frame of the orbit with origin at perihelion are then easily expressed in Fourier series of the mean anomaly, as

$\begin{matrix} Ẋ = - na (\frac{a}{r} \sin E) Ẏ = na \sqrt{1 - e^{2}} (\frac{a}{r} \cos E) . \end{matrix}$ $Mathematical equation: \begin{equation} \EQM{ \dot X = -na \left(\frac{a}{r} \sin E\right) \\ \dot Y = na \sqrt{1-e^2} \left(\frac{a}{r} \cos E\right) . } \end{equation}$ (39)

If we denote 𝒵 = Ẋ + iẎ, we then have for the spatial problem in the fixed reference frame, with the same notations as in the previous section (Laskar & Robutel 1995, Eq. (37)),

$V = ℜ (μ_{*} 𝒵 𝒵̅' e^{i (ω - ω^{'})} + ν_{*} 𝒵 𝒵^{'} e^{i (ω + ω^{'})}) .$ $Mathematical equation: \begin{equation} V= \Re\left(\tmu \cZ \bar\cZ' \e^{\i(\om-\om')} + \tnu\cZ\cZ' \e^{\i(\om+\om')}\right) . \end{equation}$ (40)

If one considers the mutual inclination J, as in Sect. 4, this expression simplifies since Ω = Ω′, s′ = 0, c′ = 1, and μ_*,ν_* are then real, with μ_* = cos²(J/2), ν_* = sin²(J/2). In the planar case (Sect. 3), μ_* = 1, ν_* = 0.

We considered here heliocentric coordinates. One could also use Jacobi coordinates for the three-body problem. In this case, the indirect part does not require additional computations as it is expressed in term of u·u′ (see Laskar 1990). It should be noted that in terms of both heliocentric coordinates or Jacobi coordinates, the indirect part does not contribute to the secular system. We have provided here the expression of the indirect part for the computation of non-secular inequalities.

6. Conclusion

We have presented a self-contained exposition of the expansion of the three-body Hamiltonian in canonical heliocentric coordinates in term of Legendre polynomials at any order in the ratio of semi-major axes α, which can also be adapted in the case of Jacobi coordinates, where only the indirect part differs. We have included here all the necessary material that allows one to write explicitly the secular Hamiltonian at order α¹⁰ for the planar case, α⁵ for the spatial case expressed in terms of the mutual inclination, and α³ for the spatial case in a fixed reference frame. With the additional computation of the required Hansen coefficients, the expressions of the Tisserand functions F_n can also be used for a straightforward computation of the expression of a non-secular inequality kM + k′M′.

As the algorithms that are presented here are very simple, we have not added to this paper any tables in electronic form. Indeed, the reader who needs to use expressions of high order that do not appear in the paper, will have no problems in exploiting the algorithms given here to derive the required expressions. The written results of the paper can then be used to check his computations for the lowest orders. For example, the computation of F₁₀₀, in the spatial case takes less than 8 min on an average 8 core desktop computer in exact rational arithmetics using TRIP, but has 2 343 926 terms, while F₅₀ needs only 7.5 s with already 164 151 terms. As the algorithms require only a few lines of code (fewer than 20 in our case), one understands that it is preferable to compute the terms when they are needed than to store the values in electronic form.

See http://exoplanet.eu

One should consult (Aksenov 1986) for a detailed discussion of the relation between the Tisserand polynomials and the inclination functions of Kaula (1962).

Appendix A: Computation of the Hansen coefficients

Table A.1

$Mathematical equation: $\tilde X_0^{-n,m}$$ is the polynomial part of the Hansen coefficients $Mathematical equation: $X_0^{-n,m}$$ .

Table A.2

Hansen coefficients $Mathematical equation: $X_0^{n,m}$$ for (0 ≤ n ≤ 10,0 ≤ m ≤ n).

The Hansen coefficients $Mathematical equation: \appendix \setcounter{section}{1} $X_k^{n,m}(e)$$ are defined as the Fourier coefficients

${(\frac{r}{a})}^{n} e^{i mv} = \sum_{k = - \infty}^{+ \infty} X_{k}^{n,m} (e) e^{i kM} .$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{equation} \left(\frac{r}{a}\right)^n \e^{\i mv} = \sum_{k=-\infty}^{+\infty} X_k^{n,m}(e) \, \e^{\i kM}. \end{equation}$ (A.1)

For arbitrary values of k ∈ Z, the Hansen coefficients $Mathematical equation: \appendix \setcounter{section}{1} $X_k^{n,m}(e)$$ can be expressed in an explicit manner as an infinite series involving Bessel functions and hypergeometric functions (Hansen 1855; Hill 1875; Tisserand 1899), or generalized Laplace coefficients (Laskar 2005), but for k = 0, $Mathematical equation: \appendix \setcounter{section}{1} $X_k^{n,m}(e)$$ reduces to a simple polynomial in e, 1/e, $Mathematical equation: \appendix \setcounter{section}{1} $\sqrt{1-e^2}$$ , and $Mathematical equation: \appendix \setcounter{section}{1} $1/\sqrt{1-e^2}$$ . The literature on the computation of the Hansen coefficients has been huge since the original work of Hansen (1855), but we do not review it here. We concentrate on the obtention of the coefficients $Mathematical equation: \appendix \setcounter{section}{1} $ X_0^{n,m}(e)$$ and $Mathematical equation: \appendix \setcounter{section}{1} $ X_0^{-n,m}(e)$$ , for n ≥ 0 and 0 ≤ m ≤ n that are required to compute the averaged planetary or lunar Hamiltonian. Using

$d M = \frac{r^{2}}{a^{2} \sqrt{1 - e^{2}}} d v = \frac{r}{a} d E$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{equation} {\rm d}M = \frac{r^2}{a^2\sqrt{1-e^2}} {\rm d}v = \frac{r}{a} {\rm d}E \end{equation}$ (A.2)

where v is the true anomaly and E the eccentric anomaly, we have for n ≥ 2

$\begin{matrix} X_{0}^{- n,m} & = & \frac{1}{\sqrt{1 - e^{2}}} \frac{1}{2 π} \int_{0}^{2 π} {(\frac{a}{r})}^{n - 2} e^{i mv} d v \\ = & \frac{1}{(1 - e^{2})^{n - 3 / 2}} \frac{1}{2 π} \int_{0}^{2 π} (1 + e \cos v)^{n - 2} e^{i mv} d v \\ = & \frac{1}{(1 - e^{2})^{n - 3 / 2}} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{eqnarray} X_0^{-n,m} &=& \frac{1}{\sqrt{1-e^2}} \frac{1}{2\pi}\int_0^{2\pi} \left(\frac{a}{r}\right)^{n-2} \e^{\i mv} {\rm d}v \nonumber \\[2mm] &=&\frac{1}{(1-e^2)^{n-3/2}} \frac{1}{2\pi}\int_0^{2\pi} (1+e\cos v)^{n-2} \e^{\i mv} {\rm d}v\nonumber \\[2mm] &=&\frac{1}{(1-e^2)^{n-3/2}} \nonumber \\[2mm] &&\times \sum_{l=0}^{[(n-2-m)/2]}\!\!\!\!\! \frac{(n-2)!}{l!\,(m+l)!\,(n-2-(m+2l))!} \left(\frac{e}{2}\right)^{m+2l} \cdot \label{eq.anomv} \end{eqnarray}$ (A.3)

We still need to consider the case n = 1 for which the expansion in true anomaly v is not suitable. We have immediately, using (A.6)

$X_{0}^{- 1, 0} = 1; X_{0}^{- 1, 1} = \frac{\sqrt{1 - e^{2}} - 1}{e} \cdot$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{equation} X_0^{-1,0}=1; \qquad X_0^{-1,1} = \frac{\sqrt{1-e^2}-1}{e} \cdot \end{equation}$ (A.4)

It is important to note that from (A.3), we have

$X_{0}^{- n,m} = 0 for n \geq 2 and m \geq n - 1.$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{equation} X_0^{-n,m} = 0 \quad \hbox{for} \quad n\geq 2 \quad \hbox{and} \quad m \geq n-1 . \end{equation}$ (A.5)

The Hansen coefficients $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{-n,m}$$ are given in Table A.1 for n ≥ 2 and m ≤ n − 2. On the other hand, the computation of $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,m}(e)$$ for n ≥ 0 is not as straightforward, as a direct expansion of

$\begin{matrix} X_{0}^{n,m} & = & \frac{1}{2 π} \int_{0}^{2 π} {(\frac{r}{a})}^{n + 1} e^{i mv} d E \\ = & \frac{1}{2 π} \int_{0}^{2 π} (1 - e \cos E)^{n + 1 - m} (\cos E - e \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{eqnarray} X_0^{n,m} &= & \frac{1}{2\pi}\int_0^{2\pi} \left(\frac{r}{a}\right)^{n+1} \e^{\i mv} {\rm d}E \nonumber \\[2mm] &=& \frac{1}{2\pi}\int_0^{2\pi} (1-e\cos E)^{n+1-m} (\cos E -e\nonumber \\[2mm] && +\, {\rm i}\sqrt{1-e^2}\sin E)^m {\rm d}E \label{eq.anomE} \end{eqnarray}$ (A.6)

in eccentric anomaly is far more complicated than the previous expression (A.3). One can still perform some formal expansion, but this would not be very useful, as it would contain expressions with many summations that cannot be easily reduced to a single summation as in (A.3). In fact, the only explicit computation of $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,m}(e)$$ available in the literature, was obtained through a complex process in (Hansen 1855; Hill 1875; Tisserand 1899), using the auxiliary variable $Mathematical equation: \appendix \setcounter{section}{1} $\beta = ({1-\sqrt{1-e^2}})/{e}$$ . In this case, it can be shown that $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,m}(e)$$ is expressed as a finite hypergeometric series in β². Using a relation among hypergeometric series due to Gauss, and the change of variable e = 2β/(1 + β²), Tisserand (1899) provides a finite expression of $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,m}(e)$$ in term of hypergeometric series of e². We present here a more direct demonstration for the obtention of an explicit expression of $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,m}(e)$$ which is a recurrence using the relation among Hansen coefficients

$X_{0}^{n,m} = \frac{1}{n + 2 - m} (\frac{2}{e} (m - 1) X_{0}^{n,m - 1} + (n + m) X_{0}^{n,m - 2}) .$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{equation} X_0^{n,m} = \frac{1}{n+2-m} \left( \frac{2}{e} (m-1) X_0^{n,m-1} + (n+m) X_0^{n,m-2} \right). \label{eq.rec} \end{equation}$ (A.7)

This relation appears in (Hughes 1981), but it was shown by Laskar (2005) that it is equivalent to a recurrence relation obtained on Laplace coefficients by Laplace (1785). The computation of $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,0}$$ for n ≥ 0 using (A.6) is straightforward and gives

$X_{0}^{n, 0} = \sum_{l = 0}^{[(n + 1) / 2]} \frac{(n + 1)!}{l! l! (n + 1 - 2 l)!} {(\frac{e}{2})}^{2 l} \cdot$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{equation} X_0^{n,0} = \sum_{l=0}^{[(n+1)/2]} \frac{(n+1)!}{l!\,l!\,(n+1-2l)!} \left(\frac{e}{2}\right)^{2l}\cdot \end{equation}$ (A.8)

The computation of $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,1}$$ can also be performed using (A.6). Changing E to − E, it is immediate to see that the part in sinE cancel, and we are left with

$X_{0}^{n, 1} = - (n + 2) \sum_{l = 0}^{[n / 2]} \frac{n!}{l! (l + 1)! (n - 2 l)!} {(\frac{e}{2})}^{2 l + 1} \cdot$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{equation} X_0^{n,1} =- (n+2)\sum_{l=0}^{[n/2]} \frac{n!}{l!\,(l+1)!\,(n-2l)!} \left(\frac{e}{2}\right)^{2l+1}\cdot \end{equation}$ (A.9)

We can now prove by recurrence using the relation (A.7) that the general form of $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,m}$$ for n ≥ 0 and 0 ≤ m ≤ n is

$\begin{matrix} X_{0}^{n,m} & = & (- 1)^{m} \frac{(n + 1 + m)!}{(n + 1)!} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{1} \begin{eqnarray} X_0^{n,m} &=&(-1)^m \frac{(n+ 1+m)!}{(n+1)!} \nonumber\\[2.5mm] &&\times \sum_{l=0}^{[(n+1-m)/2]} \frac{(n+1-m)!}{l!\,(m+l)!\,(n+1- m-2l)!} \left(\frac{e}{2}\right)^{m+2l}\cdot \label{eq.xnm} \end{eqnarray}$ (A.10)

The computation is delicate but straightforward. One can first show that the two first elements (for l = 0) of the polynomial expressions of $Mathematical equation: \appendix \setcounter{section}{1} $ X_0^{n,m-1}$$ and $Mathematical equation: \appendix \setcounter{section}{1} $ X_0^{n,m-2}$$ in (A.7) cancel. We then change the index from l to l + 1 and show that the general term of the sum in (A.7) gives the general term of $Mathematical equation: \appendix \setcounter{section}{1} $ X_0^{n,m}$$ , and that the last term of $Mathematical equation: \appendix \setcounter{section}{1} $ X_0^{n,m-2}$$ will give the last term of $Mathematical equation: \appendix \setcounter{section}{1} $ X_0^{n,m}$$ . The expression (A.10) is equivalent to the expression of (Tisserand 1899).The expressions of the Hansen coefficients $Mathematical equation: \appendix \setcounter{section}{1} $X_0^{n,m}$$ are given in Table A.2 for 0 ≤ n ≤ 10 and 0 ≤ m ≤ n.

Appendix B: Computation of Tisserand polynomials and Tisserand functions

Expansions valid for all inclinations were already given by Hansen (1855), and Tisserand (1885) in his researches on asteroidal motions, but these derivations are rather complex. Although in (B.3) we present some of Tisserand’s computations, we make first, as in the planar case some direct expansions and reordering of terms to obtain expressions that are as compact as possible.

B.1. Tisserand functions with respect to mutual inclination

With x = u − u′, y = u + u′, μ = cos²(J/2), ν = sin²(J/2), we have

$u \cdot u^{'} = μ \cos x + ν \cos y,$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} \uu\dpp\uu' = \mu \cos x + \nu \cos y , \label{eq.uuxy} \end{equation}$ (B.1)

and with a straightforward expansion in complex notation, we obtain

$(u \cdot u^{'})^{m} = \frac{m!}{2^{m}} \sum_{\begin{matrix} k_{1} + k_{2} + k_{3} + k_{4} = m \\ 0 \leq k_{i} \leq m \end{matrix}} \frac{μ^{k_{1} + k_{2}} ν^{k_{3} + k_{4}}}{k_{1}! k_{2}! k_{3}! k_{4}!} e^{i x (k_{1} - k_{2})} e^{i y (k_{3} - k_{4})},$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} (\uu\dpp\uu')^m = \frac{m!}{2^m} \!\! \!\! \sum_{\tiny \begin{array}{c}{k_1+k_2+k_3+k_4=m} \\ {0\leq k_i\leq m}\end{array}}\!\!\!\! \frac{\mu^{k_1+k_2}\nu^{k_3+k_4}}{k_1!k_2!k_3!k_4!} \e^{\i x(k_1-k_2)}\e^{\i y(k_3-k_4)} ,\; \end{equation}$ (B.2)

which we can reorder as

$(u \cdot u^{'})^{m} = \frac{m!}{2^{m}} \sum_{r = 0}^{m} \sum_{p = - r}^{m - r} 𝒞_{p,m - 2 r - p}^{(m)} e^{i xp} e^{i y (m - 2 r - p)}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} (\uu\dpp\uu')^m = \frac{m!}{2^m} \sum_{r=0}^{m} \sum_{p=-r}^{m-r} \cC_{p,m-2r-p}^{(m)}\, \e^{\i x\,p}\e^{\i y\,(m-2r-p)} \end{equation}$ (B.3)

with

$𝒞_{p,m - 2 r - p}^{(m)} = \sum_{l = \max (0, - p)}^{\min (r,m - r - p)} \frac{μ^{2 l + p} ν^{m - 2 l - p}}{l! (l + p)! (r - l)! (m - p - r - l)!} \cdot$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} \cC_{p,m-2r-p}^{(m)} = \sum_{l=\max(0,-p)}^{\min(r,m-r-p)} \frac{\mu^{2l+p}\nu^{m-2l-p}}{l!(l+p)!(r-l)!(m-p-r-l)!}\cdot \label{eq.Cmpq} \end{equation}$ (B.4)

With the expansion (8), we have then

$F_{n} = \sum_{k = 0}^{[n / 2]} q_{n,k} \sum_{r = 0}^{n - 2 k} \sum_{p = - r}^{n - 2 k - r} 𝒞_{p,n - 2 k - 2 r - p}^{(n - 2 k)} e^{i xp} e^{i y (n - 2 k - 2 r - p)}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} F_n = \sum_{k=0}^{[n/2]} q_{n,k}\, \sum_{r=0}^{n-2k} \sum_{p=-r}^{n-2k-r} \cC_{p,n-2k-2r-p}^{(n-2k)}\, \e^{\i x\,p}\e^{\i y\,(n-2k-2r-p)} \label{eq.fkrp} \end{equation}$ (B.5)

with

$q_{n,k} = p_{n,k} \frac{(n - 2 k)!}{2^{n - 2 k}} = \frac{(- 1)^{k}}{2^{(2 n - 2 k)}} \frac{(2 n - 2 k)!}{k! (n - k)!} \cdot$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} q_{n,k}=p_{n,k}\, \frac{(n-2k)!}{2^{n-2k}} = \frac{(-1)^k}{2^{(2n-2k)}} \frac{(2n-2k)!}{k! (n-k)!} \cdot \label{eq.qnk} \end{equation}$ (B.6)

We can then reorder (B.5) with s = k + r as

$F_{n} = \sum_{s = 0}^{n} \sum_{k = 0}^{\min (s,n - s)} \sum_{p = k - s}^{n - k - s} q_{n,k} 𝒞_{p,n - 2 s - p}^{(n - 2 k)} e^{i xp} e^{i y (n - 2 s - p)},$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} F_n = \sum_{s=0}^{n} \sum_{k=0}^{\min(s,n-s)}\sum_{p=k-s}^{n-k-s} q_{n,k}\, \cC_{p,n-2s-p}^{(n-2k)}\, \e^{\i x\,p}\e^{\i y\,(n-2s-p)} , \label{eq.fskp} \end{equation}$ (B.7)

exchange the sum in p and k

$F_{n} = \sum_{s = 0}^{n} \sum_{p = - s}^{n - s} \sum_{k = 0}^{\min (s + p,n - s - p,s,n - s)} q_{n,k} 𝒞_{p,n - 2 s - p}^{(n - 2 k)} e^{i xp} e^{i y (n - 2 s - p)}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} F_n = \sum_{s=0}^{n} \sum_{p= -s}^{n -s} \sum_{k=0}^{\min(s+p,n-s-p,s,n-s)}\!\!\!\! q_{n,k}\, \cC_{p,n-2s-p}^{(n-2k)}\, \e^{\i x\,p}\e^{\i y\,(n-2s-p)} \label{eq.fskpii} \end{equation}$ (B.8)

and set q = p + s which gives

$F_{n} = \sum_{s = 0}^{n} \sum_{q = 0}^{n} (\sum_{k = 0}^{\min (q,n - q,s,n - s)} q_{n,k} 𝒞_{q - s,n - q - s}^{(n - 2 k)}) e^{i x (q - s)} e^{i y (n - q - s)} .$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} F_n = \sum_{s=0}^{n} \sum_{q= 0}^{n } \left( \sum_{k=0}^{\min(q,n-q,s,n-s)}\!\!\!\! q_{n,k} \, \cC_{q-s,n-q-s}^{(n-2k)}\right) \, \e^{\i x\,(q-s)}\e^{\i y\,(n-q-s)}. \label{eq.fskpiiii} \end{equation}$ (B.9)

That is, with x = u − u′ and y = u + u′,

$F_{n} = \sum_{s = 0}^{n} \sum_{q = 0}^{n} 𝒬_{s,q}^{(n)} (μ,ν) e^{i u (n - 2 s)} e^{i u^{'} (n - 2 q)}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} F_n = \sum_{s=0}^{n} \sum_{q= 0}^{n } \cQ^{(n)}_{s,q} (\mu,\nu) \, \e^{\i u\,(n-2s)}\e^{\i u'\,(n-2q)} \label{eq.fnuuB} \end{equation}$ (B.10)

with

$𝒬_{s,q}^{(n)} (J) = \sum_{k = 0}^{\min (q,n - q,s,n - s)} q_{n,k} 𝒞_{q - s,n - q - s}^{(n - 2 k)},$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} \cQ^{(n)}_{s,q}(J) = \sum_{k=0}^{\min(q,n-q,s,n-s)} q_{n,k}\, \cC_{q-s,n-q-s}^{(n-2k)} , \label{eq.tis} \end{equation}$ (B.11)

that is

$\begin{matrix} 𝒬_{s,q}^{(n)} (J) & = & \sum_{k = 0}^{\min (q,n - q,s,n - s)} q_{n,k} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} \cQ^{(n)}_{s,q}(J) &= &\sum_{k=0}^{\min(q,n-q,s,n-s)} q_{n,k} \nonumber \\ &&\hspace*{-1.2cm}\times \sum_{l=\max(0,s-q)}^{\min(s-k,n-k-q)} \frac{\mu^{2l+q-s}\nu^{n-2k-q+s-2l }}{l!(l+q-s)!(s-k-l)!(n-k-q-l)!} \cdot \label{eq.tisi} \end{eqnarray}$ (B.12)

B.2. Tisserand functions with respect to a fixed reference frame

In the case of a fixed reference frame (Sect. 4.2), the expressions differ slightly, but are more complicated. We have (35)

$u \cdot u^{'} = a \cos x - a^{'} \sin x + b \cos y - b^{'} \sin y,$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} \uu\dpp\uu' = a\cos x -a'\sin x + b\cos y -b'\sin y , \end{equation}$ (B.13)

that is, with μ_* = a + ia′, ν_* = b + ib′,

$u \cdot u^{'} = \frac{1}{2} [μ_{*} e^{i x} + μ_{*} e^{- i x} + ν_{*} e^{i y} + ν_{*} e^{- i y}] .$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} \uu\dpp\uu' = \frac{1}{2}\left[ \tmu \e^{\i x} + \bmus \e^{-\i x}+ \tnu \e^{\i y} + \bnus \e^{-\i y}\right]. \end{equation}$ (B.14)

This is similar to the previous expression (B.1), and the computation of the Tisserand functions F_n will follow the same lines and gives the more general expression

$F_{n} = \sum_{s = 0}^{n} \sum_{q = 0}^{n} \begin{matrix} ˜ \\ 𝒬_{s,q}^{(n)} \end{matrix} e^{i u (n - 2 s)} e^{i u^{'} (n - 2 q)}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} F_n = \sum_{s=0}^{n} \sum_{q= 0}^{n } \tilde{\cQ}^{(n)}_{s,q} \, \e^{\i u\,(n-2s)}\e^{\i u'\,(n-2q)} \label{eq.fnuuii} \end{equation}$ (B.15)

with

$\begin{matrix} \begin{matrix} ˜ \\ 𝒬_{s,q}^{(n)} \end{matrix} (μ_{*}, ν_{*}) & = & \sum_{k = 0}^{\min (q,n - q,s,n - s)} q_{n,k} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} \tilde{\cQ}^{(n)}_{s,q}(\tmu,\tnu) &=& \sum_{k=0}^{\min(q,n-q,s,n-s)} q_{n,k}\, \nonumber\\ &&\hspace*{-1.8cm}\times\sum_{l=\max(0,s-q)}^{\min(s-k,n-k-q)} \frac{\tmu^{l+q-s}\,\bmus^{\ l}\,\tnu^{n-k-q-l }\,\bnus^{\ s-k-l }}{l!(l+q-s)!(s-k-l)!(n-k-q-l)!} \cdot \label{eq.tisii} \end{eqnarray}$ (B.16)

The coefficients $Mathematical equation: \appendix \setcounter{section}{2} $\cQ^{(n)}_{s,q}(J)$$ (B.12) and $Mathematical equation: \appendix \setcounter{section}{2} $\tilde{\cQ}^{(n)}_{s,q}(\tmu,\tnu)$$ (B.16) in the Tisserand functions are very similar. Indeed, using μ = 1 − ν and $Mathematical equation: \appendix \setcounter{section}{2} $\abs{\tmu}=1-\abs{\tnu}$$ (34), we have

$\begin{matrix} 𝒬_{s,q}^{(n)} (J) & = & μ^{q - s} ν^{n - q - s} A_{q - s,n - q - s}^{(n)} (ν), \\ \begin{matrix} ˜ \\ 𝒬_{s,q}^{(n)} \end{matrix} (μ_{*}, ν_{*}) & = \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} \cQ^{(n)}_{s,q}(J) &=& \mu^{q-s}\nu^{n-q-s}A^{(n)}_{q-s,n-q-s}(\nu), \nonumber\\ \tilde{\cQ}^{(n)}_{s,q}(\tmu,\tnu) &=& \tmu^{q-s}\tnu^{n-q-s}A^{(n)}_{q-s,n-q-s}(\abs{\tnu}), \label{eq.Ansq} \end{eqnarray}$ (B.17)

with $Mathematical equation: \appendix \setcounter{section}{2} $\nu=\abs{\tnu}=\sin^2(J/2)$$ and

$\begin{matrix} A_{q - s,n - q - s}^{(n)} (x) = & \sum_{k = 0}^{\min (q,n - q,s,n - s)} q_{n,k} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} A^{(n)}_{q-s,n-q-s}(x) =&& \sum_{k=0}^{\min(q,n-q,s,n-s)} q_{n,k}\nonumber\\ &&\hspace*{-2cm}\times \sum_{l=\max(0,s-q)}^{\min(s-k,n-k-q)}\!\! \frac{(1-x)^{2l}x^{-2k+2s-2l}}{l!(l+q-s)!(s-k-l)!(n-k-q-l)!} \cdot \label{Anqs} \end{eqnarray}$ (B.18)

B.3. Tisserand simplification

The above expression (B.18) is a double sum. In the mutual inclination case, Tisserand (1885) could reduce it to a single sum using hypergeometric functions. He first considers the differential equation satisfied by Legendre polynomials

$(1 - z^{2}) \frac{d^{2} P_{n}}{d z^{2}} - 2 z \frac{d P_{n}}{d z} + n (n + 1) P_{n} = 0.$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} (1-z^2)\frac{{\rm d}^2P_n}{{\rm d}z^2}-2z\frac{{\rm d}P_n}{{\rm d}z}+n(n+1)P_n = 0 . \label{eq.diffPn} \end{equation}$ (B.19)

As F_n = P_n(z) (Eq. (8)) with z = cosx + ν(cosy − cosx) (Eq. (B.1)), we have from (B.19),

$\begin{matrix} ν (1 - ν) \frac{\partial^{2} F_{n}}{\partial ν^{2}} & + & (1 - 2 ν) \frac{\partial F_{n}}{∂ν} + \frac{1}{1 - ν} \frac{\partial^{2} F_{n}}{\partial x^{2}} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} \nu(1-\nu)\Dron{^2F_n}{\nu^2} \!\!&+&\!\! (1-2\nu)\Dron{F_n}{\nu} +\frac{1}{1-\nu}\Dron{^2F_n}{x^2} \nonumber\\ && +\frac{1}{\nu}\Dron{^2F_n}{y^2} +n(n+1)F_n=0. \label{eq.diffFn} \end{eqnarray}$ (B.20)

Using expressions (B.10, B.17), one then replaces F_n in (B.20) by

$F_{n} = \sum_{k} \sum_{l} (1 - ν)^{k} ν^{l} A_{k,l}^{(n)} (ν) e^{i kx} e^{i ly},$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} F_n = \sum_k\sum_l(1-\nu)^k\nu^lA^{(n)}_{k,l}(\nu)\e^{\i kx}\e^{\i ly} , \end{equation}$ (B.21)

which leads to

$\sum_{k} \sum_{l} (1 - ν)^{k} ν^{l} 𝒜_{k,l}^{(n)} (ν) e^{i kx} e^{i ly} = 0,$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} \sum_k\sum_l(1-\nu)^k\nu^l\cA^{(n)}_{k,l}(\nu)\e^{\i kx}\e^{\i ly}=0 , \label{eq.diffFnii} \end{equation}$ (B.22)

where

$\begin{matrix} 𝒜_{k,l}^{(n)} (ν) & = & ν (1 - ν) \frac{d^{2} A_{k,l}^{(n)}}{d ν^{2}} + [1 + 2 l - 2 (k + l + 1) ν] \frac{d A_{k,l}^{(n)}}{d ν} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} \cA^{(n)}_{k,l}(\nu) &=& \nu(1-\nu)\frac{{\rm d}^2A^{(n)}_{k,l}}{{\rm d}\nu^2} +\big[1+2l-2(k+l+1)\nu\big]\frac{{\rm d}A^{(n)}_{k,l}}{{\rm d}\nu} \nonumber\\ &&+\,(n-k-l)(n+k+l+1)A^{(n)}_{k,l}. \label{eq.diffAn} \end{eqnarray}$ (B.23)

The equality (B.22) must be satisfied for all x and y, thus all the coefficients $Mathematical equation: \appendix \setcounter{section}{2} $\cA^{(n)}_{k,l}(\nu)$$ (B.23) are equal to 0. The solutions of (B.23) are

$A_{k,l}^{(n)} (ν) = K_{k,l}^{(n)} F (k + l - n,k + l + n + 1; 1 + 2 l; ν)$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} A^{(n)}_{k,l}(\nu) = K^{(n)}_{k,l}F(k+l-n,k+l+n+1;1+2l;\nu) \label{eq.solAn} \end{equation}$ (B.24)

or equivalently

$\begin{matrix} A_{q - s,n - q - s}^{(n)} (ν) & = & K_{q - s,n - q - s}^{(n)} F (- 2 s, 2 n - 2 s + 1; \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} A^{(n)}_{q-s,n-q-s}(\nu) &=& K^{(n)}_{q-s,n-q-s}F(-2s,2n-2s+1; \nonumber\\ && 2n-2q-2s+1;\nu) , \label{eq.solAnii} \end{eqnarray}$ (B.25)

where $Mathematical equation: \appendix \setcounter{section}{2} $K^{(n)}_{k,l}$$ is an unknown constant and F an hypergeometric function (e.g. Whittaker & Watson 1927). Let us assume that the quantity 2n − 2q − 2s + 1 is positive. If it is not the case, one can make the change of variable (s,q) → (n − s′,n − q′). Then 2n − 2q′ − 2s′ + 1 is positive and since F_n is real, $Mathematical equation: \appendix \setcounter{section}{2} $\cQ^{(n)}_{s,q}(\nu)=\cQ^{(n)}_{s',q'}(\nu)$$ (B.10). From (B.25), one thus has

$\begin{matrix} A_{q - s,n - q - s}^{(n)} (ν) & = & K_{q - s,n - q - s}^{(n)} \frac{(2 s)! (2 n - 2 q - 2 s)!}{(2 n - 2 s)!} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} A^{(n)}_{q-s,n-q-s}(\nu) &=& K^{(n)}_{q-s,n-q-s} \frac{(2s)!(2n-2q-2s)!}{(2n-2s)!} \nonumber\\ &&\times\sum_{k=0}^{2s}(-1)^k \frac{(2n-2s+k)!}{(2s\! -\! k)!(2n\! -\! 2q\! -\! 2s\! +\! k)!}\frac{\nu^k}{k!}\cdot \label{eq.An} \end{eqnarray}$ (B.26)

Tisserand (1885) needs then some lengthy computations to determine $Mathematical equation: \appendix \setcounter{section}{2} $K^{(n)}_{k,l}$$ from the coefficient of νⁿ in F_n. Here, we use the expression (B.12) with (B.17) and (B.26). With μ = 1 − ν, we get

$\begin{matrix} K_{q - s,n - q - s}^{(n)} & (- 1)^{q - s} \frac{(2 n - 2 q - 2 s)! (2 n)!}{(2 n - 2 s)! (2 n - 2 q)!} = (- 1)^{q - s} \frac{(2 n)!}{2^{2 n} n!} \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} K^{(n)}_{q-s,n-q-s}&&\!\!\!\!\!\! (-1)^{q-s}\frac{(2n-2q-2s)!(2n)!}{(2n-2s)!(2n-2q)!} =(-1)^{q-s}\frac{(2n)!}{2^{2n}n!} \nonumber\\ &&\times \sum_{l=\max(0,s-q)}^{\min(s,n-q)} \frac{1}{l!(l+q-s)!(s-l)!(n\!-\!q\!-\!l)!}\cdot \label{eq.Kn} \end{eqnarray}$ (B.27)

Calculating the coefficient of the term x^s in (1 + x)^n − q(1 + x)^q = (1 + x)ⁿ, one finds

$\begin{matrix} \sum_{l = \max (0,s - q)}^{\min (s,n - q)} & \frac{1}{l! (l + q - s)! (s - l)! (n - q - l)!} = \\ \frac{n!}{s! (n - s)! q! (n - q)!} \cdot \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} \sum_{l=\max(0,s-q)}^{\min(s,n-q)}&&\!\!\!\frac{1}{l!(l+q-s)!(s-l)!(n-q-l)!}= \nonumber\\ && \hspace*{3cm} \frac{n!}{s!(n-s)!q!(n-q)!}\cdot \end{eqnarray}$ (B.28)

Thus, from (B.27), we obtain

$K_{q - s,n - q - s}^{(n)} = 12^{2 n} \frac{(2 n - 2 s)! (2 n - 2 q)!}{s! (n - s)! q! (n - q)! (2 n - 2 q - 2 s)!} \cdot$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} K^{(n)}_{q-s,n-q-s} = \frac{1}{2^{2n}}\frac{(2n-2s)!(2n-2q)!}{s!(n-s)!q!(n-q)!(2n-2q-2s)!}\cdot \label{eq.Knqs} \end{equation}$ (B.29)

Finally, the most general case, where inclinations are defined with respect to a fixed reference plane, can be as well derived from (B.17) and (B.29). We have

$F_{n} = \sum_{s = 0}^{n} \sum_{q = 0}^{n} \begin{matrix} ˜ \\ 𝒬_{s,q}^{(n)} \end{matrix} (μ_{*}, ν_{*}) e^{i (n - 2 s) u} e^{i (n - 2 q) u^{'}},$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} F_n = \sum_{s=0}^n\sum_{q=0}^n\tilde{\cQ}^{(n)}_{s,q}(\tmu,\tnu) \e^{\i(n-2s)u}\e^{\i(n-2q)u'} , \end{equation}$ (B.30)

where

$\begin{matrix} ˜ \\ 𝒬_{s,q}^{(n)} \end{matrix} (μ_{*}, ν_{*}) = {\begin{matrix} μ_{*}^{q - s} ν_{*}^{n - q - s} A_{q - s,n - q - s}^{(n)} (| ν_{*} |) & if s + q \leq n, \\ μ \begin{matrix} s - q \\ * \end{matrix} ν \begin{matrix} s + q - n \\ * \end{matrix} A_{s - q,s + q - n}^{(n)} (| ν_{*} |) & else, \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{equation} \tilde{\cQ}^{(n)}_{s,q}(\tmu,\tnu) = \left\{ \begin{array}{ll} \tmu^{q-s}\,\tnu^{n-q-s}A^{(n)}_{q-s,n-q-s}(\abs{\tnu}) & \textrm{if } s+q\leq n , \\[1mm] \bmus^{\ s-q}\,\bnus^{\ s+q-n}A^{(n)}_{s-q,s+q-n}(\abs{\tnu}) & \textrm{else}, \end{array} \right. \label{eq.Qtilde} \end{equation}$ (B.31)

and

$\begin{matrix} A_{q - s,n - q - s}^{(n)} (x) & = & \frac{1}{2^{2 n}} \frac{(2 s)! (2 n - 2 q)!}{s! (n - s)! q! (n - q)!} \\ \times \sum_{k = 0}^{2 s} (- 1)^{k} \frac{(2 n - 2 s + k)!}{(2 s - k)! (2 n - 2 q - 2 s + k)!} \frac{x^{k}}{k!} \cdot \end{matrix}$ $Mathematical equation: \appendix \setcounter{section}{2} \begin{eqnarray} A^{(n)}_{q-s,n-q-s}(x)\! &=&\! \frac{1}{2^{2n}} \frac{(2s)!(2n-2q)!}{s!(n-s)!q!(n-q)!} \nonumber \\ &&\times \sum_{k=0}^{2s}(-1)^k\frac{(2n-2s+k)!}{(2s\!-\!k)!(2n\!-\!2q\!-\!2s\!+\!k)!} \frac{x^k}{k!}\cdot \end{eqnarray}$ (B.32)

In (B.31), we used the fact that $Mathematical equation: \appendix \setcounter{section}{2} $\tilde{\cQ}^{(n)}_{s,q}$$ is the complex conjugate of $Mathematical equation: \appendix \setcounter{section}{2} $\tilde{\cQ}^{(n)}_{n-s,n-q}$$ since F_n is real.

References

Abu-El-Ata, N., & Chapront, J. 1975, A&A, 38, 57 [NASA ADS] [Google Scholar]
Aksenov, E. P. 1986, SvA, 30, 221 [NASA ADS] [Google Scholar]
Blaes, O., Lee, M. H., & Socrates, A. 2002, ApJ, 578, 775 [NASA ADS] [CrossRef] [Google Scholar]
Brumberg, V. A. 1967, Bull. Inst. Theor. Astron., XI, 73 [Google Scholar]
Brumberg, V. A., Evdokimova, L. S., & Kochina, N. G. 1971, Celest. Mech., 3, 197 [Google Scholar]
Farago, F., & Laskar, J. 2010, MNRAS, 401, 1189 [NASA ADS] [CrossRef] [Google Scholar]
Ford, E. B., Kozinsky, B., & Rasio, F. A. 2000, ApJ, 535, 385 [NASA ADS] [CrossRef] [Google Scholar]
Gastineau, M., & Laskar, J. 2009, TRIP 1.0, TRIP Reference manual, IMCCE, Paris Observatory, http://www.imcce.fr/Equipes/ASD/trip/trip.html [Google Scholar]
Giacaglia, G. E. O. 1974, Celest. Mech., 9, 239 [NASA ADS] [CrossRef] [Google Scholar]
Hansen, P. 1855, Abhandl. d. K. S. Ges. d. Wissensch, IV, 182 [Google Scholar]
Hill, G. W. 1875, The Analyst, II, 161 [CrossRef] [Google Scholar]
Hughes, S. 1981, Celest. Mech., 25, 101 [NASA ADS] [CrossRef] [Google Scholar]
Kaula, W. M. 1962, AJ, 67, 300 [NASA ADS] [CrossRef] [Google Scholar]
Kozai, Y. 1959, AJ, 64, 367 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
Krymolowski, Y., & Mazeh, T. 1999, MNRAS, 304, 720 [NASA ADS] [CrossRef] [Google Scholar]
Laplace, P. 1785, in Oeuvres Complètes, T11 (Paris: Gauthier-Villars) [Google Scholar]
Laskar, J. 1990, in Les Méthodes Modernes de la Mécanique Céleste, ed. C. F. D. Benest [Google Scholar]
Laskar, J. 2005, Celest. Mech. Dynam. Astron., 91, 351 [NASA ADS] [CrossRef] [Google Scholar]
Laskar, J., & Robutel, P. 1995, Celest Mech. Dynam. Astron., 62, 193 [Google Scholar]
Lee, M. H., & Peale, S. J. 2003, ApJ, 592, 1201 [NASA ADS] [CrossRef] [Google Scholar]
Lidov, M. L., & Ziglin, S. L. 1976, Celest. Mech., 13, 471 [NASA ADS] [CrossRef] [Google Scholar]
Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [NASA ADS] [CrossRef] [Google Scholar]
Migaszewski, C., & Goździewski, K. 2008, MNRAS, 388, 789 [NASA ADS] [CrossRef] [Google Scholar]
Tisserand, F. 1885, Annales de l’Observatoire de Paris, 18, 1 [Google Scholar]
Tisserand, F. 1899, Traité de Mécanique Céleste (Paris: Gauthier-Villars) [Google Scholar]
Whittaker, E., & Watson, G. 1927, A course of modern analysis (Cambridge University Press) [Google Scholar]

All Tables

Table 1

Tisserand functions for the planar case (Eq. (8)).

In the text

Table 2

Tisserand functions for the spatial case.

In the text

Table 3

Tisserand functions for the spatial case in a fixed reference frame.

In the text

Table A.1

$Mathematical equation: $\tilde X_0^{-n,m}$$ is the polynomial part of the Hansen coefficients $Mathematical equation: $X_0^{-n,m}$$ .

In the text

Table A.2

Hansen coefficients $Mathematical equation: $X_0^{n,m}$$ for (0 ≤ n ≤ 10,0 ≤ m ≤ n).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.

[R1] Abu-El-Ata, N., & Chapront, J. 1975, A&A, 38, 57 [NASA ADS] [Google Scholar]

[R2] Aksenov, E. P. 1986, SvA, 30, 221 [NASA ADS] [Google Scholar]

[R3] Blaes, O., Lee, M. H., & Socrates, A. 2002, ApJ, 578, 775 [NASA ADS] [CrossRef] [Google Scholar]

[R4] Brumberg, V. A. 1967, Bull. Inst. Theor. Astron., XI, 73 [Google Scholar]

[R5] Brumberg, V. A., Evdokimova, L. S., & Kochina, N. G. 1971, Celest. Mech., 3, 197 [Google Scholar]

[R6] Farago, F., & Laskar, J. 2010, MNRAS, 401, 1189 [NASA ADS] [CrossRef] [Google Scholar]

[R7] Ford, E. B., Kozinsky, B., & Rasio, F. A. 2000, ApJ, 535, 385 [NASA ADS] [CrossRef] [Google Scholar]

[R8] Gastineau, M., & Laskar, J. 2009, TRIP 1.0, TRIP Reference manual, IMCCE, Paris Observatory, http://www.imcce.fr/Equipes/ASD/trip/trip.html [Google Scholar]

[R9] Giacaglia, G. E. O. 1974, Celest. Mech., 9, 239 [NASA ADS] [CrossRef] [Google Scholar]

[R10] Hansen, P. 1855, Abhandl. d. K. S. Ges. d. Wissensch, IV, 182 [Google Scholar]

[R11] Hill, G. W. 1875, The Analyst, II, 161 [CrossRef] [Google Scholar]

[R12] Hughes, S. 1981, Celest. Mech., 25, 101 [NASA ADS] [CrossRef] [Google Scholar]

[R13] Kaula, W. M. 1962, AJ, 67, 300 [NASA ADS] [CrossRef] [Google Scholar]

[R14] Kozai, Y. 1959, AJ, 64, 367 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]

[R15] Krymolowski, Y., & Mazeh, T. 1999, MNRAS, 304, 720 [NASA ADS] [CrossRef] [Google Scholar]

[R16] Laplace, P. 1785, in Oeuvres Complètes, T11 (Paris: Gauthier-Villars) [Google Scholar]

[R17] Laskar, J. 1990, in Les Méthodes Modernes de la Mécanique Céleste, ed. C. F. D. Benest [Google Scholar]

[R18] Laskar, J. 2005, Celest. Mech. Dynam. Astron., 91, 351 [NASA ADS] [CrossRef] [Google Scholar]

[R19] Laskar, J., & Robutel, P. 1995, Celest Mech. Dynam. Astron., 62, 193 [Google Scholar]

[R20] Lee, M. H., & Peale, S. J. 2003, ApJ, 592, 1201 [NASA ADS] [CrossRef] [Google Scholar]

[R21] Lidov, M. L., & Ziglin, S. L. 1976, Celest. Mech., 13, 471 [NASA ADS] [CrossRef] [Google Scholar]

[R22] Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [NASA ADS] [CrossRef] [Google Scholar]

[R23] Migaszewski, C., & Goździewski, K. 2008, MNRAS, 388, 789 [NASA ADS] [CrossRef] [Google Scholar]

[R24] Tisserand, F. 1885, Annales de l’Observatoire de Paris, 18, 1 [Google Scholar]

[R25] Tisserand, F. 1899, Traité de Mécanique Céleste (Paris: Gauthier-Villars) [Google Scholar]

[R26] Whittaker, E., & Watson, G. 1927, A course of modern analysis (Cambridge University Press) [Google Scholar]