Diese Arbeit befasst sich mit der Untersuchung und Manipulation von Halbleitern, bei denen die Spin-Bahn-Kopplung (SBK) in Kombination mit einem Bruch der strukturellen Inversionssymmetrie zu einer impulsabhängigen Spinaufspaltung der Bandstruktur führt. Von besonderem Interesse ist hierbei der Zusammenhang zwischen der spinabhängigen elektronischen Struktur und der strukturellen Geometrie. Dieser wird durch eine Kombination komplementärer, oberflächensensitiver Messmethoden - insbesondere Rastertunnelmikroskopie (STM) und Photoelektronenspektroskopie (PES) - an geeigneten Modellsystemen untersucht. Der experimentelle Fokus liegt dabei auf den polaren Halbleitern BiTeX (X =I, Br, Cl). Zusätzliche Experimente werden an dünnen Schichten der topologischen Isolatoren (TI) Bi1,1-xSb0;9+xSe3 (x = 0. . . 1,1) und Bi2Te2Se durchgeführt. Die inversionsasymmetrische Kristallstruktur in BiTeX führt zur Existenz zweier nicht-äquivalenter Oberflächen mit unterschiedlicher Terminierung (Te oder X) und invertierter atomarer Stapelfolge. STM-Aufnahmen der Oberflächen gespaltener Einkristalle belegen für BiTeI(0001) eine Koexistenz beider Terminierungen auf einer Längenskala von etwa 100 nm, die sich auf Stapelfehler im Kristallvolumen zurückführen lassen. Diese Domänen sind groß genug, um eine vollständig entwickelte Banddispersion auszubilden und erzeugen daher eine Kombination der Bandstrukturen beider Terminierungen bei räumlich integrierenden Messmethoden. BiTeBr(0001) und BiTeCl(0001) hingegen zeichnen sich durch homogene Terminierungen auf einer makroskopischen Längenskala aus. Atomar aufgelöste STM-Messungen zeigen für die drei Systeme unterschiedliche Defektdichten der einzelnen Lagen sowie verschiedene strukturelle Beeinflussungen durch die Halogene. PES-Messungen belegen einen starken Einfluss der Terminierung auf verschiedene Eigenschaften der Oberflächen, insbesondere auf die elektronische Bandstruktur, die Austrittsarbeit sowie auf die Wechselwirkung mit Adsorbaten. Die unterschiedliche Elektronegativität der Halogene resultiert in verschieden starken Ladungsübergängen innerhalb der kovalent-ionisch gebundenen BiTe+ X- Einheitszelle. Eine erweiterte Analyse der Oberflächeneigenschaften ist durch die Bedampfung mit Cs möglich, wobei eine Änderung der elektronischen Struktur durch die Wechselwirkung mit dem Alkalimetall studiert wird. Modifiziert man die Kristallstruktur sowie die chemische Zusammensetzung von BiTeI(0001) nahe der Oberfläche durch Heizen im Vakuum, bewirkt dies eine Veränderung der Bandstruktur in zwei Schritten. So führt zunächst der Verlust von Iod zum Verlust der Rashba-Aufspaltung, was vermutlich durch eine Aufhebung der Inversionsasymmetrie in der Einheitszelle verursacht wird. Anschließend bildet sich eine neue Kristallstruktur, die topologisch nichttriviale Oberflächenzustände hervorbringt. Der Umordnungsprozess betrifft allerdings nur die Kristalloberfläche - im Volumen bleibt die inversionsasymmetrische Einheitszelle erhalten. Einem derartigen Hybridsystem werden bislang unbekannte elektronische Eigenschaften vorausgesagt. Eine systematische Untersuchung von Dünnschicht-TIs, die mittels Molekularstrahlepitaxie (MBE) erzeugt wurden, zeigt eine Veränderung der Morphologie und elektronischen Struktur in Abhängigkeit von Stöchiometrie und Substrat. Der Vergleich zwischen MBE und gewachsenen Einkristallen offenbart deutliche Unterschiede. Bei einem der Dünnschichtsysteme tritt sogar eine lokal inhomogene Zustandsdichte im Bindungsenergiebereich des topologischen Oberflächenzustands auf.

Im Rahmen der vorliegenden Dissertation werden mit unterschiedlichen Analysemethoden die Korrelationen zwischen den strukturellen, elektronischen und magnetischen Eigenschaften von Selten Erd-basierten intermetallischen Oberflächenlegierungen anhand der beiden Probensysteme LaPt$_5$/Pt(111) und CePt$_5$/Pt(111) untersucht. Darüber hinaus werden die strukturellen Eigenschaften von Adsorbat-induzierten Oberflächenrekonstruktionen im sub-ML Bereich in reduzierten Dimensionen auf der Halbleiteroberfläche Si(111) anhand der beiden Materialsysteme Si(111)-(5$\times$2)-Au und Si(111)-($\sqrt{3}\times\sqrt{3}$)R30${\degree}$-Sn mit der Methode LEED-IV analysiert. Das erste experimentelle Kapitel dieser Arbeit behandelt die intermetallische Oberflächenlegierung LaPt$_5$/Pt(111), die sich ausbildet wenn La-Atome auf einem sauberen Pt(111)-Substrat abgeschieden werden und anschließend thermische Energie hinzugefügt wird. Die Dicke der gebildeten Legierung lässt sich über die zuvor angebotene Menge an La-Atomen variieren und resultiert aufgrund der Gitterfehlanpassung von Pt(111) und den obenauf liegenden LaPt$_5$-Filmen in sechs unterschiedliche Beugungsmuster im LEED, deren Überstrukturvektoren durch zwei unterschiedliche Rotationsausrichtungen in Bezug auf das Gitter des Substrats und unterschiedlichen lateralen Gitterkonstanten der Filme gekennzeichnet sind. Die atomare Struktur kann auf eine gemeinsame Kristallstruktur zurückgeführt werden, deren Stöchiometrie aus dickenabhängigen AES-Messungen zu LaPt$_5$ mit einer Pt-reichen Oberflächenabschlusslage bestimmt werden konnte. Die Ergebnisse einer durchgeführten LEED-IV Studie bestätigen das Wachstum der Filme in der CaCu$_5$-Struktur, wobei die Oberflächenterminierungslage im Vergleich zum Volumengitter ein zusätzliches Pt-Atom pro Einheitszelle aufweist, das zusätzlich um einen Wert von \unit{0.26}{\angstrom} aus der Oberfläche hervorsteht. Die La-Atome, die direkt unterhalb der Terminierungslage liegen, erfahren eine Verschiebung in entgegengesetzter Richtung, so dass im Vergleich zum Volumen der Filme eine lokal veränderte Symmetrie im oberflächennahen Bereich vorherrscht und sich auf die elektronischen Eigenschaften der LaPt$_5$-Filme auswirkt. Darüber hinaus wurden die Schwingungseigenschaften der LaPt$_5$-Filme mittels der polarisierten in situ Raman-Spektroskopie bestimmt, bei der die auftretenden Schwingungspeaks durch die Kenntnis der atomaren Struktur und mit Überlegungen aus der Gruppentheorie unterschiedlichen Tiefenbereichen der LaPt$_5$-Filme (Volumen und Oberfläche) zugewiesen werden konnten. Im zweiten experimentellen Kapitel liegt der Fokus auf der atomaren Struktur sowie auf den elektronischen und magnetischen Eigenschaften des Kondo- und Schwerfermionensystems CePt$_5$/Pt(111). In Abhängigkeit von der vor dem Legierungsprozess angebotenen Menge an Ce-Atomen auf dem Pt(111)-Substrat konnten insgesamt sieben verschiedene LEED-Phasen der CePt$_5$-Filme identifiziert werden, deren jeweilige Oberflächenrekonstruktionen durch eine unterschiedliche Rotationsausrichtung in Bezug auf das Pt(111)-Substrat gekennzeichnet sind. Zusätzlich ist die laterale Gitterkonstante einem Prozess aus Verspannung und Dehnung aufgrund der Gitterfehlanpassung von Film und Substrat ausgesetzt. Eine durchgeführte LEED-IV Analyse bestätigt das Wachstum der Filme in der CaCu$_5$-Struktur mit einer Pt-reichen Oberflächenabschlusslage, deren Pt$_3$-Kagom\'{e}-Lage im Vergleich zum Volumengitter mit einem zusätzlichen Pt-Atom pro Einheitszelle gefüllt ist. Die strukturellen Ergebnisse stimmen mit erzielten Resultaten aus früheren Arbeiten überein und verdeutlichen zudem die isostrukturellen Eigenschaften zur intermetallischen Oberflächenlegierung LaPt$_5$/Pt(111). Dies ermöglicht durch geeignete Vergleichsexperimente an LaPt$_5$/Pt(111) die induzierten Phänomene der $4f$-Elektronen bezüglich des Kondo- und Schwerfermionenverhaltens bei CePt$_5$/Pt(111) zu bestimmen, da La-Atome in ihrem atomaren Aufbau keine $4f$-Elektronen beherbergen. Mit der polarisierten in situ Raman-Spektroskopie aufgenommene Spektren anhand von unterschiedlich dicken CePt$_5$-Filmen beinhalten sowohl charakteristische Schwingungspeaks als auch elektronische Übergänge. Das spektroskopische Verhalten der Schwingungspeaks zeigt dabei nicht nur Gemeinsamkeiten zu LaPt$_5$/Pt(111) bei der Zuweisung der Schwingungsmoden zu den jeweiligen Tiefenbereichen in den CePt$_5$-Filmen, sondern es treten auch Unterschiede auf, da eine CePt$_5$-Schwingungsmode einem anormalen Temperaturverhalten unterliegt, das auf die Wechselwirkung mit den $4f$-Elektronen zurückzuführen ist. Weitere spezifische Raman-Signaturen, die elektronischen Übergängen in Form von Kristallfeldniveauaufspaltungen der $4f$-Elektronen von Ce zugewiesen werden konnten, resultieren ebenfalls aus unterschiedlichen Regionen der CePt$_5$-Filme (Oberfläche, inneres Volumen, Interface). Die magnetischen Eigenschaften der CePt$_5$-Filme wurden mit XAS und XMCD an den Ce M$_{4,5}$-Kanten in Abhängigkeit von der Temperatur, dem Einfallswinkel, der Filmdicke und der Stärke des Magnetfelds analysiert. Die markanten Übergänge zwischen unterschiedlichen Curie-Weiss-Regimen in der inversen Suszeptibilität erlauben Rückschlüsse über das Kristallfeldaufspaltungsschema, die Kondo- und die RKKY-Wechselwirkung und korrelieren mit der Ce-Valenz. Zudem konnte bei tiefen Temperaturen ein Übergang in den kohärenten Schwerfermionen-Zustand für alle untersuchten CePt$_5$-Filmdicken in dieser Arbeit nachgewiesen werden. Durch die Vorhersage eines metamagnetischen Lifshitz-Übergangs für diese Filme, der sich in der Magnetfeldabhängigkeit des magnetischen Moments äußert, konnte durch die Aufnahme von Magnetisierungskurven bei tiefen Temperaturen und hohen Magnetfeldern auf zwei weitere charakteristische Energieskalen der renormalisierten Bandstruktur zugegriffen werden. Das dritte experimentelle Kapitel widmet sich der mit LEED und LEED-IV durchgeführten Aufklärung der atomaren Struktur eines quasi-eindimensionalen Elektronensystems, bei dem sich die gebildeten Au-Nanodrähte auf der Si(111)-Oberfläche durch eine Si(111)-(5$\times$2)-Au Rekonstruktion beschreiben lassen. Die aufgenommenen LEED-Bilder mit ihren markanten Beugungsreflexen und sogenannten Streifen deuten auf drei gleichwertige Rotationsdomänen, die jeweils um einen Winkel von \unit{120}{\degree} gegeneinander gedreht sind, auf der Oberfläche hin. Zudem konnte aus einer Simulation der Beugungsbilder das Auftreten von Streifen durch drei zusätzliche Spiegeldomänen, die eine Phasenverschiebung von einem halben Überstrukturvektor einführen und bei einer sorgfältigen LEED-IV Analyse ebenfalls berücksichtigt werden sollten, erklärt werden. Aus den in der Literatur nach einer zweiten Rekalibrierung der nötigen Menge an Au-Atomen zur Ausbildung der Si(111)-(5$\times$2)-Au Rekonstruktion in den letzten Jahren heftig diskutierten Strukturmodellen gibt das von Kwon und Kang aufgestellte Geometriemodell (KK-Modell) die beobachteten energieabhängigen Intensitätsmodulationen in den experimentellen Daten beim Vergleich mit theoretisch berechneten IV-Kurven am besten wieder. Für dieses Modell nimmt der R-Faktor nach Pendry bei den unabhängig voneinander betrachteten drei Energieserien unter verschiedenen Einfallswinkeln der Elektronen auf die Probenoberfläche stets den kleinsten Wert an. Unter der expliziten Berücksichtigung von Si-Adatomen, die sich zusätzlich auf der Oberfläche befinden und in einer (5$\times$4)-Einheitszelle beschrieben werden können, bleibt das KK-Modell das zu präferierende Strukturmodell zur Beschreibung der ausgebildeten Au-Ketten und der Si-Honigwabenstruktur bei der Si(111)-(5$\times$2)-Au Oberflächenrekonstruktion. Im letzten experimentellen Kapitel wird ein zweidimensionales Elektronensystem -- die $\alpha$-Si(111)-($\sqrt{3}\times\sqrt{3}$)R30${\degree}$-Sn Oberflächenrekonstruktion, die sich bei 1/3 ML an Sn-Adsorbaten auf dem Si(111)-Substrat ausbildet -- im Hinblick auf die atomare Struktur bei Raumtemperatur mit LEED und LEED-IV untersucht. Aus den insgesamt sechs in die Analyse aufgenommenen Strukturmodellen, bei denen die Sn-Atome innerhalb der rekonstruierten ($\sqrt{3}\times\sqrt{3}$)R30${\degree}$-Einheitszelle unterschiedliche Adsorptionsplätze auf einer ideal terminierten Si(111)-Oberfläche einnehmen, konnte ein Legierungsverhalten, wie es bei der $\gamma$-Si(111)-($\sqrt{3}\times\sqrt{3}$)R30${\degree}$-Sn Phase auftritt, ausgeschlossen werden. Die Sn-Atome ordnen sich ausschließlich auf der Oberfläche neu an und führen zu einer Relaxation des darunterliegenden Substrats, deren atomare Verschiebungen sich bis in die sechste Si-Lage nachverfolgen lassen. Im Vergleich zu früheren Strukturaufklärungen an diesem Materialsystem bestätigt diese Analyse, dass sich die abgeschiedenen Sn-Atome auf T$_4$-Adsorptionsplätzen energetisch günstig anlagern, wobei die bei drei unterschiedlichen Einfallswinkeln aufgenommenen experimentellen Daten an unterschiedlichen Probenpositionen auf ein vorhandenes bzw. fehlendes Si-Atom auf einem S$_5$-Gitterplatz im darunterliegenden Si(111)-Substrat hindeuten. Außerdem konnte das theoretisch vorhergesagte dynamische Fluktuations-Modell aufgrund der sehr stark erhöhten thermischen Auslenkungen der Sn-Atome aus ihrer Gleichgewichtslage in den Modellrechnungen zur dynamischen Streutheorie nachgewiesen werden. Dies könnte neben den unregelmäßig angeordneten Si-Fehlstellen eine Ursache für das Ausbleiben des strukturell reversiblen Phasenübergangs von einer ($\sqrt{3}\times\sqrt{3}$)R30${\degree}$-Phase zu einer (3$\times$3)-Phase bei tiefen Temperaturen, wie er beispielsweise beim elektronisch vergleichbaren Adsorbatsystem Ge(111)-($\sqrt{3}\times\sqrt{3}$)R30${\degree}$-Sn auftritt, sein.

Ion-specific effects occur in a huge variety of aqueous solutions of electrolytes and larger molecules like peptides, altering properties such as viscosity, enzyme activity, protein stability, and salting-in and salting-out behavior of proteins. Typically, these type of effects are rationalized in terms of the Hofmeister series, which originally orders cations and anions according to their ability to enhance or suppress the solubility of proteins in water. This empirical order, however, is still not understood yet. Quite some effort was made to gain a molecular level understanding of this phenomenon, yet no consensus has been found about the underlying mechanisms and the determination and localization of the interaction sites. Resonant inelastic soft x-ray scattering (RIXS) combines x-ray emission (XES) and absorption spectroscopies (XAS), probing the partial local density of states of both occupied and unoccupied electronic states and is thus a promising candidate to shed more light onto the issue. The studies presented in this work are directed towards an improved understanding of the interaction between salts and peptides. In order to address this topic, the impact of different physical environments on the electronic structure of small molecules (i.e., methanol and glycine derived peptides) is investigated systematically using soft x-ray spectroscopic methods, corroborated with density functional theory (DFT) calculations. In a first step, molecules without any interactions to the surrounding are investigated, using gas-phase methanol as a model system. Thereby, the local and element specific character of RIXS is demonstrated and used to separately probe the local electronic structure of methanol’s hydroxyl and methyl group, respectively. The attribution of the observed emission features to distinct molecular orbitals is confirmed by DFT calculations, which also quantitatively explain the different relative intensities of the emission features. For resonant excitation of the O K pre-edge absorption resonance, strong isotope effects are found that are explained by dynamical processes at the hydroxyl group. This serves as an excellent example for possible consequences of a local change in the geometric structure or symmetry of a molecule on its electronic structure. In the following, the sample system is expanded to the amino acid glycine and its smallest derived peptides diglycine and triglycine. As a first step, they are studied in their crystalline form in solid state. Again, a comprehensive picture of the electronic structure is developed by measuring RIXS maps at the oxygen and nitrogen K absorption edge, corroborated by DFT calculations. Similar to the case of methanol, dynamic processes at the protonated amino group of the molecules after exciting the nitrogen atom have a strong influence on the emission spectra. Furthermore, it is shown that RIXS can be used to selectively excite the peptide nitrogen to probe the electronic structure around it. A simple building block approach for XES spectra is applied to separate the contribution of the emission attributed to transitions into core holes at the peptide and the amino nitrogen, respectively. In the aqueous solution, the surrounding water molecules slightly change the electronic structure, probably via interactions with the charged functional groups. The effects on the x-ray emission spectra, however, are rather small. Much bigger changes are observed when manipulating the protonation state of the functional groups by adjusting the pH value of the solution. A protonation of the carboxyl group at low pH values, as well as a deprotonation of the amino group at high pH values lead to striking changes in the shape of the RIXS maps. In a comprehensive study of glycine’s XES spectra at varying pH values, changes in the local electronic structure are not only observed in the immediate surrounding of the manipulated functional groups but also in more distant moieties of the molecule. Finally, the study is extended to mixed aqueous solutions of diglycine and a variety of different salts as examples for systems where Hofmeister effects are observed. To investigate the influence of different cations and anions on the electronic structure of diglycine, two series of chlorine and potassium salts are used. Ion-specific effects are identified for both cases. Some of the changes in the x-ray emission spectra of diglycine in the mixed solutions qualitatively follow the Hofmeister series as a function of the used salt. The observed trends thereby indicate an increased interaction between the electron density around the peptide oxygen with the cations, whereas anions seem to interact with the amino group of the peptide.

Im Rahmen der vorliegenden Arbeit werden mit einem Rastertunnelmikroskop (RTM) Ladungsdichtemodulationen (LDM) auf Oberflächen von drei verschiedenen Probensystemen untersucht. Bei den Proben handelt es sich um Chrom auf Wolfram(110), Iridiumditellurid (IrTe2) als Volumenmaterial und Eisen auf Rhodium(001). Es werden sowohl die Temperaturabhängigkeit der Phasenübergänge als auch die Wechselwirkung zwischen magnetischen und elektronischen Eigenschaften analysiert. Chrom (Cr) ist ein einfaches Übergangsmetall, in dem sowohl eine klassische Ladungsdichtewelle (LDW) als auch eine Spindichtewelle (SDW) auftreten. Die im Experiment betrachteten Cr-Inseln auf Wolfram(110) schlagen eine Brücke zwischen dem Volumenmaterial und ultradünnen Schichten. Dabei zeigt sich der Zusammenhang zwischen elektronischen und magnetischen Eigenschaften in der Ausbildung einer LDW-Lücke und dem gleichzeitigen Verschwinden des magnetischen Kontrastes bei lokalen Schichtdicken von dCr = 4nm. Dies kann durch eine Rotation des Spindichtewellenvektors Q erklärt werden. Für dCr < 3nm verschwindet die LDW erneut. Zusätzlich zur LDW und SDW entsteht aufgrund der unterschiedlichen Gitterparameter von Chrom und Wolfram bei lokalen Schichtdicken von dCr < 3nm eine Moiré-Überstruktur. IrTe2 ist Gegenstand zahlreicher aktueller Forschungsaktivitäten und weist eine LDM mit gleichzeitiger Transformation des atomaren Gitters auf. Ein Phasenübergang erster Ordnung erzeugt zunächst bei der Übergangstemperatur TC = 275K eine Modulation mit dem Wellenvektor q = 1/5(1, 1, 0). Mithilfe temperaturabhängiger RTM-Messungen kann das Phasendiagramm um einen weiteren Übergang erster Ordnung bei TS = 180K erweitert werden. Dabei bilden sich zunehmend Te-Dimere an der sichtbaren (001)-Oberfläche und IrTe2 wechselt in einen Grundzustand mit maximaler Dichte von Dimeren und dem Wellenvektor q = 1/6(1, 1, 0). Der Mechanismus beider Phasenübergänge wird durch die Probenqualität und die Oberflächenpräparation beeinflusst, sodass die Phasenübergänge erster Ordnung teilweise verlangsamt ablaufen. Durch eine Analyse der Oberflächendynamik am Phasenübergang kann der zugrundeliegende Mechanismus des Domänenwachstums im Realraum untersucht werden. Im letzten Teil der Arbeit werden ultradünne Eisenfilme auf Rhodium(001) betrachtet. Dabei treten auf der Doppellage Eisen (Fe) auf Rhodium (Rh) spannungsabhängige elektronische Modulationen mit senkrecht zueinander orientierten Wellenvektoren q1 = [(0, 30 ± 0, 03), 0, 0] und q2 = [0, (0, 30 ± 0, 03), 0] in Richtung [100] und [010] auf. Temperaturabhängige Messungen zeigen die stetige Verkleinerung der Modulation beim Erwärmen der Probe und somit einen Phasenübergang zweiter Ordnung. Die LDM tritt auch auf der dritten und vierten Lage Eisen mit gleichgerichteten aber kleineren Wellenvektoren q auf. Spinpolarisierte RTM-Daten zeigen einen c(2×2)-Antiferromagnetismus auf einer Monolage Eisen. Für Fe-Bedeckungen von 1ML - 5ML tritt Ferromagnetismus perpendikular zur Oberfläche auf. Diese Messungen zeigen erstmals gleichzeitiges Auftreten einer elektronischen und magnetischen Phase in einem reinen 3d-Übergangsmetall im Realraum.

In this thesis, I present a model system for carbohydrate interactions with single-crystalline Ru surfaces. Geometric and electronic properties of copper phthalocyanine (CuPc) on top of graphene on hexagonal Ru(0001), rectangular Ru(10-10) and vicinal Ru(1,1,-2,10) surfaces have been studied. First, the Fermi surfaces and band structures of the three Ru surfaces were investigated by high-resolution angle-resolved photoemission spectroscopy. The experimental data and theoretical calculations allow to derive detailed information about the momentum-resolved electronic structure. The results can be used as a reference to understand the chemical and catalytic properties of Ru surfaces. Second, graphene layers were prepared on the three different Ru surfaces. Using low-energy electron diffraction and scanning tunneling microscopy, it was found that graphene can be grown in well-ordered structures on all three surfaces, hexagonal Ru(0001), rectangular Ru(10-10) and vicinal Ru(1,1,-2,10), although they have different surface symmetries. Evidence for a strong interaction between graphene and Ru surfaces is a 1.3-1.7e V increase in the graphene pi-bands binding energy with respect to free-standing graphene sheets. This energy variation is due to the hybridization between the graphene pi bands and the Ru 4d electrons, while the lattice mismatch does not play an important role in the bonding between graphene and Ru surfaces. Finally, the geometric and electronic structures of CuPc on Ru(10-10), graphene/Ru(10-10), and graphene/Ru(0001) have been studied in detail. CuPc molecules can be grown well-ordered on Ru(10-10) but not on Ru(0001). The growth of CuPc on graphene/Ru(10-10) and Ru(0001) is dominated by the Moire pattern of graphene. CuPc molecules form well-ordered structures with rectangular unit cells on graphene/Ru(10-10) and Ru(0001). The distance of adjacent CuPc molecules is 1.5 and 1.3 nm on graphene/Ru(0001) and 1.54 and 1.37 nm on graphene/Ru(10-10). This indicates that the molecule-substrate interaction dominates over the intermolecular interaction for CuPc molecules on graphene/Ru(10-10) and graphene/Ru(0001).