OSSOS: The eccentricity and inclination distributions of the stable Neptunian Trojans

The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely em:laced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the libration amplitude of the stable Neptunian Trojan population can be well modeled as a Rayleigh distribution with a libration amplitude width σ_Aϕ of 15°. When taken as a whole, the Neptune Trojan population can be acceptably modeled with a Rayleigh eccentricity distribution of width σ_e of 0.045 and a typical sin(i) × Gaussian inclination distribution with a width σ_i of 14±2^∘; however, these distributions are only marginally acceptable. This is likely because, even after accounting for survey detection biases, the known large (H_r<8) and small (H_r≥8) Neptune Trojans appear to have markedly different eccentricities and inclinations. We propose that like the classical Kuiper belt, the stable intrinsic Neptunian Trojan population have dynamically ‘hot’ and dynamically ‘cold’ components to its eccentricity/inclination distribution, with σ_e−cold∼0.02/σ_i−cold∼6° and σ_e−hot∼0.05/σ_i−hot∼18°. In this scenario, the ‘cold’ L4 Neptunian Trojan population lacks the H_r≥8.0 members and has 13₋₆⁺¹¹ ‘cold’ Trojans with H_r<8.0. On the other hand, the ‘hot’ L4 Neptunian Trojan population has 136₋₇₅⁺⁸⁴ Trojans with H_r<10 — a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range.