This Letter introduces the so-called Quasi Time-Reversible scheme based on Grassmann extrapolation (QTR G-Ext) of density matrices for an accurate calculation of initial guesses in Born-Oppenheimer Molecular Dynamics (BOMD) simulations. The method shows excellent results on four large molecular systems that are representative of real-life production applications, ranging from 21 to 94 atoms simulated with Kohn-Sham (KS) density functional theory surrounded with a classical environment with 6k to 16k atoms. Namely, it clearly reduces the number of self-consistent field iterations while at the same time achieving energy-conserving simulations, resulting in a considerable speed-up of BOMD simulations even when tight convergence of the KS equations is required.
Pes, F.; Polack, É.; Mazzeo, P.; Dusson, G.; Stamm, B. & Lipparini, F.
J. Phys. Chem. Lett. 14, 43, 9720-9726 (2023) https://pubs.acs.org/doi/10.1021/acs.jpclett.3c02098
