Can computational approaches aid in untangling the inherent complexity of practical organic photovoltaic systems?

The correspondence between the conformational and chromophoric properties of amorphous conjugated polymers in mesoscale condensed systems

For π-conjugated polymers, the notion of spectroscopic units or "chromophores" provides illuminating insights into the experimentally observed absorption/emission spectra and the mechanisms of energy/charge transfer. To date, however, no statistical analysis has revealed a direct correspondence between chromophoric and conformational properties-with the latter being fundamental to polymer semiconductors. Herein, we propose a "persistence length" calculation to re-evaluate chain conformation over a full conjugation length. The mesoscale condensed systems of MEH-PPV and MEH-PPV/C₆₀ hybrid (system size ∼10 × 10 × 10 nm³) are utilized as two prototypical model systems, along with a full range of segmental lengths (2-20-mer) and five lowest singlet excited states to hint at the generality of the features presented. We demonstrate, for the first time, that two properly re-defined conformational factors that characterize chain folding and planarity, respectively, capture excellently the population distribution of chromophores in both systems investigated. In contrast, the conventional strategy of utilizing two adjacent monomer units to characterize (local) chain conformation results in only an inconspicuous correlation between the two, as previously reported. It is further shown that chain folding-and not chain planarity-is more relevant in capturing the associated oscillator strength for the first excited state, where the transient dipole moments are known to align with the chain conformation, although the corresponding excitation energy and exciton size seem relatively unaffected. The observed effects of C₆₀ on the MEH-PPV adsorption spectra also agree with recent experimental trends. Overall, the present findings are expected to aid future multiscale computer simulations and spectroscopy-data interpretations for polymer semiconductors and their hybrid systems.

Multiscale Molecular Simulation of Solution Processing of SMDPPEH: PCBM Small-Molecule Organic Solar Cells

Solution-processed small-molecule organic solar cells are a promising renewable energy source because of their low production cost, mechanical flexibility, and light weight relative to their pure inorganic counterparts. In this work, we developed a coarse-grained (CG) Gay-Berne ellipsoid molecular simulation model based on atomistic trajectories from all-atom molecular dynamics simulations of smaller system sizes to systematically study the nanomorphology of the SMDPPEH/PCBM/solvent ternary blend during solution processing, including the blade-coating process by applying external shear to the solution. With the significantly reduced overall system degrees of freedom and computational acceleration from GPU, we were able to go well beyond the limitation of conventional all-atom molecular simulations with a system size on the order of hundreds of nanometers with mesoscale molecular detail. Our simulations indicate that, similar to polymer solar cells, the optimal blending ratio in small-molecule organic solar cells must provide the highest specific interfacial area for efficient exciton dissociation, while retaining balanced hole/electron transport pathway percolation. We also reveal that blade-coating processes have a significant impact on nanomorphology. For given donor/acceptor blending ratios, applying an external shear force can effectively promote donor/acceptor phase segregation and stacking in the SMDPPEH domains. The present study demonstrated the capability of an ellipsoid-based coarse-grained model for studying the nanomorphology evolution of small-molecule organic solar cells during solution processing/blade-coating and provided links between fabrication protocols and device nanomorphologies.

Petascale Simulations of the Morphology and the Molecular Interface of Bulk Heterojunctions

Understanding how additives interact and segregate within bulk heterojunction (BHJ) thin films is critical for exercising control over structure at multiple length scales and delivering improvements in photovoltaic performance. The morphological evolution of poly(3-hexylthiophene) (P3HT) and phenyl-C61-butyric acid methyl ester (PCBM) blends that are commensurate with the size of a BHJ thin film is examined using petascale coarse-grained molecular dynamics simulations. Comparisons between two-component and three-component systems containing short P3HT chains as additives undergoing thermal annealing demonstrate that the short chains alter the morphology in apparently useful ways: they efficiently migrate to the P3HT/PCBM interface, increasing the P3HT domain size and interfacial area. Simulation results agree with depth profiles determined from neutron reflectometry measurements that reveal PCBM enrichment near substrate and air interfaces but a decrease in that PCBM enrichment when a small amount of short P3HT chains are integrated into the BHJ blend. Atomistic simulations of the P3HT/PCBM blend interfaces show a nonmonotonic dependence of the interfacial thickness as a function of number of repeat units in the oligomeric P3HT additive, and the thiophene rings orient parallel to the interfacial plane as they approach the PCBM domain. Using the nanoscale geometries of the P3HT oligomers, LUMO and HOMO energy levels calculated by density functional theory are found to be invariant across the donor/acceptor interface. These connections between additives, processing, and morphology at all length scales are generally useful for efforts to improve device performance.

Multiscale Simulation and Modeling of Multilayer Heteroepitactic Growth of C60 on Pentacene

We apply multiscale methods to describe the strained growth of multiple layers of C60 on a thin film of pentacene. We study this growth in the presence of a monolayer pentacene step to compare our simulations to recent experimental studies by Breuer and Witte of submonolayer growth in the presence of monolayer steps. The molecular-level details of this organic semiconductor interface have ramifications on the macroscale structural and electronic behavior of this system and allow us to describe several unexplained experimental observations for this system. The growth of a C60 thin film on a pentacene surface is complicated by the differing crystal habits of the two component species, leading to heteroepitactical growth. In order to probe this growth, we use three computational methods that offer different approaches to coarse-graining the system and differing degrees of computational efficiency. We present a new, efficient reaction-diffusion continuum model for 2D systems whose results compare well with mesoscale kinetic Monte Carlo (KMC) results for submonolayer growth. KMC extends our ability to simulate multiple layers but requires a library of predefined rates for event transitions. Coarse-grained molecular dynamics (CGMD) circumvents KMC's need for predefined lattices, allowing defects and grain boundaries to provide a more realistic thin film morphology. For multilayer growth, in this particularly suitable candidate for coarse-graining, CGMD is a preferable approach to KMC. Combining the results from these three methods, we show that the lattice strain induced by heteroepitactical growth promotes 3D growth and the creation of defects in the first monolayer. The CGMD results are consistent with experimental results on the same system by Conrad et al. and by Breuer and Witte in which C60 aggregates change from a 2D structure at low temperature to 3D clusters along the pentacene step edges at higher temperatures.

The prediction of the morphology and PCE of small molecular organic solar cells

The predicted morphology, domain size, PCE (power conversion efficiency) of Small Molecular Organic Solar Cells.

Structure and optoelectronic properties of helical pyridine–furan, pyridine–pyrrole and pyridine–thiophene oligomers

The optoelectronic properties of helical structures pyridine–furan, pyridine–pyrrole and pyridine–thiophene oligomers are explored using density functional theory methods.

Interfacial properties and design of functional energy materials

CONSPECTUS: The vital importance of energy to society continues to demand a relentless pursuit of energy responsive materials that can bridge fundamental chemical structures at the molecular level and achieve improved functionality and performance. This demand can potentially be realized by harnessing the power of self-assembly, a spontaneous process where molecules or much larger entities form ordered aggregates as a consequence of predominately noncovalent (weak) interactions. Self-assembly is the key to bottom-up design of molecular devices, because the nearly atomic-level control is very difficult to realize in a top-down, for example, lithographic, approach. However, while function in simple systems such as single crystals can often be evaluated a priori, predicting the function of the great variety of self-assembled molecular architectures is complicated by the lack of understanding and control over nanoscale interactions, mesoscale architectures, and macroscale order. To establish a foundation toward delivering practical solutions, it is critical to develop an understanding of the chemical and physical mechanisms responsible for the self-assembly of molecular and hybrid materials on various support substrates. Typical molecular self-assembly involves noncovalent intermolecular and substrate-molecule interactions. These interactions remain poorly understood, due to the combination of many-body interactions compounded by local or collective influences from the substrate atomic lattice and electronic structure. Progress toward unraveling the underlying physicochemical processes that control the structure and macroscopic physical, chemical, mechanical, electrical, and transport properties of materials increasingly requires tight integration of theory, modeling, and simulation with precision synthesis, advanced experimental characterization, and device measurements. Theory, modeling, and simulation can accelerate the process of materials understanding and design by providing atomic level understanding of the underlying physicochemical phenomena (illuminating connections between experiments). It can also provide the ability to explore new materials and conditions before they are realized in the laboratory. With tight integration and feedback with experiment, it becomes feasible to identify promising materials or processes for targeted energy applications. In this Account, we highlight recent advances and success in using an integrated approach based on electronic structure simulations and scanning probe microscopy techniques to study and design functional materials formed from the self-assembly of molecules into supramolecular or polymeric architectures on substrates.

Efficient Computational Screening of Organic Polymer Photovoltaics

There has been increasing interest in rational, computationally driven design methods for materials, including organic photovoltaics (OPVs). Our approach focuses on a screening "pipeline", using a genetic algorithm for first stage screening and multiple filtering stages for further refinement. An important step forward is to expand our diversity of candidate compounds, including both synthetic and property-based measures of diversity. For example, top monomer pairs from our screening are all donor-donor (D-D) combinations, in contrast with the typical donor-acceptor (D-A) motif used in organic photovoltaics. We also find a strong "sequence effect", in which the average HOMO-LUMO gap of tetramers changes by ∼0.2 eV as a function of monomer sequence (e.g., ABBA versus BAAB); this has rarely been explored in conjugated polymers. Beyond such optoelectronic optimization, we discuss other properties needed for high-efficiency organic solar cells, and applications of screening methods to other areas, including non-fullerene n-type materials, tandem cells, and improving charge and exciton transport.