More recently, Meek and Levine have demonstrated that interpolating the derivative coupling pairwise between adiabatic states (and then averaging the T matrix over that time step) provides a better approximation for the case of two electronic states and allows for even larger time steps. This approach can be extended to the case of more than two electronic states if we recognized that the T matrix can be calculated directly from the logarithm of the overlap matrix 23