First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature.
Rainer HeggerRobert BinderIrene BurghardtPublished in: Journal of chemical theory and computation (2020)
We report on first-principles quantum-dynamical and quantum-classical simulations of photoinduced exciton dynamics in oligothiophene chain segments, representative of intrachain exciton migration in the poly(3-hexylthiophene) (P3HT) polymer. Following up on our recent study (Binder R.; Burghardt, I. Faraday Discuss. 2020, 221, 406), multilayer multiconfiguration time-dependent Hartree calculations for a short oligothiophene segment comprising 20 monomer units (OT-20) are carried out to obtain full quantum-dynamical simulations at finite temperature. These are employed to benchmark mean-field Ehrenfest calculations, which are shown to give qualitatively correct results for the present system. Periodic boundary conditions turn out to significantly improve earlier estimates of diffusion coefficients. Using the Ehrenfest approach, a series of calculations are subsequently carried out for larger lattices (OT-40 to OT-80), leading to estimates for temperature-dependent mean-squared displacements, which are found to exhibit a near-linear dependence as a function of time. The resulting diffusion coefficient estimates are an increasing function of temperature, whose detailed functional form depends on the degree of static disorder. With a realistic static disorder parameter (σs ≃ 0.06 eV), the diffusion coefficients decrease from D ∼ 1 × 10-2 cm2 s-1 to D ∼ 1 × 10-3 cm2 s-1, in qualitative agreement with experimental data for P3HT. The dynamical scenario obtained from our simulations shows that exciton migration in P3HT-type chains is a largely adiabatic process throughout the temperature regime we investigated (i.e., T = 50-300 K). The resulting picture of exciton migration is a coherent, but not bandlike, motion of an exciton-polaron driven by fluctuations induced by low-frequency modes. This process acquires partial hopping character if static disorder becomes prominent and Anderson localization sets in.