First-principles study of the pentacene/Cu(1Â 1Â 1) interface: Adsorption states and vacuum level shifts

We have studied the interaction of pentacene with a Cu(1 1 1) surface using density functional theory (DFT) within a generalized gradient approximation (GGA) and the van der Waals density functional [vdW-DF, M. Dion, H. Rydberg, E. Schröder, D.C. Langreth, B.I. Lundqvist, Phys. Rev. Lett. 92 (2004) 246401]. The adsorption energy is accurately predicted by vdW-DF, while the equilibrium distances between pentacene and the metal substrate (ZC) are overestimated by both GGA and vdW-DF. The work function changes depend significantly on ZC. The experimental work function change can be successfully reproduced by GGA if the experimentally reported adsorption geometry is used, whereas the magnitude of the work function change is underestimated if calculated adsorption geometries are applied. We examined the IDIS model [H. Vázquez, R. Qszwaldowski, P. Pou, J. Ortega, R. Pérez, F. Flores, A. Kahn, Europhys. Lett. 65 (2004) 802] to compare it with the GGA results. The interface dipoles estimated by the IDIS model fairly agree with the GGA results, provided that the adsorption distance is large. On the other hand, they tend to deviate from the GGA results as the adsorption distance becomes smaller, where back donation from the metal surface to the adsorbate occurs. Our analysis reveals that at experimentally reported metal-organic distance, back donation is significant enough to induce polarization of pentacene molecules perpendicular to the surface, which leads to a reduction of the work function. Thus, at the experimentally reported metal-organic distance, the work function change estimated by a simple IDIS model deviates from that calculated by self-consistent GGA calculations. We also found that at the experimentally reported metal-organic distance, the transferred electrons create weak chemical bonds between pentacene and the Cu(1 1 1) surface, illustrating the reactive nature of pentacene.