Long-Range Embedding of Molecular Ions and Excitations in a Polarizable Molecular Environment

Quantum chemical package Jaguar: A survey of recent developments and unique features

This paper is dedicated to the quantum chemical package Jaguar, which is commercial software developed and distributed by Schrödinger, Inc. We discuss Jaguar's scientific features that are relevant to chemical research as well as describe those aspects of the program that are pertinent to the user interface, the organization of the computer code, and its maintenance and testing. Among the scientific topics that feature prominently in this paper are the quantum chemical methods grounded in the pseudospectral approach. A number of multistep workflows dependent on Jaguar are covered: prediction of protonation equilibria in aqueous solutions (particularly calculations of tautomeric stability and pKa), reactivity predictions based on automated transition state search, assembly of Boltzmann-averaged spectra such as vibrational and electronic circular dichroism, as well as nuclear magnetic resonance. Discussed also are quantum chemical calculations that are oriented toward materials science applications, in particular, prediction of properties of optoelectronic materials and organic semiconductors, and molecular catalyst design. The topic of treatment of conformations inevitably comes up in real world research projects and is considered as part of all the workflows mentioned above. In addition, we examine the role of machine learning methods in quantum chemical calculations performed by Jaguar, from auxiliary functions that return the approximate calculation runtime in a user interface, to prediction of actual molecular properties. The current work is second in a series of reviews of Jaguar, the first having been published more than ten years ago. Thus, this paper serves as a rare milestone on the path that is being traversed by Jaguar's development in more than thirty years of its existence.

Efficient and Accurate Estimation of Electronic Polarization Energies for Organic Semiconductors: An Embedding Charge Quantum Mechanics/Continuum Dielectric Model

Electronic polarization plays a pivotal role in determining the molecular energy levels of organic semiconductors (OSCs) in the condensed phase. However, accurate estimation of the electronic polarization energy is a challenging task due to the intricate imbalance between the precision and efficiency. In this work, we have developed an embedding charge quantum mechanics/continuum dielectric (EC-QM/CD) model, which enables quantitative evaluation of the ionization potential (IP), electron affinity (EA), and polarization energy in both crystalline and amorphous solids for OSCs. The benchmark calculations on both p-type OSCs of oligoacenes and n-type OSCs of A-D-A small-molecule acceptors show that the values of IP, EA, and polarization energy obtained by EC-QM/CD are in good accordance with the experimental measurements or the results by high-precision methods, while the computational costs are substantially reduced. Given its balance between the accuracy and efficiency, the EC-QM/CD model exhibits considerable potential to broaden the applications in the field of OSCs, for instance, high-throughput screening by using solid-state energy levels or polarization energies as critical descriptors.

New Avenues for Organic Solar Cells Using Intrinsically Charge-Generating Materials

The molecular electron acceptor material Y6 has been a key part of the most recent surge in organic solar cell sunlight-to-electricity power conversion efficiency, which is now approaching 20%. Numerous studies have sought to understand the fundamental photophysical reasons for the exceptional performance of Y6 and its growing family of structural derivatives. Though significant uncertainty about several details remains, many have concluded that initially photogenerated excited states rapidly convert into electron-hole charge pairs in the neat material. These charge pairs are characterized by location of the electron and hole on different Y6 molecules, in contrast to the Frenkel excitons that dominate the behavior of most organic semiconductor materials. Here, we summarize the current state of knowledge regarding Y6 photophysics and the key observations that have led to it. We then link this understanding to other advances, such as the role of quadrupolar fields in donor-acceptor blends, and the importance of molecular interactions and organization in providing the structural basis for Y6's properties. Finally, we turn our attention to ways of making use of the new photophysics of Y6, and suggest molecular doping, crystal structure tuning, and electric field engineering as promising avenues for future exploration.

Semitransparent Organic Photovoltaics Utilizing Intrinsic Charge Generation in Non-Fullerene Acceptors

In organic semiconductors, a donor/acceptor heterojunction is typically required for efficient dissociation of excitons. Using transient absorption spectroscopy to study the dynamics of excited states in non-fullerene acceptors (NFAs), it is shown that NFAs can generate charges without a donor/acceptor interface. This is due to the fact that dielectric solvation provides a driving force sufficient to dissociate the excited state and form the charge-transfer (CT) state. The CT state is further dissociated into free charges at interfaces between polycrystalline regions in neat NFAs. For IEICO-4F, incorporating just 9 wt% donor polymer PTB7-Th in neat films greatly boosts charge generation, enhancing efficient exciton separation into free charges. This property is utilized to fabricate donor-dilute organic photovoltaics (OPV) delivering a power conversion efficiency of 8.3% in the case of opaque devices with a metal top-electrode and an active layer average visible transmittance (AVT) of 75%. It is shown that the intrinsic charge generation in low-bandgap NFAs contributes to the overall photocurrent generation. IEICO-4F-based OPVs with limited PTB7-Th content have high thermal resilience demonstrating little drop in performance over 700 h. PTB7-Th:IEICO-4F semitransparent OPVs are leveraged to fabricate an 8-series connected semitransparent module, demonstrating light-utilization efficiency of 2.2% alongside an AVT of 63%.

Reduced bimolecular charge recombination in efficient organic solar cells comprising non-fullerene acceptors

Bimolecular charge recombination is one of the most important loss processes in organic solar cells. However, the bimolecular recombination rate in solar cells based on novel non-fullerene acceptors is mostly unclear. Moreover, the origin of the reduced-Langevin recombination rate in bulk heterojunction solar cells in general is still poorly understood. Here, we investigate the bimolecular recombination rate and charge transport in a series of high-performance organic solar cells based on non-fullerene acceptors. From steady-state dark injection measurements and drift-diffusion simulations of the current-voltage characteristics under illumination, Langevin reduction factors of up to over two orders of magnitude are observed. The reduced recombination is essential for the high fill factors of these solar cells. The Langevin reduction factors are observed to correlate with the quadrupole moment of the acceptors, which is responsible for band bending at the donor-acceptor interface, forming a barrier for charge recombination. Overall these results therefore show that suppressed bimolecular recombination is essential for the performance of organic solar cells and provide design rules for novel materials.

Theoretical Study of Charge Mobility in Crystal Porphine and a Computer Design of a Porphine-Based Semiconductive Discotic Liquid Mesophase

Organic semiconductors are the focus of numerous studies; they are used in electronic devices. Modern research involves the production of neuromorphic organic materials, including those based on liquid crystal materials. The purpose of this work involves the theoretical modeling of molecules (the "core with branches" type) to construct a discotic mesophase capable of performing the functions of a neuromorphic material. For this purpose, the conductivity of crystal porphine, which can act as the nucleus of a molecule of the "core with branches" type, was investigated. The Marcus theory charge mobility values for the hole and electron were 0.148 and 0.088 cm2/V·s, respectively (the MOO method for calculating transfer integrals), and 0.561 and 0.160 cm2/V·s (DIPRO method). Based on TD-HF (HF-3c level of theory) calculations, possible structures of molecules for the formation of a discotic mesophase are proposed.

Virtual Screening for Organic Solar Cells and Light Emitting Diodes

The field of organic semiconductors is multifaceted and the potentially suitable molecular compounds are very diverse. Representative examples include discotic liquid crystals, dye-sensitized solar cells, conjugated polymers, and graphene-based low-dimensional materials. This huge variety not only represents enormous challenges for synthesis but also for theory, which aims at a comprehensive understanding and structuring of the plethora of possible compounds. Eventually computational methods should point to new, better materials, which have not yet been synthesized. In this perspective, it is shown that the answer to this question rests upon the delicate balance between computational efficiency and accuracy of the methods used in the virtual screening. To illustrate the fundamentals of virtual screening, chemical design of non-fullerene acceptors, thermally activated delayed fluorescence emitters, and nanographenes are discussed.

Free charge photogeneration in a single component high photovoltaic efficiency organic semiconductor

Organic photovoltaics (OPVs) promise cheap and flexible solar energy. Whereas light generates free charges in silicon photovoltaics, excitons are normally formed in organic semiconductors due to their low dielectric constants, and require molecular heterojunctions to split into charges. Recent record efficiency OPVs utilise the small molecule, Y6, and its analogues, which – unlike previous organic semiconductors – have low band-gaps and high dielectric constants. We show that, in Y6 films, these factors lead to intrinsic free charge generation without a heterojunction. Intensity-dependent spectroscopy reveals that 60–90% of excitons form free charges at AM1.5 light intensity. Bimolecular recombination, and hole traps constrain single component Y6 photovoltaics to low efficiencies, but recombination is reduced by small quantities of donor. Quantum-chemical calculations reveal strong coupling between exciton and CT states, and an intermolecular polarisation pattern that drives exciton dissociation. Our results challenge how current OPVs operate, and renew the possibility of efficient single-component OPVs.

When light hits organic semiconductors, bound charge pairs, called excitons, are usually produced. Here, the authors show that in the best performing organic solar material to date, free charges, rather than excitons, are directly created by light.

Virtual Screening of TADF Emitters for Single-Layer OLEDs

Thermally-activated delayed fluorescence (TADF) is a concept which helps to harvest triplet excitations, boosting the efficiency of an organic light-emitting diode. TADF can be observed in molecules with spatially separated donor and acceptor groups with a reduced triplet-singlet energy level splitting. TADF materials with balanced electron and hole transport are attractive for realizing efficient single-layer organic light emitting diodes, greatly simplifying their manufacturing and improving their stability. Our goal here is to computationally screen such materials and provide a comprehensive database of compounds with a range of emission wavelengths, ionization energies, and electron affinities.

Band gap engineering in blended organic semiconductor films based on dielectric interactions

Blending organic molecules to tune their energy levels is currently being investigated as an approach to engineer the bulk and interfacial optoelectronic properties of organic semiconductors. It has been proven that the ionization energy and electron affinity can be equally shifted in the same direction by electrostatic effects controlled by blending similar halogenated derivatives with different energetics. Here we show that the energy gap of organic semiconductors can also be tuned by blending. We use oligothiophenes with different numbers of thiophene rings as an example and investigate their structure and electronic properties. Photoelectron spectroscopy and inverse photoelectron spectroscopy show tunability of the single-particle gap, with the optical gaps showing similar, but smaller, effects. Theoretical analysis shows that this tuning is mainly caused by a change in the dielectric constant with blend ratio. Further studies will explore the practical impact of this energy-level engineering strategy for optoelectronic devices.

Computing Charging and Polarization Energies of Small Organic Molecules Embedded into Amorphous Materials with Quantum Accuracy

The ionization potential, electron affinity, and cation/anion polarization energies (IP, EA, P⁽⁺⁾, P^(-)) of organic molecules determine injection barriers, charge carriers balance, doping efficiency, and light outcoupling in organic electronics devices, such as organic light-emitting diodes (OLEDs). Computing IP and EA of isolated molecules is a common task for quantum chemistry methods. However, once molecules are embedded in an amorphous organic matrix, IP and EA values change, and accurate predictions become challenging. Here, we present a revised quantum embedding method [Friederich et al. J. Chem. Theory Comput. 2014, 10 (9), 3720-3725] that accurately predicts the dielectric permittivity and ionization potentials in three test materials, NPB, TCTA, and C60, and allows straightforward interpretation of their nature. The method paves the way toward reliable virtual screening of amorphous organic semiconductors with targeted IP/EA, polarization energies, and relative dielectric permittivity.

Intrinsic efficiency limits in low-bandgap non-fullerene acceptor organic solar cells

In bulk heterojunction (BHJ) organic solar cells (OSCs) both the electron affinity (EA) and ionization energy (IE) offsets at the donor-acceptor interface should equally control exciton dissociation. Here, we demonstrate that in low-bandgap non-fullerene acceptor (NFA) BHJs ultrafast donor-to-acceptor energy transfer precedes hole transfer from the acceptor to the donor and thus renders the EA offset virtually unimportant. Moreover, sizeable bulk IE offsets of about 0.5 eV are needed for efficient charge transfer and high internal quantum efficiencies, since energy level bending at the donor-NFA interface caused by the acceptors' quadrupole moments prevents efficient exciton-to-charge-transfer state conversion at low IE offsets. The same bending, however, is the origin of the barrier-less charge transfer state to free charge conversion. Our results provide a comprehensive picture of the photophysics of NFA-based blends, and show that sizeable bulk IE offsets are essential to design efficient BHJ OSCs based on low-bandgap NFAs.

Barrierless Free Charge Generation in the High-Performance PM6:Y6 Bulk Heterojunction Non-Fullerene Solar Cell

Organic solar cells are currently experiencing a second golden age thanks to the development of novel non-fullerene acceptors (NFAs). Surprisingly, some of these blends exhibit high efficiencies despite a low energy offset at the heterojunction. Herein, free charge generation in the high-performance blend of the donor polymer PM6 with the NFA Y6 is thoroughly investigated as a function of internal field, temperature and excitation energy. Results show that photocurrent generation is essentially barrierless with near-unity efficiency, regardless of excitation energy. Efficient charge separation is maintained over a wide temperature range, down to 100 K, despite the small driving force for charge generation. Studies on a blend with a low concentration of the NFA, measurements of the energetic disorder, and theoretical modeling suggest that CT state dissociation is assisted by the electrostatic interfacial field which for Y6 is large enough to compensate the Coulomb dissociation barrier.

A window to trap-free charge transport in organic semiconducting thin films

Organic semiconductors, which serve as the active component in devices, such as solar cells, light-emitting diodes and field-effect transistors¹, often exhibit highly unipolar charge transport, meaning that they predominantly conduct either electrons or holes. Here, we identify an energy window inside which organic semiconductors do not experience charge trapping for device-relevant thicknesses in the range of 100 to 300 nm, leading to trap-free charge transport of both carriers. When the ionization energy of a material surpasses 6 eV, hole trapping will limit the hole transport, whereas an electron affinity lower than 3.6 eV will give rise to trap-limited electron transport. When both energy levels are within this window, trap-free bipolar charge transport occurs. Based on simulations, water clusters are proposed to be the source of hole trapping. Organic semiconductors with energy levels situated within this energy window may lead to optoelectronic devices with enhanced performance. However, for blue-emitting light-emitting diodes, which require an energy gap of 3 eV, removing or disabling charge traps will remain a challenge.

Impact of molecular quadrupole moments on the energy levels at organic heterojunctions

The functionality of organic semiconductor devices crucially depends on molecular energies, namely the ionisation energy and the electron affinity. Ionisation energy and electron affinity values of thin films are, however, sensitive to film morphology and composition, making their prediction challenging. In a combined experimental and simulation study on zinc-phthalocyanine and its fluorinated derivatives, we show that changes in ionisation energy as a function of molecular orientation in neat films or mixing ratio in blends are proportional to the molecular quadrupole component along the π-π-stacking direction. We apply these findings to organic solar cells and demonstrate how the electrostatic interactions can be tuned to optimise the energy of the charge-transfer state at the donor-acceptor interface and the dissociation barrier for free charge carrier generation. The confirmation of the correlation between interfacial energies and quadrupole moments for other materials indicates its relevance for small molecules and polymers.

Read between the Molecules: Computational Insights into Organic Semiconductors

The performance and key electronic properties of molecular organic semiconductors are dictated by the interplay between the chemistry of the molecular core and the intermolecular factors of which manipulation has inspired both experimentalists and theorists. This Perspective presents major computational challenges and modern methodological strategies to advance the field. The discussion ranges from insights and design principles at the quantum chemical level, in-depth atomistic modeling based on multiscale protocols, morphological prediction and characterization as well as energy-property maps involving data-driven analysis. A personal overview of the past achievements and future direction is also provided.

Universal strategy for Ohmic hole injection into organic semiconductors with high ionization energies

Barrier-free (Ohmic) contacts are a key requirement for efficient organic optoelectronic devices, such as organic light-emitting diodes, solar cells, and field-effect transistors. Here, we propose a simple and robust way of forming an Ohmic hole contact on organic semiconductors with a high ionization energy (IE). The injected hole current from high-work-function metal-oxide electrodes is improved by more than an order of magnitude by using an interlayer for which the sole requirement is that it has a higher IE than the organic semiconductor. Insertion of the interlayer results in electrostatic decoupling of the electrode from the semiconductor and realignment of the Fermi level with the IE of the organic semiconductor. The Ohmic-contact formation is illustrated for a number of material combinations and solves the problem of hole injection into organic semiconductors with a high IE of up to 6 eV.

Hot and Cold Charge-Transfer Mechanisms in Organic Photovoltaics: Insights into the Excited States of Donor/Acceptor Interfaces

The evolution of the excited-state manifold in organic D/A aggregates (e.g., the prototypical P3HT/PCBM) is investigated through a bottom-up approach via first-principles calculations. We show how the excited-state energies, the charge transfer (CT) states, and the electron-hole density distributions are strongly influenced by the size, the orientation, and the position (i.e., on-top versus on-edge phases) of P3HT/PCBM domains. We discuss how the structural order influences the excited-state electronic structure, providing an atomistic interpretation of the photophysics of organic blends. We show how the simultaneous presence of on-top and on-edge phases does not alter the optical absorption spectrum of the blend but does affect the photophysics. Photovoltaic processes such as (i) the simultaneous charge generation obtained from hot and cold excitations, (ii) the instantaneous and delayed charge separation, and (iii) the pump-push-probe charge generation can be interpreted based on our study.

Study of the cold charge transfer state separation at the TQ1/PC71 BM interface

Charge transfer (CT) state separation is one of the most critical processes in the functioning of an organic solar cell. In this article, we study a bilayer of TQ1 and PC₇₁ BM molecules presenting disorder at the interface, obtained by means of Molecular Dynamics. The study of the CT state splitting can be first analyzed through the CT state splitting diagram, introduced in a previous work. Through this analysis, we identify the possibility of CT state splitting within Marcus Theory in function of the electric field. Once the right range of electric fields has been identified, we perform Kinetic Monte Carlo simulations to estimate percentages and times for the CT state splitting and the free charge carriers collection. Statistical information extracted from these simulations allows us to highlight the importance of polarization and to test the limits of the predictions given by the CT state splitting diagram. © 2017 Wiley Periodicals, Inc.