To date, there are several implementations of FEH calculations using time-dependent density functional theory (TDDFT) in the linear-response limit. For instance, in Ref. 9 a diabatization scheme with a transformation matrix that reproduces local transition dipoles has been proposed. This approach was extended in Ref. 10 to include charge transfer states via minimization of the difference between the diabatized and a reference transition charge matrix.