Molecular Modeling and Simulation of Conjugated Polymer Oligomers: Ground and Excited State Chain Dynamics of PPV in the Gas Phase

Ultrafast Hot Exciton Nonadiabatic Excited-State Dynamics in Green Fluorescent Protein Chromophore Analogue

The ultrafast high-energy nonadiabatic excited-state dynamics of the benzylidenedimethylimidazolinone chromophore dimer has been investigated using an electronic structure method coupled with on-the-fly quantitative wave function analysis to gain insight into the photophysics of hot excitons in biological systems. The dynamical simulation provides a rationalization of the behavior of the exciton in a dimer after the photoabsorption of light to higher-energy states. The results suggest that hot exciton localization within the manifold of excited states is caused by the hindrance of torsional rotation due to imidazolinone (I) or phenolate (P) bonds i.e., ΦI- or ΦP-dihedral rotation, in the monomeric units of a dimer. This hindrance arises due to weak π-π stacking interaction in the dimer, resulting in an energetically uphill excited-state barrier for ΦI- and ΦP-twisted rotation, impeding the isomerization process in the chromophore. Thus, this study highlights the potential impact of the weak π-π interaction in regulating the photodynamics of the green fluorescent protein chromophore derivatives.

Optoelectronic Current through Unbiased Monolayer Amorphous Carbon Nanojunctions

Monolayer amorphous carbon (MAC) is a recently synthesized disordered 2D carbon material. An ensemble of MAC nanofragments contains diverse manifestations of lattice disorder, and because of disorder the key unifying characteristic of this ensemble is poor electronic conductance. Curiously, our computational analysis of the electronic properties of MAC nanofragments revealed an additional commonality: a robust presence of charge-transfer character for electronic transitions from the occupied to virtual orbitals. This charge-transfer property suggests possible applications in optoelectronics. In this Letter, we demonstrate computationally that a laser pulse induces directional electronic currents in unbiased MAC nanojunctions and discuss the effects of pulse intensity on the magnitude of electron transfer.

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

How to calculate charge mobility in molecular materials from surface hopping non-adiabatic molecular dynamics - beyond the hopping/band paradigm

Charge transport in high mobility organic semiconductors is in an intermediate regime between small polaron hopping and band transport limits. We have recently shown that surface hopping non-adiabatic molecular dynamics is a powerful method for prediction of charge transport mechanisms in organic materials and for near-quantitative prediction of charge mobilities at room temperature where the effects of nuclear zero-point motion and tunneling are still relatively small [S. Giannini et al., Nat. Commun., 2019, 10, 3843]. Here we assess and critically discuss the extensions to Tully's original method that have led to this success: (i) correction for missing electronic decoherence, (ii) detection of trivial crossings and (iii) removal of decoherence correction-induced spurious charge transfer. If any one of these corrections is not included, the charge mobility diverges with system size, each for different physical reasons. Yet if they are included, convergence with system size, detailed balance and good internal consistency are achieved.

Detailed balance, internal consistency, and energy conservation in fragment orbital-based surface hopping

We have recently introduced an efficient semi-empirical non-adiabatic molecular dynamics method for the simulation of charge transfer/transport in molecules and molecular materials, denoted fragment orbital-based surface hopping (FOB-SH) [J. Spencer et al., J. Chem. Phys. 145, 064102 (2016)]. In this method, the charge carrier wavefunction is expanded in a set of charge localized, diabatic electronic states and propagated in the time-dependent potential due to classical nuclear motion. Here we derive and implement an exact expression for the non-adiabatic coupling vectors between the adiabatic electronic states in terms of nuclear gradients of the diabatic electronic states. With the non-adiabatic coupling vectors (NACVs) available, we investigate how different flavours of fewest switches surface hopping affect detailed balance, internal consistency, and total energy conservation for electron hole transfer in a molecular dimer with two electronic states. We find that FOB-SH satisfies detailed balance across a wide range of diabatic electronic coupling strengths provided that the velocities are adjusted along the direction of the NACV to satisfy total energy conservation upon a surface hop. This criterion produces the right fraction of energy-forbidden (frustrated) hops, which is essential for correct population of excited states, especially when diabatic couplings are on the order of the thermal energy or larger, as in organic semiconductors and DNA. Furthermore, we find that FOB-SH is internally consistent, that is, the electronic surface population matches the average quantum amplitudes, but only in the limit of small diabatic couplings. For large diabatic couplings, inconsistencies are observed as the decrease in excited state population due to frustrated hops is not matched by a corresponding decrease in quantum amplitudes. The derivation provided here for the NACV should be generally applicable to any electronic structure approach where the electronic Hamiltonian is constructed in a diabatic electronic state basis.

Charge Transport in Molecular Materials: An Assessment of Computational Methods

The booming field of molecular electronics has fostered a surge of computational research on electronic properties of organic molecular solids. In particular, with respect to a microscopic understanding of transport and loss mechanisms, theoretical studies assume an ever-increasing role. Owing to the tremendous diversity of organic molecular materials, a great number of computational methods have been put forward to suit every possible charge transport regime, material, and need for accuracy. With this review article we aim at providing a compendium of the available methods, their theoretical foundations, and their ranges of validity. We illustrate these through applications found in the literature. The focus is on methods available for organic molecular crystals, but mention is made wherever techniques are suitable for use in other related materials such as disordered or polymeric systems.

An insight into non-emissive excited states in conjugated polymers

Conjugated polymers in the solid state usually exhibit low fluorescence quantum yields, which limit their applications in many areas such as light-emitting diodes. Despite considerable research efforts, the underlying mechanism still remains controversial and elusive. Here, the nature and properties of excited states in the archetypal polythiophene are investigated via aggregates suspended in solvents with different dielectric constants (ɛ). In relatively polar solvents (ɛ>∼ 3), the aggregates exhibit a low fluorescence quantum yield (QY) of 2–5%, similar to bulk films, however, in relatively nonpolar solvents (ɛ<∼ 3) they demonstrate much higher fluorescence QY up to 20–30%. A series of mixed quantum-classical atomistic simulations illustrate that dielectric induced stabilization of nonradiative charge-transfer (CT) type states can lead to similar drastic reduction in fluorescence QY as seen experimentally. Fluorescence lifetime measurement reveals that the CT-type states exist as a competitive channel of the formation of emissive exciton-type states.

Conjugated polymers in thin films exhibit low fluorescence quantum yields, but the mechanism is still unclear. Here, Hu et al. show the trade-off between charge transfer and emissive exciton states, whilst the former can be suppressed via dielectric-induced stabilization for large fluorescence quantum yields.

Excited state dynamics of organic semi-conducting materials

Time-resolved absorption and emission spectroscopy has been applied to investigate the dynamics of excited state processes in oligomer models for semi-conducting organic materials. Following the photo-excitation of a pentamer oligomer that is a model for the conjugated polymer MEH-PPV, an ultrafast component of a few picoseconds is observed for the decay of the initially formed transient species. Variable temperature absorption and emission spectra combined with X-ray crystallography and calculations support the assignment of this rapid relaxation process to an excited state conformational rearrangement from non-planar to more planar molecular configurations. The implications of the results for the overall photophysics of conjugated polymers are considered.

Absorption and fluorescence spectra of poly(p-phenylenevinylene) (PPV) oligomers: an ab initio simulation

The absorption and fluorescence spectra of poly(p-phenylenevinylene) (PPV) oligomers with up to seven repeat units were theoretically investigated using the algebraic diagrammatic construction method to second order, ADC(2), combined with the resolution-of-the-identity (RI) approach. The ground and first excited state geometries of the oligomers were fully optimized. Vertical excitation energies and oscillator strengths of the first four transitions were computed. The vibrational broadening of the absorption and fluorescence spectra was studied using a semiclassical nuclear ensemble method. After correcting for basis set and solvent effects, we achieved a balanced description of the absorption and fluorescence spectra by means of the ADC(2) approach. This fact is documented by the computed Stokes shift along the PPV series, which is in good agreement with the experimental values. The experimentally observed band width of the UV absorption and fluorescence spectra is well reproduced by the present simulations showing that the nuclear ensemble generated should be well suitable for consecutive surface hopping dynamics simulations.

Intermolecular interactions and charge transfer transitions in aromatic hydrocarbon–tetracyanoethylene complexes

ADC(2) calculations accurately describe charge transfer transitions in complexes of the tetracyanoethylene electron acceptor and three distinct aromatic donors.

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.

Coherent Excitation Transfer Driven by Torsional Dynamics: a Model Hamiltonian for PPV Type Systems

A model Hamiltonian is proposed for systems of poly-phenylene-vinylene (PPV) type, which describes excitation energy transfer (EET) in the presence of torsional and bond-length-alternation coordinates that couple neighboring fragments. This model is employed to describe coherent EET, starting from a partially delocalized exciton state which is initially separated from the rest of the chain by a defect due to a pronounced torsional displacement at one of the inter-fragment junctions. EET is found to be driven by the restoring force of the torsions towards the planar geometry. Quantum dynamical calculations are carried out for a minimal five-site model, using the multiconfiguration time-dependent Hartree (MCTDH) method.

Prediction of Excitation Energies for Conjugated Oligomers and Polymers from Time-Dependent Density Functional Theory

With technological advances, light-emitting conjugated oligomers and polymers have become competitive candidates in the commercial market of light-emitting diodes for display and other technologies, due to the ultralow cost, light weight, and flexibility. Prediction of excitation energies of these systems plays a crucial role in the understanding of their optical properties and device design. In this review article, we discuss the calculation of excitation energies with time-dependent density functional theory, which is one of the most successful methods in the investigation of the dynamical response of molecular systems to external perturbation, owing to its high computational efficiency.