phenacene single crystals - ACS Publications - American Chemical

Oct 10, 2018 - relaxed using the Quantum Espresso40 codes with the lattice parameters xed to the experimentally reported values.30,41. The cuto energy...
5 downloads 0 Views 1MB Size
Subscriber access provided by Gothenburg University Library

Communication

Surface structure of organic semiconductor [n]phenacene single crystals Yusuke Wakabayashi, Masanari Nakamura, Kaori Sasaki, Takahiro Maeda, Yuutaro Kishi, Hiroyuki Ishii, Nobuhiko Kobayashi, Susumu Yanagisawa, Yuma Shimo, and Yoshihiro Kubozono J. Am. Chem. Soc., Just Accepted Manuscript • DOI: 10.1021/jacs.8b08811 • Publication Date (Web): 15 Oct 2018 Downloaded from http://pubs.acs.org on October 16, 2018

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 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 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.

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 5 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

Journal of the American Chemical Society

Surface structure of organic semiconductor [n]phenacene single crystals Yusuke Wakabayashi, Hiroyuki Ishii,

†Division



∗,†

Masanari Nakamura,

Nobuhiko Kobayashi,





Kaori Sasaki,



Takahiro Maeda,



Susumu Yanagisawa, Kubozono

§

Yuma Shimo,

§



Yuutaro Kishi,



and Yoshihiro

of Materials Physics, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan

‡Department

of Applied Physics, Faculty of Pure and Applied Science, University of Tsukuba, Ibaraki

¶Department

of Physics and Earth Sciences, Faculty of Science, University of the Ryukyus, Okinawa

§Research

305-8573, Japan 903-0213, Japan

Institute for Interdisciplinary Science, Okayama University, Okayama, 700-8530, Japan Received October 10, 2018; E-mail: [email protected] u.ac.jp

Surface structure relaxation of organic semiconductors aects the properties of organic devices, although such relaxation has not been well explored. Only two examples have been experimentally reported; tetracene shows a large surface relaxation while rubrene exhibits no relaxation. Therefore, a systematic investigation of the surface relaxation is conducted on [n]phenacenes (n = 5, 7, and 9). Electron density analyses are performed based on the synchrotron surface X-ray scattering with the aid of rst-principles calculations. The results show little surface relaxation in [n]phenacenes. Abstract:

The recent development of organic semiconductor techniques allows us to design and fabricate various molecular electronic devices. Light emitting diodes and solar cells are commercially available, and laser devices have also been studied. 1 The physical properties of organic semiconductors are aected by the molecular arrangement. Such a structural dependence of the transport, optical, or mechanical properties has been studied from the point of view of polymorphism. 2,3 The distance and angle between neighboring molecules dominate the transfer integral. 46 Molecular vibrations in the lattice structure also aect the carrier conductivity owing to the electron-vibration couplings. 711 When we think about organic eld-eect transistors (OFETs), the carriers are accumulated in a few molecular layers at the surface. Therefore, the surface structure, rather than the crystal structure, is relevant to the OFET properties. It is well known that the surface structure is sometimes dierent from the bulk crystals; such structural modulation is called surface relaxation or surface reconstruction, and is one of the central issues of surface science. This can be regarded as surface-induced polymorphism. Theoretical calculations predict little surface relaxation in many cases, 1215 while the techniques for energy calculation of van der Waals crystals are in the development phase. 16,17 Experimentally, surface structures are often studied by scanning probe microscopy and grazing incidence X-ray diraction. The former, which observes the in-plane structure, reports no superstructure in many organic semiconductors. 1820 The latter is usually used on lms, and derives the lattice parameters and preferred orientation. 21,22 The molecular arrangement in the lm can also be studied by this technique, although it requires much eort and is less common. 2325 Although the best device performance is al-

ways found in single crystals, structural studies of the surface relaxation in organic semiconductor crystals are sparse. It has been reported 26,27 that a large surface relaxation is observed in a tetracene single crystal, and little relaxation is observed in a rubrene single crystal. Hence, it is important to clarify if surface relaxation is common in organic semiconductors, and determining the conditions that produce a large surface relaxation will be even more interesting. To clarify the general tendency of the surface relaxation, a systematic study of the surface structure is required. In this study, we performed a surface electron density analysis of [n]phenacenes (n = 5, 7, and 9) having herringbone structure, 28 which is similar to tetracene. [n]phenacenes exhibit a carrier mobility of 5 to 18 cm2 V−1 s−1 , 29,30 which is comparable with a rubrene OFET. The larger the value of n, the larger the mobility. The depth dependence of the crystal structure can be studied by the X-ray scattering from surfaces, 31,32 and the technique was applied also on organic crystals. 33 While there are several analyzing methods, 3436 only coherent Bragg rod analysis (COBRA) 37,38 has been successfully used to derive the depth prole of the electron density of organic semiconductor surfaces. 26,27 This method uses the scattering amplitude from the bulk structure part as the reference wave of the hologram to obtain the unknown scattering amplitude from the relaxed surfaces, which is treated as the object wave in holography. Therefore, we need the bulk crystal structure for all the samples. Unfortunately, [n]phenacenes with n ≥ 6 crystallize in the shape of a lm, which prevents us from performing single crystal X-ray structure analysis. Instead, we used the calculated structure. Surface electron density analysis based on the calculated bulk crystal structure has never been done, and is, therefore, a challenge in terms of both theoretical and analysis techniques. The rst-principles calculations were performed using the van der Waals density functional (vdW-DF) 39 with a planewave basis set. The internal degrees of freedom were fully relaxed using the Quantum Espresso 40 codes with the lattice parameters xed to the experimentally reported values. 30,41 The cuto energy for the plane wave was 120 Ry, and the Brillouin zone integration was performed with a 4 × 4 × 2 k-point set. For the vdW-DF family, rev-vdW-DF2 42 was employed. The atomic congurations of [n]phenacenes for n = 7 and 9 were calculated starting from the initial guess constructed from a similar molecular orientation as that in a n=5 crys-

ACS Paragon Plus Environment 1

Journal of the American Chemical Society

(a)

(b)

(c)

(a) Crystal structure of [5]phenacene obtained by single crystal structure analysis. 43 (b), (c) Crystal structures of [n]phenacenes (n=7 and 9) obtained by DFT calculation (see text). Figure 1.

Surface X-ray scattering proles for [n]phenecanes along the (00L) lines are presented in Fig. 2 by the open symbols. Note that the horizontal axis is expressed in reciprocal lattice unit (r.l.u.) for each crystal. As can be seen, the intensity proles were similar to each other. Sharp peaks at integer L positions are the Bragg peaks. A relatively strong signal at the anti-Bragg position was found at L = n + 1.5 for all [n]phenecanes. 1018 10

16

1014

datasets. The data treatment was the same as the rubrene and tetracene cases. 26,27 The eect of the interpolation was minor (see supporting information). The thick blue curves in Fig. 2 show the interpolated experimental results, which were used for electron density analysis. Electron density proles derived by COBRA for [n]phenacenes are presented in Fig. 3 by the black curves. The red curves represent the electron density proles calculated for the ideal surface. The similarity between them was clear. For n = 5, the topmost layer had smaller contrast in the electron density prole than the ideal surface, showing that the topmost molecules have larger positional uctuation than those in bulk. All n = 5, 7 and 9 proles show negligible amount of deviation from the ideal surface structure except for the uctuation. They showed a clear contrast with the electron density prole of tetracene presented in panel (d), which exhibited a large surface relaxation. Our experiment shows that the [n]phenacenes have little surface relaxation. Another nding is that no structural change below the second molecular layer has been observed so far, whereas the topmost layer often show larger positional uctuation or relaxation (old rubrene, 26 [5]phenacene and tetracene). The surface eect appears to be limited to the topmost layer, which may be a consequence of the weak inter-layer interaction. 0.8 -3 Electron density (Å )

tal structure. We used the experimental values of lattice parameters, 30,41 and assumed no additional symmetry for n = 7 and 9 samples. The resulting DFT-derived structures are presented in Fig. 1 together with the results of the single crystal structure analysis of [5]phenacene. 43

1012 10

10

n=9

108

n=7

106 Expt. Si subst. Ideal model Mod. bulk str. COBRA

104 10

2

n=5

100

0

1

2

3

4

5

6

7

8

9 10 11 12 13

L(r.l.u) Figure 2. Surface X-ray scattering intensity distribution of [n]phenacenes along the 00L-rod. The horizontal axis is expressed in reciprocal lattice unit for each crystal. The signal from the Si sample holder for n = 9 data is also ploted with gray symbols. The blue curves show the interpolated experimental data for COBRA, and the red curves show the calculated intensity proles based on the crystal structure (n = 5) or DFT-derived structures (n = 7 and 9). The red dashed curve shows the calculated intensity distribution from modied bulk crystal structure (see text.)

The thin red curves show the calculated intensity proles from ideal surfaces, that is, the sudden termination of the bulk structure without any distortion. The observed proles were well reproduced by this calculation for n = 5 and n = 7. For n = 9, slight modication of the bulk crystal structure, 3.5◦ rotation of the molecules from the DFT-derived structure, improved the similarity between the measured and calculated intensity (red dashed curve, see supporting information for detail). Those features showed that the degree of surface relaxation was small. Coherent Bragg rod analysis was performed using a homemade code, which has already been successfully used for rubrene and tetracene, 26,27 for interpolated experimental

Electron density (Å-3 )

Intensity(a.u.)

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

(a) n = 5

(b) n = 7

(c) n = 9

(d) tetracene

0.6 0.4 0.2 0.0 0.8

COBRA Ideal surface

0.6 0.4 0.2 0.0

-10

0

10

20 30 z(Å)

40

50

-10

0

10

20 30 z(Å)

40

50

60

Electron density prole as a function of depth for [n]phenacenes together with that of tetracene. 27 The origin of the horizontal axis is chosen as the surface. The oscillation of the electron density above the surface (negative-z region) is caused by the termination error of the Fourier transformation. Figure 3.

The surface X-ray scattering is caused by not only the surface but also the interface between the phenacene and SiO2 sample holder. Previous studies have shown that the interface is either very rough or as at as the surface. 26,27,44 The large FET mobility of phenacenes 29,30 suggests that the interface is at. The surface scattering intensity prole for each phenacene was well explained by one specic surface structure presented in Fig. 3, showing that the interface structure was the same as the surface structure. Let us examine what makes the surface relaxed. Morisaki 27 et al. pointed out that pentacene, which has a very similar structure as tetracene, has various kinds of polymorph structures, which suggests that the surface relaxation in tetracene is caused by the various structures having similar energy in the bulk form. The phenacenes reported in this study crystallized in a specic structure very stably. At this stage, we propose that the number of possible polymorph structures for the bulk crystal structure reects the possible degree of

ACS Paragon Plus Environment 2

Page 2 of 5

Page 3 of 5 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

Journal of the American Chemical Society

surface relaxation. Recently, the FET mobility, which is what we can measure experimentally, of phenacenes was calculated 45 based on the bulk crystal structure by using Marcus theory, and the values were close to the experimentally reported ones. 29,30 Although the Marcus theory has been used for ages, modern understanding of the transport in organic semiconductors does not support using it. 46,47 Theoretical examination based on modern transport theory for organic semiconductor is desired. On the other hand, using the bulk crystal structure to estimate the surface transport properties for phenacenes is supported by our present work. Surface structure relaxation of [n]phenacenes is examined by the combination of rst-principles calculations and surface X-ray scattering measurements. The phase retrieval analysis of the surface X-ray scattering using the calculated bulk structure was successfully demonstrated. As a result, [n]phenacenes exhibit little surface relaxation, which is very dierent from the tetracene case.

Acknowledgement The authors thank Mr. Y. Araya and Ms. S. Hamao for the experimental support. This work was supported by Grants-In-Aid for Scientic Research (JSPS KAKENHI, Grant Nos JP26287080, JP26105004, JP26105008, and JP26105011). The synchrotron radiation experiments at the Photon Factory were performed with the approval of the Photon Factory Program Advisory Committee (Proposal No. 2015S2-009). The computations were performed at the Supercomputer Center, ISSP, University of Tokyo, and Center for Computational Sciences, University of Tsukuba.

(11) (12) (13) (14) (15) (16) (17) (18) (19) (20)

(21) (22)

Supporting Information Available Experimental procedures and detail of the electron density analysis are provided in the supporting information. The Supporting Information is available free of charge on the ACS Publications website at DOI:

References

(23) (24) (25) (26)

(1) Samuel, I.; Turnbull, G. Organic semiconductor lasers. Chem. Rev. 2007, 107, 12721295. (2) Chung, H.; Diao, Y. Polymorphism as an emerging design strategy for high performance organic electronics. J. Mater. Chem. C 2016, 4, 39153933. (3) Jones, A.; Chattopadhyay, B.; Geerts, Y.; Resel, R. Substrateinduced and thin-lm phases: Polymorphism of organic materials on surfaces. Adv. Funct. Mater. 2016, 26, 22332255. (4) Mori, T.; Kobayashi, A.; Sasaki, Y.; Kobayashi, H.; Saito, G.; Inokuchi, H. Bull. Chem. Soc. Jpn. 1984, 57, 627633. (5) Da Silva Filho, D.; Kim, E.-G.; Brédas, J.-L. Transport properties in the rubrene crystal: Electronic coupling and vibrational reorganization energy. Adv. Mater. 2005, 17, 10721076. (6) Li, Y.; Coropceanu, V.; Brédas, J.-L. Thermal narrowing of the electronic bandwidths in organic molecular semiconductors: Impact of the crystal thermal expansion. J. Phys. Chem. Lett. 2012, 3, 33253329. (7) Valeev, E.; Coropceanu, V.; Da Silva Filho, D.; Salman, S.; Brédas, J.-L. Eect of electronic polarization on charge-transport parameters in molecular organic semiconductors. J. Am. Chem. Soc. 2006, 128, 98829886. (8) Troisi, A.; Orlandi, G. Charge-Transport Regime of Crystalline Organic Semiconductors: Diusion Limited by Thermal O-Diagonal Electronic Disorder. Phys. Rev. Lett. 2006, 96, 086601. (9) Ishii, H.; Honma, K.; Kobayashi, N.; Hirose, K. Wave-packet approach to transport properties of carrier coupled with intermolecular and intramolecular vibrations of organic semiconductors. Phys. Rev. B 2012, 85 . (10) Ishii, H.; Kobayashi, N.; Hirose, K. Charge transport calculations by a wave-packet dynamical approach using maxi-

(27) (28) (29) (30) (31) (32) (33) (34) (35) (36)

mally localized Wannier functions based on density functional theory: Application to high-mobility organic semiconductors. Phys. Rev. B 2017, 95, 035433. Fratini, S.; Ciuchi, S.; Mayou, D.; De Laissardiére, G.; Troisi, A. A map of high-mobility molecular semiconductors. Nat. Mater. 2017, 16, 9981002. Marcon, V.; Raos, G. Free energies of molecular crystal surfaces by computer simulation: Application to tetrathiophene. J. Am. Chem. Soc. 2006, 128, 14081409. Massaro, F.; Moret, M.; Bruno, M.; Rubbo, M.; Aquilano, D. Equilibrium and growth morphology of oligoacenes: Periodic bond chains (PBC) analysis of tetracene crystal. Cryst. Growth Des. 2011, 11, 46394646. Massaro, F.; Moret, M.; Bruno, M.; Aquilano, D. Equilibrium and growth morphology of oligoacenes: Periodic bond chains analysis of naphthalene, anthracene, and pentacene crystals. Cryst. Growth Des. 2012, 12, 982989. Massaro, F.; Moret, M.; Bruno, M.; Aquilano, D. Equilibrium and growth morphology of oligothiophenes: Periodic bond chain analysis of quaterthiophene and sexithiophene crystals. Cryst. Growth Des. 2013, 13, 13341341. Yanagisawa, S.; Okuma, K.; Inaoka, T.; Hamada, I. Recent progress in predicting structural and electronic properties of organic solids with the van der Waals density functional. J. Electron. Spectrosc. Relat. Phenom. 2015, 204, 159167. Hermann, J.; DiStasio, J., R.A.; Tkatchenko, A. First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications. Chem. Rev. 2017, 117, 47144758. Overney, R.; Howald, L.; Frommer, J.; Meyer, E.; Güntherodt, H.-J. Molecular surface structure of tetracene mapped by the atomic force microscope. J. Chem. Phys. 1991, 94, 8441 8443. Minato, T.; Aoki, H.; Fukidome, H.; Wagner, T.; Itaya, K. Highresolution molecular images of rubrene single crystals obtained by frequency modulation atomic force microscopy. Appl. Phys. Lett. 2009, 95, 093302. Casalegno, M.; Moret, M.; Resel, R.; Raos, G. Surface Reconstructions in Organic Crystals: Simulations of the Eect of Temperature and Defectivity on Bulk and (001) Surfaces of 2,2 ′ :6 ′ ,2 ″ -Ternaphthalene. Cryst. Growth Des. 2016, 16, 412422. Li, J.; Zhao, Y.; Tan, H.; Guo, Y.; Di, C.-A.; Yu, G.; Liu, Y.; Lin, M.; Lim, S.; Zhou, Y.; Su, H.; Ong, B. A stable solutionprocessed polymer semiconductor with record high-mobility for printed transistors. Sci. Rep. 2012, 2 . Verploegen, E.; Miller, C.; Schmidt, K.; Bao, Z.; Toney, M. Manipulating the morphology of P3HT-PCBM bulk heterojunction blends with solvent vapor annealing. Chem. Mater. 2012, 24, 39233931. Mannsfeld, S.; Virkar, A.; Reese, C.; Toney, M.; Bao, Z. Precise structure of pentacene monolayers on amorphous silicon oxide and relation to charge transport. Adv. Mater. 2009, 21, 2294 2298. Mannsfeld, S.; Tang, M.; Bao, Z. Thin lm structure of triisopropylsilylethynyl-functionalized pentacene and tetraceno[2,3-b]thiophene from grazing incidence x-ray diraction. Adv. Mater. 2011, 23, 127131. Shirasawa, T.; Yanagisawa, S.; Hatada, S.-N.; Voegeli, W.; Morikawa, Y.; Takahashi, T. A New Pentacene Polymorph Induced by Interaction with a Bi(0001) Substrate. J. Phys. Chem. C 2018, 122, 62406245. Wakabayashi, Y.; Takeya, J.; Kimura, T. Sub-Å Resolution Electron Density Analysis of the Surface of Organic Rubrene Crystals. Phys. Rev. Lett. 2010, 104, 066103. Morisaki, H.; Koretsune, T.; Hotta, C.; Takeya, J.; Kimura, T.; Wakabayashi, Y. Large surface relaxation in the organic semiconductor tetracene. Nat. Commun. 2014, 5, 5400. Desiraju, G.; Gavezzotti, A. Crystal structures of polynuclear aromatic hydrocarbons. Classication, rationalization and prediction from molecular structure. Acta Crystallogr., Sect. B 1989, 45, 473482. Kubozono, Y.; He, X.; Hamao, S.; Teranishi, K.; Goto, H.; Eguchi, R.; Kambe, T.; Gohda, S.; Nishihara, Y. Transistor application of phenacene molecules and their characteristics. Eur. J. Inorg. Chem. 2014, 2014, 38063819. Shimo, Y.; Mikami, T.; Hamao, S.; Goto, H.; Okamoto, H.; Eguchi, R.; Gohda, S.; Hayashi, Y.; Kubozono, Y. Synthesis and transistor application of the extremely extended phenacene molecule, [9]phenacene. Sci. Rep. 2016, 6 . Andrews, S.; Cowley, R. Scattering of x-rays from crystal surfaces. J. Phys. C: Solid State Phys. 1985, 18, 64276439. Robinson, I. Crystal truncation rods and surface roughness. Phys. Rev. B 1986, 33, 38303836. Gidalevitz, D.; Feidenhans'l, R.; Smilgies, D.; Leiserowitz, L. X-ray scattering from surfaces of organic crystals. Surf. Rev. Lett. 1997, 4, 721732. Fienup, J. Phase retrieval algorithms: A comparison. Appl. Opt. 1982, 21, 27582769. Saldin, D.; Harder, R.; Shneerson, V.; Moritz, W. Phase retrieval methods for surface x-ray diraction. J. Phys.: Condens. Matter 2001, 13, 1068910707. Björck, M.; Schlepütz, C.; Pauli, S.; Martoccia, D.; Herger, R.; Willmott, P. Atomic imaging of thin lms with surface x-ray diraction: Introducing DCAF. J. Phys.: Condens. Matter

ACS Paragon Plus Environment 3

Journal of the American Chemical Society 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

2008, 20, 445006. (37) Yacoby, Y.; Sowwan, M.; Stern, E.; Cross, J.; Brewe, D.; Pindak, R.; Pitney, J.; Dufresne, E.; Clarke, R. Direct determination of epitaxial interface structure in Gd2 O3 passivation of GaAs. Nat. Mater. 2002, 1, 99101. (38) Sowwan, M.; Yacoby, Y.; Pitney, J.; MacHarrie, R.; Hong, M.; Cross, J.; Walko, D.; Clarke, R.; Pindak, R.; Stern, E. Direct atomic structure determination of epitaxially grown lms: Gd2 O3 on GaAs(100). Phys. Rev. B 2002, 66, 2053111 20531112. (39) Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I. Van der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92 . (40) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.; Cococcioni, M.; Dabo, I.; Dal Corso, A.; De Gironcoli, S.; Fabris, S.; Fratesi, G.; Gebauer, R.; Gerstmann, U.; Gougoussis, C.; Kokalj, A.; Lazzeri, M.; Martin-Samos, L.; Marzari, N.; Mauri, F.; Mazzarello, R.; Paolini, S.; Pasquarello, A.; Paulatto, L.; Sbraccia, C.; Scandolo, S.; Sclauzero, G.; Seitsonen, A.; Smogunov, A.; Umari, P.; Wentzcovitch, R. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J. Phys.: Condens. Matter 2009, 21 . (41) He, X.; Eguchi, R.; Goto, H.; Uesugi, E.; Hamao, S.; Takabayashi, Y.; Kubozono, Y. Fabrication of single crystal eldeect transistors with phenacene-type molecules and their excellent transistor characteristics. Org. Electron. 2013, 14, 1673  1682. (42) Hamada, I. van der Waals density functional made accurate. Phys. Rev. B 2014, 89, 121103. (43) De, A.; Ghosh, R.; Roychowdhury, S.; Roychowdhury, P. Structural analysis of picene, C2 2H1 4. Acta Crystallogr., Sect. C 1985, 41, 907909. (44) Wakabayashi, Y.; Takeya, J.; Kimura, T. Interface structure of the rubrene crystal eld eect transistor. J. Appl. Phys. 2011, 110 . (45) Nguyen, T.; Roy, P.; Shim, J. Remarkable charge-transfer mobility from [6] to [10]phenacene as a high performance p-type organic semiconductor. Phys. Chem. Chem. Phys. 2018, 20, 86588667. (46) Fratini, S.; Mayou, D.; Ciuchi, S. The transient localization scenario for charge transport in crystalline organic materials. Adv. Funct. Mater. 2016, 26, 22922315. (47) Oberhofer, H.; Reuter, K.; Blumberger, J. Charge Transport in Molecular Materials: An Assessment of Computational Methods. Chem. Rev. 2017, 117, 1031910357.

ACS Paragon Plus Environment 4

Page 4 of 5

Page 5 of 5

Graphical TOC Entry

108

0.7 0.6

Exp. Ideal COBRA

107

n=5

106

0.5

105

Intensity

Electron density(Å-3)

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

Journal of the American Chemical Society

0.4 0.3

104 103 102

0.2

101

0.1

100

0.0

10-1

-13.52

0.00

13.52

27.04

z(Å)

40.56

54.08

0

2

4

6

00L

ACS Paragon Plus Environment 5

8

10