Robust implicit adaptive low rank time-stepping methods for matrix differential equations