Tkatchenko-Scheffler does not work with variable cell optimization

See attached file: works nicely for calculation='relax', oscillates between two energies for calculation='vc-relax'. This is the typical signature of a mismatch between energies and forces. ts+vcrelax.in