Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 46561.
Published online 2017 April 20. doi:  10.1038/srep46561
PMCID: PMC5397854

Model potential for the description of metal/organic interface states


We present an analytical one-dimensional model potential for the description of electronic interface states that form at the interface between a metal surface and flat-lying adlayers of π-conjugated organic molecules. The model utilizes graphene as a universal representation of these organic adlayers. It predicts the energy position of the interface state as well as the overlap of its wave function with the bulk metal without free fitting parameters. We show that the energy of the interface state depends systematically on the bond distance between the carbon backbone of the adayers and the metal. The general applicability and robustness of the model is demonstrated by a comparison of the calculated energies with numerous experimental results for a number of flat-lying organic molecules on different closed-packed metal surfaces that cover a large range of bond distances.

The charge transfer at the interface between a metal and a layer of organic molecules plays a decisive role in the functionality of organic semiconductor devices and for future applications of molecular electronics. It depends crucially on the energy alignment and the wave function overlap of electronic states at such interface, which also governs the binding and even the growth of the molecular layer1,2,3. These interface states can either originate from localized molecular orbitals of the organic layer or from delocalized electronic states of the metal. The latter is characteristic for π-conjugated organic molecular layers as became first apparent for perylene-tetracarboxylic-acid-dianhydride(PTCDA) on Ag(111) for which an unoccupied, strongly dispersive interface state has been found4,5. A similar interface states has been also observed for the naphthalene-based variant NTCDA on Ag(111)6, for PTCDA on Ag(100)7 and meanwhile also for a number of other systems8,9,10,11,12,13,14,15,16,17,18,19,20. Time- and angle-resolved two-photon photoemission(2PPE) experiments on PTCDA/Ag(111) concluded from the dispersion and the rather short inelastic lifetime of this state that it must originate from the Shockley surface state of the bare Ag(111) substrate which is upshifted from below the metallic Fermi level by as much as 0.7 eV due to the interaction with the molecular layer5,21. This interpretation was subsequently confirmed by density-functional theory(DFT) calculations6,22,23,24,25,26, which showed that the hybridization of molecular and metallic states is rather small in the region of the projected band gap of the metal.

It turned out, however, that a realistic description of organic molecules on metal surfaces by DFT is challenging although it is one of the most widely used approaches for the determination of the geometric and electronic structure at surfaces and interfaces27. The large size of the organic molecules does not only require a large supercell within a slab model in the lateral direction. A reasonable description of the intrinsic Shockley surface state of the metal makes it necessary to consider also a large number of metal layers24. Both make such calculations very time-consuming. Moreover, metal-organic interfaces require tailored calculations methods26, because conventional DFT neither correctly accounts for van der Waals forces which have an important contribution to the interaction between organic molecules and metal surfaces nor for the correct long-range interaction in front of metal surfaces.

In order to highlight the main physical mechanism for the formation of the delocalized interface state at organic/metal interfaces without the help of complex DFT calculations, we propose in this letter a one-dimensional description by an analytical model potential. The choice of the potential is inspired by previous work on surface states of clean metals28. We recently used it for the description of image-potential states in graphene/metal systems29. Unlike other, more adsorbate-specific model potentials that have been used to describe interfacial electronic states,30,31 its main parameter is simply the distance of the carbon plane from the metal substrate. We will show that the same model potential not only predicts the energy of the interface state in various graphene/metal systems, but can be applied to large class of flat lying molecular layers with a similar π-π interaction as in graphene.

Our model calculation clearly illustrates how the interface state develops from the former Shockley-type surface state of the bare metal substrate with increasing interaction between the molecular film and the metal. By comparing our model results with available experimental data for different organic molecules, we show that our model is able to describe the systematic dependence of the interface state’s energy on the bond distance between the carbon backbone and the metal with predictive power if this distance and the work function are known. Moreover, the model reveals how the wave function overlap of the interface state with both, the bulk metal and the molecular overlayer depends on the carbon-metal distance.

Results and Discussion

Model potential

Our one-dimensional model potential is based on the nearly-free-electron model for the bulk. In this approximation, the electronic states within the metal are given by the solution of the Schrödinger equation under the influence of a weak periodic pseudopotential. Already within this simple model, the formation of the Shockley surface state can be described by introducing a potential barrier at the surface. For illustration of the basic properties of this surface state, we first recall the textbook example of a step barrier at a distance z0 from the position of the topmost atom at z = 032,33. In the direction perpendicular to the surface, the potential is in this case given by

An external file that holds a picture, illustration, etc.
Object name is srep46561-m1.jpg

with the reciprocal lattice vector g = 2π/a and the distance between lattice planes a. With an appropriate choice of the inner potential V0 and the corrugation Vg, this model provides a good approximation of the energy-momentum dispersion E(k) of electron states for the chosen direction in simple metals such as Al, Ag or Cu which are derived from sp electrons34. The model predicts an energy gap 2Vg in which bulk electronic states are forbidden. The bulk eigenfunctions(Bloch states) at the bottom and at the top of this gap correspond to the k-value at the Brillouin zone boundary An external file that holds a picture, illustration, etc.
Object name is srep46561-m2.jpg, and are simply given by An external file that holds a picture, illustration, etc.
Object name is srep46561-m3.jpg. In case of a repulsive potential(Vg > 0), An external file that holds a picture, illustration, etc.
Object name is srep46561-m4.jpg corresponds to the upper and An external file that holds a picture, illustration, etc.
Object name is srep46561-m5.jpg to the lower band. The corresponding eigenvalues are An external file that holds a picture, illustration, etc.
Object name is srep46561-m6.jpg.

The potential barrier at the surface leads to an additional solution with an energy within the bulk band gap. The corresponding eigenfunction is located at the surface and consists of an exponentially decaying cosine Bloch function in the bulk and a simple decaying exponential on the vacuum side,

An external file that holds a picture, illustration, etc.
Object name is srep46561-m7.jpg

with An external file that holds a picture, illustration, etc.
Object name is srep46561-m8.jpg and An external file that holds a picture, illustration, etc.
Object name is srep46561-m9.jpg. The phase shift δ of the bulk cosine function is directly related to the energy E of the surface state. It varies from δ = π/2 at the bottom to δ = 0 at the top of the bulk band gap. The wavefunctions in the bulk and in the vacuum can only be matched for one specific value of the phase δ or the respective energy33. The surface state’s energy depends on both, the height V0 and the position z0 of the surface barrier. It is, however, much more sensitive on the latter. Figure 1 illustrates how the energy of the surface state increases for increasing V0 and decreasing z0. With parameters for Ag(111)(a = 2.36 Å, Vg = 2.15 eV, V0 = 9.56 eV) and z0 = a/2, for example, it turns out that changing the surface state’s energy by 1 eV requires a change of the barrier height as large as 4 eV, but only a change of z0 by ~a/10. The key point for the following discussion is that adlayers of organic molecules modify the distance as well as the height of the barrier, but shift the energy of the surface state predominantly according to their adsorption distance.

Figure 1
Solutions Ψ(z) of the Schrödinger equation for the Shockley surface state(red solid lines) for two different positions z0 and heights V0 of the surface barrier (blue solid lines).

For our model potential we use a more realistic description of the potential barrier at metal surfaces as was introduced by Chulkov et al.28. This approach accounts for the long-range image-potential which is matched to the periodic bulk potential such that the model potential and its derivative is continuous in space. By fitting the matching parameters, not only the work function and the energies of the image-potential states, but also the energy of the Shockley surface state on a number of simple and noble metal surfaces can be quantitatively well reproduced28. Our model potential combines this metal potential with a potential for the molecular adlayer.

Recently, we have proposed such one-dimensional model potential for a description of the image-potential states of graphene(g) on metal substrates29. We could show that the energy of the image-potential states as well the coupling of their wave functions to the metal bulk systematically varies as a function of the carbon-metal distance dC. This potential is composed of four parts

An external file that holds a picture, illustration, etc.
Object name is srep46561-m10.jpg

where Vm(z) denotes the metal potential and Vg(z) the potential of the π-conjugated graphene layer. The latter is a parameterized analytic expression of the potential proposed by Silkin et al.35. dC is the distance of the carbon atoms in the graphene layer with respect to the position of the outermost metal atoms located at z = 0. VΦ(z) and δV(z) are corrections that consider the difference in work function between the bare and the graphene covered metal and the influence of higher-order image-charges, respectively. Beside the distance and change of the workfunction, V(z) is fixed by the separate properties of the metal and the adlayer and does not contain further free fitting parameters. In particular, the metal potential quantitatively reproduce not only the image-potential states, but also the Shockley surface state on the(111) noble metal surfaces.

We apply this model potential at first to a graphene layer on Ag(111) and then show that the results for this system can be directly related to adlayers of flat-lying organic molecules containing carbon rings that have a similar π-π interaction as in graphene as long as the corresponding carbon-metal distance and work function is considered. For the Ag(111) substrate, we use the parameters given in ref. 28. For the work function of the combined system, we use Φ = 4.24 eV at dC = 3.33 Å as reported in ref. 36. Vg(z), VΦ(z) and δV(z) are determined as described in ref. 29. Figure 2a) shows the combined potential for an exemplary carbon-metal distance of dC = 5 Å. The wave functions Ψ and energies E of the Shockley surface state(SS) for bare Ag(111) and the interface state(IS) for graphene covered Ag(111) have been calculated at the center of the surface Brillouin zone(An external file that holds a picture, illustration, etc.
Object name is srep46561-m11.jpg-point) by solving the one-dimensional Schrödinger equation numerically by using Numerov’s method.

Figure 2
(a) One-dimensional model potential V(z)(blue solid line) for a carbon layer on Ag(111) at an exemplary metal-C distance of dC = 5.0 Å. The positions of the uppermost Ag atomic layer and the carbon layer are depicted by ...

We characterize the coupling strength of the interface state to the metal bulk by the the fraction p of its probability density for z < z0 = a/2. For bare Ag(111) at the An external file that holds a picture, illustration, etc.
Object name is srep46561-m12.jpg-point, we calculate p = 76.23% and ESS = −59 meV relative to the Fermi level.


In the following we discuss the transition of the Shockley state of the bare metal into the interface state upon approaching the carbon layer to the Ag(111) surface. For this purpose, we have calculated the wavefunction and energy of the former surface state for different carbon-metal distance dC. In real systems, dC varies in the range of 2.2–3.7 Å(compare Table 1) reflecting strong and weak interaction, respectively. We start with Fig. 2b at a much larger distance of dC = 8 Å, where the surface potential of the bare Ag(111) substrate(dotted line) is substantially modified only at distances that are larger than the extension of the Shockley surface state into the vacuum. At this distance, the probability density of the interface state(red solid line) is still basically identical with that of the Shockley surface state of bare Ag(111), but its energy relative to the Fermi level of EIS = −136 meV is slightly reduced. This is caused by the reduction of the barrier in the region between the metal and the carbon layer as can be seen by comparing the dashed and the solid line in Fig. 2b. For smaller distances dC, there is an interplay between a further reduction of this barrier and the approach of the barrier between the carbon layer and the vacuum closer to the metal. The first leads to a decrease of the interface state’s energy, the latter to an increase. At a distance of dC = 4 Å(Fig. 2c), which is still larger than in real systems, one can already observe the transition from the Shockley surface state of the metal to the actual interface state. Its energy of EIS = −309 meV is, however, even further reduced. Its probability density leaks further into the vacuum and develops two small maxima around the position of the carbon layer. The overlap with the metal bulk decreases to p = 73.83%. For a further reduction of dC, the approach of the barrier between the carbon layer and the vacuum starts to dominate the effective barrier and results in an upshift of the interface state’s energy with decreasing dC. This case is illustrated in Fig. 2d, which shows the calculated probability density for dC = 2.86 Å. This is just the experimental determined distance between the carbon backbone of PTCDA and the metal surface when adsorbed as a flat-lying layer on Ag(111)37. At this distance, the calculated energy of the interface state is now shifted substantially above the Fermi level to EIS = +331 meV. The interface state becomes therefore an unoccupied state. The reallocation of the probability density from the metal surface into the interface and the vacuum region is even stronger as compared to dC = 4.00 Å and the penetration into the metal further decreases to p = 64.70%. The blue solid line shows for comparison the laterally averaged result of a DFT calculation that has been performed for PTCDA on a nine-layer thick Ag(111) slab24. The agreement between our results for the one-dimensional model potential and the DFT-result is remarkably good, in particular within in metal and between the metal and the carbon plane. Above the carbon plane, however, the DFT result extends much less into the vacuum. A part of this difference might be connected to the fact that the DFT calculation does not account for the long-range image-potential that is well described by our model potential. On the other hand, our one-dimensional model does not consider lateral variations of the electronic structure which are most important within the molecular layer.

Table 1
Experimental data on the carbon-metal distance d C and the energy shift ΔE IS between the interface state and the surface state for single layers of flat-lying organic molecules and graphene layers on metal surfaces.

For a comparison with experimental results, we list in Table 1 available data on the energy difference ΔEIS = EISESS between the interface state and the surface state as well as on the carbon-metal distance dC for a variety of adlayers of organic molecules on metal surfaces. In order to emphasize the prototypic character of a graphene layer on a metal surface for these systems, we added also data of such systems. If available, dC has been taken from x-ray standing wave experiments, otherwise from DFT, low energy electron diffraction or surface x-ray diffraction. Fig. 3 depicts the correlation between ΔEIS and dC for the experimental and the calculated results. The red solid line shows the calculated results for a work function of the combined graphene/Ag(111) system of Φ = 4.24 eV at dC = 3.33 Å36. Since the work function difference between the covered and bare metal surfaces varies substantially for the different systems even at comparable dC, we depict with the gray areas the variation of the calculated results when changing the work function by ±1 eV. As can be clearly seen, the work function, i.e. the height of the potential barrier, has only a minor influence on the energy of the interface state. The latter depends instead much more sensitive on dC and shows a strong increase for distances below ~3.25 Å. Our model calculation can well reproduce this trend for the majority of the experimental data. We can not confirm the predicted downshift of the interface state at large dC because this regime is not covered by real systems. The agreement between model calculation and experiment described here is rather good, even for systems like PTCDA and NTCDA on Ag(111) that are chemisorbed. A comparison with the results of DFT calculations for these two cases shows that not only the energies but also the laterally averaged wavefunctions are described well by the model24. The good agreement obtained for these chemisorbed organic layers is surprising only at first glance. Experiments utilizing so-called orbital tomography show that the main effect of chemisorption is simply a filling or partial filling of the former LUMO orbital38,39. Chemisorption does not result in a major distortion of the π system which would cause the model to break down.

Figure 3
Energy shift ΔEIS of the interface state with respect to the energy of the former surface state on the bare metal as a function of the carbon-metal distance dC.

Systems with stronger bonds of the functional groups of the organic molecules typically show a pronounced bending of the molecules towards the surface which reduces the carbon-metal distance at the edges40. Because we relate dC to the center of the carbon backbone, this might explain the larger deviation between model and experiment for PTCDA/Ag(100). For this system, also DFT calculations7 clearly underestimate ΔEIS(c.f. Table 1). Even if our simple one-dimensional model gives a reasonable explanation of the main physical mechanism for the interface state formation in a number of systems, it can, however, not account for more complex chemical interactions between the substrate and the molecules. This might explain the more pronounced deviations between model and experiment for the phthalocyanines ZnPc/Cu(111) and F16ZnPc/Cu(111) that are subject to a stronger corrugation or distortion of the molecular layer on the substrate41. With the exception of these two phthalocyanines molecules, the model calculations slightly underestimate the IS energy for most systems. These small, but systematic deviations are most likely caused by molecule-metal interactions that are not taken into account by the model. Adsorption modifies the charge distribution of both the metal surface and the molecular layer. One well-known consequence is the formation of a surface dipole by the so-called cushion effect42 which our model takes into account via the work function in VΦ(z). However, also the metal as well as the graphen part of the potential will deviate slightly from the functional form of the isolated systems Vm(z) and Vg(z  dC) used in our model. One expects Vg(z  dC) to become asymmetric with respect to the carbon plane and, most importantly, the position of the image-plane in the Chulkov potential Vm(z) to move closer the substrate thus leading to an additional upshift of the IS. In order to keep the model potential parameter free we refrained from correcting for this overall small effect. A further possibility for the underestimation of the model to the experimental data is a shift of the last atomic layer of the substrate after molecular adsorption compared to the bulk values. Such a shift is not implemented in the model but it is known that the energetic position of the Shockley state depends on the lattice constant of the substrate43.

Finally, we note that the present one-dimensional model can, in principle, be used to also predict the effective mass of interface states, i.e. the dispersion An external file that holds a picture, illustration, etc.
Object name is srep46561-m13.jpg near the An external file that holds a picture, illustration, etc.
Object name is srep46561-m14.jpg point k|| = 0. For this purpose, one simply has to solve the Schrödinger equation not only for one potential V(z) containing the metallic part at the An external file that holds a picture, illustration, etc.
Object name is srep46561-m15.jpg point An external file that holds a picture, illustration, etc.
Object name is srep46561-m16.jpg, as done above, but for a whole set of An external file that holds a picture, illustration, etc.
Object name is srep46561-m17.jpg. This set of potentials is readily available for all the metals included in Table 1 as their band structure is well known. One should, however, keep in mind that this approach neglects any lateral corrugation of the molecular part Vg(z). It should thus be restricted to the vicinity of the An external file that holds a picture, illustration, etc.
Object name is srep46561-m18.jpg point.

Additional Information

How to cite this article: Armbrust, N. et al. Model potential for the description of metal/organic interface states. Sci. Rep. 7, 46561; doi: 10.1038/srep46561(2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We gratefully acknowledge funding by the Deutsche Forschungsgemeinschaft through SFB 1083 Structure and Dynamics of Internal Interfaces, Project B6. F. S. acknowledges support by the CSIC through the I-Link 1065 program.


The authors declare no competing financial interests.

Author Contributions U.H. conceived the model, N.A. carried out the calculations. All authors discussed, wrote, and reviewed the manuscript.


  • Koch N., Ueno N. & Wee A. T. S.(eds.) The Molecule-Metal Interface(Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, 2013).
  • Lindstrom C. D. & Zhu X. Y. Photoinduced Electron Transfer at Molecule-Metal Interfaces. Chem. Rev. 106, 4281–300(2006). [PubMed]
  • Hwang J., Wan A. & Kahn A. Energetics of Metal-Organic Interfaces: New Experiments and Assessment of the Field. Mat. Sci. Eng. R. 64, 1–31(2009).
  • Temirov R., Soubatch S., Luican A. & Tautz F. S. Free-Electron-Like Dispersion in an Organic Monolayer Film on a Metal Substrate. Nature 444, 350–53(2006). [PubMed]
  • Schwalb C. H. et al. . Electron Lifetime in a Shockley-Type Metal-Organic Interface State. Phys. Rev. Lett. 101, 146801(2008). [PubMed]
  • Marks M. et al. . Energy shift and wave function overlap of metal-organic interface states. Phys. Rev. B 84, 081301(2011).
  • Galbraith M. C. E., Marks M., Tonner R. & Höfer U. Formation of an Organic/Metal Interface State from a Shockley Resonance. J. Phys. Chem. Lett. 5, 50(2014). [PubMed]
  • Tamai A. et al. . Electronic Structure at the C60/Metal Interface: An Angle-Resolved Photoemission and First-Principles Study. Phys. Rev. B 77, 075134(2008).
  • Yamane H., Kanai K., Ouchi Y., Ueno N. & Seki K. Impact of Interface Geometric Structure on Organic-Metal Interface Energetics and Subsequent Films Electronic Structure. J. Electron Spectrosc. 174, 28–34(2009).
  • Scheybal A. et al. . Modification of the Cu(110) Shockley Surface State by an Adsorbed Pentacene Monolayer. Phys. Rev. B 79, 115406(2009).
  • Ziroff J., Gold P., Bendounan A., Forster F. & Reinert F. Adsorption Energy and Geometry of Physisorbed Organic Molecules on Au(111) Probed by Surface-State Photoemission. Surf. Sci. 603, 354–58(2009).
  • Andrews K. M. & Pearl T. P. Modification of Ag(111) Surface Electronic Structure Via Weak Molecular Adsorption of Adenine Measured with Low Temperature Scanning Tunneling Microscopy and Spectroscopy. J. Chem. Phys. 132, 214701(2010). [PubMed]
  • Heinrich B. W. et al. . Dispersion and Localization of Electronic States at a Ferrocene/Cu(111) Interface. Phys. Rev. Lett. 107, 216801(2011). [PubMed]
  • Schmid M., Zirzlmeier J., Steinrück H. P. & Gottfried J. M. Interfacial Interactions of Iron(II) Tetrapyrrole Complexes on Au(111). J. Phys. Chem. C 115, 17028–17035(2011).
  • Hong S. Y., Yeh P. C., Dadap J. I. & Osgood R. M. Interfacial Dipole Formation and Surface-Electron Confinement in Low-Coverage Self-Assembled Thiol Layers: Thiophenol and p-Fluorothiophenol on Cu(111). Acs Nano 6, 10622–10631(2012). [PubMed]
  • Faraggi M. N. et al. . Bonding and Charge Transfer in Metal–Organic Coordination Networks on Au111) with Strong Acceptor Molecules. J. Phys. Chem. C 116, 24558–24565(2012).
  • Wiessner M. et al. . Different views on the electronic structure of nanoscale graphene: aromatic molecule versus quantum dot. New J. Phys. 14, 113008(2012).
  • Umbach T. R., Fernandez-Torrente I., Ladenthin J. N., Pascual J. I. & Franke K. J. Enhanced charge transfer in a monolayer of the organic charge transfer complex TTF–TNAP on Au(111). J. Phys.: Condens. Matter 24, 354003(2012). [PubMed]
  • Ilyas N., Kelly L. L. & Monti O. L. A. Electronic structure and dynamics of quasi-2D states of vanadyl naphthalocyanine on Au(111). Mol. Phys. 111, 2175–2188(2013).
  • Heidorn S. et al. . Influence of Substrate Surface-Induced Defects on the Interface State between NaCl(100) and Ag(111). Journal of Physical Chemistry C 117, 16095–16103(2013).
  • Sachs S. et al. . Electronic structure at the perylene-tetracarboxylic acid dianhydride/Ag(111) interface studied with two-photon photoelectron spectroscopy. J. Chem. Phys. 131, 144701(2009). [PubMed]
  • Dyer M. S. & Persson M. The nature of the observed free-electron-like state in a PTCDA monolayer on Ag(111). New J. Phys. 12, 063014(2010).
  • Zaitsev N. L., Nechaev I. A. & Chulkov E. V. Change in surface states of Ag(111) thin films upon adsorption of a monolayer of PTCDA organic molecules. J. Exp. Theor. Phys. 110, 114–20(2010).
  • Zaitsev N. L., Nechaev I. A., Echenique P. M. & Chulkov E. V. Transformation of the Ag(111) Surface State Due to Molecule-Surface Interaction with Ordered Organic Molecular Monolayers. Phys. Rev. B 85, 115301(2012).
  • Tsirkin S. S. et al. . Inelastic decay of electrons in Shockley-type metal-organic interface states. Phys. Rev. B 92, 235434(2015).
  • Jakob P. et al. . Adsorption geometry and interface states: Relaxed and compressed phases of NTCDA/Ag(111). Phys. Rev. B 94, 125436(2016).
  • Pilania G., Zhu H. & Ramprasad R. It Applications of Modern Density Functional Theory to Surfaces and Interfaces. 271–312(John Wiley & Sons, Inc., Hoboken, NJ, USA, 2012).
  • Chulkov E. V., Silkin V. M. & Echenique P. M. Image Potential States on Metal Surfaces: Binding Energies and Wave Functions. Surf. Sci. 437, 330–52(1999).
  • Armbrust N., Güdde J. & Höfer U. Formation of image-potential states at the graphene/metal interface. New J. Phys. 17, 103043(2015).
  • Berthold W., Rebentrost F., Feulner P. & Höfer U. Influence of Ar, Kr, and Xe layers on the energies and lifetimes of image-potential states on Cu(100). Appl. Phys. A 78, 131(2004).
  • Zhao J., Feng M., Dougherty D. B., Sun H. & Petek H. Molecular Electronic Level Alignment at Weakly Coupled Organic Film/Metal Interfaces. ACS Nano 8, 10988–97(2014). [PubMed]
  • Maue A. W. Die Oberflächenwellen in der Elektronentheorie der Metalle. Z. Phys. A-Hadron. Nucl. 94, 717–41(1935).
  • Zangwill A. Physics at Surfaces(Cambridge University Press, Cambridge, 1988).
  • Ashcroft N. W. & Mermin N. D. Solid State Physics(Holt, Rinehart and Winston, New York, 1976).
  • Silkin V. M. et al. . Image Potential States in Graphene. Phys. Rev. B 80, 121408(2009).
  • Giovannetti G. et al. . Doping Graphene with Metal Contacts. Phys. Rev. Lett. 101, 026803(2008). [PubMed]
  • Hauschild A. et al. . Molecular Distortions and Chemical Bonding of a Large π-Conjugated Molecule on a Metal Surface. Phys. Rev. Lett. 94, 036106(2005). [PubMed]
  • Wießner M. et al. . Complete determination of molecular orbitals by measurement of phase symmetry and electron density. Nature Communications 5, 4156(2014). [PubMed]
  • Weiß S. et al. . Exploring three-dimensional orbital imaging with energy-dependent photoemission tomography. Nature Communications 6, 8287(2015). [PMC free article] [PubMed]
  • Bauer O. et al. . Role of Functional Groups in Surface Bonding of Planar π-Conjugated Molecules. Phys. Rev. B 86, 235431(2012).
  • Yamane H. et al. . Site-Specific Geometric and Electronic Relaxations at Organic-Metal Interfaces. Phys. Rev. Lett. 105, 046103(2010). [PubMed]
  • Witte G., Lukas S., Bagus P. S. & Wöll C. Vacuum Level Alignment at Organic/Metal Junctions: “Cushion” Effect and the Interface Dipole. Appl. Phys. Lett. 87, 263502(2005).
  • Zaitsev N. L., Nechaev I. A., Höfer U. & Chulkov E. V. Adsorption geometry and electronic properties of flat-lying monolayers of tetracene on the Ag(111) surface. Phys. Rev. B 94, 155452(2016).
  • Hauschild A. et al. . Normal-incidence x-ray standing-wave determination of the adsorption geometry of PTCDA on Ag(111): Comparison of the ordered room-temperature and disordered low-temperature phases. Phys. Rev. B 81, 125432(2010).
  • Stadtmüller B. et al. . Commensurate Registry and Chemisorption at a Hetero-Organic Interface. Phys. Rev. Lett. 108, 106103(2012). [PubMed]
  • To be published.
  • Stadler C. et al. . Molecular distortion of NTCDA upon adsorption on Ag(111): a normal incidence x-ray standing wave study. New J. Phys. 9, 50(2007).
  • Kröger I., Stadtmüller B. & Kumpf C. Submonolayer and multilayer growth of titaniumoxide-phthalocyanine on Ag(111). New J. Phys. 18, 111002(2016).
  • Kröger I. et al. . Submonolayer Growth of H2-Phthalocyanine on Ag(111). Phys. Rev. B 86, 195412(2012).
  • Caplins B. W., Suich D. E., Shearer A. J. & Harris C. B. Metal/Phthalocyanine Hybrid Interface States on Ag(111). J. Phys. Chem. Lett. 5, 1679–84(2014). [PubMed]
  • Soubatch S., Kröger I., Kumpf C. & Tautz F. S. Structure and Growth of Tetracene on Ag(111). Phys. Rev. B 84, 195440(2011).
  • Duhm S. et al. . Influence of Intramolecular Polar Bonds on Interface Energetics in Perfluoro-Pentacene on Ag(111). Phys. Rev. B 81, 045418(2010).
  • Henze S. K. M., Bauer O., Lee T. L., Sokolowski M. & Tautz F. S. Vertical bonding distances of PTCDA on Au(111) and Ag(111): Relation to the bonding type. Surf. Sci. 601, 1566–73(2007).
  • Armbrust N., Güdde J., Jakob P. & Höfer U. Time-Resolved Two-Photon Photoemission of Unoccupied Electronic States of Periodically Rippled Graphene on Ru(0001). Phys. Rev. Lett. 108, 056801(2012). [PubMed]
  • Wang B., Bocquet M. L., Güenther S. & Wintterlin J. Comment on “Periodically Rippled Graphene: Growth and Spatially Resolved Electronic Structure”. Phys. Rev. Lett. 101, 099703(2008). [PubMed]
  • Wang B., Bocquet M. L., Marchini S., Günther S. & Wintterlin J. Chemical Origin of a Graphene Moire Overlayer on Ru(0001). Phys. Chem. Chem. Phys. 10, 3530–34(2008). [PubMed]
  • Moritz W. et al. . Structure Determination of the Coincidence Phase of Graphene on Ru(0001). Phys. Rev. Lett. 104, 136102(2010). [PubMed]
  • Pagliara S. et al. . Nature of the Surface States at the Single-Layer Graphene/Cu(111) and Graphene/Polycrystalline-Cu Interfaces. Phys. Rev. B 91, 195440(2015).
  • Jeon C. et al. . Rotated Domains in Chemical Vapor Deposition-Grown Monolayer Graphene on Cu(111): An Angle-Resolved Photoemission Study. Nanoscale 5, 8210–14(2013). [PubMed]
  • Jolie W., Craes F. & Busse C. Graphene on Weakly Interacting Metals: Dirac States Versus Surface States. Phys. Rev. B 91, 115419(2015).
  • Jean F. et al. . Topography of the Graphene/Ir(111) Moire Studied by Surface X-Ray Diffraction. Phys. Rev. B 91, 245424(2015).
  • Niesner D. et al. . Trapping Surface Electrons on Graphene Layers and Islands. Phys. Rev. B 85, 081402(2012).
  • Reinert F., Nicolay G., Schmidt S., Ehm D. & Hüfner S. Direct Measurements of the L-Gap Surface States on the(111) Face of Noble Metals by Photoelectron Spectroscopy. Phys. Rev. B 63, 115415(2001).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group