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.

The impact of film deposition and annealing on the nanostructure and dielectric constant of organic semiconductor thin films.

The strategy of using a bulk-heterojunction light-absorbing layer has led to the most efficient organic solar cells. However, optimising the blend morphology to maximise light absorption, charge generation and extraction can be challenging. Homojunction devices containing a single component have the potential to overcome the challenges associated with bulk heterojunction films. A strategy towards this goal is to increase the dielectric constant of the organic semiconductor to ≈10, which in principle would lead to free charge carrier generation upon photoexcitation. However, the factors that affect the thin film dielectric constants are still not well understood. In this work we report an organic semiconductor material that can be solution processed or vacuum evaporated to form good quality thin films to explore the effect of chromophore structure and film morphology on the dielectric constant and other optoelectronic properties. 2,2'-[(4,4,4',4'-Tetrakis{2-[2-methoxyethoxy]ethyl}-4H,4'H-{2,2'-bi[cyclo-penta[2,1-b:3,4-b']dithiophene]}-6,6'-diyl)bis(methaneylylidene)]dimalononitrile [D(CPDT-DCV)] was designed to have high electron-affinity end groups and low ionisation-potential central moieties. It can be processed from solution or be thermally evaporated, with the film morphology changing from face-on to a herringbone arrangement upon solvent or thermal annealing. The glycol solubilising groups led to the static dielectric constant (taken from capacitance measurements) of the films to be between 6 and 7 (independent of processing conditions), while the optical frequency dielectric constant depended on the processing conditions. The less ordered solution processed film was found to have the lowest optical frequency dielectric constant of 3.6 at 2.0 × 1014 Hz, which did not change upon annealing. In contrast, the more ordered evaporated film had an optical frequency dielectric constant 20% higher at 4.2 and thermal annealing further increased it to 4.5, which is amongst the highest reported for an organic semiconductor at that frequency. Finally, the more ordered evaporated films had more balanced charge transport, which did not change upon annealing.

Understanding the performance differences between solution and vacuum deposited OLEDs: A computational approach.

Solution-processing of organic light-emitting diode films has potential advantages in terms of cost and scalability over vacuum-deposition for large area applications. However, solution processed small molecule films can have lower overall device performance. Here, novel molecular dynamics techniques are developed to enable faster simulation of solvent evaporation that occurs during solution processing and give films of thicknesses relevant to real devices. All-atom molecular dynamics simulations are then used in combination with kinetic Monte Carlo transport modeling to examine how differences in morphology stemming from solution or vacuum film deposition affect charge transport and exciton dynamics in films consisting of light-emitting bis(2-phenylpyridine)(acetylacetonate)iridium(III) [Ir(ppy)2(acac)] guest molecules in a 4,4'-bis(N-carbazolyl)biphenyl host. While the structures of the films deposited from vacuum and solution were found to differ, critically, only minor variations in the transport properties were predicted by the simulations even if trapped solvent was present.

Unraveling exciton processes in Ir(ppy)3:CBP OLED films upon photoexcitation.

Emissive layers in phosphorescent organic light-emitting diodes commonly make use of guest-host blends such as Ir(ppy)3:CBP to achieve high external quantum efficiencies. However, while the Ir(ppy)3:CBP blend has been studied experimentally, crucial questions remain regarding how exciton diffusion is dependent on the distribution of the guest in the host, which can currently only be addressed at the atomic level via computational modeling. In this work, kinetic Monte Carlo simulations are utilized to gain insight into exciton diffusion in Ir(ppy)3:CBP blend films. The effects of both guest concentration and exciton density on various system properties are analyzed, including the probability of singlet excitons being converted to triplets, and the probability of those triplets decaying radiatively. Significantly, these simulations suggest that triplet diffusion occurs almost exclusively via guest-guest Dexter transfer and that concentration quenching of triplets induced by guest-guest intermolecular dipole-dipole interactions has a negligible effect at high exciton densities due to the prevalence of triplet-triplet annihilation. Furthermore, results for vacuum deposited morphologies derived from molecular dynamics simulations are compared to the results obtained using a simple cubic lattice approximation with randomly distributed guest molecules. We show that while differences in host-based processes such as singlet diffusion are observed, overall, the results on the fate of the excitons are in good agreement for the two morphology types, particularly for guest-based processes at low guest concentrations where guest clustering is limited.

Pyrrolo[3,2-b]pyrrole-1,4-dione (IsoDPP) End Capped with Napthalimide or Phthalimide: Novel Small Molecular Acceptors for Organic Solar Cells.

We introduce two novel solution-processable electron acceptors based on an isomeric core of the much explored diketopyrrolopyrrole (DPP) moiety, namely pyrrolo[3,2-b]pyrrole-1,4-dione (IsoDPP). The newly designed and synthesized compounds, 6,6'-[(1,4-bis{4-decylphenyl}-2,5-dioxo-1,2,4,5-tetrahydropyrrolo[3,2-b]pyrrole-3,6-diyl)bis(thiophene-5,2-diyl)]bis[2-(2-butyloctyl)-1H-benzo[de]isoquinoline-1,3(2H)-dione] (NAI-IsoDPP-NAI) and 5,5'-[(1,4-bis{4-decylphenyl}-2,5-dioxo-1,2,4,5-tetrahydropyrrolo[3,2-b]pyrrole-3,6-diyl)bis(thiophene-5,2-diyl)]bis[2-(2-butyloctyl)isoindoline-1,3-dione] (PI-IsoDPP-PI) have been synthesized via Suzuki couplings using IsoDPP as a central building block and napthalimide or phthalimide as end-capping groups. The materials both exhibit good solubility in a wide range of organic solvents including chloroform (CF), dichloromethane (DCM), and tetrahydrofuran (THF), and have a high thermal stability. The new materials absorb in the wavelength range of 300-600 nm and both compounds have similar electron affinities, with the electron affinities that are compatible with their use as acceptors in donor-acceptor bulk heterojunction (BHJ) organic solar cells. BHJ devices comprising the NAI-IsoDPP-NAI acceptor with poly(3-n-hexylthiophene) (P3HT) as the donor were found to have a better performance than the PI-IsoDPP-PI containing cells, with the best device having a VOC of 0.92 V, a JSC of 1.7 mAcm-2, a FF of 63%, and a PCE of 0.97%.

Perdeuterated Conjugated Polymers for Ultralow-Frequency Magnetic Resonance of OLEDs.

The formation of excitons in OLEDs is spin dependent and can be controlled by electron-paramagnetic resonance, affecting device resistance and electroluminescence yield. We explore electrically detected magnetic resonance in the regime of very low magnetic fields (<1 mT). A pronounced feature emerges at zero field in addition to the conventional spin- 1 / 2 Zeeman resonance for which the Larmor frequency matches that of the incident radiation. By comparing a conventional π-conjugated polymer as the active material to a perdeuterated analogue, we demonstrate the interplay between the zero-field feature and local hyperfine fields. The zero-field peak results from a quasistatic magnetic-field effect of the RF radiation for periods comparable to the carrier-pair lifetime. Zeeman resonances are resolved down to 3.2 MHz, approximately twice the Larmor frequency of an electron in Earth's field. However, since reducing hyperfine fields sharpens the Zeeman peak at the cost of an increased zero-field peak, we suggest that this result may constitute a fundamental low-field limit of magnetic resonance in carrier-pair-based systems. OLEDs offer an alternative solid-state platform to investigate the radical-pair mechanism of magnetic-field effects in photochemical reactions, allowing models of biological magnetoreception to be tested by measuring spin decoherence directly in the time domain by pulsed experiments.

Calculating transition dipole moments of phosphorescent emitters for efficient organic light-emitting diodes.

The out-coupling of light from an organic light-emitting diode, and thus its efficiency, strongly depends on the orientation of the transition dipole moment (TDM) of the emitting molecules with respect to the substrate surface. Despite the importance of this quantity, theoretical investigations of the direction of the TDM of phosphorescent emitters based on iridium(iii) complexes remain limited. One challenge is to find an appropriate level of theory able to accurately predict the direction of the TDM. Here, we report relativistic time-dependent density functional theory (TDDFT) calculations of the TDM, emission energies and lifetimes for both the ground-state (S0) and triplet (T1) excited-state geometries of fac-tris(2-phenylpyridyl)iridium(iii) (Ir(ppy)3), using the two-component zero-order regular approximation (ZORA) or including spin-orbit coupling (SOC) perturbatively using the simpler one-component (scalar) formulation. We show that the one- and two-component approaches give similar emission energies and overall radiative lifetimes for each individual geometry. Use of the S0 geometry leads to two of the excited triplet substates being degenerate, with the degeneracy lifted for the T1 geometry, with the latter matching experiment. Two-component calculations using the T1 geometry give results for the direction of the TDM more consistent with experiment than calculations using the S0 geometry. Finally, we show that adding a dielectric medium does not affect the direction of TDM significantly, but leads to better agreement with the experimentally measured radiative lifetime.

Understanding charge transport in Ir(ppy)3:CBP OLED films.

Ir(ppy)3:CBP blends have been widely studied as the emissive layer in organic light emitting diodes (OLEDs), yet crucial questions about charge transport within the layer remain unaddressed. Recent molecular dynamics simulations show that the Ir(ppy)3 molecules are not isolated from each other, but at concentrations of as low as 5 wt. % can be part of connected pathways. Such connectivity raises the question of how the iridium(iii) complexes contribute to long-range charge transport in the blend. We implement a kinetic Monte Carlo transport model to probe the guest concentration dependence of charge mobility and show that distinct minima appear at approximately 10 wt. % Ir(ppy)3 due to an increased number of trap states that can include interconnected complexes within the blend film. The depth of the minima is shown to be dependent on the electric field and to vary between electrons and holes due to their different trapping depths arising from the different ionization potentials and electron affinities of the guest and host molecules. Typical guest-host OLEDs use a guest concentration below 10 wt. % to avoid triplet-triplet annihilation, so these results suggest that optimal device performance is achieved when there is significant charge trapping on the iridium(iii) complex guest molecules and minimum interactions of the emissive chromophores that can lead to triplet-triplet annihilation.

Charge Generation Pathways in Organic Solar Cells: Assessing the Contribution from the Electron Acceptor.

Photocurrent generation in organic bulk heterojunction (BHJ) solar cells is most commonly understood as a process which predominantly involves photoexcitation of the lower ionization potential species (donor) followed by electron transfer to the higher electron affinity material (acceptor) [i.e., photoinduced electron transfer (PET), which we term Channel I]. A mirror process also occurs in which photocurrent is generated through photoexcitation of the acceptor followed by hole transfer to the nonexcited donor or photoinduced hole transfer (PHT), which we term Channel II. The role of Channel II photocurrent generation has often been neglected due to overlap of the individual absorption spectra of the donor and acceptor materials that are commonly used. More recently Channel II charge generation has been explored for several reasons. First, many of the new high-efficiency polymeric donors are used as the minority component in bulk heterojunction blends, and therefore, the acceptor absorption is a significant fraction of the total; second, nonfullerene acceptors have been prepared, which through careful design, allow for spectral separation from the donor material, facilitating fundamental studies on charge generation. In this article, we review the methodologies for investigating the two charge generation channels. We also discuss the factors that affect charge generation via Channel I and II pathways, including energy levels of the materials involved, exciton diffusion, and other considerations. Finally, we take a comprehensive look at the nonfullerene acceptor literature and discuss what information about Channel I and Channel II can be obtained from the experiments conducted and what other experiments could be undertaken to provide further information about the operational efficiencies of Channels I and II.

The Molecular Origin of Anisotropic Emission in an Organic Light-Emitting Diode.

Atomistic nonequilibrium molecular dynamics simulations have been used to model the induction of molecular orientation anisotropy within the emission layer of an organic light-emitting diode (OLED) formed by vapor deposition. Two emitter species were compared: racemic fac-tris(2-phenylpyridine)iridium(III) (Ir(ppy)3) and trans-bis(2-phenylpyridine)(acetylacetonate)iridium(III) (Ir(ppy)2(acac)). The simulations show that the molecular symmetry axes of both emitters preferentially align perpendicular to the surface during deposition. The molecular arrangement formed on deposition combined with consideration of the transition dipole moments provides insight into experimental reports that Ir(ppy)3 shows isotropic emission, while Ir(ppy)2(acac) displays improved efficiency due to an apparent preferential alignment of the transition dipole vectors parallel to the substrate. The simulations indicate that this difference is not due to differences in the extent of emitter alignment, but rather differences in the direction of the transition dipoles within the two complexes.

Organohalide Perovskites for Solar Energy Conversion.

Lead-based organohalide perovskites have recently emerged as arguably the most promising of all next generation thin film solar cell technologies. Power conversion efficiencies have reached 20% in less than 5 years, and their application to other optoelectronic device platforms such as photodetectors and light emitting diodes is being increasingly reported. Organohalide perovskites can be solution processed or evaporated at low temperatures to form simple thin film photojunctions, thus delivering the potential for the holy grail of high efficiency, low embedded energy, and low cost photovoltaics. The initial device-driven "perovskite fever" has more recently given way to efforts to better understand how these materials work in solar cells, and deeper elucidation of their structure-property relationships. In this Account, we focus on this element of organohalide perovskite chemistry and physics in particular examining critical electro-optical, morphological, and architectural phenomena. We first examine basic crystal and chemical structure, and how this impacts important solar-cell related properties such as the optical gap. We then turn to deeper electronic phenomena such as carrier mobilities, trap densities, and recombination dynamics, as well as examining ionic and dielectric properties and how these two types of physics impact each other. The issue of whether organohalide perovskites are predominantly nonexcitonic at room temperature is currently a matter of some debate, and we summarize the evidence for what appears to be the emerging field consensus: an exciton binding energy of order 10 meV. Having discussed the important basic chemistry and physics we turn to more device-related considerations including processing, morphology, architecture, thin film electro-optics and interfacial energetics. These phenomena directly impact solar cell performance parameters such as open circuit voltage, short circuit current density, internal and external quantum efficiency, fill factor, and ultimately the all-important power conversion efficiency. Finally, we address the key challenges pertinent to actually delivering a new and viable solar cell technology. These include long-term cell stability, scaling to the module level, and the toxicity associated with lead. Organohalide perovskites not only offer exciting possibilities for next generation optoelectronics and photovoltaics, but are an intriguing class of material crossing the boundaries of molecular solids and banded inorganic semiconductors. This is a potential area of rich new chemistry, materials science, and physics.

The structural impact of water sorption on device-quality melanin thin films.

The melanins are a class of pigmentary bio-macromolecules ubiquitous in the biosphere. They possess an intriguing set of physico-chemical properties and have been shown to exhibit hybrid protonic-electronic electrical conductivity, a feature derived from a process termed chemical self-doping driven by the sorption of water. Although the mechanism underlying the electrical conduction has been established, how the sorbed water interacts with the melanin structure at the physical level has not. Herein we use neutron reflectometry to study changes in the structure of synthetic melanin thin films as a function of H2O and D2O vapour pressure. Water is found to be taken up evenly throughout the films, and by employing the contrast effect, the existence of labile protons through reversible deuterium exchange is demonstrated. Finally, we determine a sorption isotherm to enable quantification of the melanin-water interactions.

Effect of n-propyl substituents on the emission properties of blue phosphorescent iridium(iii) complexes.

Ligand substitution is often used for tuning the emission color of phosphorescent iridium(iii) complexes that are used in organic light-emitting diodes. However, in addition to tuning the emission color, the substituents can also affect the radiative and non-radiative decay rates of the excited state and hence the photoluminescence quantum yield. Understanding the substituent effect is therefore important for the design of new iridium(iii) complexes with specific emission properties. Using (time dependent) density functional methods, we investigate the substituent effect of n-propyl groups on the structure, emission color, and emission efficiency of fac-tris(1-methyl-5-phenyl-[1,2,4]triazolyl)iridium(iii) based phosphorescent complexes by comparing the calculated results for structural models with and without the n-propyl substituents. We find that attachment of the n-propyl groups increases the length of three Ir-N bonds, and although the emission color does not change significantly, the radiative and non-radiative rates do, leading to a prediction of enhanced blue phosphorescence emission efficiency. Furthermore, the calculations show that the attachment of the n-propyl groups leads to a larger activation energy to degradation and the formation of dark states.

Elucidating the Spatial Arrangement of Emitter Molecules in Organic Light-Emitting Diode Films.

The effect of varying the emitter concentration on the structural properties of an archetypal phosphorescent blend consisting of 4,4'-bis(N-carbazolyl)biphenyl and tris(2-phenylpyridyl)iridium(III) has been investigated using non-equilibrium molecular dynamics (MD) simulations that mimic the process of vacuum deposition. By comparison with reflectometry measurements, we show that the simulations provide an accurate model of the average density of such films. The emitter molecules were found not to be evenly distributed throughout film, but rather they can form networks that provide charge and/or energy migration pathways, even at emitter concentrations as low as ≈5 weight percent. At slightly higher concentrations, percolated networks form that span the entire system. While such networks would give improved charge transport, they could also lead to more non-radiative pathways for the emissive state and a resultant loss of efficiency.

Detection of Explosive Vapors: The Roles of Exciton and Molecular Diffusion in Real-Time Sensing.

Time-resolved quartz crystal microbalance with in situ fluorescence measurements are used to monitor the sorption of the nitroaromatic (explosive) vapor, 2,4-dinitrotoluene (DNT) into a porous pentiptycene-containing poly(phenyleneethynylene) sensing film. Correlation of the nitroaromatic mass uptake with fluorescence quenching shows that the analyte diffusion follows the Case-II transport model, a film-swelling-limited process, in which a sharp diffusional front propagates at a constant velocity through the film. At a low vapor pressure of DNT of ≈16 ppb, the analyte concentration in the front is sufficiently high to give an average fluorophore-analyte separation of ≈1.5 nm. Hence, a long exciton diffusion length is not required for real-time sensing in the solid state. Rather the diffusion behavior of the analyte and the strength of the binding interaction between the analyte and the polymer play first-order roles in the fluorescence quenching process.

Bond Fission and Non-Radiative Decay in Iridium(III) Complexes.

We investigate the role of metal-ligand bond fission in the nonradiative decay of excited states in iridium(III) complexes with applications in blue organic light-emitting diodes (OLEDs). We report density functional theory (DFT) calculations of the potential energy surfaces upon lengthening an iridium-nitrogen (Ir-N) bond. In all cases we find that for bond lengths comparable to those of the ground state the lowest energy state is a triplet with significant metal-to-ligand change transfer character ((3)MLCT). But, as the Ir-N bond is lengthened there is a sudden transition to a regime where the lowest excited state is a triplet with significant metal centered character ((3)MC). Time-dependent DFT relativistic calculations including spin-orbit coupling perturbatively show that the radiative decay rate from the (3)MC state is orders of magnitude slower than that from the (3)MLCT state. The calculated barrier height between the (3)MLCT and (3)MC regimes is clearly correlated with previously measured nonradiative decay rates, suggesting that thermal population of the (3)MC state is the dominant nonradiative decay process at ambient temperature. In particular, fluorination both drives the emission of these complexes to a deeper blue color and lowers the (3)MLCT-(3)MC barrier. If the Ir-N bond is shortened in the (3)MC state another N atom is pushed away from the Ir, resulting in the breaking of this bond, suggesting that once the Ir-N bond breaks the damage to the complex is permanent-this will have important implications for the lifetimes of devices using this type of complex as the active material. The consequences of these results for the design of more efficient blue phosphors for OLED applications are discussed.

Photophysics of detection of explosive vapours via luminescence quenching of thin films: impact of inter-molecular interactions.

Fluorescence-based detection of explosive analytes requires an understanding of the nature of the excited state responsible for the luminescence response of a sensing material. Many measurements are carried out to elucidate the fundamental photophysical properties of an emissive material in solution. However, simple transfer of the understanding gained from the solution measurements to the solid-state can lead to errors. This is in part due to the absence of inter-molecular interactions of the chromophores in solution, which are present in the solid-state. To understand the role of inter-molecular interactions on the detection of explosive analytes we have chosen dendrimers from two different families, D1 and D2, which allow facile control of the inter-molecular interactions through the choice of dendrons and emissive chromophores. Using ultrafast transient absorption spectroscopy we find that the solution photoinduced absorption (PA) for both materials can be explained in terms of the generation of singlet excitons, which decay to the ground state, or intersystem cross (ISC) to form a triplet exciton. In neat films however, we observe different photophysical behaviours; first, ISC to the triplet state does not occur, and second, depending on the chromophore, charge transfer and charge separated states are formed. Furthermore, we find that when either dendrimer is interfaced with analyte vapour, the singlet state is strongly quenched, generating a charge transfer state that undergoes geminate recombination.

Tuning the optoelectronic properties of nonfullerene electron acceptors.

Broad spectral coverage over the solar spectrum is necessary for photovoltaic technologies and is a focus for organic solar cells. We report a series of small-molecule, nonfullerene electron acceptors containing the [(benzo[c][1,2,5]thiadiazol-4-yl)methylene]malononitrile unit as a high electron affinity component. The optoelectronic properties of these molecules were fine-tuned with the objective of attaining strong absorption at longer wavelengths by changing the low-ionization-potential moiety. The electron-accepting function of these materials was investigated with poly(3-n-hexylthiophene) (P3HT) as a standard electron donor. Significant photocurrent generation in the near infrared region, with an external quantum yield reaching as high as 22 % at 700 nm and an onset >800 nm was achieved. The results support efficient hole transfer to P3HT taking place after light absorption by the acceptor molecules. A Channel II-dominated power conversion efficiency of up to 1.5 % was, thus, achieved.

Interplay of Zero-Field Splitting and Excited State Geometry Relaxation in fac-Ir(ppy)3.

The lowest energy triplet state, T1, of organometallic complexes based on iridium(III) is of fundamental interest, as the behavior of molecules in this state determines the suitability of the complex for use in many applications, e.g., organic light-emitting diodes. Previous characterization of T1 in fac-Ir(ppy)3 suggests that the trigonal symmetry of the complex is weakly broken in the excited state. Here we report relativistic time dependent density functional calculations of the zero-field splitting (ZFS) of fac-Ir(ppy)3 in the ground state (S0) and lowest energy triplet (T1) geometries and at intermediate geometries. We show that the energy scale of the geometry relaxation in the T1 state is large compared to the ZFS. Thus, the natural analysis of the ZFS and the radiative decay rates, based on the assumption that the structural distortion is a small perturbation, fails dramatically. In contrast, our calculations of these quantities are in good agreement with experiment.

Planar silver nanowire, carbon nanotube and PEDOT:PSS nanocomposite transparent electrodes.

Highly conductive, transparent and flexible planar electrodes were fabricated using interwoven silver nanowires and single-walled carbon nanotubes (AgNW:SWCNT) in a PEDOT:PSS matrix via an epoxy transfer method from a silicon template. The planar electrodes achieved a sheet resistance of 6.6 ± 0.0 Ω/□ and an average transmission of 86% between 400 and 800 nm. A high figure of merit of 367 Ω-1 is reported for the electrodes, which is much higher than that measured for indium tin oxide and reported for other AgNW composites. The AgNW:SWCNT:PEDOT:PSS electrode was used to fabricate low temperature (annealing free) devices demonstrating their potential to function with a range of organic semiconducting polymer:fullerene bulk heterojunction blend systems.