Is the dynamics of Jupiter family comets amenable to Monte Carlo modelling?
Authors: Fouchard M.1; Froeschlé C.1; Valsecchi G.B.2
Source: Monthly Notices of the Royal Astronomical Society, Volume 344, Number 4, October 2003 , pp. 1283-1295(13)
Publisher: Blackwell Publishing
- Free Content
- New Content
- Subscribed Content
- Free Trial Content
Abstract:
In order to explore the range of applicability of Monte Carlo modelling, we reconsider in depth the Monte Carlo simulations of Froeschlé & Rickman aimed at the study of the orbital evolution of Jupiter family comets. We first analyse the computation of perturbation distributions, and find that to get statistically significant results, the sample size of these distributions must be larger than several times 105. We then show the necessity of introducing a procedure, referred to as the Spurious Perturbations Removal (SPR) procedure, to clean the perturbation distribution of any spurious contribution due to comets undergoing a close encounter with Jupiter at the beginning or/and the end of the perturbation computations. The partition of element space into 30 boxes in the (Q, q) plane ( 4
Q 13 au and 0 q 7 au ) does not lead to dynamically homogeneous boxes, as is required for Monte Carlo simulations; these inhomogeneities are due to various reasons. Different ways to make such a partition are discussed. The strong influence of mean motion resonances on the dynamics of Jupiter family comets implies that resonances must be taken into account in any Monte Carlo simulation.
Keywords: methods: numerical; methods: statistical; celestial mechanics; comets: general
Document Type: Research article
DOI: 10.1046/j.1365-8711.2003.06904.x
Affiliations: 1: Observatoire de la Côte d'Azure, CNRS/UMR 6529, Bv. de l'Observatoire, B.P. 4229, F-06304 Nice cedex 4, France 2: IASF-CNR, via Fosso del Cavaliere, I-00133 Roma, Italy
- Free Content
- New Content
- Subscribed Content
- Free Trial Content
Click here for Page Help