Energetics and Hierarchical Interactions of Metal–Phthalocyanines

Jul 23, 2013 - The adsorption of metal–phthalocyanine (MPc) layers (M = Fe, Co, Cu) assembled on graphene/Ir(111) is studied by means of temperature...
1 downloads 15 Views 1MB Size
Subscriber access provided by UNIV OF YORK

Article

Energetics and Hierarchical Interactions of MetalPhthalocyanine Adsorbed on Graphene/Ir(111) Mattia Scardamaglia, Claudia Struzzi, Silvano Lizzit, Matteo Dalmiglio, Paolo Lacovig, Alessandro Baraldi, Carlo Mariani, and Maria Grazia Betti Langmuir, Just Accepted Manuscript • DOI: 10.1021/la401850v • Publication Date (Web): 23 Jul 2013 Downloaded from http://pubs.acs.org on July 29, 2013

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Langmuir is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Energetics and Hierarchical Interactions of Metal-Phthalocyanines Adsorbed on Graphene/Ir(111) Mattia Scardamaglia,† Claudia Struzzi,‡ Silvano Lizzit,‡ Matteo Dalmiglio,‡ Paolo Lacovig,‡ Alessandro Baraldi,¶,§ Carlo Mariani,∗,† and Maria Grazia Betti∥ Dipartimento di Fisica, CNISM, CNIS, Università di Roma La Sapienza, Piazzale Aldo Moro 2, I - 00185 Roma, Italy, Elettra - Sincrotrone Trieste, Str. St. 14 km 163.5, I-34149 Trieste, Italy, Physics Department and Center of Excellence for Nanostructured Materials, University of Trieste, Via Valerio 2, I-34127 Trieste, Italy, IOM-CNR Laboratorio TASC, AREA Science Park, S.S.14 Km 163.5, I-34012 Trieste, Italy, and Dipartimento di Fisica, Università di Roma La Sapienza, Piazzale Aldo Moro 2, I - 00185 Roma, Italy E-mail: [email protected]

Abstract The adsorption of Metal-Phthalocyanine (MPc) layers (M = Fe, Co, Cu) assembled on graphene/Ir(111) is studied by means of temperature programmed X-ray Photoemission Spec∗ To

whom correspondence should be addressed di Roma La Sapienza, CNISM, CNIS ‡ Sincrotrone Trieste ¶ Università di Trieste § IOM-CNR, Trieste ∥ Università di Roma La Sapienza † Università

1 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

troscopy (XPS) and near-edge x-ray absorption fine structure (NEXAFS). The balance between interaction forces among the organometallic molecules and the underlying graphene gives rise to flat-lying molecular layers, weakly interacting with the underlying graphene. Further MPc layers pile up face-on onto the first layer, up to a few nm thickness, as deduced by NEXAFS. The FePc, CoPc and CuPc multilayers present comparable desorption temperatures, compatible with molecule-molecule interactions dominated by van der Waals forces between the π -conjugated macrocycles. The MPc single layers desorb from graphene/Ir at higher temperatures. The CuPc single-layer desorbs at lower temperature than the FePc and CoPc single-layers, suggesting a higher adsorption energy of the FePc and CoPc single-layers on graphene/Ir with respect to CuPc, with increasing molecule-substrate interaction in the order ECuPc < EFePc ∼ ECoPc .

Introduction Aromatic oligomers can be used as building blocks for nanometer-sized devices on graphene, exploiting the transport properties of π -conjugation of the macrocycle electrons, by tailoring the band structure, by doping and/or charge transfer, for potential nanoelectronic applications. 1–4 Organic molecule deposition on patterned surfaces is a viable method for inducing self-assembling and for tuning the charge transport, depending on the specific molecule-substrate interaction. MetalPhthalocyanines (M-C32 H16 N8 , MPcs) are square-shaped planar molecules with a central metal atom, four pyrrole rings and four benzene rings, 5 with the π -orbitals extended on the macrocycles ensuring conduction on substrates. 6 MPcs, embedding a metal atom in the organic cage, also constitute an useful shuttle for organizing individual metal/magnetic units when they are adsorbed on crystalline surfaces. 7–12 Planar arrays of MPcs with different motifs can actually form on the graphene moiré pattern, 13–18 as it has been recently shown by scanning tunneling microscopy measurements. 19–22 In fact, graphene acts as a buffer layer decoupling the metal substrate from the molecular layer, 23 inducing different configurations depending on the molecule-graphene interaction. 2 ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

An important parameter determining the specific molecular adsorption is the interaction strength between the molecule and graphene. A high graphene moiré corrugation can drive long-range ordered superstructures of MPcs, like the kagome super-lattice observed for FePc assembled on graphene prepared on Ru. 19,21 On the other hand, when the moiré corrugation plays a minor role, the molecule-molecule interaction dominates over the molecule-graphene interaction, like in the case of MPcs grown on graphene prepared on Ir, where square closely-packed long-range ordered molecular structures are formed. 21,22 We have recently demonstrated that FePc grown on graphene/Ir(111) lies flat from the single-layer (SL) 24 up to a multilayer coverage, 23 as determined by near-edge X-ray absorption fine-structure spectroscopy (NEXAFS). The FePc molecules of the SL arrange in a planar molecular structure, without modifying the substrate periodicity, and induce a slight electron doping to the Ir-supported graphene, while at the same time graphene acts as a buffer layer efficiently decoupling the molecules from the metallic substrate. 23 A strong mixing of the d molecular orbitals with metallic substrate have been recently reported for the adsorption of FePc and CoPc on metallic substrates, due to a rehybridization of the electronic states localized on the central atoms, breaking the 4-fold symmetry of the molecular orbitals of the MPc molecules. The molecular adsorption is controlled by a symmetry-determined mixing between the electronic states of the MPc metal center and of the metal substrate, as recently observed for MPcs adsorbed on Au(110) 8,25 Despite the weak influence of these π -conjugated molecules on the unperturbed graphene electronic properties, different central metal atoms in the MPc organic cage result in different adsorption energies on the same substrate. The adsorption energy of a single molecular layer on graphene is an important parameter for determining the actual molecule-graphene interaction, and temperature-programmed desorption experiments represent an useful way to investigate the interaction strength between a molecule and a substrate. It is possible to follow in-situ and in real-time the evolution of the molecular layer adsorbed on the substrate, by performing Temperature-Programmed XPS 26 measurements on the relevant core-levels. In this work, we present a systematic fast-XPS study of three MPcs adsorbed on graphene/Ir(111): FePc, CoPc and CuPc. These three molecules differ for the occupation of the

3 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

molecular orbitals involving the central metal states, thus for the symmetry of the partially filled orbitals involved in the interaction with the substrate. Graphene grown on the Ir surface presents a lower interaction than on other surfaces, 27 thus the graphene/Ir substrate can act as an efficient buffer layer for molecular adsorption. Nevertheless, we bring to light a higher interaction of the first FePc and CoPc layers in contact with graphene/Ir than the first CuPc SL, thus suggesting a role of the metal-derived states reinforcing the π -π van der Waals interaction, in the binding to graphene.

Experimental Synchrotron-radiation XPS data have been measured at the SuperESCA beamline of the Elettra synchrotron radiation facility in Trieste, Italy. 28 Photon energy of 130 eV, 500 eV and 400 eV were used for the valence band (VB), the N 1s and the C 1s core-levels, respectively. Photoelectrons were collected at normal emission, with a Phoibos electron energy analyzer equipped with a home-made delay line detection system. Photoemission data have been normalized to the beam intensity and the BE scale was calibrated using the Fermi edge of the Ir substrate. NEXAFS measurements across the C and N K-edges were taken at SuperESCA in Auger yield mode, with linearly polarized radiation (horizontal), with the electric field vector either parallel (incidence angle γ =0◦ , transverse electric field) or almost normal (γ =70◦ , transverse magnetic field) to the surface plane, by appropriately rotating the sample. The base pressure in the ultra-high-vacuum (UHV) chamber of the experimental station was better than 1×10−10 mbar. The UHV chamber is equipped with all necessary ancillary facilities for sample cleaning and characterization. The Ir(111) surface was cleaned by cycles of sputtering (Ar+ at E = 2 keV and T = 300 K), followed by annealing and oxygen/hydrogen-exposures cycles, to quench any carbon contamination. The graphene single-layer was prepared by repeated cycles of C2 H4 dosing followed by annealing up to 1300 K. With this procedure, we achieved a long-range ordered layer, 13–15 whose quality and cleanness were checked by means of XPS and

4 ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

low-energy electron diffraction (LEED). The MPc molecules were sublimated in UHV by means of home-made resistively heated quartz crucibles. The nominal thickness of the molecular layer was measured by using an oscillating quartz microbalance, and a single-layer is defined as a flat and compact molecular layer, 29 corresponding to about 3.4 Å nominal coverage. During MPc evaporation the sample was kept at room temperature (RT) and the overall pressure during deposition was in the low 10−9 mbar range. Fast-XPS measurements were taken after deposition of a few nm-thick MPc molecular layer. We linearly increased the sample temperature with a rate of 0.3 K/s, following in-situ and in realtime the desorption process. With this procedure, we followed the temperature evolution of the molecular layers adsorbed on the substrate, by analyzing the N 1s core-level intensity, binding energy and lineshape. The N 1s core-level represents a clear signature of the MPc molecule, because the molecular C 1s component is difficult to deconvolve from the graphene signal. 24 Finally, we estimate the desorption energy of the molecular film deposited on the graphene/Ir surface, by appropriately treating the data within the Redhead equation. 30

Results and Discussion The Metal-Phthalocyanine layers have been grown on freshly prepared graphene on Ir(111), by subsequent deposition steps at room-temperature. The adsorption geometry of CoPc and CuPc layers adsorbed on graphene can be investigated by the dichroic response of the C and N absorption edges, as a function of the incident linearly polarized radiation. The N K-edge near-edge absorption spectra for the CoPc and CuPc SL and thin-films prepared on graphene/Ir(111) are shown in Fig. Figure 1. Data show a definite dichroic response for both systems, due to the welldefined molecular orientation of both MPcs. In particular, the N K-edge spectra for both systems presents a first resonance at about 398.4 eV photon energy, constituted by two close structures due to photoexcitation from the 1s level of the two nitrogen atoms in the pyrrole ring to the lowestunoccupied molecular-orbital (LUMO) of π ∗ -symmetry. The higher-energy structures at 400.6,

5 ACS Paragon Plus Environment

Langmuir

402.5, and 404.0 eV photon energies, respectively, are due to transitions to higher-lying unoccupied states mainly involving the pyrrole ring, 31 while the resonances above 406 eV photon energy are associated to transition to unoccupied orbitals with σ ∗ symmetry, in agreement with the attribution given for FePc and CoPc planar thin-film (TF) grown on other crystalline substrates. 9 The clear dichroic response of the signal is fingerprint of flat-lying molecules, from the SL up to the formation of a thin-film on the graphene-Ir substrate. Furthermore, the absorption resonances maintain the same energy position and lineshape in the different coverage phases, indicating a negligible interaction of the molecules with graphene and confirming the decoupling role played by graphene with respect to the bare Ir surface. 23

Intensity (arb. units)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 27

N-K edge

N-K edge

TF CoPc

TF CuPc

E

g

SL CoPc

400

405

410

2SL CuPc

400

405

410

Photon Energy (eV) Figure 1: NEXAFS across the N K-edge for CoPc SL and TF (left) and CuPc (≃ 2 SL and TF (right) adsorbed on graphene/Ir(111). Red curves correspond to γ =0◦ and blue curves correspond to γ =70◦ , where γ is the incidence angle between the incident radiation and the surface normal, as reported in the sketch in the inset. The NEXAFS C K-edge data at the SL coverage for CoPc and CuPc on graphene/Ir are shown 6 ACS Paragon Plus Environment

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Intensity (arb. units)

Page 7 of 27

C-K edge

C-K edge

SL CoPc

2SL CuPc

E

g

graphene/Ir

285

290

295

graphene/Ir

300

285

290

295

300

Photon Energy (eV) Figure 2: NEXAFS across the C K-edge for a CoPc SL (left) and 2 SL of CuPc (right) adsorbed on graphene/Ir(111), and for graphene/Ir(111) for comparison (bottom curves). Red curves correspond to γ =0◦ and blue curves correspond to γ =70◦ , where γ is the incidence angle between the incident radiation and the surface normal, as reported in the sketch in the inset.

7 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

in Fig. Figure 2, compared with the C K-edge from clean graphene/Ir. The absorption structures in clean graphene are fingerprint of formation of a good single graphene sheet. 23,32,33 The adsorption of the first CoPc and CuPc layer determines the emergence of new C-related features characteristic of metal-phthalocyanines, 23 superimposed onto the intense graphene signal, maintaining the dichroic signal among the two polarizations, characteristic of flat layers. These NEXAFS results on CoPc and CuPc grown on graphene/Ir, are analogous to the data recently taken and discussed for FePc grown on graphene/Ir, building-up planar layers from the SL to a multi-layer coverage, while maintaining unperturbed resonance lineshape. 23 Thus, MPc growth on graphene/Ir is characterized by a flat molecular orientation and weak interaction of the organic macrocycles with the underlying graphene. The N 1s core-levels for CoPc and CuPc thin films grown on graphene/Ir, taken at room temperature, are reported in the middle panels of Fig. Figure 3. Their lineshape is due to the presence of two non-equivalent nitrogen atoms in the phthalocyanine molecule, the isoindole nitrogen linked to the central atom, and the azomethine nitrogen bridging two carbon atoms in the pyrrole rings, respectively, as shown in the fit superimposed to the experimental data, in agreement with previous determinations in MPc thin-films. 6,34,35 The slightly different N 1s binding energy measured for CoPc (5 nm-thick film) and CuPc (2.4 nm-thick film) is due to a small band bending variation in the semiconducting molecular films, depending on the TF thickness. In order to follow the N 1s core-level intensity evolution as a function of temperature, we fit the data in the following as a single peak. 36 The intensity evolution of the N 1s core levels as a function of temperature is shown through two-dimensional (2D) projected plots of the fast-XPS spectra, as shown in the upper panels of Fig. Figure 3. The single spectra of the time-lapsed sequence have been fitted with single Gaussian curves and the results of the fitting procedure, full-width at half-maximum (FWHM) and binding energy (BE) for both molecular films, are reported in the bottom panels of Fig. Figure 3. The N 1s peak intensity stays constant up to about 590 K for CoPc and 570 K for CuPc on graphene/Ir. Above these temperatures, we observe a drop of the intensity, a shift to lower binding energy of the N 1s peak position, and no effects on the peak width. We attribute this intensity drop

8 ACS Paragon Plus Environment

Page 8 of 27

Page 9 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 3: N 1s core-level of CoPc (left) and CuPc (right) thin-films adsorbed on graphene/Ir(111). Top panels: fast-XPS data as 2D trace of the spectral intensity as a function of BE and temperature (higher to lower intensity signal represented in false colors from red to blue). Middle panels: N 1s core-level data (dotted lines) taken at the initial temperature, RT; experimental data fitted with two 9 N atoms. Bottom panels: BE (blue squares) and components associated to the twoACS non-equivalent Paragon Plus Environment FWHM (green squares) evolution of the N 1s core level as a function of temperature.

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 27

to the desorption of the flat molecular thin-film adsorbed onto the graphene/Ir(111) surface. The desorption temperature for these thin-films is comparable among different MPcs and basically due to the weak intermolecular forces. An analogous temperature was found for NiPc TF desorption from highly-oriented pyrolitic graphite (HOPG), 37 and for FePc TF desorption from from a variety of substrates, from Au to HOPG. 38 It is thus compatible with the multilayer desorption of physisorbed molecular MPc multi-layers dominated by van der Waals interactions. Increasing the temperature above 620 K, the N 1s core-level intensity for the MPcs reaches a plateau, while the BE reaches the same value obtained for a single-layer as-deposited on graphene/Ir at RT. The desorption of the MPc thin films leaves a stable SL adsorbed on the graphene/Ir surface. 23 We report in Fig. Figure 4 the fast-XPS data for FePc (left), CoPc (center) and CuPc (right) on graphene/Ir, as a function of temperature, in order to determine the desorption temperature of the first molecular layer in direct contact with the graphene/Ir substrate. The CuPc SL (right panel) presents a fast and total desorption at about 670 K, corresponding to the disappearance of the N 1s peak. The N 1s core-level shape evolution for the CuPc SL is consistent with the full desorption of the single-layer in contact with graphene. It has be fitted with a single component 36 which starts to decrease above 670 K temperature, that we consider as the desorption temperature. The N 1s peak lineshape evolution for FePc and CoPc as a function of temperature shows an intensity decrease above 820 K, taken as desorption temperature for the SLs, while maintaining a lower decreasing rate with respect to the case of CuPc. We notice that the molecular layer resulting from annealing a multi-layer at 775 K, just below the desorption temperature of the SL, presents the same spectroscopic signatures (binding energy, absorption resonances and spectral density of states in the valence band) of the SL as deposited at RT (data presented in the Supplementary Information section). After the first desorption stage, the N 1s cannot be fitted with a single component above 900 K and it has been fitted with two components, as reported in Figure 5 and Figure 6 for FePc and CoPc, respectively. In fact, there is an apparent shift of the N 1s core-level to lower binding energies actually due to the growth of a second N 1s component at lower BE for both FePc and

10 ACS Paragon Plus Environment

Page 11 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 4: Top panels: 2D plot of the fast-XPS data of the N 1s core-level of FePc (left), CoPc (center) and CuPc (right) single-layers adsorbed on graphene/Ir(111). Higher to lower intensity signal represented in false colors from red to blue. Bottom panels: BE (blue squares) and FWHM (green squares) evolution of the N 1s core level as a function of temperature.

11 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 27

CoPc systems, which becomes dominant at very high temperature (above 1100 K). This lower BE component appears above the SL desorption temperature, it grows in intensity till the highest temperature reached in the experiment, and is present only at high temperatures and only under high-brilliance and focussed x-ray beam. Above 900 K, also the spectroscopic signature of the C 1s core-level presents a modified shape with respect to the typical graphene layer (see Supplementary Information). The lineshape modification of C 1s and the emerging of the new N 1s component can be attributed to doping of the graphene layer with pyridinic nitrogen 39,40 deriving from the FePc ad CoPc molecular fragmentation. The determination of the desoption temperature is independent from the presence of the doping N defects, as explained in the Supplementary information section, and the majority of molecules of the SLs desorb in the temperature range where the radiationinduced molecular fragmentation is not yet efficient, while at higher temperature the cooperation of temperature and radiation damage induces the uptake of doping N. The undesired effect of beam damage, introducing N defects in graphene, may be exploited as a possible method for doping and patterning, once performed in a controlled way. The desorption temperatures in each phase, shown in Table Table 1, can be correlated to the activation energies for desorption, thus to the adsorption energy on the surface. We can estimate the activation energy EA for the desorption of the MPc TF and of the MPc SL on graphene/Ir(111), by using the Redhead equation: 30 ] [ ν Tmax EA EA = KB Tmax ln − ln ; β KB Tmax

(1)

where KB is the Boltzmann constant, ν a frequency factor, Tmax the desorption temperature,

β the heating rate. In our experiment, the heating rate β is 0.3 K/s, the quantity ln(EA /KB Tmax ) is a small constant value compared to the first term in brackets. These approximations are valid for first order kinetic processes that are non dissociative desorptions without multiple process of re-adsorption. In these processes, the desorption temperature does not change with the initial layer thickness, as in the present case, discussed in the Supplementary Information section. The frequency factor for large molecules with many internal degrees of freedom can be much higher 12 ACS Paragon Plus Environment

Page 13 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 5: Fast-XPS N 1s core-level spectra of FePc on graphene/Ir(111), as a function of temperature. Left panel: single energy distribution curve data, vertically stacked for clarity. Central panel: fit of significant spectra at different temperatures. Right panel: intensity evolution of the two N 1s components as a function of temperature (molecular N, red lines, doping N, blue lines); N 1s peak intensity normalized to the SL intensity.

13 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6: Fast-XPS N 1s core-level spectra of CoPc on graphene/Ir(111), as a function of temperature. Left panel: single energy distribution curve data, vertically stacked for clarity. Central panel: fit of significant spectra at different temperatures. Right panel: intensity evolution of the two N 1s components as a function of temperature (molecular N, red lines, intercalated N, blue lines); N 1s peak intensity normalized to the SL intensity.

14 ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

than in small molecules, 41 and we can assume ν to be 1018 s−1 . The choice of this frequency factor has been taken in order to obtain an activation energy for the CuPc thin-film grown on graphene-Ir in agreement with the measured sublimation enthalpy of CuPc thin-films. 42 With this assumption on the frequency factor, we estimated activation energies for the desorption of the thin-films and of the SL of three MPcs adsorbed on graphene/Ir, resumed in Table Table 1. The thin-films of all MPc molecular species present very close activation energies of about 2.2-2.4 eV, compatible with the van der Waals interaction between the π -organic macrocycles of the metalphthalocyanines 42,43 producing the molecular stacking into molecular planes, and common to all three molecules, basically independent on the specific central metal atom. Table 1: Activation energies for desorption, calculated with the Redhead equation (1) for the adsorption of FePc, CoPc and CuPc on graphene prepared on Ir(111): thin-film and single-layer activation energies. Temperature error ±15 K, energy error ±0.1 eV. (a) from Ref. [ 38 ]. TF SL

FePc 2.4 eV (600 K(a) ) 3.2 eV (820 K)

CoPc CuPc 2.3 eV (590 K) 2.2 eV (570 K) 3.2 eV (820 K) 2.6 eV (670 K)

On the other hand, the activation energy of the SL on graphene/Ir is different among the three molecules, i.e. 2.6 eV for CuPc, and 3.2 eV for FePc and CoPc. The activation energy estimated for of the CuPc SL (2.6 eV) is comparable to that of similar organic oligomer adsorbed on metal surfaces 44,45 or on organic buffer layer onto a metal, 46 and it is consistent with an estimation of activation energy basically exerted between aromatic species, driven by π − π interaction. 47 The activation energy of the FePc and CoPc SLs (3.2 eV), instead, presents a stronger interaction than the CuPc SL, and is dominated by stronger molecule-substrate forces. Comparably adsorption energies are found in molecular systems on surfaces with stronger interaction due to a particular surface morphology like vicinal surfaces, 44,48 where the surface step morphology creates preferentially and energetically favorable adsorption site, or for different large organic molecules adsorbed on graphene/Ru(0001), 49 where the lateral variations in the molecule-substrate interactions play an important role. The frequency factor ν for large molecules can also change due to different growth morphology (planar, up-right geometry). 49 Even considering a modification by two orders 15 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of magnitude of the pre-factor, associated to a freezing of the modes, in contact with the substrate, the absolute energy value would change by less than 10% However, it does not affect the relative energy difference between the various MPcs studied here in each coverage phase, because the three MPc molecular layers on graphene/Ir maintain the planar geometry from sub-SL to the multi-layer coverage, and they present the same organic cage embedding the central metal atom. The changes in the absolute energy estimation due to the specific choice of a frequency factor, do not modify the definite energy difference observed among FePc and CuPc SL on graphene/Ir. We have assumed the same frequency pre-factor in the Redhead equation to evaluate the adsorption energy for the multilayer and the single layers. The MPc layers, in contact with the graphene, can have degree of freedom frozen by the adsorption geometry and this can induce increasing values for the frequency factor. Even assuming a frequency factor by orders of magnitude higher than the average value used for this estimation, which affects the absolute energy, the relative increase in activation energy between FePc and CoPc, and CuPc SLs, is maintained. In general, for aromatic molecules adsorbed on graphene, both graphene-like (CC) and benzene-like (CH) carbon atoms contribute to the bonding, which also depends on the competition between dipolar, nonlocal correlation and electrostatic energy contributions to the total energy, varying as a function of the molecule-to-graphene distance. 47 Recently, the lowest-energy linear assembly of coplanar unsupported FePc chains has been predicted by means of density functional theory (DFT) calculations by computing the interaction energy as a function of the intermolecular distance, and stable long-range ordered chain configurations have been obtained, without incurring any energy barrier. 50 The intermolecular attractive electrostatic potential is the major driving force to stabilize the molecular assembly, and the molecule-molecule distance is dictated by the repulsion barrier. Long-range ordered configurations, mainly driven by the MPc-MPc mutual electrostatic interactions, have been observed for CoPc assembled on graphene/Ir(111) 22 and graphite. 51 State-of-the-art first principle DFT calculations evaluate a negligible interaction of MPc with graphene if van der Waals forces are not included. 51,52 The inclusion of van der Waals functional in the DFT calculations of CuPc on

16 ACS Paragon Plus Environment

Page 16 of 27

Page 17 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

graphene induces an energy minimum for flat-lying molecules with a binding energy of 3.37 eV. 53 Very recent preliminary calculations within the DFT approach, including different van der Waals functionals, obtain basically undistorted flat-lying MPc molecules on graphene, whose minimum energy configuration is a hollow site for the central metal atom on the carbon network, with a weak dependence on a specific adsorption site, 52 and where the main contribution to the adsorption energy is due to the organic macrocycles without the metallic center (2 eV). 52 These considerations can rationalize the formation of stable ordered flat MPc layers on graphene. As far as it concerns the role of the central metal atom in the adsorption on graphene, there is not yet any theoretical analysis enlightening a clear hybridization and/or interaction process explaining the observed adsorption energies. On the other hand, theoretical and experimental data have rationalized the interaction of MPc on metallic substrates, where the energy landscape is strongly influenced from the MPc central metal atoms orbitals mixed and hybridized with the underlying metallic states, as for example on Au and Ag surfaces. 50,54,55 Previous experimental analysis of FePc adsorption on graphene-Ir enlightened the absence of molecular orbital occupation, hybridization or significant charge transfer between graphene and the molecules. 23 Thus, we suggest that very slight deformation of the molecules with subsequent broken symmetry of the crystal field might play a role in the case of FePc and CoPc on graphene-Ir, claiming for refined theoretical analysis.

Conclusions Temperature-Programmed XPS was used to follow the evolution of MPc (M = Fe, Co, Cu) molecular layers adsorbed onto the graphene/Ir(111) surface to determine the desorption temperatures. Different phases have been singled out, by following the N 1s core-level intensity, binding-energy and lineshape evolution, upon linearly increasing the temperature, and estimating the adsorption energy. Very close values of activation energy (≃ 2.2-2.4 eV) is determined for all MPcs for the desorption of the flat multi-layer film, while a higher activation energy is obtained for the SL of FePc and CoPc with graphene/Ir (3.2 eV), with respect to SL of CuPc (2.6 eV). The flat-lying MPc

17 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 27

single-layers are stabilized by cooperative molecule-molecule interactions, van der Waals forces and dipolar interaction with the underlying graphene. The higher interaction strength experienced by FePc and CoPc SLs with graphene/Ir can be related to tiny molecular deformation and crystal field effects, also depending on the metal-associated d-orbital occupancy. The energetics and the role of competing interaction in graphene functionalized with organic molecules is a key challenge in graphene research, to move from model structures to more complex systems. A comprehensive theoretical contribution in the determination of adsorption energy in van der Waals systems with competing and hierarchical interactions is demanding, but highly wished.

Acknowledgements We kindly thank the Elettra synchrotron radiation facility for the experimental assistance. Marco Buongiorno Nardelli, Arrigo Calzolari and Andrea Pelissetto are gratefully acknowledged for useful scientific discussions. Work funded by PRIN grant 20105ZZTSE "GRAF" of the Italian Ministery for Research (MIUR), by Roma "La Sapienza" University funds, and by the Eletra user support.

Supporting Information Available Valence band, core-level photoemission and NEXAFS data for the FePc SL as-deposited on graphene/Ir(111) at RT, and as obtained after a thermal treatment of a FePc multi-layer, together with core-level data taken at very high temperatures, and temperature programmed desorption mass spectra, are available in the SI. This information is available free of charge via the Internet at http://pubs.acs.org.

References 1. Huang, H.; Chen, S.; Gao, X.; Chen, W.; Wee, A. T. S. Structural and Electronic Properties of PTCDA Thin Films on Epitaxial Graphene. ACS Nano 2009, 3, 3431–3436.

18 ACS Paragon Plus Environment

Page 19 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

2. Wang, X.; Xu, J.-B.; Xie, W.; Du, J. Quantitative Analysis of Graphene Doping by Organic Molecular Charge Transfer. The Journal of Physical Chemistry C 2011, 115, 7596–7602. 3. Zhang, W.; Lin, C.-T.; Liu, K.-K.; Tite, T.; Su, C.-Y.; Chang, C.-H.; Lee, Y.-H.; Chu, C.-W.; Wei, K.-H.; Kuo, J.-L. et al. Opening an electrical band gap of bilayer graphene with molecular doping. ACS Nano 2011, 5, 7517–24. 4. Mao, H. Y.; Wang, Y.; Niu, T. C.; Zhon, J. Q.; Huang, M. Y.; Qi, D. C.; Loh, K. P.; Wee, A. T. S.; Chen, W. Chemical vapor deposition graphene as structural template to control interfacial molecular orientation of chloroaluminium phthalocyanine. Applied Physics Letters 2011, 99, 093301. 5. Forrest, S. R. Ultrathin Organic Films Grown by Organic Molecular Beam Deposition and Related Techniques. Chemical Reviews 1997, 97, 1793 – 1896. 6. Evangelista, F.; Ruocco, A.; Gotter, R.; Cossaro, A.; Floreano, L.; Morgante, A.; Crispoldi, F.; Betti, M. G.; Mariani, C. Electronic states of CuPc chains on the Au(110) surface. The Journal of Chemical Physics 2009, 131, 174710. 7. Cossaro, A.; Cvetko, D.; Bavdek, G.; Floreano, L.; Gotter, R.; Morgante, A. CopperPhthalocyanine Induced Reconstruction of Au (110). Journal of Physical Chemistry B 2004, 108, 14671–14676. 8. Gargiani, P.; Angelucci, M.; Mariani, C.; Betti, M. G. Metal-phthalocyanine chains on the Au(110) surface : Interaction states versus d-metal states occupancy. Physical Review B 2010, 81, 085412. 9. Betti, M. G.; Gargiani, P.; Frisenda, R.; Biagi, R.; Cossaro, A.; Verdini, A.; Floreano, L.; Mariani, C. Localized and dispersive electronic states at ordered FePc and CoPc chains on Au(110). Journal of Physical Chemistry C 2010, 114, 21638–21644.

19 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

10. Kröger, I.; Stadtmüller, B.; Stadler, C.; Ziroff, J.; Kochler, M.; Stahl, A.; Pollinger, F.; Lee, T.L.; Zegenhagen, J.; Reinert, F. et al. Submonolayer growth of copper-phthalocyanine on Ag(111). New Journal of Physics 2010, 12, 083038. 11. Annese, E.; Fujii, J.; Vobornik, I.; Rossi, G. Structure and Electron States of Cophthalocyanine Interacting With the Cu(111) Surface. The Journal of Physical Chemistry C 2011, 115, 17409–17416. 12. Petraki, F.; Peisert, H.; Latteyer, F.; Ayg, U.; Vollmer, A.; Chass, T. Impact of the 3d Electronic States of Cobalt and Manganese Phthalocyanines on the Electronic Structure at the Interface to Ag(111). The Journal of Physical Chemistry C 2011, 115, 21334–21340. 13. Lacovig, P.; Pozzo, M.; Alfè, D.; Vilmercati, P.; Baraldi, A.; Lizzit, S. Growth of DomeShaped Carbon Nanoislands on Ir(111): The Intermediate between Carbidic Clusters and Quasi-Free-Standing Graphene. Physical Review Letters 2009, 103, 166101. 14. N’Diaye, A. T.; Coraux, J.; Plasa, T. N.; Busse, C.; Michely, T. Structure of epitaxial graphene on Ir(111). New Journal of Physics 2008, 10, 043033. 15. Pletikosi´c, I.; Kralj, M.; Pervan, P.; Brako, R.; Coraux, J.; Diaye, A. T. N.; Busse, C.; Michely, T. Dirac Cones and Minigaps for Graphene on Ir(111). Physical Review Letters 2009, 102, 056808. 16. Preobrajenski, A. B.; Ng, M. L.; Vinogradov, A. S.; Mårtensson, N. Controlling graphene corrugation on lattice-mismatched substrates. Physical Review B 2008, 78, 073401. 17. Murata, Y.; Petrova, V.; Kappes, B. B.; Ebnonnasir, A.; Petrov, I.; Xie, Y.-H.; Ciobanu, C. V.; Kodambaka, S. Moiré Superstructures of Graphene on Faceted Nickel Islands. ACS Nano 2010, 4, 6509 – 5514. 18. Merino, P.; Svec, M.; Pinardi, A. L.; Otero, G.; Martín-Gago, J. A. Strain-Driven Moiré Su-

20 ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

perstructures of Epitaxial Graphene on Transition Metal Surfaces. ACS Nano 2011, 5, 5627 – 5634. 19. Mao, J.; Zhang, H.; Jiang, Y.; Pan, Y.; Gao, M.; Xiao, W.; Gao, H.-J. Tunability of supramolecular Kagome lattices of magnetic phthalocyanines using graphene-based moire patterns as templates. Journal of the American Chemical Society 2009, 131, 14136 – 14137. 20. Zhang, H.; Sun, J.; Low, T.; Zhang, L.; Pan, Y.; Liu, Q.; Mao, J.; Zhou, H.; Guo, H.; Du, S. et al. Assembly of iron phthalocyanine and pentacene molecules on a graphene monolayer grown on Ru(0001). Physical Review B 2011, 84, 245436. 21. Yang, K.; Xiao, W.; Jiang, Y.; Zhang, H.; Liu, L.; Mao, J.; Zhou, H.; Du, S.; Gao, H.-J. Molecule-Substrate Coupling between Metal Phthalocyanines and Epitaxial Graphene Grown on Ru(0001) and Pt(111). The Journal of Physical Chemistry C 2012, 116, 14052 – 14056. 22. Hämäläinen, S. K.; Stepanova, M.; Drost, R.; Liljeroth, P.; Lahtinen, J.; Sainio, J. SelfAssembly of Cobalt-Phthalocyanine Molecules on Epitaxial Graphene on Ir(111). The Journal of Physical Chemistry C 2012, 116, 20433 – 20437. 23. Scardamaglia, M.; Lisi, S.; Lizzit, S.; Baraldi, A.; Larciprete, R.; Mariani, C.; Betti, M. G. Graphene Induced Substrate Decoupling and Ideal Doping of a Self-Assembled IronPhthalocyanine Single Layer. Journal of Physical Chemisty C 2013, 116, 3019 – 3027. 24. Scardamaglia, M.; Forte, G.; Lizzit, S.; Baraldi, A.; Lacovig, P.; Larciprete, R.; Mariani, C.; Betti, M. G. Metal-phthalocyanine array on the moiré pattern of a graphene sheet. Journal of Nanoparticle Research 2011, 13, 6013–6020. 25. Gargiani, P.; Rossi, G.; Biagi, R.; Pedio, M.; Fortuna, S.; Calzolari, A.; Fabris, S.; Cezar, J. C.; Brookes, N.; Betti, M. G. Spin and orbital configuration of Metal Phthalocyanine chains assembled on the Au(110) surface. Physical Review B 2013, 87, 165407.

21 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

26. Baraldi, A.; Comelli, G.; Lizzit, S.; Cocco, D.; Paolucci, G.; Rosei, R. Temperature programmed X-ray photoelectron spectroscopy: A new technique for the study of surface kinetics. Surface Science 1996, 367, L67–L72. 27. Batzill, M. The surface science of graphene: Metal interfaces, CVD synthesis, nanoribbons, chemical modifications, and defects. Surface Science Reports 2012, 67, 83–115. 28. Baraldi, A.; Comelli, G.; Lizzit, S.; Kiskinova, M.; Paolucci, G. Real-time X-ray photoelectron spectroscopy of surface reactions. Surface Science Reports 2003, 49, 169–224. 29. Betti, M. G.; Gargiani, P.; Mariani, C.; Biagi, R.; Fujii, J.; Rossi, G.; Resta, A.; Fabris, S.; Fortuna, S.; Torrelles, X. et al. Structural Phases of Ordered FePc-Nanochains Self-Assembled on Au(110). Langmuir 2012, 28, 13232 – 13240. 30. Redhead, P. Thermal desorption of gases. Vacuum 1962, 12, 203 – 211. 31. Calabrese, A.; Floreano, L.; Verdini, A.; Mariani, C.; Betti, M. G. Filling empty states in a CuPc single layer on the Au(110) surface via electron injection. Physical Review B 2009, 79, 115446. 32. Pacilé, D.; Papagno, M.; Rodríguez, A. F.; Grioni, M.; Papagno, L. Near-Edge X-Ray Absorption Fine-Structure Investigation of Graphene. Physical Review Letters 2008, 101, 066806. 33. Papagno, M.; Rodríguez, A. F.; Girit, c. O.; Meyer, J. C.; Zettl, A.; Pacilé, D. Polarizationdependent C K near-edge X-ray absorption fine-structure of graphene. Chemical Physics Letters 2009, 475, 269–271. 34. Papageorgiou, N.; Ferro, Y.; Salomon, E.; Allouche, A.; Layet, J. M.; Giovanelli, L.; Lay, G. L. Geometry and electronic structure of lead phthalocyanine: Quantum calculations via densityfunctional theory and photoemission measurements. Physical Review B 2003, 68, 235105. 35. Papageorgiou, N.; Salomon, E.; Angot, T.; Layet, J.-M.; Giovanelli, L.; Lay, G. L. Physics

22 ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

of ultra-thin phthalocyanine films on semiconductors. Progress in Surface Science 2004, 77, 139–170. 36. In this fit, the double-component N peak associated to the intact molecule has been treated as a single peak. 37. Ottaviano, L.; Di Nardo, S.; Lozzi, L.; Passacantando, M.; Picozzi, P.; Santucci, S. Thin and ultra-thin films of nickel phthalocyanine grown on highly oriented pyrolitic graphite: An XPS, UHV-AFM and air tapping-mode AFM study. Surface Science 1997, 373, 318 – 322. 38. Struzzi, C.; Scardamaglia, M.; Angelucci, M.; Massimi, L.; Mariani, C.; Betti, M. G. Adsorption Energy of Iron-Phthalocyanine on Crystal Surfaces. Il Nuovo Cimento C 2013, 36, 51 – 56. 39. Usachov, D.; Vilkov, O.; Grüneis, A.; Haberer, D.; Fedorov, A.; Adamchuk, V.; Preobrajenski, A.; Dudin, P.; Barinov, A.; Oehzelt, M. et al. Nitrogen-Doped Graphene: Efficient Growth, Structure, and Electronic Properties. Nano Letters 2011, 11, 5401 – 5407. 40. Kim, K.-j.; Yang, S.; Park, Y.; Myungjin, L.; Kim, B.; Lee, H. Annealing Effects after Nitrogen Ion Casting on Monolayer and Multilayer Graphene. The Journal of Physical Chemistry C 2013, 117, 2129 – 2134. 41. Zacharia, R.; Ulbricht, H.; Hertel, T. Interlayer cohesive energy of graphite from thermal desorption of polyaromatic hydrocarbons. Physical Review B 2004, 69, 155406. 42. Basova, T.; Semyannikov, P.; Plyashkevich, V.; Hassan, A.; Igumenov, I. Volatile Phthalocyanines: Vapor Pressure and Thermodynamics. Critical Reviews in Solid State and Materials Sciences 2009, 34, 180–189. 43. Bonderman, D. P.; Cater, E. D.; Bennett, W. E. Vapor pressures, mass spectra, magnetic susceptibilities, and thermodynamics of some phthalocyanine compounds. Journal of Chemical Engineering Data 1970, 15, 396–400. 23 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

44. Baldacchini, C.; Mariani, C.; Betti, M. G. Adsorption of pentacene on filled d-band metal surfaces: long-range ordering and adsorption energy. The Journal of Chemical Physics 2006, 124, 154702. 45. Kanjilal, A.; Ottaviano, L.; Di Castro, V.; Beccari, M.; Betti, M. G.; Mariani, C. Pentacene Grown on Self-Assembled Monolayer: Adsorption Energy, Interface Dipole, and Electronic Properties. Journal of Physical Chemistry C 2007, 111, 286 – 293. 46. Schuerlein, T.; R.M., A. Formation and characterization of epitaxial phthalocyanine and perylene monolayers and bilayers on Cu(100): Low-energy electron diffraction and thermal desorption mass spectrometry studies. Journal of Vacuum Science and Technology A 1994, 12, 1992 – 1997. 47. Björk, J.; Hanke, F.; Palma, C.-A.; Samori, P.; Cecchini, M.; Persson, M. Adsorption of Aromatic and Anti-Aromatic Systems on Graphene through π -π Stacking. Journal of Physical Chemistry Letters 2010, 1, 3407 – 3412. 48. Angelucci, M.; Gargiani, P.; Mariani, C.; Betti, M. G. Potassium-doped FePc thin-film on metal surfaces: observation of different empty state occupation. Journal of Nanoparticle Research 2011, 13, 5967–5973. 49. Roos, M.; Künzel, D.; Uhl, B.; Huang, H.-H.; Alves, O. B.; Hoster, H. E.; Gross, A.; Behm, R. J. Hierarchical Interactions and Their Influence upon the Adsorption of Organic Molecules on a Graphene Film. Journal of the American Chemicl Society 2011, 133, 9208 – 9211. 50. Fortuna, S.; Gargiani, P.; Betti, M. G.; Mariani, C.; Calzolari, A.; Modesti, S.; Fabris, S. Molecule-Driven Substrate Reconstruction in the Two-Dimensional Self-Organization of FePhthalocyanines on Au(110). The Journal of Physical Chemistry C 2012, 116, 6251 – 6258. 51. Scheffler, M.; Smykalla, L.; Baumann, D.; Schlegel, R.; HÃd’nke, T.; Toader, M.; Büchner, B.;

24 ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Hietschold, M.; Hess, C. Structural study of monolayer cobalt phthalocyanine adsorbed on graphite. Surface Science 2013, 608, 55 – 60. 52. Marco Buongiorno Nardelli and Arrigo Calzolari, forthcoming paper to be published. 2013, 53. Xiao, K.; Deng, W.; Keum, J. K.; Yoon, M.; Vlassiouk, I. V.; Clark, K. W.; Li, A.-P.; Kravchenko, I. I.; Gu, G.; Payzant, E. A. et al. Surface-Induced Orientation Control of CuPc Molecules for the Epitaxial Growth of Highly Ordered Organic Crystals on Graphene. Journal of the American Chemical Society 2013, 135, 3680–3687. 54. Betti, M. G.; Gargiani, P.; Mariani, C.; Turchini, S.; Zema, N.; Fortuna, S.; Calzolari, A.; Fabris, S. Formation of Hybrid Electronic States in FePc Chains Mediated by the Au(110) Surface. The Journal of Physical Chemistry C 2012, 116, 8657–8663. 55. Petraki, F.; Peisert, H.; Aygül, F., U. Latteyer; Uihlein, J.; Vollmer, A.; Chassé, T. Electronic Structure of FePc and Interface Properties on Ag(111) and Au(100). The Journal of Physical Chemistry C 2011, 116, 11110–11116.

25 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

Langmuir

ACS Paragon Plus Environment