First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. I. Generalized Frenkel–Holstein Hamiltonian

A generalized Frenkel-Holstein Hamiltonian is constructed to describe exciton migration in oligo(para-phenylene vinylene) chains, based on excited state electronic structure data for an oligomer comprising 20 monomer units (OPV-20). Time-dependent density functional theory calculations using the ωB97XD hybrid functional are employed in conjunction with a transition density analysis to study the low-lying singlet excitations and demonstrate that these can be characterized to a good approximation as a Frenkel exciton manifold. Based on these findings, we employ the analytic mapping procedure of Binder et al. [J. Chem. Phys. 141, 014101 (2014)] to translate one-dimensional (1D) and two-dimensional (2D) potential energy surface (PES) scans to a fully anharmonic, generalized Frenkel-Holstein (FH) Hamiltonian. A 1D PES scan is carried out for intra-ring quinoid distortion modes, while 2D PES scans are performed for the anharmonically coupled inter-monomer torsional and vinylene bridge bond length alternation modes. The kinetic energy is constructed in curvilinear coordinates by an exact numerical procedure, using the TNUM Fortran code. As a result, a fully molecular-based, generalized FH Hamiltonian is obtained, which is subsequently employed for quantum exciton dynamics simulations, as shown in Paper II [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)].