The pre-computed precession-tilt solutions PT-ZB18a(Ed,Td) span the past 100 Myr. Available solutions include all combinations of dynamical ellipticity values between 0.9950 and 1.0050 in increments of 0.0010 and tidal dissipation values between 0.000 and 1.2000 in increments of 0.1.
Source
All astronomical solutions by Zeebe can be found on http://www.soest.hawaii.edu/oceanography/faculty/zeebe_files/Astro.html.
They can be loaded into R quickly, using get_solution()
.
References
Zeebe, R. E. and Lourens, L. J. (2022). Geologically constrained astronomical solutions for the Cenozoic era. Earth and Planetary Science Letters. doi:10.1016/j.epsl.2022.117595