Monte Carlo simulation based on dynamic disorder model in organic semiconductors: From coherent to incoherent transport

Abstract

Quantum transport in disordered systems is studied using a polaron-based master equation. The polaron approach is capable of bridging the results from the coherent band-like transport regime governed by the Redfield equation to incoherent hopping transport in the classical regime. A non-monotonic dependence of the diffusion coefficient is observed both as a function of temperature and system-phonon coupling strength. In the band-like transport regime, the diffusion coefficient is shown to be linearly proportional to the system-phonon coupling strength and vanishes at zero coupling due to Anderson localization. In the opposite classical hopping regime, we correctly recover the dynamics described by the Fermi's Golden Rule and establish that the scaling of the diffusion coefficient depends on the phonon bath relaxation time. In both the hopping and band-like transport regimes, it is demonstrated that at low temperature, the zero-point fluctuations of the bath lead to non-zero transport rates and hence a finite diffusion constant. Application to rubrene and other organic semiconductor materials shows a good agreement with experimental mobility data. C 2015 AIP Publishing LLC. [http://dx.

The Transient Localization Scenario for Charge Transport in Crystalline Organic Materials

“…For the electron-phonon couplings typical of organics crystals, it is rather the spread ∆X 2 of the wavefunction for a single initial set of molecular displacements which grows with time, leading to a fully delocalized wavefunction in the long time limit (see Ref. [103] for an approach which includes decoherence in order to solve this issue). The interpretation of the time snapshots of the electronic density should therefore taken with care.…”

Drift of charge carriers in crystalline organic semiconductors

“…[33,34] It is recently applied to the crystalline organic semiconductors, resulting in a diffusion with localized states and band-like behavior. [20] Besides, the near-equilibrium Boltzmann distribution is maintained for the electronic system. [19] Here we generalize this method to account for the drift motion in finite electric fields, which is implemented in site representation with Miller-Abrahams type energydependent reweighing factor.…”

“…[17,18] Very recently, two of the present authors proposed an instantaneous decoherence correction (IDC) approach with energy-dependent reweighing factors to account for the decoherence and energy relaxation processes. [19,20] It is shown that this method is able to maintain the near-equilibrium distribution of electronic states in the diffusion dynamics, which is a sign of proper treatment of energy relaxation.…”

A variational surface hopping algorithm for the sub-Ohmic spin-boson model

“…It is thus clear that our surface hopping algorithm can simultaneously include negative and positive temperature dependencies. This points to possible important applications of our algorithm to complicated transport problems in molecular systems, where both the coherent and incoherent mechanisms coexist [39].…”

“…(4) The roles of the bosons are twofold: to initiate the thermally assisted hopping process, and to act as scattering centers for the wave packet propagation in a bandlike process [39]. As opposed to usual schemes for including the temperature effects which often incorporate only one of the two roles, in our algorithm we take both into consideration.…”