“…Despite the limitations of PBE, it has been used to describe charge transfer systems, such as perylene-3,4,9,10-tetracarboxylicdiimide (PTCDI), 31 perylene-3,4,9,10tetracarboxylic dianhydride (PTCDA) 32 and C 60 33 on ZnO, as well as a donor-acceptor complex of tetrafluorotetracyanoquinodimethane (F 4 -TCNQ) with t-butyl carbazole-phosphonic acid modified ITO. 34 In these earlier works, the work function changes and charge-transfer characters have been calculated at the PBE level and good agreement with experiment has been reported. [31][32][33] In the gold calculations, a 2 × 2 × 1 Monkhorst-Pack k-point grids was used for geometry optimizations for all unit cells, while 6 × 6 × 1, 3 × 6 × 1, and 3 × 3 × 1 Monkhorst-Pack k-point grids were used for self-consistent total-energy calculations for the smallest, remarkably similar, with differences on the order of 0.01 |e| (see Table S6 in the Supporting Information; note that the Gaussian smearing method was used for the isolated systems and only the Γ-point was employed in such calculations).…”