Steric control of the donor/acceptor interface: implications in organic photovoltaic charge generation.

The performance of organic photovoltaic (OPV) devices is currently limited by modest short-circuit current densities. Approaches toward improving this output parameter may provide new avenues to advance OPV technologies and the basic science of charge transfer in organic semiconductors. This work highlights how steric control of the charge separation interface can be effectively tuned in OPV devices. By introducing an octylphenyl substituent onto the investigated polymer backbones, the thermally relaxed charge-transfer state, and potentially excited charge-transfer states, can be raised in energy. This decreases the barrier to charge separation and results in increased photocurrent generation. This finding is of particular significance for nonfullerene OPVs, which have many potential advantages such as tunable energy levels and spectral breadth, but are prone to poor exciton separation efficiencies. Computational, spectroscopic, and synthetic methods were combined to develop a structure-property relationship that correlates polymer substituents with charge-transfer state energies and, ultimately, device efficiencies.

Molecular understanding of organic solar cells: the challenges.

Our objective in this Account is 3-fold. First, we provide an overview of the optical and electronic processes that take place in a solid-state organic solar cell, which we define as a cell in which the semiconducting materials between the electrodes are organic, be them polymers, oligomers, or small molecules; this discussion is also meant to set the conceptual framework in which many of the contributions to this Special Issue on Photovoltaics can be viewed. We successively turn our attention to (i) optical absorption and exciton formation, (ii) exciton migration to the donor-acceptor interface, (iii) exciton dissociation into charge carriers, resulting in the appearance of holes in the donor and electrons in the acceptor, (iv) charge-carrier mobility, and (v) charge collection at the electrodes. For each of these processes, we also describe the theoretical challenges that need to be overcome to gain a comprehensive understanding at the molecular level. Finally, we highlight recent theoretical advances, in particular regarding the determination of the energetics and dynamics at organic-organic interfaces, and underline that the right balance needs to be found for the optimization of material parameters that often result in opposite effects on the photovoltaic performance.

Polarization energies in oligoacene semiconductor crystals.

Characterization of the electronically polarized environment and the nuclear relaxation that accompanies charge carriers is fundamental to charge transport in crystalline, polycrystalline, and amorphous organic solids. To study the polarization effects of localized charged carriers, we use quantum/classical QM/MM approaches with charge redistribution and polarizable force field schemes and apply them to crystals of naphthalene through pentacene. We describe the results of a comprehensive investigation of the electronic polarization energies in molecular crystal structures of these oligoacenes and discuss as well the evolution of the nuclear relaxation energies calculated for model oligoacene systems.

Experimental and theoretical study of temperature dependent exciton delocalization and relaxation in anthracene thin films.

The spectroscopy of solid anthracene is examined both experimentally and theoretically. To avoid experimental complications such as self-absorption and polariton effects, ultrathin polycrystalline films deposited on transparent substrates are studied. To separate the contributions from different emitting species, the emission is resolved in both time and wavelength. The spectroscopic data are interpreted in terms of a three-state kinetic model, where two excited states, a high energy state 1 and a low energy state 2, both contribute to the luminescence and are kinetically coupled. Using this model, we analyze the spectral lineshape, relative quantum yield, and relaxation rates as a function of temperature. For state 1, we find that the ratio of the 0-0 vibronic peak to the 0-1 peak is enhanced by roughly a factor of 3.5 at low temperature, while the quantum yield and decay rates also increase by a similar factor. These observations are explained using a theoretical model previously developed for herringbone polyacene crystals. The early-time emission lineshape is consistent with that expected for a linear aggregate corresponding to an edge-dislocation defect. The results of experiment and theory are quantitatively compared at different temperatures in order to estimate that the singlet exciton in our polycrystalline films is delocalized over about ten molecules. Within these domains, the exciton's coherence length steadily increases as the temperature drops, until it reaches the limits of the domain, whereupon it saturates and remains constant as the temperature is lowered further. While the theoretical modeling correctly reproduces the temperature dependence of the fluorescence spectral lineshape, the decay of the singlet exciton appears to be determined by a trapping process that becomes more rapid as the temperature is lowered. This more rapid decay is consistent with accelerated trapping due to increased delocalization of the exciton at lower temperatures. These observations suggest that exciton coherence can play an important role in both radiative and nonradiative decay channels in these materials. Our results show that the spectroscopy of polyacene solids can be analyzed in a self-consistent fashion to obtain information about electronic delocalization and domain sizes.

Theoretical characterization of titanyl phthalocyanine as a p-type organic semiconductor: short intermolecular pi-pi interactions yield large electronic couplings and hole transport bandwidths.

The charge-transport properties of the triclinic phase II crystal of titanyl phthalocyanine (alpha-TiOPc) are explored within both a hopping and bandlike regime. Electronic coupling elements in convex- and concave-type dimers are calculated using density functional theory, and the relationship between molecular structure and crystal packing structure in model dimer configurations is considered. Hole transport bandwidths derived from crystal structure dimers are compared to those obtained from electronic band structure calculations; very good agreement between the two approaches is found. The calculations predict large hole bandwidths, on the order of 0.4 eV, and correspondingly very low hole reorganization energies.

On the mechanism of peripentacene formation from pentacene: computational studies of a prototype for graphene formation from smaller acenes.

The formation of peripentacene during the high-temperature vacuum sublimation of pentacene (P) in the presence of trace amounts of 6,13-dihydropentacene (DHP) has been studied computationally with density functional theory. Computational and kinetic analyses indicate that competing mechanisms involving a series of H atom transfers initiated by hydrogen transfer from DHP to P can account for the formation of peripentacene. The overall reaction is predicted to proceed with a free energy barrier of 36.1 kcal/mol and to be autocatalytic. Kinetic modeling supports the proposed mechanism.

Why 6-methylpentacene deconjugates but avoids the thermally allowed unimolecular mechanism.

[reaction: see text] Tautomeric equilibria involving hydrogen migration in methylacenes were explored computationally using DFT methods. As the aromatic system becomes more extended, the methylene isomer is predicted to become favored. Reasonable-looking pericyclic sigmatropic hydrogen shifts are found to be energetically prohibitive, and bimolecular mechanisms involving radical pair intermediates are energetically feasible.

Lone-pair orbital interactions in polythiaadamantanes.

The electronic structures of a series of polythiaadamantanes from thiaadamantane through 2,4,6,8,9,10-hexathiaadamantane (HTA) have been analyzed using density functional theory calculations in conjunction with Hückel and natural bond orbital analysis. The effects of multiple sulfur p-type lone-pair orbital interactions on ionization potentials, hole mobilities, and electronic coupling have been determined. An overall increase in the average energy of the lone-pair orbitals as the number of sulfur atoms increases is predicted, with the exact positioning of the HOMO depending on specific lone-pair interactions. Separation of through-bond (TB) and through-space (TS) interactions between intramolecular sulfur atoms has been performed using localized molecular orbitals and model systems based on interacting hydrogen sulfide molecules. TB interations were found to reduce orbital splitting, while TS interactions were found to increase orbital splitting. TS interactions were more or less constant from one polythiaadamantane to the next, and the contributions of TB effects to individual orbital energies vary depending on the relative orientation of sulfur atoms as determined by the sigma molecular framework. Electronic coupling between intermolecular sulfur lone-pair orbitals was determined by investigating unique dimer pairs observed in the crystal structure of HTA. Electronic coupling is not as strong as expected given the short intermolecular S-S distances observed in the crystal structure. In general, B3LYP/6-31G(d) and B3LYP/6-311+G(d,p) give very similar orbital energies and splittings.

Theoretical studies of quantum amplified isomerizations for imaging systems involving hexamethyl Dewar benzene and related systems.

The ring-opening reactions of the radical cations of hexamethyl Dewar benzene (1) and Dewar benzene have been studied using density functional theory (DFT) and complete active-space self-consistent field (CASSCF) calculations. Compound 1 is known to undergo photoinitiated ring opening by a radical cation chain mechanism, termed "quantum amplified isomerization" (QAI), which is due to the high quantum yield. Why QAI is efficient for 1 but not other reactions is explained computationally. Two radical cation minima of 1 and transition states located near avoided crossings are identified. The state crossings are characterized by conical intersections corresponding to degeneracy between doublet surfaces. Ring opening occurs by formation of the radical cation followed by a decrease in the flap dihedral angle. A rate-limiting Cs transition state leads to a second stable radical cation with an elongated transannular C-C bond and an increased flap dihedral. This structure proceeds through a conrotatory-like pathway of Cs symmetry to give the benzene radical cation. The role of electron transfer was investigated by evaluating oxidation of various systems using adiabatic ionization energies and electron affinities calculated from neutral and cation geometries. Electron-transfer theory was applied to 1 to investigate the limiting effects of back-electron transfer as it is related to the unusual stability of the two radical cations. Expected changes in optical properties between reactants and products of Dewar benzene compounds and other systems known to undergo QAI were characterized by computing frequency-dependent indices of refraction from isotropic polarizabilities. In particular, the reaction of 1 shows greater contrast in index of refraction than that of the Dewar benzene parent system.

Electronic structures and properties of twisted polyacenes.

The effects of twisting on the electronic structures and properties of polyacenes were studied computationally using DFT methods. Singlet-triplet and HOMO-LUMO gaps and vertical S0-S1 transition energies are marginally affected as a function of end-to-end twist angle. The large twist induced by the phenyl substituents in 9,10,11,20,21,22-hexaphenyltetrabenzo[a,c,l,n]pentacene has little influence on its electronic structure.

DFT study of pericyclic reaction cascades in the synthesis of antibiotic TAN-1085.

DFT calculations show that aromatic and bis-methoxy substituent effects in a synthetic precursor of TAN-1085 strongly favor a [1,7] sigmatropic hydrogen shift over the 6 pi electrocyclic pathway, rendering the latter unfavorable in synthesis.

Impact of multiple cation-pi interactions upon calix[4]arene substrate binding and specificity.

The cation-pi interaction influence on the conformation and binding of calix[4]arenes to alkali-metal cations has been studied using a dehydroxylated model. The model allows for the separation of cooperative cation-pi and electrostatic forces commonly found in the binding motifs found in calixarene complexes. Starting from the four well-known calix[4]arene conformations, six conformers for this dehydroxylated model (cone, partial cone, flattened cone, chair, 1,2-alternate, and 1,3-alternate) have been characterized by geometry optimization and frequency analysis using the Becke three-parameter exchange functional with the nonlocal correlation functional of Lee, Yang, and Parr and the 6-31G(d) basis set. Without the stabilization provided by the hydroxyl hydrogen bonds in calix[4]arene, neither the cone nor the 1,2-alternate conformation is computed to be a ground-state structure. The partial cone, flattened cone, chair, and 1,3-alternate conformers have been identified as ground-state structures in a vacuum, with the partial cone and the 1,3-alternate as the lowest energy minima in the aromatic model. The C(4)(v)() cone conformation is found to be a transition structure separating the flattened cone (C(2)(v)()) conformers. The energetic and structural preferences of the calix[4]arene model change dramatically when it is bound to Li(+), Na(+), and K(+). The number of pi-faces, the positioning of these pi-faces with respect to the cations, and the nature of the cation were studied as factors in the binding strength. A detailed study of the distances and angles between the aromatic ring centroids and the cations reveals the energetic advantages of multiple weak cation-pi interactions. The geometries are often far from the optimal cation-pi interaction in which the cation approaches in a perpendicular path the aromatic ring center, where the quadrupole moment is strongest. The results reveal that multiple weaker nonoptimal cation-pi interactions contribute significantly to the overall binding strength. This theoretical analysis underscores the importance of neighboring aromatic faces and provides new insight into the significance of cation-pi binding, not only for calix[4]arenes, but also for other supramolecular and biological systems.

Electrostatic interactions between substituents as regioselectivity control elements in Diels-Alder cycloadditions. A DFT study of cycloadditions of 1-methoxy-4-trimethylsiloxy dienes with acrylonitrile.

The regioselectivities of Diels-Alder reactions of 1-methoxy-4-trimethylsiloxy-1,3-butadiene and the corresponding o-xylylene with acrylonitrile were explored with density functional theory. The transition state of the reaction of the diene with acrylonitrile was studied in both the gas phase and in a low dielectric (epsilon = 2.247) solvent. The regioselectivities of these reactions are predicted to be controlled by the direct electrostatic interactions between the diene and dienophile substituents. The electron-donating capacities of methoxy and trimethylsiloxy substituents are shown to be very similar and to not contribute to regioselectivity control. A prediction is made that the cycloadditions of the o-xylylene will be unselective, while 1-methoxy-4-trimethylsiloxy-1,3-butadiene will give a small preference ( approximately 5:1) in favor of the proximal methoxy and cyano groups. The thermochromic behavior of trans-disubstituted disiloxy benzocyclobutene was also explored.