Exciting H2 Molecules for Graphene Functionalization - ACS Nano

Dec 18, 2017 - Hydrogen functionalization of graphene by exposure to vibrationally excited H2 molecules is investigated by combined scanning tunneling...
0 downloads 6 Views 1MB Size
Article

Exciting H2 Molecules for ACS Nano is published by the American Graphene Functionalization Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American ChemicalLibrary Subscriber access provided by Gothenburg University Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced

Line Kyhl, Regis Bisson, Richard Balog, Michael N. Groves, Esben Leonhard Kolsbjerg, Andrew Martin Cassidy, Jakob Holm Jørgensen, ACS is published by the American Susanne Halkjær, Jill A.Nano Miwa, Antonija Chemical Society. 1155 Sixteenth Grubiši# #abo, Thierry Angot, Philip Hofmann, Street N.W., Washington, DC 20036 Published by American ChemicalLibrary Subscriber access provided by Gothenburg University Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced

Mohammad Alif Arman, Samuli Urpelainen, Paolo Lacovig, Luca Bignardi, Hendrik Bluhm, Jan Knudsen, Bjørk Hammer, and Liv Hornekaer ACS Nano is published by the American

ACS Nano, Just Accepted Manuscript • DOI: Chemical Society. 1155 10.1021/ Sixteenth acsnano.7b07079 • Publication Date (Web): 18 Dec 2017 Street N.W., Washington, DC 20036

Published by American ChemicalLibrary Subscriber access provided by Gothenburg University Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced

Downloaded from http:// pubs.acs.org on December 20, 2017

Just

ACS Nano is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Accepted Published by American ChemicalLibrary Subscriber access provided by Gothenburg University Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced

“Just Accepted” manuscripts have been peer-reviewed and acc online prior to technical editing, formatting for publication and a Society provides “Just Accepted” as a free service to the ACS Nano is published by the American dissemination of scientific material as soon as possible after ac Chemical Society. 1155 Sixteenth N.W., Washington, DC HTML 20036 abstract. appear in full in PDF formatStreet accompanied by an Published by American ChemicalLibrary Subscriber access provided by Gothenburg University Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced

fully peer reviewed, but should not be considered the official vers readers and citable by the Digital Object Identifier (DOI®). “Just to authors. Therefore, the “Just Accepted” Web site may not in in the journal. After a manuscript technically ACS Nanois is published by theedited Americanand form Chemical Society. Sixteenth Accepted” Web site and published as an1155 ASAP article. Note tha Street N.W., Washington, DC 20036 changes to theaccess manuscript text and/or graphics which could a Published by American ChemicalLibrary Subscriber provided by Gothenburg University Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced

and ethical guidelines that apply to the journal pertain. ACS or consequences arising from the use of information contained ACS Nano is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American ChemicalLibrary Subscriber access provided by Gothenburg University Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced

Page 1 of 25 1 2 3 4 5 6 7 8

ACS Nano

ACS Paragon Plus Environment

20 Å

ACS Nano 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 2 of 25

Exciting H2 Molecules for Graphene Functionalization Line Kyhl,a Régis Bisson,b Richard Balog,c Michael N. Groves,c Esben Leonhard Kolsbjerg,a Andrew Martin Cassidy,c Jakob Holm Jørgensen,a Susanne Halkjær,c Jill A. Miwa,a,c Antonija Grubišić Čabo,a Thierry Angot,b Philip Hofmann,c Mohammad Alif Arman,d Samuli Urpelainen,d Paolo Lacovig,e Luca Bignardi,e Hendrik Bluhm,f Jan Knudsen,d, g Bjørk Hammer,a, c Liv Hornekaera, c a

iNANO, University of Aarhus, DK-8000 Aarhus C, Denmark b

c

Aix-Marseille Univ, CNRS, PIIM, Marseille, France

Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, Denmark d

e

The MAX IV Laboratory, Lund University, Sweden

Elettra - Sincrotrone Trieste S.C.p.A., S.S. 14 km 163.5, 34012 Trieste, Italy

f

Chemical Sciences Division and Advanced Light Source, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 g

Division of Synchrotron Radiation Research, Lund University, Sweden

Corresponding author: *Email: [email protected] (L.H) KEYWORDS: graphene, vibrational excitation, nanostructured functionalization, band gap engineering, molecular hydrogen, catalysis

1 ACS Paragon Plus Environment

Page 3 of 25 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 Nano

Abstract Hydrogen functionalization of graphene by exposure to vibrationally excited H2 molecules is investigated by combined scanning tunneling microscopy, high resolution electron energy loss spectroscopy, x-ray photoelectron spectroscopy measurements and density functional theory calculations. The measurements reveal that vibrationally excited H2 molecules dissociatively adsorb on graphene on Ir(111) resulting in nano-patterned hydrogen functionalization structures. Calculations demonstrate that the presence of the Ir surface below the graphene lowers the H2 dissociative adsorption barrier and allows for the adsorption reaction at energies well below the dissociation threshold of the H-H bond. The first reacting H2 molecule must contain considerable vibrational energy to overcome the dissociative adsorption barrier. However, this initial adsorption further activates the surface resulting in reduced barriers for dissociative adsorption of subsequent H2 molecules. This enables functionalization by H2 molecules with lower vibrational energy, yielding an avalanche effect for the hydrogenation reaction. These results provide an example of a catalytically active graphene coated surface and additionally set the stage for a reinterpretation of previous experimental work involving elevated H2 background gas pressures in the presence of hot filaments.

2 ACS Paragon Plus Environment

ACS Nano 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 4 of 25

The controlled chemical functionalization of graphene is of broad interest for modifying its chemical, optical, mechanical and electronic properties. The adsorption of a single H atom onto the graphene lattice induces magnetic moments in the sheet1 enabling the development of spintronic devices. At a higher hydrogen coverage, manipulation of the global electronic structure is possible and in some systems a tunable band gap can be obtained, which is of interest for the use of graphene in electronic applications.2-3 A complete hydrogenation of graphene by H atoms yields the two-dimensional electronically insulating sp3 structure called graphane,4 which has been proposed as a hydrogen storage material.5-6 The interaction between graphene and its substrate is a crucial factor for the reactivity of graphene and thus for its chemical functionalization.7 Several experimental and theoretical studies have shown that hot atomic hydrogen can chemisorb onto the basal plane of graphite, free-standing graphene and graphene on metal and semiconductor substrates.2-4, 8-16 For graphene on Ir(111) (gr/Ir(111)) the reactivity varies over the surface, since the graphene-Ir interaction is modulated with a ~25 Å periodicity, reflecting the moiré pattern that originates from the ~10% lattice mismatch between graphene and Ir(111).17-18 In so-called HCP and FCC areas of the moiré structure the position of every second carbon atom coincides with the position of an Ir atom below.2-3 In these HCP and FCC areas, the graphene lattice can distort upon hydrogen functionalization and form a configuration where every second C atom binds to the underlying Ir, while neighboring C atoms bind to H atoms on top.3, 10 The configuration becomes thermodynamically stable at a certain H cluster size and adsorption at the HCP site is found to be slightly more favorable than the FCC site. 2 A selective functionalization of HCP areas by hot H atoms, has been achieved by elevating the sample temperature during the atomic exposure. Such selective functionalization leads to highly

ACS Paragon Plus Environment

3

Page 5 of 25 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 Nano

ordered hydrogenation structures, including the opening of a gap in the electronic band structure.2

While functionalization of graphene with atomic H is well established, functionalization using H2 molecules is less explored19-21 and is in general energetically unfavorable due to the high stability of the H2 molecule. For free-standing graphene the barrier for dissociative adsorption is estimated to 3.1-3.4 eV.8, 22-25 It has been predicted, however, that the barrier can be lowered in the presence of a perpendicular electric field26-27 or at defects and edges.24, 28-29 For non-free-standing graphene, the barrier can be substantially lower as shown in this report. A few experimental studies report partial hydrogenation of graphene at high pressure H2 (PH2 > 1 bar).19-21 Additionally, it has been proposed that nanoparticles can dissociate H2 and transfer the resulting hydrogen atoms to the graphene through spillover effects.30-32 The validity of this spillover phenomenon was however questioned in a recent theoretical study.33 The H2 adsorption energy and activation barrier is consistently reported to depend on the graphene curvature.9, 34-35

Here we present combined scanning tunneling microscopy (STM), high resolution electron energy loss spectroscopy (HREELS) and x-ray photoemission spectroscopy (XPS) data, that demonstrate that highly vibrationally excited H2 molecules can functionalize gr/Ir(111), yielding highly ordered nano-patterned graphene, via dissociative adsorption reactions. Ambient pressure XPS (APXPS) excludes the possibility of graphene hydrogenation through direct dissociation of H2 molecules on the Ir(111) surface. The experimental observations are supported by density functional (DFT) calculations showing that the dissociation barrier is lowered substantially by the graphene-mediated catalytic effect of the Ir substrate. Moreover, the barrier to adsorption

ACS Paragon Plus Environment

4

ACS Nano 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 25

continues to decrease for further addition of H2 at sites adjacent to adsorbed molecules due to local distortion effects. Consequently, an avalanche of H2 dissociation and H-addition reactions becomes possible, providing an efficient route towards hydrogenation.

Results and discussion Highly vibrationally excited H2 is produced via Eley-Rideal abstraction reactions between gas phase H atoms and H atoms adsorbed on metal surfaces. Such reactions have been shown to yield vibrationally excited H2 molecules exhibiting super-thermal vibrational energy level populations up to at least ν=9.36-37 The exact population distribution depends on the metal surface on which the recombination occurs.36 While gaseous H atoms are expected to either adsorb or undergo abstraction reactions upon collision with metal surfaces, the vibrationally excited H2 molecules are found to survive hundreds of collisions before thermalizing.36 A high yield of highly vibrationally excited molecules was obtained by letting H2 or D2 gas pass through a ~2300 K ThO2Ir filament (method i) or a hot W capillary with temperature T W (method ii). The filament and W capillary are both sources of hot atomic H. Thus, to avoid the influence of these hot H atoms in our experiments, the sample was always placed out of line of sight and far from any atomic H source. This yields a complicated trajectory for hot H atoms between the source and the sample surface. In their trajectory, the gas phase H atoms will inevitably collide with metallic elements. During such collisions, they have a high probability to recombine with adsorbed H atoms from the metallic surfaces to yield gas phase vibrationally excited H2.36 To exclude the possibility that any hot H atoms reach the sample surface in our experiments, a highly oriented pyrolytic graphite (HOPG) sample was exposed to vibrationally excited H 2 as described above. Graphite becomes hydrogen functionalized when exposed to hot H atoms.8 In

ACS Paragon Plus Environment

5

Page 7 of 25 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 Nano

our experiments with vibrationally excited H2 the graphite was, however, completely inert, as shown by HREELS (see Supporting Information) and STM measurements, even for exposure times up to ~12 hours, and the influence of hot H atoms can thus be excluded.

STM results Fig. 1a-d display STM images of gr/Ir(111) following exposure to vibrationally excited H 2 or D2 molecules. Fig. 1a depicts gr/Ir(111) after a ~32 min. exposure (method i) at a background pressure of 2×10-5 mbar H2. Three types of bright features related to the H2 exposure are identified and indicated by arrows: Disc shaped (type I), donut shaped (type II) and merged (type III). Type I and II features are outlined in Fig. 1b and a type III feature is outlined in Fig. 1c. Type III features are represented by the merging of neighboring type I and/or II features. It is observed (not shown here) that during a continuous exposure to vibrationally excited H2, type I features evolve into type II, which may merge at higher exposure into type III. Type II and III features sometimes appear instantaneously, which we ascribe to the limited time resolution of the STM imaging. The image in Fig. 1a is the last in a series of STM images (see Supporting Information, S1) obtained during a continuous exposure to vibrationally excited H2. The accumulated number of nucleation sites was obtained from these images at different exposure times and plotted in Fig. 1e (red dots). The plot shows a linear trend. The fractional apparent functionalized area vs. time is plotted in Fig. 1e (blue triangles) displaying a square dependency on the fluence. These findings indicate, that the dissociative adsorption of excited H2 is more likely to lead to an expansion of existing clusters rather than the initiation of new ones. The absolute values in this plot are not a direct measure of the hydrogen coverage since the appearance of features in STM varies with the bias voltage applied. From the STM data in Fig. 1

ACS Paragon Plus Environment

6

ACS Nano 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 8 of 25

it is evident that features of type I and II appear almost exclusively on one area of the moiré structure. In Fig. 1b this area is seen to be the brightest part of the moiré unit cell, which is identified as the HCP site.2, 18, 38 When other moiré areas are occupied, it is observed to almost exclusively occur in a type III structure where type II features have grown together to occupy neighboring moiré areas. Inspection of Fig. 1a reveals a strong tendency for H-functionalized areas to cluster in neighboring moiré unit cells. In the inset in Fig. 1f, a Fourier transform of a larger scale STM image (see Supporting Information, S2), obtained on the same area as the image presented in Fig. 1d, is shown. The hexagonal pattern observed in the Fourier transform underlines, that hydrogen functionalization structures appear preferentially on one site on the moiré unit cell conferring a high degree of ordering to the functionalized surface. A line profile along the line indicated on the Fourier transform is presented. The separation between the spots in the Fourier transform, indicated by arrows on the line profile, corresponds to a real space distance of ~21.5 Å. This equals ~25 Å × cos(30°) and thus demonstrates the preferred occupancy of one site in the moiré super lattice, identified in the STM images as the HCP site. The exposure to vibrationally excited H2, thus leads to a highly ordered global patterning of the graphene, with moiré periodicity. No obvious isotopic dependence was observed in the STM experiments. An entirely clean gr/Ir(111) surface can be regained by annealing the sample above ~720 K.

ACS Paragon Plus Environment

7

Page 9 of 25 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 Nano

Figure 1: a)-d) STM images of gr/Ir(111) exposed to vibrationally excited H2 or D2: a) H2 with method i) at P=2×10-5 mbar for 32 min. b) D2 with method ii) at P=5×10-7 mbar for 20 min, TW > 2000 K. Features of types I and II are indicated with circles. The rhombus outlines a moiré unit cell, the “*” and the “+” denote the FCC and HCP areas, respectively and the corners of the rhombus mark the atop regions. c) D2 with method ii) at P=5×10-7 mbar for 60 min, TW > 2000 K. A type III feature is outlined. d) D2 with method ii) at P=5×10-7 mbar for 60 min, TW > 2000 K. e) Time evolution of hydrogenation structures during exposure to vibrationally excited H 2, using method i), at P=1×10-5 mbar. Red dots: total number of nucleation sites, normalized to the imaged area, plotted against time. Error bars represent the intrinsic error assuming a Poisson distribution. Blue triangles: Apparent hydrogenated area, in % of total image area, plotted against time. f) A line profile through the Fourier transform shown in the inset, along the line indicated.

ACS Paragon Plus Environment

8

ACS Nano 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 25

The Fourier transform is performed on a larger scale STM image (see Supporting Information, S2) in the same area as shown in the STM image in d). The Fourier transform illustrates a high degree of ordering on the surface after the exposure to vibrationally excited H2. The separation of the peaks indicated by arrows in the line profile, corresponds to a real space separation of ~21.5 Å which equals ~25 Å × cos(30°), confirming a global hydrogen induced patterning with moiré superlattice periodicity. Imaging parameters for a-d): a) Vt=478.2 mV, It=0.790 nA. b) Vt=67.1 mV, It=1.090 nA. c) Vt= -351.9 mV, It= -0.310 nA. d) Vt= -351.9 mV, It= -0.320 nA.

HREELS results Fig. 2 depicts HREEL spectra of functionalized graphene after exposure to D atoms (red crosses) and after exposure to vibrationally excited D2 molecules (black circles), both at saturation coverage. The HREEL spectrum obtained after exposure to highly excited D2 molecules depicts all characteristic features found for gr/Ir(111) functionalized with atomic D.10 The two most intense vibrational features in the spectra are positioned at 108 meV and 245 meV corresponding to the C-D bending (δC-D) and stretching (νC-D) modes, respectively.10 The presence of these features in both spectra, prove that regardless of deuterium precursor, atomic D or vibrationally excited D2, a formation of C-D covalent bonds takes place yielding chemically functionalized gr/Ir(111). The two spectra differ only in intensity (reflecting different saturation coverages) and by the peak at 195 meV, corresponding to the clean graphene TO/LO phonon, which is stronger in the spectrum for the D2 exposed sample. The persistence of this phonon after vibrationally excited D2 saturation exposure indicates that larger graphene areas are unfunctionalized, consistent with the lower intensity of the C-D vibrational features.10

ACS Paragon Plus Environment

9

Page 11 of 25 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 Nano

Figure 2: HREEL spectra of gr/Ir(111) exposed to D atoms (flux 3.5×1013 D/cm2s, fluence 1.0×1016 D/cm2) (red) and to vibrationally excited D2 using method i) at 2×10-4 mbar for 45 min. (black), both up to saturation. Incident energy E0 = 5 eV and momentum transfer of 0.07-0.1 Å-1 over the presented energy range. XPS results XPS measurements were performed on gr/Ir(111) as a function of the exposure time to vibrationally excited H2 produced by method ii) at a background pressure of 4.5×10-8 mbar. A spectrum at saturation coverage is shown in Fig. 3a. The spectrum was fit following the procedure used to analyse gr/Ir(111) exposed to hot H atoms,11 and is composed of an sp2 carbon component (Cc) and three C-H related components at positions +0.81 eV (Ca), +0.44 eV (Cb) and -0.27 eV (Cd) relative to the position of Cc. The C-H related components have been suggested11 to originate from H clusters in FCC and/or HCP regions (Cb), vacancies in the clusters as well as C atoms neighboring C-H bonds (Cd) and finally C-H in atop regions (Ca). Fig. 3b shows the evolution of the relative intensities of the four components during the continuous exposure of gr/Ir(111) to vibrationally excited H2. The spectra shift to higher binding energy as the degree of functionalization increases, which we attribute to doping, and only the

ACS Paragon Plus Environment

10

ACS Nano 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 25

relative positions of the components were fixed during the fitting procedure (see Methods). At saturation coverage, the loss of sp2 carbon (Cc component) was ~62%, which is lower than what has been reported for the functionalization using atomic H, where a ~70% loss in the Cc component was observed.11 This is consistent with HREELS in Fig. 2. Besides the lower saturation coverage, the main difference for vibrationally excited H2 compared to H atom functionalization, is a relatively lower contribution from the Ca component in the case of the former. To verify the described mechanism for production of vibrationally excited H2, the temperature of the W capillary (TW) was varied, keeping a constant pressure of 5×10-7 mbar and exposure time of 20 min. The sample was flashed to T = 900 K between different exposures. The relative component intensities plotted against TW are shown in Fig. 3c. It can be seen that functionalization is initiated at an onset of TW > 1800 K, which corresponds to the onset of atomic H production,39 in agreement with the proposed mechanism for the generation of highly vibrationally excited H2 described above.

Figure 3: a) XPS (hν = 400 eV) C1s core level spectrum of gr/Ir(111) after ~1.5 h exposure to vibrationally excited H2 at P=4.5×10−8 mbar, produced by method ii), TW = ~2018 K (saturation coverage). The spectrum is the last spectrum from the uptake series shown in b). b) Uptake of vibrationally excited H2 shown as the relative intensities of components Ca, Cb, Cc and Cd over

ACS Paragon Plus Environment

11

Page 13 of 25 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 Nano

time during exposure at same parameters as in a). c) The relative intensities of components Ca, Cb, Cc and Cd in the XPS (hν = 390 eV) C1s core level spectrum for a series of experiments where gr/Ir(111) was exposed for 20 min P=5×10−

7

mbar of H2, method ii), with a variable

temperature TW. DFT results Fig. 4a displays calculated reaction pathways for the dissociative adsorption of H2 on the HCP areas of the moiré unit cell for dissociation and adsorption into an ortho-(blue), para-(red) and meta-(green) configurations. Full lines represent calculations where the graphene and the Ir (except the bottom Ir layer) are allowed to fully relax. The dashed lines represent the same calculations for free-standing graphene with no underlying Ir substrate. All energies are referenced to the relaxed graphene and gas phase H2 molecules. Consistent with previous reports, the H2 adsorption configuration with the lowest barrier, on free-standing graphene, is found to be the para-configuration (3.22 eV, dashed red line). This barrier, as well as the transition and final state, is unchanged in the presence of an underlying Ir surface (3.19 eV, full red line). However, for the adsorption into an ortho- or meta-configuration the barrier when Ir is present is lowered by 0.40 eV to 3.86 eV and by 1.14 eV to 2.85 eV respectively (compare dashed to full lines). The presence of the Ir substrate thus changes the most likely adsorption configuration (lowest barrier) from the para- to the meta-position. Additionally, meta-adsorption gives the most thermodynamically stable product due to the possibility for the center C atom to bind down to the Ir surface, leading to a local sp2 to sp3 rehybridization of the graphene lattice. The internal energy needed for a gas-phase H2 molecule to overcome the 2.85 eV barrier for adsorption into the meta-configuration requires vibrational excitation to the ν≥7 level (~3.0 eV).36

ACS Paragon Plus Environment

12

ACS Nano 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 14 of 25

Figure 4: a) Reaction pathways for one H2 molecule to adsorb onto ortho-, meta-, and para-sites of an HCP region of graphene over an Ir(111) surface. Solid lines refer to when the graphene sheet and Ir (except for the bottom layer) are allowed to fully relax. Dashed lines reflect the calculations for free-standing graphene. Calculated transition states are plotted (side-view) as well as the final adsorption into the ortho- para- and meta-configuration (top-view). The pathway for the ortho-site is for a semi-constrained reaction. b) The reaction pathway for a second H2 molecule adsorbing in the vicinity of an adsorbed meta-dimer. c) The reaction pathway for a third H2 molecule adsorbing in the vicinity of the first two meta-adsorbed H2 molecules with a top-view of the final configuration. d) The potential energy surface around the transition state into the meta-configuration expanded in the H2 height, Z, above the graphene and the internal H-H bond length, b. The table gives the calculated barrier heights for the three dimer configurations for: free-standing graphene (no Ir); in the case where the Ir atoms are not allowed to move (Ir con); and for the fully relaxed system (Ir) displayed in the figure.

ACS Paragon Plus Environment

13

Page 15 of 25 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 Nano

The barrier for adsorption into the para-configuration would require the population of the ν≥8 level (~3.3 eV).36 If an H2 molecule adsorbs into the para-configuration on gr/Ir(111) (red full line) it is thermodynamically favorable for it to diffuse (~0.8 eV diffusion barrier) to the more stable meta-configuration. No H atom diffusion is expected to occur, once dissociation into the stable meta-configuration is established. Adsorption into the ortho-configuration is unlikely since an excitation to the ν≥11 level is required to overcome the ~4 eV barrier.36 While exhibiting different adsorption barriers, the reaction coordinates for adsorption into the para- or meta-configurations also undergo very different transition states, see insets in Fig. 4. For the adsorption in para-configuration, two C atoms in the graphene lattice move slightly upwards to accommodate H2 adsorption in a para-dimer configuration. In the case of metaadsorption, an Ir atom moves slightly out of the Ir(111) plane towards the graphene lattice and several C atoms move toward the Ir surface simultaneously as the H2 approaches. After binding to a C atom, the Ir atom relaxes back into the Ir(111) plane, while the adjacent C atoms bind to the H atoms from the dissociated H2. Thus, a stabilized sp3 hybridized state is produced locally in the graphene sheet. The rearrangement of atoms in the described transition state for meta-site adsorption, especially the movement of the Ir atom in and out of the plane, happens on a timescale much slower than the vibrational timescale of the excited H2 molecule. Hence it is unclear whether the incoming excited H2 molecule will experience the fully relaxed potential energy surface. The potential energy barriers for hydrogen adsorption were explored, by keeping the Ir lattice constrained in its undisturbed state (“Ir con.” in the table in Fig. 4), since the movement of only the graphene sheet is more likely to match the relevant timescales. The barrier is unchanged for ortho- and para-adsorption, as expected, and increases only by 0.09 eV for meta-adsorption. With such a small increase it seems that the movement of the Ir atom is not the

ACS Paragon Plus Environment

14

ACS Nano 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 16 of 25

determining factor for the reaction mechanism, while the mere presence of the Ir surface has a huge impact on the barriers, through stabilization of the final state. The DFT calculations, see Fig. 4b-c, also reveal that the adsorption of a subsequent impinging H2 molecule, into the vicinity of an H2 that has already dissociated and adsorbed into a metaconfiguration on gr/Ir(111), has a barrier that is 0.76 eV lower and the reaction is exothermic. Furthermore, the barrier is lowered by an additional 0.19 eV for the third impinging H2 molecule, which then only needs vibrational excitation to the ν≥4 level to overcome a 1.90 eV adsorption barrier. This is consistent with our STM observations in Fig. 1e, which show a faster increase of the total area of adsorbed hydrogen (triangles in Fig. 1e) compared to the increase in number of nucleation sites (circles in Fig. 1e). The barriers are expected to further decrease for subsequent H2 adsorption until the H-C-Ir cluster becomes thermodynamically stable at a critical size.2 Thus, the adsorption of the first H2 molecule triggers an avalanche reaction for the formation of clusters mainly in HCP sites as observed by STM. These clusters may extend to include adsorption into FCC sites via the perimeter of the HCP-site cluster as the distortion of the graphene makes more sites energetically available for dissociation of the highly excited molecules. This is consistent with the STM observations for higher coverages (Fig. 1c-f). Fig. 4d shows an elbow plot of the potential energy surface (PES) in the vicinity of the transition state for adsorbing H2 into the meta-site obtained by moving the H atoms only. The transition state is located in the curved region of the PES. Such energy landscapes have been reported to lead to strong vibrational coupling.40 The adsorption scheme for hydrogen at the atop areas of the gr/Ir(111) moiré unit cell should be similar to free-standing graphene and graphite. The dissociation of H2 on graphite predominantly occurs in a para-adsorption configuration.8 The related adsorption barrier is

ACS Paragon Plus Environment

15

Page 17 of 25 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 Nano

expected to be similar to that for free-standing graphene and is significantly higher than the barriers found for dissociative adsorption of hydrogen on gr/Ir(111). Hence, only hot H atoms can adsorb onto graphite as well as on the atop regions of gr/Ir(111), which is consistent with HREELS (Fig. 2) and XPS measurements (Fig. 3), showing that saturation coverage for exposure to vibrationally excited H2 is lower than for exposure to atomic H. However, even though the barrier for para-site adsorption is higher, vibrationally excited H2 in ν=8 or ν=9, should have enough energy to overcome it. This contradicts our observations, indicating that the energy loss to the surface, needed to stabilize the reaction products in the dissociative adsorption process, is less efficient for free-standing graphene. To corroborate the presented interpretation of our results the viability of an alternative scenario involving H2 dissociation on bare patches of Ir and subsequent reaction with the overlying graphene in a spillover-like scheme30-32 was investigated by APXPS. The APXPS measurements (see Supporting Information, S4) show that non-vibrationally excited H2 molecules, even at highpressure exposures above 0.01 mbar, do not lead to chemisorption of H onto gr/Ir(111). The experiments yield identical negative results for the monolayer as well as sub-monolayer (Ir exposed) gr/Ir(111) samples.

Conclusions To conclude, we have demonstrated a route to graphene functionalization using vibrationally excited H2 molecules. The functionalization occurs in a highly ordered manner and is therefore a viable tool for band gap engineering of graphene. The functionalization route exploits the finding that graphene can mediate the catalytic activity of an underlying metal substrate resulting in reduced dissociative adsorption barriers for molecular hydrogen on graphene. Furthermore, the

ACS Paragon Plus Environment

16

ACS Nano 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 25

described mechanism exhibits an avalanche effect where the first dissociative adsorption event leads to reduced barriers for subsequent dissociative adsorption. Not only is this study an important proof of concept, but the essential insight provided by the mechanism, where nucleation is the most energetically expensive step for the graphene hydrogenation, provides an approach towards the detailed control of graphene and its properties, through chemical functionalization. The presented results further demonstrate that at high chamber pressures of H2 an atomic source, like an ion gauge filament, may produce vibrationally excited H2 molecules, with potential influence on the outcome of the experiment. In the absence of a source of atomic hydrogen, vibrationally excited H2 molecules may be produced via VUV or UV excitation processes.41

Methods

Experimental details

A high-quality gr/Ir(111) sample was prepared by chemical vapor deposition procedures and/or temperature programmed growth.17 STM was performed using home-built Aarhus-type STMs42 at Aarhus University, Aarhus, Denmark. HREELS was performed at Aix-Marseille Université, Marseille, France in a setup described elsewhere.10 APXPS was performed at beamline 11.0.2 at the Advanced Light Source, Berkeley, California and at SPECIES43-45 at the MAX II storage ring of MaxIV (formerly Max lab), Lund, Sweden, time-resolved XPS was performed at the SuperESCA beamline at Elettra, Trieste, Italy and conventional high-resolution XPS was performed at beamline I311 at the MAX II storage ring of MaxIV (formerly Max lab), Lund, Sweden. All experiments were performed under ultrahigh vacuum (UHV) conditions, unless stated otherwise.

ACS Paragon Plus Environment

17

Page 19 of 25 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 Nano

Hydrogen background pressures used in the experiments varied from 1×10-6 - 2×10-4 mbar (method i) and 4.5×10-8 - 7×10-7 mbar (method ii). The mean free path of H2 molecules at 300 K at these pressures corresponds to 0.5 m - 2400 m, exceeding the geometry of the chamber setups.

Computational methods

Density functional theory (DFT) was used to calculate adsorption potential energies and reaction barriers for H2 on graphene supported by Ir(111). The real-space, grid-based, projector augmented wave method GPAW46-47 was employed and supported by the Atomic Simulation Environment (ASE).48 In order to include the dispersion related effects which govern the interaction between graphene and the Ir surface, the vdW enabled functional optB88-vdW was used.49 The parameters used for the presented calculations, which reproduce experimental observations2, are as follows: an 8 × 8 graphene lattice over a 7 × 7 Ir(111) surface simulated the surface for the adsorption of H2. The Ir slab was composed of three layers where the bottom layer was held fixed. The dimensions of the cell parallel to the surface are defined by the optB88vdW calculated lattice constant for graphene (agraphene,optB88−vdW = 2.465). This functional uses the vdW-DF framework proposed by Dion et. al.50 with a reparameterized version of B88 exchange51 fitted to the S22 dataset.52 It expands the Ir lattice by 1.6% relative to the optB88-vdW calculated value of 3.921 Å. At least 7.5 Å separates any atom from the top or bottom of the cell. A grid spacing of 0.1779 Å defines the length of each grid point. Finally, a 2 × 2 k-point mesh was sampled and all forces were relaxed until they were below 0.02 eV/Å. Reaction pathways are determined using the AutoNEB code found in ASE with the CI-EB method.53 As it is not possible to make the direct hydrogenation of the ortho-configuration we have used a constraint CI-EB relaxation for completeness - to ensure that a direct barrier was found. A free CI-EB relaxation forces the atoms first into the meta-configuration before moving into the orthoconfiguration. In the constrained case the transition state has been allowed to only move in the plane

ACS Paragon Plus Environment

18

ACS Nano 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 20 of 25

crossing the right of the two attacked atoms in the insert in Fig. 4 with a normal vector pointing in the direction of the left attacked atom.

Data Analysis

STM images were lightly filtered from noise and flattened to compensate for piezo creep. HREELS spectra were normalized and calibrated to the intensities and positions respectively of their elastic peak. No further data handling was performed for the presented spectra. XPS data (obtained from SuperESCA, Elettra, Trieste (Fig. 3a,b), from I311 Maxlab, Lund (Fig. 3c, S3) and beamline 11.0.2 at the Advanced Light Source, Berkeley, California (Fig. S4) was fitted using Doniach-Sunjic profiles convoluted with a Gaussian profile (DS°G). Clean gr/Ir were fit using a Lorentzian full width at half maximum (LFWHM) of 0.135 eV and a Gaussian full width at half maximum (GFWHM) of 0.177 eV for data obtained at SuperESCA, and LFWHM = 0.158 eV, GFWHM = 0.129 eV for data obtained at I311. The GFWHM was allowed to broaden slightly during the vibrationally excited H2 exposure, yielding a maximum of GFWHM = 0.279 eV for the last spectrum in the uptake series presented in Fig. 3a and 3b. For the T W = 2018 K exposure, presented in Fig. S3, the Cc component has GFWHM = 0.217 eV. The C1s core levels change in binding energy upon formation of CH bonds due to doping effects.11 The component related to clean sp2 graphene (Cc) was therefore allowed to shift its binding energy position, when fitting spectra that represent vibrationally excited H2 exposed samples, and the positions of the C-H related components (Ca, Cb and Cd) were fixed relative to the position of Cc by +0.81, +0.44 eV and -0.27 eV respectively. This procedure is equivalent to the fitting scheme for the data presented by Balog et al.11 for gr/Ir(111) exposed to hot H atoms. Peak assignment was also done according to this previous work. APXPS data (obtained at beamline 11.0.2 at the Advanced Light Source, Berkeley) shown in Fig. S4 (left inset) were fit with a LFWHM = 0.135 and GFWHM =

ACS Paragon Plus Environment

19

Page 21 of 25 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 Nano

0.35 eV. The binding energy scale for all XPS data was calibrated to the Fermi edge. Linear backgrounds were used for the fits.

ACKNOWLEDGMENT We acknowledge financial support from The Danish Council for Independent Research (grant no. 0602-02566B and grant no. 0602-02265B), Innovation Fund Denmark (NIAGRA), The European Research Council (CoG GRANN), Villum Centre of Excellence for Dirac Materials. H.B. acknowledges support by the Director, Office of Science, Office of Basic Energy Sciences, and by the Division of Chemical Sciences, Geosciences and Biosciences of the U.S. Department of Energy at LBNL under Contract No. DE-AC02-05CH11231. The Advanced Light Source is supported by the Director, Office of Science, Office of Basic Energy Sciences of the U.S. Department of Energy at LBNL under Contract No. DE-AC02-05CH11231. ASSOCIATED CONTENT The authors declare no competing financial interest. Supporting Information Available:

See Online Supporting Information for supplemental

STM, HREELS and XPS data. AUTHOR INFORMATION Corresponding Author *Email: [email protected] (L.H) Author Contributions

ACS Paragon Plus Environment

20

ACS Nano 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 22 of 25

L.K., R. Bisson, R. Balog, and L.H. designed the experiments. L.K., S.H., and J.M., performed STM experiments. L.K., R. Bisson, R. Balog and T.A. performed HREELS experiments. L.K., A.C, J.J., R. Balog, A. Č., M.A., J.K., P.L. and L.B. performed XPS measurements. L.K., H.B., S.U., and J.K. performed APXPS measurements. M.G., E.K. and B.H. performed DFT calculations. A. G. Č. helped prepare samples. L.K., R. Bisson, R. Balog and L.H. analyzed the data and wrote the manuscript. All authors contributed to discussions.

1. González-Herrero, H.; Gómez-Rodríguez, J. M.; Mallet, P.; Moaied, M.; Palacios, J. J.; Salgado, C.; Ugeda, M. M.; Veuillen, J.-Y.; Yndurain, F.; Brihuega, I., Atomic-Scale Control of Graphene Magnetism by Using Hydrogen Atoms. Science 2016, 352, 437-441. 2. Jørgensen, J. H.; Čabo, A. G. i.; Balog, R.; Kyhl, L.; Groves, M. N.; Cassidy, A. M.; Bruix, A.; Bianchi, M.; Dendzik, M.; Arman, M. A., Symmetry-Driven Band Gap Engineering in Hydrogen Functionalized Graphene. ACS Nano 2016, 10, 10798-10807. 3. Balog, R.; Jorgensen, B.; Nilsson, L.; Andersen, M.; Rienks, E.; Bianchi, M.; Fanetti, M.; Laegsgaard, E.; Baraldi, A.; Lizzit, S.; Sljivancanin, Z.; Besenbacher, F.; Hammer, B.; Pedersen, T. G.; Hofmann, P.; Hornekaer, L., Bandgap Opening in Graphene Induced by Patterned Hydrogen Adsorption. Nat. Mater. 2010, 9, 315-319. 4. Elias, D.; Nair, R.; Mohiuddin, T.; Morozov, S.; Blake, P.; Halsall, M.; Ferrari, A.; Boukhvalov, D.; Katsnelson, M.; Geim, A., Control of Graphene's Properties by Reversible Hydrogenation: Evidence for Graphane. Science 2009, 323, 610-613. Tozzini, V.; Pellegrini, V., Prospects for Hydrogen Storage in Graphene. Phys. Chem. 5. Chem. Phys. 2013, 15, 80-89. 6. Camiola, V. D.; Farchioni, R.; Cavallucci, T.; Rossi, A.; Pellegrini, V.; Tozzini, V., Hydrogen Storage in Rippled Graphene: Perspectives from Multi-Scale Simulations. Front. Mater. 2015, 2, 3. Trotochaud, L.; Head, A. R.; Karslıoğlu, O.; Kyhl, L.; Bluhm, H., Ambient Pressure 7. Photoelectron Spectroscopy: Practical Considerations and Experimental Frontiers. J. Phys.: Condens. Matter 2016, 29, 053002. 8. Hornekaer, L.; Sljivancanin, Z.; Xu, W.; Otero, R.; Rauls, E.; Stensgaard, I.; Laegsgaard, E.; Hammer, B.; Besenbacher, F., Metastable Structures and Recombination Pathways for Atomic Hydrogen on the Graphite (0001) Surface. Phys. Rev. Lett. 2006, 96, 156104. 9. Balog, R.; Jørgensen, B.; Wells, J.; Lægsgaard, E.; Hofmann, P.; Besenbacher, F.; Hornekær, L., Atomic Hydrogen Adsorbate Structures on Graphene. J. Am. Chem. Soc. 2009, 131, 8744-8745. 10. Kyhl, L.; Balog, R.; Angot, T.; Hornekaer, L.; Bisson, R., Hydrogenated Graphene on Ir (111): A High-Resolution Electron Energy Loss Spectroscopy Study of the Vibrational Spectrum. Phys. Rev. B 2016, 93, 115403.

ACS Paragon Plus Environment

21

Page 23 of 25 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 Nano

11. Balog, R.; Andersen, M.; Jorgensen, B.; Sljivancanin, Z.; Hammer, B.; Baraldi, A.; Larciprete, R.; Hofmann, P.; Hornekaer, L.; Lizzit, S., Controlling Hydrogenation of Graphene on Ir(111). ACS Nano 2013, 7, 3823-3832. 12. Bocquet, F. C.; Bisson, R.; Themlin, J. M.; Layet, J. M.; Angot, T., Deuterium Adsorption on (and Desorption from) Sic(0001)-(3 X 3), (Root 3 X Root 3) R30 Degrees, (6 Root 3 X 6 Root 3) R30 Degrees and Quasi-Free-Standing Graphene Obtained by Hydrogen Intercalation. J. Phys. D: Appl. Phys. 2014, 47, 094014. 13. Hyunil, K.; Balgar, T.; Hasselbrink, E., Is There Sp 3-Bound H on Epitaxial Graphene? Evidence for Adsorption on Both Sides of the Sheet. Chem. Phys. Lett. 2012, 546, 12-17. 14. Hyunil, K.; Balgar, T.; Hasselbrink, E., The Stretching Vibration of Hydrogen Adsorbed on Epitaxial Graphene Studied by Sum-Frequency Generation Spectroscopy. Chem. Phys. Lett. 2011, 508, 1-5. 15. Luo, Z.; Yu, T.; Kim, K.-j.; Ni, Z.; You, Y.; Lim, S.; Shen, Z.; Wang, S.; Lin, J., Thickness-Dependent Reversible Hydrogenation of Graphene Layers. ACS Nano 2009, 3, 17811788. Zhao, W.; Gebhardt, J.; Späth, F.; Gotterbarm, K.; Gleichweit, C.; Steinrück, H. P.; 16. Görling, A.; Papp, C., Reversible Hydrogenation of Graphene on Ni (111)—Synthesis of “Graphone”. Chem. – Eur. J. 2015, 21, 3347-3358. 17. Coraux, J.; N'Diaye, A. T.; Engler, M.; Busse, C.; Wall, D.; Buckanie, N.; Heringdorf, F.; van Gastel, R.; Poelsema, B.; Michely, T., Growth of Graphene on Ir(111). New J. Phys. 2009, 11, 023006. 18. N'Diaye, A. T.; Coraux, J.; Plasa, T. N.; Busse, C.; Michely, T., Structure of Epitaxial Graphene on Ir(111). New J. Phys. 2008, 10, 043033. 19. Kim, B. H.; Hong, S. J.; Baek, S. J.; Jeong, H. Y.; Park, N.; Lee, M.; Lee, S. W.; Park, M.; Chu, S. W.; Shin, H. S.; Lim, J.; Lee, J. C.; Jun, Y.; Park, Y. W., N-Type Graphene Induced by Dissociative H2 Adsorption at Room Temperature. Sci. Rep. 2012, 2, 690. 20. Smith, D.; Howie, R. T.; Crowe, I. F.; Simionescu, C. L.; Muryn, C.; Vishnyakov, V.; Novoselov, K. S.; Kim, Y.-J.; Halsall, M. P.; Gregoryanz, E., Hydrogenation of Graphene by Reaction at High Pressure and High Temperature. ACS Nano 2015, 9, 8279-8283. 21. Shin, D. S.; Kim, Y. B.; Kim, D. Y.; Choi, T. H.; Kim, B. H., Direct Measurement of the Amount of Dissociated Hydrogen Atoms Attached on Graphene. Synth. Met. 2015, 200, 80-84. 22. Miura, Y.; Kasai, H.; Dino, W.; Nakanishi, H.; Sugimoto, T., First Principles Studies for the Dissociative Adsorption of H2 on Graphene. J. Appl. Phys. 2003, 93, 3395-3400. 23. McKay, H.; Wales, D. J.; Jenkins, S. J.; Verges, J. A.; deAndres, P. L., Hydrogen on Graphene under Stress: Molecular Dissociation and Gap Opening. Phys. Rev. B 2010, 81, 075425. 24. Wang, N.; Wang, L.; Tan, Q.; Pan, Y.-X., Effects of Hydroxyl Group on H 2 Dissociation on Graphene: A Density Functional Theory Study. J. Energy Chem. 2013, 22, 493497. 25. Doi, K.; Onishi, I.; Kawano, S., Dissociative Adsorption of H2 Molecules on Steric Graphene Surface: Ab Initio MD Study Based on DFT. Comput. Theor. Chem. 2012, 994, 54-64. 26. Ao, Z.; Peeters, F., Electric Field Activated Hydrogen Dissociative Adsorption to Nitrogen-Doped Graphene. J. Phys. Chem. C 2010, 114, 14503-14509. 27. Ao, Z.; Peeters, F., Electric Field: A Catalyst for Hydrogenation of Graphene. Appl. Phys. Lett. 2010, 96, 253106.

ACS Paragon Plus Environment

22

ACS Nano 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 24 of 25

28. Allouche, A.; Ferro, Y., Dissociative Adsorption of Small Molecules at Vacancies on the Graphite (0001) Surface. Carbon 2006, 44, 3320-3327. 29. Diño, W. A.; Nakanishi, H.; Kasai, H.; Sugimoto, T.; Kondo, T., H2 Dissociative Adsorption at the Zigzag Edges of Graphite. e-J. Surf. Sci. Nanotechnol. 2004, 2, 77-80. 30. Divya, P.; Ramaprabhu, S., Hydrogen Storage in Platinum Decorated Hydrogen Exfoliated Graphene Sheets by Spillover Mechanism. Phys. Chem. Chem. Phys. 2014, 16, 26725-26729. 31. Parambhath, V. B.; Nagar, R.; Sethupathi, K.; Ramaprabhu, S., Investigation of Spillover Mechanism in Palladium Decorated Hydrogen Exfoliated Functionalized Graphene. J. Phys. Chem. C 2011, 115, 15679-15685. 32. Gao, Y.; Zhao, N.; Li, J.; Liu, E.; He, C.; Shi, C., Hydrogen Spillover Storage on CaDecorated Graphene. Int. J. Hydrogen Energy 2012, 37, 11835-11841. 33. Blanco-Rey, M.; Juaristi, J. I. a.; Alducin, M.; Lopez, M. J.; Alonso, J. A., Is Spillover Relevant for Hydrogen Adsorption and Storage in Porous Carbons Doped with Palladium Nanoparticles? J. Phys. Chem. C 2016, 120, 17357-17364. Tozzini, V.; Pellegrini, V., Reversible Hydrogen Storage by Controlled Buckling of 34. Graphene Layers. J. Phys. Chem. C 2011, 115, 25523-25528. 35. Goler, S.; Coletti, C.; Tozzini, V.; Piazza, V.; Mashoff, T.; Beltram, F.; Pellegrini, V.; Heun, S., Influence of Graphene Curvature on Hydrogen Adsorption: Toward Hydrogen Storage Devices. J. Phys. Chem. C 2013, 117, 11506-11513. Hall, R. I.; Cadez, I.; Landau, M.; Pichou, F.; Schermann, C., Vibrational Excitation of 36. Hydrogen via Recombinative Desorption of Atomic Hydrogen Gas on a Metal Surface. Phys. Rev. Lett. 1988, 60, 337. 37. Eenshuistra, P. J.; Bonnie, J. H. M.; Los, J.; Hopman, H. J., Observation of Exceptionally High Vibrational Excitation of Hydrogen Molecules Formed by Wall Recombination. Phys. Rev. Lett. 1988, 60, 341. 38. Gerber, T.; Busse, C.; Mysliveček, J.; Coraux, J.; Michely, T., A Versatile Fabrication Method for Cluster Superlattices. New J. Phys. 2009, 11, 103045. 39. Tschersich, K., Intensity of a Source of Atomic Hydrogen Based on a Hot Capillary. J. Appl. Phys. 2000, 87, 2565-2573. Gross, A., Theoretical Surface Science. Springer: 2014; Vol. 1. 40. 41. Chelkowski, S.; Bandrauk, A., Raman Chirped Adiabatic Passage: A New Method for Selective Excitation of High Vibrational States. J. Raman spectrosc. 1997, 28, 459-466. 42. Laegsgaard, E.; Besenbacher, F.; Mortensen, K.; Stensgaard, I., A Fully Automated, Thimble-Size Scanning Tunnelling Microscope. J. Microsc. (Oxford, U.K.) 1988, 152, 663-669. 43. Knudsen, J.; Andersen, J. N.; Schnadt, J., A Versatile Instrument for Ambient Pressure X-Ray Photoelectron Spectroscopy: The Lund Cell Approach. Surf. Sci. 2016, 646, 160-169. 44. Schnadt, J.; Knudsen, J.; Andersen, J. N.; Siegbahn, H.; Pietzsch, A.; Hennies, F.; Johansson, N.; Mårtensson, N.; Öhrwall, G.; Bahr, S., The New Ambient-Pressure X-Ray Photoelectron Spectroscopy Instrument at Max-Lab. J. Synchrotron Radiat. 2012, 19, 701-704. 45. Urpelainen, S.; Såthe, C.; Grizolli, W.; Agåker, M.; Head, A. R.; Andersson, M.; Huang, S.-W.; Jensen, B. N.; Wallén, E.; Tarawneh, H., The Species Beamline at the Max Iv Laboratory: A Facility for Soft X-Ray Rixs and Apxps. J. Synchrotron Radiat. 2017, 24, 344-353. 46. Mortensen, J. J.; Hansen, L. B.; Jacobsen, K. W., Real-Space Grid Implementation of the Projector Augmented Wave Method. Phys. Rev. B 2005, 71, 035109.

ACS Paragon Plus Environment

23

Page 25 of 25 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 Nano

47. Enkovaara, J.; Rostgaard, C.; Mortensen, J. J.; Chen, J.; Dułak, M.; Ferrighi, L.; Gavnholt, J.; Glinsvad, C.; Haikola, V.; Hansen, H., Electronic Structure Calculations with Gpaw: A Real-Space Implementation of the Projector Augmented-Wave Method. J. Phys.: Condens. Matter 2010, 22, 253202. 48. Larsen, A.; Mortensen, J.; Blomqvist, J.; Castelli, I.; Christensen, R.; Dulak, M.; Friis, J.; Groves, M.; Hammer, B.; Hargus, C., The Atomic Simulation Environment—a Python Library for Working with Atoms. J. Phys.: Condens. Matter 2017, 273002. 49. Klimeš, J.; Bowler, D. R.; Michaelides, A., A Critical Assessment of Theoretical Methods for Finding Reaction Pathways and Transition States of Surface Processes. J. Phys.: Condens. Matter 2010, 22, 074203. 50. Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I., Van Der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. 51. Becke, A. D., Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098. 52. Jurečka, P.; Šponer, J.; Černý, J.; Hobza, P., Benchmark Database of Accurate (Mp2 and Ccsd (T) Complete Basis Set Limit) Interaction Energies of Small Model Complexes, DNA Base Pairs, and Amino Acid Pairs. Phys. Chem. Chem. Phys. 2006, 8, 1985-1993. 53. Kolsbjerg, E. L.; Groves, M. N.; Hammer, B., An Automated Nudged Elastic Band Method. J. Chem. Phys. 2016, 145, 094107.

ACS Paragon Plus Environment

24