 Poynting–Robertson drag invariant

The orbit-averaged Poynting–Robertson drag equations preserve $C=a(1-e^2)e^{-4/5}$ for $0<e<1$. This relates shrinking semi-major axis to decreasing orbital eccentricity; a circular orbit is treated by its separate limiting solution.