Normalized to: Pérez-Hernández, J.
[1]
oai:arXiv.org:1904.09342 [pdf] - 1895594
On the dynamics of Comet 1P/Halley: Lyapunov and power spectra
Submitted: 2019-04-19
Using a purely Newtonian model for the Solar System, we investigate the
dynamics of comet 1P/Halley considering in particular the Lyapunov and power
spectra of its orbit, using the nominal initial conditions of JPL's Horizons
system. We carry out precise numerical integrations of the $(N+1)$-restricted
problem and the first variational equations, considering a time span of
$2\times10^5$~yr. The power spectra are dominated by a broadband component,
with peaks located at the current planetary frequencies, including
contributions from Jupiter, Venus, the Earth and Saturn, as well as the $1:6$
resonance among Halley and Jupiter and higher harmonics. From the average value
of the maximum Lyapunov exponent we estimate the Lyapunov time of the comet's
nominal orbit, obtaining $\tau_L \simeq 562$~yr; the remaining independent
Lyapunov exponents (not related by time-reversal symmetry) tend asymptotically
to zero as $t^{-1/2}$. Yet, our results do not display convergence of the
maximum Lyapunov exponent. We argue that the lack of convergence of the maximum
Lyapunov exponent is a signature of transient chaos which will lead to an
eventual ejection of the comet from the Solar System.