Adsorption of Dichlorobenzene on Au and Pt Stepped Surfaces Using

Sep 4, 2012 - The adsorption of dichlorobenzene on flat (111) and stepped (332) Au and Pt surfaces was studied using density functional theory with bo...
0 downloads 0 Views 1MB Size
Article pubs.acs.org/JPCC

Adsorption of Dichlorobenzene on Au and Pt Stepped Surfaces Using van der Waals Density Functional Theory Rengin Peköz,* Karen Johnston, and Davide Donadio Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany ABSTRACT: The adsorption of dichlorobenzene on flat (111) and stepped (332) Au and Pt surfaces was studied using density functional theory with both a conventional generalized gradient approximation (GGA) and a fully nonlocal van der Waals density functional (vdWDF). The equilibrium geometries and adsorption energies were computed for several different adsorption configurations. The two functionals yielded qualitatively different results, with the GGA functional predicting only weak binding compared to vdW-DF, demonstrating the importance of including nonlocal dispersion. By analyzing the electronic density and projected density of states, it was found that the interaction of dichlorobenzene with the two surfaces caused a charge redistribution, especially for the stepped surfaces. Moreover, adsorption on the step edge on Au(332) was dominated by nonlocal dispersion, whereas adsorption on the Pt(332) step was dominated by chemical bonding.

1. INTRODUCTION Graphene nanoribbons (GNRs) have a nonzero band gap, and the band structure can be tuned by changing the width and shape of the edges, which makes them very appealing for electronic applications. To tailor defect-free, atomically precise GNRs with desired properties, a well-controlled growth mechanism is required. Among the different synthesis methods recently applied, a bottom-up approach using surface-supported synthesis is a powerful technique because of its ease in controlling the quality of the nanoribbons.1 Surface-supported GNR fabrication converts halogenated aromatic precursor molecules into linear polyphenylenes, which are then transformed into extended fully aromatic structures by cyclodehydrogenation. Many experimental and theoretical studies have focused on understanding the details of these chemical processes. Theoretical studies, in particular, have investigated the cyclodehydrogenation of individual molecules and the selfassembly of molecules on surfaces as well as the effect of the substrate on such processes.2−7 Much less is known about the adsorption, diffusion, and self-assembly of molecular precursors into polymer chains. Improved control of the growth setup and of the final products stems from a fundamental understanding of the initial stages of GNR fabrication and, in particular, on the interaction between molecular precursors and metal surfaces. Corrugated surfaces (such as high-index vicinal surfaces,8−10 or superstructures formed due to lattice mismatch11) have been proposed to be more suitable than flat surfaces for building ordered 2D molecular structures12−15 since, especially at low coverage, step edges are preferential adsorption sites.16−18 Using low-temperature scanning tunneling microscopy (STM), Pascual and co-workers17 found that benzene molecules preferentially adsorb on step edges, while terraces remain free of molecules for low coverage. Increasing the coverage resulted in single rows of molecules adsorbed at the step edges. Very recently, Zhang et al.18 investigated the adsorption properties of iron-phthalocyanine on the 22 × √3 reconstructed Au(111) © 2012 American Chemical Society

surface with monatomic steps using a combination of lowtemperature STM and density functional theory (DFT) calculations, and they also reported adsorption at the step edges. A study by Treier et al.12 of the large organic molecule hexa-peri-hexabenzocoronene on the Au(11 12 12) surface used a combination of STM, synchrotron-based X-ray photoelectron diffraction (XPD), and DFT calculations and found that the stable adsorption structure is tilted and situated on kink sites. Although there are many studies on the adsorption behavior of aromatic molecules, such as benzene, pyridine, pyrazine, thiophene, cyclopentandienyl, and n-butane,19−22 on different metal surfaces, only few are concerned with halogenated molecules. Despite the fact that it is not possible to grow GNRs from dichlorobenzene (DCB), DCB is a simple prototype for larger and more complex GNR precursors. For example, Milosevic et al.23 investigated experimentally the constitutional meta−para isomerization of DCB adsorbed on Ag(111) and Cu(111) surfaces, but gave no details about the adsorption behavior. Fagan and co-workers24,25 explored the interaction of 1,2-DCB with semiconducting and metallic single-walled carbon nanotubes (SWCNTs) based on both DFT calculations and experiments. However, none of these studies revealed any information about the adsorption properties of DCB molecules on metal surfaces. In this work, we perform first-principles calculations to investigate the role of van der Waals (vdW) interactions on the adsorption of the 1,4-DCB molecule on flat (111) and stepped Au and Pt surfaces. Even though high-index metallic surfaces exhibit complex reconstructions,26 we here consider (332) ascleaved metal surfaces as model step-terrace systems, with a relatively small size and a well-defined geometry that permits comparisons between different metals at an affordable Received: June 20, 2012 Revised: September 3, 2012 Published: September 4, 2012 20409

dx.doi.org/10.1021/jp306062s | J. Phys. Chem. C 2012, 116, 20409−20416

The Journal of Physical Chemistry C

Article

3. RESULTS AND DISCUSSION 3.1. Adsorption Structure and Energy. The adsorption of DCB on Au(111) and Pt(111) surfaces is investigated for flat-lying orientations, in which the center of the molecule is positioned on the top (T), bridge (B), or fcc hollow (H) sites. For each site, two orientations are considered, with an angle of θ = 0° or 30° between the Cl−Cl vector and the (100) axis, as shown in Figure 1.

computational cost. We employ the vdW density functional (vdW-DF) for general geometries proposed by Dion et al.,27 using the implementation described by Gulans et al.28 Adsorption sites, geometries, energies, and electronic densities are explored using vdW-DF, and the results are compared to those obtained using a conventional GGA functional. Furthermore, the effect of the substrate is also explored by comparing the adsorption properties of DCB on Au and Pt surfaces, which have different electronic structures, that is, fully occupied d band (Au) versus partially occupied d band (Pt).

2. THEORY AND METHODS All calculations were performed within the framework of DFT as implemented in the VASP 4.6 package.29,30 The electronic wave functions were expanded in a plane-wave basis set, and frozen-core all-electron projector augmented wave (PAW)31 potentials were used to describe the core electrons. The VASP code was adapted to include vdW interactions via a selfconsistent implementation of the vdW-DF functional.27,28 Both the PBE version32,33 of the generalized gradient approximation (GGA) and the vdW-DF functional27,28 with PBE exchange were used. This fully nonlocal vdW functional has been successfully applied to a broad spectrum of systems, including layered structures, such as graphite34 and boron nitride,35 dimers of atoms and small molecules,28 polycyclic aromatic hydrocarbons (PAHs),36 silicon37,38 and metals,21,39−41 polymers,42 and carbon nanotubes (CNTs).43 Four different surfaces were used as substrates: the flat Au(111) and Pt(111) surfaces and the stepped Au(332) and Pt(332) surfaces. The Au(Pt) flat surface has a hexagonal unit cell with dimensions a = 2.95(2.81) Å for PBE and a = 2.99(2.85) Å for vdW-DF. Both the flat and the stepped surfaces have four atomic layers with the two bottom layers fixed during the simulations. In the stepped surfaces, the terraces have a (111) orientation and the steps are spaced just over five atomic rows apart. A vacuum region of 30 Å separates adjacent slabs, and the cutoff energy for the plane waves was 400 eV. The Brillouin zone integration was performed using the Monkhorst−Pack scheme with 4 × 4 × 1 and 2 × 3 × 1 mesh for the flat and stepped surfaces, respectively. All structures were optimized with a convergence criterion of 10 meV/Å for the forces and 0.1 meV for the energy. Several test calculations for the convergence of results with the number of fixed bottom layers and vacuum space were also carried out. The binding energy for DCB on Pt decreases by 0.02 eV from fixing the bottom two layers to fixing only the bottom layer. Increasing the vacuum region from 20 to 40 Å increased the binding energy by 0.001 eV. For the adsorption calculations on (111) surfaces, a 4 × 4 unit cell was constructed, corresponding to a coverage of 0.0625 monolayers (ML), where one ML is defined here as one molecule adsorbed on one surface unit cell. The adsorption energy, Eads, at a given vertical distance z of the molecular plane from the surface, is defined by Eads = Emolecule + Esurface − Etotal

Figure 1. Au(111) surface showing a dichlorobenzene molecule on the fcc hollow site, H, oriented with angle θ = 30°. The T, B, and H labels show where the center of the molecule would be on top, bridge, and fcc hollow sites, respectively. Atoms are colored to distinguish between the layers. Top, second, and third layers are colored from light to dark yellow.

Figure 2 shows the dependence of the binding energy of DCB on the top and hollow sites of Au(111) and Pt(111) as a

(1)

where Etotal, Emolecule, and Esurface are the total energies of the adsorbed system, the isolated DCB molecule, and the clean metal surface, respectively. A positive value of Eads means that the adsorbed system is energetically favorable relative to the isolated state. The adsorption energy and distance from the surface at equilibrium is denoted by Eeq ads and zeq, respectively.

Figure 2. Calculated adsorption energy of a dicholorobenzene molecule on (a) Au(111) and (b) Pt(111) surfaces with PBE and vdW-DF functionals as a function of the distance between the surface and the center of the molecule. 20410

dx.doi.org/10.1021/jp306062s | J. Phys. Chem. C 2012, 116, 20409−20416

The Journal of Physical Chemistry C

Article

the three high symmetry sites (hollow, top, and bridge). Therefore, an upper limit for the diffusion barrier can be estimated as the difference between the adsorption energies of the highest- and lowest-energy configurations, which is ∼0.05 eV for Au and ∼0.13 eV for Pt. This low diffusion barrier is consistent with previous results regarding benzene diffusion on metal surfaces. STM images showed that, even at 4 K with low coverage, benzene molecules were mobile over the Au(111) surfaces with diffusion barriers ∼ 0.15 eV.46 While benzene molecules on Pt(111) were frozen at 4 K, they were able to diffuse at room temperature,47 and in the 280−520 K temperature range, the diffusion barrier was measured to be ∼0.43 eV.48,49 Recent vdW-DF calculations reported by Berland et al.50 slightly underestimated the diffusion barrier for benzene on the Cu(111) surface. We now consider the adsorption behavior of DCB on the (332) stepped surfaces. The lower symmetry of the (332) surfaces yields several inequivalent adsorption sites along the step edge, several of which are shown in Figure 3. Note that, for

function of the perpendicular distance (z) between the metal surface and the center of the molecule. For both Au and Pt, the PBE functional predicts a negligible equilibrium binding energy of about 0.1 eV, whereas vdW-DF yields equilibrium binding energies of 1.0−1.1 eV. This shows that, in both systems, dispersion forces are the dominant interaction, as previously demonstrated for benzene on gold.22 Furthermore, inclusion of the vdW interaction brings the molecule closer to the surface and decreases zeq by 0.1−0.2 Å. The equilibrium distances, zeq, and adsorption energies, Eeq ads, of DCB on the high symmetry sites of Au(111) and Pt(111) using vdW-DF are displayed in Table 1. We found that the Table 1. Adsorption Sites (See Figure 1), Angles (θ), Energies (Eeq ads), Equilibrium Distances (zeq), and Charge Displacement (ΔQ) for Dichlorobenzene on Au(111) and Pt(111) Surfaces from vdW-DF Calculations Au(111)

Pt(111)

site

θ [°]

Eeq ads [eV]

zeq [Å]

ΔQ [e]

Eeq ads [eV]

zeq [Å]

ΔQ [e]

T T H H B B

0 30 0 30 0 30

1.00 0.97 1.00 1.02 1.01 1.01

3.51 3.51 3.50 3.46 3.49 3.48

0.26 0.22 0.25 0.29 0.25 0.26

1.12 1.06 1.14 1.19 1.17 1.19

3.41 3.48 3.34 3.29 3.31 3.23

0.46 0.30 0.44 0.60 0.54 0.62

most favorable adsorption site for Au(111) is H-30 with Eeq ads = 1.02 eV and zeq = 3.46 Å. The difference in adsorption energy between this site and the least-favorable one (T-30) is only 0.05 eV. The energy differences between the various configurations of DCB on Au(111) are at the limit of the accuracy of DFTbased approaches, and thus, our calculations cannot be conclusive about the optimal geometry but, instead, indicate that the potential energy surface of DCB on Au(111) displays several nearly equivalent minima. For Pt(111), the energy differences between different adsorption sites are also small (∼0.1 eV), but the adsorption of DCB on Pt(111) is slightly stronger than on Au. The adsorption energy of the optimal configurations (H-30 and B-30) is 1.19 eV, and the equilibrium distance is shorter by 0.2−0.3 Å than on the Au(111) surface. For the B-30 configuration, both Cl atoms lie on low-symmetry surface sites, and the molecule is found to lie flat on both Au and Pt surfaces. In contrast, for the H-30 configuration on Au(Pt), one Cl atom lies on a top site and the other on a hollow hcp site, with a coordination of one and three metal atoms, respectively. This leads to a tilted adsorption structure, where the Cl atom on the top site is closer to the surface by 0.08(0.33) Å than the Cl atom on the hcp site. The corresponding tilt angle of the H-30 configuration on Au and Pt is 1° and 3.5°, respectively. The adsorption differences between Au and Pt stems from the stronger coupling between the p bands of the DCB to the d bands of Pt, as will be discussed in depth in the following sections. While it is expected that DCB interaction with Au is almost exclusively dispersive, it is surprising that the different position and electronic occupation of the d band of Pt, with respect to Au, produces only minor effects in the adsorption of DCB on the (111) surfaces. Diffusion barriers for the (111) surfaces were calculated using the nudged elastic band (NEB) method44,45 implemented in the VASP code. Only configurations with θ = 30° were considered, and it was found that there are no barriers between

Figure 3. Schematic diagram of the Au(332) surface showing DCB with various initial adsorption configurations. Upper and lower panels show the top and side views, respectively, and the black lines indicate the unit cell of Au(332). The step-edge Au atoms are shown with white spheres. The details of the configurations are given in Table 2. The optimized (c) and (e) configurations have a tilted DCB structure similar to, but less tilted than, configuration (g).

the (332) surfaces, we redefine θ to be the angle between the Cl−Cl vector and the step edge. Eleven different configurations of DCB on Au(332) and four on Pt(332) were calculated with the vdW-DF functional, and the results are reported in Table 2. Since the structure of the (332) terrace is the same as that of the (111) surface, we expect similar adsorption properties at 20411

dx.doi.org/10.1021/jp306062s | J. Phys. Chem. C 2012, 116, 20409−20416

The Journal of Physical Chemistry C

Article

Table 2. Adsorption Configurations (See Figure 3), Angles (θ), Energies (Eeq ads), Equilibrium Distances (zeq), and Charge Displacement (ΔQ) for Dichlorobenzene on Au(332) and Pt(332) Surfaces Using the vdW-DF Functionala configuration (a) (b) (c) (d) (e) (f) (g) (h)

S S S Tr S S S Tr Tr Tr Tr

Tshift T B H H Ttilt Ttilt T T T H

Au(332)

Pt(332)

θ [°]

Eeq ads [eV]

zeq [Å]

ΔQ [e]

0 0 0 0 0 90 0 90 30 60 30

0.92 0.91 0.89 0.98 0.90 0.99 1.10 0.96 0.96 0.98 0.99

3.24 3.54 3.25 3.41 3.43 3.42 2.95 3.40 3.48 3.53 3.40

0.38 0.48 0.39 0.34 0.33 0.61 0.60 0.31 0.25 0.30 0.36

Eeq ads [eV]

zeq [Å]

ΔQ [e]

1.84

2.27

2.72

1.76 1.60 1.06

2.29 2.40 3.49

2.20 1.81 0.40

a Abbreviations: S, step; Tr, terrace; T, top; B, bridge; H, hollow. θ is the angle that the Cl−Cl vector makes with the step edge. zeq is the perpendicular distance between the geometric center of the molecule and the average height of the terrace atoms, except for configuration (g), which is tilted on the step, where it is the perpendicular distance between the plane of the geometric center of dichlorobenzene and the step edge atoms.

Figure 4. PDOS of DCB molecule adsorbed on (a) Au(111) with configuration H-30, (b) Pt(111) with configuration H-30, (c) Au(332) with configuration (f), and (d) Pt(332) with configuration (f). The green and black lines in the upper panels indicate the surface atom PDOS when DCB is adsorbed and for the clean surface, respectively. The Fermi level is set to 0 eV. Electronic states are broadened by 0.2 eV.

energy of Eeq ads = 1.84 eV compared to 1.19 eV (B-30 site) on the flat surface. Differences between the binding energies on flat surfaces and at the step edge are very large also at the PBE level. The PBE energy of the S-Ttilt-90 configuration is Eeq ads = 1.14 eV, which is ∼1 eV more attractive than the configurations on Pt(111). By comparing vdW-DF versus PBE calculations, flat versus stepped, and Au versus Pt surfaces, it is clear that the electronic properties of the substrate have an important effect on the binding properties; that is, Eeq ads of DCB on Au is smaller than that of Pt for both functionals and surface types. The higher adsorption energies of the configurations on the step edges are due to the different atomic environment, such as the coordination of the step atoms, which leads to different electronic properties and a higher reactivity of the stepped surface (Smoluchowski effect).51 Moreover, for the Au stepped surfaces, the dominant interaction still arises from vdW forces,

terrace sites (denoted as Tr in Table 2). This is indeed the case for both Au and Pt: Eads and zeq are nearly the same on the (111) surface and (332) terrace. On Au(332), only configuration (g), in which DCB lies on the step with the Cl atoms aligned with the step edge, has a larger equilibrium adsorption energy than the terrace sites, with Eeq ads = 1.10 eV. This adsorption energy is only 0.08 eV larger than that on the flat (111) surface, indicating that there is no substantial difference in the chemical interactions between the two situations and that, for both flat and stepped Au surfaces, the nature of the interaction is predominantly dispersive. This interpretation is supported by the calculation of the S-Ttilt-0 configuration at the PBE level, which yields Eeq ads = 1.12 eV, which is only ∼0.05 eV larger than that of the configurations on the flat surface. On Pt(332), the chemical behavior is rather different, and the binding energy on the step is much larger. For example, the S-B-0, shown in Figure 3c, has an adsorption 20412

dx.doi.org/10.1021/jp306062s | J. Phys. Chem. C 2012, 116, 20409−20416

The Journal of Physical Chemistry C

Article

Figure 5. 1D planar-averaged total electron density, ρtotal(z), and electronic density difference, Δρav(z), for DCB and benzene (BZ) in the H-30 configuration adsorbed on (a) Au(111) and (d) Pt(111). Δρav(z) is multiplied by 400. 3D isosurface of the electronic density differences with an isovalue of ±0.003 e/Å3 of DCB (b, e) and benzene (c, f) on Au(111) and Pt(111) surfaces are shown. Blue and white regions in 3D plots represent charge accumulation and depletion, respectively.

molecule, leading to a stronger metal/adsorbate interaction.52 However, the effect on the d band of Au is smaller than that on Pt. The difference between the adsorption of DCB on flat and stepped surfaces is further analyzed by examining the electronic charge density. The charge density difference or deformation due to the adsorption of the molecule on the surface is calculated as

whereas in the case of the Pt step, the main contribution to binding energy is due to chemical bonding. 3.2. Electronic Structure. To further interpret the chemistry of the adsorption processes, we analyzed charge density and partial density of states (PDOS). For both Au and Pt, we consider the H-30 structure on the (111) surfaces and configuration (f) on the (332) surfaces. In Figure 4a,c, we show the PDOS for DCB on Au(111) and Au(332), respectively. The top panel in each case shows the Au 5d states, and the bottom panel shows the 2s(3s) and 2p(3p) states of the C(Cl) atoms of DCB. When DCB adsorbs on Au(111) with the H-30 configuration, neither the PDOS of the gold surface nor the molecule changes significantly, except for small differences in the C 2p and Cl 3p valence states just below the Fermi level, as seen in Figure 4a. The adsorption of DCB on Au(332) has a more pronounced effect on the molecular states, which shift and decrease in intensity, indicating a mixing between the electronic states of DCB and the d band of Au, as seen in Figure 4c. The higher adsorption energies of DCB on Pt surfaces imply that there should be a stronger hybridization between the molecular and the surface states compared to the Au case. The PDOSs of DCB on Pt(111) and Pt(332) are shown in Figure 4b,d, respectively. The change in the molecular PDOS of DCB on Pt is clear even for the (111) surface, where the three highest occupied molecular states and the two lowest unoccupied molecular states are slightly shifted and their intensity is reduced by the interaction with the surface. For DCB on Pt(332), the two highest occupied states and the unoccupied C 2p states are almost bleached, as seen in Figure 4d. Furthermore, while there is no significant effect of DCB adsorption on the d states of the flat surfaces, the d states of both Au(332) and Pt(332) are coupled to the p states of the

Δρ = ρtotal − ρsurface − ρmolecule

(2)

where ρtotal is the electronic charge density of DCB adsorbed on the surface, ρsurface is the charge density of the Au or Pt slab in the adsorbed configuration without the presence of DCB, and ρmolecule is the charge density of the molecule in the adsorbed configuration in vacuum. For a more quantitative analysis, the average charge density difference profile, Δρav(z), is calculated by averaging over the plane perpendicular to the (111) surface for both the flat and the stepped surfaces Δρav (z) =

∫ Δρ dx dy

(3)

The total charge densities and electronic density differences of DCB with configuration H-30 (and benzene, which will be discussed below) on Au(111) and Pt(111) are shown in Figure 5. When DCB is adsorbed on the Au or Pt surface, it has an overall positive charge and there are regions of electronic density accumulation (Δρav > 0) and depletion (Δρav < 0) in the metal. For DCB on Au(111), the charge accumulation (blue regions in Figure 5b) is localized above the top layer of atoms (Figure 5a,b), whereas for Pt(111), the charge displacement is larger and penetrates into the subsurface layer (Figure 5d,e). In addition to the charge accumulation at the surface, there is also a small area of charge depletion between 20413

dx.doi.org/10.1021/jp306062s | J. Phys. Chem. C 2012, 116, 20409−20416

The Journal of Physical Chemistry C

Article

Figure 6. Average electronic density difference, Δρav(z), of DCB on (a) Au(111) with the H-30 configuration, (b) Pt(111) with the H-30 configuration, (c) Au(332) with configuration (f), and (d) Pt(332) with configuration (f). Δρav(z) is multiplied by 200. The positions of the surface and subsurface atoms for Au and Pt and for DCB on the (111) surface are indicated by vertical lines. The shaded gray area represents the tilted DCB on the (332) surface. 3D electronic density differences of DCB on Au(332) and Pt(332) with an isovalue of ±0.003 e/Å3 are shown in (e) and (f) with side and top views of the unit cell. The blue and white regions represent charge accumulation and depletion, respectively.

Pt(111) with the H-30 configuration using PBE was found to be 0.40 eV, which is only 0.1 eV larger than the result of Morin et al.58 The PBE calculations together with experiments59 show that benzene is chemisorbed on Pt(111). This is qualitatively different from the case of DCB on Pt(111), which has a PBE adsorption energy of around 0.1 eV, showing that DCB is not chemically bonded to the surface. The different binding properties of benzene and DCB on Pt(111) arise because the large Cl atoms keep the aromatic ring structure farther from the surface, thus preventing chemical bonding. The electronic density differences of benzene on Au(111) and Pt(111) are shown in Figure 5. For both surfaces, the shape of Δρav(z) for benzene and DCB is very similar, although, for benzene, the amplitude of the charge fluctuations is larger, especially on the Pt surface. This is due to a larger overlap of the electronic states of benzene and the d band of Pt. The 3D electronic density differences, Δρav, of benzene on Au and Pt, shown in Figure 5c,f, respectively, provide more information about the electronic properties of these structures. For benzene in the H-30 configuration, all the hydrogen atoms are on lowsymmetry sites, which preserves the symmetry of the system, resulting in a completely flat molecule lying on the surface with a symmetric charge distribution (see Figure 5c,f), in contrast to the asymmetric and tilted DCB molecule with the same configuration (see Figure 5b,e). The charge reorganization due to the adsorption of the molecule on the surface can be analyzed quantitatively by integrating the charge density difference profile. Following ref 60, we define the charge displacement as

the surface and the Cl atom on the top site. As mentioned previously, this Cl atom is closer to the surface than the other Cl atom by 0.33 Å, which corresponds to the observed tilt of 3.5° between the DCB plane and the Pt surface. This alternating charge density pattern between the surface and the Cl atom induces a dipole. The averaged electronic density difference profile, Δρav(z), reveals that there is almost no dipole formation for the lower layers of the Au(111) surface, whereas for Pt(111), the charge difference oscillates down to the lowest layers. For the lowest two layers, the oscillation is very small and decreases toward the bottom surface layer, which can also be seen in the 3D electronic density difference (Δρ; see Figure 5b,e). To elucidate the specific role of the chlorine atoms, we compare the electronic structure of DCB on metal surfaces with benzene on metal surfaces. The latter system has been widely investigated both experimentally and theoretically. Theoretical calculations,53 which used the same vdW-DF and exchange functionals as the present study, predicted an adsorption energy of ∼0.8 eV for benzene on Au(111), in reasonable agreement with experiment54 and in good agreement with other calculations.22,55 The experimental adsorption energy of benzene on Pt(111) was measured in the range of 1.84−2.25 eV56 for low coverage, and 1.21−1.34 eV49,57 for higher coverages. DFT calculations using PBE, which does not include dispersion, predicted that, for benzene on Pt(111), the adsorption energy is in the range of 0.30−0.90 eV depending on the adsorption site,58 which is consistent with our PBE calculations. For example, the binding energy of benzene/ 20414

dx.doi.org/10.1021/jp306062s | J. Phys. Chem. C 2012, 116, 20409−20416

The Journal of Physical Chemistry C ΔQ =

∫ |Δρ(r)| d3r

Article

calculations with a conventional GGA-PBE and a vdW-DF functional. We found that adsorption on both Au(111) and Pt(111) surfaces is governed by dispersion interactions. The slightly stronger binding energy of DCB on Pt(111) compared to Au(111) stems from the larger interaction between the d orbitals of Pt and the highest occupied molecular orbitals of DCB. For the stepped surfaces, the step edges were found to be preferential adsorption sites; however, large differences are observed between Pt and Au. DCB on the Au(332) surface binds only slightly more strongly than on the flat surface with a difference in adsorption energy of around 0.1 eV, and the interaction remains mostly dispersive. In contrast, on the Pt(332) surface, DCB is chemisorbed at the edge of the steps with an adsorption energy 0.4−0.6 eV larger than that on the terrace or flat surface. In view of exploiting stepped metallic surfaces as supports to grow GNR, we can conclude that, for Au, the very small differences in adsorption energies between steps and terraces may be insufficient to induce a significant alignment of graphene precursors along the edge of the steps at room temperature. The diffusion barriers of DCB on (111) surfaces were found to be very low, for both Au and Pt surfaces. Therefore, on the stepped surface, molecules will diffuse freely along the terraces, but in the case of Pt, they would get trapped at step edges. This investigation of adsorption and diffusion of small precursor molecules should provide valuable information to improve the synthesis of GNRs. However, further work is needed to elucidate the kinetics of molecular diffusion on the surfaces on larger scales. Future work will focus on the coupling of precursor molecules and multiscale modeling of the diffusion and dynamics of molecules on gold surfaces.

(4)

The values of ΔQ are reported in Tables 1 and 2. ΔQ correlates well with binding energies and is inversely proportional to adsorption distances. This general trend is not completely fulfilled for the case of dispersive adsorption of DCB on Au; however, adsorption on step sites leads to larger ΔQ than on flat surfaces. A stronger relation between binding energy and ΔQ is found for DCB on Pt, as different charge displacements reflect different types of bonding, from physisorption to chemisorption. For DCB on Au(111), ΔQ is about 0.29 e, and for benzene on Au(111), ΔQ = 0.49 e. On Pt(111), ΔQ is 0.62 e for DCB and 1.18 e for benzene. Consistent with the 3D plots, it is clear that benzene has a larger charge displacement on both surfaces than DCB and that the Pt surface has a large charge displacement than Au. The work function of a surface is defined as Φmetal = Evac − EF, where Evac and EF are the electrostatic potential in the center of the vacuum region and Fermi energy of the systems, respectively. The deformation dipole due to the charge rearrangement at the interface mentioned above is thus partly responsible for the decrease of the work function upon adsorption. The calculated ΦAu(111) and ΦPt(111) are 5.47 and 5.98 eV, which are in good agreement with the experimental results of 5.3161,62 and 5.70 eV,62 respectively. The adsorption of DCB(benzene) on flat Au and Pt surfaces decreases the work function by 0.15(0.23) and 0.25(0.41) eV, respectively. This is rather small in the case of DCB and indicates that there are no significant dipoles formed between DCB and the surfaces. For the stepped surfaces, we explore the electronic density differences taking configuration (f) as an example. The electronic density difference profiles, Δρav(z), are rather different from those of the flat surfaces, as shown in Figure 6. For DCB on Au, the region between the surface and the subsurface layers has a small accumulation of charge (see Figure 6c). For DCB on Pt, the fluctuations are much larger and the region between the surface and the subsurface layers has both accumulation and depletion regions (see Figure 6d). A small charge depletion can be seen in the region between the surface and the lower part of the DCB molecule. Within the range of the DCB molecule, there is a large accumulation layer, followed by a charge depletion region, the magnitude of which is larger than the accumulation part. Furthermore, for the stepped surfaces, there is a weak polarization of the electronic states, which penetrates only into the Au subsurface layer but continues down to the bottom layer for Pt. From the 3D representation of electronic density difference, shown in Figure 6e,f, it is observed that the dipole is formed along the direction perpendicular to the DCB molecule through the step edge. For configuration (f) on Au(332), shown in Figure 6c, the charge displacement is equal to 0.60 e. Configuration (f) on Pt(332) has ΔQ = 2.72 e. This demonstrates the large difference in the type of bonding of DCB to Au and Pt stepped surfaces. The charge rearrangement between the aromatic regions of DCB adsorbed on both Au and Pt steps is more pronounced than on flat surfaces. This is supported by both the 3D electronic density difference plots and the total integrated charge displacement data.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We thank Björn Baumeier for a critical reading of the manuscript, and Xinliang Feng and Alexandre Tkatchenko for valuable discussions. We acknowledge the RZG of the Max Planck Society for computational resources. R.P. and D.D. acknowledge funding from the MPRG program of the Max Planck Society. Funding for K.J. was provided by the German Research Foundation Priority Program DFG SPP 1369.



REFERENCES

(1) Cai, J.; Ruffieux, P.; Jaafar, R.; Bieri, M.; Braun, T.; Blankenburg, S.; Muoth, M.; Seitsonen, A. P.; Saleh, M.; Feng, X.; et al. Nature 2010, 466, 470−473. (2) Bieri, M.; Nguyen, M.-T.; Gröning, O.; Cai, J.; Treier, M.; AïtMansour, K.; Ruffieux, P.; Pignedoli, C. A.; Passerone, D.; Kastler, M.; et al. J. Am. Chem. Soc. 2010, 132, 16669−16676. (3) Blankenburg, S.; Rauls, E.; Schmidt, W. G. J. Phys. Chem. Lett. 2010, 1, 3266−3270. (4) Nguyen, M. T.; Pignedoli, C. A.; Passerone, D. Phys. Chem. Chem. Phys. 2011, 13, 154−160. (5) Björk, J.; Stafström, S.; Hanke, F. J. Am. Chem. Soc. 2011, 133, 14884−14887. (6) Treier, M.; Pignedoli, C. A.; Laino, T.; Rieger, R.; Müllen, K.; Passerone, D.; Fasel, R. Nat. Chem. 2011, 3, 61−67.

4. CONCLUSIONS We investigated the adsorption properties of dichlorobenzene on flat (111) and stepped (332) Au and Pt surfaces using DFT 20415

dx.doi.org/10.1021/jp306062s | J. Phys. Chem. C 2012, 116, 20409−20416

The Journal of Physical Chemistry C

Article

(7) Peköz, R.; Feng, X.; Donadio, D. J. Phys.: Condens. Matter 2012, 24, 104023−104031. (8) Rousset, S.; Repain, V.; Baudot, G.; Garreau, Y.; Lecoeur, J. J. Phys.: Condens. Matter 2003, 15, S3363−S3392. (9) Mugarza, A.; Ortega, J. E. J. Phys.: Condens. Matter 2003, 15, S3281−S3310. (10) Kröger, J.; Jensen, H.; Néel, N.; Berndt, R. Surf. Sci. 2007, 601, 4180−4184. (11) Corso, M.; Auwärter, W.; Muntwiler, M.; Tamai, A.; Greber, T.; Osterwalder, J. Science 2004, 303, 217−220. (12) Treier, M.; Ruffieux, P.; Schillinger, R.; Greber, T.; Müllen, K.; Fasel, R. Surf. Sci. 2008, 602, L84−L88. (13) Canas-Ventura, M. E.; Xiao, W.; Wasserfallen, D.; Müllen, K.; Brune, H.; Barth, J. V.; Fasel, R. Angew. Chem., Int. Ed. 2007, 46, 1814−1818. (14) Vladimirova, M.; Stengel, M.; Vita, A. D.; Baldereschi, A.; Bö hringer, M.; Morgenstern, K.; Berndt, R.; Schneider, W.-D. Europhys. Lett. 2001, 56, 254−260. (15) Jiang, Y. H.; Xiao, W. D.; Liu, L. W.; Zhang, L. Z.; Lian, J. C.; Yang, K.; Du, S. X.; Gao, H.-J. J. Phys. Chem. C 2011, 115, 21750− 21754. (16) Zambelli, T.; Jiang, P.; Lagoute, J.; Grillo, S. E.; Gauthier, S.; Gourdon, A.; Joachim, C. Phys. Rev. B 2002, 66, 075410. (17) Pascual, J. I.; Jackiw, J. J.; Kely, K. F.; Conrad, H.; Rust, H.-P.; Weiss, P. S. Phys. Rev. B 2000, 62, 12632−12635. (18) Zhang, L.; Cheng, Z.; Huan, Q.; He, X.; Lin, X.; Gao, L.; Deng, Z.; Jiang, N.; Liu, Q.; Du, S.; Guo, H.; Gao, H.-J. J. Phys. Chem. C 2011, 115, 10791−10796. (19) Atodiresei, N.; Caciuc, V.; Lazic, P.; Blügel, S. Phys. Rev. Lett. 2009, 102, 136809. (20) Lee, K.; Morikawa, Y.; Langreth, D. C. Phys. Rev. B 2010, 82, 155461. (21) Sony, P.; Puschnig, P.; Nabok, D.; Ambrosch-Draxl, C. Phys. Rev. Lett. 2007, 99, 176401. (22) Tonigold, K.; Gross, A. J. Chem. Phys. 2010, 132, 224701− 224710. (23) Simic-Milosevic, V.; Meyer, J.; Morgenstern, K. Phys. Chem. Chem. Phys. 2008, 10, 1916−1920. (24) Fagan, S. B.; Filho, A. G. S.; Lima, J. O. G.; Filho, J. M.; Ferreira, O. P.; Mazali, I. O.; Alves, O. L.; Dresselhaus, M. S. Nano Lett. 2004, 4, 1285−1288. (25) Fagan, S. B.; Girao, E. C.; Filho, J. M.; Filho, A. G. S. Int. J. Quantum Chem. 2006, 106, 2558−2563. (26) Gaspari, R.; Pignedoli, C. A.; Fasel, R.; Treier, M.; Passerone, D. Phys. Rev. B 2010, 82, 041408. (27) Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I. Phys. Rev. Lett. 2004, 92, 246401. (28) Gulans, A.; Puska, M. J.; Nieminen, R. M. Phys. Rev. B 2009, 79, 201105. (29) Kresse, G.; Furthmüller, J. Comput. Mater. Sci. 1996, 6, 15−20. (30) Kresse, G.; Furthmüller, J. Phys. Rev. B 1996, 54, 11169−11186. (31) Blöchl, P. E. Phys. Rev. B 1994, 50, 17953−17979. (32) Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1996, 77, 3865−3868. (33) Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1997, 78, 1396. (34) Chakarova-Käck, S. D.; Schröder, E.; Lundqvist, B. I.; Langreth, D. C. Phys. Rev. Lett. 2006, 96, 146107. (35) Rydberg, H.; Dion, M.; Jacobson, N.; Schröder, E.; Hyldgaard, P.; Simak, S. I.; Langreth, D. C.; Lundqvist, B. I. Phys. Rev. Lett. 2003, 91, 126402. (36) Chakarova-Käck, S. D.; Vojvodic, A.; Kleis, J.; Hyldgaard, P.; Schröder, E. New J. Phys. 2010, 12, 013017−013032. (37) Johnston, K.; Kleis, J.; Lundqvist, B. I.; Nieminen, R. M. Phys. Rev. B 2008, 77, 121404. (38) Johnston, K.; Gulans, A.; Verho, T.; Puska, M. J. Phys. Rev. B 2010, 81, 235428. (39) Toyoda, K.; Nakano, Y.; Hamada, I.; Lee, K.; Yanagisawa, S.; Morikawa, Y. Surf. Sci. 2009, 603, 2912−2922.

(40) Yanagisawa, S.; Lee, K.; Morikawa, Y. J. Chem. Phys. 2008, 128, 244704−144716. (41) Hamada, I.; Lee, K.; Morikawa, Y. Phys. Rev. B 2010, 81, 115452. (42) Kleis, J.; Lundqvist, B. I.; Langreth, D. C.; Schröder, E. Phys. Rev. B 2007, 76, 100201. (43) Kleis, J.; Schröder, E.; Hyldgaard, P. Phys. Rev. B 2008, 77, 205422. (44) Mills, G.; Jónsson, H.; Schenter, G. K. Surf. Sci. 1995, 324, 305− 337. (45) Henkelman, G.; Jónsson, H. J. Chem. Phys. 2000, 113, 9978− 9985. (46) Sykes, E. C. H.; Mantooth, B. A.; Han, P.; Donhauser, Z. J.; Weiss, P. S. J. Am. Chem. Soc. 2005, 127, 7255−7260. (47) Weiss, P. S.; Eigler, D. M. Phys. Rev. Lett. 1993, 71, 3139−3142. (48) Tsai, M.-C.; Muetterties, E. L. J. Am. Chem. Soc. 1982, 104, 2534−2539. (49) Campbell, J. M.; Seimanides, S.; Campbell, C. T. J. Phys. Chem. 1989, 93, 815−826. (50) Berland, K.; Einstein, T. L.; Hyldgaard, P. Phys. Rev. B 2009, 80, 155431. (51) Smoluchowski, R. Phys. Rev. 1941, 60, 661−674. (52) Nørskov, J. K.; Abild-Pedersen, F.; Studt, F.; Bligaard, T. Proc. Natl. Acad. Sci. U.S.A. 2011, 108, 937−943. (53) Johnston, K.; Harmandaris, V. J. Phys. Chem. C 2011, 115, 14707−14717. (54) Syomin, D.; Kim, J.; Koel, B. E.; Ellison, G. B. J. Phys. Chem. B 2001, 105, 8387−8394. (55) McNellis, E. R.; Meyer, J.; Reuter, K. Phys. Rev. B 2009, 80, 205414. (56) Ihm, H.; Ajo, H. M.; Gottfried, J. M.; Bera, P.; Campbell, C. T. J. Phys. Chem. B 2004, 108, 14627−14633. (57) Xu, C.; Tsai, Y.-L.; Koel, B. E. J. Phys. Chem. 1994, 98, 585−593. (58) Morin, C.; Simon, D.; Sautet, P. J. Phys. Chem. B 2003, 107, 2995−3002. (59) Wander, A.; Held, G.; Hwang, R. Q.; Blackman, G. S.; Xu, M. L.; de Andres, P.; Van Hove, M. A.; Somorjai, G. A. Surf. Sci. 1991, 249, 21−34. (60) Delle Site, L.; Sebastiani, D. Phys. Rev. B 2004, 70, 115401. (61) Paniago, R.; Matzdorf, R.; Meister, G.; Goldmann, A. Surf. Sci. 1995, 336, 113−122. (62) Skriver, H. L.; Rosengaard, N. M. Phys. Rev. B 1992, 46, 7157− 7168.

20416

dx.doi.org/10.1021/jp306062s | J. Phys. Chem. C 2012, 116, 20409−20416