A smooth time-step selection formula for the midpoint method is derived which minimize deviations in the Hamiltonian function along piecewise-linear phase space trajectories of autonomous Hamiltonian systems. The time-step formula is implemented in a second order predictor/corrector scheme and applied to Kepler's problem. The formula significantly improves energy conservation as well as the accuracy of the configuration space trajectory. Peak errors in position and momentum coordinates are not significantly reduced, but the time behavior of the errors is markedly more regular.
Shibberu, Yosi, "A Variable Time-Step Midpoint Scheme for Hamiltonian Systems" (1995). Mathematical Sciences Technical Reports (MSTR). 120.