Quantum localization and delocalization of charge carriers in organic semiconducting crystals

Exciton Transport in the Nonfullerene Acceptor O-IDTBR from Nonadiabatic Molecular Dynamics

Theory, computation, and experiment have given strong evidence that charge carriers in organic molecular crystals form partially delocalized quantum objects that diffuse very efficiently via a mechanism termed transient delocalization. It is currently unclear how prevalent this mechanism is for exciton transport. Here we carry out simulation of singlet Frenkel excitons (FE) in a molecular organic semiconductor that belongs to the class of nonfullerene acceptors, O-IDTBR, using the recently introduced FE surface hopping nonadiabatic molecular dynamics method. We find that FE are, on average, localized on a single molecule in the crystal due to sizable reorganization energy and moderate excitonic couplings. Yet, our simulations suggest that the diffusion mechanism is more complex than simple local hopping; in addition to hopping, we observe frequent transient delocalization events where the exciton wave function expands over 10 or more molecules for a short period of time in response to thermal excitations within the excitonic band, followed by de-excitation and contraction onto a single molecule. The transient delocalization events lead to an increase in the diffusion constant by a factor of 3-4, depending on the crystallographic direction as compared to the situation where only local hopping events are considered. Intriguingly, O-IDTBR appears to be a moderately anisotropic 3D "conductor" for excitons but a highly anisotropic 2D conductor for electrons. Taken together with previous simulation results, two trends seem to emerge for molecular organic crystals: excitons tend to be more localized and slower than charge carriers due to higher internal reorganization energy, while exciton transport tends to be more isotropic than charge transport due to the weaker distance dependence of excitonic versus electronic coupling.

Non-adiabatic molecular dynamics simulations provide new insights into the exciton transfer in the Fenna-Matthews-Olson complex

A trajectory surface hopping approach, which uses machine learning to speed up the most time-consuming steps, has been adopted to investigate the exciton transfer in light-harvesting systems. The present neural networks achieve high accuracy in predicting both Coulomb couplings and excitation energies. The latter are predicted taking into account the environment of the pigments. Direct simulation of exciton dynamics through light-harvesting complexes on significant time scales is usually challenging due to the coupled motion of nuclear and electronic degrees of freedom in these rather large systems containing several relatively large pigments. In the present approach, however, we are able to evaluate a statistically significant number of non-adiabatic molecular dynamics trajectories with respect to exciton delocalization and exciton paths. The formalism is applied to the Fenna-Matthews-Olson complex of green sulfur bacteria, which transfers energy from the light-harvesting chlorosome to the reaction center with astonishing efficiency. The system has been studied experimentally and theoretically for decades. In total, we were able to simulate non-adiabatically more than 30 ns, sampling also the relevant space of parameters within their uncertainty. Our simulations show that the driving force supplied by the energy landscape resulting from electrostatic tuning is sufficient to funnel the energy towards site 3, from where it can be transferred to the reaction center. However, the high efficiency of transfer within a picosecond timescale can be attributed to the rather unusual properties of the BChl a molecules, resulting in very low inner and outer-sphere reorganization energies, not matched by any other organic molecule, e.g., used in organic electronics. A comparison with electron and exciton transfer in organic materials is particularly illuminating, suggesting a mechanism to explain the comparably high transfer efficiency.

Extending non-adiabatic rate theory to strong electronic couplings in the Marcus inverted regime

Electron transfer reactions play an essential role in many chemical and biological processes. Fermi's golden rule (GR), which assumes that the coupling between electronic states is small, has formed the foundation of electron transfer rate theory; however, in short range electron/energy transfer reactions, this coupling can become very large, and, therefore, Fermi's GR fails to make even qualitatively accurate rate predictions. In this paper, I present a simple modified GR theory to describe electron transfer in the Marcus inverted regime at arbitrarily large electronic coupling strengths. This theory is based on an optimal global rotation of the diabatic states, which makes it compatible with existing methods for calculating GR rates that can account for nuclear quantum effects with anharmonic potentials. Furthermore, the optimal GR (OGR) theory can also be combined with analytic theories for non-adiabatic rates, such as Marcus theory and Marcus-Levich-Jortner theory, offering clear physical insights into strong electronic coupling effects in non-adiabatic processes. OGR theory is also tested on a large set of spin-boson models and an anharmonic model against exact quantum dynamics calculations, where it performs well, correctly predicting rate turnover at large coupling strengths. Finally, an example application to a boron-dipyrromethane-anthracene photosensitizer reveals that strong coupling effects inhibit excited state charge recombination in this system, reducing the rate of this process by a factor of 4. Overall, OGR theory offers a new approach to calculating electron transfer rates at strong couplings, offering new physical insights into a range of non-adiabatic processes.

Transient Localization from the Interaction with Quantum Bosons

We carefully revisit the electron-boson scattering problem, going beyond weak-coupling expansions and popular semiclassical treatments. By providing numerically exact results valid at finite temperatures, we demonstrate the existence of a broad regime of electron-boson scattering where quantum localization processes become relevant despite the absence of extrinsic disorder. Localization in the Anderson sense is caused by the dynamical randomness resulting from a large thermal boson population, being, however, effective only at transient times before diffusion can set in. Compelling evidence of this transient localization phenomenon is provided by the observation of a distinctive displaced Drude peak in the optical absorption and the ensuing suppression of conductivity. Our findings identify a general route for anomalous metallic behavior that can broadly apply in interacting quantum matter.

Accessing the electronic structure of liquid crystalline semiconductors with bottom-up electronic coarse-graining

Understanding the relationship between multiscale morphology and electronic structure is a grand challenge for semiconducting soft materials. Computational studies aimed at characterizing these relationships require the complex integration of quantum-chemical (QC) calculations, all-atom and coarse-grained (CG) molecular dynamics simulations, and back-mapping approaches. However, these methods pose substantial computational challenges that limit their application to the requisite length scales of soft material morphologies. Here, we demonstrate the bottom-up electronic coarse-graining (ECG) of morphology-dependent electronic structure in the liquid-crystal-forming semiconductor, 2-(4-methoxyphenyl)-7-octyl-benzothienobenzothiophene (BTBT). ECG is applied to construct density functional theory (DFT)-accurate valence band Hamiltonians of the isotropic and smectic liquid crystal (LC) phases using only the CG representation of BTBT. By bypassing the atomistic resolution and its prohibitive computational costs, ECG enables the first calculations of the morphology dependence of the electronic structure of charge carriers across LC phases at the ∼20 nm length scale, with robust statistical sampling. Kinetic Monte Carlo (kMC) simulations reveal a strong morphology dependence on zero-field charge mobility among different LC phases as well as the presence of two-molecule charge carriers that act as traps and hinder charge transport. We leverage these results to further evaluate the feasibility of developing mesoscopic, field-based ECG models in future works. The fully CG approach to electronic property predictions in LC semiconductors opens a new computational direction for designing electronic processes in soft materials at their characteristic length scales.

The organo-metal-like nature of long-range conduction in cable bacteria

Cable bacteria are filamentous, multicellular microorganisms that display an exceptional form of biological electron transport across centimeter-scale distances. Currents are guided through a network of nickel-containing protein fibers within the cell envelope. Still, the mechanism of long-range conduction remains unresolved. Here, we characterize the conductance of the fiber network under dry and wet, physiologically relevant, conditions. Our data reveal that the fiber conductivity is high (median value: 27 S cm-1; range: 2 to 564 S cm-1), does not show any redox signature, has a low thermal activation energy (Ea = 69 ± 23 meV), and is not affected by humidity or the presence of ions. These features set the nickel-based conduction mechanism in cable bacteria apart from other known forms of biological electron transport. As such, conduction resembles that of an organic semi-metal with a high charge carrier density. Our observation that biochemistry can synthesize an organo-metal-like structure opens the way for novel bio-based electronic technologies.

MesoHOPS: Size-invariant scaling calculations of multi-excitation open quantum systems

The photoexcitation dynamics of molecular materials on the 10-100 nm length scale depend on complex interactions between electronic and vibrational degrees of freedom, rendering exact calculations difficult or intractable. The adaptive Hierarchy of Pure States (adHOPS) is a formally exact method that leverages the locality imposed by interactions between thermal environments and electronic excitations to achieve size-invariant scaling calculations for single-excitation processes in systems described by a Frenkel-Holstein Hamiltonian. Here, we extend adHOPS to account for arbitrary couplings between thermal environments and vertical excitation energies, enabling formally exact, size-invariant calculations that involve multiple excitations or states with shared thermal environments. In addition, we introduce a low-temperature correction and an effective integration of the noise to reduce the computational expense of including ultrafast vibrational relaxation in Hierarchy of Pure States (HOPS) simulations. We present these advances in the latest version of the open-source MesoHOPS library and use MesoHOPS to characterize charge separation at a one-dimensional organic heterojunction when both the electron and hole are mobile.

Unified Entropy-Ruled Einstein's Relation for Bulk and Low-Dimensional Molecular-Material Systems: A Hopping-to-Band Shift Paradigm

We present a unified paradigm on entropy-ruled Einstein's diffusion-mobility relation (μ/D ratio) for 1D, 2D, and 3D free-electron solid state systems. The localization transport in the extended molecules is well approximated by the continuum time-delayed hopping factor within our unified entropy-ruled transport method of noninteracting quantum systems. Moreover, we generalize an entropy-dependent diffusion relation for 1D, 2D, and 3D systems as defined by D d , h e f f = D d , h e f f = 0 ⁡ exp ( ( d - 1 ) h e f f d + 2 ) , where heff and d are the effective entropy and dimension (d = 1, 2, 3), respectively. This generalized relation is valid for both equilibrium and nonequilibrium transport systems since the parameter heff is closely connected with the nonequilibrium fluctuation theorem-based entropy production rule. Importantly, we herein revisit the Boltzmann approach using an entropy-ruled method for mobility calculation for the universal quantum materials that is expressed as μ d = [ ( d d + 2 ) q d h e f f d η ] v F 2 τ 2 , where v F 2 τ 2 is the diffusion constant for band transport systems and η is the chemical potential. According to our entropy-ruled μ/D relation, the Navamani-Shockley diode equation is transformed.

On the application of Marcus-Hush theory to small polaron chemical dynamics in oxides: its relationship to the Holstein model and the importance of lattice-orbital symmetries

The chemical dynamics of small polaron hopping within oxides is often interpreted through two-site variations on Marcus-Hush theory, while from a physics perspective small polaron hopping is more often approached from Holstein's solid-state formalism. Here we seek to provide a chemically oriented viewpoint, focusing on small polaron hopping in oxides, concerning these two phenomenological frameworks by employing both tight-binding modelling and first-principles calculations. First, within a semiclassical approach the Marcus-Hush relations are overviewed as a two-site reduction of Holstein's model. Within the single-band regime, similarities and differences between Holstein derived small polaron hopping and the Marcus-Hush model are also discussed. In this context the emergence of adiabaticity (or, conversely, diabaticity) is also explored within each framework both analytically and by directly evolving the system wavefunction. Then, through first-principles calculations of select oxides we explore how coupled lattice and orbital symmetries can impact on hopping properties - in a manner that is quite distinct typical chemical applications of Marcus-Hush theory. These results are then related back to the Holstein model to explore the relative applicability of the two frameworks towards interpreting small polaron hopping properties, where it is emphasized that the Holstein model offers an increasingly more appealing physicochemical interpretation of hopping processes as band and/or coupling interactions increase. Overall, this work aims to strengthen the physically oriented exploration of small polarons and their physicochemical properties in the growing oxide chemistry community.

A model analysis of centimeter-long electron transport in cable bacteria

The recent discovery of cable bacteria has greatly expanded the known length scale of biological electron transport, as these multi-cellular bacteria are capable of mediating electrical currents across centimeter-scale distances. To enable such long-range conduction, cable bacteria embed a network of regularly spaced, parallel protein fibers in their cell envelope. These fibers exhibit extraordinary electrical properties for a biological material, including an electrical conductivity that can exceed 100 S cm-1. Traditionally, long-range electron transport through proteins is described as a multi-step hopping process, in which the individual hopping steps are described by Marcus electron transport theory. Here, we investigate to what extent such a classical hopping model can explain the conductance data recorded for individual cable bacterium filaments. To this end, the conductive fiber network in cable bacteria is modelled as a set of parallel one-dimensional hopping chains. Comparison of model simulated and experimental current(I)/voltage(V) curves, reveals that the charge transport is field-driven rather than concentration-driven, and there is no significant injection barrier between electrodes and filaments. However, the observed high conductivity levels (>100 S cm-1) can only be reproduced, if we include much longer hopping distances (a > 10 nm) and lower reorganisation energies (λ < 0.2 eV) than conventionally used in electron relay models of protein structures. Overall, our model analysis suggests that the conduction mechanism in cable bacteria is markedly distinct from other known forms of long-range biological electron transport, such as in multi-heme cytochromes.

Spin relaxation of electron and hole polarons in ambipolar conjugated polymers

The charge-transport properties of conjugated polymers have been studied extensively for opto-electronic device applications. Some polymer semiconductors not only support the ambipolar transport of electrons and holes, but do so with comparable carrier mobilities. This opens the possibility of gaining deeper insight into the charge-transport physics of these complex materials via comparison between electron and hole dynamics while keeping other factors, such as polymer microstructure, equal. Here, we use field-induced electron spin resonance spectroscopy to compare the spin relaxation behavior of electron and hole polarons in three ambipolar conjugated polymers. Our experiments show unique relaxation regimes as a function of temperature for electrons and holes, whereby at lower temperatures electrons relax slower than holes, but at higher temperatures, in the so-called spin-shuttling regime, the trend is reversed. On the basis of theoretical simulations, we attribute this to differences in the delocalization of electron and hole wavefunctions and show that spin relaxation in the spin shuttling regimes provides a sensitive probe of the intimate coupling between charge and structural dynamics.

Interlayer Charge Transport in 2D Lead Halide Perovskites from First Principles

We report on the implementation of a versatile projection-operator diabatization approach to calculate electronic coupling integrals in layered periodic systems. The approach is applied to model charge transport across the saturated organic spacers in two-dimensional (2D) lead halide perovskites. The calculations yield out-of-plane charge transfer rates that decay exponentially with the increasing length of the alkyl chain, range from a few nanoseconds to milliseconds, and are supportive of a hopping mechanism. Most importantly, we show that the charge carriers strongly couple to distortions of the Pb-I framework and that accounting for the associated nonlocal dynamic disorder increases the thermally averaged interlayer rates by a few orders of magnitude compared to the frozen-ion 0 K-optimized structure. Our formalism provides the first comprehensive insight into the role of the organic spacer cation on vertical transport in 2D lead halide perovskites and can be readily extended to functional π-conjugated spacers, where we expect the improved energy alignment with the inorganic layout to speed up the charge transfer between the semiconducting layers.

Stabilized coupled trajectory mixed quantum-classical algorithm with improved energy conservation: CTMQC-EDI

Coupled trajectory mixed quantum-classical (CTMQC) dynamics is a rigorous approach to trajectory-based non-adiabatic dynamics, which has recently seen an improvement to energy conservation via the introduction of the CTMQC-E algorithm. Despite this, the method's two key quantities distinguishing it from Ehrenfest dynamics, the modified Born-Oppenheimer momentum and the quantum momentum, require regularization procedures in certain circumstances. Such procedures in the latter can cause instabilities, leading to undesirable effects, such as energy drift and spurious population transfer, which is expected to become increasingly prevalent when the system gets larger as such events would happen more frequently. We propose a further modification to CTMQC-E, which includes a redefinition of the quantum momentum, CTMQC-EDI (double intercept), such that it has no formal divergences. We then show for Tully models I-III and the double arch model that the algorithm has greatly improved total energy conservation and negligible spurious population transfer at all times, in particular in regions of strong non-adiabatic coupling. CTMQC-EDI, therefore, shows promise as a numerically robust non-adiabatic dynamics technique that accounts for decoherence from first principles and that is scalable to large molecular systems and materials.

Large Coherent States Formed from Disordered k-Regular Random Graphs

The present work is motivated by the need for robust, large-scale coherent states that can play possible roles as quantum resources. A challenge is that large, complex systems tend to be fragile. However, emergent phenomena in classical systems tend to become more robust with scale. Do these classical systems inspire ways to think about robust quantum networks? This question is studied by characterizing the complex quantum states produced by mapping interactions between a set of qubits from structure in graphs. We focus on maps based on k-regular random graphs where many edges were randomly deleted. We ask how many edge deletions can be tolerated. Surprisingly, it was found that the emergent coherent state characteristic of these graphs was robust to a substantial number of edge deletions. The analysis considers the possible role of the expander property of k-regular random graphs.

Transiently delocalized states enhance hole mobility in organic molecular semiconductors

Evidence shows that charge carriers in organic semiconductors self-localize because of dynamic disorder. Nevertheless, some organic semiconductors feature reduced mobility at increasing temperature, a hallmark for delocalized band transport. Here we present the temperature-dependent mobility in two record-mobility organic semiconductors: dinaphtho[2,3-b:2',3'-f]thieno[3,2-b]-thiophene (DNTT) and its alkylated derivative, C8-DNTT-C8. By combining terahertz photoconductivity measurements with atomistic non-adiabatic molecular dynamics simulations, we show that while both crystals display a power-law decrease of the mobility (μ) with temperature (T) following μ ∝ T -n, the exponent n differs substantially. Modelling reveals that the differences between the two chemically similar semiconductors can be traced to the delocalization of the different states that are thermally accessible by charge carriers, which in turn depends on their specific electronic band structure. The emerging picture is that of holes surfing on a dynamic manifold of vibrationally dressed extended states with a temperature-dependent mobility that provides a sensitive fingerprint for the underlying density of states.

Charge carrier dynamics in conducting polymer PEDOT using ab initio molecular dynamics simulations

As conducting polymers become increasingly important in electronic devices, understanding their charge transport is essential for material and device development. Various semi-empirical approaches have been used to describe temporal charge carrier dynamics in these materials, but there have yet to be any theoretical approaches utilizing ab initio molecular dynamics. In this work, we develop a computational technique based on ab initio Car-Parrinello molecular dynamics to trace charge carrier temporal motion in archetypical conducting polymer poly(3,4-ethylenedioxythiophene) (PEDOT). Particularly, we analyze charge dynamics in a single PEDOT chain and in two coupled chains with different degrees of coupling and study the effect of temperature. In our model we first initiate a positively charged polaron (compensated by a negative counterion) at one end of the chain, and subsequently displace the counterion to the other end of the chain and trace polaron dynamics in the system by monitoring bond length alternation in the PEDOT backbone and charge density distribution. We find that at low temperature (T = 1 K) the polaron distortion gradually disappears from its initial location and reappears near the new position of the counterion. At the room temperature (T = 300 K), we find that the distortions induced by polaron, and atomic vibrations are of the same magnitude, which makes tracking the polaron distortion challenging because it is hidden behind the temperature-induced vibrations. The novel approach developed in this work can be used to study polaron mobility along and between the chains, investigate charge transport in highly doped polymers, and explore other flexible polymers, including n-doped ones.

Density-functional theory (DFT) and time-dependent DFT study of the chemical and physical origins of key photoproperties of end-group derivatives of a nonfullerene acceptor molecule for bulk heterojunction organic solar cells

As emphasized in a recent review article (Chem. Rev. 2022, 122, 14180), organic solar cell (OSC) photoconversion efficiency has been rapidly evolving with results increasingly comparable to those of traditional inorganic solar cells. Historically, OSC performance improvement focused first on the morphology of P3HT:PC 61 BM solar cells then went through different stages to shift lately interest towards nonfullerene acceptors (NFAs) as a replacement ofPC 61 BM acceptor (ACC) molecule. Here, we use density-functional theory (DFT) and time-dependent DFT to investigate four novel NFAs of A-D-A (acceptor-donor-acceptor) form derived from the recently synthesized IDIC-4Cl (Dyes Pigm. 2019, 166, 196). Our level of theory is carefully evaluated for IDIC-4Cl and then applied to the four novel NFAs in order to understand how chemical modifications lead to physical changes in cyclic voltammetry (CV) frontier molecular orbital energies and absorption spectra in solution. Finally we design and apply a new type of Scharber plot for NFAs based upon some simple but we think reasonable assumptions. Unlike the original Scharber plots where a larger DON band gap favors a larger PCE, our modified Scharber plot reflects the fact that a smaller ACC band gap may favor PCE by filling in gaps in the DON acceptor spectrum. We predict that only the candidate molecule with the least good acceptor A, with the highest frontier molecular orbital energies, and one of the larger CV lowest unoccupied molecular orbital (LUMO) - highest unoccupied molecular orbital (HOMO) gaps, will yield a PM6:ACC PCE exceeding that of the parent IDIC-4Cl ACC. This candidate also shows the largest oscillator strength for the primary1 (HOMO, LUMO) charge- transfer transition and the largest degree of delocalization of charge transfer of any of the ACC molecules investigated here.

Promoting Solar-driven Hydrogen Peroxide Production over Thiazole-based Conjugated Polymers via Generating and Converting Singlet Oxygen

Solar-driven synthesis of hydrogen peroxide (H2 O2 ) from water and air provides a low-cost and eco-friendly alternative route to the traditional anthraquinone method. Herein, four thiazole-based conjugated polymers (Tz-CPs: TTz, BTz, TBTz and BBTz) are synthesized via aldimine condensation. BBTz exhibits the highest H2 O2 production rate of 7274 μmol g-1  h-1 in pure water. Further, the reaction path is analyzed by electron paramagnetic resonance (EPR), in situ diffuse reflectance infrared Fourier transform (DRIFT) and theoretical calculation, highlighting the prominent role of singlet oxygen (1 O2 ). The generation of 1 O2 occurs through the oxidation of superoxide radical (⋅O2 - ) and subsequent conversion into endoperoxides via [4+2] cycloaddition over BBTz, which promotes charge separation and reduces the barrier for H2 O2 production. This work provides new insight into the mechanism of photocatalytic O2 reduction and the molecular design of superior single-polymer photocatalysts.

Localization measures of parity adapted U(D)-spin coherent states applied to the phase space analysis of the D-level Lipkin-Meshkov-Glick model

We study phase space properties of critical, parity symmetric, N-qudit systems undergoing a quantum phase transition (QPT) in the thermodynamic N→∞ limit. The D=3 level (qutrit) Lipkin-Meshkov-Glick model is eventually examined as a particular example. For this purpose, we consider U(D)-spin coherent states (DSCS), generalizing the standard D=2 atomic coherent states, to define the coherent state representation Q_{ψ} (Husimi function) of a symmetric N-qudit state |ψ〉 in the phase space CP^{D-1} (complex projective manifold). DSCS are good variational approximations to the ground state of an N-qudit system, especially in the N→∞ limit, where the discrete parity symmetry Z_{2}^{D-1} is spontaneously broken. For finite N, parity can be restored by projecting DSCS onto 2^{D-1} different parity invariant subspaces, which define generalized "Schrödinger cat states" reproducing quite faithfully low-lying Hamiltonian eigenstates obtained by numerical diagonalization. Precursors of the QPT are then visualized for finite N by plotting the Husimi function of these parity projected DSCS in phase space, together with their Husimi moments and Wehrl entropy, in the neighborhood of the critical points. These are good localization measures and markers of the QPT.

Ambipolar Charge Transport in Organic Semiconductors: How Intramolecular Reorganization Energy Is Controlled by Diradical Character

The charged forms of π-conjugated chromophores are relevant in the field of organic electronics as charge carriers in optoelectronic devices, but also as energy storage substrates in organic batteries. In this context, intramolecular reorganization energy plays an important role in controlling material efficiency. In this work, we investigate how the diradical character influences the reorganization energies of holes and electrons by considering a library of diradicaloid chromophores. We determine the reorganization energies with the four-point adiabatic potential method using quantum-chemical calculations at density functional theory (DFT) level. To assess the role of diradical character, we compare the results obtained, assuming both closed-shell and open-shell representations of the neutral species. The study shows how the diradical character impacts the geometrical and electronic structure of neutral species, which in turn control the magnitude of reorganization energies for both charge carriers. Based on computed geometries of neutral and charged species, we propose a simple scheme to rationalize the small, computed reorganization energies for both n-type and p-type charge transport. The study is supplemented with the calculation of intermolecular electronic couplings governing charge transport for selected diradicals, further supporting the ambipolar character of the investigated diradicals.

Jumping Kinetic Monte Carlo: Fast and Accurate Simulations of Partially Delocalized Charge Transport in Organic Semiconductors

Developing devices using disordered organic semiconductors requires accurate and practical models of charge transport. In these materials, charge transport occurs through partially delocalized states in an intermediate regime between localized hopping and delocalized band conduction. Partial delocalization can increase mobilities by orders of magnitude compared to those with conventional hopping, making it important for the design of materials and devices. Although delocalization, disorder, and polaron formation can be described using delocalized kinetic Monte Carlo (dKMC), it is a computationally expensive method. Here, we develop jumping kinetic Monte Carlo (jKMC), a model that approaches the accuracy of dKMC for modest amounts of delocalization (such as those found in disordered organic semiconductors), with a computational cost comparable to that of conventional hopping. jKMC achieves its computational performance by modeling conduction using identical spherical polarons, yielding a simple delocalization correction to the Marcus hopping rate that allows polarons to jump over their nearest neighbors. jKMC can be used in regimes of partial delocalization inaccessible to dKMC to show that modest delocalization can increase mobilities by as much as 2 orders of magnitude.

Charge Localization in Acene Crystals from Ab Initio Electronic Structure

The performance of Koopmans-compliant hybrid functionals in reproducing the electronic structure of organic crystals is tested for a series of acene crystals. The calculated band gaps are found to be consistent with those achieved with the GW method at a fraction of the computational cost and in excellent accord with the experimental results at room temperature, when including the thermal renormalization. The energetics of excess holes and electrons reveals a struggle between polaronic localization and band-like delocalization. The consequences of these results on the transport properties of acene crystals are discussed.

Orbital-resolved observation of singlet fission

Singlet fission1-13 may boost photovoltaic efficiency14-16 by transforming a singlet exciton into two triplet excitons and thereby doubling the number of excited charge carriers. The primary step of singlet fission is the ultrafast creation of the correlated triplet pair17. Whereas several mechanisms have been proposed to explain this step, none has emerged as a consensus. The challenge lies in tracking the transient excitonic states. Here we use time- and angle-resolved photoemission spectroscopy to observe the primary step of singlet fission in crystalline pentacene. Our results indicate a charge-transfer mediated mechanism with a hybridization of Frenkel and charge-transfer states in the lowest bright singlet exciton. We gained intimate knowledge about the localization and the orbital character of the exciton wave functions recorded in momentum maps. This allowed us to directly compare the localization of singlet and bitriplet excitons and decompose energetically overlapping states on the basis of their orbital character. Orbital- and localization-resolved many-body dynamics promise deep insights into the mechanics governing molecular systems18-20 and topological materials21-23.

Highly Ordered Supramolecular Materials of Phase-Separated Block Molecules for Long-Range Exciton Transport

Efficient energy transport over long distances is essential for optoelectronic and light-harvesting devices. Although self-assembled nanofibers of organic molecules are shown to exhibit long exciton diffusion lengths, alignment of these nanofibers into films with large, organized domains with similar properties remains a challenge. Here, it is shown how the functionalization of C₃ -symmetric carbonyl-bridged triarylamine trisamide (CBT) with oligodimethylsiloxane (oDMS) side chains of discrete length leads to fully covered surfaces with aligned domains up to 125 × 70 µm² in which long-range exciton transport takes place. The nanoscale morphology within the domains consists of highly ordered nanofibers with discrete intercolumnar spacings within a soft amorphous oDMS matrix. The oDMS prevents bundling of the CBT fibers, reducing the number of defects within the CBT columns. As a result, the columns have a high degree of coherence, leading to exciton diffusion lengths of a few hundred nanometers with exciton diffusivities (≈0.05 cm² s^-1 ) that are comparable to those of a crystalline tetracene. These findings represent the next step toward fully covered surfaces of highly aligned nanofibers through functionalization with oDMS.

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.

Mechanism of Delocalization-Enhanced Exciton Transport in Disordered Organic Semiconductors

Large exciton diffusion lengths generally improve the performance of organic semiconductor devices, because they enable energy to be transported farther during the exciton lifetime. However, the physics of exciton motion in disordered organic materials is not fully understood, and modeling the transport of quantum-mechanically delocalized excitons in disordered organic semiconductors is a computational challenge. Here, we describe delocalized kinetic Monte Carlo (dKMC), the first model of three-dimensional exciton transport in organic semiconductors that includes delocalization, disorder, and polaron formation. We find that delocalization can dramatically increase exciton transport; for example, delocalization across less than two molecules in each direction can increase the exciton diffusion coefficient by over an order of magnitude. The mechanism for the enhancement is 2-fold: delocalization enables excitons to hop both more frequently and further in each hop. We also quantify the effect of transient delocalization (short-lived periods where excitons become highly delocalized) and show that it depends strongly upon the disorder and transition dipole moments.

High Electron Mobility Hot-Exciton Induced Delayed Fluorescent Organic Semiconductors

The development of high mobility emissive organic semiconductors is of great significance for the fabrication of miniaturized optoelectronic devices, such as organic light emitting transistors. However, great challenge exists in designing key materials, especially those who integrates triplet exciton utilization ability. Herein, dinaphthylanthracene diimides (DNADIs), with 2,6-extended anthracene donor, and 3'- or 4'-substituted naphthalene monoimide acceptors were designed and synthesized. By introducing acceptor-donor-acceptor structure, both materials show high electron mobility. Moreover, by fine-tuning of substitution sites, good integration with high solid state photoluminescence quantum yield of 26 %, high electron mobility of 0.02 cm²  V^-1  s^-1 , and the feature of hot-exciton induced delayed fluorescence were obtained in 4'-DNADI. This work opens a new avenue for developing high electron mobility emissive organic semiconductors with efficient utilization of triplet excitons.

Phonon-Induced Localization of Excitons in Molecular Crystals from First Principles

The spatial extent of excitons in molecular systems underpins their photophysics and utility for optoelectronic applications. Phonons are reported to lead to both exciton localization and delocalization. However, a microscopic understanding of phonon-induced (de)localization is lacking, in particular, how localized states form, the role of specific vibrations, and the relative importance of quantum and thermal nuclear fluctuations. Here, we present a first-principles study of these phenomena in solid pentacene, a prototypical molecular crystal, capturing the formation of bound excitons, exciton-phonon coupling to all orders, and phonon anharmonicity, using density functional theory, the ab initio GW-Bethe-Salpeter equation approach, finite-difference, and path integral techniques. We find that for pentacene zero-point nuclear motion causes uniformly strong localization, with thermal motion providing additional localization only for Wannier-Mott-like excitons. Anharmonic effects drive temperature-dependent localization, and, while such effects prevent the emergence of highly delocalized excitons, we explore the conditions under which these might be realized.

Efficiently predicting directional carrier mobilities in organic materials with the Boltzmann transport equation

Describing charge carrier anisotropy in crystalline organic semiconductors with ab initio methods is challenging because of the weak intermolecular interactions that lead to both localized and delocalized charge transport mechanisms. Small polaron hopping models (localized) are generally used to describe materials with small charge carrier mobilities, while periodic band models (delocalized) are used to describe materials with high charge carrier mobilities. Here, we prove the advantage of applying the constant relaxation time approximation of the Boltzmann transport equation (BTE) to efficiently predict the anisotropic hole mobilities of several unsubstituted (anthracene, tetracene, pentacene, and hexacene) and substituted (2,6-diphenylanthracene, rubrene, and TIPS-pentacene) high-mobility n-acene single crystals. Several density functionals are used to optimize the crystals, and the composite density functional PBEsol0-3c/sol-def2-mSVP predicts the most experimentally similar geometries, adequate indirect bandgaps, and the theoretically consistent n-acene charge transport mobility trend. Similarities between BTE and Marcus mobilities are presented for each crystal. BTE and Marcus charge carrier mobilities computed at the same geometry result in similar mobility trends, differing mostly in materials with more substitutions or structurally complex substituents. By using a reduced number of calculations, BTE is able to predict anisotropic carrier mobilities efficiently and effectively for a range of high-mobility organic semiconductors.

Dynamics of Phonon-Assisted Holes Trapping and Transport over Chemical Defects in Polyethylene

Charge trapping and transport over chemical defects in polyethylene have significant impacts on its electrical and dielectric properties. However, the dynamics of this phenomenon and its underlying mechanisms remain unclear. To understand this fundamental aspect, we conducted a time-domain ab initio nonadiabatic molecular dynamics study of phonon-assisted holes dynamics in polyethylene over C═O and C-OH defect states. Our results suggest that the hole transfer and energy fluctuations substantially depend on temperature and local morphology. When the temperature decreases from 300 to 100 K, the hole transfer efficiency and the energy fluctuations are severely suppressed due to the weakened interactions between holes and phonons. Furthermore, amorphous polyethylene exhibits a severe suppression of the hole transfer process compared to crystalline polyethylene. An explanation for the influence of morphology on the hole transfer process can be found in the differences in the hole-phonon coupling and the electronic coupling between two chemical defect states in crystalline and amorphous polyethylene. Advancing the fundamental understanding of the dynamics of hole transfer over chemical effects in polymers is a key to improving their insulating properties for the next-generation high-voltage cables.

Nonpolymer Organic Solar Cells: Microscopic Phonon Control to Suppress Nonradiative Voltage Loss via Charge-Separated State

Recent remarkable developments on nonfullerene solar cells have reached a photoelectric conversion efficiency (PCE) of 18% by tuning the band energy levels in small molecular acceptors. In this regard, understanding the impact of small donor molecules on nonpolymer solar cells is essential. Here, we systematically investigated mechanisms of solar cell performance using diketopyrrolopyrrole (DPP)-tetrabenzoporphyrin (BP) conjugates of C4-DPP-H₂BP and C4-DPP-ZnBP, where C4 represents the butyl group substituted at the DPP unit as small p-type molecules, while an acceptor of [6,6]-phenyl-C₆₁-buthylic acid methyl ester is employed. We clarified the microscopic origins of the photocarrier caused by phonon-assisted one-dimensional (1D) electron-hole dissociations at the donor-acceptor interface. Using a time-resolved electron paramagnetic resonance, we have characterized controlled charge-recombination by manipulating disorders in π-π donor stacking. This ensures carrier transport through stacking molecular conformations to suppress nonradiative voltage loss capturing specific interfacial radical pairs separated by 1.8 nm in bulk-heterojunction solar cells. We show that, while disordered lattice motions by the π-π stackings via zinc ligation are essential to enhance the entropy for charge dissociations at the interface, too much ordered crystallinity causes the backscattering phonon to reduce the open-circuit voltage by geminate charge-recombination.

Ultrafast and Sensitive Self-Powered Photodetector Based on Graphene/Pentacene Single Crystal Heterostructure with Weak Light Detection Capacity

Organic materials exhibit efficient light absorption and low-temperature, large-scale processability, and have stimulated enormous research efforts for next-generation optoelectronics. While, high-performance organic devices with fast speed and high responsivity still face intractable challenges, due to their intrinsic limitations including finite carrier mobility and high exciton binding energy. Here an ultrafast and highly sensitive broadband phototransistor is demonstrated by integrating high-quality pentacene single crystal with monolayer graphene. Encouragingly, the -3 dB bandwidth can reach up to 26 kHz, which is a record-speed for such sensitized organic phototransistors. Enormous absorption, long exciton diffusion length of pentacene crystal, and efficient interfacial charge transfer enable a high responsivity of >10⁵  A W^-1  and specific detectivity of >10¹¹  Jones. Moreover, self-powered weak-light detection is realized using a simple asymmetric configuration, and the obvious zero-bias photoresponses can be displayed even under 750 nW cm^-2  light intensity. Excellent response speed and photoresponsivity enable high-speed image sensor capability in UV-Vis ranges.  The results offer a practical strategy for constructing high-performance self-powered organic hybrid photodetectors, with strong applicability in wireless, weak-light detection, and video-frame-rate imaging applications.

Direct observation of ultrafast singlet exciton fission in three dimensions

We present quantitative ultrafast interferometric pump-probe microscopy capable of tracking of photoexcitations with sub-10 nm spatial precision in three dimensions with 15 fs temporal resolution, through retrieval of the full transient photoinduced complex refractive index. We use this methodology to study the spatiotemporal dynamics of the quantum coherent photophysical process of ultrafast singlet exciton fission. Measurements on microcrystalline pentacene films grown on glass (SiO₂) and boron nitride (hBN) reveal a 25 nm, 70 fs expansion of the joint-density-of-states along the crystal a,c-axes accompanied by a 6 nm, 115 fs change in the exciton density along the crystal b-axis. We propose that photogenerated singlet excitons expand along the direction of maximal orbital π-overlap in the crystal a,c-plane to form correlated triplet pairs, which subsequently electronically decouples into free triplets along the crystal b-axis due to molecular sliding motion of neighbouring pentacene molecules. Our methodology lays the foundation for the study of three dimensional transport on ultrafast timescales.

Computational Method for Evaluating the Thermoelectric Power Factor for Organic Materials Modeled by the Holstein Model: A Time-Dependent Density Matrix Renormalization Group Formalism

Organic/polymeric materials are of emerging importance for thermoelectric conversion. The soft nature of these materials implies strong electron-phonon coupling, often leading to carrier localization. This poses great challenges for the conventional Boltzmann transport description based on relaxation time approximation and band structure calculations. In this work, combining the Kubo formula with the finite-temperature time-dependent density matrix renormalization group (FT-TD-DMRG) in the grand canonical ensemble, we developed a nearly exact algorithm to calculate the thermoelectric power factor PF = α² σ, where α is the Seebeck coefficient and σ is the electrical conductivity, and apply the algorithm to Holstein Hamiltonian with electron-phonon coupling to model organic materials. Our algorithm can provide a unified description covering the weak coupling limit described by the bandlike Boltzmann transport to the strong coupling hopping limit.

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.

Quantitative Hole Mobility Simulation and Validation in Substituted Acenes

Knowledge of the full phonon spectrum is essential to accurately calculate the dynamic disorder (σ) and hole mobility (μ_h) in organic semiconductors (OSCs). However, most vibrational spectroscopy techniques under-measure the phonons, thus limiting the phonon validation. Here, we measure and model the full phonon spectrum using multiple spectroscopic techniques and predict μ_h using σ from only the Γ-point and the full Brillouin zone (FBZ). We find that only inelastic neutron scattering (INS) provides validation of all phonon modes, and that σ in a set of small molecule semiconductors can be miscalculated by up to 28% when comparing Γ-point against FBZ calculations. A subsequent mode analysis shows that many modes contribute to σ and that no single mode dominates. Our results demonstrate the importance of a thoroughly validated phonon calculation, and a need to develop design rules considering the full spectrum of phonon modes.

Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package

Constrained density functional theory (CDFT) is a powerful tool for the prediction of electron transfer parameters in condensed phase simulations at a reasonable computational cost. In this work we present an extension to CDFT in the popular mixed Gaussian/plane wave electronic structure package CP2K, implementing the additional force terms arising from a constraint based on Hirshfeld charge partitioning. This improves upon the existing Becke partitioning scheme, which is prone to give unphysical atomic charges. We verify this implementation for a variety of systems: electron transfer in (H₂O)₂⁺ in a vacuum, electron tunnelling between oxygen vacancy centers in solid MgO, and electron self-exchange in aqueous Ru²⁺–Ru³⁺. We find good agreement with previous plane-wave CDFT results for the same systems, but at a significantly lower computational cost, and we discuss the general reliability of condensed phase CDFT calculations.

Exciton transport in molecular organic semiconductors boosted by transient quantum delocalization

Designing molecular materials with very large exciton diffusion lengths would remove some of the intrinsic limitations of present-day organic optoelectronic devices. Yet, the nature of excitons in these materials is still not sufficiently well understood. Here we present Frenkel exciton surface hopping, an efficient method to propagate excitons through truly nano-scale materials by solving the time-dependent Schrödinger equation coupled to nuclear motion. We find a clear correlation between diffusion constant and quantum delocalization of the exciton. In materials featuring some of the highest diffusion lengths to date, e.g. the non-fullerene acceptor Y6, the exciton propagates via a transient delocalization mechanism, reminiscent to what was recently proposed for charge transport. Yet, the extent of delocalization is rather modest, even in Y6, and found to be limited by the relatively large exciton reorganization energy. On this basis we chart out a path for rationally improving exciton transport in organic optoelectronic materials.

Charge Transport in Organic Semiconductors: The Perspective from Nonadiabatic Molecular Dynamics

Organic semiconductors (OSs) are an exciting class of materials that have enabled disruptive technologies in this century including large-area electronics, flexible displays, and inexpensive solar cells. All of these technologies rely on the motion of electrical charges within the material and the diffusivity of these charges critically determines their performance. In this respect, it is remarkable that the nature of the charge transport in these materials has puzzled the community for so many years, even for apparently simple systems such as molecular single crystals: some experiments would better fit an interpretation in terms of a localized particle picture, akin to molecular or biological electron transfer, while others are in better agreement with a wave-like interpretation, more akin to band transport in metals.

Exciting recent progress in the theory and simulation of charge carrier transport in OSs has now led to a unified understanding of these disparate findings, and this Account will review one of these tools developed in our laboratory in some detail: direct charge carrier propagation by quantum-classical nonadiabatic molecular dynamics. One finds that even in defect-free crystals the charge carrier can either localize on a single molecule or substantially delocalize over a large number of molecules depending on the relative strength of electronic couplings between the molecules, reorganization, or charge trapping energy of the molecule and thermal fluctuations of electronic couplings and site energies, also known as electron–phonon couplings.

Our simulations predict that in molecular OSs exhibiting some of the highest measured charge mobilities to date, the charge carrier forms “flickering” polarons, objects that are delocalized over 10–20 molecules on average and that constantly change their shape and extension under the influence of thermal disorder. The flickering polarons propagate through the OS by short (≈10 fs long) bursts of the wave function that lead to an expansion of the polaron to about twice its size, resulting in spatial displacement, carrier diffusion, charge mobility, and electrical conductivity. Arguably best termed “transient delocalization”, this mechanistic scenario is very similar to the one assumed in transient localization theory and supports its assertions. We also review recent applications of our methodology to charge transport in disordered and nanocrystalline samples, which allows us to understand the influence of defects and grain boundaries on the charge propagation.

Unfortunately, the energetically favorable packing structures of typical OSs, whether molecular or polymeric, places fundamental constraints on charge mobilities/electronic conductivity compared to inorganic semiconductors, which limits their range of applications. In this Account, we review the design rules that could pave the way for new very high-mobility OS materials and we argue that 2D covalent organic frameworks are one of the most promising candidates to satisfy them.

We conclude that our nonadiabatic dynamics method is a powerful approach for predicting charge carrier transport in crystalline and disordered materials. We close with a brief outlook on extensions of the method to exciton transport, dissociation, and recombination. This will bring us a step closer to an understanding of the birth, survival, and annihiliation of charges at interfaces of optoelectronic devices.

Thermally activated intra-chain charge transport in high charge-carrier mobility copolymers

Disordered or even seemingly amorphous, donor–acceptor type, conjugated copolymers with high charge-carrier mobility have emerged as a new class of functional materials, where transport along the conjugated backbone is key. Here, we report on non-adiabatic molecular dynamics simulations of charge-carrier transport along chains of poly (indacenodithiophene-co-benzothiadiazole), within a model Hamiltonian parameterized against first-principles calculations. We predict thermally activated charge transport associated with a slightly twisted ground-state conformation, on par with experimental results. Our results also demonstrate that the energy mismatch between the hole on the donor vs the acceptor units of the copolymer drives localization of the charge carriers and limits the intra-chain charge-carrier mobility. We predict that room-temperature mobility values in excess of 10 cm2 V−1 s−1 can be achieved through proper chemical tuning of the component monomer units.

Efficient Surface Hopping Approach for Modeling Charge Transport in Organic Semiconductors

The trajectory surface hopping (TSH) method is nowadays widely applied to study the charge/exciton transport process in organic semiconductors (OSCs). In the present study, we systematically examine the performance of two approximations in the fewest switched surface hopping (FSSH) simulations for charge transport (CT) in several representative OSCs. These approximations include (i) the substitution of the nuclear velocity scaling along the nonadiabatic coupling vector (NCV) by rescaling the hopping probability with the Boltzmann factor (Boltzmann correction (BC)) and (ii) a phenomenological approach to treat the quantum feedback from the electronic system to the nuclear system (implicit charge relaxation (IR)) in the OSCs. We find that charge mobilities computed by FSSH-BC-IR are in very good agreement with the mobilities obtained by standard FSSH simulations with explicit charge relaxation (FSSH-ER), however, at reduced computational cost. A key parameter determining the charge carrier mobility is the reorganization energy, which is sensitively dependent on DFT functionals applied. By employing the IR approximation, the FSSH method allows systematic investigation of the effect of the reorganization energies obtained by different DFT functionals like B3LYP or ωB97XD on CT in OSCs. In comparison to the experiments, FSSH-BC-IR using ωB97XD reorganization energy underestimates mobilities in the low-coupling regime, which may indicate the lack of nuclear quantum effects (e.g., zero point energy (ZPE)) in the simulations. The mobilities obtained by FSSH-BC-IR using the B3LYP reorganization energy agree well with experimental values in 3 orders of magnitude. The accidental agreement may be the consequence of the underestimation of the reorganization energy by the B3LYP functional, which compensates for the neglect of nuclear ZPE in the simulations.

All-Atom Nonadiabatic Dynamics Simulation of Hybrid Graphene Nanoribbons Based on Wannier Analysis and Machine Learning

Trajectory surface hopping combined with ab initio electronic structure calculations is a popular and powerful approach for on-the-fly nonadiabatic dynamics simulations. For large systems, however, this remains a significant challenge because of the unaffordable computational cost of large-scale electronic structure calculations. Here, we present an efficient divide-and-conquer approach to construct the system Hamiltonian based on Wannier analysis and machine learning. In detail, the large system under investigation is first decomposed into small building blocks, and then all possible segments formed by building blocks within a cutoff distance are found out. Ab initio molecular dynamics is carried out to generate a sequence of geometries for each equivalent segment with periodicity. The Hamiltonian matrices in the maximum localized Wannier function (MLWF) basis are obtained for all geometries and utilized to train artificial neural networks (ANNs) for the structure-dependent Hamiltonian elements. Taking advantage of the orthogonality and spatial locality of MLWFs, the one-electron Hamiltonian of a large system at arbitrary geometry can be directly constructed by the trained ANNs. As demonstrations, we study charge transport in a zigzag graphene nanoribbon (GNR), a coved GNR, and a series of hybrid GNRs with a state-of-the-art surface hopping method. The interplay between delocalized and localized states is found to determine the electron dynamics in hybrid GNRs. Our approach has successfully studied GNRs with >10 000 atoms, paving the way for efficient and reliable all-atom nonadiabatic dynamics simulation of general systems.

Temperature-Dependent Structural Phase Transition in Rubrene Single Crystals: The Missing Piece from the Charge Mobility Puzzle?

Accurate structural models for rubrene, the benchmark organic semiconductor, derived from synchrotron X-ray data in the temperature range of 100-300 K, show that its cofacially stacked tetracene backbone units remain blocked with respect to each other upon cooling to 200 K and start to slip below that temperature. The release of the blocked slippage occurs at approximately the same temperature as the hole mobility crossover. The blocking between 200 and 300 K is caused by a negative correlation between the relatively small thermal expansion along the crystallographic b-axis and the relatively large widening of the angle between herringbone-stacked tetracene units. DFT calculations reveal that this blocked slippage is accompanied by a discontinuity in the variation with temperature of the electronic couplings associated with hole transport between cofacially stacked tetracene backbones.

Ultrafast estimation of electronic couplings for electron transfer between pi-conjugated organic molecules. II

The development of highly efficient methods for the calculation of electronic coupling matrix elements between the electron donor and acceptor is an important goal in theoretical organic semiconductor research. In Paper I [F. Gajdos, S. Valner, F. Hoffmann, J. Spencer, M. Breuer, A. Kubas, M. Dupuis, and J. Blumberger, J. Chem. Theory Comput. 10, 4653 (2014)], we introduced the analytic overlap method (AOM) for this purpose, which is an ultrafast electronic coupling estimator parameterized to and orders of magnitude faster than density functional theory (DFT) calculations at a reasonably small loss in accuracy. In this work, we reparameterize and extend the AOM to molecules containing nitrogen, oxygen, fluorine, and sulfur heteroatoms using 921 dimer configurations from the recently introduced HAB79 dataset. We find again a very good linear correlation between the frontier orbital overlap, calculated ultrafast in an optimized minimum Slater basis, and DFT reference electronic couplings. The new parameterization scheme is shown to be transferable to sulfur-containing polyaromatic hydrocarbons in experimentally resolved dimeric configurations. Our extension of the AOM enables high-throughput screening of very large databases of chemically diverse organic crystal structures and the application of computationally intense non-adiabatic molecular dynamics methods to charge transport in state-of-the-art organic semiconductors, e.g., non-fullerene acceptors.

HAB79: A new molecular dataset for benchmarking DFT and DFTB electronic couplings against high-level ab initio calculations

A new molecular dataset called HAB79 is introduced to provide ab initio reference values for electronic couplings (transfer integrals) and to benchmark density functional theory (DFT) and density functional tight-binding (DFTB) calculations. The HAB79 dataset is composed of 79 planar heterocyclic polyaromatic hydrocarbon molecules frequently encountered in organic (opto)electronics, arranged to 921 structurally diverse dimer configurations. We show that CASSCF/NEVPT2 with a minimal active space provides a robust reference method that can be applied to the relatively large molecules of the dataset. Electronic couplings are largest for cofacial dimers, in particular, sulfur-containing polyaromatic hydrocarbons, with values in excess of 0.5 eV, followed by parallel displaced cofacial dimers. V-shaped dimer motifs, often encountered in the herringbone layers of organic crystals, exhibit medium-sized couplings, whereas T-shaped dimers have the lowest couplings. DFT values obtained from the projector operator-based diabatization (POD) method are initially benchmarked against the smaller databases HAB11 (HAB7-) and found to systematically improve when climbing Jacob's ladder, giving mean relative unsigned errors (MRUEs) of 27.7% (26.3%) for the generalized gradient approximation (GGA) functional BLYP, 20.7% (15.8%) for hybrid functional B3LYP, and 5.2% (7.5%) for the long-range corrected hybrid functional omega-B97X. Cost-effective POD in combination with a GGA functional and very efficient DFTB calculations on the dimers of the HAB79 database give a good linear correlation with the CASSCF/NEVPT2 reference data, which, after scaling with a multiplicative constant, gives reasonably small MRUEs of 17.9% and 40.1%, respectively, bearing in mind that couplings in HAB79 vary over 4 orders of magnitude. The ab initio reference data reported here are expected to be useful for benchmarking other DFT or semi-empirical approaches for electronic coupling calculations.

Impact of Nanoscale Morphology on Charge Carrier Delocalization and Mobility in an Organic Semiconductor

A central challenge of organic semiconductor research is to make cheap, disordered materials that exhibit high electrical conductivity. Unfortunately, this endeavor is hampered by the poor fundamental understanding of the relationship between molecular packing structure and charge carrier mobility. Here a novel computational methodology is presented that fills this gap. Using a melt-quench procedure it is shown that amorphous pentacene spontaneously self-assembles to nanocrystalline structures that, at long quench times, form the characteristic herringbone layer of the single crystal. Quantum dynamical simulations of electron hole transport show a clear correlation between the crystallinity of the sample, the quantum delocalization, and the mobility of the charge carrier. Surprisingly, the long-held belief that charge carriers form relatively localized polarons in disordered OS is only valid for fully amorphous structures-for nanocrystalline and crystalline samples, significant charge carrier delocalization over several nanometers occurs that underpins their improved conductivities. The good agreement with experimentally available data makes the presented methodology a robust computational tool for the predictive engineering of disordered organic materials.

High-throughput virtual screening for organic electronics: a comparative study of alternative strategies

We present a review of the field of high-throughput virtual screening for organic electronics materials focusing on the sequence of methodological choices that determine each virtual screening protocol. These choices are present in all high-throughput virtual screenings and addressing them systematically will lead to optimised workflows and improve their applicability. We consider the range of properties that can be computed and illustrate how their accuracy can be determined depending on the quality and size of the experimental datasets. The approaches to generate candidates for virtual screening are also extremely varied and their relative strengths and weaknesses are discussed. The analysis of high-throughput virtual screening is almost never limited to the identification of top candidates and often new patterns and structure-property relations are the most interesting findings of such searches. The review reveals a very dynamic field constantly adapting to match an evolving landscape of applications, methodologies and datasets.

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.