Search tips
Search criteria 


Logo of nanoreslettspringer open web sitethis articlemanuscript submissionregistrationjournal front pagespringer open web site
Nanoscale Res Lett. 2012; 7(1): 248.
Published online 2012 May 14. doi:  10.1186/1556-276X-7-248
PMCID: PMC3479028

Different interface orientations of pentacene and PTCDA induce different degrees of disorder


Organic polymers or crystals are commonly used in manufacturing of today‘s electronically functional devices (OLEDs, organic solar cells, etc). Understanding their morphology in general and at the interface in particular is of paramount importance. Proper knowledge of molecular orientation at interfaces is essential for predicting optoelectronic properties such as exciton diffusion length, charge carrier mobility, and molecular quadrupole moments. Two promising candidates are pentacene and 3,4:9,10-perylenetetracarboxylic dianhydride (PTCDA). Different orientations of pentacene on PTCDA have been investigated using an atomistic molecular dynamics approach. Here, we show that the degree of disorder at the interface depends largely on the crystal orientation and that more ordered interfaces generally suffer from large vacancy formation.

Keywords: Organic interfaces, Organic electronic devices, Interface disorder, Molecular dynamics, PTCDA, Pentacene


Organic light emitting diodes (OLEDs), organic solar cells, organic thin films transistors, etc. are made of organic polymers or crystals [1-3]. The effect of the disorder in organic devices on optoelectronic properties was analyzed by Rim et al. [4]. They showed an increased photocurrent generation with improved molecular order. It occurs due to the influence of the stacking on the exciton diffusion length. Hu et al. measured a strong dependence of the conductance across highly oriented pentacene nanocrystals on the packing orientation [5]. The influence of packing on charge transport in organic solids was also analyzed using Monte Carlo methods [6]. Kwiatkowski et al. [6] were able to predict the mobilities of electron and holes for ordered and disordered Alq3. Different functional organic materials were reviewed by Ishii et al. [7]. They highlighted the energy level alignment and electronic structures at organic/inorganic and organic/organic interfaces of, for example, Alq3, 3,4:9,10-perylenetetracarboxylic dianhydride (PTCDA) and 1,4,5,8-tetrathiafulvalene (TTF).

In our work, the morphology of interfaces between pentacene [8] and PTCDA [9] was analyzed (Figure (Figure1a).1a). Both molecules form different crystal modifications. Pentacene is known to have a high temperature (HT) and a low temperature (LT) polymorph. Yoneya et al. [8] showed that the LT polymorph is destabilized by substrates and transforms into HT polymorph. Therefore, the HT polymorph was used as the base for simulations. For PTCDA, the αpolymorph [9] was used.

Figure 1
Pentacene and PTCDA: Chemical formulas and interface formation. (a) Chemical formulas of PTCDA (top) and pentacene (bottom) are presented. (b) Example of the realistic interface formed by PTCDA (212) on pentacene (100). After full MD relaxation cycle, ...

Molecular orientation at interfaces is decisive for predicting optoelectronic properties such as exciton diffusion length [10], charge carrier mobility [11], and molecular quadrupole moments [12]. Verlaak et al. analyzed the impact of the molecular quadrupole moments, influenced by e. g., material and crystal orientation on the interface energetics. An insight on models of electronic processes across organic interfaces is given by Beljonne et al. [13], while a review of the corresponding theoretical approaches is presented by Brédas [14].

Our study of organic-organic pentacene/PDCDA interfaces is organized as follows: after a brief introduction presented above, we proceed with the presentation of the methods followed by the results and some conclusive remarks.


The molecular dynamic (MD) simulations of the interfaces between PTCDA and pentacene have been performed with the atomistic molecular dynamics package GROMACS (Stockholm Center for Biomembrane Research, Stockholm, Sweden and Biomedical Centre, Uppsala, Sweden) [15] using the generalized amber force field (GAFF) parameterization [16] for organic molecules, having Yoneya et al.’s work [8] in mind, and ESP charges [17] calculated with the semi-empirical quantum chemistry package MOPAC (Stewart Computational Chemistry, Colorado Springs, CO, USA) [18]. The parameter conversion from amber to GROMACS was done with the help of Antechamber python parser interface (ACPYPE) [19], the recommended tool for using GAFF with GROMACS, cf [8,20-22]. After simulation, a check of basic molecule parameters was done and the results for the example of pentacene are presented in Table Table1.1. A more detailed report on relative errors in energy, dehidrals, etc can be found in the ACPYPE wiki [23].

Table 1
Comparison of calculated and experimental relevant parameters

The systems were simulated with a step size of 0.5 fs for more than 3 ns at a temperature of 300 K using a Berendsen thermostat [26] for temperature control. The van der Waals cut-off was set to 1.2 nm, the Coulomb cut-off to 5 nm and the relative permittivity was set to four which was taken from Wang et al. [27]. No periodic boundary conditions were used owing to the different crystal lattices.

Three surfaces were chosen and combined. For pentacene the surfaces (100), (010), and (001) were used and for PTCDA the surfaces used are (102), (-221), and (212) as defined by Miller indices. The combination of these surfaces led to nine different interface facets, e. g., (212) on (010) and (-221) on (001), as depicted in Figure Figure1b,c1b,c showing their relaxed structures, leaving rotation and translation as degrees of freedom. An optimal relative orientation within each of these nine facets was found by performing four simulations each with relative orientations from being twisted against each other. After a total energy comparison, the structure with the lowest mean energy per molecule of the fully relaxed systems was chosen. As an example, the energy-evolution for the interface facet (-221)//(100) is shown in Figure Figure2.2. The set of simulations were done on systems arranged to fill a 10 × 10 × 10 nm3 cube with each crystal type, filling half the space.

Figure 2
Time-evolution of mean energy per molecule for the four interfaces of (-221) PTCDA and (100) pentacene. The triangles mark the mean energy at subsequent time steps where each relative orientation is represented by a different color. After few hundred ...

Results and discussion

In order to quantify the disorder at the pentacene/PTCDA interface, we used distribution of ϕ, defined as the angle between the molecular and the interface plane (or rather their respective normals) as shown in Figure Figure3.3. Owing to the fact that the molecules will start to relax, they will start to deviate from the bulk values. The more molecules have different ϕ, the more disordered is the structure.

Figure 3
Definition of the angle defining the molecular orientation along with a distribution at one interface. (a) ϕis the angle between the normal of the molecule plane and the normal of the interface plane (z-axis), i. e. the angle between the molecule ...

In the histograms of Figure Figure4,4, the y-axis was defined as distance in Å from the (ideal) interface in z-direction, while the x-axis shows the angle distribution. Light blue regions mark the disordered regions. Two clear patterns can be observed: 1) size of the disordered region can vary from 2 to 16 Å, and 2) the disorder seems to spread asymmetrically from the ideal interface, clearly preferring pentacene-rich regions. The first pattern can be explained as having two competing effects at the interface, one being the optimization of the intermolecular distance/interaction and the other being the conservation of bulk properties. The second pattern can be understood in the light of much stronger ΠΠ stacking of the PTCDA molecules, leading to a stronger intermolecular interactions, and greater energies are required to disrupt these molecules from their bulk positions when compared to pentacene bulk.

Figure 4
Angular distribution of angleϕas function of distance from the interface. The 2d-histograms (a–i) show the angle ϕdepending on the distance to the interface at z=0, where ϕis the angle between the molecule plane and the ...


Analysis of PTCDA/pentacene interfaces was performed with two emerging messages: there seems to be two competing effects, one coming from intermolecular interaction, which leads to disordered interfaces, while the other coming from the preservation of bulk properties results in large interfacial vacancies. Both of the effects would lead to dramatically diminished transport properties. Namely, increased disorder would cause greater energy disorder of the interfacial hopping sites, while interfacial vacancies would lead to diminished intermolecular overlaps, or hopping matrix elements. Whether which of the competing effects is influencing more the hopping transport properties is the focus of our ongoing research. Our second observation is that pentacene seems to be, in general, a more flexible material, which can be observed from the fact that the disordered regions are predominantly pentacene-rich.

Competing interests

The authors declare that they have no competing interests.

Author’s contributions

AP carried out the molecular dynamics calculations, the setup of the initial system and helped in drafting of the manuscript, and revisions. VM helped in analysis and interpretation of data, and drafted the manuscript and revisions. RM provided the calculation of the partial charges. WW participated in the design of the study, formulated the original scientific question and helped in analysis and interpretation of data. All authors read and approved the final manuscript.

Authors’ information

AP is Ph.D. student, VM and RM have Ph.D. degree in physics, and WW is an associate professor at Karlsruhe Institute of Technology.


This work was supported by bwGRiD resources and the FP7 project MINOTOR. bwGRiD is the central collection of computing power within the state of Baden-Wuerttemberg operated by eight universities, providing access for local researchers. Further thanks go to Ivan Kondov from SCC/KIT.


  • Yun C, Cho H, Kang H, Lee Y. Electron injection via pentacene thin films for efficient inverted organic light-emitting diodes. Appl Phys Lett. 2009;95:053301. doi: 10.1063/1.3192361. [Cross Ref]
  • Roncaliu J. Molecular bulk heterojunctions: an emerging approach to organic solar cells. Acc Chem Res. 2009;42:1719–1730. doi: 10.1021/ar900041b. [PubMed] [Cross Ref]
  • Reese Bao. Organic single-crystal field-effect transistors. Mat Today. 2007;10:20–27.
  • Rim S, Fink R, Schöneboom J, Erk P, Peumans P. Effect of molecular packing on the exciton diffusion length in organic solar cells. Appl Phys Lett. 2007;91:173504. doi: 10.1063/1.2783202. [Cross Ref]
  • Hu W-S, Tao Y-T, Chen Y-F, Chang C-S. Orientation-dependent conductance study of pentacene nanocrystals by conductive atomic force microscopy. Appl Phys Lett. 2008;93:053304. doi: 10.1063/1.2960343. [Cross Ref]
  • Kwiatkowski J, Nelson J, Li J, Brédas J, Wenzel W, Lennartz C. Simulating charge transport in tris (8-hydroxyquinoline) aluminium (Alq3) Phys Chem Chem Phys. 2008;10:1852–1858. [PubMed]
  • Ishii H, Sugiyama K, Ito E, Seki K. Energy level alignment and interfacial electronic structures at organic/metal and organic/organic Interfaces. Adv Mater. 1999;11:8.
  • Yoneya M, Kawasaki M, Ando M. Molecular dynamics simulations of pentacene thin films: The effect of surface on polymorph selection. J Mater Chem. 2010;20:10397–10402. doi: 10.1039/c0jm01577f. [Cross Ref]
  • Ogawa T, Kuwamoto K, Isoda S, Kobayashi T, Karl N. 3,4:9,10-Perylenetetracarboxylic dianhydride (PTCDA) by electron crystallography. Acta Cryst B. 1999;55:123–130. doi: 10.1107/S0108768198009872. [PubMed] [Cross Ref]
  • Najafov H, Lee B, Zhou Q, Feldman L, Podzorov V. Observation of long-range exciton diffusion in highly ordered organic semiconductors. Nat Mater. 2010;9:938–943. doi: 10.1038/nmat2872. [PubMed] [Cross Ref]
  • Vehoff T, Baumeier B, Troisi A, Andrienko D. Charge transport in organic crystals: role of disorder and topological connectivity. J Am Chem Soc. 2010;13:11702–11708. [PubMed]
  • Verlaak S, Beljonne D, Cheyns D, Rolin C, Linares M, Castet F, Cornil J, Heremans P. Electronic structure and geminate pair energetics at organic–organic interfaces: the case of pentacene/C60 heterojunctions. Adv Func Mat. 2009;19:3809–3814. doi: 10.1002/adfm.200901233. [Cross Ref]
  • Beljonne D, Cornil J, Muccioli L, Zannoni CJ, Castet F. Electronic processes at organic-organic interfaces: insight from modeling and implications for opto-electronic devices. Chem Mater. 2011;23:591–609. doi: 10.1021/cm1023426. [Cross Ref]
  • Brédas J, Norton J, Cornil J, Coropceanu V. Molecular understanding of organic solar cells: the challenges. Acc Chem Res. 2009;42:1691–1699. doi: 10.1021/ar900099h. [PubMed] [Cross Ref]
  • Lindahl E, Hess B, van der Spoel D. GROMACS 3.0: a package for molecular simulation and trajectory analysis. J of Mol Model. 2001;7:306–317.
  • Wang J, Wolf R, Caldwell J, Kollman P, Case D. Development and testing of a general amber force field. J Comput Chem. 2004;25:1157. doi: 10.1002/jcc.20035. [PubMed] [Cross Ref]
  • Singh UC, Kollman PA. An approach to computing electrostatic charges for molecules. J Comp Chem. 1983;5:129–145.
  • Stewart JJP. Optimization of parameters for semiempirical methods V: modification of NDDO approximations and application to 70 elements. J Mol Modeling. 2007;13:1173–1213. doi: 10.1007/s00894-007-0233-4. [PMC free article] [PubMed] [Cross Ref]
  • Sousa da Silva A, Vranken W, Laue E. ACPYPE - Antechamber python parser interface. [] [PMC free article] [PubMed]
  • Ramalho TC, França TC, Cortopassi WA, Gonçalves AS, da Silva AW, da Cunha EF. Topology and dynamics of the interaction between 5-nitroimidazole radiosensitizers and duplex DNA studied by a combination of docking, molecular dynamic simulations and NMR spectroscopy. J Mol Struc. 2011;992:65–71. doi: 10.1016/j.molstruc.2011.02.042. [Cross Ref]
  • Punkvang A, Saparpakorn P, Hannongbuam S, Wolschann P, Beyer A, Pungpo P. Investigating the structural basis of arylamides to improve potency against M. tuberculosis strain through molecular dynamics simulations. Europ J Med Chem. 2010;45:5585–5593. doi: 10.1016/j.ejmech.2010.09.008. [PubMed] [Cross Ref]
  • Balajee R, Rajan MSD. Molecular docking and simulation studies of farnesyl trasnferase with the potential inhibitor theflavin. J Appl Pharm Sci. 2011;8:141–148.
  • ACPYPE Wiki: Testing ACPYPE amb2gmx function. []
  • Campbell RB, Robertson MJ, Trotter J. The crystal and molecular structure of pentacene. Acta Cryst. 1961;14:705. doi: 10.1107/S0365110X61002163. [Cross Ref]
  • Endres RG, Fong CY, Yang LH, Witte G, Wöll C. Structural and electronic properties of pentacene molecule and molecular pentacene solid. Comp Mat Sci. 2004;29:362–370. doi: 10.1016/j.commatsci.2003.09.006. [Cross Ref]
  • Berendsen HJC, Postma JPM, van Gunsteren W F, DiNola A, Haak JR. Molecular dynamics with coupling to an external bath. J Chem Phys. 1984;81:3684. doi: 10.1063/1.448118. [Cross Ref]
  • Wang Y, Chengb H, Wanga Y, Hub T, Hob J, Leeb C, Leia T, Yeha C. Influence of measuring environment on the electrical characteristics of pentacene-based thin film transistors. Thin Solid Films. 2004;467:215–219. doi: 10.1016/j.tsf.2004.04.001. [Cross Ref]

Articles from Nanoscale Research Letters are provided here courtesy of Springer