Exciton Dynamics in Disordered Poly(p-phenylenevinylene). 1. Ultrafast Interconversion and Dynamical Localization

The Role of Double Excitations in Exciton Dynamics of Multiazobenzenes: Trisazobenzenophane as a Test Case

Molecular exciton dynamics underlie energy and charge transfer processes in organic multichromophoric systems. A particularly interesting class of the latter is multiphotochromic systems made of molecules capable of photochemical transformations. Exciton dynamics in assemblies of photoswitches have been recently investigated using either the molecular exciton model or supermolecular configuration interaction (CI) singles, both approaches being based on a semiempirical Hamiltonian and combined with surface hopping molecular dynamics. Here, we study how inclusion of double excitations in nonadiabatic dynamics simulations affects exciton dynamics of multiazobenzenes, using trisazobenzenophane as an example. We find that both CI singles and CI singles and doubles yield virtually the same time scale of dynamical exciton localization, ∼50 fs for the studied multiazobenzene. However, inclusion of double excitations considerably affects the excited state lifetimes and isomerization quantum yields.

Ultrafast Fluorescence Depolarization in Conjugated Polymers

We report on large-scale simulations of intrachain exciton dynamics in poly(para-phenylenevinylene). Our theoretical model describes Frenkel exciton coupling to both fast, quantized C-C bond vibrations and slow, classical torsional modes. We also incorporate system-bath interactions. The dynamics is simulated using the time evolution block decimation method, which avoids the failures of the Ehrenfest approximation to describe decoherence processes and nonadiabatic interstate conversion. System-bath interactions are modeled using quantum trajectories and Lindblad quantum jump operators. We find that following photoexcitation, the quantum mechanical entanglement of the exciton and C-C bond phonons causes exciton-site decoherence. Next, system-bath interactions cause the stochastic collapse of high-energy delocalized excitons into chromophores. Finally, torsional relaxation causes additional exciton-density localization. We relate these dynamical processes to the predicted fluorescence depolarization, extract the time scales corresponding to them, and thus interpret the observed sub-ps fluorescence depolarization.

First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature

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 cm² s^-1 to D ∼ 1 × 10^-3 cm² 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.

Structural and Photophysical Templating of Conjugated Polyelectrolytes with Single-Stranded DNA

A promising approach to influence and control the photophysical properties of conjugated polymers is directing their molecular conformation by templating. We explore here the templating effect of single-stranded DNA oligomers (ssDNAs) on cationic polythiophenes with the goal to uncover the intermolecular interactions that direct the polymer backbone conformation. We have comprehensively characterized the optical behavior and structure of the polythiophenes in conformationally distinct complexes depending on the sequence of nucleic bases and addressed the effect on the ultrafast excited-state relaxation. This, in combination with molecular dynamics simulations, allowed us a detailed atomistic-level understanding of the structure-property correlations. We find that electrostatic and other noncovalent interactions direct the assembly with the polymer, and we identify that optimal templating is achieved with (ideally 10-20) consecutive cytosine bases through numerous π-stacking interactions with the thiophene rings and side groups of the polymer, leading to a rigid assembly with ssDNA, with highly ordered chains and unique optical signatures. Our insights are an important step forward in an effective approach to structural templating and optoelectronic control of conjugated polymers and organic materials in general.

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)].

First-principles quantum simulations of exciton diffusion on a minimal oligothiophene chain at finite temperature

High-dimensional multiconfigurational quantum dynamics simulations are carried out at finite temperature to simulate exciton diffusion on an oligothiophene chain, representative of a segment of the poly(3-hexylthiophene) (P3HT) polymer. The ab initio parametrized site-based Hamiltonian of Binder et al. [Phys. Rev. Lett., 2018, 120, 227401] is employed to model a 20-site system, including intra-ring and inter-ring high-frequency modes as well as torsional modes which undergo thermal fluctuations induced by an explicit harmonic oscillator bath. The system-bath dynamics is treated within the setting of a stochastic mean-field Schrödinger equation. For the 20-site excitonic system, a total of 20 Frenkel states and 248 modes are propagated using the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method. The resulting dynamics can be interpreted in terms of the coherent motion of an exciton-polaron quasi-particle stochastically driven by torsional fluctuations. This dynamics yields a near-linear mean squared displacement (MSD) as a function of time, from which a diffusion coefficient can be deduced which increases with temperature, up to 5.7 × 10-3 cm2 s-1 at T = 300 K.

Conformational Dynamics Guides Coherent Exciton Migration in Conjugated Polymer Materials: First-Principles Quantum Dynamical Study

We report on high-dimensional quantum dynamical simulations of photoinduced exciton migration in a single-chain oligothiophene segment, in view of elucidating the controversial nature of the elementary exciton transport steps in semiconducting polymers. A novel first-principles parametrized Frenkel J aggregate Hamiltonian is employed that goes significantly beyond the standard Frenkel-Holstein Hamiltonian. Departing from a nonequilibrium state created by photoexcitation, these simulations provide evidence of an ultrafast two-timescale process at low temperatures, involving exciton-polaron formation within tens of femtoseconds (fs), followed by torsional relaxation on an ∼400  fs timescale. The second step is the driving force for exciton migration, as initial conjugation breaks are removed by dynamical planarization. The quantum coherent nature of the elementary exciton migration step is consistent with experimental observations highlighting the correlated and vibrationally coherent nature of the dynamics on ultrafast timescales.

Excited-state structural relaxation and exciton delocalization dynamics in linear and cyclic π-conjugated oligothiophenes

π-Conjugated oligothiophene is considered a chain segment of its polymeric counterpart with simper excited-state dynamics and spectral signatures.

Non-adiabatic excited state molecular dynamics of phenylene ethynylene dendrimer using a multiconfigurational Ehrenfest approach

Photoinduced dynamics of electronic and vibrational unidirectional energy transfer between meta-linked building blocks in a phenylene ethynylene dendrimer is simulated using a multiconfigurational Ehrenfest in time-dependent diabatic basis (MCE-TDDB) method.

Simulations of singlet exciton diffusion in organic semiconductors: a review

Recent advances in exciton diffusion simulations in conjugated materials are presented in this review.

Strong Photophysical Similarities between Conjugated Polymers and J-aggregates

The photophysical properties of emissive conjugated polymer (CP) chains are compared to those of linear J-aggregates. The two systems share many properties in common, including a red-shifted absorption spectrum with increasing chain/aggregate length, enhanced radiative decay rates (superradiance) relative to a single monomer/molecule, and several vibronic signatures involving the vinyl-stretching mode common to many conjugated molecules. In particular, the scaling of the 0-0/0-1 photoluminescence ratio and radiative decay rate with the inverse square root of temperature in red-phase polydiacetylene is also characteristic of linear, disorder-free J-aggregates. The strong photophysical resemblance is traced to the excitonic band structure; in one-dimensional direct band gap semiconductors as well as J-aggregates, the exciton band curvature is positive at the gamma point (k = 0).

Coarse-grained simulations of the solution-phase self-assembly of poly(3-hexylthiophene) nanostructures

Under certain conditions the conjugated polymer poly(3-hexylthiophene) (P3HT) self-assembles into high-aspect-ratio nanostructures (known as nanofibres, nanowires, or nanoribbons) when cooled below its solubility limit in a marginal solvent such as anisole. Such nanostructures are potentially beneficial for organic photovoltaic device performance. In this work, Langevin dynamics simulations of a coarse-grained model of P3HT in implicit anisole solvent are used to study the self-assembly of P3HT nanostructures for polymer chain lengths and concentrations used experimentally to prepare P3HT nanofibres. The coarse-grained model is parametrised to match the local structure and dynamics of an atomistic model with explicit solvent. Nanofibres are also prepared experimentally and characterised by atomic force microscopy and UV-vis spectroscopy. The simulations match the experimental phase behaviour of P3HT in anisole, showing aggregation of P3HT at 293 and 308 K but not at 323 or 353 K. Single-chain simulations at 293 K reveal two distinct nano-scale aggregate morphologies: hairpins and helices. Hairpin aggregates, which are the precursors of nanofibres, are slightly favoured energetically at 293 K for nuclei of the critical size of ≈80 monomers for aggregation. Consequently, chains in multi-chain aggregates adopt the hairpin morphology exclusively in simulations at experimental concentrations at 293 K. The simulated aggregate sizes match experimentally measured nanofibre widths. An estimate of the shift in UV-vis absorption of P3HT due to the change in conjugation length with aggregation in the simulations agrees reasonably well with experiment and shows that most of the spectral red shift that occurs with nanofibre formation is due to increased planarisation of the P3HT chains. In addition to providing insight into the mechanisms of nanofibre formation, the simulations resolve details of the molecular-level organisation of chains in P3HT nanofibres hitherto inaccessible by experiment.

Electronically excited states in poly(p-phenylenevinylene): vertical excitations and torsional potentials from high-level ab initio calculations

Ab initio second-order algebraic diagrammatic construction (ADC(2)) calculations using the resolution of the identity (RI) method have been performed on poly-(p-phenylenevinylene) (PPV) oligomers with chain lengths up to eight phenyl rings. Vertical excitation energies for the four lowest π–π* excitations and geometry relaxation effects for the lowest excited state (S₁) are reported. Extrapolation to infinite chain length shows good agreement with analogous data derived from experiment. Analysis of the bond length alternation (BLA) based on the optimized S₁ geometry provides conclusive evidence for the localization of the defect in the center of the oligomer chain. Torsional potentials have been computed for the four excited states investigated and the transition densities divided into fragment contributions have been used to identify excitonic interactions. The present investigation provides benchmark results, which can be used (i) as reference for lower level methods and (ii) give the possibility to parametrize an effective Frenkel exciton Hamiltonian for quantum dynamical simulations of ultrafast exciton transfer dynamics in PPV type systems.