Orbital evolution of Phaethon under different dynamical models (left column panels) and their differences (right panels). The evolutions of the semimajor axis, eccentricity, and inclination are presented in the left column panels, where the gray lines stand for the evolution without the Yarkovsky force, the blue lines represent the evolution with a linear approximation of the Yarkovsky force, and the red lines are for the evolution with the Yarkovsky force produced from DeepONet.