Simulation of vapor-phase deposition and growth of a pentacene thin film on C60 (001)

Kinetic Trapping of Charge-Transfer Molecules at Metal Interfaces

Despite the common expectation that conjugated organic molecules on metals adsorb in a flat-lying layer, several recent studies have found coverage-dependent transitions to upright-standing phases, which exhibit notably different physical properties. In this work, we argue that from an energetic perspective, thermodynamically stable upright-standing phases may be more common than hitherto thought. However, for kinetic reasons, this phase may often not be observed experimentally. Using first-principles kinetic Monte Carlo simulations, we find that the structure with lower molecular density is (almost) always formed first, reminiscent of Ostwald's rule of stages. The phase transitions to the upright-standing phase are likely to be kinetically hindered under the conditions typically used in surface science. The simulation results are experimentally confirmed for the adsorption of tetracyanoethylene on Cu(111) using infrared and X-ray photoemission spectroscopy. Investigating both the role of the growth conditions and the energetics of the interface, we find that the time for the phase transition is determined mostly by the deposition rate and, thus, is mostly independent of the nature of the molecule.

PyThinFilm: Automated Molecular Dynamics Simulation Protocols for the Generation of Thin Film Morphologies

The performance of organic optoelectronic devices, such as organic light-emitting diodes (OLEDs) and organic solar cells (OSCs), is intrinsically related to the molecular-scale morphology of the thin films from which they are composed. However, the experimental characterization of morphology at the molecular level is challenging due to the often amorphous or at best semicrystalline nature of these films. Classical molecular modeling techniques, such as molecular dynamics (MD) simulation, are increasingly used to understand the relationship between morphology and the properties of thin-film devices. PyThinFilm (github.com/ATB-UQ/PyThinFilm) is an open-source Python package which allows fully automated MD simulations of thin film growth to be performed using vacuum and/or solution deposition processes. PyThinFilm utilizes the GROMACS simulation package in combination with interaction parameters from the Automated Topology Builder (atb.uq.edu.au). Here, PyThinFilm is described along with an overview of applications in which PyThinFilm has been used to study the thin films of organic semiconductor materials typically used in OLEDs and OSCs.

Interplay of orientational order and roughness in simulated thin film growth of anisotropically interacting particles

Roughness and orientational order in thin films of anisotropic particles are investigated using kinetic Monte Carlo simulations on a cubic lattice. Anisotropic next-neighbor interactions between the lattice particles were chosen to mimic the effects of shape anisotropy in the interactions of disk- or rodlike molecules with van der Waals attractions. Increasing anisotropy leads first to a preferred orientation in the film (which is close to the corresponding equilibrium transition) while the qualitative mode of roughness evolution (known from isotropic systems) does not change. At strong anisotropies, an effective step-edge (Ehrlich-Schwoebel) barrier appears and a nonequilibrium roughening effect is found, accompanied by reordering in the film which can be interpreted as the nucleation and growth of domains of lying-down disks or rods. The information on order and roughness is combined into a diagram of dynamic growth modes.

A simple predictor of interface orientation of fluids of disk-like anisotropic particles and its implications for organic semiconductors

From classical molecular dynamics simulations, we identify a simple and general predictor of molecular orientation at solid and vapour interfaces of isotropic fluids of disk-like anisotropic particles based on their shape and interaction anisotropy. For a wide variety of inter-particle interactions, temperatures, and substrate types within the range of typical organic semiconductors and their processing conditions, we find remarkable universal scaling of the orientation at the interface with the free energy calculated from pair interactions between close-packed nearest neighbours and an empirically derived universal relationship between the entropy and the shape anisotropy and bulk volume fraction of the fluid particles. The face-on orientation of fluid particles at the solid interface is generally predicted to be the equilibrium structure, although the alignment can be controlled by tuning the particle shape and substrate type, while changing the strength of fluid-fluid interactions is likely to play a less effective role. At the vapour interface, only the side-on structure is predicted, and conditions for which the face-on structure may be preferred, such as low temperature, low interaction anisotropy, or low shape anisotropy, are likely to result in little orientation preference (due to the low anisotropy) or be associated with a phase transition to an anisotropic bulk phase for systems with interactions in the range of typical organic semiconductors. Based on these results, we propose a set of guidelines for the rational design and processing of organic semiconductors to achieve a target orientation at a solid or vapour interface.

Toward Targeted Kinetic Trapping of Organic–Inorganic Interfaces: A Computational Case Study

Properties of inorganic–organic interfaces, such as their interface dipole, strongly depend on the structural arrangements of the organic molecules. A prime example is tetracyanoethylene (TCNE) on Cu(111), which shows two different phases with significantly different work functions. However, the thermodynamically preferred phase is not always the one that is best suited for a given application. Rather, it may be desirable to selectively grow a kinetically trapped structure. In this work, we employ density functional theory and transition state theory to discuss under which conditions such a kinetic trapping might be possible for the model system of TCNE on Cu. Specifically, we want to trap the molecules in the first layer in a flat-lying orientation. This requires temperatures that are sufficiently low to suppress the reorientation of the molecules, which is thermodynamically more favorable for high dosages, but still high enough to enable ordered growth through diffusion of molecules. On the basis of the temperature-dependent diffusion and reorientation rates, we propose a temperature range at which the reorientation can be successfully suppressed.

Topics in the mathematical design of materials

We present a perspective on several current research directions relevant to the mathematical design of new materials. We discuss: (i) design problems for phase-transforming and shape-morphing materials, (ii) epitaxy as an approach of central importance in the design of advanced semiconductor materials, (iii) selected design problems in soft matter, (iv) mathematical problems in magnetic materials, (v) some open problems in liquid crystals and soft materials and (vi) mathematical problems on liquid crystal colloids. The presentation combines topics from soft and hard condensed matter, with specific focus on those design themes where mathematical approaches could possibly lead to exciting progress. This article is part of the theme issue 'Topics in mathematical design of complex materials'.

Lattice gas study of thin-film growth scenarios and transitions between them: Role of substrate

Thin-film growth is investigated in two types of lattice gas models where substrate and film particles are different, expressed by unequal interaction energy parameters. The first is of solid-on-solid type, whereas the second additionally incorporates desorption, diffusion in the gas phase above the film and readsorption at the film (appropriate for growth in colloidal systems). In both models, the difference between particle-substrate and particle-particle interactions plays a central role for the evolution of the film morphology at intermediate times. The models exhibit a dynamic layering transition which occurs at generally lower substrate attraction strengths than the equilibrium layering transition. A second, flattening transition is found where initial island growth transforms to layer-by-layer growth at intermediate deposition times. Combined with the known roughening behavior in such models for very large deposition times, we present four global growth scenarios, charting out the possible types of roughness evolution.

Fate of Low-Lying Charge-Transfer Excited States in a Donor:Acceptor Blend with a Large Energy Offset

In an effort to gain a comprehensive picture of the interfacial states in bulk heterojunction solar cells, we provide a combined experimental-theoretical analysis of the energetics and dynamics of low-lying electronic charge-transfer (CT) states in donor:acceptor blends with a large frontier orbital energy offset. By varying the blend composition and temperature, we unravel the static and dynamic contributions to the disordered density of states (DOS) of the CT-state manifold and assess their recombination to the ground state. Namely, we find that static disorder (conformational and electrostatic) shapes the CT DOS and that fast nonradiative recombination crops the low-energy tail of the distribution probed by external quantum efficiency (EQE) measurements (thereby largely contributing to voltage losses). Our results then question the standard practice of extracting microscopic parameters such as exciton energy and energetic disorder from EQE.

Evolution and Morphology of Thin Films Formed by Solvent Evaporation: An Organic Semiconductor Case Study

The crucial role played by the solution-vapor interface in determining the growth and morphology of an organic semiconductor thin film formed by solvent evaporation has been examined in atomic detail. Specifically, how the loss of individual solvent molecules from the surface of the solution induces solute assembly has been studied using molecular dynamics simulations. The system consisted of bis(2-phenylpyridine) (acetylacetonate)iridium(III) [Ir(ppy)₂(acac)] and 4,4'-bis(N-carbazolyl)-1,1'-biphenyl (CBP) in chloroform at 310 K. The simulations clearly indicate that (a) the system does not undergo uniform phase separation (spinodal decomposition), (b) solute aggregation initiates at the solution-vapor interface, (c) the distribution of solvent in the film is nonhomogeneous, (d) this nonhomogeneous distribution can induce preferential alignment of host molecules, and (e) a portion of the solvent likely remains trapped within the film. The work not only demonstrates the ability to directly model evaporation in atomic detail on the relevant length scales but also shows that atomistic simulations have the potential to shed new light on morphological properties of a wide range of organic semiconductor devices manufactured using solution-processing methods.

Impact of Nanomorphology on Surface Doping of Organic Semiconductors: The Pentacene-C60F48 Interface

Establishing the rather complex correlation between the structure and the charge transfer in organic-organic heterostructures is of utmost importance for organic electronics and requires spatially resolved structural, chemical, and electronic details. Insight into this issue is provided here by combining atomic force microscopy, Kelvin probe force microscopy, photoemission electron microscopy, and low-energy electron microscopy for investigating a case study. We select the interface formed by pentacene (PEN), benchmark among the donor organic semiconductors, and a p-type dopant from the family of fluorinated fullerenes. As for Buckminsterfullerene (C₆₀), the growth of its fluorinated derivative C₆₀F₄₈ is influenced by the thickness and crystallinity of the PEN buffer layer, but the behavior is markedly different. We provide a microscopic description of the C₆₀F₄₈/PEN interface formation and analyze the consequences in the electronic properties of the final heterostructure. For just one single layer of PEN, a laterally complete but noncompact C₆₀F₄₈/PEN interface is created, importantly affecting the surface work function. Nonetheless, from the very beginning of the second layer formation, the presence of epitaxial and nonepitaxial PEN domains dramatically influences the growth dynamics and extremely well packed two-dimensional C₆₀F₄₈ islands develop. Insightful elemental maps of the C₆₀F₄₈/PEN surface spatially resolve the nonuniform distribution of the dopant molecules, which leads to a heterogeneous work function landscape.

The interplay of interfaces, supramolecular assembly, and electronics in organic semiconductors

Organic semiconductors, which include a diverse range of carbon-based small molecules and polymers with interesting optoelectronic properties, offer many advantages over conventional inorganic semiconductors such as silicon and are growing in importance in electronic applications. Although these materials are now the basis of a lucrative industry in electronic displays, many promising applications such as photovoltaics remain largely untapped. One major impediment to more rapid development and widespread adoption of organic semiconductor technologies is that device performance is not easily predicted from the chemical structure of the constituent molecules. Fundamentally, this is because organic semiconductor molecules, unlike inorganic materials, interact by weak non-covalent forces, resulting in significant structural disorder that can strongly impact electronic properties. Nevertheless, directional forces between generally anisotropic organic-semiconductor molecules, combined with translational symmetry breaking at interfaces, can be exploited to control supramolecular order and consequent electronic properties in these materials. This review surveys recent advances in understanding of supramolecular assembly at organic-semiconductor interfaces and its impact on device properties in a number of applications, including transistors, light-emitting diodes, and photovoltaics. Recent progress and challenges in computer simulations of supramolecular assembly and orientational anisotropy at these interfaces is also addressed.

Controlling Interface Morphology and Layer Crystallinity in Organic Heterostructures: Microscopic View on C60 Island Formation on Pentacene Buffer Layers

Controlling the crystallinity of organic thin films is an important aspect in the improvement of organic electronic devices. However, because of high molecular mass, structural anisotropy, and weak intermolecular van der Waals bonding, crystalline ordering is not easily accomplished. While film preparation at elevated substrate temperature often improves the crystalline quality, this approach cannot be applied to temperature-sensitive materials such as plastic foils used as substrates for flexible electronics. Here, we examine in detail a low-temperature approach to improve film crystallinity by using ultrathin pentacene (PEN) buffer layers that allow crystalline growth of buckminsterfullerene (C₆₀) thin films while without such buffer layers, only amorphous fullerene films are formed upon room-temperature deposition on various support substrates. Remarkably, this effect depends critically on the thickness of the PEN buffer and requires a thickness of at least two monolayers to induce crystalline growth, whereas a buffer layer consisting of a monolayer of PEN again yields amorphous C₆₀ films. Combining crystallographic investigations by X-ray diffraction and atomic force microscopy measurements, we determine distinct nucleation sites on buffer layers of different thickness, which are correlated to the amorphous, respectively crystalline C₆₀ islands. Our microscopic analysis reveals distinct differences for the nucleation and diffusivity of fullerenes on the PEN monolayer and on thicker buffer layers, which are attributed to the molecular arrangement in the PEN monolayer. Finally, we show that the crystalline C₆₀ films are exclusively (111)-oriented and the fullerene islands are even heteroepitaxially aligned on the PEN buffer.

Two temperature regimes of triplet transfer in the dissociation of the correlated triplet pair after singlet fission

The ability to undergo spin-allowed exciton multiplication makes singlet fission materials promising for photovoltaic applications. Here, we examine the separation of correlated triplet pairs, 1(T…T), in polycrystalline pentacene films via temperature-dependent transient absorption spectroscopy. Single wavelength analysis reveals a profound delay in 1(T…T) dynamics. Moreover, the dynamics of 1(T…T) exhibit temperature dependence, whereas other features show no discernable temperature dependence. Previous literatures have suggested that correlated triplet separation is mediated by a thermally activated hopping process. Surprisingly, we found that the time constants governing triplet pair separation display two distinct temperature-dependent regimes of triplet transport. The high temperature regime follows a thermally activated hopping mechanism. The experimentally derived reorganization energy and electronic coupling is verified by density matrix renormalization group quantum chemical calculations. In addition, we evaluated the low temperature regime and show that the trend can be modelled by a Miller–Abrahams-type model that incorporates the effects of energetic disorder. We conclude that the correlated triplet pair separation is mediated by thermally activated hopping or a disorder driven Miller–Abrahams-type mechanism at high and low temperature, respectively. We observe that crossover between two regimes occurs ∼226 K. We find the time constant for triplet–triplet energy transfer to be 1.8 ps at ambient temperature and 21 ps at 77 K.

Toward Design of Novel Materials for Organic Electronics

Materials for organic electronics are presently used in prominent applications, such as displays in mobile devices, while being intensely researched for other purposes, such as organic photovoltaics, large-area devices, and thin-film transistors. Many of the challenges to improve and optimize these applications are material related and there is a nearly infinite chemical space that needs to be explored to identify the most suitable material candidates. Established experimental approaches struggle with the size and complexity of this chemical space. Herein, the development of simulation methods is addressed, with a particular emphasis on predictive multiscale protocols, to complement experimental research in the identification of novel materials and illustrate the potential of these methods with a few prominent recent applications. Finally, the potential of machine learning and methods based on artificial intelligence is discussed to further accelerate the search for new materials.

Diels-Alder adduct formation at solid interfaces between fullerenes and acenes

Understanding organic-organic interfaces is rather challenging due to their large complexity regarding morphology, molecular orientation at the interface, interdiffusion, and energetics. One additional important but often neglected aspect are chemical reactions occuring at such interfaces. For solid interfaces between pentacene and Buckminster-Fullerene (C₆₀) recently very efficient Diels-Alder (D-A) adduct formation has been reported. Considering the importance of pentacene/C₆₀ as prototypical donor-acceptor combination to study fundamental processes in organic photovoltaics, understanding this effect is essential. In this work, we provide detailed NEXAFS-based investigations with respect to the temperature-dependence and reaction zone depth of this effect. Moreover, we widely vary the interface morphology and observe that the D-A adduct formation is most efficient for bulk heterojunctions of pentacene and C₆₀. By also investigating further material combinations such as PEN/C₆₀-PCBM and interfaces between C₆₀ and functionalized acenes, we observe trends for the occurrence of the D-A adduct formation correlated with the different chemical properties of the involved compounds.

A molecular dynamics study on the interface morphology of vapor-deposited amorphous organic thin films

The interfaces between amorphous organic layers play an important role in the efficiency and lifetime of organic light emitting diodes (OLEDs). However, an atomistic understanding of the interface morphology is still poor. In this study, we theoretically investigate the interfacial structure of amorphous organic films using molecular dynamics simulations that mimic vapor-deposition processes. We find that molecularly sharp interfaces are formed by the vapor-deposition process as the interface thickness spans only a mono- or double-layer in terms of lie-down geometry. Interestingly, the interface is more diffusive into the upper layer due to asymmetric interdiffusion during the vapor-deposition process, which is well described by a simple random-walk model. Additionally, we investigate the change in the molecular orientation of interdiffused molecules, which is crucial for device performance.

Pentacene Crystal Growth on Silica and Layer-Dependent Step-Edge Barrier from Atomistic Simulations

Understanding and controlling the growth of organic crystals deposited from the vapor phase is important for fundamental materials science and necessary for applications in pharmaceutical and organic electronics industries. Here, this process is studied for the paradigmatic case of pentacene on silica by means of a specifically tailored computational approach inspired by the experimental vapor deposition process. This scheme is able to reproduce the early stages of the thin-film formation, characterized by a quasi layer-by-layer growth, thus showcasing its potential as a tool complementary to experimental techniques for investigating organic crystals. Crystalline islands of standing molecules are formed at a critical coverage, as a result of a collective reorientation of disordered aggregates of flat-lying molecules. The growth then proceeds by sequential attachment of molecules at the cluster and then terrace edges. Free-energy calculations allowed us to characterize the step-edge barrier for descending the terraces, a fundamental parameter for growth models for which only indirect experimental measurements are available. The barrier is found to be layer-dependent (approximately 1 kcal/mol for the first monolayer on silica, 2 kcal/mol for the second monolayer) and to extend over a distance comparable with the molecular length.

Computational Design of Thermally Activated Delayed Fluorescence Materials: The Challenges Ahead

Thermally activated delayed fluorescence (TADF) offers promise for all-organic light-emitting diodes with quantum efficiencies competing with those of transition-metal-based phosphorescent devices. While computational efforts have so far largely focused on gas-phase calculations of singlet and triplet excitation energies, the design of TADF materials requires multiple methodological developments targeting among others a quantitative description of electronic excitation energetics, fully accounting for environmental electrostatics and molecular conformational effects, the accurate assessment of the quantum mechanical interactions that trigger the elementary electronic processes involved in TADF, and a robust picture for the dynamics of these fundamental processes. In this Perspective, we describe some recent progress along those lines and highlight the main challenges ahead for modeling, which we hope will be useful to the whole TADF community.

Crystal alignment of caffeine deposited onto single crystal surfaces via hot-wall epitaxy

Crystal growth of caffeine on single crystalline surfaces yields needle or bird-like shaped crystals depending on surface chemistry and symmetry.

Defined crystal growth is highly demanded for technological applications but also fundamental research. Within this work, the crystal growth of the asymmetric molecule caffeine was studied on single crystalline surfaces of muscovite mica, sodium chloride and potassium chloride. While elongated needle-like crystals grow on muscovite mica and sodium chloride, smaller individual “bird-like” structures were observed on potassium chloride. Depending on the surface type and temperature, the disk-shaped caffeine molecules prefer either an edge-on or flat-on orientation with respect to the surface, but in each case, a defined crystallographic relation between the surface and caffeine crystallites was determined by using the X-ray pole figure technique. On muscovite mica and sodium chloride, needle-like crystallites with edge-on oriented molecules aligned mainly with the unit cell c-axis (which coincides with the long needle axis) along the [1–10]_mica, [100]_mica, [110]_mica and [110]_NaCl, [1–10]_NaCl directions, respectively. Crystals consisting of flat-on oriented molecules on KCl showed also defined alignments with respect to the substrate, but due to the altered molecule–substrate contact, the b-axis aligned along [110]_KCl and [1–10]_KCl. Growth at elevated temperatures enabled changes in the crystal growth whereby more defined structures formed on NaCl. On KCl, the bird-like structures remained very similar, while caffeine on the mica surface at elevated temperatures resulted in even additional texture forming with the caffeine molecules now also favoring a flat-on orientation with respect to the surface. The systematic variation of various system parameters demonstrates how sensitive the growth behavior of caffeine on this variety of substrates is.

Predicting the Conditions for Homeotropic Anchoring of Liquid Crystals at a Soft Surface. 4-n-Pentyl-4'-cyanobiphenyl on Alkylsilane Self-Assembled Monolayers

We have studied, using atomistic molecular dynamics simulations, the alignment of the nematic liquid-crystal 4-n-pentyl-4'-cyanobiphenyl (5CB) on self-assembled monolayers (SAMs) formed from octadecyl- and/or hexyltrichlorosilane (OTS and HTS) attached to glassy silica. We find a planar alignment on OTS at full coverage and an intermediate situation at partial OTS coverage because of the penetration of 5CB molecules into the monolayer, which also removes the tilt of the OTS SAM. Binary mixtures of HTS and OTS SAMs instead induce homeotropic (i.e., perpendicular) alignment. A comparison with the existing experimental literature is provided.

Theoretical Study of the Charge-Transfer State Separation within Marcus Theory: The C60-Anthracene Case Study

We study, within Marcus theory, the possibility of the charge-transfer (CT) state splitting at organic interfaces and a subsequent transport of the free charge carriers to the electrodes. As a case study we analyze model anthracene-C60 interfaces. Kinetic Monte Carlo (KMC) simulations on the cold CT state were performed at a range of applied electric fields, and with the fields applied at a range of angles to the interface to simulate the action of the electric field in a bulk heterojunction (BHJ) interface. The results show that the inclusion of polarization in our model increases CT state dissociation and charge collection. The effect of the electric field on CT state splitting and free charge carrier conduction is analyzed in detail with and without polarization. Also, depending on the relative orientation of the anthracene and C60 molecules at the interface, CT state splitting shows different behavior with respect to both applied field strength and applied field angle. The importance of the hot CT in helping the charge carrier dissociation is also analyzed in our scheme.

Molecular Dynamics Simulations and their Application to Thin-film Devices

The performance of devices based on organic semiconductors strongly depends on the molecular organisation in thin films. Due to the intrinsic complexity of these systems, a combination of theoretical modelling and experimental techniques is often the key to achieve a full understanding of their inner working. Here, we introduce the modelling of organic semiconductors by means of molecular dynamics simulations. We describe the basic theoretical framework of the technique and review the most popular class of force fields used to model organic materials, paying particular attention to the peculiarities of confined systems like nano-thick films. Representative studies of the organisation of organic functional materials in thin film phases are also reviewed.

Precise, Self-Limited Epitaxy of Ultrathin Organic Semiconductors and Heterojunctions Tailored by van der Waals Interactions

Precise assembly of semiconductor heterojunctions is the key to realize many optoelectronic devices. By exploiting the strong and tunable van der Waals (vdW) forces between graphene and organic small molecules, we demonstrate layer-by-layer epitaxy of ultrathin organic semiconductors and heterostructures with unprecedented precision with well-defined number of layers and self-limited characteristics. We further demonstrate organic p-n heterojunctions with molecularly flat interface, which exhibit excellent rectifying behavior and photovoltaic responses. The self-limited organic molecular beam epitaxy (SLOMBE) is generically applicable for many layered small-molecule semiconductors and may lead to advanced organic optoelectronic devices beyond bulk heterojunctions.

Surface Reconstructions in Organic Crystals: Simulations of the Effect of Temperature and Defectivity on Bulk and (001) Surfaces of 2,2':6',2″-Ternaphthalene

2,2':6',2″-Ternaphthalene (NNN) is a novel, blue-emitting material, suitable for preparation of organic light-emitting diodes. Its crystal structure has been solved recently, but its thermal behavior and surface properties have not yet been explored, partly due to the difficulty in obtaining high quality crystals. In the present study we use classical molecular dynamics to investigate the thermal behavior of bulk and (001) surfaces of NNN. Our bulk simulations indicate the occurrence of a phase transition at about 600 K. The transition is facilitated by the presence of a free (001) surface, since a reconstruction leading to a very similar structure occurs around 550 K in our surface models. This holds for both ideal and defective surface models, containing a small number of vacancies (one or two missing molecules in the outermost layer). In all cases, the process is triggered by thermal motion and involves the reorientation of the molecules with respect to the (001) plane. Both the bulk and surface phases share the monoclinic space group P21/a with a herringbone disposition of molecules. These findings and their implications for the use of NNN in organic electronics are discussed.

The role of charge transfer in the energy level alignment at the pentacene/C60 interface

Understanding the mechanism of energy level alignment at organic-organic interfaces is a crucial line of research to optimize applications in organic electronics. We address this problem for the C60-pentacene interface by performing local-orbital Density Functional Theory (DFT) calculations, including the effect of the charging energies on the energy gap of both organic materials. The results are analyzed within the induced density of interface states (IDIS) model. We find that the induced interface potential is in the range of 0.06-0.10 eV, in good agreement with the experimental evidence, and that such potential is mainly induced by the small, but non-negligible, charge transfer between the two compounds and the multipolar contribution associated with pentacene. We also suggest that an appropriate external intercompound potential could create an insulator-metal transition at the interface.

Diffusion-controlled growth of molecular heterostructures: fabrication of two-, one-, and zero-dimensional C(60) nanostructures on pentacene substrates

A variety of low dimensional C60 structures has been grown on supporting pentacene multilayers. By choice of substrate temperature during growth the effective diffusion length of evaporated fullerenes and their nucleation at terraces or step edges can be precisely controlled. AFM and SEM measurements show that this enables the fabrication of either 2D adlayers or solely 1D chains decorating substrate steps, while at elevated growth temperature continuous wetting of step edges is prohibited and instead the formation of separated C60 clusters pinned at the pentacene step edges occurs. Remarkably, all structures remain thermally stable at room temperature once they are formed. In addition the various fullerene structures have been overgrown by an additional pentacene capping layer. Utilizing the different probe depth of XRD and NEXAFS, we found that no contiguous pentacene film is formed on the 2D C60 structure, whereas an encapsulation of the 1D and 0D structures with uniformly upright oriented pentacene is achieved, hence allowing the fabrication of low dimensional buried organic heterostructures.

Singlet exciton fission in polycrystalline pentacene: from photophysics toward devices

Singlet exciton fission is the process in conjugated organic molecules bywhich a photogenerated singlet exciton couples to a nearby chromophore in the ground state, creating a pair of triplet excitons. Researchers first reported this phenomenon in the 1960s, an event that sparked further studies in the following decade. These investigations used fluorescence spectroscopy to establish that exciton fission occurred in single crystals of several acenes. However, research interest has been recently rekindled by the possibility that singlet fission could be used as a carrier multiplication technique to enhance the efficiency of photovoltaic cells. The most successful architecture to-date involves sensitizing a red-absorbing photoactive layer with a blue-absorbing material that undergoes fission, thereby generating additional photocurrent from higher-energy photons. The quest for improved solar cells has spurred a drive to better understand the fission process, which has received timely aid from modern techniques for time-resolved spectroscopy, quantum chemistry, and small-molecule device fabrication. However, the consensus interpretation of the initial studies using ultrafast transient absorption spectroscopy was that exciton fission was suppressed in polycrystalline thin films of pentacene, a material that would be otherwise expected to be an ideal model system, as well as a viable candidate for fission-sensitized photovoltaic devices. In this Account, we review the results of our recent transient absorption and device-based studies of polycrystalline pentacene. We address the controversy surrounding the assignment of spectroscopic features in transient absorption data, and illustrate how a consistent interpretation is possible. This work underpins our conclusion that singlet fission in pentacene is extraordinarily rapid (∼80 fs) and is thus the dominant decay channel for the photoexcited singlet exciton. Further, we discuss our demonstration that triplet excitons generated via singlet fission in pentacene can be dissociated at an interface with a suitable electron acceptor, such as fullerenes and infrared-absorbing inorganic semiconducting quantum dots. We highlight our recent reports of a pentacene/PbSe hybrid solar cell with a power conversion efficiency of 4.7% and of a pentacene/PbSe/amorphous silicon photovoltaic device. Although substantive challenges remain, both to better our understanding of the mechanism of singlet exciton fission and to optimize device performance, this realization of a solar cell where photocurrent is simultaneously contributed from a blue-absorbing fission-capable material and an infrared-absorbing conventional cell is an important step towards a dual-bandgap, single-junction, fission-enhanced photovoltaic device, which could one day surpass the Shockley-Queisser limit.

Exploring the energy landscape of the charge transport levels in organic semiconductors at the molecular scale

The extraordinary semiconducting properties of conjugated organic materials continue to attract attention across disciplines including materials science, engineering, chemistry, and physics, particularly with application to organic electronics. Such materials are used as active components in light-emitting diodes, field-effect transistors, or photovoltaic cells, as a substitute for (mostly Si-based) inorganic semiconducting materials. Many strategies developed for inorganic semiconductor device building (doping, p-n junctions, etc.) have been attempted, often successfully, with organics, even though the key electronic and photophysical properties of organic thin films are fundamentally different from those of their bulk inorganic counterparts. In particular, organic materials consist of individual units (molecules or conjugated segments) that are coupled by weak intermolecular forces. The flexibility of organic synthesis has allowed the development of more efficient opto-electronic devices including impressive improvements in quantum yields for charge generation in organic solar cells and in light emission in electroluminescent displays. Nonetheless, a number of fundamental questions regarding the working principles of these devices remain that preclude their full optimization. For example, the role of intermolecular interactions in driving the geometric and electronic structures of solid-state conjugated materials, though ubiquitous in organic electronic devices, has long been overlooked, especially when it comes to these interfaces with other (in)organic materials or metals. Because they are soft and in most cases disordered, conjugated organic materials support localized electrons or holes associated with local geometric distortions, also known as polarons, as primary charge carriers. The spatial localization of excess charges in organics together with low dielectric constant (ε) entails very large electrostatic effects. It is therefore not obvious how these strongly interacting electron-hole pairs can potentially escape from their Coulomb well, a process that is at the heart of photoconversion or molecular doping. Yet they do, with near-quantitative yield in some cases. Limited screening by the low dielectric medium in organic materials leads to subtle static and dynamic electronic polarization effects that strongly impact the energy landscape for charges, which offers a rationale for this apparent inconsistency. In this Account, we use different theoretical approaches to predict the energy landscape of charge carriers at the molecular level and review a few case studies highlighting the role of electrostatic interactions in conjugated organic molecules. We describe the pros and cons of different theoretical approaches that provide access to the energy landscape defining the motion of charge carriers. We illustrate the applications of these approaches through selected examples involving OFETs, OLEDs, and solar cells. The three selected examples collectively show that energetic disorder governs device performances and highlights the relevance of theoretical tools to probe energy landscapes in molecular assemblies.

Molecular thin films on solid surfaces: mechanisms of melting

We use molecular dynamics simulations to study the melting of pentane and hexane monolayers adsorbed on the basal plane of graphite. For both of these systems, the temperature-dependent structures and the melting temperatures agree well with experiment. A detailed analysis reveals that a mechanism involving the promotion of molecules to the second layer underlies melting in these systems. In the second-layer promotion mechanism, a small fraction of molecules transition into the second layer around the melting temperature, leaving vacant space in the first layer to facilitate disordering. The second-layer promotion mechanism arises because of the weaker molecule-surface interaction in our study than that in previous studies. The weaker molecule-surface interaction is consistent with experimental temperature-programmed desorption studies.

Enhanced charge transfer by phenyl groups at a rubrene/C60 interface

Exciton dynamics at an interface between an electron donor, rubrene, and a C(60) acceptor is studied by nonadiabatic quantum molecular dynamics simulation. Simulation results reveal an essential role of the phenyl groups in rubrene in increasing the charge-transfer rate by an order-of-magnitude. The atomistic mechanism of the enhanced charge transfer is found to be the amplification of aromatic breathing modes by the phenyl groups, which causes large fluctuations of electronic excitation energies. These findings provide insight into molecular structure design for efficient solar cells, while explaining recent experimental observations.

Anisotropic behavior of organic molecules on prepatterned surfaces

The nucleation of organic molecules on surfaces, prepatterned with stripes, is investigated with emphasis on anisotropy effects. Representing the molecules as ellipsoids, the related particle-particle interaction is modeled by means of a generalized Gay-Berne potential for similar biaxial particles. The orientation behavior of these ellipsoidal molecules induced by the stripe pattern is studied for the first monolayer by performing kinetic Monte Carlo simulations. It is shown how the properties of the particle alignment depend on energy scales, temperature, and flux. Based on the fact the particles strictly arrange in rows, it is furthermore instructive to analyze the orientation behavior within the different rows. Finally, the transfer of orientation from a preset row of molecules with fixed orientation to other nucleating particles is examined.