In Eq. 10, we assume that 𝜆 𝑖𝑗 is the sum of a classical component (𝜆 𝑖𝑗 (𝐶) = (1 -𝑓 𝑄 )𝜆 𝑖𝑗 ) due to low-frequency nuclear modes and a quantum component (𝜆 𝑖𝑗 (𝑄) = 𝑓 𝑄 𝜆 𝑖𝑗 ) due to one effective high-frequency mode ℏ𝜔 ≈ 1450 cm -1 . 𝑓 𝑄 = 0.4 is the reorganization energy partition factor which is close to the value emerges from electronic structure calculations and has been utilized in a previous study. In addition, 𝑤 and 𝑤 ′ are the initial and final vibrational states along the high-frequency mode, 𝑃(𝑤) is the Boltzmann population of the initial vibrational state, and 𝐹𝐶 𝑖𝑗,𝑤𝑤 ′ is the Franck-Condon overlap integral between the vibronic states |𝑖, 𝑤⟩ and |𝑗, 𝑤 ′ ⟩. Employing these calculated hopping rates, the intra-chain charge-carrier mobility (𝜇) is calculated within low charge-carrier and weak electric field approximation within a standard master equation approach. We assume that the electric field modifies the electrostatic potential linearly along the chain (rod-like chain model ).