Carmen
Coppola
ab,
Adriana
Pecoraro
c,
Ana B.
Muñoz-García
c,
Rossella
Infantino
de,
Alessio
Dessì
e,
Gianna
Reginato
e,
Riccardo
Basosi
abe,
Adalgisa
Sinicropi
*abe and
Michele
Pavone
*f
aDipartimento di Biotecnologie, Chimica e Farmacia, R2ES Lab, Università degli Studi di Siena, via Aldo Moro 2, 53100, Siena, Italy. E-mail: adalgisa.sinicropi@unisi.it
bCSGI, Consorzio per lo Sviluppo dei Sistemi a Grande Interfase, via della Lastruccia 3, 50019, Sesto Fiorentino, Italy
cDipartimento di Fisica “Ettore Pancini”, Università degli Studi di Napoli Federico II, Comp. Univ. Monte Sant’Angelo, Via Cintia 21, 80126, Napoli, Italy
dDipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100, Siena, Italy
eIstituto di Chimica dei Composti OrganoMetallici (ICCOM), Consiglio Nazionale delle Ricerche (CNR), Via Madonna del Piano 10, 50019, Sesto Fiorentino, Italy
fDipartimento di Scienze Chimiche, Università degli Studi di Napoli Federico II, Comp. Univ. Monte Sant’Angelo, Via Cintia 21, 80126, Napoli, Italy. E-mail: michele.pavone@unina.it
First published on 24th May 2022
Recently, great research efforts have been devoted to perovskite solar cells (PSCs) leading to sunlight-to-power conversion efficiencies above 25%. However, several barriers still hinder the full deployment of these devices. Critical issues are related to PCE stability and device lifetimes, which could be improved by targeted engineering of the hole transport material (HTM). Indeed, the HTM is not only responsible for transporting holes and preventing direct contact between the photo-active perovskite and the charge collector layer, but it plays important structural and protective roles too. As alternatives to the widely used yet expensive and unstable Spiro-OMeTAD, organic HTMs based on triphenylamine (TPA) and phenothiazine (PTZ) moieties have been proposed. However, their performances in PSC devices, and in particular their interfacial properties with the most popular methylammonium lead iodide perovskite (MAPI) still need investigations to be fully determined. In this framework, here we report a first-principles study on the structural and the electronic properties of a recently designed TPA and PTZ-based HTM (HTM1) and its interface with the MAPI (001) surface, considering both the PbI2- and the MAI-terminations. We also addressed already known HTM molecular systems to allow for a direct comparison with the recently proposed HTM1: we characterized the molecular parameters and the MAPI/HTM interfacial properties for Spiro-OMeTAD, PTZ1, and PTZ2. Our results suggest that good adhesion properties do not ensure effective and efficient MAPI-HTM hole injection. Despite the theoretical good alignment between HTM1 HOMO and MAPI valence band edge, our results for the mutually polarized interface point out the lack of a sufficient driving force for hole transport. While the hole mobility of HTM1 outperforms those of the other HTM molecules, for this HTM molecule, our findings suggest the application of lead halide perovskite compositions other than MAPI, with substituents that lower its valence band maximum potential value.
To overcome the drawbacks of Spiro-OMeTAD, a wide library of novel doped or dopant-free organic HTMs has been proposed and developed.13,14 In particular, HTMs based on triphenylamines (TPA) and phenothiazine (PTZ) moieties (Scheme S1, ESI†) have attracted great research interest. The TPA moieties feature a non-planar geometry that confers unique hole transport properties, while the PTZ unit has an electron-donating butterfly structure, suitable for several kinds of functionalization, which allows the formation of stable radical cations.14–26 HTMs designed from a combination of TPA and PTZ moieties, e.g.PTZ2 reported by Grisorio et al.,26 have shown efficiencies comparable to those obtained with Spiro-OMeTAD. With this aim, we recently designed in silico four novel TPA and PTZ-based HTMs, HTM1-4, and investigated their electronic properties by means of density functional theory (DFT) and time dependent density functional theory (TDDFT) methods; the results obtained for these novel HTMs were comparable to those obtained for Spiro-OMeTAD and PTZ2.27 Among them, HTM1 showed the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) positions at −4.88 eV and −1.35 eV, respectively, hence higher in energy than the VBM and CBM of MAPI, respectively. Moreover, the wavefunction plot analysis revealed the HOMO of neutral HTM1, as well as the singly occupied molecular orbital (SOMO) of (HTM1)˙+, are fully delocalized over the molecule, while the LUMO is confined on its central core. Such correct energy levels alignment and the slight HOMO LUMO overlap suggest that HTM1 could possess good intermolecular hole-transport ability.25 Moreover, the calculated absorption maximum of neutral HTM1 was found to be 337 nm (3.68 eV), which minimizes competition with the perovskite absorption in the visible range, and it has shown good chemical hardness, higher than that of Spiro-OMeTAD.27 Further synthesis and spectroscopic characterization ratified HTM1 (also named BPT-1) as a valid HTM in PSCs.28
Although the design process of highly efficient PSCs largely benefits from physicochemical studies on isolated HTMs, it is also crucial to investigate their physicochemical interactions with the perovskite, as they strongly influence the ultimate performance of the devices.13 Nevertheless, due to the large size of state-of-the-art HTMs, perovskite/organic HTM interfaces have been rarely studied from the computational point of view, and, thus, HTM adsorption features and interface electronic properties at the microscopic level are yet to be fully understood. To the best of our knowledge, the few computational studies available have explored the interfacial behavior of Spiro-OMeTAD on different MAPI surfaces.29–32 The early work by Torres et al.30 investigated the binding modes of methoxybenzene, i.e.Spiro-OMeTAD anchoring groups, on the MAPI (001) surface to clarify the HTM/perovskite attachment mechanism. They found that interaction mainly occurs when methoxy groups occupy MA+ sites, while a repulsive behavior is observed with the surface iodine atoms.30 Saliba et al.,29 instead, theoretically investigated the whole Spiro-OMeTAD/MAPI interface. They focused on the interaction between the perovskite and the molecule and confirmed that it does not involve the core of the molecule but only the protruding methoxy groups.29 The entire Spiro-OMeTAD/MAPI system has also been studied by Yin and co-workers.31 This work points out the importance of the perovskite exposed facet in promoting the hole transfer. The most recent work devoted to a full Spiro-OMeTAD/MAPI interface by Wang et al.,32 explores the role of interface termination, namely PbI2 and MAI, in determining both the strength of the interaction with the MAPI (110) surface and the ability to ensure the charge transfer across the interface.32 Besides the ones with Spiro-OMeTAD, other interfaces of the MAPI with TPA-based molecules have been theoretically investigated, both in terms of the interaction and driving force for the hole injection.33–36 Motivated by the inspiring preliminary results obtained for isolated HTM1,27 in this work we aim to assess the potential of the TPA and PTZ-based HTM1 molecule as HTM in PSCs involving MAPI as an absorber by first-principles DFT-calculations. We first evaluate key electronic properties of isolated HTM1 (Scheme 1(a)) related to the hole mobility and band alignment with the perovskite, such as the reorganization energies and redox potentials. In order to assess perovskite/HTM mutual polarization effects, we then investigate MAPI/HTM1 interfaces, with HTM1 in contact with either the PbI2- or MAI-terminated MAPI (001) surface. Binding energies are rationalized in terms of structural modifications and bond formation, which are the main features leading to new channels for charge transport. We also assess how the mutual perovskite-HTM1 polarization affects the HOMO-VBM alignment, providing results that contradict the results of isolated molecules, but allow us to understand available experimental outcomes, with direct implications on the final device performances. To allow a direct comparison with state-of-the-art and well-known TPA and PTZ-based HTMs, we extend our interfacial study to Spiro-OMeTAD (Scheme 1(b)), PTZ1 (Scheme 1(c)) and PTZ2 (Scheme 1(d)).26
![]() | ||
Scheme 1 Molecular structure drawings of molecules investigated in this work: (a) HTM1 (b) Spiro-OMeTAD (c) PTZ1 and (d) PTZ2. |
MOLECULE | METHODS | λ (eV) | G ox (eV) |
---|---|---|---|
HTM1 | PBE/6-31G(d,p) | 0.064 | −4.48 |
B3LYP/6-31G(d,p) | 0.092 | −4.99 | |
Spiro-OMeTAD | PBE/6-31G(d,p) | 0.087 | −4.61 |
B3LYP/6-31G(d,p) | 0.132 | −5.00 | |
PTZ1 | PBE/6-31G(d,p) | 0.233 | −4.66 |
B3LYP/6-31G(d,p) | 0.381 | −4.95 | |
PTZ2 | PBE/6-31G(d,p) | 0.139 | −4.71 |
B3LYP/6-31G(d,p) | 0.217 | −5.10 |
The reorganization energy values present the trend HTM1 < Spiro-OMeTAD < PTZ2 < PTZ1, at both levels of theory, with the λ of Spiro-OMeTAD in accordance with the B3LYP results of Chi et al. in dichloromethane (0.131 eV).45 This outcome reports that the TPA and PTZ-based dimeric HTM1 has a lower reorganization energy than those of the monomeric PTZ1 and PTZ2. Additionally, it also has lower reorganization energy than that of Spiro-OMeTAD, hence we predicted good hole mobility and potentially minor need for further additive upon device assembly. The higher reorganization energies of Spiro-OMeTAD, PTZ2 and PTZ1 could be among the reasons behind the need for p-type doping, as reported in the literature.26,46
The low reorganization energy for HTM1 can be associated with a slightly more rigid structure than those of the other molecules. Indeed, moving from its neutral to its oxidized form, we observed only a little average decrease of ∼0.01 Å along main bond lengths (C–C) and ∼3.1° in main dihedral angles (C–C–C–C) among the PTZ and the TPA moieties. Similarly, upon oxidation, the average bond lengths (C–N) and dihedral angles (C–C–N–C) among the 9,9′-spirobifluorene core and the arylamine moieties in Spiro-OMeTAD decreased by ∼0.01 Å and ∼4.6°, respectively. On the other hand, the decrease in the main bond lengths and dihedral angles among the PTZ and the TPA moieties is ∼0.2 Å (C–C) and ∼7° (C–C–C–C) moving from neutral to oxidized PTZ2 and of ∼0.2 Å (C–N) and ∼12.7° (C–C–N–C) moving from neutral to oxidized PTZ1. Hence, by a difference of ca. 3.9°, in PTZ2 the torsional freedom of the C–C–C–C dihedral angles among the PTZ and the TPA moieties is higher than that observed for HTM1, while the C–C–N–C dihedral angles of PTZ1 have a torsional freedom higher than that of PTZ2 of ca 5.7°. These increased values correspond to higher reorganization energy values with respect to those of HTM1.
Oxidation potentials (Gox) have been calculated as the difference between the free energies of the neutral and the oxidized states of each molecule (G0 and G+, respectively), i.e. Gox = G0 − G+, and are listed in Table 1. The computed Gox values follow the expected trend for both DFT functionals. Indeed, a recent article by Wang et al.47 reports oxidation potentials of 50 compounds, theoretically determined employing different exchange correlation functionals and compared with experimental values: pure functionals as PBE underestimate Gox, while the presence of Hartree–Fock-like exchange in hybrid functionals increases the Gox values.47 In particular, the calculated Gox for Spiro-OMeTAD with B3LYP/6-31G(d,p) in this work (−5.00 eV) is in good agreement with the one reported by Wang et al. at B3LYP/6-31G* (−4.98 eV) and with Spiro-OMeTAD experimental oxidation potential in CH2Cl2 (-5.09 eV).32 Calculated Gox of PTZ1 (-4.95 eV) and PTZ2 (-5.10 eV) are in fair agreement with the experimental oxidation potentials obtained at −4.77 eV and −5.15 eV by Grisorio et al.,26 respectively. All lying in a narrow range, the Gox of HTM1 calculated and reported here for the first time is very similar to those of the other three analysed HTMs and it is higher in energy than the experimental value of the MAPI VBM potential (∼−5.4 eV).9,38
A slab supercell with different periodicity within the xy plane has been considered for each termination and each HTM, according to their size, i.e. 7 × 7 for HTM1, 4 × 4 for Spiro-OMeTAD and PTZ1 and 6 × 6 for PTZ2. In order to account for the molecule's complex structure and to obtain reliable energy minima for each termination, different adsorption modes of HTM1 have been considered. HTM1 has been placed on the MAPI surface either exposing the PTZ or the TPA moieties, which have been labeled as down and up orientations, respectively (Fig. S2, ESI†). The same strategy and labeling have been adopted also for PTZ1 and PTZ2 (Fig. S3 and S4, respectively, ESI†). In the case of Spiro-OMeTAD, it has been previously reported that it interacts with the MAPI surface exposing the methoxybenzene of the arylamine moieties,29–32 thus the up and down orientations refer to the methoxybenzene of the two different facets of the 9,9′-spirobifluorene core. The minimum-energy interfaces from preliminary rigid scans have been allowed to relax without symmetry constraints. After relaxation, HTMs/MAPI:PbI2 and HTMs/MAPI:MAI interface binding energies have been computed as follows:
Ebind = Einterface − EMAPI − Emolecule |
HTM1 preferably adopts the up orientation on both PbI2 and MAI terminations, which causes some distortions in the overall surfaces (Fig. 1(a) and 2(a)) due to the interaction of methoxybenzene of the TPA moieties with the MAPI. In particular, HTM1/MAPI:MAI interaction is stronger than that calculated for HTM1/MAPI:PbI2, due to the formation of some hydrogen bonds between the oxygen atoms of HTM1 methoxybenzene groups and the hydrogen atoms of MAPI methylammonium cations. Differently, when HTM1 adopts the down orientation on both MAPI:PbI2 and MAPI:MAI, the resultant interaction with the MAPI (001) surface is just a localized event, which is initiated by the sulfur atoms of the PTZ moieties (Fig. S2, ESI†).
Despite the molecular similarities between HTM1 and PTZ2, the latter gets stabilized by adopting the down orientation on both MAPI (001) terminations (Fig. 1(b) and 2(b)). In these cases, PTZ2 interacts with MAPI by exposing the sulfur atom of the PTZ moiety and additional interactions occur through the methoxybenzene of the TPA moieties and from the methoxybenzene directly linked to the PTZ nitrogen. These additional interactions are unlikely to occur in HTM1, since the benzene ring linking the two PTZ moieties has limited rotational freedom and the molecule is overall more rigid, which limits the interactions between the TPA moieties and the MAPI in the down orientation (Fig. S6, ESI†). Regarding PTZ1/MAPI, the up orientations are preferred, where the HTM exposes the metoxybenzene of arylamine moieties. More convenient binding energies for PTZ1/MAPI:MAI than for PTZ1/MAPI:PbI2 can be attributed to the formation of hydrogen bonds between the oxygen of methoxybenzene and the hydrogens of methylammonium cations (Fig. 1(c) and 2(c)), as it also occurs for HTM1/MAPI:MAI (up). Indeed, differently from PTZ2, when PTZ1 adopts the down orientation, the sulfur atom of the PTZ moiety is less involved in the interaction with MAPI, because of a more pronounced torsion of the methoxybenzene directly linked to the PTZ nitrogen, which has an almost parallel position during the interaction in PTZ2 instead (Fig. S7, ESI†). As for Spiro-OMeTAD, our calculations are in close agreement with previous results29–32 reporting the down orientation as the preferred one on MAPI surfaces (Fig. 1(d) and 2(d)). Spiro-OMeTAD, indeed, interacts with the MAPI through the methoxybenzene of arylamine moieties. In particular, when Spiro-OMeTAD adopts the down orientation, methoxybenzene has a preferential arrangement with respect to the MAPI surface, which allows a closer interaction between the molecule and the surface, thus contributing to stabilizing the binding energy values. Hence, in the case of Spiro-OMeTAD/MAPI:PbI2 the oxygens of methoxybenzene interact with the Pb atoms of the MAPI (Fig. 1(d)), while in the case of Spiro-OMeTAD/MAPI:MAI the oxygens of methoxybenzene interact with the hydrogens of methylammonium cations (Fig. 2(d)). While in the first case the interaction is stronger than those of HTM1/MAPI:PbI2 and PTZ1/MAPI:PbI2, Spiro-OMeTAD/MAPI:MAI has lower binding energy with respect to HTM1/MAPI:MAI and PTZ1/MAPI:MAI, which are more stabilized thanks to the formation of hydrogen bonds between the molecules and the MAPI surface. Indeed, while in HTM1/MAPI:MAI and PTZ1/MAPI:MAI the distances among the oxygen atoms of the methoxybenzene and some hydrogen atoms of the methylammonium cations recall the typical distances of an intermolecular hydrogen bond (e.g. dO7–H828 = 1.74 Å in HTM1/MAPI, dO1–H275 = 1.73Å in PTZ1/MAPI), major distances are found in Spiro-OMeTAD/MAPI (e.g. dO2–H143 = 2.17 Å).
![]() | ||
Fig. 3 pDOS plots of (a) HTM1/MAPI:PbI2 (up), (b) PTZ2/MAPI:PbI2 (down), (c) PTZ1/MAPI:PbI2 (up) and (d) Spiro-OMeTAD/MAPI:PbI2 (down), computed at GGA-PBE, including D3BJ dispersion forces. |
![]() | ||
Fig. 4 pDOS plots of (a) HTM1/MAPI:MAI (up), (b) PTZ2/MAPI:MAI (down), (c) PTZ1/MAPI:MAI (up) and (d) Spiro-OMeTAD/MAPI:MAI (down) computed at GGA-PBE, including D3BJ dispersion forces. |
Regarding PbI2-terminated surfaces (Fig. 3), the pDOS reveal that the VBM of MAPI slightly lies at a higher energy than the HOMO level of HTM1 and PTZ2 (Fig. 3(a and b)). Thus for the standard MAPI perovskite, there is no adequate driving force for the hole-injection process and undesired charge-recombination events can occur. Differently, in the cases of Spiro-OMeTAD and PTZ1, the pDOS highlight a slight but convenient driving force for the hole injection at the MAPI:PbI2 interfaces (see also the zoom-in of the pDOS in Fig. S8 in the ESI†). Indeed, the HOMO levels of Spiro-OMeTAD and PTZ1 are slightly higher in energy than the MAPI VB edge, thus granting a convenient extraction of the hole from the photo-excited MAPI.
As for HTMs/MAPI:MAI interfaces (Fig. 4), the VBM of MAPI is higher in energy than the HOMO levels of all HTMs considered. Thus, these interfaces, although energetically very stable, do not provide a sufficient driving force for the holes to be injected into the HTMs.
Considering these results, the prototypical MAPI perovskite has resulted in an inefficient charge transport process when combined with TPA and PTZ-based molecules. For this reason, if it is to be adopted for the construction of PSC devices, the chemical doping of HTM1 could be considered as a valid approach to increase its HOMO energy level. Still, since isolated HTM1 possesses optimal properties to be used as a dopant-free HTM and it is well known that the doping could lead to perovskite degradation,3,4,9,13–15 an alternative and valid strategy could consist in the exploration of different kinds of perovskite compositions to understand the potential of this molecule in PSC devices. As a matter of fact, the properties of dopant-free HTM1 (also named BPT-1) have been recently investigated in a methylammonium-free p–i–n PSC, adopting Cs0.17FA0.83Pb(I0.9Br0.1)3 as the perovskite and achieving 17.26% PCE, while that obtained employing the standard PTAA in the same experimental conditions is found to be 17.96%.28 This result further corroborates the suitability of HTM1 as HTM, as well as the influence of the perovskite formulation on the hole transport process. Additionally, the promising performances of PSC devices built with PTZ2 as HTM have been reported together with the use of dopants but most of all considering the employment of the mixed halide (bromine and iodine) and cations (methylammonium, formamidinium and cesium) perovskite.26 Hence, this is an additional proof that the interaction between the HTMs and the perovskite is not the only parameter influencing the charge transport process. Indeed, a good binding motif does not always provide sufficient driving force for the hole injection/extraction, but also the perovskite composition plays a fundamental role in determining an effective charge transport process.
The calculations revealed that HTM1 has the lowest reorganization energy, hence it could possess good hole mobility, while its Gox is well aligned with those of Spiro-OMeTAD, PTZ1, and PTZ2, which would allow its application as HTM in common PSCs. The investigation of the structural and electronic properties of HTM1/MAPI has shown that HTM1 has a good interaction with both MAPI:PbI2 and MAPI:MAI, when it exposes the methoxybenzene of TPA moieties. In particular, the interaction with MAPI:MAI results in a more stabilized binding energy, due to the formation of hydrogen bonds between the oxygen atoms of methoxybenzene and the hydrogen atoms of methylammonium cations. Nevertheless, from the pDOS plots it has been found that the HOMO level of HTM1 has lower energy than the VBM of MAPI, thus we found that sufficient driving force is not available for the hole injection process to occur. An inadequate driving force for the hole injection/extraction has also been found for PTZ2/MAPI, even if the latter presents the most stabilized binding energies thanks to the presence of additional interactions. In the end, Spiro-OMeTAD and PTZ1 have shown a convenient but low driving force for the hole extraction from the MAPI:PbI2 configuration. Hence, from this study, it emerged that the intrinsic HTM properties and the binding motif are not the only parameters influencing the charge transport process, while the perovskite composition can also strongly affect it. Besides MAPI perovskite is often taken as a reference for the investigation of electronic properties at the interface with HTMs, it has resulted to be inadequate for TPA and PTZ-based molecules. This finding supports the limited employment of MAPI for the real fabrication of devices. For this reason and to shed further light on interfacial processes, additional computational studies involving a different kind of perovskite composition should be performed. Moreover, more TPA and PTZ-based HTM conformations on MAPI surface could be explored to further analyze the dynamics of the charge transport processes.
HTMs reorganization energies in vacuo of compounds HTM1, Spiro-OMeTAD, PTZ1 and PTZ2 have been calculated at the PBE/6-31G(d,p) and B3LYP/6-31G(d,p) levels of theory, on the lowest energy neutral and oxidized geometries optimized at the PBE/6-31G(d,p) and B3LYP/6-31G(d,p) levels of theory, respectively.
To investigate HTMs/MAPI structural and electronic properties, the cubic structure of the MAPI bulk and the four considered compounds have been optimized with Periodic Boundary Conditions (PBC) via DFT at PBE-GGA, using the DZP basis set, including 5d10 semicore electrons for Pb atoms,60 along with Troullier–Martins norm-conserving pseudopotentials,61,62 as implemented in the SIESTA 4.1 program package.63 A mesh cutoff of 400 Ry and 8 × 8 × 8 and 1 × 1 × 1 Γ-centered Monkhorst–Pack k-point grid for sampling the Brillouin zone have been considered respectively for the optimization of the MAPI bulk and for the optimization of HTM1, Spiro-OMeTAD, PTZ1 and PTZ2.
In the same way, the HTMs/MAPI energy minima have been optimized at PBE-GGA, using the DZP basis set, including 5d10 semicore electrons for Pb atoms, along with Troullier–Martins norm-conserving pseudopotentials. A mesh cutoff of 400 Ry and a 1 × 1 × 1 Γ-centered Monkhorst–Pack k-point grid have been considered. SCF energy threshold and total force threshold have been set at 10−5 eV and 0.01 eV Å−1, respectively, for each geometry optimization. The binding energies have been determined through single point calculations on the optimized systems at 8 × 8 × 1 Γ-centered Monkhorst–Pack k-point grid. The energies of all investigated systems have been corrected including the dispersion force calculated by the D3BJ damping scheme. We must note that more accurate quantitative predictions of the electronic properties would require a perturbative approach (e.g. GW) and inclusion of the spin orbit coupling (SOC) effects, which influence both the optical and electronic properties of compounds involving heavy atoms such as Pb. However, such an approach is more computationally expensive than DFT and would not be feasible for our interface systems. On the other hand, the combination of DFT+SOC predicts bandgaps very different from experiments. In particular, SOC addition affects to a larger extent the conduction band (made of Pb, p states) while retaining the position of the iodine (p) states forming the valence band. This results in a strong underestimation of the experimental bandgap.11 Moreover, we expect that the SOC would equally affect all the considered interfaces, so neglecting it in all calculations does not prevent us to make a qualitative comparison of the different HTMs, that is the purpose of this work.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2cp01270g |
This journal is © the Owner Societies 2022 |