--- title: Democratic Heliocentric Coordinates url: https://www.emergentmind.com/topics/democratic-heliocentric-coordinates-dhc type: topic --- # Democratic Heliocentric Coordinates Democratic heliocentric coordinates (DHC) are a canonical coordinate system employed in symplectic N-body integrations, most notably within Wisdom–Holman (WH) integrators, to model the long-term evolution of planetary systems. DHC are constructed to treat all planets equivalently with respect to the heliocentric frame, and enable operator splitting of the Hamiltonian into analytically tractable Keplerian, interplanetary, and so-called “jump” (solar) terms. Recent numerical experiments demonstrate that DHC introduce an eccentricity-dependent artificial precession in planetary orbits, significantly affecting the rates of secular instabilities, such as those involving Mercury in the Solar System, unless the integration timestep is reduced to values much smaller than are typically used [2601.08019]. ## 1. Canonical Variables and DHC Transformations The DHC construction proceeds from standard barycentric N-body variables—positions $r_i$, momenta $p_i$, masses $m_0$ (Sun), $m_1,\ldots,m_N$ (planets), and total mass $M = \sum_{i=0}^N m_i$. The canonical change of variables consists of: - Center-of-mass coordinate and momentum: - $Q_0 = \frac{1}{M} \sum_{i=0}^N m_i r_i$ - $P_0 = \sum_{i=0}^N p_i$ - Heliocentric “relative” coordinates and their conjugate momenta for $i=1,\ldots,N$: - $Q_i = r_i - r_0$ - $P_i = p_i - \frac{m_i}{M} \sum_{j=0}^N p_j$ These transformed variables are canonical, i.e., $\{Q_\alpha, P_\beta\} = \delta_{\alpha\beta}$, and their inverses express barycentric positions and momenta in terms of $Q_i, P_i$, and the center of mass variables. The total momentum $P_0$ is conserved, so the barycentric motion trivially decouples from the internal planetary evolution [2601.08019]. ## 2. Wisdom–Holman Hamiltonian Splitting in DHC In DHC, the full N-body Hamiltonian is split exactly into three terms: - $H_K$: Kepler Hamiltonian for each planet’s heliocentric motion, $$ H_K = \sum_{i=1}^N \left( \frac{\|P_i\|^2}{2m_i} - G \frac{m_0 m_i}{\|Q_i\|} \right) $$ - $H_I^{\rm pp}$: Interplanetary potential, $$ H_I^{\rm pp} = - G \sum_{1 \leq i < j \leq N} \frac{m_i m_j}{\|Q_i - Q_j\|} $$ - $H_I^{\rm Sun}$: Solar “jump” term, $$ H_I^{\rm Sun} = \frac{1}{2m_0} \left\| \sum_{i=1}^N P_i \right\|^2 $$ The standard second-order WH map advances the system over timestep $\Delta t$ via a sequence ("kick–kick–drift–kick–kick"): 1. **Kick**: update $Q_i$ by $H_I^{\rm Sun}$, 2. **Kick**: update $P_i$ by $H_I^{\rm pp}$, 3. **Drift**: exactly solve each heliocentric two-body problem for $\Delta t$, 4. **Kick**: repeat step 2, 5. **Kick**: repeat step 1. Each component of the Hamiltonian is solved exactly within its substep; only their noncommutativity induces integration errors [2601.08019]. ## 3. Eccentricity-Dependent Artificial Precession in DHC Due to the Baker–Campbell–Hausdorff expansion, second-order WH splitting in DHC generates leading errors that couple $H_K$, $H_I^{\rm pp}$, and notably $H_I^{\rm Sun}$. This manifests as an $\mathcal{O}(\Delta t^2)$ secular term in each heliocentric Kepler subsystem, generating an artificial apsidal precession. Numerical results show: - The artificial precession rate $\Delta g$ is proportional to $\Delta t^2$ and increases rapidly with orbital eccentricity $e$. - This precession exceeds the physical general-relativistic precession $g_{\rm GR} \approx 43''$ cy$^{-1}$ when $\Delta t \gtrsim T_f/4$, with $$ T_f = \frac{2\pi}{n_1} \frac{(1-e_1)^2}{\sqrt{1-e_1^2}} $$ where $n_1$ is the planet’s mean motion. - For $e_1 = 0.6$–$0.8$, the DHC-induced precession grows rapidly at $\Delta t \sim T_f$, undermining accuracy at typical symplectic timesteps; Jacobi-coordinate splits remain accurate for $\Delta t$ up to $10$ days or greater, even at high $e_1$ [2601.08019]. ## 4. Numerical Integration Outcomes in the Solar System Long-term ensemble integrations serve to illustrate the practical consequences of DHC’s artificial precession. Key findings from $5$ Gyr Solar System runs are: | Integration Setting | Δt | Instability Rate (%) | Notes | |--------------------------|----------------|----------------------|-----------------------------------------| | DHC | ≈ 6 days | 0 | Mercury instabilities artificially suppressed | | DHC | ≈ 0.6 days | ≈ 1 | Agrees with Jacobi results | | Jacobi (with GR) | 6 days | 100 | Instabilities present (converged) | | DHC (restart, high $e$) | 6 days | 0 | Artificially stabilized | | DHC (restart, high $e$) | 0.6 days | 100 | Converged (instabilities appear) | These results demonstrate that with standard DHC timesteps (several days), secular increases in eccentricity (e.g., Mercury’s $e_1$ surge via $g_1$–$g_5$ resonance) are suppressed by the numerical artifact, leading to an unrealistically low instability rate. Reducing the timestep by an order of magnitude restores physical instability rates and agreement with Jacobi splits [2601.08019]. ## 5. Comparison with Jacobi Coordinate Wisdom–Holman Integrators Jacobi coordinates employ a hierarchical splitting wherein planet $i$ orbits the cumulative mass $M_{i-1} = \sum_{k=0}^{i-1} m_k$. The Kepler and interaction Hamiltonians are: - $H_K^{\rm Jacobi} = \sum_{i=1}^N \left( \frac{\|P_i\|^2}{2\mu_i} - G \frac{M_{i-1} m_i}{\|Q_i\|} \right)$, with reduced mass $\mu_i = m_i M_{i-1}/(m_i + M_{i-1})$, - $H_I^{\rm Jacobi} = -G \sum_{i