On this page you can find interesting graphs, comparsions and features of Yarkovksy effects, which were not covered by previous graphs.
The effect of nonzero albedo and nonunity emissivity is clearly visible in these comparsion graphs da/dt. The yellow line is plotted for albedo = 0.2, infrared emissivity 0.9. The drift rates are approximately 50 % smaller if you compare them with drift rates for albedo = 0.0 and emissivity 0.9.
High eccentricity (in this simulation e = 0.9) can cause a very important increase of da/dt and also de/dt rates. These two graphs are plotted for diurnal variant of Yarkovsky force acting on regolith covered bodies.
We have observed an increase of semimajor axis for fragments smaller than 0.3 m and eccentricity 0.9, when only seasonal effect was working (diurnal effect was supressed by the orientation of fragment's spin axis). The y-axis of the plot is linear, not logarithmic as previously.
Drift of semimajor axis for wide range of heliocentric distances for small asteroid fragments with radius R = 1 m. Notice that diurnal variant of thermal effect is not maximal for the smallest initial semimajor axis, however changes of semimajor axes such close to the Sun are not linear for the integration time 0.1 Myr.