David M. Hernandez

Wisdom-Holman (WH) methods are algorithms used as a basis for a wide range of codes used to solve problems in solar system and planetary dynamics. The problems range from the growth and migration of planets to the stability of the solar system. In many cases, these codes work with Democratic Heliocentric Coordinates (DHC) which offer some advantages. However, it has been noted these coordinates affect the dynamics of solar system bodies in simulations, in particular Mercury's, and introduce artificial precession which affects solar system stability. In this work, we analytically derive the two-body artificial precession induced by DHC. We show the effect is small for solar system bodies, but the artificial effect on Jupiter is $242$ times larger than on Mercury. In a two-body Mercury-Sun system with general relativity (GR), artificial precession is negligible compared to GR precession, even with extreme timesteps that amplify the numerical effects. In a full solar system model, numerical effects are amplified further.