Osman Barış
Malcıoğlu‡
*ab,
Irene
Bechis
a and
Michel
Bockstedte§
*ab
aChemistry and Physics of Materials, University of Salzburg, Jakob-Haringer-Str. 2a, A-5020 Salzburg, Austria. E-mail: mbaris@metu.edu.tr
bSolid state theory, Friedrich-Alexander-University Erlangen-Nürnberg, Staudtstr. 7B2, 91058 Erlangen, Germany. E-mail: michel.bockstedte@fau.de
First published on 21st January 2020
Thin porphyrin films as employed in modern optical devices or photovoltaic applications show deviating electronic and optical properties from the gasphase species. Any understanding of the physical origin may pave way to a specific engineering of these properties via ligand or substituent control. Here we investigate the impact of crystallization of prototypical porphyrins on the electronic levels and optical properties in the framework of density functional theory and many-body perturbation theory. Crystallization substantially shrinks the HOMO–LUMO gap based on polarization effects. We find a shift of the HOMO to higher energy is consistent with recent experiment of MgTPP multilayer film on Ag (100) [A. Classen et al., Phys. Rev. B, 2017, 95, 115414]. Calculated excitation spectra demonstrate a significant redshift of excitation bands except for the Q bands. These lowest excitation bands, in stark contrast to the strong HOMO–LUMO gap renormalization, remain essentially the same as in the gas phase. Our work underlines the possibility of band-gap engineering via ligand-controlled modification of the polarizability.
Archetypal free-base tetraphenylporphyrin (H2TPP) and Magnesium meso-tetraphenylporphyrin (MgTPP) represent naturally occurring porphyrins frequently encountered in biology. A recent experiment reveals that when H2TPP and MgTPP form thick films on Ag(100) substrate, the quasiparticle band gap shrinks compared to gas phase.18,19 These experiments identify a ∼0.5 eV upshift of the HOMO in the multilayer films with respect to the molecules in the gas phase independent of the film thickness.
In fact, level shifts and a corresponding band-gap renormalization were observed for a wide range molecular crystals in photo emission experiments20 and was attributed to the polarization energy due to the ionization of molecules in the crystal. Recent state of the art calculations using high-level methods such as many body perturbation theory or DFT with optimal tuning of range separated hybrid functionals were able to access these effects,21,22 whereas using semi-local and hybrid exchange–correlation energy functionals failed to capture level-shifts and band gap renormalization. How these effects manifest themselves in the electronic and optical properties of the porphyrin crystals is, however, unclear.
In this article we address the impact of the supramolecular packing on the electronic and optical spectra of archetypical porphyrins H2TPP and MgTPP. We focus on the experimentally reported molecular crystal17 and various other packings of the molecules in the crystalline phase. Besides the fundamental interest in the properties of the porphyrin crystal, the results are also applicable to the modeling of a thick film on a substrate. Indeed, substrate–film interaction effects couple weakly to the bulk part of the film owing to the VdW bonding between the molecules. We demonstrate that the deviation between the gas phase data and the multilayer films in experiment is due to intermolecular interactions. The quasiparticle DOS features for the experimental crystalline packing align well with the photoemission spectra on thick films,19,23 in particular regarding the HOMO position. In the experimental crystalline packing, the vertical transitions of the visible part of the spectrum remains unchanged, whereas in the UV region redshifts by more than 1 eV do occur.
The electronic properties of the isolated molecules and crystals are obtained within the sc-GW027 level of theory using LDA orbitals as a starting point (cf. ESI,† Section S2, for details). This level of theory reproduces NIST calibration data28 for the molecules in the gas phase. In the crystal, the band structure is obtained by interpolating the quasiparticle levels along the k-path trajectory using the maximally localized Wannier functions (MLWFs) method of Marzari and Vanderbilt.29 The details of the parameters and procedures can be found in the ESI,† Section S3.
The optical properties in the visible and UV regions are calculated by solving the Bethe–Salpether equation30 (BSE). Due to the complex nature of the excited state potential energy surface of porphyrins, we restrict the investigation to vertical transitions.
In addition to the P triclinic crystal structure, we investigate three additional possible scenarios for thick MgTPP films following the strong and weak interaction idea from Perepogu and Bangal:32 (1) The simple stack where porphyrin molecules are directly on top of each other. (2) Rotated stack, where porphyrin layers are rotated by 90° with respect to each other (3) Herringbone where each porphyrin is coordinated with the phenyl rings of another porphyrin perpendicular to it. The binding energies of the models are reported in the ESI,† Table S5. The experimentally reported crystal packing, which can be considered as an offset form of the simple stack, is the most favorable energetically.
In the following we focus on the experimentally observed P triclinic crystal structure. The experimental and calculated unit cell parameters are reported in the ESI,† Table S1. The porphine skeleton of H2TPP is significantly distorted in the crystal compared to the isolated molecule and assumes a saddle shape. The point group becomes Ci whereas the gas phase H2TPP has the point group D2h and the porphine macrocycle has D4h. Our calculations identify the strong VdW interaction both between macrocycles and between phenyl ring substituents as the major driving force behind this distortion. In the ESI,† Fig. S1 and Table S2, we report key interatomic distances within the phenyl, and pyrrole sections of the molecule as well as the orientation of the phenyl rings with respect to the macrocycle plane, and the overall macrocycle size. The distorted molecule in the crystal (cf. Table S3, ESI†) shows a significant change in the orientation of the phenyl rings of 10° with respect the isolated molecule and only small changes in bond distances.
In comparison with the experimental structural parameters, the calculated values for the volume and inter molecular distances using the D3 or MDB@rSC dispersion correction variants are slightly smaller: the crystal volume is by 6% lower and the distance between adjacent macro cycles is by 4–5% shorter. On the other hand, the calculated intramolecular bond distance result slightly larger than the corresponding experimental values (cf. ESI† Tables S2). For the latter, the different VdW schemes have little effect and the results follow the common tendency of the PBE exchange–correlation functional. The other structural parameters hinge on the flavor of the VdW scheme (cf. ESI† Tables S1 and S2). Note, however, that the structure optimisation does not account for zero point vibrations, effects of anharmonic vibrations or other thermal effects that are contained in the experimental data.
An interesting aspect is the energy cost for cleaving the layered material (cleavage energy) which we determine based on DFT-VdW calculations. The calculated cleavage energy of the H2TPP and the MgTPP crystal is 0.198 J m−2 and 0.207 J m−2 respectively (cf.Fig. 2). Using the same level of approximation yields a cleavage energy of 0.324 J m−2 for graphite (0.36 J m−2 experimentally). Although the calculated values are smaller than the value for graphite, they are still within the range of other commercially interesting layered materials.
Fig. 2 Calculated cleavage energy of the H2TPP and MgTPP crystals, and graphite as a function of the layer separation (symbols as obtained with DFT-D3, lines guide the eye only). |
Now, we turn our attention to the electronic properties of the P triclinic crystal. The comparison of the gas phase quasiparticle levels with the DOS in crystalline phase (along with relevant experimental values) is presented in Fig. 3. The band structure of the crystal is presented in Fig. 4. All levels are aligned with respect to the vacuum level as a common reference potential. Note that our approach replicates the available NIST data28 for all similar porphyrins in the gas phase as indicated in Fig. 3. The crystalline porphyrin H2TPP is an indirect band-gap semiconductor with a band gap of 2.84 eV, which is more than 1.6 eV smaller compared to the gas phase. The band gap of the crystal at the Γ-point is slightly larger than the indirect gap (+0.3 eV). At the PBE-VdW level we obtain an indirect band gap of the H2TPP crystal (cf. Table S6, ESI†) that is only 0.3 eV smaller than that of the gas phase molecule. This difference is comparable to the dispersion of the highest occupied and lowest unoccupied bands obtained at this level of theory (cf. Table S6, ESI†). The failure to reproduce the band-gap renormalization or level shifts at this DFT-level of theory is not unexpected and is in par with earlier findings.21,33
Fig. 3 Calculated molecular levels and crystal DOS (sc-GW0) aligned to the vacuum level. For the H2TPP molecule, levels are reported both in the gas phase geometry, and the distorted geometry found in the crystal. Gas phase experiments: (A) ref. 36 (B) ref. 37 (C) ref. 38. Multi-layer film on Ag(100) from ref. 18 (D) HOMO, (E) HOMO−1, and (F) LUMO; the reported peak positions are broadened using Gaussians according to the corresponding peak width given in ref. 18. Peak height is adjusted for clarity. |
The HOMO level of the gas phase H2TPP is shifted upwards by ∼0.9 eV and the LUMO level is shifted downwards by ∼0.7 eV in the crystal. Similary the HOMO level of gas phase MgTPP upshifts by ∼0.9 eV and the LUMO level downshifts by ∼0.9 eV. The bands in the crystal are slightly dispersive (for details cf. ESI† Section S2, Table S6 and Fig. S2). This manifests itself as an uneven broadening in the DOS. There are a number of level crossings between LUMO and LUMO+1 bands which are degenerate orbitals in the gas phase. The band structure of the MgTPP crystal is quite similar to that of the H2TPP crystal, with a band gap of 2.95 eV.
Note that a very similar band-gap renormalization is observed in the simple stacking scheme of MgTPP, where the HOMO–LUMO gap reduces to 3.0 eV. The packing density seems to be more relevant to the band-gap renormalization than the exact molecular packing. This suggests that ligand controlled modification of the packing density is a viable strategy for band-gap engineering.
In order to understand the large band-gap renormalization of 36%, we considered the influence of the significant distortion of the molecule in the crystal and intermolecular interaction. A comparison of the energy levels of the isolated H2TPP molecule in its fully relaxed configuration and the distorted geometry in the crystal showed that the distortion into a saddle shape and the rotation of the phenyl rings have negligble effect (cf.Fig. 3). Our finding of a weak dispersion at the GW- and DFT-levels of theory furthermore show that weak interaction between the molecules that causes this band dispersion does not explain the band-gap renormalization or the distinct shifts of the HOMO and LUMO bands. This is in par with similar calculations for other molecular crystals and is explained there by electronic polarization effects.34,35
The calculated HOMO band energy position of both H2TPP and MgTPP crystal fit with photoemission experiments on thick films on Ag(100)18 (cf.Fig. 3). In these experiments the LUMO states were investigated using two-photon photoemission spectroscopy (2PPE).18 As shown in Fig. 3 the broad LUMO peaks are centered at somewhat higher energies than the calculated values, however, the discrepancy is less than the peak width. The LUMO peaks are broadened most likely due to unresolved vibrational modes given their Gaussian shape.18 Besides the high level of vibrational excitation in the LUMO peak, the unknown film structure may also contribute to the deviation. Nevertheless, a reduction of the HOMO–LUMO gap of 1.0 eV and 0.5 eV is also reported18 for the thick MgTPP and H2TPP films, respectively, which although being smaller, confirms the trend found in our calculations for the pure crystalline structures.
Now we focus on the optical properties. The porphyrins share similar features in their optical spectra, namely two peaks in the visible region (Q bands), a very intense peak (B or Soret band), a shoulder on the B band (N band), and two other small peaks, the L and M bands. Q bands are often considered of critical importance for energy harvesting and lighting applications. The origin of these typical features are well captured by the Gouterman model16 which involves the two uppermost occupied orbitals and the two lowest unoccupied orbitals of the porphine skeleton. The effect of the substituents on these features are limited and often indirect.
We calculate the vertical optical transitions using BSE presented in Fig. 5. We observe that, the important Q band vertical transitions are not effected by the formation of the crystal, apart from an additional small feature at 2.5 eV. On the contrary, the effect of crystallization is clearly observed in the Soret band. It redshifts from UV by 0.5 eV and moves into the visible range. Furthermore the reduced symmetry and anisotropic screening in the layered crystal lead, in the case of MgTPP, to an anisotropic splitting of the Soret band. It is interesting to note that the shift in the Soret band aligns with the purple appearance of crystalline H2TPP.
Fig. 5 Imaginary part of the dielectric response, obtained from Bethe–Salpether equation: (a) H2TPP gas phase and crystal and (b) MgTPP gas phase and crystal. |
In recent experiments the optical properties of H2TPP and MgTPP films were investigated.39,40 Note that a direct comparison of the optical spectra presented in this and the experimental works requires incorporation of the vibronic effects, since vibronic enhancement of orders of magnitude in Q bands is well known in porphyrins. Furthermore, the excited states are dynamic Jahn–Teller active.41 Given the flexibility of the macrocycle and the low impact of the distortion on the quasiparticle levels, we expect nevertheless that the trends we describe are robust. Indeed, also in these experiments39,40 only a very small shift in the Q bands (H2TPP films:39 <70 meV) and a much larger shift in the Soret band (MgTPP films:40 ∼200 meV) were observed in the mostly amorphous films. Although the effect in the Soret band of the amorphous MgTPP film is by a factor of two smaller than predicted for the crystal, most likely due to the more isotropic screening in the amorphous film, the experiments match our findings.
The impact of crystallization on electronic levels and optical transitions are qualitatively different. Despite of the large band-gap renormalization, optical Q band vertical transitions remain at the same energy with an additional low intensity transition. The Soret band, however, redshifts towards blue in the crystal, which makes the purple/blue spectral region accessible in solar applications.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9cp06040e |
‡ Present address: METU Orta Dogu Teknik Universitesi, Fizik Bolumu, Ankara, Turkey. |
§ Present address: Institute for Theoretical Physics, Johannes Keppler University, Linz, Austria. |
This journal is © the Owner Societies 2020 |