Concurrent Effects of Delocalization and Internal Conversion Tune Charge Separation at Regioregular Polythiophene-Fullerene Heterojunctions

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

Polariton enhanced free charge carrier generation in donor-acceptor cavity systems by a second-hybridization mechanism

Cavity quantum electrodynamics has been studied as a potential approach to modify free charge carrier generation in donor-acceptor heterojunctions because of the delocalization and controllable energy level properties of hybridized light-matter states known as polaritons. However, in many experimental systems, cavity coupling decreases charge separation. Here, we theoretically study the quantum dynamics of a coherent and dissipative donor-acceptor cavity system, to investigate the dynamical mechanism and further discover the conditions under which polaritons may enhance free charge carrier generation. We use open quantum system methods based on single-pulse pumping to find that polaritons have the potential to connect excitonic states and charge separated states, further enhancing free charge generation on an ultrafast timescale of several hundred femtoseconds. The mechanism involves polaritons with optimal energy levels that allow the exciton to overcome the high Coulomb barrier induced by electron-hole attraction. Moreover, we propose that a second-hybridization between a polariton state and dark states with similar energy enables the formation of the hybrid charge separated states that are optically active. These two mechanisms lead to a maximum of 50% enhancement of free charge carrier generation on a short timescale. However, our simulation reveals that on the longer timescale of picoseconds, internal conversion and cavity loss dominate and suppress free charge carrier generation, reproducing the experimental results. Thus, our work shows that polaritons can affect the charge separation mechanism and promote free charge carrier generation efficiency, but predominantly on a short timescale after photoexcitation.

Even a little delocalization produces large kinetic enhancements of charge-separation efficiency in organic photovoltaics

In organic photovoltaics, charges can separate efficiently even if their Coulomb attraction is an order of magnitude greater than the available thermal energy. Delocalization has been suggested to explain this fact, because it could increase the initial separation of charges in the charge-transfer (CT) state, reducing their attraction. However, understanding the mechanism requires a kinetic model of delocalized charge separation, which has proven difficult because it involves tracking the correlated quantum-mechanical motion of the electron and the hole in large simulation boxes required for disordered materials. Here, we report the first three-dimensional simulations of charge-separation dynamics in the presence of disorder, delocalization, and polaron formation, finding that even slight delocalization, across less than two molecules, can substantially enhance the charge-separation efficiency, even starting with thermalized CT states. Delocalization does not enhance efficiency by reducing the Coulomb attraction; instead, the enhancement is a kinetic effect produced by the increased overlap of electronic states.

Exciton Dissociation in a Model Organic Interface: Excitonic State-Based Surface Hopping versus Multiconfigurational Time-Dependent Hartree

Quantum dynamical simulations are essential for a molecular-level understanding of light-induced processes in optoelectronic materials, but they tend to be computationally demanding. We introduce an efficient mixed quantum-classical nonadiabatic molecular dynamics method termed eXcitonic state-based Surface Hopping (X-SH), which propagates the electronic Schrödinger equation in the space of local excitonic and charge-transfer electronic states, coupled to the thermal motion of the nuclear degrees of freedom. The method is applied to exciton decay in a 1D model of a fullerene-oligothiophene junction, and the results are compared to the ones from a fully quantum dynamical treatment at the level of the Multilayer Multiconfigurational Time-Dependent Hartree (ML-MCTDH) approach. Both methods predict that charge-separated states are formed on the 10-100 fs time scale via multiple "hot-exciton dissociation" pathways. The results demonstrate that X-SH is a promising tool advancing the simulation of photoexcited processes from the molecular to the true nanomaterials scale.

A New Frontier in Exciton Transport: Transient Delocalization

Efficient exciton transport is crucial to the application of organic semiconductors (OSCs) in light-harvesting devices. While the physics of exciton transport in highly disordered media is well-explored, the description of transport in structurally and energetically ordered OSCs is less established, despite such materials being favorable for devices. In this Perspective we describe and highlight recent research pointing toward a highly efficient exciton transport mechanism which occurs in ordered OSCs, transient delocalization. Here, exciton-phonon couplings play a critical role in allowing localized exciton states to temporarily access higher-energy delocalized states whereupon they move large distances. The mechanism shows great promise for facilitating long-range exciton transport and may allow for improved device efficiencies and new device architectures. However, many fundamental questions on transient delocalization remain to be answered. These questions and suggested next steps are summarized.

The Ultrafast Quantum Dynamics of Photoexcited Adenine-Thymine Basepair Investigated with a Fragment-based Diabatization and a Linear Vibronic Coupling Model

In this contribution we present a quantum dynamical study of the photoexcited hydrogen bonded base pair adenine–thymine (AT) in a Watson–Crick arrangement. To that end, we parametrize Linear Vibronic Coupling (LVC) models with Time-Dependent Density Functional Theory (TD-DFT) calculations, exploiting a fragment diabatization scheme (FrD) we have developed to define diabatic states on the basis of individual chromophores in a multichromophoric system. Wavepacket propagations were run with the multilayer extension of the Multiconfiguration Time-Dependent Hartree method. We considered excitations to the three lowest bright states, a ππ* state of thymine and two ππ* states (L_a and L_b) of adenine, and we found that on the 100 fs time scale the main decay pathways involve intramonomer population transfers toward nπ* states of the same nucleobase. In AT this transfer is less effective than in the isolated nucleobases, because hydrogen bonding destabilizes the nπ* states. The population transfer to the A → T charge transfer state is negligible, making the ultrafast (femtosecond) decay through the proton coupled electron transfer mechanism unlikely, in line with experimental results in apolar solvents. The excitation energy transfer is also very small. We carefully compare the predictions of LVC Hamiltonians obtained with different sets of diabatic states, defined so to match either local states of the two separated monomers or the base pair adiabatic states in the Franck–Condon region. To that end we also extend the flexibility of the FrD-LVC approach, introducing a new strategy to define fragments diabatic states that account for the effect of the rest of the multichromohoric system through a Molecular Mechanics potential.

Long-Range Interactions Boost Singlet Exciton Diffusion in Nanofibers of π-Extended Polymer Chains

Raising the distance covered by singlet excitons during their lifetimes to values maximizing light absorption (a few hundred nm) would solve the exciton diffusion bottleneck issue and lift the constraint for fine (∼10 nm) phase segregation in bulk heterojunction organic solar cells. In that context, the recent report of highly ordered conjugated polymer nanofibers featuring singlet exciton diffusion length, L_D, in excess of 300 nm is both appealing and intriguing [Jin, X.; et al. Science 2018, 360 (6391), 897-900]. Here, on the basis of nonadiabatic molecular dynamics simulations, we demonstrate that singlet exciton diffusion in poly(3-hexylthiophene) (P3HT) fibers is highly sensitive to the interplay between delocalization along the polymer chains and long-range interactions along the stacks. Remarkably, the diffusion coefficient is predicted to rocket by 3 orders of magnitude when going beyond nearest-neighbor intermolecular interactions in fibers of extended (30-mer) polymer chains and to be resilient to interchain energetic and positional disorders.

Ultrafast Charge Transfer and Relaxation at a Donor-Acceptor Interface

The efficiency of charge separation in organic photovoltaic materials is crucially determined by the underlying dynamics of the charge transfer (CT) excitons and their dissociation into free electrons and holes. To unravel the main principles of the underlying mechanism on a molecular level, we construct a toy model of electronically coupled donors interacting with a manifold of CT exciton states. In particular, we set up a ladder of CT site energies to model the exciton dissociation. To mimic the complexity of the exciton dynamics at the donor-acceptor interface, the electronic CT manifold is designed to include two vibrational modes that are vibronically coupled to the excitons. We examine the impact of the electronic and vibrational coherences and the structure of the vibronic manifold on the transfer efficiency and charge recombination. Optimal configurations of the vibronic CT manifold are revealed. In particular, the rate of charge recombination can be minimized when the transient dynamics are carefully explored. Such a toy model can be used as a guide for the design of organic materials for efficient photovoltaic devices.

Fragment Diabatization Linear Vibronic Coupling Model for Quantum Dynamics of Multichromophoric Systems: Population of the Charge-Transfer State in the Photoexcited Guanine-Cytosine Pair

We introduce a method (FrD-LVC) based on a fragment diabatization (FrD) for the parametrization of a linear vibronic coupling (LVC) model suitable for studying the photophysics of multichromophore systems. In combination with effective quantum dynamics (QD) propagations with multilayer multiconfigurational time-dependent Hartree (ML-MCTDH), the FrD-LVC approach gives access to the study of the competition between intrachromophore decays, like those at conical intersections, and interchromophore processes, like exciton localization/delocalization and the involvement of charge-transfer (CT) states. We used FrD-LVC parametrized with time-dependent density functional theory (TD-DFT) calculations, adopting either CAM-B3LYP or ωB97X-D functionals, to study the ultrafast photoexcited QD of a guanine-cytosine (GC) hydrogen-bonded pair, within a Watson-Crick arrangement, considering up to 12 coupled diabatic electronic states and the effect of all of the 99 vibrational coordinates. The bright excited states localized on C and, especially, on G are predicted to be strongly coupled to the G → C CT state, which is efficiently and quickly populated after an excitation to any of the four lowest energy bright local excited states. Our QD simulations show that more than 80% of the excited population on G and ∼50% of that on C decay to this CT state in less than 50 fs. We investigate the role of vibronic effects in the population of the CT state and show that it depends mainly on its large reorganization energy so that it can occur even when it is significantly less stable than the bright states in the Franck-Condon region. At the same time, we document that the formation of the GC pair almost suppresses the involvement of dark nπ* excited states in the photoactivated dynamics.

Modeling and Simulating the Excited-State Dynamics of a System with Condensed Phases: A Machine Learning Approach

Simulating the irreversible quantum dynamics of exciton- and electron-transfer problems poses a nontrivial challenge. Because the irreversibility of the system dynamics is a result of quantum thermal activation and dissipation caused by the surrounding environment, it is necessary to include infinite environmental degrees of freedom in the simulation. Because the capabilities of full quantum dynamics simulations that include the surrounding molecular degrees of freedom are limited, employing a system-bath model is a practical approach. In such a model, the dynamics of excitons or electrons are described by a system Hamiltonian, while the other degrees of freedom that arise from the environmental molecules are described by a harmonic oscillator bath (HOB) and system-bath interaction parameters. By extending on a previous study of molecular liquids [ J. Chem. Theory Comput. 2020, 16, 2099], here, we construct a system-bath model for exciton- and electron-transfer problems by means of a machine learning approach. We determine both the system and system-bath interaction parameters, including the spectral distribution of the bath, using the electronic excitation energies obtained from a quantum mechanics/molecular mechanics (QM/MM) simulation that is conducted as a function of time. Using the analytical expressions of optical response functions, we calculate linear and two-dimensional electronic spectra (2DES) for indocarbocyanine dimers in methanol. From these results, we demonstrate the capability of our approach to elucidate the nonequilibrium exciton dynamics of a quantum system in a nonintuitive manner.

Vibrationally Resolved Absorption Spectra and Exciton Dynamics in Zinc Phthalocyanine Aggregates: Effects of Aggregation Lengths and Remote Exciton Transfer

The vibrationally resolved absorption spectra and exciton dynamics in the α-zinc phthalocyanine aggregates are theoretically investigated by using a non-Markovian stochastic Schrödinger equation. The model Hamiltonian adopted for spectral and dynamic simulations explicitly includes the couplings for both nearest-neighbor and remote exciton transfer, and it is parametrized from first-principles calculations. The results indicate that aggregation lengths and remote exciton transfer significantly influence the relative energy alignment between delocalized Frenkel exciton (FE) and charge transfer (CT) states, which in turn strongly affects the relative intensities of the two absorption peaks in the Q-band region. Analytical formulas are derived to establish quantitative structure-spectra relationships in aggregates, and they offer simple patterns to extract electronic-state properties directly from absorption spectra. The dynamics simulations reveal that the light absorption can directly generate mixed states with both FE and CT features, but it is hard for the photoexcitation from the Q-band region to generate free carriers due to the high energies of charge-separated states.

Quantum Dynamics of Exciton Transport and Dissociation in Multichromophoric Systems

Due to the subtle interplay of site-to-site electronic couplings, exciton delocalization, nonadiabatic effects, and vibronic couplings, quantum dynamical studies are needed to elucidate the details of ultrafast photoinduced energy and charge transfer events in organic multichromophoric systems. In this vein, we review an approach that combines first-principles parameterized lattice Hamiltonians with accurate quantum dynamical simulations using advanced multiconfigurational methods. Focusing on the elementary transfer steps in organic functional materials, we address coherent exciton migration and creation of charge transfer excitons in homopolymers, notably representative of the poly(3-hexylthiophene) material, as well as exciton dissociation at polymer:fullerene heterojunctions. We emphasize the role of coherent transfer, trapping effects due to high-frequency phonon modes, and thermal activation due to low-frequency soft modes that drive a diffusive dynamics.

Exciton Spatial Dynamics and Self-Trapping in Carbon Nanocages

Three-dimensional cage-shaped molecules formed from chainlike structures hold potential as unique optoelectronic materials and host compounds. Their optical, structural, and dynamical features are tunable by changes in shape and size. We perform a comparison of these properties for three sizes of strained conjugated [n.n.n]carbon nanocages composed of three paraphenylene chains (bridges) of length n = 4, 5, or 6. The exciton intramolecular redistribution occurring during nonradiative relaxation has been explored using nonadiabatic excited-state molecular dynamics. Our results provide atomistic insight into the conformational features associated with the observed red- and blue-shift trends in the absorption and fluorescence spectra, respectively, with increasing nanocage size. Their internal conversion processes involve intramolecular energy transfer that leads to exciton self-trapping on a few phenylene units at the center of a single bridge. The dependence of these dynamical features on the size of the nanocage can be used to tune their host-guest chemical properties and their use for organic electronics and catenane-like applications.

Delocalised kinetic Monte Carlo for simulating delocalisation-enhanced charge and exciton transport in disordered materials

Charge transport is well understood in both highly ordered materials (band conduction) or highly disordered ones (hopping conduction). In moderately disordered materials—including many organic semiconductors—the approximations valid in either extreme break down, making it difficult to accurately model the conduction. In particular, describing wavefunction delocalisation requires a quantum treatment, which is difficult in disordered materials that lack periodicity. Here, we present the first three-dimensional model of partially delocalised charge and exciton transport in materials in the intermediate disorder regime. Our approach is based on polaron-transformed Redfield theory, but overcomes several computational roadblocks by mapping the quantum-mechanical techniques onto kinetic Monte Carlo. Our theory, delocalised kinetic Monte Carlo (dKMC), shows that the fundamental physics of transport in moderately disordered materials is that of charges hopping between partially delocalised electronic states. Our results reveal why standard kinetic Monte Carlo can dramatically underestimate mobilities even in disordered organic semiconductors, where even a little delocalisation can substantially enhance mobilities, as well as showing that three-dimensional calculations capture important delocalisation effects neglected in lower-dimensional approximations.

The first three-dimensional model of transport in moderately disordered materials shows that a little delocalisation can dramatically enhance mobilities.

Applying Marcus theory to describe the carrier transports in organic semiconductors: Limitations and beyond

Marcus theory has been successfully applied to molecular design for organic semiconductors with the aid of quantum chemistry calculations for the molecular parameters: the intermolecular electronic coupling V and the intramolecular charge reorganization energy λ. The assumption behind this is the localized nature of the electronic state for representing the charge carriers, being holes or electrons. As far as the quantitative description of carrier mobility is concerned, the direct application of Marcus semiclassical theory usually led to underestimation of the experimental data. A number of effects going beyond such a semiclassical description will be introduced here, including the quantum nuclear effect, dynamic disorder, and delocalization effects. The recently developed quantum dynamics simulation at the time-dependent density matrix renormalization group theory is briefly discussed. The latter was shown to be a quickly emerging efficient quantum dynamics method for the complex system.

Exciton dissociation and charge separation at donor–acceptor interfaces from quantum-classical dynamics simulations

Nonadiabatic dynamics simulations based on the quantum-classical Liouville equation are employed to study the real-time dynamics of exciton dissociation and charge separation at a model donor–acceptor interface.

Coherent Real-Space Charge Transport Across a Donor-Acceptor Interface Mediated by Vibronic Couplings

There is growing experimental and theoretical evidence that vibronic couplings, couplings between electronic and nuclear degrees of freedom, play a fundamental role in ultrafast excited-state dynamics in organic donor-acceptor hybrids. Whereas vibronic coupling has been shown to support charge separation at donor-acceptor interfaces, so far, little is known about its role in the real-space transport of charges in such systems. Here we theoretically study charge transport in thiophene:fullerene stacks using time-dependent density functional tight-binding theory combined with Ehrenfest molecular dynamics for open systems. Our results reveal coherent oscillations of the charge density between neighboring donor sites, persisting for ∼200 fs and promoting charge transport within the polymer stacks. At the donor-acceptor interface, vibronic wave packets are launched, propagating coherently over distances of more than 3 nm into the acceptor region. This supports previous experimental observations of long-range ballistic charge-carrier motion in organic photovoltaic systems and highlights the importance of vibronic coupling engineering as a concept for tailoring the functionality of hybrid organic devices.

Bridge-Mediated Charge Separation in Isomeric N-Annulated Perylene Diimide Dimers

The possibility and rate of charge separation (CS) in donor-bridge-acceptor molecules mainly depend on two factors: electronic coupling and solvent effects. The question of how CS occurred in two identical chromophores is fundamental, as it is particularly interesting for potential molecular electronics applications and the photosynthetic reaction centers (RCs). Conjugated bridge definitely plays a crucial role in electronic coupling. To determine the bridge-mediated charge separation dynamics between the two identical chromophores, the isomeric N-annulated perylene diimide dimers (para-BDNP and meta-BDNP) with different conjugated bridge structures have been comparatively investigated in different solvents using femtosecond transient absorption spectra (fs-TA). It is found that the charge separation is disfavored in weak polar solvent, whereas direct spectroscopic signatures of radicals are observed in polar solvents, and the rate of charge separation increases as the solvent polarity increasing. To our surprise, the rate of charge separation in m-BDNP is more than an order of magnitude slower than that in p-BDNP, although there is a larger negative ΔG_CS in m-BDNP. The slow CS rate that occurred in m-BDNP mainly results from the intrinsic destructive interference of the wave function through the meta-substituted bridge. The roles of solvent effects in free energy and electronic coupling for charge separation are further identified with quantum calculations.

Coherent Charge Transfer Exciton Formation in Regioregular P3HT: A Quantum Dynamical Study

The ultrafast formation of charge transfer excitons (CTXs) in regioregular poly(3-hexyl thiophene) (rrP3HT) domains is elucidated by electronic structure and quantum dynamical studies of an aggregate model system comprising five stacked quaterthiophene units. Using a multistate vibronic coupling Hamiltonian parametrized by TDDFT calculations for 13 electronic states of Frenkel and CTX type, along with 78 vibrational modes, quantum dynamical simulations are carried out using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method. In line with time-resolved spectroscopic results [ De Sio , A. ; et al. Nat. Commun. 2016 , 7 , 13742 ], it is found that CTX formation occurs immediately upon photoexcitation, accompanied by sustained regular oscillations with a ∼22 fs periodicity. These coherent features, whose presence may seem surprising in a high-dimensional aggregate or thin film material, can be traced back to a dominant vibronic signature of CC stretch-type high-frequency modes. These vibrational signatures are found to be enhanced due to a collective vibronic response that is prompted by the initial generation of a delocalized bright exciton and its subsequent relaxation, by internal conversion, to a polaronic local exciton ground state.

The relationship between the coherent size, binding energy and dissociation dynamics of charge transfer excitons at organic interfaces

At organic semiconductor interfaces, an electron and a hole can be bound Coulombically to form an interfacial charge transfer (CT) exciton. It is still under debate how a CT exciton can overcome its strong binding and dissociate into free carriers. Experimentally, capturing the evolution of the CT exciton on time (fs-ps) and length scales (nm) in which the dissociation process occurs is challenging. To overcome this challenge, time-resolved two photon photoemission spectroscopy is used to measure the binding energies and electronic coherent sizes of a series of CT states at organic interfaces, and capture the temporal dynamics of these CT excitons after their excitation. Using zinc phthalocyanine (ZnPc)/fullerene (C₆₀) interface as a model system, it is shown that the interfacial CT process first populates a hot CT state with a coherent size of ~4 nm. Hot and delocalized CT excitons subsequently relax into CT excitons with lower energies and smaller coherent sizes. To correlate the CT exciton properties with the dissociation efficiency, we develop a method that exploits graphene field effect transistors to probe the rate and yield of free carrier generation at the interface. Our results show that exciton dissociation can be more efficient if one can extract electrons from the hot and delocalized CT state. We propose a cascade structure that would serve this purpose.

Charge separation at an organic/inorganic nano-hybrid interface: atomistic simulations of a para-sexiphenyl ZnO system

A prototypical organic/inorganic interface is considered which is formed by vertical stacking of 20 para-sexiphenyl molecules physisorbed on a ZnO nano-cluster of 3903 atoms. Charge separation kinetics at the interface are investigated for their dependence on ultrafast optical excitation. In order to analyze the spatio-temporal evolution of the Frenkel exciton in the organic part and the formation of charge separated states a first principles parameterized Hamiltonian is introduced and the related time-dependent Schroedinger equation is solved. By determining the interface absorption spectrum the optically addressable states can be uncovered. The work continues our previous studies of J. Phys. Chem. Lett., 2018, 9, 209, but with a changed type of surface passivation. This prevents trapping of electrons close to the surface. Charge separated states are formed by direct optical excitation and also by exciton decay at the interface. Electron migration away from the interface into bulk regions becomes possible. The hole stays close to the interface for all excitation scenarios. Finally, it is demonstrated that energetic disorder is of minor influence.

Charge Transfer Driven Structural Relaxation in a Push-Pull Azobenzene Dye-Semiconductor Complex

Photoexcited structural dynamics in azo-compounds may differ fundamentally whether the push-pull photochromic azo-compound is isolated or forms a heterogeneous charge transfer complex, due to a sudden oxidation of the chromophore. Herein, we use a quantum-classical self-consistent approach that incorporates nonadiabatic excited-state electronic quantum dynamics into molecular mechanics to study the photoexcited dynamics of the push-pull azo-compound para-Methyl Red in the gas phase and sensitizing the (101) anatase surface of TiO₂. We find that the photoinduced S₂/S₀ trans-to- cis isomerization of para-Methyl Red in the gas phase occurs through a pedal-like torsion around the ϕ_CNNC dihedral angle, without evidence to support the inversion mechanism, likewise in the parent azobenzene molecule. However, the photoexcited structural relaxation of the charge transfer complex para-Methyl Red/TiO₂ contrasts essentially with the isolated azo-compounds. Immediately after photoexcitation, the excited electron flows into the TiO₂ conduction band, with an injection time constant of ≃5 fs, and no indication of isomerization is observed during the 1.5 ps simulations. Instead, a strong vibronic relaxation occurs that excites the NN stretching mode of the azo-group, which is ultimately ascribed to the NA relaxation, and delocalization, of the hole wavepacket.

Charge Separation in Donor-C60 Complexes with Real-Time Green Functions: The Importance of Nonlocal Correlations

We use the nonequilibrium Green function (NEGF) method to perform real-time simulations of the ultrafast electron dynamics of photoexcited donor-C₆₀ complexes modeled by a Pariser-Parr-Pople Hamiltonian. The NEGF results are compared to mean-field Hartree-Fock (HF) calculations to disentangle the role of correlations. Initial benchmarking against numerically highly accurate time-dependent density matrix renormalization group calculations verifies the accuracy of NEGF. We then find that charge-transfer (CT) excitons partially decay into charge separated (CS) states if dynamical nonlocal correlation corrections are included. This CS process occurs in ∼10 fs after photoexcitation. In contrast, the probability of exciton recombination is almost 100% in HF simulations. These results are largely unaffected by nuclear vibrations; the latter become however essential whenever level misalignment hinders the CT process. The robust nature of our findings indicates that ultrafast CS driven by correlation-induced decoherence may occur in many organic nanoscale systems, but it will only be correctly predicted by theoretical treatments that include time-nonlocal correlations.

Atomistic Simulations of Charge Separation at a Nanohybrid Interface: Relevance of Photoinduced Initial State Preparation

Charge separation kinetics at a nanohybrid interface are investigated in their dependence on ultrafast optical excitation. A prototypical organic/inorganic interface is considered. It is formed by a vertical stacking of 20 para-sexiphenyl molecules physisorbed on a ZnO nanocluster of 3783 atoms. A first principle parametrized Hamiltonian is employed, and the photoinduced subpicosecond evolution of Frenkel-excitons in the organic part is analyzed besides the formation of charge separated states across the interface. The interface absorption spectrum is calculated. Together, the data indicate that the charge separation is based on the direct excitation of the charge separated states but also on the migration of created Frenkel excitons to the interface with subsequent decay. Further, the photoinduced interface dynamics are compared with data resulting from direct set-ups of an initially excited state. Mostly such set-ups lead to substantially different charge separation processes.

Exciton-Coupled Electron Transfer Process Controlled by Non-Markovian Environments

We theoretically investigate an exciton-coupled electron transfer (XCET) process that is conversion of an exciton into a charge transfer state. This conversion happens in an exciton transfer (XT) process, and the electron moves away in an electron transfer (ET) process in multiple environments (baths). This XCET process plays an essential role in the harvesting of solar energy in biological and photovoltaic materials. We develop a practical theoretical model to study the efficiency of the XCET process that occurs either in consecutive or concerted processes under the influence of non-Markovian baths. The role of quantum coherence in the XT-ET system and the baths is investigated using reduced hierarchal equations of motion (HEOM). This model includes independent baths for each XT and ET state, in addition to a XCET bath for the conversion process. We found that, while quantum system-bath coherence is important in the XT and ET processes, coherence between the XT and ET processes must be suppressed in order to realize that an efficient irreversible XCET process through the weak off-diagonal interaction between the XT and ET bridge sites arises from an XCET bath.

Charge and energy migration in molecular clusters: A stochastic Schrödinger equation approach

The performance of stochastic Schrödinger equations for simulating dynamic phenomena in large scale open quantum systems is studied. Going beyond small system sizes, commonly used master equation approaches become inadequate. In this regime, wave function based methods profit from their inherent scaling benefit and present a promising tool to study, for example, exciton and charge carrier dynamics in huge and complex molecular structures. In the first part of this work, a strict analytic derivation is presented. It starts with the finite temperature reduced density operator expanded in coherent reservoir states and ends up with two linear stochastic Schrödinger equations. Both equations are valid in the weak and intermediate coupling limit and can be properly related to two existing approaches in literature. In the second part, we focus on the numerical solution of these equations. The main issue is the missing norm conservation of the wave function propagation which may lead to numerical discrepancies. To illustrate this, we simulate the exciton dynamics in the Fenna-Matthews-Olson complex in direct comparison with the data from literature. Subsequently a strategy for the proper computational handling of the linear stochastic Schrödinger equation is exposed particularly with regard to large systems. Here, we study charge carrier transfer kinetics in realistic hybrid organic/inorganic para-sexiphenyl/ZnO systems of different extension.

Vibronic coupling in organic semiconductors for photovoltaics

Ultrafast two-dimensional electronic spectroscopy reveals vibronically-assisted coherent charge transport and separation in organic materials and opens up new perspectives for artificial light-to-current conversion.

Ab initio molecular dynamics of thiophene: the interplay of internal conversion and intersystem crossing

Ab initio on-the-fly molecular dynamics reveals that excited thiophene decays via low lying conical intersections and via intersystem crossing. Open-ring structures are responsible for the observed long life times.

Ultrafast decoherence dynamics govern photocarrier generation efficiencies in polymer solar cells

All-organic-based photovoltaic solar cells have attracted considerable attention because of their low-cost processing and short energy payback time. In such systems the primary dissociation of an optical excitation into a pair of photocarriers has been recently shown to be extremely rapid and efficient, but the physical reason for this remains unclear. Here, two-dimensional photocurrent excitation spectroscopy, a novel non-linear optical spectroscopy, is used to probe the ultrafast coherent decay of photoexcitations into charge-producing states in a polymer:fullerene based solar cell. The two-dimensional photocurrent spectra are interpreted by introducing a theoretical model for the description of the coupling of the electronic states of the system to an external environment and to the applied laser fields. The experimental data show no cross-peaks in the twodimensional photocurrent spectra, as predicted by the model for coherence times between the exciton and the photocurrent producing states of 20 fs or less.

Molecular Packing Determines Charge Separation in a Liquid Crystalline Bisthiophene-Perylene Diimide Donor-Acceptor Material

Combined electronic structure and quantum dynamical calculations are employed to investigate charge separation in a novel class of covalently bound bisthiophene-perylene diimide type donor-acceptor (DA) co-oligomer aggregates. In an earlier spectroscopic study of this DA system in a smectic liquid crystalline (LC) film, efficient and ultrafast (subpicosecond) initial charge separation was found to be followed by rapid recombination. By comparison, the same DA system in solution exhibits ultrafast resonant energy transfer followed by slower (picosecond scale) charge separation. The present first-principles study explains these contrasting observations, highlighting the role of an efficient intermolecular charge-transfer pathway that results from the molecular packing in the LC phase. Despite the efficiency of this primary charge-transfer step, long-range charge separation is impeded by a comparatively high Coulomb barrier in conjunction with small electron- and hole-transfer integrals. Quantum dynamical calculations are carried out for a fragment-based model Hamiltonian, parametrized by ab initio second-order Algebraic Diagrammatic Construction (ADC(2)) and Time-Dependent Density Functional Theory (TDDFT) electronic structure calculations. Simulations of coherent vibronic quantum dynamics for up to 156 electronic states and 48 modes are performed using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method. Excellent agreement with experimentally determined charge separation time scales is obtained, and the spatially coherent nature of the dynamics is analyzed.

Charge Separation and Recombination at Polymer-Fullerene Heterojunctions: Delocalization and Hybridization Effects

We address charge separation and recombination in polymer/fullerene solar cells with a multiscale modeling built from accurate atomistic inputs and accounting for disorder, interface electrostatics and genuine quantum effects on equal footings. Our results show that bound localized charge transfer states at the interface coexist with a large majority of thermally accessible delocalized space-separated states that can be also reached by direct photoexcitation, thanks to their strong hybridization with singlet polymer excitons. These findings reconcile the recent experimental reports of ultrafast exciton separation ("hot" process) with the evidence that high quantum yields do not require excess electronic or vibrational energy ("cold" process), and show that delocalization, by shifting the density of charge transfer states toward larger effective electron-hole radii, may reduce energy losses through charge recombination.