Here S ab ′ ij ′ is a square matrix of size N o where the i-th orbital is replaced by the virtual orbital a for the geometry at time t and the j ′ -th orbital is replaced by the virtual orbital b ′ for the geometry at time t + dt.