Skip to main content
  • Research Article
  • Open access
  • Published:

Structure-dependent optoelectronic properties of perylene, di-indenoperylene (DIP) isolated molecule and DIP molecular crystal


Theoretical simulations were designed by first principles approach of density functional theory to investigate the structural and optoelectronic properties of different structural classes of perylene; isolated perylene, diindeno[1,2,3-cd:1′,2′,3′-lm]perylene (DIP) molecule and DIP molecular crystal. The presence of molecular interactions in DIP crystal proved its structure-dependent behaviours. The herringbone molecular arrangement of DIP crystal has influenced the electronic properties by triggering the intermolecular interactions that reduced the energy gaps between HOMO and LUMO of the crystal. Strong hybridization resulting from dense charges population near zero Fermi energy has pushed valence band maxima in the density of states of all perylene structures to higher energies. Under small energy input, charges are transferred continuously as observed in the spectra of conductivity and dielectric. The existence of strong absorption intensities are consistent with the former works and supported by the obtained polarized reflectivity and loss spectra.


Decades ago, significant research efforts were diverted from inorganic to organic electronic materials due to their interesting properties and promising applications in multi-functional organic semiconductor devices such as solar cells [1,2,3,4], sensors [5], nanolasers [6], field-effect transistors (OFETs) [7, 8] and organic light-emitting diode (OLEDs) [9, 10]. Indirectly, these materials have demonstrated a strong stability, high carrier mobility and better efficiency in terms of lifetime, manufacturing cost as well as the outcome performance of organic electronic materials. Organic electronic materials might be divided into two major classes: single molecules and molecular crystals that were studied thoroughly to tune into desirable properties before they were prepared in several forms as nanosheets [11,12,13], nanorod [3, 4], nanotube [12, 13] and etcetera to fulfil the demands of the global community.

Perylene is one of the small model compounds that have been widely used in organic photoelectric and photoconductive fields [11]. The potential implementations using perylene derivatives are the recognition of the practical use of perylene, which comparable to inorganic compounds [29, 32]. Its polycyclic hydrocarbon constructed with the conjunction of single and double bonds proposed the presence of high charges population. The dense charges population have attributed to the presence of weak π–π interaction and/or Van der Waals forces of perylene molecular crystal structure. Both of these weak interactions have a huge responsibility in determining the system properties. However they are difficult to compute either in the gas phase or in periodic solid structures experimentally, thus their existence in perylene molecular crystal is easier to be carried out through theoretical work. The non-covalently bonded perylene molecular crystal was found to be crystallized into two crystal phases: α form and β form. The perylene crystal structure was first determined by Donaldson et al. [14]. However, the accuracy of the two bonds connecting the two naphthalene nuclei was not well resolved. A couple of years later, exact values of the two connector bonds were finally found to be longer (1.50 Å) from the actual value of the expected molecular symmetry [15]. Since then, researchers earnestly investigate the proposed β-perylene crystal than α-perylene crystal, as the results, numerous inventions and findings was breakthrough up to twentieth century.

In 2007, Heinrich et al. [16] have proposed a new perylene crystal structure of diindeno[1,2,3-cd:1′,2′,3′-lm]perylene also known as DIP that hit the market with its remarkable properties. DIP is a small molecular crystal complex that diverted the researcher’s attention from the former perylene molecular crystal. For instance, DIP crystal structure compromised with a different arrangement of molecules that lead to different behaviour of photons and electrons and thus generating different level of performance efficiency. As reported by Liao et al. [17], most of the organic molecular crystals are structure-dependent which later affect the performance of electronic applications especially. Although it is important, the relationship between the building block of molecules and the formation of the organic molecular crystal is remaining elusive.

In the present study, investigations were made to theoretically interpret the structural, electronic and optical properties of perylene in two phases of the isolated molecule and molecular crystal structure. The relationship between the molecule arrangements in DIP organic crystal from the isolated DIP molecule which originated from isolated perylene molecule is highlighted in parallel to their optoelectronic properties. We employed a theoretical approach based on first principles calculation of density functional theory implemented within DMol3 code and employed different exchange and correlation functionals Exc that governed by two parameterizations of LDA, seven parameterization types of GGA, BLYP, and B3LYP. A basis set of DNP 4.4 together with pseudopotential scheme is used to provide high accuracy in the calculations especially on the fairly large compound as DIP crystal.

Computational methodology

All the calculations were carried out by the developed density functional theory program of DMol3 package [18, 19] that mounted in the Materials Studio package. Expansion of wavefunction is dealt by applying numerical basis approach of double numerical plus polarization (DNP) including a polarization p-function on all hydrogen atoms to grasp all the interactions that may occur in both isolated perylene molecule, DIP molecule and DIP molecular crystal structures. Basis file version 4.4 is applied to deliver slightly improved heats of formation. The accuracy level of DNP 4.4 basis set is comparable to the 6-31G* Gaussian orbital basis set. The optimization of perylene and DIP molecule with convergence tolerance of maximum force 0.004 Ha/Å, maximum displacement 0.005 Å, and energy 0.00002 Ha is used and the structures were considered isolated to avoid outside interactions that may affect the result. The quality of integration accuracy and orbital cutoff were set at fine which approximately ~ 0.1 eV/atom. In order to gain insight into the formation mechanism from the isolated perylene molecule up to the perylene molecule extension of DIP molecule, then perylene DIP crystal, unrestricted spin polarization has been chosen. By applying the pseudopotential [20] approach, the core orbitals were treated as frozen. A thermal smearing of 0.0005–0.001 Ha is used for the phase transition effect and to enhance the speed of convergence process on the orbital occupation. All perylene structures were suppressed until the SCF convergence reach 0.0001 eV/atom, and then the structures were allowed to fully relaxed.

In the beginning, exchange–correlation functionals Exc that cover from local spin density approximations (LSDAs) [21], generalized gradient approximations (GGAs) [22] and hybrid functional of B3LYP [23, 24] were employed to achieve optimum accuracy in the calculations. In LSDA, the Exc is dependent on the electron density at every point of the system that is written as

$$ E_{\text{XC}}^{\text{LDA}} \left[ {\rho \left( \varvec{r} \right)} \right] = \mathop \int \nolimits \rho \left( \varvec{r} \right) \varepsilon_{\text{XC}}^{\text{Hom}} \left( {\rho \left( \varvec{r} \right)} \right){\text{d}}\varvec{\tau} $$

where \( \varepsilon_{\text{XC}}^{\text{Hom}} \left( {\rho \left( \varvec{r} \right)} \right) \) is the exchange correlation energy density of an interacting homogeneous electron gas at the electron density \( \rho \left( \varvec{r} \right) \) and \( \varvec{r} \) is the position vector. \( \varepsilon_{\text{XC}}^{\text{Hom}} \left( {\rho \left( \varvec{r} \right)} \right) \) is a function of the local density that can be separated into the function of exchange \( \varepsilon_{X} \left( {\rho \left( \varvec{r} \right)} \right) \) and correlation \( \varepsilon_{C} \left( {\rho \left( \varvec{r} \right)} \right) \) parts.

Whereas in GGA, every point is included by the gradient of electron density \( \nabla \rho \left( \varvec{r} \right) \) into the Exc functional. The Exc of GGA can be written as

$$ E_{\text{XC}}^{\text{GGA}} \left[ {\rho \left( \varvec{r} \right)} \right] = \mathop \int \nolimits \rho \left( \varvec{r} \right) F_{{\varepsilon_{\text{XC}}^{\text{Hom}} }} \left( {\rho \left( \varvec{r} \right), \;\nabla \rho \left( \varvec{r} \right)} \right){\text{d}}\varvec{\tau} $$

where \( F_{{\varepsilon_{\text{XC}}^{\text{Hom}} }} \) is referred to the functional of generalized gradient approximations.

Indifferent with LDA and GGA functional, the hybrid functional B3LYP is constructed by an approach where the Hartree–Fock exact exchange functional is linearly combined with any number of Exc of density functional from LDA or GGA. The Exc of B3LYP can be expressed as

$$ E_{\text{XC}}^{\text{B3LYP}} =\; E_{X}^{\text{LDA}} +\; a_{0} \left( {E_{X}^{\text{HF}} - E_{X}^{\text{LDA}} } \right) +\; a_{X} \left( {E_{X}^{\text{GGA}} - E_{X}^{\text{LDA}} } \right) +\; E_{C}^{\text{LDA}} + \;a_{c} \left( {E_{C}^{\text{GGA}} - E_{C}^{\text{LDA}} } \right) $$

where a 0 = 0.20, a X  = 0.72, a c  = 0.81 are the functional parameter, \( E_{X}^{\text{HF}} \) is the Hatree–Fock exchange energy, \( E_{X}^{\text{GGA}} \) is the GGA exchange energy, \( E_{C}^{\text{LDA}} \) is the LDA correlation energy and \( E_{C}^{\text{GGA}} \) is the GGA correlation energy. These functionals were included to provide extensive study on the electronic properties and an overview on the energy gaps.

As stated previously, calculations by LDAs, GGAs and B3LYP were only performed in the beginning. Later, only GGA formalism with parameterization of Perdew-Burke-Ernzerhof (PBE) [22] is introduced since it was able to deliver reasonable results. Bond length, angle, and volume are some part of the structural properties that are being discussed. Also, the electronic properties are governed by DOS, HOMO–LUMO, and electronic energy gap results will be detailed. And lastly, conductivity, dielectric function, absorption, reflectivity and loss function are discussed for optical properties part. To provide an overview in terms of the material performance, the possibility of loss occurrence was calculated in percentage by analytical through Eq. 4 below;

$$ {\text{Possibility of loss occurence}} = \frac{\text{loss intensity}}{\text{absorption intensity}} \times 100\% $$

Results and discussion

A molecule with definite structure provides a model to investigate the basic interactions of molecules within a single crystal. In the present study, the influence of definite molecule structure plays an important role. Thus, an isolated perylene molecule has been optimized and graphically was presented in Fig. 1a. The obtained lengths of C–H bonds are ± 1.090 Å, single bonds of C2–C3, C3–C4, C17–C18, C18–C19 are 1.415 Å, single bonds of C1–C6, C5–C10, C11–C16, C15–C20 are 1.404 Å, single bonds of C7–C8, C8–C9, C12–C13, C13–C14 are 1.431 Å, single bonds of C7–C12, C9–C14 is 1.470 Å and double bonds of C1=C2, C4=C5, C16=C17, C19=C20 are 1.379 Å, double bonds of C6=C7, C9=C10, C11=C12, C14=C15 are 1.396 Å and last but not least double bonds of C3=C8, C13=C18 are 1.436 Å. The structural geometry of perylene molecule that involved with the values of bond length and the angle is found to be agreed well with the previous experimental [14, 15, 25] and theoretical works [26, 27].

Fig. 1
figure 1

Graphical view of a isolated perylene molecule, b DIP molecule and c DIP single crystal where grey and white balls represent carbon and hydrogen atoms

In order to investigate the molecular interactions effect on a single DIP crystal, the perylene molecule extension of isolated DIP molecule was optimized. As seen in Fig. 1b, DIP molecule is a planar-perylene based molecule where an indeno molecule was attached on both ends of perylene molecule resulting in the formation of a covalent bond; two carbon atoms are shared between perylene molecule and indeno molecule. This has summed up the structure of DIP molecule that contained 32 carbon atoms and 16 hydrogen atoms. As observed, the obtained bond lengths of DIP molecule are slightly different as compared to perylene molecule. The C–H bonds are the same ± 1.090 Å, while the other bonds as in the perylene molecule are slightly increased except C3–C8 and C13–C18 where the bond lengths are slightly reduced. For the attached indeno molecule of DIP molecule, all of the bond lengths are approximately ± 1.400 Å.

From the optimized DIP molecule, DIP crystal structure that consists of 64 carbon atoms and 32 hydrogen atoms was optimized, refer Table 1. With cell dimensions of a = 7.1709 Å, b = 8.5496 Å, c = 16.7981 Å and an angle of β = 92.416°, DIP crystallizes in the primitive-centered monoclinic lattice. The space group of P21/C with a unit cell volume of 1028.95. These values are found to agree well with the previously reported works [16]. In Fig. 1c we displayed the graphical view of the optimized DIP structure. As observed, partially four molecules that were arranged in the ‘herringbone’ molecular stacking structure were taken as the building block of DIP crystal. These molecules arrangement are expected to be the major influence in finding the electronic and optical properties later.

Table 1 Structural properties of the optimized DIP molecular crystal

As known, usual research questions are always prompted by the need for better system efficiency which commonly involved with carrier mobility, charge transport, opening window, and energy gap in order to achieve excellent practical applications. Despite, the solutions are basically relied on understanding the fundamental of its electronic behaviours. Hence to investigate the transition of molecule arrangement from an isolated perylene molecule to a single crystal DIP, the spin-polarized electronic structures were determined in the present study.

For electronic structure, various Exc were taken into account to provide an exclusive comparison in terms of different adopted methods. Tabulated energy gap values of both perylene and DIP molecule and DIP crystal were presented in Tables 2, 3 and 4, and several of the energy gaps were illustrated in graph forms respective to their HOMO and LUMO levels as shown in Fig. 2. A pair of HOMO–LUMO located at the outermost boundaries of the electrons in the molecule is strongly interacting with each other since they lie the closest in the energy of any pair of orbitals in the molecular structure. One may relate the differences in energy between HOMO and LUMO gives the energy gap, Eg = ELUMO − EHOMO. As shown in Fig. 2, similar pattern from all structures was observed where HOMOs–LUMOs are lying at highest energy level when LDA functional was adopted, then following by GGA functional and lastly BLYP functional. Although HOMO–LUMO by LDA functional is lying at the highest energy level, the energy gap obtained is relatively low due to the discrepancy of LDA functional in underestimating the gap values.

Table 2 The energy gaps of isolated perylene molecule using various Exc functional
Table 3 The energy gaps of isolated DIP molecule using various Exc functional
Table 4 The energy gaps of DIP molecular crystal using various Exc functional
Fig. 2
figure 2

Graph of energy gaps and HOMO–LUMO of a perylene molecule, b DIP molecule and c DIP crystal

By referring to Fig. 2, HOMO and LUMO of perylene and DIP molecule are evenly distributed with positive and negative charges. For HOMOperylene_mol, a cloud of charges majorly concentrated on the double bonds, while charges on LUMOperylene_mol are more focused on the single bonds. This supports the existence of localized charges in HOMOperylene_mol contributed from the formation of π bonding. Similar isosurface pattern is also observed on HOMODIP_mol for the original part of perylene molecule. The presence of indeno part of DIP molecule does not affect much on the HOMODIP_mol as compared to LUMODIP_mol. However, by referring Table 3, the HOMODIP_mol clearly lied on higher energy than HOMOperylene_mol.

In LUMOperylene_mol, the cloud charges are smaller in size than in HOMOperylene_mol and mostly are formed on the single bonds. Whereas the presence of indeno part of DIP molecule is more stressed than the type of bonding since it clearly elevates the population of charges on LUMODIP_mol. As observed on the shared carbon atoms of DIP molecule, a larger cloud of charges is formed on LUMODIP_mol. By comparison, less concentration of charges is carried by the single bond than the double bond. However, the compatibility of HOMO–LUMOperylene_mol and HOMO–LUMODIP_mol is answered with the reasonable energy gaps produced. The π* antibonding of both of the LUMO molecules proved its capability in attracting the charges from the HOMO molecules side well. By overcoming small energy of ~ 1.909 and ~ 1.490 eV (mean energy gap values of perylene and DIP molecule, exclude BLYP and B3LYP), the charges from HOMO are transferred to LUMO for perylene [28] and DIP molecule respectively.

On the other hand, the energy gaps of DIP crystal are clearly producing much smaller value compared to the molecules (refer to Tables 2, 3, 4). Typical discrepancy involving LDA [21] functional can be concluded in this study since the gap value has been underestimated either in molecules or crystal form of perylene. The existence of four herringbone DIP molecules in DIP crystal has initiated small intermolecular interaction that reduced the energy gap by ~ 0.092 eV, where the mean energy gap values of DIP is ~ 1.398 eV (exclude BLYP and B3LYP). As displayed in Fig. 2, the larger lobes formation of charges cloud on HOMODIP_crystal that were attracted towards each other also was partially influenced by the molecule arrangement. Indifferently, with the LUMODIP_crystal, the free charges were seen lingering among the lobes and they were able to move freely. Hence the possibility for transferring the charges from HOMODIP_crystal to LUMODIP_crystal is increased because the free charges will attract the opposite charges on the HOMODIP_crystal side stronger.

Since DIP is a crystal, the band-gap calculation is performed within GGA–PBE Exc functional and the electronic bandstructure was presented in Fig. 3. No change in momentum has been observed. It reveals that the energy required for an electron to jump from conduction to valence band is relatively small. It was proved that perylene is one of the donors known due to its highly-conducting organic crystal characteristic [29]. The minimum valence band and maximum conduction band were shown to be direct in nature along the Γ symmetry point as displayed in Fig. 3. The conduction band was fully occupied up to the zero Fermi energy level, while the minimum valence band is lying at 0.046 Ha = 1.252 eV.

Fig. 3
figure 3

DIP crystal band structures along the high-symmetry points in Brillouin zone

From the density of states (DOS) of all perylene structures, the population of charges is becoming denser in the order from the less charges tabulation of isolated perylene molecule followed by the isolated DIP molecule and the most highly tabulated is DIP crystal (as seen from Fig. 4). This suggests that the existence of indeno molecule, free charges and the weak intermolecular forces within DIP crystal clearly affected the charges population compared to perylene and DIP molecules. As observed in all DOSs, the charges population are highly concentrated in the valence region, which indicates that all structures possess high tendency to behave as n-type semiconductor materials. Besides, the presence of strong hybridization has elevated the charges population near the zero Fermi energy level and consequently, the valence band maxima of the molecules and DIP crystal are pushed to higher energies. Since charges population are dense near zero Fermi energy, the possibility to succeed in transmitting the charges from valence region to conduction region is increased as well. As we compared, the gaps produced in DOS are in parallel with the previously found on the HOMO–LUMO energy gaps and band gap values for all the structures.

Fig. 4
figure 4

Computed density of states (DOS) of perylene molecule, DIP molecule and DIP molecular crystal

Besides electronic structure, the behaviour regarding the optical characteristics is also one of the most interesting topics to dig in for an optoelectronic material. The contribution of efficient optical properties, to capture and transmit the falling photons and turned the respective energy into electricity, is the key to the dream of optoelectronic practical applications. In this study, the frequency-dependent optical properties that covered with conductivity, dielectric function, absorption, reflectivity and loss function of perylene molecule and DIP crystal were calculated as displayed in Figs. 5 and 6. As mentioned earlier, LDA functional consistently give lower energy gap, whilst BLYP and B3LYP approaches are known to be expensive. Thereby experience GGA-PBE functional was trusted to perform these calculations [22]. In addition, standard GGA functionals were found to yield reasonable results [30].

Fig. 5
figure 5

Real and imaginary parts of conductivity and dielectric functions spectra of perylene molecule (1st row), DIP molecule (2nd row) and DIP crystal (3rd row)

Fig. 6
figure 6

a Absorption, b reflectivity and c loss function spectra of perylene molecule (red spectra), DIP molecule (blue spectra) and DIP crystal (black spectra)

Figure 5 shows the conductivity and dielectric function curves for both molecules and crystal of perylene. As seen, values in all of the properties are higher in DIP crystal structure than in perylene molecules. The first peak in the real conductivity spectra was led by DIP crystal at 1.6450 eV with a value of 0.8625, compared to perylene and DIP molecules where it was found at 3.6875 eV with a value of 0.0419 and 4.6875 eV with a value of 0.3499 respectively. Stronger optical intensity produced by DIP crystal in conductivity spectrum is continued up to 17.500 eV where the last peak with the intensity of 0.675 is observed at 13.6575 eV, whereas no visible peak is observed in the energy range 10–17.5 eV for perylene molecule and only a small peak is observed at 13.8541 eV with intensity 0.0807 for DIP molecule. Herringbone molecule arrangement in DIP crystal proved to carry extra strength in the conducting ability since charges were able to transfer under small energy ≥ 1.645 eV contributes from the interaction between molecules.

Dielectric coefficients are closely related to the performance of conductivity; where better dielectric properties will give a hand in conducting the electrons. Thus, to support the findings in conductivity, calculations of dielectric functions were been executed as well, and the results are displayed in Fig. 5. As seen, the first peak of perylene molecule with εreal(perylene_mol) = 1.0803 at 2.9583 eV was slightly underperformed compared to the same class of its molecule structure. While DIP molecule performing better within a moderate range where the first peak is found at 2.5417 eV with εreal(DIP_mol) = 1.5732. This totally opposite on DIP crystal where its properties were found to perform significantly excellent even at the extremely low presence of energy, 0 eV. DIP crystal conquered in transferring the electrons from HOMO to LUMO with high intensity of εreal(DIP_crystal); 0 eV = 5.4556 at 0 eV and εreal(DIP_crystal); 0.79 eV = 6.1111 at 0.7917 eV which represents the highest peak. By comparison, no plateau is seen on DIP crystal and molecule dielectric spectrum which contradicts to perylene molecule where it has reached a plateau at 6.548 eV.

Experimental measurements of molecular crystal absorption usually encounter some difficulties because the observed absorption spectrum frequently suffers from reflection on the crystal surfaces. However, for the same case, theoretical approaches could provide reliable results without the presence of external disturbances. As displayed in Fig. 6a, the absorption spectra k of DIP molecular crystal and both perylene molecules were successfully obtained from the theoretical approach. The absorption spectrum of DIP crystal shows an almost similar pattern with a previous study on the former β-perylene crystal [31] and DIP molecule [27] structure. By comparison, DIP crystal possess higher intensity than the obtained perylene molecule, DIP molecule, and the previous study. Highest peak was observed at 4.972 eV with an intensity of k DIP_crystal; 4.97 eV = 85,563, while in perylene and DIP molecule the highest peak intensity are k perylene_mol; 3.72 eV = 2738 and k DIP_mol; 4.94eV = 23,901 respectively. This suggesting there is no presence of reflectivity due to the disturbances in the calculations.

To validate these findings, reflectivity calculations were performed as well. Figure 6b depicts reflectivity spectra of all three perylene structures. All structures present small reflectivity intensities, especially in perylene molecule spectrum. Since reflectivity intensities are small, directly absorption ability is enhanced since the photon falls are transmitted into the material easily with the reflection that rarely happened. Obviously, two visible peaks are observed in all structures, where the highest peak intensity (x, y) is (3.3889, 0.0006295) for perylene molecule, (4.9722, 0.02325) for DIP molecule and (5.8889, 0.4211) for DIP crystal. All reflectivity spectra have a background in the energy region from 0 to 10 eV, the intensity of which increases gradually with increasing energies, then gradually decreases after the highest peaks are reached. This indicates the presence of continuous reflectivity corresponds to the absorption spectra calculation; however with a very weak/small reflectivity intensities, it will not affect the obtained absorption spectra.

By considering the presence of reflectivity, loss function L calculations were executed in order to count probability of loss occurrence. Similar spectra of loss function that take place near to the highest peaks of reflectivity spectra are observed in all perylene structures have partially supported the existence of reflectivity (as seen in Fig. 6c). Moreover, the loss occurred in DIP crystal that is higher than perylene and DIP molecules is in line with the finding of higher reflectivity intensity in DIP crystal. Although there were loss occurrences in the structures, the maximum loss intensities were considered small since maximum L perylene_mol; 3.718 eV = 0.0885, L DIP_mol; 5.395 eV = 0.6291 and L DIP_crystal; 5.379 eV = 5.4026. In order to illustrate the possibility of loss occurrence, Eq. 4 was analytically calculated upon energy that produced maximum loss intensities, and the absorption intensity was extracted on the same energy. Hence, the maximum possibility of loss to occur respective to the selective energy was only 0.003210% for perylene molecule, 0.002919% for DIP molecule and 0.006693% for DIP crystal. Thus, the absorption disturbances which were majorly contributed from reflectivity can be neglected since the loss percentages were extremely small.


First principles of DFT approach calculations have been successfully executed with various exchange–correlation functionals (E xc ) in order to investigate the influence of molecules arrangement of the perylene-based molecule in DIP molecular crystal. In addition, comparative study of optoelectronic properties between those perylene classes of isolated perylene molecule, isolated DIP molecule and DIP molecular crystal were performed as well. Our investigations have revealed the effect of indeno molecule and the presence of molecular interactions from the herringbone molecular arrangement of DIP crystal. The indeno molecule has elevated the population of charges on both DIP molecule and crystal. The molecular interactions of DIP molecules within DIP crystal have affected the obtained energy gap to be reduced since the HOMODIP_crystal and the LUMODIP_crystal were attracted stronger to each other. The small energy gap resulted from the molecular attraction clearly contributes to the intense spectra of conductivity and dielectric function. Directly, this has increased the attracting and conducting abilities of DIP crystal. Other optical parameters such as absorption, reflectivity, and loss function were evaluated over same energy range as conductivity and dielectric function and our results are agreed well with the previous studies. Besides, the maximum loss that may occur is extremely low with 0.003210% for perylene molecule, 0.002919% for DIP molecule and 0.006693% for DIP crystal.


  1. Li G, Zhu R, Yang Y (2012) High-performance multiple-donor bulk heterojunction solar cells. Nat Photonics 6(3):153–161

    Article  CAS  Google Scholar 

  2. Raj MR, Ramkumar S, Anandan S (2013) Photovoltaic studies on perylene diimide-based copolymers containing electronic push–pull chromophores. RSC Adv 2013(3):5108–5120

    Article  Google Scholar 

  3. Huynh WU, Dittmer JJ, Alivisatos AP (2002) Hybrid nanorod-polymer solar cells. Science 295(5564):2425–2427

    Article  CAS  Google Scholar 

  4. Menke SM, Holmes RJ (2014) Exciton diffusion in organic photovoltaic cells. Energy Environ Sci 2014(7):499–512

    Article  Google Scholar 

  5. Harsányi G (2000) Polymer films in sensor applications: a review of present uses and future possibilities. Sens Rev 20(2):98–105

    Article  Google Scholar 

  6. Wang H, Xie Z, Ma Y, Shen J (2007) Progress on the optoelectronic functional organic crystals. Sci China Ser B 50(4):433–452

    Article  CAS  Google Scholar 

  7. Gao Y, Ding H, Wang H, Yan D (2007) Electronic structure of interfaces between copper-hexadecafluoro-phthalocyanine and 2,5-bis(4-biphenylyl) bithiophene. Appl Phys Lett 91(14):142112

    Article  Google Scholar 

  8. Beaujuge PM, Fréchet JM (2011) Molecular design and ordering effects in π-functional materials for transistor and solar cell applications. J Am Chem Soc 133(50):20009–20029

    Article  CAS  Google Scholar 

  9. Kulkarni AP, Tonzola CJ, Babel A, Jenekhe SA (2004) Electron transport materials for organic light-emitting diodes. Chem Mater 16(23):4556–4573

    Article  CAS  Google Scholar 

  10. Burroughes J, Bradley D, Brown A, Marks R, Mackay K, Friend R, Burns P, Holmes A (1990) Light-emitting diodes based on conjugated polymers. Nature 347(6293):539–541

    Article  CAS  Google Scholar 

  11. Aldahhak H, Matencio S, Barrena E, Ocal C, Schmidta WG, Rauls E (2015) Structure formation in diindenoperylene thin films on copper(111). Phys Chem Chem Phys 17:8776–8783

    Article  CAS  Google Scholar 

  12. Lai Y, Li H, Pan J, Guo J, Kang L, Cao Z (2015) Synthesis of ultrathin nanosheets of perylene. Cryst Growth Des 15(3):1011–1016

    Article  CAS  Google Scholar 

  13. Zhang X, Zhang X, Zou K, Lee C-S, Lee S-T (2007) Single-crystal nanoribbons, nanotubes, and nanowires from intramolecular charge-transfer organic molecules. J Am Chem Soc 129(12):3527–3532

    Article  CAS  Google Scholar 

  14. Donaldson D, Robertson J, White J (1953) The crystal and molecular structure of perylene. In: Proceedings of the Royal Society of London A: mathematical, physical and engineering sciences (unpublished)

  15. Camerman AT, Trotter J (1964) The crystal and molecular structure of perylene. In: Proceedings of the Royal Society of London A: mathematical, physical and engineering sciences (unpublished)

  16. Heinrich MA, Pflaum J, Tripathi AK, Frey W, Steigerwald ML, Siegrist T (2007) Enantiotropic polymorphism in di-indenoperylene. J Phys Chem C 111(51):18878–18881

    Article  CAS  Google Scholar 

  17. Liao Q, Zhang H, Zhu W, Hu K, Fu H (2014) Perylene crystals: tuning optoelectronic properties by dimensional-controlled synthesis. J Mater Chem C 45:9695–9700

    Article  Google Scholar 

  18. Delley B (1990) An all-electron numerical method for solving the local density functional for polyatomic molecules. J Chem Phys 92(1):508–517

    Article  CAS  Google Scholar 

  19. Delley B (2000) From molecules to solids with the DMol3 approach. J Chem Phys 113(18):7756–7764

    Article  CAS  Google Scholar 

  20. Delley B (2002) Hardness conserving semilocal pseudopotentials. Phys Rev B 66(15):155125

    Article  Google Scholar 

  21. Vosko S, Wilk L, Nusair M (1980) Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis. Can J Phys 58(8):1200–1211

    Article  CAS  Google Scholar 

  22. Perdew JP, Burke K, Ernzerhof M (1996) Generalized gradient approximation made simple. Phys Rev Lett 77(18):3865

    Article  CAS  Google Scholar 

  23. Lee C, Yang W, Parr RG (1988) Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys Rev B 37(2):785–789

    Article  CAS  Google Scholar 

  24. Becke AD (1993) Density-functional thermochemistry. III. The role of exact exchange. J Chem Phys 98(7):5648–5652

    Article  CAS  Google Scholar 

  25. Goodwin TH (1960) Calculated bond lengths in some cyclic compounds. Part VII. The series naphthalene, perylene, terrylene, quaterrylene, and the lengths of low-order bonds. J Chem Soc (Resumed) 4851–4856

  26. Fedorov I, Zhuravlev Y, Berveno V (2013) Structural and electronic properties of perylene from first principles calculations. J Chem Phys 138(9):094509

    Article  CAS  Google Scholar 

  27. Ramaniah LM, Boero M (2006) Structural, electronic, and optical properties of the diindenoperylene molecule from first-principles density-functional theory. Phys Rev A 74(4):042505

    Article  Google Scholar 

  28. Ferraris J, Cowan D (1973) Electron transfer in a new highly conducting donor–acceptor complex. J Am Chem Soc 95:948

    Article  CAS  Google Scholar 

  29. Kimball D, Munns R, Fisher SP, Bartolo M, Valdez J, Teat SJ, Tsui S, Reinheimer EW (2015) Utilizing perylene in new organic donor–acceptor materials: highlighting the synthesis, structure and physical properties of perylene-pDNB and perylene-TCNP. J Chem Crystallogr 45(4):169–177

    Article  CAS  Google Scholar 

  30. Liu W, Tkatchenko A, Scheffler M (2014) Modeling adsorption and reactions of organic molecules at metal surfaces. Acc Chem Res 47(11):3369–3377

    Article  CAS  Google Scholar 

  31. Fuke K, Kaya K, Kajiwara T, Nagakura S (1976) The polarized reflection and absorption spectra of perylene crystals in monomeric and dimeric forms. J Mol Spectrosc 63(1):98–107

    Article  CAS  Google Scholar 

  32. Roy UN, Cui Y, Hawrami R, Burger A, Orona L, Goldstein JT (2006) AgGaSe 2: a highly photoconductive material. Solid State Commun 139(10):527–530

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

The main idea for the work was thought up by SGS and RA. Theoretical calculations were performed by MM. The analysis of the theoretical data was carried out by MM, RA and SGS. All authors have participated in writing the manuscript. All authors read and approved the final manuscript.


Authors from UTM acknowledge Ministry of Education (MOE)–Malaysia and Universiti Teknologi Malaysia (UTM) for financial support/facilities for this research through Grant Numbers; Q.J130000.2526.04H14 and Q.J130000.2526.06H15 and R.J130000.7826.4F508. Financial support to SGS for this study is acknowledged in the form of an internal Grant IRG 16417 from office of research at Alfaisal University.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Souraya Goumri-Said.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mohamad, M., Ahmed, R., Shaari, A. et al. Structure-dependent optoelectronic properties of perylene, di-indenoperylene (DIP) isolated molecule and DIP molecular crystal. Chemistry Central Journal 11, 125 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: