In summary, we developed a numerically nearly exact method for calculating the power factor of organic thermoelectric materials modeled by electron-phonon (Holstein) model via combining Kubo formula with finite-temperature time-dependent density matrix renormalization group (FT-TD-DMRG) in the grand canonical ensemble. We demonstrated the validity and computational feasibility for a broad parameter range.