The charge-carrier dynamics model presented below is an adaptation of that introduced in Ref. and adopted also in Refs. and . Within this model, charge-carrier's motion is usually considered as a series of hopping events between localized eigenstates ({|𝑖⟩}) which are eigenstates of the zeroth-order Hamiltonian given in Eq. 3. It is also assumed that the process follows a first-order reaction kinetics, and the coupling of the electronic degrees of freedom with the nuclear vibrations equilibrates the charge-carrier between successive hopping events. From that point of view, we consider a Holstein-type local electron-phonon coupling term (Eq. 4) as the first-order perturbation to the zeroth-order Hamiltonian given in Eq. 3.