“…In the initial stages of geometry or reaction path optimization involving crude searches, the less accurate but more rapid method is sufficient, such as the L€ owdin-Inverse-L€ owdin (LIL) density matrix extrapolation, [11,34,40] relying only overlap matrix calculations at any new nuclear geometry, and after a reasonably reduced neighborhood of the desired conformation is found in the configuration space by the crude method, the fullfledged ADMA approach is used for a more accurate determination. [47,48,50,51,53,54] These two methods, ADMA and LIL, are, indeed, fully compatible, as they both use, as input as well as output, density matrices of the same dimension and involving a common AO basis set.…”