### Abstract

Monte-Carlo codes generally treat planestesimal-planet encounters using the two-body scattering approximation, which can be inaccurate when relative velocities are low; however, Monte-Carlo codes using the two-body approximation frequently produce results consistent with more accurate codes using numerical integration. To better understand why this breakdown occurs at low velocities, and to test a hypothesis from Greenberg et al. (1988, Icarus 75, 1-29) that may explain the unexpected accuracy of Monte-Carlo codes, we numerically integrate test body trajectories using a unique set of orbital elements defined by the geometry of the two-body approximation. This new coordinate system is ideal for examining the effects of distant planetary perturbations on particle trajectories all the way to encounter with the planet. Our results show that the failure of the two-body approximation is caused by distant planetary perturbations modifying the approach geometry of the test bodies; behavior at encounter follows two-body scattering even at very low relative velocities. By testing particle swarms encountering a planet, we found that some test bodies, whose approach orbits were shifted by distant planetary perturbations, were then replaced by similarly shifted nearby test bodies. The "particle replacement" mechanism explains why Monte-Carlo codes frequently yield outcome results comparable to numerical integration results. Moreover, we found that the relative velocity of a test body at encounter is not the critical parameter in determining the "breakdown" of two-body scattering outcome statistics; instead, we found that the semimajor axis of the test body relative to the size of the planet's Hill sphere (or the synodic period of the test body when mass is included) is much more diagnostic. Thus, our results verify that Monte-Carlo models can yield statistically accurate results, even if individual particles do not behave as assumed in those codes.

Original language | English (US) |
---|---|

Pages (from-to) | 288-301 |

Number of pages | 14 |

Journal | Icarus |

Volume | 125 |

Issue number | 2 |

State | Published - Feb 1997 |

### Fingerprint

### ASJC Scopus subject areas

- Space and Planetary Science
- Astronomy and Astrophysics

### Cite this

*Icarus*,

*125*(2), 288-301.

**Mapping the effects of distant perturbations on particle-planet interactions.** / Bottke, William F.; Greenberg, Richard J.; Carusi, Andrea; Valsecchi, Giovanni B.

Research output: Contribution to journal › Article

*Icarus*, vol. 125, no. 2, pp. 288-301.

}

TY - JOUR

T1 - Mapping the effects of distant perturbations on particle-planet interactions

AU - Bottke, William F.

AU - Greenberg, Richard J.

AU - Carusi, Andrea

AU - Valsecchi, Giovanni B.

PY - 1997/2

Y1 - 1997/2

N2 - Monte-Carlo codes generally treat planestesimal-planet encounters using the two-body scattering approximation, which can be inaccurate when relative velocities are low; however, Monte-Carlo codes using the two-body approximation frequently produce results consistent with more accurate codes using numerical integration. To better understand why this breakdown occurs at low velocities, and to test a hypothesis from Greenberg et al. (1988, Icarus 75, 1-29) that may explain the unexpected accuracy of Monte-Carlo codes, we numerically integrate test body trajectories using a unique set of orbital elements defined by the geometry of the two-body approximation. This new coordinate system is ideal for examining the effects of distant planetary perturbations on particle trajectories all the way to encounter with the planet. Our results show that the failure of the two-body approximation is caused by distant planetary perturbations modifying the approach geometry of the test bodies; behavior at encounter follows two-body scattering even at very low relative velocities. By testing particle swarms encountering a planet, we found that some test bodies, whose approach orbits were shifted by distant planetary perturbations, were then replaced by similarly shifted nearby test bodies. The "particle replacement" mechanism explains why Monte-Carlo codes frequently yield outcome results comparable to numerical integration results. Moreover, we found that the relative velocity of a test body at encounter is not the critical parameter in determining the "breakdown" of two-body scattering outcome statistics; instead, we found that the semimajor axis of the test body relative to the size of the planet's Hill sphere (or the synodic period of the test body when mass is included) is much more diagnostic. Thus, our results verify that Monte-Carlo models can yield statistically accurate results, even if individual particles do not behave as assumed in those codes.

AB - Monte-Carlo codes generally treat planestesimal-planet encounters using the two-body scattering approximation, which can be inaccurate when relative velocities are low; however, Monte-Carlo codes using the two-body approximation frequently produce results consistent with more accurate codes using numerical integration. To better understand why this breakdown occurs at low velocities, and to test a hypothesis from Greenberg et al. (1988, Icarus 75, 1-29) that may explain the unexpected accuracy of Monte-Carlo codes, we numerically integrate test body trajectories using a unique set of orbital elements defined by the geometry of the two-body approximation. This new coordinate system is ideal for examining the effects of distant planetary perturbations on particle trajectories all the way to encounter with the planet. Our results show that the failure of the two-body approximation is caused by distant planetary perturbations modifying the approach geometry of the test bodies; behavior at encounter follows two-body scattering even at very low relative velocities. By testing particle swarms encountering a planet, we found that some test bodies, whose approach orbits were shifted by distant planetary perturbations, were then replaced by similarly shifted nearby test bodies. The "particle replacement" mechanism explains why Monte-Carlo codes frequently yield outcome results comparable to numerical integration results. Moreover, we found that the relative velocity of a test body at encounter is not the critical parameter in determining the "breakdown" of two-body scattering outcome statistics; instead, we found that the semimajor axis of the test body relative to the size of the planet's Hill sphere (or the synodic period of the test body when mass is included) is much more diagnostic. Thus, our results verify that Monte-Carlo models can yield statistically accurate results, even if individual particles do not behave as assumed in those codes.

UR - http://www.scopus.com/inward/record.url?scp=0031071165&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0031071165&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0031071165

VL - 125

SP - 288

EP - 301

JO - Icarus

JF - Icarus

SN - 0019-1035

IS - 2

ER -