Orbital stability of Earth Trojans

The only discovery of Earth Trojan 2010 TK7 and the subsequent launch of OSIRIS-REx have motived us to investigate the stability around the triangular Lagrange points of the Earth, L-4 and L-5. In this paper we present detailed dynamical maps on the (a(0); i(0)) plane with the spectral number (SN) indicating the stability. Two main stability regions, separated by a chaotic region arising from the nu(3) and nu(4) secular resonances, are found at low (i(0) <= 15 degrees) and moderate (24 <= i(0) <= 37 degrees) inclinations, respectively. The most stable orbits reside below i(0) = 10 degrees and they can survive the age of the solar system. The nodal secular resonance nu(13) could vary the inclinations from 0 degrees to similar to 10 degrees according to their initial values, while nu(14) could pump up the inclinations to similar to 20 degrees and upwards. The fine structures in the dynamical maps are related to higher degree secular resonances, of which different types dominate different areas. The dynamical behaviour of the tadpole and horseshoe orbits, reflected in their secular precession, show great differences in the frequency space. The secular resonances involving the tadpole orbits are more sensitive to the frequency drift of the inner planets, thus the instabilities could sweep across the phase space, leading to the clearance of tadpole orbits. We are more likely to find terrestrial companions on horseshoe orbits. The Yarkovsky effect could destabilize Earth Trojans in varying degrees. We numerically obtain the formula describing the stabilities affected by the Yarkovsky effect and find the asymmetry between the prograde and retrograde rotating Earth Trojans. The existence of small primordial Earth Trojans that avoid being detected but survive the Yarkovsky effect for 4.5 Gyr is substantially ruled out.