# Tidal evolution of the Keplerian elements

## TL;DR

This paper derives corrected formulas for the rates of Keplerian orbital elements under tidal perturbations, addressing errors in previous literature and providing expressions applicable to systems with comparable masses.

## Contribution

It presents a re-derivation of the tidal evolution equations, correcting previous formulas and clarifying the contributions to orbital inclination changes for binary systems.

## Key findings

- Derived new expressions for da/dt, de/dt, di/dt including higher order terms.
- Identified and corrected errors in Kaula's classical formulas for systems with comparable masses.
- Highlighted the importance of angular momentum conservation in orbital inclination evolution.

## Abstract

We address the expressions for the rates of the Keplerian orbital elements within a two-body problem perturbed by the tides in both partners. The formulae for these rates have appeared in the literature in various forms, at times with errors. We reconsider, from scratch, the derivation of these rates and arrive at the Lagrange-type equations which, in some details, differ from the corresponding equations obtained previously by Kaula (1964). We also write down detailed expressions for $da/dt$, $de/dt$ and $di/dt$, to order $e^4$. They differ from Kaula's expressions which contain a redundant factor of $M/(M+M^{\prime}),$ with $M$ and $M^{\prime}$ being the masses of the primary and the secondary. As Kaula was interested in the Earth-Moon system, this redundant factor was close to unity and was unimportant in his developments. This factor, however, must be reinstated when Kaula's theory is applied to a binary composed of partners of comparable masses. We have found that, while it is legitimate to simply sum the primary's and secondary's inputs in $da/dt$ or $de/dt$, this is not the case for $di/dt$. So our expression for $di/dt$ differs from that of Kaula in two regards. First, the contribution due to the dissipation in the secondary averages out when the apsidal precession is uniform. Second, we have obtained an additional term which emerges owing to the conservation of the angular momentum: a change in the inclination of the orbit causes a change of the primary's plane of equator.