Virtual Screening for High Carrier Mobility in Organic Semiconductors

Sep 23, 2016 - Thus we readily find many acclaimed compounds as well as a number of most promising materials that have not yet been considered for an ...
0 downloads 10 Views 1MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Letter

Virtual Screening for High Carrier Mobility in Organic Semiconductors Christoph Schober, Karsten Reuter, and Harald Oberhofer J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.6b01657 • Publication Date (Web): 23 Sep 2016 Downloaded from http://pubs.acs.org on September 26, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Virtual Screening for High Carrier Mobility in Organic Semiconductors Christoph Schober, Karsten Reuter, and Harald Oberhofer∗ Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universi¨at M¨ unchen, Lichtenbergstr. 4, D-85747 Garching, Germany E-mail: [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 2 of 13

Page 3 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Electronic devices with components chiefly made from organic semiconductors raise great expectations, thanks to their novel material properties and favorable economic footprint. 1–6 Substantial experimental and theoretical resources have correspondingly been allocated to understand and identify organic semiconductors that are suitable for everyday application. 7–9 One of the key figures of merit for this is the mobility of charge carriers, which determines the efficiency of charge transport through the material and thus its applicability as light emitter, transistor or solar cell. Unfortunately, this intrinsic material property is difficult to access experimentally. Measured mobilities depend heavily on the employed device architecture 10,11 or device fabrication process, 12–14 and orders of magnitude improvements in carrier mobility have been reported for the same material simply by optimizing these factors. 15–17 While this highlights the relevance and prospects of an eventual device optimization, the limited access to the intrinsic mobility prevents the identification of suitable candidate materials for which such optimization efforts are to be initiated. In this situation highthroughput virtual screening approaches promise an intriguing and increasingly employed alternative. 18,19 These approaches assess the fitness of a (ideally huge) number of candidate materials by calculating computationally undemanding key quantities. The lower the computational cost to determine such descriptor(s), the more materials can be screened. The more representative the descriptor(s) for the actual figure of merit, the higher the likelihood to identify promising candidates. As the quantitative calculation of explicit charge transfer rates in real organic semiconductors is still demanding, 20 prominent descriptors for the mobility of charge carriers 21,22 are the transfer integral Hab and the reorganization energy λ. 23 Roughly speaking, the former mostly reflects the quantum mechanical overlap of the frontier orbitals involved in the charge transport, while the latter quantity accounts for the response of all other charges in the systems to the local change in charge state. While both quantities are intuitively understood in the context of hopping type transport, they are also meaningful descriptors beyond this regime. 24 The electronic coupling, for example, directly enters the tight-binding description of the carrier effective mass in the band-transport pic-

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ture 25,26 and λ can—albeit more indirectly—be viewed as a measure of the electron/phonon coupling strength 27 which e.g. plays a role in the relaxation time approximation of the band transport mobility. 28 Finally, another factor potentially influencing transport properties in organic semiconductors is the degree of disorder present in the system, 29 leading for example to a variation in site energies of localised charges on inequivalent sites. However, this effect is mostly of importance in amorphous systems, 17,30 and can thus for the sake of brevity be omitted for the crystalline semiconductors considered here (see supporting information for an analysis of site energy variations). In a pioneering screening study, Sokolov et al. 24 employed λ to assess the suitability of seven modifications to the well-known dinaphtho-thienothiophene (DNTT) organic semiconductor. 31 The choice for the reorganization energy was thereby not arbitrary, but motivated by the fact that λ can with some limitations be reduced to a single molecule property due to the comparatively low contribution of the external reorganisation energy. 32–34 This allows to screen candidate materials without the prerequisite of having to know their crystal structure. In comparison, the transfer integral Hab is solely defined in condensed phases and therefore requires knowledge about the crystal environment of the system. Moreover, prediction of charge transfer necessitates the determination of charge carrier percolation pathways with large couplings, 35 which is only possible given a full knowledge of the crystal geometry. Employing Hab in a virtual screening approach of unknown materials would therefore have to include a crystal structure prediction routine, which in itself is a huge and computationally generally involved challenge. 36,37 Of course, this disadvantage does not apply to already synthesized and structurally characterized materials. The Cambridge Structural Database (CSD) 38 alone contains for instance about 750.000 molecular crystal structures (CSD Summary Statistic, 16. February 2015), the by far largest part of which has never been investigated for a suitability in organic electronic devices. In this work we correspondingly devise a computational protocol to screen this compound library with a twofold objective. On the one hand, computing Hab and λ for

4

ACS Paragon Plus Environment

Page 4 of 13

Page 5 of 13

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

rise to a vast number of Hab calculations, which we efficiently handle on the level of fragmentorbital density-functional theory (DFT) 39,40 using the FHI-aims package 41 and organizing the computation tasks and workflows within the Fireworks-framework. 42 All compounds exhibiting at least one transfer integral that exceeds a minimal threshold value of 50 meV are subsequently recalculated with improved computational accuracy. If the minimum Hab value is confirmed, the reorganization energy of the system is calculated. This calculation relies on a QM/MM embedding scheme to account for differences between gas-phase and solid state structures, which can severely bias the predicted reorganization energies. 32 With further details provided in the SI, we thus arrive at an extensive (Hab , λ) descriptor database comprising a total of 10214 systems. Figure 2 depicts the 2D histogram of this database, together with the respective 1D distributions for both Hab and λ. Here we thereby concentrate on the maximum Hab value obtained for each crystal, viewing the largest coupling as performance critical for the system, while in the subsequent detailed analysis we also take percolation pathways into account. Ideal candidate systems promising a high intrinsic carrier mobility will exhibit a low reorganization energy λ and a large electronic coupling Hab . With increasing molecular size it is expected for Hab and λ to decrease due to the increasing delocalization of the wave function, 43 which is indeed confirmed for our dataset (see Fig. S5). A low reorganization energy of a system is therefore very likely to be nullified by a low transfer integral, while a large value of the transfer integral becomes meaningless if the associated reorganization energy is too high. The latter is less likely to occur due to the more favorable distribution of the reorganization energies with a clear maximum around 200 meV. The almost exponential decay of the Hab distribution instead shows that a low reorganization energy is very likely connected with a low transfer integral, rendering a virtual screening for low λs alone a less suitable approach to identify promising organic semiconductor materials. To identify the most promising candidates, we thus need to look for outliers to this trend, highlighting the need for an inclusive study as presented here.

6

ACS Paragon Plus Environment

Page 6 of 13

Page 7 of 13

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 13

Table 1: Selected compounds identified in Fig. 2 as promising p-type organic semiconductors. Compiled are the appealing (Hab , λ)-descriptor values, the crystal space group, as well as the low number of steps required to synthesize each compound from commercially available chemicals according to the information provided in the CSD database.

chem. formula Hab / meV λ / meV space group synthesis steps

S1

S2

S3

S4

C8 H9 NS2 190 169 P -1 0

C16 H12 Cl2 131 152 P 21/a 0

C22 H14 O2 120 167 P 21 21 21 1

C20 H14 Br2 85 121 P 21/n 3

however, also many more compounds in this area that exhibit similarly, if not even more promising descriptor values and that have apparently not yet been considered for an application in organic electronics. Deferring a more in depth exploration of this data to a forthcoming publication we present in Fig. 2 and Table 1 a small selection of corresponding materials. This selection has been based on the unique meta-data available for every compound in the CSD database, 38 like the ease of synthesis from commercially available chemicals, straightforward crystallization, or good stability. All selected structures also show at least one continuous percolation pathway with Hab values above 50 meV (see SI for details). Beyond such immediate leads, we nevertheless expect even larger insight and prospects from a subsequent mining of the generated large descriptor database. Analysis of trends and correlations, e.g. between crystal structure, percolation pathways, chemical backbone, type of functionalization, closest end-to-end distance in the crystal packaging, will clearly accelerate the theoretical design and discovery of future high mobility organic semiconductors. Computational Methods All DFT calculations were performed with the FHI-aims package. 41,45 Electronic wave functions were expanded in a tier-2 numeric atomic orbital (NAO) basis and tight integration grids, if not indicated otherwise. Total energies, electronic couplings and forces are calculated using the Becke exchange functional in combination with the 8

ACS Paragon Plus Environment

Page 9 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

correlation functional by Lee, Yang, and Parr (BLYP). 46,47 Scalar relativistic effects were considered on the level of the zero-order regular approximation (ZORA). Electronic coupling Hab : After constructing the nearest neighbour dimers using an inhouse code we compute electronic couplings for each dimer using the fragment molecular orbital approach (FO-DFT) 39 with the H 2n−1 @D+ A Reorganisation energy λ: A non-periodic QM/MM embedding scheme with a shell of neighbouring molecules is used to mimic the solid state environment during the geometry optimisation of a single (central) molecule. Following an ONIOM-scheme, 48 the MM-forces of the central molecule are subtracted from the full MM region. All atomic positions in the neighbour cell are constrained, and only the coordinates of the central molecule are allowed to change. In the MM calculations using the ”Universal Force Field” (UFF) 49 and the Qeq charge equilibration scheme 50 within the LAMMPS 51 framework, only the two-body interactions are considered. Ackknowledgement The authors gratefully acknowledge support from the Solar Technologies Go Hybrid initiative of the State of Bavaria and the Leibniz Supercomputing Centre for the use of the SuperMUC high performance computing facility. Supporting Information Available Detailed computational methods. Selection criteria and additional statistics for the presented dataset. Details and percolation pathways of suggested structures S1-S4.

References (1) Forrest, S. R. The Path to Ubiquitous and Low-Cost Organic Electronic Appliances on Plastic. Nature 2004, 428, 911–918. (2) Gundlach, D. J. Organic Electronics: Low Power, High Impact. Nat. Mater. 2007, 6, 173–174. (3) Minemawari, H.; Yamada, T.; Matsui, H.; Tsutsumi, J.; Haas, S.; Chiba, R.; Kumai, R.; Hasegawa, T. Inkjet Printing of Single-Crystal Films. Nature 2011, 475, 364–367.

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(4) Sun, Y.; Welch, G. C.; Leong, W. L.; Takacs, C. J.; Bazan, G. C.; Heeger, A. J. Solution-Processed Small-Molecule Solar Cells with 6.7% Efficiency. Nat. Mater. 2012, 11, 44–48. (5) Anikeeva, P.; Koppes, R. A. Restoring the Sense of Touch. Science 2015, 350, 274–275. (6) Stavrinidou, E.; Gabrielsson, R.; Gomez, E.; Crispin, X.; Nilsson, O.; Simon, D. T.; Berggren, M. Electronic Plants. Sci. Adv. 2015, 1, e1501136. (7) Anthony, J. E. Organic Electronics: Addressing Challenges. Nat. Mater. 2014, 13, 773–775. (8) Venkateshvaran, D.; Nikolka, M.; Sadhanala, A.; Lemaur, V.; Zelazny, M.; Kepa, M.; Hurhangee, M.; Kronemeijer, A. J.; Pecunia, V.; Nasrallah, I. et al. Approaching Disorder-Free Transport in High-Mobility Conjugated Polymers. Nature 2014, 515, 384–388. (9) Kim, G.; Kang, S.-J.; Dutta, G. K.; Han, Y.-K.; Shin, T. J.; Noh, Y.-Y.; Yang, C. A Thienoisoindigo-Naphthalene Polymer with Ultrahigh Mobility of 14.4 cm2/V·s That Substantially Exceeds Benchmark Values for Amorphous Silicon Semiconductors. J. Am. Chem. Soc. 2014, 136, 9477–9483. (10) Uno, M.; Uemura, T.; Kanaoka, Y.; Chen, Z.; Facchetti, A.; Takeya, J. High-Speed Organic Single-Crystal Transistors Gated with Short-Channel Air Gaps: Efficient Hole and Electron Injection in Organic Semiconductor Crystals. Org. Electron. 2013, 14, 1656–1662. (11) Orgiu, E.; George, J.; Hutchison, J. A.; Devaux, E.; Dayen, J. F.; Doudin, B.; Stellacci, F.; Genet, C.; Schachenmayer, J.; Genes, C. et al. Conductivity in Organic Semiconductors Hybridized with the Vacuum Field. Nat. Mater. 2015, 14, 1123–1129. (12) Reese, C.; Bao, Z. Organic Single-Crystal Field-Effect Transistors. Mater. Today 2007, 10, 20–27. (13) Zhou, Y.; Fuentes-Hernandez, C.; Shim, J.; Meyer, J.; Giordano, A. J.; Li, H.; Winget, P.; Papadopoulos, T.; Cheun, H.; Kim, J. et al. A Universal Method to Produce Low–Work Function Electrodes for Organic Electronics. Science 2012, 336, 327–332. (14) Yuan, Y.; Giri, G.; Ayzner, A. L.; Zoombelt, A. P.; Mannsfeld, S. C. B.; Chen, J.; Nordlund, D.; Toney, M. F.; Huang, J.; Bao, Z. Ultra-High Mobility Transparent Organic Thin Film Transistors Grown by an Off-Centre Spin-Coating Method. Nat. Commun. 2014, 5, 3005. (15) Goldmann, C.; Haas, S.; Krellner, C.; Pernstich, K. P.; Gundlach, D. J.; Batlogg, B. Hole Mobility in Organic Single Crystals Measured by a “Flip-Crystal” Field-Effect Technique. J. Appl. Phys. 2004, 96, 2080–2086. (16) Dong, H.; Fu, X.; Liu, J.; Wang, Z.; Hu, W. 25th Anniversary Article: Key Points for High-Mobility Organic Field-Effect Transistors. Adv. Mater. 2013, 25, 6158–6183. 10

ACS Paragon Plus Environment

Page 10 of 13

Page 11 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(17) Yuan, Y.; Giri, G.; Ayzner, A. L.; Zoombelt, A. P.; Mannsfeld, S. C. B.; Chen, J.; Nordlund, D.; Toney, M. F.; Huang, J.; Bao, Z. Ultra-high Mobility Transparent Organic Thin Film Transistors Grown by an Off-Centre Spin-Coating Method. Nat. Commun. 2014, 5, 3005. (18) Shoichet, B. K. Virtual Screening of Chemical Libraries. Nature 2004, 432, 862–865. (19) Pyzer-Knapp, E. O.; Suh, C.; G´omez-Bombarelli, R.; Aguilera-Iparraguirre, J.; AspuruGuzik, A. What Is High-Throughput Virtual Screening? A Perspective from Organic Materials Discovery. Annu. Rev. Mater. Res. 2015, 45, 195–216. (20) Coropceanu, V.; Li, H.; Winget, P.; Zhu, L.; Br´edas, J.-L. Electronic-Structure Theory of Organic Semiconductors: Charge-Transport Parameters and Metal/Organic Interfaces. Annu. Rev. Mater. Res. 2013, 43, 63–87. (21) Br´edas, J.-L.; Beljonne, D.; Coropceanu, V.; Cornil, J. Charge-Transfer and EnergyTransfer Processes in π-Conjugated Oligomers and Polymers: A Molecular Picture. Chem. Rev. 2004, 104, 4971–5004. (22) Sirringhaus, H. 25th Anniversary Article: Organic Field-Effect Transistors: The Path Beyond Amorphous Silicon. Adv. Mater. 2014, 26, 1319–1335. (23) Marcus, R. A. Electron Transfer Reactions In Chemistry. Theory and Experiment. Rev. Mod. Phys. 1993, 65, 599–610. (24) Sokolov, A. N.; Atahan-Evrenk, S.; Mondal, R.; Akkerman, H. B.; S´anchezCarrera, R. S.; Granados-Focil, S.; Schrier, J.; Mannsfeld, S. C. B.; Zoombelt, A. P.; Bao, Z. et al. From Computational Discovery to Experimental Characterization of a High Hole Mobility Organic Crystal. Nat. Commun. 2011, 2, 437. (25) Cheng, Y. C.; Silbey, R. J.; Filho, D. A. d. S.; Calbert, J. P.; Cornil, J.; Br´edas, J. L. Three-Dimensional Band Structure and Bandlike Mobility in Oligoacene Single Crystals: A Theoretical Investigation. J. Chem. Phys. 2003, 118, 3764–3774. (26) Troisi, A.; Orlandi, G. Band Structure of the Four Pentacene Polymorphs and Effect on the Hole Mobility at Low Temperature. J. Phys. Chem. B 2005, 109, 1849–1856. (27) Malagoli, M.; Coropceanu, V.; Filho, D. A. d. S.; Br´edas, J. L. A Multimode Analysis of the Gas-Phase Photoelectron Spectra in Oligoacenes. J. Chem. Phys. 2004, 120, 7490–7496. (28) Xi, J.; Long, M.; Tang, L.; Wang, D.; Shuai, Z. First-Principles Prediction of Charge Mobility in Carbon and Organic Nanomaterials. Nanoscale 2012, 4, 4348–4369. (29) B¨assler, H. Charge Transport in Disordered Organic Photoconductors a Monte Carlo Simulation Study. Phys. Status Solidi B 1993, 175, 15–56.

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(30) Hoffmann, S. T.; Jaiser, F.; Hayer, A.; B¨assler, H.; Unger, T.; Athanasopoulos, S.; Neher, D.; K¨ohler, A. How Do Disorder, Reorganization, and Localization Influence the Hole Mobility in Conjugated Copolymers? J. Am. Chem. Soc 2013, 135, 1772– 1782. (31) Yamamoto, T.; Takimiya, K. Facile Synthesis of Highly π-Extended Heteroarenes, Dinaphtho[2,3-b:2‘,3‘-f]chalcogenopheno[3,2-b]chalcogenophenes, and their Application to Field-Effect Transistors. J. Am. Chem. Soc. 2007, 129, 2224–2225. (32) Mas-Torrent, M.; Hadley, P.; Bromley, S. T.; Ribas, X.; Tarr´es, J.; Mas, M.; Molins, E.; Veciana, J.; Rovira, C. Correlation Between Crystal Structure and Mobility in Organic Field-Effect Transistors Based on Single Crystals of Tetrathiafulvalene Derivatives. J. Am. Chem. Soc. 2004, 126, 8546–8553. (33) Norton, J. E.; Br´edas, J.-L. Polarization Energies in Oligoacene Semiconductor Crystals. J. Am. Chem. Soc. 2008, 130, 12377–12384. (34) Kwiatkowski, J. J.; Frost, J. M.; Nelson, J. The Effect of Morphology on Electron Field-Effect Mobility in Disordered C60 Thin Films. Nano Lett. 2009, 9, 1085–1090. (35) Vehoff, T.; Baumeier, B.; Troisi, A.; Andrienko, D. Charge Transport in Organic Crystals: Role of Disorder and Topological Connectivity. J. Am. Chem. Soc. 2010, 132, 11702–11708. (36) Woodley, S.; Catlow, C. Crystal Structure Prediction From First Principles. Nat. Mater. 2008, 7, 937–946. (37) Tiekink, E. R. T. Molecular Crystals by Design? Chem. Commun. 2014, 50, 11079– 11082. (38) Allen, F. H. The Cambridge Structural Database: A Quarter of a Million Crystal Structures and Rising. Acta Crystallogr., Sect. B: Struct. Sci., Cryst. Eng. Mater. 2002, 58, 380–388. (39) Senthilkumar, K.; Grozema, F. C.; Bickelhaupt, F. M.; Siebbeles, L. D. A. Charge Transport in Columnar Stacked Triphenylenes: Effects of Conformational Fluctuations on Charge Transfer Integrals and Site Energies. J. Chem. Phys. 2003, 119, 9809–9817. (40) Schober, C.; Reuter, K.; Oberhofer, H. Critical Analysis of Fragment-Orbital DFT Schemes for the Calculation of Electronic Coupling Values. J. Chem. Phys. 2016, 144, 054103. (41) Blum, V.; Gehrke, R.; Hanke, F.; Havu, P.; Havu, V.; Ren, X.; Reuter, K.; Scheffler, M. Ab Initio Molecular Simulations with Numeric Atom-Centered Orbitals. Comput. Phys. Commun. 2009, 180, 2175–2196. (42) Jain, A.; Ong, S. P.; Chen, W.; Medasani, B.; Qu, X.; Kocher, M.; Brafman, M.; Petretto, G.; Rignanese, G.-M.; Hautier, G. et al. FireWorks: A Dynamic Workflow System Designed for High-Throughput Applications. Concurr. Comput. Pract. Exp. 2015, 5037–5059. 12

ACS Paragon Plus Environment

Page 12 of 13

Page 13 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(43) Zade, S. S.; Zamoshchik, N.; Bendikov, M. From Short Conjugated Oligomers to Conjugated Polymers. Lessons from Studies on Long Conjugated Oligomers. Acc. Chem. Res. 2011, 44, 14–24. (44) Takahashi, Y.; Hasegawa, T.; Horiuchi, S.; Kumai, R.; Tokura, Y.; Saito, G. High Mobility Organic Field-Effect Transistor Based on Hexamethylenetetrathiafulvalene with Organic Metal Electrodes. Chem. Mater. 2007, 19, 6382–6384. (45) Zhang, I. Y.; Ren, X.; Rinke, P.; Blum, V.; Scheffler, M. Numeric Atom-CenteredOrbital Basis Sets with Valence-Correlation Consistency from H to Ar. New J. Phys. 2013, 15, 123033. (46) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098–3100. (47) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula Into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785–789. (48) Svensson, M.; Humbel, S.; Froese, R. D. J.; Matsubara, T.; Sieber, S.; Morokuma, K. ONIOM: A Multilayered Integrated MO + MM Method for Geometry Optimizations and Single Point Energy Predictions. A Test for Diels-Alder Reactions and Pt(P(tBu)3)2 + H2 Oxidative Addition. J. Phys. Chem. 1996, 100, 19357–19363. (49) Rappe, A. K.; Casewit, C. J.; Colwell, K. S.; Goddard, W. A.; Skiff, W. M. UFF, a Full Periodic Table Force Field for Molecular Mechanics and Molecular Dynamics Simulations. J. Am. Chem. Soc. 1992, 114, 10024–10035. (50) Rappe, A. K.; Goddard, W. A. Charge Equilibration for Molecular Dynamics Simulations. J. Phys. Chem. 1991, 95, 3358–3363. (51) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19.

13

ACS Paragon Plus Environment