Past dynamical evolution of Phaethon’s spin axis direction s: left column of panels for the P1 solution (λ,β) = (319°,−39°), right column of panels for the P2 solution (λ,β) = (84°,−39°). The upper two panels show the spin evolution over the whole integrated interval of 1 Myr, the bottom two panels zoom at the past 50 kyr. Each time the nominal solution is shown by the black line, while the gray lines are for the 50 spin clones (see the text). Each of the four sections shows (i) the obliquity ε at the top (angle between s and normal to the osculating orbital plane); and (ii) the angle α between s and the direction to the Sun at pericenter.