High-order regularised symplectic integrator for collisional planetary systems

Research output: Contribution to journal › Journal article › Research › peer-review

Standard

High-order regularised symplectic integrator for collisional planetary systems. / Petit, Antoine C.; Laskar, Jacques; Boué, Gwenaël; Gastineau, Mickaël.

In: Astronomy and Astrophysics, Vol. 628, A32, 2019.

Research output: Contribution to journal › Journal article › Research › peer-review

Harvard

Petit, AC, Laskar, J, Boué, G & Gastineau, M 2019, 'High-order regularised symplectic integrator for collisional planetary systems', Astronomy and Astrophysics, vol. 628, A32. https://doi.org/10.1051/0004-6361/201935786

APA

Petit, A. C., Laskar, J., Boué, G., & Gastineau, M. (2019). High-order regularised symplectic integrator for collisional planetary systems. Astronomy and Astrophysics, 628, [A32]. https://doi.org/10.1051/0004-6361/201935786

Vancouver

Petit AC, Laskar J, Boué G, Gastineau M. High-order regularised symplectic integrator for collisional planetary systems. Astronomy and Astrophysics. 2019;628. A32. https://doi.org/10.1051/0004-6361/201935786

Author

Petit, Antoine C. ; Laskar, Jacques ; Boué, Gwenaël ; Gastineau, Mickaël. / High-order regularised symplectic integrator for collisional planetary systems. In: Astronomy and Astrophysics. 2019 ; Vol. 628.

Bibtex

@article{f96d502ec5814660868f2eee8eac2bce,

title = "High-order regularised symplectic integrator for collisional planetary systems",

abstract = "We present a new mixed variable symplectic (MVS) integrator for planetary systems that fully resolves close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators while also reducing the effective step size when two planets encounter. We used a high-order MVS scheme so that it was possible to integrate with large time-steps far away from close encounters. We show that this algorithm is able to resolve almost exact collisions (i.e. with a mutual separation of a fraction of the physical radius) while using the same time-step as in a weakly perturbed problem such as the solar system. We demonstrate the long-term behaviour in systems of six super-Earths that experience strong scattering for 50 kyr. We compare our algorithm to hybrid methods such as MERCURY and show that for an equivalent cost, we obtain better energy conservation.",

keywords = "Celestial mechanics, Methods: Numerical, Planets and satellites: Dynamical evolution and stability",

author = "Petit, {Antoine C.} and Jacques Laskar and Gwena{\"e}l Bou{\'e} and Micka{\"e}l Gastineau",

note = "Publisher Copyright: {\textcopyright} 2019 A. C. Petit et al.",

year = "2019",

doi = "10.1051/0004-6361/201935786",

language = "English",

volume = "628",

journal = "Astronomy & Astrophysics",

issn = "0004-6361",

publisher = "E D P Sciences",

}

RIS

TY - JOUR

T1 - High-order regularised symplectic integrator for collisional planetary systems

AU - Petit, Antoine C.

AU - Laskar, Jacques

AU - Boué, Gwenaël

AU - Gastineau, Mickaël

PY - 2019

Y1 - 2019

N2 - We present a new mixed variable symplectic (MVS) integrator for planetary systems that fully resolves close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators while also reducing the effective step size when two planets encounter. We used a high-order MVS scheme so that it was possible to integrate with large time-steps far away from close encounters. We show that this algorithm is able to resolve almost exact collisions (i.e. with a mutual separation of a fraction of the physical radius) while using the same time-step as in a weakly perturbed problem such as the solar system. We demonstrate the long-term behaviour in systems of six super-Earths that experience strong scattering for 50 kyr. We compare our algorithm to hybrid methods such as MERCURY and show that for an equivalent cost, we obtain better energy conservation.

AB - We present a new mixed variable symplectic (MVS) integrator for planetary systems that fully resolves close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators while also reducing the effective step size when two planets encounter. We used a high-order MVS scheme so that it was possible to integrate with large time-steps far away from close encounters. We show that this algorithm is able to resolve almost exact collisions (i.e. with a mutual separation of a fraction of the physical radius) while using the same time-step as in a weakly perturbed problem such as the solar system. We demonstrate the long-term behaviour in systems of six super-Earths that experience strong scattering for 50 kyr. We compare our algorithm to hybrid methods such as MERCURY and show that for an equivalent cost, we obtain better energy conservation.

KW - Celestial mechanics

KW - Methods: Numerical

KW - Planets and satellites: Dynamical evolution and stability

U2 - 10.1051/0004-6361/201935786

DO - 10.1051/0004-6361/201935786

M3 - Journal article

AN - SCOPUS:85070237838

VL - 628

JO - Astronomy & Astrophysics

JF - Astronomy & Astrophysics

SN - 0004-6361

M1 - A32

ER -

ID: 327138388

Globe Institute