Mesoscale molecular network formation in amorphous organic materials

Self-assembly prediction of architecture-controlled bottlebrush copolymers in solution using graph convolutional networks

The investigation of bottlebrush copolymer self-assembly in solution involves a comprehensive approach integrating simulation and experimental research, due to their unique physical characteristics. However, the intricate architecture of bottlebrush copolymers and the diverse solvent conditions introduce a wide range of parameter spaces. In this study, we investigated the solution self-assembly behavior of bottlebrush copolymers by combining dissipative particle dynamics (DPD) simulation results and machine learning (ML) including graph convolutional networks (GCNs). The architecture of bottlebrush copolymers is encoded by graphs including connectivity, side chain length, bead types, and interaction parameters of DPD simulation. Using GCN, we accurately predicted the single chain properties of bottlebrush copolymers with over 95% accuracy. Furthermore, phase behavior was precisely predicted using these single chain properties. Shapley additive explanations (SHAP) values of single chain properties to the various self-assembly morphologies were calculated to investigate the correlation between single chain properties and morphologies. In addition, we analyzed single chain properties and phase behavior as a function of DPD interaction parameters, extracting relevant physical properties for vesicle morphology formation. This work paves the way for tailored design in solution of self-assembled nanostructures of bottlebrush copolymers, offering a GCN framework for precise prediction of self-assembly morphologies under various chain architectures and solvent conditions.

The effects of side chain engineering on the morphology and charge transport of the A-DA1D-A type of non-fullerene acceptor: a multiscale study

The appearance of the A-DA1D-A type of non-fullerene acceptor Y6 and its derivatives significantly improves the power conversion efficiency of organic solar cells. However, the effects of the modulation of the side chains of Y6 on its morphology and charge transport in organic thin films are still not well understood. In this work, we have systematically studied the effects of symmetric modifications of the length of alkyl side chains and the types, such as branched or straight alkyl chains, and the introduction of heteroatoms to side chains on these properties. A multiscale study, including density functional theory and classical molecular dynamics simulations, has been used to answer this open question. We find that face-on configurations are generally dominant for the AA, A1A1, and DD stacking of molecular pairs. With respect to prototype Y6, the introduction of oxygen atoms to outer alkyl side chains could enhance AA stacking but worsen the electrical network and enlarge the reorganization energy during electron transfer, and changing outer side straight alkyl chains to branched chains ruins π-π stacking of all units significantly. Finally, we discover that shortening outer alkyl side chains appropriately or changing inner branched chains to straight chains with the same number of carbon atoms is a good strategy to improve the molecular π-π stacking and electron mobility of Y6 while changing outer straight side chains to branched chains or introducing oxygen atoms to outer straight chains is the opposite. This study provides a new insight into the relationship between morphology and electron mobility and will be helpful for the design of future high-performance non-fullerene acceptors.

Controlling swelling in mixed transport polymers through alkyl side-chain physical cross-linking

Semiconducting conjugated polymers bearing glycol side chains can simultaneously transport both electronic and ionic charges with high charge mobilities, making them ideal electrode materials for a range of bioelectronic devices. However, heavily glycolated conjugated polymer films have been observed to swell irreversibly when subjected to an electrochemical bias in an aqueous electrolyte. The excessive swelling can lead to the degradation of their microstructure, and subsequently reduced device performance. An effective strategy to control polymer film swelling is to copolymerize glycolated repeat units with a fraction of monomers bearing alkyl side chains, although the microscopic mechanism that constrains swelling is unknown. Here we investigate, experimentally and computationally, a series of archetypal mixed transporting copolymers with varying ratios of glycolated and alkylated repeat units. Experimentally we observe that exchanging 10% of the glycol side chains for alkyl leads to significantly reduced film swelling and an increase in electrochemical stability. Through molecular dynamics simulation of the amorphous phase of the materials, we observe the formation of polymer networks mediated by alkyl side-chain interactions. When in the presence of water, the network becomes increasingly connected, counteracting the volumetric expansion of the polymer film.

Electron Transport in Organic Photovoltaic Acceptor Materials: Improving the Carrier Mobilities by Intramolecular and Intermolecular Modulations

High carrier mobility is beneficial to increase the active-layer thickness while maintaining a high fill factor, which is crucial to further improve the light harvesting and organic photovoltaic efficiency. The aim of this Perspective is to elucidate the electron transport mechanisms in prototypical non-fullerene (NF) acceptors through our recent theoretical studies. The electron transport in A-D-A small-molecule acceptors (SMAs), e.g., ITIC and Y6, is mainly determined by end-group π-π stacking. Relative to ITIC, the angular backbone along with more flexible side chains leads to Y6 having a closer stacking and enhanced intermolecular electronic connectivity. For polymerized rylene diimide acceptors, to achieve high electron mobilities, they need to simultaneously increase intramolecular and intermolecular connectivity. Finally, finely tuning the π-bridge modes to enhance intramolecular superexchange coupling is essential to develop novel polymerized A-D-A SMAs.

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.

Design of a Fully Non-Fused Bulk Heterojunction toward Efficient and Low-Cost Organic Photovoltaics

To modulate the miscibility between donor and acceptor materials both possessing fully non-fused ring structures, a series of electron acceptors (A4T-16, A4T-31 and A4T-32) with different polar functional substituents were synthesized and investigated. The three acceptors show good planarity, high conformational stability, complementary absorption and energy levels with the non-fused polymer donor (PTVT-BT). Among them, A4T-32 possesses the strongest polar functional group and shows the highest surface energy, which facilitates morphological modulation in the bulk heterojunction (BHJ) blend. Benefiting from the proper morphology control method, an impressive power conversion efficiency (PCE) of approaching 16.0 % and a superior fill factor over 0.795 are achieved in the PTVT-BT : A4T-32-based organic photovoltaic cells with superior photoactive materials price advantage, which represent the highest value for the cells based on the non-fused blend films. Notably, this cell maintains ≈84 % of its initial PCE after nearly 2000 h under the continuous simulated 1-sun-illumination. In addition, the flexible PTVT-BT : A4T-32-based cells were fabricated and delivered a decent PCE of 14.6 %. This work provides an effective molecular design strategy for the non-fused non-fullerene acceptors (NFAs) from the aspect of bulk morphology control in fully non-fused BHJ layers, which is crucial for their practical applications.

Controlling Aggregation-Induced Two-Photon Absorption Enhancement via Intermolecular Interactions

Historically, two-photon absorption (2PA) cross sections reported in the literature have been derived from solution-phase measurements. However, such techniques fail to grasp the implications of how these cross sections can be impacted by varying degrees of aggregation or in the condensed phase as bulk solids or thin films. For a precise determination of how aggregation impacts 2PA at a molecular level, computational methods present themselves as ideal. Herein, a series of quadrupolar π-conjugated dyes were simulated by molecular dynamics (MD) in the gas phase and condensed phase. In the condensed phase, their intermolecular interactions and electronic coupling behavior were fully characterized, both quantitatively and qualitatively. Using quadratic-response time-dependent density functional theory, 2PA cross sections of structures derived from MD trajectories were calculated. Comparisons are made between gas-phase and condensed-phase results, and enhancement factors are defined to show how certain dyes may experience changes in their respective 2PA cross sections as a function of aggregation. It was found that these cross sections depend heavily on conformational locking in the condensed phase and relative stacking arrangements. J-aggregates were associated with enhanced 2PA and H-aggregates with quenched 2PA activity. However, in a highly disordered aggregate, the effects of these stacking arrangements are averaged out of the bulk result, and the effects of conformational locking dominate.

Systematic Merging of Nonfullerene Acceptor π-Extension and Tetrafluorination Strategies Affords Polymer Solar Cells with >16% Efficiency

The end-capping group (EG) is the essential electron-withdrawing component of nonfullerene acceptors (NFAs) in bulk heterojunction (BHJ) organic solar cells (OSCs). To systematically probe the impact of two frequent EG functionalization strategies, π-extension and halogenation, in A-DAD-A type NFAs, we synthesized and characterized four such NFAs: BT-BIC, LIC, L4F, and BO-L4F. To assess the relative importance of these strategies, we contrast these NFAs with the baseline acceptors, Y5 and Y6. Up to 16.6% power conversion efficiency (PCE) in binary inverted OSCs with BT-BO-L4F combining π-extension and halogenation was achieved. When these two factors are combined, the effect on optical absorption is cumulative. Single-crystal π-π stacking distances are similar for the EG strategies of π-extension. Increasing the alkyl substituent length from BT-L4F to BT-BO-L4F significantly alters the packing motif and eliminates the EG core interactions of BT-L4F. Electronic structure computations reveal some of the largest NFA π-π electronic couplings observed to date, 103.8 meV in BT-L4F and 47.5 meV in BT-BO-L4F. Computed electronic reorganization energies, 132 and 133 meV for BT-L4F and BT-BO-L4F, respectively, are also lower than Y6 (150 meV). BHJ blends show preferential π-face-on orientation, and both fluorination and π-extension increase NFA crystallinity. Femto/nanosecond transient absorption spectroscopy (fs/nsTA) and integrated photocurrent device analysis (IPDA) indicate that π-extension modifies the phase separation to enhance film ordering and carrier mobility, while fluorination suppresses unimolecular recombination. This systematic study highlights the synergistic effects of NFA π-extension and fluorination in affording efficient OSCs and provides insights into designing next-generation materials.

Dynamics of the sub-ambient gelation and shearing of solutions of P3HT and P3HT blends towards active layer formation in bulk heterojunction organic solar cells

Organic solar cells (OSCs) containing an active layer consisting of a nanostructured blend of a conjugated polymer like poly(3-hexylthiophene) (P3HT) and an electron acceptor have the potential of competing against silicon-based photovoltaic panels. However, this potential is largely unfulfilled first due to interrelated production and stability issues of organic solar cells and second due to the unscalable nature of the generally employed spin coating process used for the fabrication of organic solar cells. Furthermore, alternatives to spin coating, especially relying on continuous polymer processing methods like extrusion and coating, cannot be readily applied due to the typically low shear viscosity and elasticity of polymer solutions making up the active layer. Recently, He et al. have reported that the gelation of P3HT with [6,6]-phenyl-C61-butyric acid methyl ester (PC60BM) under sub-ambient conditions can provide a new route to the processing of the active layers of bulk heterojunction solar cells. Furthermore, increases in power conversion efficiencies (PCEs) of the P3HT/PC60BM active layer were determined to be possible under certain shearing and thermal histories of the P3HT/PC60BM gels. Here oscillatory and steady torsional flows were used to investigate the gel formation dynamics of P3HT with a recently proposed non-fullerene acceptor o-IDTBR under sub-ambient conditions and compared with the gelation behavior of P3HT/PC60BM blends. The rheological material functions as well as the gel strengths defined on the basis of linear viscoelastic material functions, characterized via small-amplitude oscillatory shearing, were observed to be functions of the P3HT and o-IDTBR concentrations, the solvent used and the shearing conditions. Overall, the P3HT gels which formed upon quenching to sub-zero temperatures were found to be stable during small-amplitude oscillatory shear (linear viscoelastic range) but broke down even at the relatively low shear rates associated with steady torsional flows, suggesting that the shearing conditions used during the processing of gels of P3HT and blends of P3HT with small molecule acceptors can alter the gel structure, possibly leading to changes in the resulting active layer performance.

Nanotechnology for catalysis and solar energy conversion

This roadmap on Nanotechnology for Catalysis and Solar Energy Conversion focuses on the application of nanotechnology in addressing the current challenges of energy conversion: 'high efficiency, stability, safety, and the potential for low-cost/scalable manufacturing' to quote from the contributed article by Nathan Lewis. This roadmap focuses on solar-to-fuel conversion, solar water splitting, solar photovoltaics and bio-catalysis. It includes dye-sensitized solar cells (DSSCs), perovskite solar cells, and organic photovoltaics. Smart engineering of colloidal quantum materials and nanostructured electrodes will improve solar-to-fuel conversion efficiency, as described in the articles by Waiskopf and Banin and Meyer. Semiconductor nanoparticles will also improve solar energy conversion efficiency, as discussed by Boschloo et al in their article on DSSCs. Perovskite solar cells have advanced rapidly in recent years, including new ideas on 2D and 3D hybrid halide perovskites, as described by Spanopoulos et al 'Next generation' solar cells using multiple exciton generation (MEG) from hot carriers, described in the article by Nozik and Beard, could lead to remarkable improvement in photovoltaic efficiency by using quantization effects in semiconductor nanostructures (quantum dots, wires or wells). These challenges will not be met without simultaneous improvement in nanoscale characterization methods. Terahertz spectroscopy, discussed in the article by Milot et al is one example of a method that is overcoming the difficulties associated with nanoscale materials characterization by avoiding electrical contacts to nanoparticles, allowing characterization during device operation, and enabling characterization of a single nanoparticle. Besides experimental advances, computational science is also meeting the challenges of nanomaterials synthesis. The article by Kohlstedt and Schatz discusses the computational frameworks being used to predict structure-property relationships in materials and devices, including machine learning methods, with an emphasis on organic photovoltaics. The contribution by Megarity and Armstrong presents the 'electrochemical leaf' for improvements in electrochemistry and beyond. In addition, biohybrid approaches can take advantage of efficient and specific enzyme catalysts. These articles present the nanoscience and technology at the forefront of renewable energy development that will have significant benefits to society.

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

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

Photophysical implications of ring fusion, linker length, and twisting angle in a series of perylenediimide-thienoacene dimers

Perylenediimide (PDI) derivatives have been widely studied as electron acceptor alternatives to fullerenes in organic photovoltaics (OPVs) because of their tunable absorption in the visible range, inexpensive synthesis, and photochemical stability. A common motif for improving device efficiency involves joining multiple PDIs together through electron-rich linkers to form a twisted acceptor-donor-acceptor molecule. Molecular features such as ring fusion are further employed to modify the structure locally and in films. These synthetic efforts have greatly enhanced OPV device efficiencies, however it remains unclear how the increasingly elaborate structural modifications affect the photophysical processes integral to efficient photon-to-charge conversion. Here we carry out a systematic study of a series of PDI dimers with thienoacene linkers in which the twist angle, linker length, and degree of ring fusion are varied to investigate the effects of these structural features on the molecular excited states and exciton recombination dynamics. Spectroscopic characterization of the dimers suggest that ring fusion causes greater coupling between the donor and acceptor components and greatly enhances the lifetime of a thienoacene to PDI charge transfer state. The lifetime of this CT state also correlates well with the linker-PDI dihedral angle, with smaller dihedral angle resulting in longer lifetime. DFT and two-photon absorption TDDFT calculations were developed in-house to model the ground state and excited transitions, providing theoretical insight into the reasons for the observed photophysical properties and identifying the charge transfer state in the excited state absorption spectra. These results highlight how the longevity of the excited state species, important for the efficient conversion of excitons to free carriers in OPV devices, can be chemically tuned by controlling ring fusion and by using steric effects to control the relative orientations of the molecular fragments. The results provide a successful rationalization of the behavior of solar cells involving these acceptor molecules.

Three dimensional nanoscale analysis reveals aperiodic mesopores in a covalent organic framework and conjugated microporous polymer

The integrated analytical approach developed in this study offers a powerful methodology for the structural characterisation of complex molecular nanomaterials. Structures of a covalent organic framework based on boronate esters (COF-5) and a conjugated microporous polymer (Aza-CMP) have been investigated by a combination of several electron microscopy techniques elucidating the three-dimensional topology of the complex polycrystalline (COF) and non-crystalline (CMP) materials. Unexpected, aperiodic mesoporous channels of 20-50 nm in diameter were found to be penetrating the COF and CMP particles, which cannot be detected by X-ray diffraction techniques. The mesopores appear to be stable under a range of different conditions and accessible to gas molecules, exhibiting a particular bonding capability with CO2 in the case of the CMP. The mesoporosity is unrelated to the intrinsic chemical structures of the COF or CMP but rather it reflects the mechanisms of polymer particle formation in a polycondensation reaction. The mesopores may be templated by clusters of solvent molecules during the COF or CMP synthesis, leaving cavities within the polymer particles. The unexpected mesoporosity discovered in COF and CMP materials begs for re-assessment of the nature of framework materials and may open new opportunities for applications of these molecular materials in gas sorption or catalysis.

Eigenvector centrality for characterization of protein allosteric pathways

Determining the principal energy-transfer pathways responsible for allosteric communication in biomolecules remains challenging, partially due to the intrinsic complexity of the systems and the lack of effective characterization methods. In this work, we introduce the eigenvector centrality metric based on mutual information to elucidate allosteric mechanisms that regulate enzymatic activity. Moreover, we propose a strategy to characterize the range of correlations that underlie the allosteric processes. We use the V-type allosteric enzyme imidazole glycerol phosphate synthase (IGPS) to test the proposed methodology. The eigenvector centrality method identifies key amino acid residues of IGPS with high susceptibility to effector binding. The findings are validated by solution NMR measurements yielding important biological insights, including direct experimental evidence for interdomain motion, the central role played by helix h[Formula: see text], and the short-range nature of correlations responsible for the allosteric mechanism. Beyond insights on IGPS allosteric pathways and the nature of residues that could be targeted by therapeutic drugs or site-directed mutagenesis, the reported findings demonstrate the eigenvector centrality analysis as a general cost-effective methodology to gain fundamental understanding of allosteric mechanisms at the molecular level.

Structural Correlations and Percolation in Twisted Perylene Diimides Using a Simple Anisotropic Coarse-Grained Model

Large, twisted, and fused conjugated molecular architectures have begun to appear more prominently in the organic semiconductor literature. From a modeling perspective, such structures present a challenge to conventional simulation techniques; atomistic resolutions are computationally inefficient, while traditional isotropic coarse-grained models do not capture the inherent anisotropies of the molecules. In this work, we develop a simple coarse-grained model that explicitly incorporates the anisotropy of these molecular architectures, thereby providing a route toward analyzing π-stacking, and thus qualitative electronic structure, at a computationally efficient coarse-grained resolution. Our simple coarse-grained model maintains relative orientations of conjugated rings, as well as inter-ring dihedrals, that are critical for understanding electronic and excitonic transport in bulk systems. We apply this model to understand structural correlations in several recently synthesized perylene diimide (PDI)-based organic semiconductors. Twisted and nonplanar molecular architectures are found to promote amorphous morphologies while maintaining local π-stacking. A graph theoretical network analysis demonstrates that these twisted molecules are more likely to form percolating three-dimensional pathways for charge motion than strictly planar molecules, which show connectivity in only one dimension.

Polaron States in Fullerene Adducts Modeled by Coarse-Grained Molecular Dynamics and Tight Binding

Strong electron-phonon coupling leads to polaron localization in molecular semiconductor materials and influences charge transport, but it is expensive to calculate atomistically. Here, we propose a simple and efficient model to determine the energy and spatial extent of polaron states within a coarse-grained representation of a disordered molecular film. We calculate the electronic structure of the molecular assembly using a tight-binding Hamiltonian and determine the polaron state self-consistently by perturbing the site energies by the dielectric response of the surrounding medium to the charge. When applied to fullerene derivatives, the method shows that polarons extend over multiple molecules in C60 but localize on single molecules in higher adducts of phenyl-C61-butyric-acid-methyl-ester (PCBM) because of packing disorder and the polar side chains. In PCBM, polarons localize on single molecules only when energetic disorder is included or when the fullerene is dispersed in a blend. The method helps to establish the conditions under which a hopping transport model is justified.

Substrate-Dependent Morphology and Its Effect on Electrical Mobility of Doped Poly(3,4-ethylenedioxythiophene) (PEDOT) Thin Films

Deposition dynamics, crystallization, molecular packing, and electronic mobility of poly(3,4-ethylenedioxythiophene) (PEDOT) thin films are affected by the nature of the substrate. Computational microscopy has been carried out to reveal the morphology-substrate dependence for PEDOT thin films doped with molecular tosylate deposited on different substrates including graphite, Si₃N₄, silicon, and amorphous SiO₂. It is shown that the substrate is instrumental in formation of the lamellar structure. PEDOT films on the ordered substrates (graphite, Si₃N₄, and silicon) exhibit preferential face-on orientation, with graphite showing the most ordered and pronounced face-on packing. In contrast, PEDOT on amorphous SiO₂ exhibits the dominant edge-on orientation, except in the dry state where both packings are equally presented. The role of water and the porosity of the substrate in formation of the edge-on structure on SiO₂ is outlined. On the basis of the calculated morphology, the multiscale calculations of the electronic transport and percolative analysis are performed outlining how the character of the substrate affects the electron mobility. It is demonstrated that good crystallinity (PEDOT on graphite substrate) and high content of edge-on (PEDOT on SiO₂ substrate) are not enough to achieve the highest electrical in-plane mobility. Instead, the least ordered material with lower degree of the edge-on content (PEDOT on silicon substrate) provides the highest mobility because it exhibits an efficient network of π-π stacked chain extending throughout the entire sample.

Closely packed, low reorganization energy π-extended postfullerene acceptors for efficient polymer solar cells

New organic semiconductors are essential for developing inexpensive, high-efficiency, solution-processable polymer solar cells (PSCs). PSC photoactive layers are typically fabricated by film-casting a donor polymer and a fullerene acceptor blend, with ensuing solvent evaporation and phase separation creating discrete conduits for photogenerated holes and electrons. Until recently, n-type fullerene acceptors dominated the PSC literature; however, indacenodithienothiophene (IDTT)-based acceptors have recently enabled remarkable PSC performance metrics, for reasons that are not entirely obvious. We report two isomeric IDTT-based acceptors 3,9-bis(2-methylene-(3-(1,1-dicyanomethylene)-benz-(5, 6)indanone))-5,5,11,11-tetrakis(4-nonylphenyl)-dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']di-thiophene (ITN-C9) and 3,9-bis(2-methylene-(3-(1,1-dicyanomethylene)-benz(6,7)indanone))-5,5,11,11-tetrakis(4-nonylphenyl)-dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']dithiophene (ITzN-C9) that shed light on the exceptional IDTT properties vis-à-vis fullerenes. The neat acceptors and blends with fluoropolymer donor poly{[4,8-bis[5-(2- ethylhexyl)-4-fluoro-2-thienyl]benzo[1,2-b:4,5-b']dithiophene2,6-diyl]-alt-[2,5-thiophenediyl[5,7-bis(2-ethylhexyl)-4,8-dioxo4H,8H-benzo[1,2-c:4,5-c']dithiophene-1,3-diyl]]} (PBDB-TF) are investigated by optical spectroscopy, cyclic voltammetry, thermogravimetric analysis, differential scanning calorimetry, single-crystal X-ray diffraction, photovoltaic response, space-charge-limited current transport, atomic force microscopy, grazing incidence wide-angle X-ray scattering, and density functional theory-level quantum chemical analysis. The data reveal that ITN-C9 and ITzN-C9 organize such that the lowest unoccupied molecular orbital-rich end groups have intermolecular π-π distances as close as 3.31(1) Å, with electronic coupling integrals as large as 38 meV, and internal reorganization energies as small as 0.133 eV, comparable to or superior to those in fullerenes. ITN-C9 and ITzN-C9 have broad solar-relevant optical absorption, and, when blended with PBDB-TF, afford devices with power conversion efficiencies near 10%. Performance differences between ITN-C9 and ITzN-C9 are understandable in terms of molecular and electronic structure distinctions via the influences on molecular packing and orientation with respect to the electrode.

Electrothermal Transistor Effect and Cyclic Electronic Currents in Multithermal Charge Transfer Networks

A theory is developed to describe the coupled transport of energy and charge in networks of electron donor-acceptor sites which are seated in a thermally heterogeneous environment, where the transfer kinetics are dominated by Marcus-type hopping rates. It is found that the coupling of heat and charge transfer in such systems gives rise to exotic transport phenomena which are absent in thermally homogeneous systems and cannot be described by standard thermoelectric relations. Specifically, the directionality and extent of thermal transistor amplification and cyclical electronic currents in a given network can be controlled by tuning the underlying temperature gradient in the system. The application of these findings toward the optimal control of multithermal currents is illustrated on a paradigmatic nanostructure.

Visualizing excitations at buried heterojunctions in organic semiconductor blends

Interfaces play a crucial role in semiconductor devices, but in many device architectures they are nanostructured, disordered and buried away from the surface of the sample. Conventional optical, X-ray and photoelectron probes often fail to provide interface-specific information in such systems. Here we develop an all-optical time-resolved method to probe the local energetic landscape and electronic dynamics at such interfaces, based on the Stark effect caused by electron-hole pairs photo-generated across the interface. Using this method, we found that the electronically active sites at the polymer/fullerene interfaces in model bulk-heterojunction blends fall within the low-energy tail of the absorption spectrum. This suggests that these sites are highly ordered compared with the bulk of the polymer film, leading to large wavefunction delocalization and low site energies. We also detected a 100 fs migration of holes from higher- to lower-energy sites, consistent with these charges moving ballistically into more ordered polymer regions. This ultrafast charge motion may be key to separating electron-hole pairs into free charges against the Coulomb interaction.

Simulating charge transport in organic semiconductors and devices: a review

Charge transport simulation can be a valuable tool to better understand, optimise and design organic transistors (OTFTs), photovoltaics (OPVs), and light-emitting diodes (OLEDs). This review presents an overview of common charge transport and device models; namely drift-diffusion, master equation, mesoscale kinetic Monte Carlo and quantum chemical Monte Carlo, and a discussion of the relative merits of each. This is followed by a review of the application of these models as applied to charge transport in organic semiconductors and devices, highlighting in particular the insights made possible by modelling. The review concludes with an outlook for charge transport modelling in organic electronics.

Improved Scaling of Molecular Network Calculations: The Emergence of Molecular Domains

The design of materials needed for the storage, delivery, and conversion of (re)useable energy is still hindered by the lack of new, hierarchical molecular screening methodologies that encode information on more than one length scale. Using a molecular network theory as a foundation, we show that to describe charge transport in disordered materials the network methodology must be scaled-up. We detail the scale-up through the use of adjacency lists and depth first search algorithms for during operations on the adjacency matrix. We consider two types of electronic acceptors, perylenediimide (PDI) and the fullerene derivative phenyl-C61-butyric acid methyl ester (PCBM), and we demonstrate that the method is scalable to length scales relevant to grain boundary and trap formations. Such boundaries lead to a decrease in the percolation ratio of PDI with system size, while the ratio for PCBM remains constant, further quantifying the stable, diverse transport pathways of PCBM and its success as a charge-accepting material.

Charge transport network dynamics in molecular aggregates

Due to the nonperiodic nature of charge transport in disordered systems, generating insight into static charge transport networks, as well as analyzing the network dynamics, can be challenging. Here, we apply time-dependent network analysis to scrutinize the charge transport networks of two representative molecular semiconductors: a rigid n-type molecule, perylenediimide, and a flexible p-type molecule, [Formula: see text] Simulations reveal the relevant timescale for local transfer integral decorrelation to be [Formula: see text]100 fs, which is shown to be faster than that of a crystalline morphology of the same molecule. Using a simple graph metric, global network changes are observed over timescales competitive with charge carrier lifetimes. These insights demonstrate that static charge transport networks are qualitatively inadequate, whereas average networks often overestimate network connectivity. Finally, a simple methodology for tracking dynamic charge transport properties is proposed.

Influence of Coherent Tunneling and Incoherent Hopping on the Charge Transfer Mechanism in Linear Donor-Bridge-Acceptor Systems

The mechanism of charge transfer has been observed to change from tunneling to hopping with increasing numbers of DNA base pairs in polynucleotides and with the length of molecular wires. The aim of this paper is to investigate this transition by examining the population dynamics using a tight-binding Hamiltonian with model parameters to describe a linear donor-bridge-acceptor (D-B-A) system. The model includes a primary vibration and an electron-vibration coupling at each site. A further coupling of the primary vibration with a secondary phonon bath allows the system to dissipate energy to the environment and reach a steady state. We apply the quantum master equation (QME) approach, based on second-order perturbation theory in a quantum dissipative system, to examine the dynamical processes involved in charge-transfer and follow the population transfer rate at the acceptor, ka, to shed light on the transition from tunneling to hopping. With a small tunneling parameter, V, the on-site population tends to localize and form polarons, and the hopping mechanism dominates the transfer process. With increasing V, the population tends to be delocalized and the tunneling mechanism dominates. The competition between incoherent hopping and coherent tunneling governs the mechanism of charge transfer. By varying V and the total number of sites, we also examine the onset of the transition from tunneling to hopping with increasing length.

Molecular helices as electron acceptors in high-performance bulk heterojunction solar cells

Despite numerous organic semiconducting materials synthesized for organic photovoltaics in the past decade, fullerenes are widely used as electron acceptors in highly efficient bulk-heterojunction solar cells. None of the non-fullerene bulk heterojunction solar cells have achieved efficiencies as high as fullerene-based solar cells. Design principles for fullerene-free acceptors remain unclear in the field. Here we report examples of helical molecular semiconductors as electron acceptors that are on par with fullerene derivatives in efficient solar cells. We achieved an 8.3% power conversion efficiency in a solar cell, which is a record high for non-fullerene bulk heterojunctions. Femtosecond transient absorption spectroscopy revealed both electron and hole transfer processes at the donor-acceptor interfaces. Atomic force microscopy reveals a mesh-like network of acceptors with pores that are tens of nanometres in diameter for efficient exciton separation and charge transport. This study describes a new motif for designing highly efficient acceptors for organic solar cells.

A Simple Index for Characterizing Charge Transport in Molecular Materials

While advances in quantum chemistry have rendered the accurate prediction of band alignment relatively straightforward, the ability to forecast a noncrystalline, multimolecule system's conductivity possesses no simple computational form. Adapting the theory of classical resistor networks, we develop an index for quantifying charge transport in bulk molecular materials, without the requirement of crystallinity. The basic behavior of this index is illustrated through its application to simple lattices and clusters of common organic photovoltaic molecules, where it is shown to reproduce experimentally known performances for these materials. This development provides a quantitative computational means for determining a priori the bulk charge transport properties of molecular materials.

A Comparison of Charge Separation Dynamics in Organic Blend Films Employing Fullerene and Perylene Diimide Electron Acceptors

We report a comparison of charge carrier dynamics and device performance for low band gap polymer PBDTTT-CT in blends with the fullerene acceptor PC71BM and a PDI derivative with similar electron affinities. Charge separation and recombination dynamics are found to be remarkably similar for these two acceptors, with both blends exhibiting efficient, ultrafast charge separation (time constants of 1.6 and 1.4 ps, respectively). The lower device performance for the PDI acceptor (1.75% compared to 3.5% for the equivalent PC71BM device) is shown to result from slower charge transport, increasing nongeminate recombination losses during charge collection.