Scanning Tunneling Microscopic Observation of Adatom-Mediated

Oct 12, 2009 - Self-Assembled Monolayers at High Coverage ... UniVersity of Denmark, DK-2800 Lyngby, Denmark, and School of Molecular and Microbial ...
0 downloads 0 Views 3MB Size
J. Phys. Chem. C 2009, 113, 19601–19608

19601

Scanning Tunneling Microscopic Observation of Adatom-Mediated Motifs on Gold-Thiol Self-Assembled Monolayers at High Coverage Yun Wang,† Qijin Chi,‡ Noel S. Hush,†,§ Jeffrey R. Reimers,*,† Jingdong Zhang,‡ and Jens Ulstrup*,‡ School of Chemistry, The UniVersity of Sydney, Australia, Department of Chemistry and NanoDTU, Technical UniVersity of Denmark, DK-2800 Lyngby, Denmark, and School of Molecular and Microbial Biosciences, The UniVersity of Sydney, Australia ReceiVed: July 2, 2009; ReVised Manuscript ReceiVed: September 9, 2009

Self-assembled monolayers (SAMs) formed by chemisorption of a branched-chain alkanethiol, 2-methyl-1propanethiol, on Au(111) surfaces were studied by in situ scanning tunneling microscopy (STM) under electrochemical potential control and analyzed using extensive density functional theory (DFT) calculations. The SAM forms in the unusual (8 × 3)-4 superlattice, producing a very complex STM image. Seventy possible structures were considered for the SAM, with the calculated lowest-energy configuration in fact predicting the details of the unusual observed STM image. The most stable structure involves two R-S-Au-S-R adatom-mediated motifs per surface cell, with steric-induced variations in the adsorbate alignment inducing the observed STM image contrasts. Observed pits covering 5.6 ( 0.5% of the SAM surface are consistent with this structure. These results provide the missing link from the structural motifs observed on surfaces at low coverage and on gold nanoparticles to the observed spectroscopic properties of high-coverage SAMs formed by methanethiol. However, the significant role attributed to intermolecular steric packing effects suggests a lack of generality for the adatom-mediated motif at high coverage. 1. Introduction Self-assembled monolayers (SAMs) have been of great interest both experimentally and theoretically as potential applications range from molecular devices to biorecognition to nanotechnology.1 Among all the SAMs studied, alkanethiol monolayers on Au(111) are the most comprehensively investigated due to their strong S-Au interaction and ease of production.2,3 Although numerous studies on thiol-gold SAMs have been performed, some fundamental details of electronic structures and formation processes remain to be clarified, such as the precise atomic-level details of structures adopted at the thiol-gold interface and the factors that control them, factors including both the headgroup interaction and the role of the aliphatic spacer.3 The sulfur-gold (S-Au) junctions that connect the substrate surface to the organic molecules are among the most important aspects of the SAM, with, e.g., the conductance of alkanethiol molecules sandwiched between two gold electrodes depending strongly on the structure.4 While thiol-gold SAMs are usually produced by exposing gold surfaces to alkanethiol or dialkanedisulfide solutions, many studies have demonstrated that the active adsorbates are in fact alkanethiyl radicals formed following loss of the thiol protons or cleavage of disulfide bonds.5-8 The purely physical interactions between alkanethiol molecules and the gold surface are too weak to form stable monolayers,9,10 while the gold surface possesses insufficient reducing power to drive thiolate formation.11 * To whom correspondence should be addressed. E-mail: [email protected] (J.U.); [email protected] (J.R.R.). † School of Chemistry, The University of Sydney. ‡ Technical University of Denmark. § School of Molecular and Microbial Biosciences, The University of Sydney.

Early insight into the nature of the interface between Au(111) and chemisorbed thiol monolayers came from theoretical studies of binding to regular bulk-like (111) surfaces; these suggest that sulfur atoms preferentially locate in the region between the classic fcc and bridge adsorption sites.5,9 However, preferential adsorption to irregular surfaces containing vacancies and/or adatoms has also been predicted.12-17 Much progress has been made experimentally in recent years concerning the structure of SAMs made from alkanethiols. Groundbreaking detailed structural studies on methanethiol SAMs using normal incidence X-ray standing wave (NIXSW) analysis by Roper et al., as well as the scanned-energy and scanned-angle S 2p photoelectron diffraction experiments by Kondoh et al., show that sulfur atoms occupy atop sites,18,19 but insufficient information was obtained to fully identify the interface structure. Some studies proposed that sulfur atoms may be associated with gold adatoms,12-14,20,21 and later this motif was observed in X-ray structures of gold nanoparticles.22 A recent scanning tunneling microscopy (STM) study on thiol-gold SAMs before and after reaction with atomic hydrogen by Kautz et al. also demonstrated that there is one additional gold adatom for every two alkanethiol adsorbate molecules,23 consistent with the adatom-mediated model, and it has also been shown that these adatoms can be liberated to form islands.21 The success of these studies raises the possibility that the adatom-mediated motif is produced in general during SAM formation between gold and alkanethiols. A very similar adatom-mediated motif has been observed for SAMs formed from the bidentate ligand 1,10-phenanthroline,24 but small changes in conditions for these SAMs lead to monolayers above the regular flat surface, and the structure is strongly coverage dependent. Direct observation of the adatom-mediated motif within SAMs formed from thiols has been achieved only at low thiol coverage,13 leaving partially answered the important question of the nature of the interface on surfaces at high

10.1021/jp906216k CCC: $40.75  2009 American Chemical Society Published on Web 10/12/2009

19602

J. Phys. Chem. C, Vol. 113, No. 45, 2009

coverage and hence the interconnection between the lowcoverage, high-coverage, and nanoparticle-coverage scenarios. The gross structural properties at high coverage of gold SAMs formed using linear alkanethiols are well known, these forming in the (3 × 23)-4 superlattice (also termed the c(4 × 2) superlattice), The (3 × 23)-4 superlattice means that there are four adsorbates in each surface cell, with lattice lengths of 3 or 23 times to the shortest Au-Au distance in bulk, respectively. Wang et al. proposed that the four adsorbate molecules in this unit cell attach with two connected in an adatom-mediated motif while two more attach directly to the surface gold layer.16 Alternatively, Gro¨nbeck et al. propose that adatom-bound motifs alone can account for the structure, possibly with or without local surface gold vacancies.17 Finally, Cossaro et al. suggested that the structure of the c(4 × 2) superlattice is quite complex involving directly attached adsorbate molecules, local vacancies, and one-dimensional adatommediated -(S-Au)n- zigzag chains.15 In principle, the observed density of pits on the gold surface3 provides significant information that could be used to assess these differing models, but their concentration is difficult to quantify. Also, while observed STM images do provide critical information, their interpretation is not always unambiguous. Many specific questions thus remain unanswered, and new probes are required. Progress has been made by considering the analogous properties of branched-chain alkanethiol monolayers. Indeed, a recent scanning tunnelling microscopy (STM) study found that a branched alkanethiol chemisorbate on the Au(111) surface formed quite different structural patterns compared to those for linear alkanethiols.25 In the report by Chi et al., the selfassembled monolayer formed from 2-methyl-2-propanethiol is found to adopt the (27 × 7)R19.1°-2 lattice, with two inequivalent radicals in a (27 × 7)R19.1° surface cell.11,25 This structure is strikingly different from the (3 × 23)-4 and similar (e.g., (3 × 3)R30°) structures found for most SAMs formed from straight-chain alkanethiols.3 In particular, there is no surface pitting involved in the SAM formation, and DFT calculations indicate that no adatom-mediated motifs are involved in the binding at all, with the appearance of the unusual double-occupied surface cell attributed to steric repulsions between the aliphatic chains.11 Hence, small changes in the steric interactions between adsorbate molecules can control SAM structure. In this work, we employ electrochemical STM technology (in situ STM) and plane-wave-based density-functional theory (DFT) to investigate the SAMs formed from a systematically chosen molecule, 2-methyl-1-propanethiol, on Au(111). This isomer of 2-methyl-2-propanethiol contains a longer chain section and a secondary carbon rather than the fully branched tertiary isomer studied previously and is intermediary with respect to the well-studied straight-chain alkanethiols. This modification illustrates the extent to which further interchain steric interactions could modify the SAM structure. 2. Methods 2.1. Experimental Measurements. 2-Methyl-1-propanethiol ((CH3)2CHCH2SH, 98%) was obtained from Sigma. The NH4Ac buffer (5 mM, pH 4.6) was prepared from 5 M stock solution (Fluka, ultrapure) and solution pH was adjusted by acetic acid (Aldrich, 99.7%). NaOH (0.1 or 0.5 M) solutions were prepared from 30% stock solution (Merck, ultrapure) and used for voltammetric measurements of reductive desorption. Absolute ethanol (ultrapure) was from Merck. Milli-Q purified water (18.2 MΩ cm) was used throughout. The Au(111) electrodes used in

Wang et al. both electrochemical and STM measurements were homemade and pretreated before use as described.23 The SAMs were prepared by soaking freshly quenched Au(111) electrodes in ethanol solution containing 2-methyl-1propanethiol (1-5 mM) overnight, followed by thorough rinsing with pure ethanol and Milli-Q water before measurements. To observe the dynamic processes of the SAM formation, 2-methyl1-propanethiol was directly added to the electrochemical STM cell. The details were similar to that for 1-propanethiol.24 Electrochemical measurements were performed using an Autolab system (Eco Chemie, Netherlands), with a three-electrode system containing a Pt counter electrode, a reversible hydrogen electrode (RHE) as the reference electrode, and an Au(111) working electrode. STM imaging was carried out in the electrochemical in situ mode using a PicoSPM system (Molecular Imaging Co.) equipped with a bipotentiostat for potential control of both the substrate and tip. The STM tips were prepared from Pt/Ir (80:20) wire by electrochemical etching and covered with Apiezon wax to reduce or eliminate Faradaic currents. 2.2. Computational Modeling. All the computations were conducted using periodic density functional theory (DFT) methods with the Vienna ab initio simulation package (VASP).26 For the electron-electron exchange and correlation interactions, the functional of Perdew and Wang (PW91),27 a form of the general gradient approximation (GGA), was used throughout. Electron-ion interactions were described using the optimized relativistic Vanderbilt-type ultrasoft pseudopotentials.28,29 A plane-wave basis set was employed to expand the smooth part of the wave functions with a kinetic energy cutoff of 300 eV, required to properly describe the carbon atoms. VASP utilized an iterative scheme to solve self-consistently the Kohn-Sham equations using residuum-minimization techniques. The geometric structure was optimized with the conjugated-gradient method. The broadening approach proposed by Methfessel and Paxton with an electronic temperature of 0.2 eV was used for the calculation of orbital occupancies. The Au(111) surface was modeled by a supercell consisting of a four-layer slab separated by a vacuum region equivalent to a six-layer thickness. When the geometry was optimized, the top two atomic layers and the adsorbates were allowed to relax, while the lower two layers were fixed at the ideal bulk-like position (using the lattice constant a0 ) 4.18 Å obtained at the same level of theory30). It has been demonstrated that the fourlayer slab is thick enough to obtain reliable results.11 The methods proposed by Neugebauer and Makov et al. were used to correct the error for the surface with a large dipole moment,31,32 when the energy converges slowly with respect to the z axis of the supercell. All the calculations of the alkanethiyl/ Au(111) systems were performed by ignoring spin polarization as studies33 on similar systems indicate that the adsorbatesubstrate interaction is strong and the effects of spin polarization are small. The Brillouin-zone integrations were performed using Monkhorst-Pack grids of special points, with (1 × 6 × 1) and (4 × 4 × 1) k-points meshes used for the (8 × 3) and (3 × 23) surface cells, respectively. These k-points meshes have been demonstrated to be sufficient for our purposes.33 STM images were simulated using the Tersoff-Hamann approximation34 with a bias voltage of 1.0 V. While the clean Au(111) surface undergoes a (22 × 3) reconstruction,35,36 thiol chemisorption lifts the reconstruction,3,12 and in these equilibrium studies binding only to the unreconstructed surface is considered. Even though experimental results for this system have been obtained in aqueous solution, we

Adatom-Mediated Motifs on Gold-Thiol SAMs

Figure 1. In situ STM images of the 2-methyl-1-propanethiol SAM assembled on a Au(111) surface obtained in 5 mM NH4Ac (pH 4.6), with the substrate lattice vectors (black) and (3 × 23)-4 substrate lattice (blue) indicated. (a) 100 nm × 100 nm, (b) 36 nm × 36 nm, (c) 7.3 nm × 7.3 nm; the blue rectangle also indicates the part of the image that is reproduced in Figure 6. STM imaging conditions: tunnel current (It) 0.25 nA, bias voltage (Vb) -0.60 V, and substrate potential (Ew) 0.26 V (vs SCE).

model only monolayers under vacuum conditions. Since the alkanethiol tail groups are hydrophobic, solvent interactions are most likely of only minor importance, in contrast to the SAMs containing ionic tail groups (such as cysteamine) in which solvation of the charged tail groups is critical.37 The thiyl rather than thiolate nature of the headgroup means that the actual surface charge transfer is small, and so once again solvent effects are not of great significance.11 3. Results and Discussion 3.1. Experimental Observations. In situ STM images observed over different scanned areas are shown in Figure 1. The surface unit cell is rectangular with dimensions of 23.3 Å × 4.8 Å (Figure 1c). Observation of the reconstruction lines of the bare single-crystal gold surface before addition of the

J. Phys. Chem. C, Vol. 113, No. 45, 2009 19603 adsorbate allows the absolute orientation of the surface-cell vectors to be determined, indicating that the SAM forms in the (8 × 3) lattice and contains 16 gold atoms per regular surface layer in the unit cell. Electrochemical measurements of the surface coverage using reductive desorption (Supporting Information) give a density of 5.9 ( 0.4 × 10-10 mol cm-2, indicating a surface coverage (denoted in the number of adsorbate molecules per Au(111) regular-surface gold atom) of Θ ) 0.26 ( 0.02 ∼ 1/4. As expected, this coverage is intermediate between those found for the 2-methyl-2-propanethiol and 1-propanethiol SAMs of 0.14 and 0.33, respectively. Four adsorbate molecules thus occupy each surface unit cell, possibly with four different internal structures. Indeed, four regions in each unit cell are observed in the in situ STM image (Figure 1c) with strongly varying contrast. Some of these features are widely spaced while others are spaced less than the minimumpossible S-S intra-adsorbate distance, suggesting that each adsorbate can produce more than one STM feature. Hence, very many adsorbate configurations could, in principle, generate this image. The most striking features of the image are two bright spots very close to each other and two dimer spots separated by a large dark region; also, the spots are not collinear. 3.2. Computational Modeling of the Monolayer Structures. Full details of all 70 computed structures for Au(111) SAMs made from 2-methyl-1-propanethiol are provided in the Supporting Information, including calculated STM images for 11 structures. Owing to the multitude of possible structures, we perform systematic investigations, starting with the nature of adsorption of noninteracting or partially interacting molecules on the flat surface (i.e., low to medium coverage). Then we increase the coverage and investigate SAM formation on the regular flat surface Au(111). Structures are then considered in which the gold atoms can rearrange within the unit cell, forming adatoms and vacancies, but transport of gold atoms to or from the unit cell is not allowed. Finally, the number of gold atoms in each surface unit cell is varied, allowing for gold adatoms above regular surfaces and possible associated pit formation. In this way, a range of low-energy structures can be identified and their suitability is further characterized in terms of their relative energy and ensuing STM images. 3.2.1. General Features of Chemisorptions on a Flat Au(111) Surface with Low to Medium CoWerage. On a regular Au(111) surface, there are four kinds of high-symmetry adsorption sites: atop, bridge, fcc, and hcp sites. Previous computational studies have demonstrated that sulfur atoms of chemisorbed alkanethiols prefer low-symmetry locations either between the bridge and fcc sites, termed FB, or between the bridge and hcp sites, termed HB.11 Hence, we calculate the adsorption energy at atop, HB, and FB sites only. Initially we consider adsorbates with similar alkane-chain conformations at sufficiently low coverage, Θ ) 1/12 ) 0.083, to avoid intermolecular steric interactions on the (3 × 23) supersurface lattice. For comparison we also consider an alternate lattice, (8 × 3), that gives a slightly lower coverage of Θ ) 1/16 but, because of its short unit-cell length in the [112j] direction, does actually display significant intermolecular steric interactions. The energy of adsorption (∆E) is determined from

∆E ) Etot - Esurf - Ead

(1)

where Etot and Esurf are the optimized total energies of the systems with and without adsorbate, respectively, and Ead represents the energy of the optimized isolated 2-methyl-1propanethiyl radical calculated allowing for spin polarization.

19604

J. Phys. Chem. C, Vol. 113, No. 45, 2009

Wang et al.

Figure 2. Elevation view of four DFT-optimized configurations A-D for the 2-methyl-1-propanethiyl radical adsorbed on the regular Au(111) surface obtained at the low coverage of Θ ) 0.083 using the (3 × 23) substrate lattice: gold, Au; green, S; gray, C; white, H.

TABLE 1: Adsorption Energies ∆E and Geometrical Properties for Single 2-Methyl-1-propanethiyl Radicals Adsorbed on (3 × 23)R30o or (8 × 3)R30o Au(111) Supersurface Lattices in Different Configurations (see Figure 2)a surface cell structure ∆E (eV)

dS-Au (Å)

(3 × 23)

2.46, 2.47, 2.45, 2.47, 2.47, 2.46, 2.48, 2.47,

(8 × 3)

A B C D A B C D

-1.70 -1.65 -1.71 -1.71 -1.62 -1.58 -1.78 -1.72

2.59 2.52 2.56 2.53 2.59 2.56 2.58 2.62

dS-C (Å) θ (degree) 1.87 1.86 1.86 1.87 1.88 1.86 1.87 1.87

21 7 9 21 23 10 7 19

a dS-Au is the shortest distance from a sulfur atom to the surface gold atom, ds-c is the C-S bond length, and θ is the tilt angle of S-C bond to the surface normal.

The most stable site is predicted to be the FB site, this being 0.16 eV more stable than the HB site and 0.34 eV more stable than the atop site, similar to the results obtained in previous theoretical studies of chemisorbed alkanethiol monolayers on regular Au(111) surfaces.9,12,38 A range of possible conformations are available at each adsorption site of the flat Au(111) surface for the adsorbate molecules, originating from the degrees of freedom afforded by the rotation and tilt angles of the aliphatic chain. However, the energy differences associated with these conformational variations are calculated to be much smaller than those associated with the variations in the surface binding sites. Four lowenergy structures, A-D, are shown in Figure 2 above the FB site on the (3 × 23) and (8 × 3) surface cells; their energetic and structural properties are listed in Table 1, while full details are provided in the Supporting Information. For the (3 × 23) surface cell, structures A, C, and D differ in energy by only 0.01 eV at ∆E ≈ -1.71 eV per adsorbate molecule, while structure B is just 0.06 eV (2 kT at room temperature) less stable. Hence, we see that at low coverage the substrate surface to headgroup interaction does not fully govern the conformations of the adsorbate molecules on the flat surface, and so other factors such as intermolecular interactions and solvation effects should affect the SAM conformations more significantly. The results provided for the lower-coverage (3 × 23) surface cell

confirm this analysis, with configurations C and D showing enhanced binding (∆E ) -1.78 eV per adsorbate molecule for C and -1.72 eV for D) owing to improved substrate relaxation. Configurations A and B are destabilized by steric interactions in the short [112j] direction (∆E ) -1.62 eV per adsorbate molecule for A and -1.58 eV for B). For 2-methyl-1propanethiol chemisorbates, steric interactions also destabilize the (3 × 23)-4 and similar lattices in which the SAMs of most straight-chain alkanethiols form, with the calculated energy of adsorption on e.g. the (3 × 3)R30° lattice being calculated to be just -1.32 eV. These results show the importance of substrate and steric interactions at low to medium coverage, effects that are expected to become more significant when four adsorbate molecules are crammed into the (8 × 3) surface cell at the high coverage (Θ ) 0.25) actually observed for stable SAMs. Considering the calculated structural properties in detail (Table 1), we see that those for 2-methyl-1-propanethiol chemisorbed monolayers directly parallel those for related molecules such as methanethiol and 2-methyl-2-propanethiol.11 After adsorption, the Au-Au bond lengths at the adsorption site increase from their optimized bulk value of 2.96 Å to over 3.50 Å while the calculated S-Au bond lengths of 2.47 Å also indicate a strong metal-molecule interaction. Applying the Bader method which decomposes the space charge density into atomic contributions to our system,39,40 it is found that only 0.07e is transferred from the substrate to the adsorbate when 2-methyl1-propanethiol chemisorbs at the FB site, where e is the magnitude of the charge on the electron. Changes of ca. 0.2e have been found in experimental studies41-43 of thiol chemisorption. While in general the significant band lineup error that occurs in DFT for molecular adsorbates on metal surfaces results in an underestimation of calculated charge flows,33 the results again clearly indicate that ionic thiolate states are not involved in the chemisorption. Thus, the adsorbate forms as a thiyl species identified as the 2-methyl-1-propanethiyl radical. 3.2.2. High-CoWerage Adsorption on the Flat Au(111) Surface. Here we consider the possibility that dense SAMs made from 2-methyl-1-propanethiol form on the flat Au(111) surface, as has been observed by STM for related SAMs.2 The most significant aspect of the observed STM images of 2-methyl-1propanethiol SAMs (Figure 1) is the starkly different separations and contrasts of the features. A similar scenario was also found for chemisorbed monolayers of 2-methyl-2-propanethiol on Au(111), a SAM for which more than one molecule per surface cell was observed to form on the flat Au(111) surface.25 For that molecule, variations in the tilt angles of S-C bonds on the surface were found to account for the STM-image topology;11 our initial computational studies of the high-coverage monolayer test the hypothesis that similar variations also account for the observed STM images for 2-methyl-1-propanethiol SAMs. A wide variety of DFT geometry optimizations were performed starting with conformations similar to C or D from Figure 2, the only conformations shown feasible in the low-density studies on the observed (8 × 3) lattice; full details of the results are provided in the Supporting Information. All of the stable structures found are of similar energy, with ∆E varying by less than 0.05 eV. However, the energy of adsorption per radical is -1.46 eV for the most stable case, this being 0.32 eV less stable than that at the lower coverage reported in Table 1. The repulsion between the adsorbates at high coverage is therefore rather strong and could be manifest through either direct steric interactions or significant distortions induced in the substrate lattice.

Adatom-Mediated Motifs on Gold-Thiol SAMs

Figure 3. DFT-optimized structure for the SAM of 2-methyl-1propanethiol chemisorbed to the regular Au(111) surface (configuration b from Figure 4) at the observed coverage of Θ ) 0.25 in the (8 × 3)-4 substrate lattice. Top, elevation view highlighting the configuration (see Figure 2) of each adsorbate molecule (gold, Au; green, S; gray, C; white, H); middle, plan view with vectors indicating the regular Au(111) surface with the location of the adsorbate S atoms (green) superimposed; bottom, simulated STM image.

The simulated STM image of the most stable configuration found is shown in Figure 3; of its four adsorbate molecules, one adopts configuration C while the others adopt configuration D. All the bright spots in the calculated STM image originate from the vertically aligned terminal methyl groups; thermally induced changes in the conformational angles will induce some variation in the calculated image, but the main features are expected to be preserved. Since the tilt angles of the S-C bond are different, the heights of the terminal methyl groups vary (see Figure 3), and therefore, the STM image contrast is modulated. However, the simulated STM image is significantly different from the observed in situ image shown in Figure 1c, with just three resolved features per cell and no large intense closely spaced image pair. Furthermore, owing to the steric interactions associated with the dimensions of the branched chain, no alternate configuration can be envisaged in which vertically aligned chains from two molecules are separated by the observed short distance between the spots of the bright pair. Thus, the hypothesis that four adsorbate molecules assemble on a regular Au(111) surface, akin to the binding observed following chemisorptions of the isomer 2-methyl-2-propanethiol, cannot be supported. 3.2.3. High-CoWerage Adsorption with Adatom-MediatedBonding Motifs aboWe the Au(111) Surface. Recently, several groups have found, both experimentally and theoretically, that the chemisorption of linear alkanethiols occurs at supersurface gold adatoms12-17 that may or may not be associated with surface-layer vacancies; in this binding motif, two adsorbate molecules are always attached to one gold adatom. As four adsorbate chains occupy the unit cell for 2-methyl-1-propanethiol SAMs, either zero, one, or two of these motifs may be formed, and as there exists many sites for each adatom and possible vacancy, a large number of configurations are possible. As a starting point we consider only the structures in which the number of adatoms equals the number of vacancies; for these structures, the number of gold atoms per simulation cell is the same as that for adsorption on the flat Au(111) surface considered previously and so the results may be directly compared. We have examined 18 possible new configurations,

J. Phys. Chem. C, Vol. 113, No. 45, 2009 19605

Figure 4. DFT-optimized structures and relative energies for the SAM of 2-methyl-1-propanethiol chemisorbed to the regular and irregular Au(111) surfaces at the observed coverage of Θ ) 0.25 in the (8 × 3)-4 substrate lattice. The lines indicate interatomic vectors on the clean (111) surface, while white circles indicate surface vacancies, gold indicate supersurface gold adatoms, and green indicates the S atoms of the adsorbate molecules. For configuration j, some adsorbate features are reproduced in the neighboring unit cell to accentuate the linearity of the structure; only the lower-energy all-trans adatom-mediated structures are shown. All energies are expressed with respect to that for configuration a (total interaction energy ∆E ) -7.15 eV per cell).

optimizing the energy of each structure using DFT; the results for 10 illustrative examples, named configurations a-j, are sketched in Figure 4, while all results are reported in the Supporting Information. Configuration a is the lowest-energy structure found, and the energies for the remaining (8 × 3)-4 lattices are expressed relative to its energy. This configuration has two adatom motifs and therefore two vacancies that are found to be located nearly underneath the adatoms. Some variants of this structure are shown in configurations d-i containing different locations for the adatom sites and different arrangements of the vacancies. As some of these configurations are less than 0.3 eV higher in energy than configuration a, the adsorbed molecules and the vacancies may be quite mobile on the surface. The lowest energy configuration involving no adatoms or vacancies (previously detailed in Figure 3) is also sketched in Figure 4 for comparison with configuration b. Its average energy of adsorption per adsorbed molecule is -1.46 eV, much less than the single configuration energies reported in Table 1, owing to the significant steric interactions between the four adsorbate molecules. This unit-cell configuration is 1.33 eV per cell less stable than configuration a, indicating that adsorption on a flat Au(111) surface will spontaneously induce rearrangement of the surface gold atoms. Much of the strain energy apparent in b is relieved through the formation of just a single adatomvacancy pair in the surface cell, with configuration c being only 0.24 eV per cell less stable than the double adatom-vacancy structure, configuration a. The possibility that mixed configurations of this type are also important for linear alkanethiol SAMs has been raised by Wang et al.16 Recently, Cossaro et al. proposed that one-dimensional -S-Au-S-Au-S- chains form the SAMs of long-chain linear alkanethiol SAMs on gold with the (3 × 23) surface cell.15 This precise structure is not feasible on the (8 × 3) surface cell owing to the considerably different unit-cell

19606

J. Phys. Chem. C, Vol. 113, No. 45, 2009

Wang et al.

Figure 5. Atomic structures for two adatom-mediated-bonding motifs adsorbed on flat Au(111) in the configuration a from Figure 4: all-trans isomer; right, all-cis isomer.

parameters, but a somewhat analogous structure is reported as configuration j. This configuration involves three adatom/ vacancy pairs per unit cell but is 1.44 eV per cell less stable than configuration a and hence should not be favored within 2-methyl-1-propanethiol SAMs. In general, an important feature of the adatom-bound adsorbate structure is whether the S-C bonds are located cis or trans with respect to each other. In a previous study of this adatom motif within methanethiol SAMs, the cis structure has been found to have the lowest energy.17 However, in the present case all structures reported in Figure 4 have lower-energy trans configurations; Figure 5 compares the atomic positions and relative energies for the cis and trans variants of configuration a. The alternative result for 2-methyl-1-propanethiyl adsorbate could arise from the different interadduct steric interactions; the trans configuration was also observed in the adatom-mediated motifs in a p-mercaptobenzoic acid (p-MBA)-protected gold nanoparticle22 and hence may in general be preferred for bulkier adsorbate molecules. The adatom-mediated structures have their Au-S bonds orientated nearly parallel to the surface, with the S-adatom-S angle being 173° for the lowest-energy structure, configuration a. However, the C-S-Au angles show significant variations between 109° and 116° for configuration a. Hence the C-S vectors rise at 52°-62° to the surface normal, quite different to the nearly parallel and nearly vertical orientations apparent in Figure 2 for molecules bound directly to Au(111). The variation in the angles is caused by the steric interactions which force the two adatom-mediated motifs to adopt slightly different structures, giving rise to a variation in height of terminal carbon atoms of 0.25 Å. Figure 6 shows a simulated STM image for this structure where it is compared to an extract of the observed in situ STM image from Figure 1c: all major qualitative features of the experimental image are reproduced. Also shown in Figure 6 are side and top views of configuration a, with the adsorbate chains numbered i-iv. This reveals that all of the bright spots originate from the terminal methyl groups, with the major factor controlling their relative brightness being the height of the carbon atoms above the surface. As each adsorbed molecule orients with two methyl groups upward, a total of eight spots of this type could in principle be found in the simulated STM images. The large central bright region of the computed STM image is assigned to the two methyl groups from molecule (ii), with the close proximity of the two groups emphasizing their significance. Two nearby resolved features arise from methyl groups of molecules (iii) and (i), while a further weak feature arises from molecule (iv). Together, these results account for all of the features in the observed in situ STM image. In addition, STM images for five alternative configurations have been simulated. Of these only the image for configuration c is in realistic agreement with the observed STM image (see the Supporting Information). Further, as the relative energy of

Figure 6. DFT-optimized structures (top panels, configuration a plan and elevation) and STM images (lower panels, configurations a, a′, and a′′ and experiment extracted from Figure 1c) for the SAM of 2-methyl-1-propanethiol chemisorbed to the regular and irregular Au(111) surfaces at the observed coverage of Θ ) 0.25 in the (8 × 3)-4 substrate lattice. The lines on the plan view indicate interatomic vectors on the clean (111) surface while white circles indicate surface vacancies, gold indicate supersurface gold adatoms, and green indicates the S atoms of the adsorbate molecules. Configurations a′ and a′′ differ from a only in that either one (a′) or two (a′′) of the surface gold vacancies are filled from an external atom source.

this structure is only 0.24 eV more than that for configuration a, it warrants further consideration. The lowest-energy alternative configuration, configuration f generates a much poorer STM image but is also retained at this stage as a possible option. All the structures considered so far contain the same number of gold atoms in each unit cell as does a flat Au(111) surface. However, Figure 4 indicates that the gold atoms may be quite mobile, requiring further refinements to be made to the optimized structures. Indeed, the clean gold surface to which the adsorbate molecules initially bind is not flat Au(111) but instead is a reconstructed surface in which 46 surface gold atoms (in two

Adatom-Mediated Motifs on Gold-Thiol SAMs rows of 23 atoms along the [110] direction) occupy only 44 (rather than 46) bulk-like positions.36,44 After formation of the SAM, this reconstruction is lifted, releasing 4.5% of the atoms in a gold layer. These excess gold atoms must either be consumed as adatoms or else merge generating new surface terraces. While terrace growth is difficult to positively identify, the opposing process of pit formation is often observed as a common feature in many thiol-gold SAMs including the present case.3 Specifically, the coverage of pits on the surface for the SAMs made from 2-methyl-1-propanethiol is estimated to be 5.6% ( 0.5% of the total surface. This value is slightly larger than that observed for linear 1-propanethiol SAMs, 4.0% ( 0.4%.35 If the gold atoms liberated by lifting of the reconstruction combine with those atoms mined from the pits, then the number of total gold atoms available as adatoms above a vacancy-free surface would increase to 10.1% ( 0.5% for the 2-methyl-1-propanethiol case, close to the adatom coverage of configuration a, 12.5%. In addition, gold atoms mined from terrace edges could also form adatoms, but this effect is difficult to quantify. We have also calculated an additional 43 configurations containing either one or no vacancies in the topmost gold layer. They include all possible variants of configurations a, c, and f, i.e., the configurations previously identified as the ones of greatest interest. All results are given in the Supporting Information, while the four lowest-energy structures found are shown in Figure 4; these four structures are those obtained by adding one Au atom (configuration a′) or two Au atoms (configuration a′′) to configuration a, one Au atom (configuration c′) to configuration c, and two Au atoms (configuration f′′) to configuration f. In these calculations, the average energy of a surface gold atom taken from a pit is estimated as one-quarter of the energy difference between a flat surface and that obtained by removing two adjacent rows (4 atoms) from the original 8 rows in the (8 × 3) surface layer, as described in the Supporting Information; the cost is expected to be slightly less if the filling atom originates from the lifting of the gold reconstruction rather than from surface mining. Filling the vacant atoms is thus found to be an exothermic process, with configurations a′′ being 0.10 eV, more stable than configuration a. Shown in Figure 6 are the simulated STM images for configurations a′ and a′′, the images that are almost identical to that for configuration a considered earlier. STM images are thus insensitive to local pitting under the SAM, and the predicted image for configuration a′′ is in excellent agreement with experimental observation. The calculated image for the nextlowest energy structure, configuration f′′, is in poor agreement with the observed image (see Supporting Information). 4. Conclusions The structure of the SAM formed by chemisorption of 2-methyl-1-propanethiol on the Au(111) surface has been investigated comprehensively both by electrochemical in situ STM measurements and by a priori DFT calculations. In addition to providing visualization of the SAM at the singlemolecule level in situ, STM establishes a means for selection from among the large number of surface-molecule and intermolecular adsorbate interactions conceived by the computations. This complementary approach resembles the complementary information offered jointly by in situ STM and electrochemistry at single-crystal electrode surfaces, with, for example, the adsorbate coverage measured by electrochemical reductive desorption significantly supporting the interpretation of the in situ STM image.37,45

J. Phys. Chem. C, Vol. 113, No. 45, 2009 19607 For 2-methyl-1-propanethiol on the Au(111) surface, an adatom-mediated binding motif is deduced. This adatommediated motif is now well established for SAMs at low coverage12-14 and for nanoparticles,22 but this is in fact the first positive identification of such a structure on a flat surface at high coverage. As key aspects of this structure have also been observed for methanethiol SAMs,18,19 the adatom motif may indeed be common. However, our results also show clearly that a wide range of alternate structures are possible in the sense of showing only very small energy differences. Steric interactions between molecules and adsorbate relaxation thus provide significant influences. STM cannot resolve variations of the basic motif involving surface-layer vacancies, but the density of surface pits can be mapped through careful measurements. This may also call for reconsideration of the effects of the solvent environment even though these effects are likely to be small for the alkanethiol systems. In any case, we see that the monolayer structure adopts the common adatom-mediated motif not because of some inherent stability but rather as a result of competition between many opposing forces. Taken together with the comprehensive STM/DFT/electrochemical approach used in our previous studies,25,35,37,45-47 a highly diverse pattern in both Au-S binding modes and surface packing of short linear and branched similar-size alkanethiols actually emerges. The straight-chain alkanethiol, 1-propanethiol, was shown to bind on a flat Au(111)-electrode surface and pack into a (3 × 3)R30° surface lattice. Extensive Au-atom surface dynamics and pitting accompanies this surface adlayer formation.35 The fully branched 2-methyl-2-propanethiol (tertbutanethiol) instead adopt a strikingly different (27 × 7)R19.1° surface cell with no surface pitting. Finally, 2-methyl-1-propanethiol (secondary butanethiol) addressed here adopts a third thiol binding and packing mode at high coverage, i.e., a (8 × 3)-4 surface lattice with two gold-atom-mediated surface entities per unit cell as a dominating binding motif. Chemically very similar adsorbates can thus produce strikingly different SAMs. Acknowledgment. We thank the Australian Research Council for funding this work under the Discovery Grant program. The use of the supercomputer facilities at the Australia Partnership for Advanced Computing (APAC) is gratefully acknowledged. Financial support from the Danish Research Council for Technology and Production Sciences (Contract No. 274-070272) and the NABIIT programme of the Danish Strategic Research Council (Contract No. 2106-04-0039) is acknowledged. Supporting Information Available: Optimized coordinates for 70 SAMs, 4 substrate lattices, the adsorbate radical along with 11 simulated STM images of SMAs, and a linear voltammogram to show reductive desorption of 2-methyl-1propanethiol. This material is available free of charge via the Internet at http://pubs.acs.org. References and Notes (1) Love, J. C.; Estroff, L. A.; Kriebel, J. K.; Nuzzo, R. G.; Whitesides, G. M. Chem. ReV. 2005, 105, 1103. (2) Vericat, C.; Vela, M. E.; Benitez, G. A.; Gago, J. A. M.; Torrelles, X.; Salvarezza, R. C. J. Phys.: Condens. Matter 2006, 18, R867. (3) Vericat, C.; Vela, M. E.; Salvarezza, R. C. Phys. Chem. Chem. Phys. 2005, 7, 3258. (4) Muller, K. H. Phys. ReV. B 2006, 73, 045403. (5) Vargas, M. C.; Giannozzi, P.; Selloni, A.; Scoles, G. J. Phys. Chem. B 2001, 105, 9509. (6) Kodama, C.; Hayashi, T.; Nozoye, H. Appl. Surf. Sci. 2001, 169, 264.

19608

J. Phys. Chem. C, Vol. 113, No. 45, 2009

(7) Hayashi, T.; Kodama, C.; Nozoye, H. Appl. Surf. Sci. 2001, 169, 100. (8) Liu, G.; Rodriguez, J. A.; Dvorak, J.; Hrbek, J.; Jirsak, T. Surf. Sci. 2002, 505, 295. (9) Yourdshahyan, Y.; Rappe, M. A. J. Chem. Phys. 2002, 117, 825. (10) Zhou, J. G.; Hagelberg, F. Phys. ReV. Lett. 2006, 97, 45505. (11) Wang, Y.; Hush, N. S.; Reimers, J. R. J. Phys. Chem. C 2007, 111, 10878. (12) Wang, Y.; Hush, N. S.; Reimers, J. R. J. Am. Chem. Soc. 2007, 129, 14532. (13) Maksymovych, P.; Sorescu, D. C.; Yates, J. T., Jr. Phys. ReV. Lett. 2006, 97, 146103. (14) Mazzarello, R.; Cossaro, A.; Verdini, A.; Rousseau, R.; Casalis, L.; Danisman, M. F.; Floreano, L.; Scandolo, S.; Morgante, A.; Scoles, G. Phys. ReV. Lett. 2007, 98, 16102. (15) Cossaro, A.; Mazzarello, R.; Rousseau, R.; Casalis, L.; Verdini, A.; Kohlmeyer, A.; Floreano, L.; Scandolo, S.; Morgante, A.; Klein, M. L.; Scoles, G. Science 2008, 321, 943. (16) Wang, J. G.; Selloni, A. J. Phys. Chem. C 2007, 111, 12149. (17) Gro¨nbeck, H.; Ha¨kkinen, H.; Whetten, R. L. J. Phys. Chem. C 2008, 112, 15940. (18) Roper, M. G.; Skegg, M. P.; Fisher, C. J.; Lee, J. J.; Dhanak, V. R.; Woodruff, D. P.; Jones, R. G. Chem. Phys. Lett. 2004, 389, 87. (19) Kondoh, H.; Iwasaki, M.; Shimada, T.; Amemiya, K.; Yokoyama, T.; Ohta, T. Phys. ReV. Lett. 2003, 90, 066102. (20) Yu, M.; Bovet, N.; Satterley, C. J.; Bengio, S.; Lovelock, K. R. J.; Milligan, P. K.; Jones, R. G.; Woodruff, D. P.; Dhanak, V. Phys. ReV. Lett. 2006, 97, 166102. (21) Li, F. S.; Zhou, W. C.; Guo, Q. M. Phys. ReV. B 2009, 79, 113412. (22) Jadzinsky, P. D.; Calero, G.; Ackerson, C. J.; Bushnell, D. A.; Kornberg, R. D. Science 2007, 318, 430. (23) Kautz, N. A.; Kandel, S. A. J. Am. Chem. Soc. 2008, 130, 6908. (24) Cafe, P. F.; Larsen, A. G.; Yang, W.; Bilic, A.; Blake, I. M.; Crossley, M. J.; Zhang, J.; Wackerbarth, H.; Ulstrup, J.; Reimers, J. R. J. Phys. Chem. C 2007, 111, 17285. (25) Chi, Q.; Zhang, J.; Ulstrup, J. J. Phys. Chem. B 2006, 110, 1102. (26) Kresse, G.; Hafner, J. Phys. ReV. B 1993, 47, RC558.

Wang et al. (27) Perdew, J. P.; Wang, Y. Phys. ReV. B 1992, 45, 13244. (28) Vanderbilt, D. Phys. ReV. B 1990, 41, 7892. (29) Kresse, G.; Hafner, J. J. Phys.: Condens. Matter 1994, 6, 8245. (30) Bilic, A.; Reimers, J. R.; Hush, N. S.; Hafner, J. J. Chem. Phys. 2002, 116, 8981. (31) Neugebauer, J.; Scheffler, M. Phys. ReV. B 1992, 46, 16067. (32) Makov, G.; Payne, M. C. Phys. ReV. B 1995, 51, 4014. (33) Bilic, A.; Reimers, J. R.; Hush, N. S. J. Chem. Phys. 2005, 122, 094708. (34) Tersoff, J.; Hamann, D. R. Phys. ReV. Lett. 1983, 50, 1998. (35) Zhang, J.; Chi, Q.; Ulstrup, J. Langmuir 2006, 22, 6203. (36) Wang, Y.; Hush, N. S.; Reimers, J. R. Phys. ReV. B 2007, 75, 233416. (37) Zhang, J.; Bilic, A.; Reimers, J. R.; Hush, N. S.; Ulstrup, J. J. Phys. Chem. B 2005, 109, 15355. (38) Maksymovych, P.; Sorescu, D. C.; Yates, J. T. J. Phys. Chem. B 2006, 110, 21161. (39) Henkelman, G.; Arnaldsson, A.; Jonsson, H. Comput. Mater. Sci. 2006, 36, 354. (40) Bader, R. Atoms in molecules: A Quantum Theory; Oxford University Press: New York, 1990. (41) Azzam, W.; Wehner, B. I.; Fischer, R. A.; Terfort, A.; Wo¨ll, C. Langmuir 2002, 18, 7766. (42) Bourg, M.-C.; Badia, A.; Lennox, R. B. J. Phys. Chem. B 2000, 104, 6562. (43) Duwez, A.-S. J. Elect. Spectrosc. Rel. Phen. 2004, 134, 97. (44) Min, B. K.; Alemozafar, A. R.; Biener, M. M.; Biener, J.; Friend, C. M. Top. Catal. 2005, 36, 77. (45) Nazmutdinov, R. R.; Zhang, J.; Zinkicheva, T. T.; Manyurov, I. R.; Ulstrup, J. Langmuir 2006, 22, 7556. (46) Wang, Y.; Hush, N. S.; Reimers, J. R. J. Am. Chem. Soc. 2007, 129, 14532. (47) Zhang, J.; Chi, Q.; Nazmutdinov, R. R.; Zinkicheva, T. T.; Bronshtein, M. D. Langmuir 2009, 25, 2232.

JP906216K