Charge carrier transport in organic semiconductors is at the heart of many revolutionary technologies ranging from organic transistors, light-emitting diodes, flexible displays and photovoltaic cells. Yet, the nature of charge carriers and their transport mechanism in these materials is still unclear. Here we show that by solving the time-dependent electronic Schrödinger equation coupled to nuclear motion for eight organic molecular crystals, the excess charge carrier forms a polaron delocalized over up to 10–20 molecules in the most conductive crystals. The polaron propagates through the crystal by diffusive jumps over several lattice spacings at a time during which it expands more than twice its size. Computed values for polaron size and charge mobility are in excellent agreement with experimental estimates and correlate very well with the recently proposed transient localization theory.
Progress in the design of high‐mobility organic semiconductors has been hampered by an incomplete fundamental understanding of the elusive charge carrier dynamics mediating electrical current in these materials. To address this problem, a novel fully atomistic non‐adiabatic molecular dynamics approach termed fragment orbital‐based surface hopping (FOB‐SH) that propagates the electron‐nuclear motion has been further improved and, for the first time, used to calculate the full 2D charge mobility tensor for the conductive planes of six structurally well characterized organic single crystals, in good agreement with available experimental data. The nature of the charge carrier in these materials is best described as a flickering polaron constantly changing shape and extensions under the influence of thermal disorder. Thermal intra‐band excitations from modestly delocalized band edge states (up to 5 nm or 10–20 molecules) to highly delocalized tail states (up to 10 nm or 40–60 molecules in the most conductive materials) give rise to short, ≈ 10 fs‐long bursts of the charge carrier wavefunction that drives the spatial displacement of the polaron, resulting in carrier diffusion and mobility. This study implies that key to the design of high‐mobility materials is a high density of strongly delocalized and thermally accessible tail states.
Germanane (GeH) and silicane (SiH) are the fully hydrogenated forms of germanene and silicene, the Ge- and Si-analogues of graphene. Here we use density-functional theory calculations to probe the properties of GeH and SiH sheets and their dependence on applied uni-axial compression. We find that GeH polymorphs with distinct hydrogen arrangements have markedly different energy band gaps. We also show that, when compressed, GeH forms superstructures with parts in low- and wide-gap geometries, enabling the creation of alternating polymorph nano-ribbons. An alternative to superstructure formation is the adoption by GeH of a corrugated form with extreme bending. Silicane shows similar behavior under compression, with either high corrugation, or successive parts with different geometries. Finally, we demonstrate that interaction with a substrate can influence the relative stability of GeH overlayer polymorphs. Overall, the results reveal ways to enhance the functionalities of these two-dimensional materials through the formation of superstructures with sizeable quantum well effects or outstanding mechanical response.
The search for conductive soft matter materials with significant charge mobility under ambient conditions has been a major priority in organic electronics (OE) research. Alkylated tetracenes are promising cost-effective candidate molecules that can be synthesized using wet chemistry methods, resulting in columnar single crystals with pronounced structural stability at and above room temperature. A remarkable characteristic of these materials is the capability of tuning the tetracene core intracolumnar stacking pattern and the crystal melting point via the side chain length and type modifications. In this study, we examine the performance of a series of alkylated tetracenes as hole conducting materials using a novel atomistic simulation technique that allows us to predict both the charge transport mechanism and mobilities. Our simulations demonstrate that molecular wires of alkylated tetracenes are capable of polaronic hole conduction at room temperature, with mobility values ranging up to 21 cm 2 V À1 s À1 , thus rendering such materials a highly promising choice for flexible OE applications. As regards the charge transfer robustness, two promising tetracene derivatives are identified with the capability of seamless inter-wire polaron delocalization, alleviating possible transfer bottlenecks due to local molecular defects. Our findings suggest that alkylated tetracenes offer an attractive route towards flexible columnar OE materials with unprecedented hole mobilities. † Electronic supplementary information (ESI) available: Description of the FOBSH methodology; systems under study; utilized force field and molecular dynamics simulations details; reorganization energy and charge transfer integral calculations details; classical force field validation results; inverse participation ratio time series for the TMT molecular crystal. See
Paper published as part of the special topic on 2D Materials Note: This paper is part of the JCP Special Topic on 2D Materials.
A new molecular dataset called HAB79 is introduced to provide ab initio reference values for electronic couplings (transfer integrals) and to benchmark density functional theory (DFT) and density functional tight-binding (DFTB) calculations. The HAB79 dataset is composed of 79 planar heterocyclic polyaromatic hydrocarbon molecules frequently encountered in organic (opto)electronics, arranged to 921 structurally diverse dimer configurations. We show that CASSCF/NEVPT2 with a minimal active space provides a robust reference method that can be applied to the relatively large molecules of the dataset. Electronic couplings are largest for cofacial dimers, in particular, sulfur-containing polyaromatic hydrocarbons, with values in excess of 0.5 eV, followed by parallel displaced cofacial dimers. V-shaped dimer motifs, often encountered in the herringbone layers of organic crystals, exhibit medium-sized couplings, whereas T-shaped dimers have the lowest couplings. DFT values obtained from the projector operator-based diabatization (POD) method are initially benchmarked against the smaller databases HAB11 (HAB7-) and found to systematically improve when climbing Jacob's ladder, giving mean relative unsigned errors (MRUEs) of 27.7% (26.3%) for the generalized gradient approximation (GGA) functional BLYP, 20.7% (15.8%) for hybrid functional B3LYP, and 5.2% (7.5%) for the long-range corrected hybrid functional omega-B97X. Cost-effective POD in combination with a GGA functional and very efficient DFTB calculations on the dimers of the HAB79 database give a good linear correlation with the CASSCF/NEVPT2 reference data, which, after scaling with a multiplicative constant, gives reasonably small MRUEs of 17.9% and 40.1%, respectively, bearing in mind that couplings in HAB79 vary over 4 orders of magnitude. The ab initio reference data reported here are expected to be useful for benchmarking other DFT or semi-empirical approaches for electronic coupling calculations.
Hydrogenated graphene, also known as graphane, is a wide band gap semiconductor in its various possible conformations. Here we show with first-principles calculations that application of uniaxial pressure can give rise to stripes of the so-called washboard conformation between segments of the most stable chair geometry. Compressed washboard-chair superlattices follow distinct zigzag patterns with interchanging chair and washboard parts. Hybrid systems of this type could be employed in various applications, such as in the templated growth of nanosystems, or as parts of nanomechanical systems.
The development of highly efficient methods for the calculation of electronic coupling matrix elements between the electron donor and acceptor is an important goal in theoretical organic semiconductor research. In Paper I [F. Gajdos, S. Valner, F. Hoffmann, J. Spencer, M. Breuer, A. Kubas, M. Dupuis, and J. Blumberger, J. Chem. Theory Comput. 10, 4653 (2014)], we introduced the analytic overlap method (AOM) for this purpose, which is an ultrafast electronic coupling estimator parameterized to and orders of magnitude faster than density functional theory (DFT) calculations at a reasonably small loss in accuracy. In this work, we reparameterize and extend the AOM to molecules containing nitrogen, oxygen, fluorine, and sulfur heteroatoms using 921 dimer configurations from the recently introduced HAB79 dataset. We find again a very good linear correlation between the frontier orbital overlap, calculated ultrafast in an optimized minimum Slater basis, and DFT reference electronic couplings. The new parameterization scheme is shown to be transferable to sulfur-containing polyaromatic hydrocarbons in experimentally resolved dimeric configurations. Our extension of the AOM enables high-throughput screening of very large databases of chemically diverse organic crystal structures and the application of computationally intense non-adiabatic molecular dynamics methods to charge transport in state-of-the-art organic semiconductors, e.g., non-fullerene acceptors.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.