Quantum dynamics of ultrafast charge transfer at an oligothiophene-fullerene heterojunction

Kylin-V: An open-source package calculating the dynamic and spectroscopic properties of large systems

Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples.

Ultrafast Excited-State Energy Transfer in Phenylene Ethynylene Dendrimer: Quantum Dynamics with the Tensor Network Method

Photoinduced excited-state energy transfer (EET) processes play an important role in solar energy conversions. Owing to their excellent photoharvesting and exciton-transport properties, phenylene ethynylene (PE) dendrimers display great potential for improving the efficiency of solar cells. In this work, we investigated the intramolecular EET dynamics in a dendrimer composed of two linear PE units (2-ring and 3-ring) using a fully quantum description based on the tensor network method. We first constructed a diabatic model Hamiltonian based on the electronic structure calculations. Using this diabatic vibronic coupling model, we tried to obtain the main features of the EET dynamics in terms of the several diabatic models with different numbers of vibrational modes (from 4 modes to 129 modes) and to explore the corresponding vibronic coupling interactions. The results show that the EET in this PE dendrimer is ultrafast. Four modes of A' symmetry play dominant roles in the dynamics; the remaining 86 modes of A' symmetry can dampen the electronic coherence; and the modes of A″ symmetry do not exhibit significant influence on the EET process. Overall, the first-order intrastate vibronic coupling terms show the dominant role in the EET dynamics, while the second-order intrastate vibronic coupling terms cause damping of the electronic coherence and slow down the overall EET process. This work provides a microscopic understanding of the EET dynamics in PE dendrimers.

Managing temperature in open quantum systems strongly coupled with structured environments

In non-perturbative non-Markovian open quantum systems, reaching either low temperatures with the hierarchical equations of motion (HEOM) or high temperatures with the Thermalized Time Evolving Density Operator with Orthogonal Polynomials Algorithm (T-TEDOPA) formalism in Hilbert space remains challenging. We compare different ways of modeling the environment. Sampling the Fourier transform of the bath correlation function, also called temperature dependent spectral density, proves to be very effective. T-TEDOPA [Tamascelli et al., Phys. Rev. Lett. 123, 090402 (2019)] uses a linear chain of oscillators with positive and negative frequencies, while HEOM is based on the complex poles of an optimized rational decomposition of the temperature dependent spectral density [Xu et al., Phys. Rev. Lett. 129, 230601 (2022)]. Resorting to the poles of the temperature independent spectral density and of the Bose function separately is an alternative when the problem due to the huge number of Bose poles at low temperatures is circumvented. Two examples illustrate the effectiveness of the HEOM and T-TEDOPA approaches: a benchmark pure dephasing case and a two-bath model simulating the dynamics of excited electronic states coupled through a conical intersection. We show the efficiency of T-TEDOPA to simulate dynamics at a finite temperature by using either continuous spectral densities or only all the intramolecular oscillators of a linear vibronic model calibrated from ab initio data of a phenylene ethynylene dimer.

Coherent Vibrations Promote Charge-Transfer across a Graphene-Based Interface

Discerning the impact of the coherent motion of the nuclei on the timing and efficiency of charge transfer at the donor-acceptor interface is essential for designing performance-enhanced optoelectronic devices. Here, we employ an experimental approach using photocurrent detection in coherent multidimensional spectroscopy to excite a donor aromatic macrocycle and collect the charge transferred to a 2D acceptor layer. For this purpose, we prepared a cobalt phthalocyanine-graphene (CoPc-Gr) interface. Unlike blends, the well-ordered architecture achieved through the physical separation of the two layers allows us to unambiguously collect the electrical signal from graphene alone and associate it with a microscopic understanding of the whole process. The CoPc-Gr interface exhibits an ultrafast electron-transfer signal, stemming from an interlayer mechanism. Remarkably, the signal presents an oscillating time evolution modulated by coherent vibrations originating from the laser-excited CoPc states. By performing Fourier analysis on the beatings and correlating it with the Raman features, along with a comprehensive first-principles characterization of the vibrational coupling in the CoPc excited states, we successfully identify both the orbitals and molecular vibrations that promote the charge transfer at the interface.

Exact Non-Markovian Quantum Dynamics on the NISQ Device Using Kraus Operators

The theory of open quantum systems has many applications ranging from simulating quantum dynamics in condensed phases to better understanding quantum-enabled technologies. At the center of theoretical chemistry are the developments of methodologies and computational tools for simulating charge and excitation energy transfer in solutions, biomolecules, and molecular aggregates. As a variety of these processes display non-Markovian behavior, classical computer simulation can be challenging due to exponential scaling with existing methods. With quantum computers holding the promise of efficient quantum simulations, in this paper, we present a new quantum algorithm based on Kraus operators that capture the exact non-Markovian effect at a finite temperature. The implementation of the Kraus operators on the quantum machine uses a combination of singular value decomposition (SVD) and optimal Walsh operators that result in shallow circuits. We demonstrate the feasibility of the algorithm by simulating the spin-boson dynamics and the exciton transfer in the Fenna-Matthews-Olson (FMO) complex. The NISQ results show very good agreement with the exact ones.

Controlling the nonadiabatic dynamics of the charge-transfer process with chirped pulses: Insights from a double-pump time-resolved fluorescence spectroscopy scheme

The manipulation of the ultrafast quantum dynamics of a molecular system can be achieved through the application of tailored light fields. This has been done in many ways in the past. In our present investigation, we show that it is possible to exert specific control over the nonadiabatic dynamics of a generic model system describing ultrafast charge-transfer within a condensed dissipative environment by using frequency-chirped pulses. By adjusting the external photoexcitation conditions, such as the chirp parameter, we show that the final population of the excitonic and charge-transfer states can be significantly altered, thereby influencing the elementary steps controlling the transfer process. In addition, we introduce an excitation scheme based on double-pump time-resolved fluorescence spectroscopy using chirped-pulse excitations. Here, our findings reveal that chirped excitations enhance the vibrational system dynamics as evidenced by the simulated spectra, where a substantial signal intensity dependence on the chirp is observed. Our simulations show that chirped pulses are a promising tool for steering the dynamics of the charge-transfer process toward a desired target outcome.

Simulations of molecular photodynamics in long timescales

Nonadiabatic dynamics simulations in the long timescale (much longer than 10 ps) are the next challenge in computational photochemistry. This paper delimits the scope of what we expect from methods to run such simulations: they should work in full nuclear dimensionality, be general enough to tackle any type of molecule and not require unrealistic computational resources. We examine the main methodological challenges we should venture to advance the field, including the computational costs of the electronic structure calculations, stability of the integration methods, accuracy of the nonadiabatic dynamics algorithms and software optimization. Based on simulations designed to shed light on each of these issues, we show how machine learning may be a crucial element for long time-scale dynamics, either as a surrogate for electronic structure calculations or aiding the parameterization of model Hamiltonians. We show that conventional methods for integrating classical equations should be adequate to extended simulations up to 1 ns and that surface hopping agrees semiquantitatively with wave packet propagation in the weak-coupling regime. We also describe our optimization of the Newton-X program to reduce computational overheads in data processing and storage. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

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.

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.

Funneling dynamics in a phenylacetylene trimer: Coherent excitation of donor excitonic states and their superposition

Funneling dynamics in conjugated dendrimers has raised great interest in the context of artificial light-harvesting processes. Photoinduced relaxation has been explored by time-resolved spectroscopy and simulations, mainly by semiclassical approaches or referring to open quantum systems methods, within the Redfield approximation. Here, we take the benefit of an ab initio investigation of a phenylacetylene trimer, and in the spirit of a divide-and-conquer approach, we focus on the early dynamics of the hierarchy of interactions. We build a simplified but realistic model by retaining only bright electronic states and selecting the vibrational domain expected to play the dominant role for timescales shorter than 500 fs. We specifically analyze the role of the in-plane high-frequency skeletal vibrational modes involving the triple bonds. Open quantum system non-adiabatic dynamics involving conical intersections is conducted by separating the electronic subsystem from the high-frequency tuning and coupling vibrational baths. This partition is implemented within a robust non-perturbative and non-Markovian method, here the hierarchical equations of motion. We will more precisely analyze the coherent preparation of donor states or of their superposition by short laser pulses with different polarizations. In particular, we extend the π-pulse strategy for the creation of a superposition to a V-type system. We study the relaxation induced by the high-frequency vibrational collective modes and the transitory dissymmetry, which results from the creation of a superposition of electronic donor states.

Vibration-mediated resonant charge separation across the donor-acceptor interface in an organic photovoltaic device

Examination of a recent open-system Ehrenfest dynamics simulation suggests that a vibration-mediate resonance may play a pivotal role in the charge transfer across a donor-acceptor interface in an organic solar cell. Based on this, a concise dissipative two-level electronic system coupled to a molecular vibrational mode is proposed and solved quantum mechanically. It is found that the charge transfer is enhanced substantially when the vibrational energy quanta is equal to the electronic energy loss across the interface. This vibration-mediate resonant charge transfer process is ultrafast, occurring within 100 fs, comparable to experimental findings. The open-system Ehrenfest dynamics simulation of the two-level model is carried out, and similar results are obtained, which confirms further that the earlier open-system Ehrenfest dynamics simulation indeed correctly predicted the occurrence of the resonant charge transfer across the donor-acceptor interface.

Multi-layer Gaussian-based multi-configuration time-dependent Hartree (ML-GMCTDH) simulations of ultrafast charge separation in a donor-acceptor complex

We report on first applications of the Multi-Layer Gaussian-based Multi-Configuration Time-Dependent Hartree (ML-GMCTDH) method [Römer et al., J. Chem. Phys. 138, 064106 (2013)] beyond its basic two-layer variant. The ML-GMCTDH scheme provides an embedding of a variationally evolving Gaussian wavepacket basis into a hierarchical tensor representation of the wavefunction. A first-principles parameterized model Hamiltonian for ultrafast non-adiabatic dynamics in an oligothiophene-fullerene charge transfer complex is employed, relying on a two-state linear vibronic coupling model that combines a distribution of tuning type modes with an intermolecular coordinate that also modulates the electronic coupling. Efficient ML-GMCTDH simulations are carried out for up to 300 vibrational modes using an implementation within the QUANTICS program. Excellent agreement with reference ML-MCTDH calculations is obtained.

Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis

The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.

On the role of non-diagonal system-environment interactions in bridge-mediated electron transfer

Bridge-mediated electron transfer (ET) between a donor and an acceptor is prototypical for the description of numerous most important ET scenarios. While multi-step ET and the interplay of sequential and direct superexchange transfer pathways in the donor-bridge-acceptor (D-B-A) model are increasingly understood, the influence of off-diagonal system-bath interactions on the transfer dynamics is less explored. Off-diagonal interactions account for the dependence of the ET coupling elements on nuclear coordinates (non-Condon effects) and are typically neglected. Here, we numerically investigate with quasi-adiabatic propagator path integral simulations the impact of off-diagonal system-environment interactions on the transfer dynamics for a wide range of scenarios in the D-B-A model. We demonstrate that off-diagonal system-environment interactions can have profound impact on the bridge-mediated ET dynamics. In the considered scenarios, the dynamics itself does not allow for a rigorous assignment of the underlying transfer mechanism. Furthermore, we demonstrate how off-diagonal system-environment interaction mediates anomalous localization by preventing long-time depopulation of the bridge B and how coherent transfer dynamics between donor D and acceptor A can be facilitated. The arising non-exponential short-time dynamics and coherent oscillations are interpreted within an equivalent Hamiltonian representation of a primary reaction coordinate model that reveals how the complex vibronic interplay of vibrational and electronic degrees of freedom underlying the non-Condon effects can impose donor-to-acceptor coherence transfer on short timescales.

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.

Dissipation-Assisted Matrix Product Factorization

Charge and energy transfer in biological and synthetic organic materials are strongly influenced by the coupling of electronic states to a highly structured dissipative environment. Nonperturbative simulations of these systems require a substantial computational effort, and current methods can only be applied to large systems if environmental structures are severely coarse grained. Time evolution methods based on tensor networks are fundamentally limited by the times that can be reached due to the buildup of entanglement in time, which quickly increases the size of the tensor representation, i.e., the bond dimension. In this Letter, we introduce a dissipation-assisted matrix product factorization (DAMPF) method that combines a tensor network representation of the vibronic state within a pseudomode description of the environment where a continuous bosonic environment is mapped into a few harmonic oscillators under Lindblad damping. This framework is particularly suitable for a tensor network representation, since damping suppresses the entanglement growth among oscillators and significantly reduces the bond dimension required to achieve a desired accuracy. We show that dissipation removes the "time-wall" limitation of existing methods, enabling the long-time simulation of large vibronic systems consisting of 10-50 sites coupled to 100-1000 underdamped modes in total and for a wide range of parameter regimes. For these reasons, we believe that our formalism will facilitate the investigation of spatially extended systems with applications to quantum biology, organic photovoltaics, and quantum thermodynamics.

Non-equilibrium relaxation of hot states in organic semiconductors: Impact of mode-selective excitation on charge transfer

The theoretical study of open quantum systems strongly coupled to a vibrational environment remains computationally challenging due to the strongly non-Markovian characteristics of the dynamics. We study this problem in the case of a molecular dimer of the organic semiconductor tetracene, the exciton states of which are strongly coupled to a few hundreds of molecular vibrations. To do so, we employ a previously developed tensor network approach, based on the formalism of matrix product states. By analyzing the entanglement structure of the system wavefunction, we can expand it in a tree tensor network state, which allows us to perform a fully quantum mechanical time evolution of the exciton-vibrational system, including the effect of 156 molecular vibrations. We simulate the dynamics of hot states, i.e., states resulting from excess energy photoexcitation, by constructing various initial bath states, and show that the exciton system indeed has a memory of those initial configurations. In particular, the specific pathway of vibrational relaxation is shown to strongly affect the quantum coherence between exciton states in time scales relevant for the ultrafast dynamics of application-relevant processes such as charge transfer. The preferential excitation of low-frequency modes leads to a limited number of relaxation pathways, thus "protecting" quantum coherence and leading to a significant increase in the charge transfer yield in the dimer structure.

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.

Effect of an underdamped vibration with both diagonal and off-diagonal exciton-phonon interactions on excitation energy transfer

A numerically exact approach, named as the hierarchical stochastic Schrödinger equation, is employed to investigate the resonant vibration-assisted excitation energy transfer in a dimer system, where an underdamped vibration with both diagonal and off-diagonal exciton-phonon interactions is incorporated. From a large parameter space over the site-energy difference, excitonic coupling, and reorganization energy, it is found that the promotion effect of the underdamped vibration is significant only when the excitonic coupling is smaller than the site-energy difference. Under the circumstance, there is an optimal strength ratio between diagonal and off-diagonal exciton-phonon interactions for the resonant vibration-assisted excitation energy transfer as the site-energy difference is greater than the reorganization energy, whereas in the opposite situation the most efficient energy transfer occurs as the exciton-phonon interaction is totally off-diagonal. © 2018 Wiley Periodicals, Inc.

From Fundamental Theories to Quantum Coherences in Electron Transfer

Photoinduced electron transfer (ET) is a cornerstone of energy transduction from light to chemistry. The past decade has seen tremendous advances in the possible role of quantum coherent effects in the light-initiated energy and ET processes in chemical, biological, and materials systems. The prevalence of such coherence effects holds a promise to increase the efficiency and robustness of transport even in the face of energetic or structural disorder. A primary motive of this Perspective is to work out how to think about "coherence" in ET reactions. We will discuss how the interplay of basic parameters governing ET reactions-like electronic coupling, interactions with the environment, and intramolecular high-frequency quantum vibrations-impact coherences. This includes revisiting the insights from the seminal work on the theory of ET and time-resolved measurements on coherent dynamics to explore the role of coherences in ET reactions. We conclude by suggesting that in addition to optical spectroscopies, validating the functional role of coherences would require simultaneous mapping of correlated electron motion and atomically resolved nuclear structure.

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.

Two-Dimensional Spectroscopy Is Being Used to Address Core Scientific Questions in Biology and Materials Science

Two-dimensional spectroscopy is a powerful tool for extracting structural and dynamic information from a wide range of chemical systems. We provide a brief overview of the ways in which two-dimensional visible and infrared spectroscopies are being applied to elucidate fundamental details of important processes in biological and materials science. The topics covered include amyloid proteins, photosynthetic complexes, ion channels, photovoltaics, batteries, as well as a variety of promising new methods in two-dimensional spectroscopy.

Full Quantum Dynamics Simulation of a Realistic Molecular System Using the Adaptive Time-Dependent Density Matrix Renormalization Group Method

The accurate theoretical interpretation of ultrafast time-resolved spectroscopy experiments relies on full quantum dynamics simulations for the investigated system, which is nevertheless computationally prohibitive for realistic molecular systems with a large number of electronic and/or vibrational degrees of freedom. In this work, we propose a unitary transformation approach for realistic vibronic Hamiltonians, which can be coped with using the adaptive time-dependent density matrix renormalization group (t-DMRG) method to efficiently evolve the nonadiabatic dynamics of a large molecular system. We demonstrate the accuracy and efficiency of this approach with an example of simulating the exciton dissociation process within an oligothiophene/fullerene heterojunction, indicating that t-DMRG can be a promising method for full quantum dynamics simulation in large chemical systems. Moreover, it is also shown that the proper vibronic features in the ultrafast electronic process can be obtained by simulating the two-dimensional (2D) electronic spectrum by virtue of the high computational efficiency of the t-DMRG method.

Charge-Transfer Dynamics in the Lowest Excited State of a Pentacene-Fullerene Complex: Implications for Organic Solar Cells

We characterize the dynamic nature of the lowest excited state in a pentacene/C₆₀ complex on the femtosecond time scale, via a combination of ab initio molecular dynamics and time-dependent density functional theory. We analyze the correlations between the molecular vibrations of the complex and the oscillations in the electron-transfer character of its lowest excited state, which point to vibration-induced coherences between the (pentacene-based) local-excitation (LE) state and the complex charge-transfer (CT) state. We discuss the implications of our results on this model system for the exciton-dissociation process in organic solar cells.

Effect of Off-Diagonal Exciton-Phonon Coupling on Intramolecular Singlet Fission

Intramolecular singlet fission (iSF) materials provide remarkable advantages in terms of tunable electronic structures, and quantum chemistry studies have indicated strong electronic coupling modulation by high frequency phonon modes. In this work, we formulate a microscopic model of iSF with simultaneous diagonal and off-diagonal coupling to high-frequency modes. A nonperturbative treatment, the Dirac-Frenkel time-dependent variational approach is adopted using the multiple Davydov trial states. It is shown that both diagonal and off-diagonal coupling can aid efficient singlet fission if excitonic coupling is weak, and fission is only facilitated by diagonal coupling if excitonic coupling is strong. In the presence of off-diagonal coupling, it is found that high frequency modes create additional fission channels for rapid iSF. Results presented here may help provide guiding principles for design of efficient singlet fission materials by directly tuning singlet-triplet interstate coupling.

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.

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.

Photovoltaic concepts inspired by coherence effects in photosynthetic systems

The past decade has seen rapid advances in our understanding of how coherent and vibronic phenomena in biological photosynthetic systems aid in the efficient transport of energy from light-harvesting antennas to photosynthetic reaction centres. Such coherence effects suggest strategies to increase transport lengths even in the presence of structural disorder. Here we explore how these principles could be exploited in making improved solar cells. We investigate in depth the case of organic materials, systems in which energy and charge transport stand to be improved by overcoming challenges that arise from the effects of static and dynamic disorder - structural and energetic - and from inherently strong electron-vibration couplings. We discuss how solar-cell device architectures can evolve to use coherence-exploiting materials, and we speculate as to the prospects for a coherent energy conversion system. We conclude with a survey of the impacts of coherence and bioinspiration on diverse solar-energy harvesting solutions, including artificial photosynthetic systems.

Tracking the coherent generation of polaron pairs in conjugated polymers

The optical excitation of organic semiconductors not only generates charge-neutral electron-hole pairs (excitons), but also charge-separated polaron pairs with high yield. The microscopic mechanisms underlying this charge separation have been debated for many years. Here we use ultrafast two-dimensional electronic spectroscopy to study the dynamics of polaron pair formation in a prototypical polymer thin film on a sub-20-fs time scale. We observe multi-period peak oscillations persisting for up to about 1 ps as distinct signatures of vibronic quantum coherence at room temperature. The measured two-dimensional spectra show pronounced peak splittings revealing that the elementary optical excitations of this polymer are hybridized exciton-polaron-pairs, strongly coupled to a dominant underdamped vibrational mode. Coherent vibronic coupling induces ultrafast polaron pair formation, accelerates the charge separation dynamics and makes it insensitive to disorder. These findings open up new perspectives for tailoring light-to-current conversion in organic materials.

Understanding of charge transfer dynamics is essential to the design of high-performance organic semiconductors for optoelectronic applications. Here, the authors show that excitons, polaron pairs and a long-lived vibrational mode are strongly coupled to each other up to 1 picosecond in polythiophene.

Slow Intramolecular Vibrational Relaxation Leads to Long-Lived Excited-State Wavepackets

Broadband optical pump and compressed white light continuum probe were used to measure the transient excited-state absorption, ground-state bleach, and stimulated emission signals of cresyl violet solution in methanol. Amplitude oscillations caused by wavepacket motion in the ground and excited electronic states were analyzed. It was found that vibrational coherences in the excited state persist for more than the experimental waiting time window of 6 ps, and the strongest mode had a dephasing time constant of 2.4 ps. We hypothesize the dephasing of the wavepacket in the excited state is predominantly caused by intramolecular vibrational relaxation (IVR). Slow IVR indicates weak mode-mode coupling and therefore weak anharmonicity of the potential of this vibration. Thus, the initially prepared vibrational wavepacket in the excited state is not significantly perturbed by nonadiabatic coupling to other electronic states, and hence the diabatic and adiabatic representations of the system are essentially identical within the Born-Oppenheimer approximation. The wavepacket therefore evolves with time in an almost harmonic potential, slowly dephased by IVR and the pure vibrational decoherence. The consistency in the position of node (phase change in the wavepacket) in the excited-state absorption and stimulated emission signals without undergoing any frequency shift until the wavepacket is completely dephased conforms to the absence of any reactive internal conversion.

Photon-absorbing charge-bridging states in organic bulk heterojunctions consisting of diketopyrrolopyrrole derivatives and PCBM

We have investigated the photo- and electrochemical properties of five diketopyrrolopyrrole (DPP) derivatives both experimentally and theoretically. In the experimental study, we found that a blend of a DPP derivative named D2 and phenyl-C61-butyric acid methyl ester (PCBM) exhibits the highest internal quantum efficiency (IQE) and power convergence efficiency (PCE) among the five derivatives investigated. In the theoretical study, we found that the open-circuit voltage can be estimated from the difference between the energy gap of frontier orbitals and the voltage loss and that the latter is suppressed when the IQE is large. Then, to investigate the factors that influence the IQE, investigations on charge recombination, hole transfer, and charge transfer induced by photoabsorption were conducted for the complexes of each DPP derivative and PCBM. It was found that D2/PCBM exhibits the largest charge-bridging upon photoabsorption, which leads to the highest IQE and PCE among the five DPP derivatives.

Full-dimensional multilayer multiconfigurational time-dependent Hartree study of electron transfer dynamics in the anthracene/C60 complex

Electron transfer at the donor-acceptor heterojunctions plays a critical role in the photoinduced process during the solar energy conversion in organic photovoltaic materials. We theoretically investigate the electron transfer process in the anthracene/C60 donor-acceptor complex by using quantum dynamics calculations. The electron-transfer model Hamiltonian with full dimensionality was built by quantum-chemical calculations. The quantum dynamics calculations were performed using the multiconfigurational time-dependent Hartree (MCTDH) theory and multilayer (ML) MCTDH methods. The latter approach (ML-MCTDH) allows us to conduct the comprehensive study on the quantum evolution of the full-dimensional electron-transfer model including 4 electronic states and 246 vibrational degrees of freedom. Our quantum dynamics calculations exhibit the ultrafast anthracene → C60 charge transfer process because of the strong coupling between excitonic and charge transfer states. This work demonstrates that the ML-MCTDH is a very powerful method to treat the quantum evolution of complex systems.

Electron transfer within a reaction path model calibrated by constrained DFT calculations: application to mixed-valence organic compounds

The quantum dynamics of electron transfer in mixed-valence organic compounds is investigated using a reaction path model calibrated by constrained density functional theory (cDFT). Constrained DFT is used to define diabatic states relevant for describing the electron transfer, to obtain equilibrium structures for each of these states and to estimate the electronic coupling between them. The harmonic analysis at the diabatic minima yields normal modes forming the dissipative bath coupled to the electronic states. In order to decrease the system-bath coupling, an effective one dimensional vibronic Hamiltonian is constructed by partitioning the modes into a linear reaction path which connects both equilibrium positions and a set of secondary vibrational modes, coupled to this reaction coordinate. Using this vibronic model Hamiltonian, dissipative quantum dynamics is carried out using Redfield theory, based on a spectral density which is determined from the cDFT results. In a first benchmark case, the model is applied to a series of mixed-valence organic compounds formed by two 1,4-dimethoxy-3-methylphenylene fragments linked by an increasing number of phenylene bridges. This allows us to examine the coherent electron transfer in extreme situations leading to a ground adiabatic state with or without a barrier and therefore to the trapping of the charge or to an easy delocalization.