Hot and Cold Charge-Transfer Mechanisms in Organic Photovoltaics: Insights into the Excited States of Donor/Acceptor Interfaces

Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications

Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

Nanoscale and Real-Time Nuclear-Electronic Dynamics Simulation Study of Charge Transfer at the Donor-Acceptor Interface in Organic Photovoltaics

Charge-transfer (CT) processes in donor-acceptor interfaces of organic photovoltaics have been challenging targets for computational chemistry owing to their nanoscale and ultrafast nature. Herein, we report real-time nuclear-electronic dynamics simulations of CT processes in a nanometer-scale donor-acceptor interface model composed of a donor poly(3-hexylthiophene-2,5-diyl) crystal and an acceptor [6,6]-phenyl-C₆₁-butyric acid methyl ester aggregate. The simulations were realized using our original reduced-scaling computational technique, namely, patchwork-approximation-based Ehrenfest dynamics. The results illustrated the CT pathway with atomic resolution, thereby rationalizing the observed excitation-energy dependence of the quantity of CT. Further, nuclear motion, which is affected by the electronic dynamics, was observed to play a significant role in the CT process by modulating molecular orbital energies. The present study suggests that microscopic CT processes strongly depend on local structures of disordered donor-acceptor interfaces as well as coupling between nuclear and electronic dynamics.

Solvent effects on excited-state relaxation dynamics of paddle-wheel BODIPY-Hexaoxatriphenylene conjugates: Insights from non-adiabatic dynamics simulations

Understanding the excited state dynamics of donor-acceptor (D-A) complexes is of fundamental importance both experimentally and theoretically. Herein, we have first explored the photoinduced dynamics of a recently synthesized paddle-wheel BODIPY-hexaoxatriphenylene (BODIPY is the abbreviation for BF2-chelated dipyrromethenes) conjugates D-A complexes with the combination of both electronic structure calculations and non-adiabatic dynamics simulations. On the basis of computational results, we concluded that the BODIPY-hexaoxatriphenylene (BH) conjugates will be promoted to the local excited (LE) states of the BODIPY fragments upon excitation, which is followed by the ultrafast exciton transfer from LE state to charge transfer (CT). Instead of the photoinduced electron transfer process proposed in previous experimental work, such a exciton transfer process is accompanied with the photoinduced hole transfer from BODIPY to hexaoxatriphenylene. Additionally, solvent effects are found to play an important role in the photoinduced dynamics. Specifically, the hole transfer dynamics is accelerated by the acetonitrile solvent, which can be ascribed to significant influences of the solvents on the charge transfer states, i.e. the energy gaps between LE and CT excitons are reduced greatly and the non-adiabatic couplings are increased in the meantime. Our present work not only provides valuable insights into the underlying photoinduced mechanism of BH, but also can be helpful for the future design of novel donor-acceptor conjugates with better optoelectronic performance.

Free Charge Carriers in Homo-Sorted π-Stacks of Donor-Acceptor Conjugates

Inspired by the high photoconversion efficiency observed in natural light-harvesting systems, the hierarchical organization of molecular building blocks has gained impetus in the past few decades. Particularly, the molecular arrangement and packing in the active layer of organic solar cells (OSCs) have garnered significant attention due to the decisive role of the nature of donor/acceptor (D/A) heterojunctions in charge carrier generation and ultimately the power conversion efficiency. This review focuses on the recent developments in emergent optoelectronic properties exhibited by self-sorted donor-on-donor/acceptor-on-acceptor arrangement of covalently linked D-A systems, highlighting the ultrafast excited state dynamics of charge transfer and transport. Segregated organization of donors and acceptors promotes the delocalization of photoinduced charges among the stacks, engendering an enhanced charge separation lifetime and percolation pathways with ambipolar conductivity and charge carrier yield. Covalently linking donors and acceptors ensure a sufficient D-A interface and interchromophoric electronic coupling as required for faster charge separation while providing better control over their supramolecular assemblies. The design strategies to attain D-A conjugate assemblies with optimal charge carrier generation efficiency, the scope of their application compared to state-of-the-art OSCs, current challenges, and future opportunities are discussed in the review. An integrated overview of rational design approaches derived from the comprehension of underlying photoinduced processes can pave the way toward superior optoelectronic devices and bring in new possibilities to the avenue of functional supramolecular architectures.

Chemical Bonding as a New Avenue for Controlling Excited-State Properties and Excitation Energy-Transfer Processes in Zinc Phthalocyanine-Fullerene Dyads

Whether chemical bonding can regulate the excited-state and optoelectronic properties of donor-acceptor dyads has been largely elusive. In this work, we used electronic structure and nonadiabatic dynamics methods to explore the excited-state properties of covalently bonded zinc phthalocyanine (ZnPc)-fullerene (C₆₀ ) dyads with a 6-6 (or 5-6) bonding configuration in which ZnPc is bonded to two carbon atoms shared by the two hexagonal rings (or a pentagonal and a hexagonal ring) in C₆₀ . In both cases, the locally excited (LE) states on ZnPc are spectroscopically bright. However, their different chemical bonding differentiates the electronic interactions between ZnPc and C₆₀ . In the 5-6 bonding configuration, the LE states on ZnPc are much higher in energy than the LE states on C₆₀ . Thus, the excitation energy transfer from ZnPc to C₆₀ is thermodynamically favorable. On the other hand, in the 6-6 bonding configuration, such a process is inhibited because the LE states on ZnPc are the lowest ones. More detailed mechanisms are elucidated from nonadiabatic dynamics simulations. In the 6-6 bonding configuration, no excitation energy transfer was observed. In contrast, in the 5-6 bonding configuration, several LE and charge-transfer (CT) excitons were shown to participate in the energy-transfer process. Further analysis reveals that the photoinduced energy transfer is mediated by a CT exciton, such that electron- and hole-transfer processes take place in a concerted but asynchronous manner in the excitation energy transfer. It is also found that high-level electronic structure methods including exciton effects are indispensable to accurately describe photoinduced energy- and electron-transfer processes. Furthermore, this work opens up new avenues for regulating the excited-state properties of molecular donor-acceptor dyads by means of chemical bonding.

Ultrafast channel I and channel II charge generation processes at a nonfullerene donor-acceptor PTB7:PDI interface is crucial for its excellent photovoltaic performance

Nonfullerene organic solar cells have received much attention in recent years due to their low cost, high absorption coefficient and excellent synthetic flexibility. However, the microscopic photoinduced dynamics at corresponding donor-acceptor interfaces remains unclear. In this work, we have firstly employed state-of-the-art TDDFT-based nonadiabatic dynamics simulations in combination with static electronic structure calculations to explore the ultrafast photoinduced dynamics at a typical nonfullerene donor-acceptor PTB7:PDI interface using a minimal model system (172 atoms). Upon excitation with specific wavelength of light, both PTB7 and PDI can be locally excited to generate |PTB7* and |PDI* excitons due to their high absorption ability and significant overlap in absorption spectrum. After that, these localized excitons gradually convert to charge transfer exciton |PTB7+PDI-, while another |PTB7-PDI+ charge transfer exciton is not involved in the whole process. Along with the exciton conversion, electron transfer from PTB7 to PDI (channel I charge generation) and the hole transfer from PDI to PTB7 (channel II charge generation) occurs simultaneously with time constants of 643 fs and 549 fs respectively. In the same time, D index that measures the centroid distance of electron and hole increases from 1.0 Å to 4.0 Å, which clearly reflects a charge transfer process at the interface. Our present work provides solid evidence that both channel I and channel II charge generation processes play important roles at PTB7:PDI interface, which could be helpful for the design of novel nonfullerene solar cells with better photovoltaic performance.

Charge-transfer electronic states in organic solar cells: a TDDFT study

The prediction of new organic photovoltaic materials in organic solar cells (OSCs) must include a precise description of charge-transfer states because they are involved in electron-transfer processes such as charge separation and charge recombination which govern the device efficiency. Also, as the experimental performance of an optoelectronic device is measured for nonequilibrium nanostructures, computational approaches need models that can incorporate morphology effects. Usually, this aspect is treated by molecular dynamics simulation (MDS) methodologies; however, methodologies and formalisms to calculate the electron-transfer processes are still controversial and sometimes do not connect their information with the phase morphologies. In this work we propose a simple and fast characterization of electron-transfer processes to find the rate constants by analysing the distribution of vertical excitation energies of both local excitation (LE) and charge-transfer (CT) states using TD-DFT calculations in the donor-acceptor pair structures which were extracted from MDS. This proposal assumes that conformational changes are prevented and equilibria are not achieved while the electron-transfer events take effect, and thus the only pathway that connects the LE and CT states is their surface crossing point where an ideal distribution might exist. Different density functionals and dialectric models were tested. The results indicate a close relationship between the proposal and experimental data for electron-transfer events, suggesting the application of this method in the rational design of new photovoltaic materials.

Chirality of a rhodamine heterodimer linked to a DNA scaffold: an experimental and computational study

The chiroptical properties of multi-chromophoric systems are governed by the intermolecular arrangement of the monomeric units. We report on a computational and experimental study of the linear optical properties and supramolecular structure of a rhodamine heterodimer assembled on a DNA scaffold. The experimental absorption and circular dichroism (CD) profiles confirm the dimer formation. Computationally, starting from low-cost DFT/TDDFT simulations of the bare dimer we attribute the measured -/+ CD sign sequence of the S1/S2 bands to a specific chiral conformation of the heterodimer. In the monomers, as typical for rhodamine dyes, the electric transition dipole of the lowest π-π* transition is parallel to the long axis of the xanthene planes. We show that in the heterodimer the sign sequence of the two CD bands is related to the orientation of these long axes. To account explicitly for environment effects, we use molecular dynamics (MD) simulations for characterizing the supramolecular structure of the two optical isomers tethered on DNA. Average absorption and CD-profiles were modeled using ab initio TDDFT calculations at the geometries sampled along a few nanosecond MD run. The absorption profiles computed for both optical isomers are in good agreement with the experimental absorption spectrum and do not allow one to discriminate between them. The computed averaged CD profiles provide the orientation of monomers in the enantiomer that is dominant under the experimental conditions.

Local Excitation/Charge-Transfer Hybridization Simultaneously Promotes Charge Generation and Reduces Nonradiative Voltage Loss in Nonfullerene Organic Solar Cells

High power conversion efficiencies in state-of-the-art nonfullerene organic solar cells (NF OSCs) call for elucidation of the underlying working mechanisms of both high photocurrent densities and low nonradiative voltage losses under small energy offsets. Here, to address this fundamental issue, we have assessed the nature of interfacial charge-transfer (CT) states in a representative small-molecule NF OSC (DRTB-T:IT-4F) by time-dependent density functional theory calculations. The calculated results point to the fact that the CT states can borrow considerable oscillator strengths from the energy-close local excitation (LE) states or be fully hybridized with these LE states by molecular aggregation at the donor-acceptor interfaces. The LE/CT hybridization can promote charge generation by direct population of thermalized CT or LE/CT states under illumination. At the same time, the increased oscillator strengths of the lowest CT state will improve the luminescence quantum efficiencies and thus reduce nonradiative voltage losses. Our work suggests that it is crucial to tune the LE/CT hybridization by optimization of the donor and acceptor molecular and interfacial structures to further improve the NF OSC performance.

A novel empirical method for quickly estimating the charge-transfer state of fullerene-donor derivatives

A new convenient equivalent conjugated electron number technique is put forward to predict the charge-transfer state in fullerene-donor derivatives.

Optimum driving energy for achieving balanced open-circuit voltage and short-circuit current density in organic bulk heterojunction solar cells

Organic bulk heterojunction solar cells generally suffer from a trade-off between the open circuit voltage (Voc) and the short circuit current density (Jsc) under a given donor/acceptor (D/A) interfacial energetic offset (or the so-called driving force). Here we theoretically investigate the optimum driving energy required for achieving the balanced Jsc and Voc simultaneously. To this end, the Jscversus the driving force ΔE curves are calculated under two different charge separation mechanisms by employing the drift-diffusion method. For the Marcus incoherent mechanism, the curve features a high plateau in a broad range of ΔE starting from 0.2 eV, which is due to the accumulation of undissociated excitons within their lifetime and signifies the possibility of obtaining a sizable Jsc under a ΔE value much smaller than the reorganization energy. After incorporating both the electron and hole transfer pathways into the device model, the calculated J-V curves are comparable to experimentally measured ones foractual blended systems of different driving forces. For the coherent mechanism, it is demonstrated that the maximum Jsc can also be achieved under the ΔE of 0.2 eV if a large proportion of the high-lying delocalized states are harvested through tuning the density of states for the charge transfer excitons to reduce the sub-gap states. This theoretical work revealed quantitatively the relationship between the interfacial energy offsets and device performance, and also provides some guidelines for identifying the macroscopic features of the actual charge separation mechanisms in bulk heterojunction solar cells.

Effects of fluorination and thermal annealing on charge recombination processes in polymer bulk-heterojunction solar cells

We investigate the effect of fluorination on the photovoltaic properties of an alternating conjugated polymer composed of 4,8-di-2-thienylbenzo[1,2-b:4,5-b']dithiophene (BDT) and 4,7-bis([2,2'-bithiophen]-5-yl)-benzo-2-1-3-thiadiazole (4TBT) units in bulk-heterojunction solar cells. The unsubstituted and fluorinated polymers afford very similar open-circuit voltages and fill factor values, but the fluorinated polymer performed better due to enhanced aggregation which provides a higher photocurrent. The photovoltaic performance of both materials improved upon thermal annealing at 150-200 °C as a result of a significantly increased fill factor and open-circuit voltage, counteracted by a slight loss in photocurrent. Detailed studies of the morphology, light intensity dependence, external quantum efficiency and electroluminescence allowed the exploration of the effects of fluorination and thermal annealing on the charge recombination and the nature of the donor-acceptor interfacial charge transfer states in these films.

Dielectric Screening Meets Optimally Tuned Density Functionals

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.

Theoretical study on mesoscopic-size impurity effects in the charge separation process of organic photocells

A bulk-heterojunction structure is often employed to develop high-performance organic photocells, in which the donor and acceptor regions are complexly intertwined.